A wake model is pursued for potential flow past a submerged, finite-length plate that is perpendicular to a uniform, horizontal stream bounded above by a free surface. The effects of gravity are included along the free surface. The approach is to adopt an open-wake model such that the wake boundaries become parallel to the undisturbed stream at some (unknown) point downstream. Boundary integral equations are formed and then discretised along the wake boundaries and free surface in order to obtain a solution numerically. In terms of the dependency of the solution on various parameters, the problem will be formulated in two ways. First, for a given Froude number and ratio of the length of the vertical plate to the draft (the depth of the bottom of the vertical plate relative to the undisturbed free surface), the effect of the wake underpressure coefficient on the size of the wake will be considered. Then, the problem will be discussed where we instead (more naturally) fix the Froude number, draft and length of the vertical, submerged plate. The dependencies of the solution on these parameters will be analysed regarding the effects on several factors, including the size of the wake, the relative lengths of the upper and lower wake boundaries, and the resulting wake underpressure coefficient.

This article offers an advanced and novel investigation into the intricate propagation dynamics of the Belousov–Zhabotinsky system with non-local delayed interaction, which exhibits dynamical transition structure from bistable to monostable. We first solved the enduring open problem concerning the existence, uniqueness and the speed sign of the bistable travelling waves. In the monostable case, we developed and derived new results for the minimal wave speed selection, which, as an application, further improved the existing investigations on pushed and pulled wavefronts. Our results can provide new estimate to the minimal speed as well as to the determinacy of the transition parameters. Moreover, these results can be directly applied to standard localised models and delayed reaction diffusion models by choosing appropriate kernel functions.

In this work, we study early warning signs for stochastic partial differential equations (SPDEs), where the linearisation around a steady state is characterised by continuous spectrum. The studied warning sign takes the form of qualitative changes in the variance as a deterministic bifurcation threshold is approached via parameter variation. Specifically, we focus on the scaling law of the variance near the transition. Since we are dealing here, in contrast to previous studies, with the case of continuous spectrum and quantitative scaling laws, it is natural to start with linearisations of the drift operator that are multiplication operators defined by analytic functions. For a one-dimensional spatial domain, we obtain precise rates of divergence. In the case of the two- and three-dimensional domains, an upper bound to the rate of the early warning sign is proven. These results are cross-validated by numerical simulations. Our theory can be generically useful for several applications, where stochastic and spatial aspects are important in combination with continuous spectrum bifurcations.

We study the community detection problem on a Gaussian mixture model, in which vertices are divided into $k\geq 2$ distinct communities. The major difference in our model is that the intensities for Gaussian perturbations are different for different entries in the observation matrix, and we do not assume that every community has the same number of vertices. We explicitly find the necessary and sufficient conditions for the exact recovery of the maximum likelihood estimation, which can give a sharp phase transition for the exact recovery even though the Gaussian perturbations are not identically distributed; see Section 7. Applications include the community detection on hypergraphs.

The purpose of the paper is twofold. First, we show that partial-data transmission eigenfunctions associated with a conductive boundary condition vanish locally around a polyhedral or conic corner in $\mathbb{R}^n$, $n=2,3$. Second, we apply the spectral property to the geometrical inverse scattering problem of determining the shape as well as its boundary impedance parameter of a conductive scatterer, independent of its medium content, by a single far-field measurement. We establish several new unique recovery results. The results extend the relevant ones in [26] in two directions: first, we consider a more general geometric setup where both polyhedral and conic corners are investigated, whereas in [26] only polyhedral corners are concerned; second, we significantly relax the regularity assumptions in [26] which is particularly useful for the geometrical inverse problem mentioned above. We develop novel technical strategies to achieve these new results.

We address the problem of optimal transport with a quadratic cost functional and a constraint on the flux through a constriction along the path. The constriction, conceptually represented by a toll station, limits the flow rate across. We provide a precise formulation which, in addition, is amenable to generalization in higher dimensions. We work out in detail the case of transport in one dimension by proving existence and uniqueness of solution. Under suitable regularity assumptions, we give an explicit construction of the transport plan. Generalization of flux constraints to higher dimensions and possible extensions of the theory are discussed.

In a smoothly bounded domain $\Omega \subset \mathbb{R}^n$, $n\ge 1$, this manuscript considers the homogeneous Neumann boundary problem for the chemotaxis system

\begin{eqnarray*} \left \{ \begin{array}{l} u_t = \Delta u - \nabla \cdot (u\nabla v), \\[5pt] v_t = \Delta v + u - \alpha uv, \end{array} \right . \end{eqnarray*}

$\alpha \gt 0$

It is shown that there exists $\delta _\star =\delta _\star (n)\gt 0$ such that for any given $\alpha \ge \frac{1}{\delta _\star }$ and for any suitably regular initial data satisfying $v(\cdot, 0)\le \delta _\star$, this problem admits a unique classical solution that stabilizes to the constant equilibrium $(\frac{1}{|\Omega |}\int _\Omega u(\cdot, 0), \, \frac{1}{\alpha })$ in the large time limit.

This work describes a hyperbolic model for cell-cell repulsion with population dynamics. We consider the pressure produced by a population of cells to describe their motion. We assume that cells try to avoid crowded areas and prefer locally empty spaces far away from the carrying capacity. Here, our main goal is to prove the existence of travelling waves with continuous profiles. This article complements our previous results about sharp travelling waves. We conclude the paper with numerical simulations of the PDE problem, illustrating such a result. An application to wound healing also illustrates the importance of travelling waves with a continuous and discontinuous profile.

Given a set of standard binary patterns and a defective pattern, the pattern retrieval task is to find the closest pattern to the defective one among these standard patterns. The Hebbian network of Kuramoto oscillators with second-order coupling provides a dynamical model for this task, and the mutual orthogonality in memorised patterns enables us to distinguish these memorised patterns from most others in terms of stability. For the sake of error-free retrieval for general problems lacking orthogonality, a unified approach was proposed which transforms the problem into a series of subproblems with orthogonality using the orthogonal lift for two patterns. In this work, we propose the least orthogonal lift for three patterns, which evidently reduces the time of solving subproblems and even the dimensions of subproblems. Furthermore, we provide an estimate for the critical strength for stability/instability of binary patterns, which is convenient in practical use. Simulation results are presented to illustrate the effectiveness of the proposed approach.

We introduce a new non-abelian quantum synchronisation model over the unitary group, represented as a gradient flow, where state matrices asymptotically converge to a common one up to phase translation. We provide a sufficient framework leading to quantum synchronisation based on Riccati-type differential inequalities. In addition, uniform time-delayed interaction is considered for modelling realistic communication, and we demonstrate that quantum synchronisation is persistent when a small time delay is allowed. Finally, numerical simulation is performed to visualise qualitative behaviours and support theoretical results.

In this paper, we present a sufficient framework to exhibit the sample path-wise asymptotic flocking dynamics of the Cucker–Smale model with unit-speed constraint and the randomly switching network topology. We employ a matrix formulation of the given equation, which allows us to evaluate the diameter of velocities with respect to the adjacency matrix of the network. Unlike the previous result on the randomly switching Cucker–Smale model, the unit-speed constraint disallows the system to be considered as a nonautonomous linear ordinary differential equation on velocity vector, which forces us to get a weaker form of the flocking estimate than the result for the original Cucker–Smale model.

We consider self-propelled rigid bodies interacting through local body-attitude alignment modelled by stochastic differential equations. We derive a hydrodynamic model of this system at large spatio-temporal scales and particle numbers in any dimension $n \geq 3$. This goal was already achieved in dimension $n=3$ or in any dimension $n \geq 3$ for a different system involving jump processes. However, the present work corresponds to huge conceptual and technical gaps compared with earlier ones. The key difficulty is to determine an auxiliary but essential object, the generalised collision invariant. We achieve this aim by using the geometrical structure of the rotation group, namely its maximal torus, Cartan subalgebra and Weyl group as well as other concepts of representation theory and Weyl’s integration formula. The resulting hydrodynamic model appears as a hyperbolic system whose coefficients depend on the generalised collision invariant.

We propose a new fractional Laplacian for bounded domains, expressed as a conservation law and thus particularly suited to finite-volume schemes. Our approach permits the direct prescription of no-flux boundary conditions. We first show the well-posedness theory for the fractional heat equation. We also develop a numerical scheme, which correctly captures the action of the fractional Laplacian and its anomalous diffusion effect. We benchmark numerical solutions for the Lévy–Fokker–Planck equation against known analytical solutions. We conclude by numerically exploring properties of these equations with respect to their stationary states and long-time asymptotics.

Consensus-based optimisation (CBO) is a versatile multi-particle metaheuristic optimisation method suitable for performing non-convex and non-smooth global optimisations in high dimensions. It has proven effective in various applications while at the same time being amenable to a theoretical convergence analysis. In this paper, we explore a variant of CBO, which incorporates truncated noise in order to enhance the well-behavedness of the statistics of the law of the dynamics. By introducing this additional truncation in the noise term of the CBO dynamics, we achieve that, in contrast to the original version, higher moments of the law of the particle system can be effectively bounded. As a result, our proposed variant exhibits enhanced convergence performance, allowing in particular for wider flexibility in choosing the noise parameter of the method as we confirm experimentally. By analysing the time evolution of the Wasserstein-$2$ distance between the empirical measure of the interacting particle system and the global minimiser of the objective function, we rigorously prove convergence in expectation of the proposed CBO variant requiring only minimal assumptions on the objective function and on the initialisation. Numerical evidences demonstrate the benefit of truncating the noise in CBO.

Measuring inequalities in a multidimensional framework is a challenging problem, which is common to most field of science and engineering. Nevertheless, despite the enormous amount of researches illustrating the fields of application of inequality indices, and of the Gini index in particular, very few consider the case of a multidimensional variable. In this paper, we consider in some details a new inequality index, based on the Fourier transform, that can be fruitfully applied to measure the degree of inhomogeneity of multivariate probability distributions. This index exhibits a number of interesting properties that make it very promising in quantifying the degree of inequality in datasets of complex and multifaceted social phenomena.

We formulate haptotaxis models of cancer invasion wherein the infiltrating cancer cells can occupy a spectrum of states in phenotype space, ranging from ‘fully mesenchymal’ to ‘fully epithelial’. The more mesenchymal cells are those that display stronger haptotaxis responses and have greater capacity to modify the extracellular matrix (ECM) through enhanced secretion of matrix-degrading enzymes (MDEs). However, as a trade-off, they have lower proliferative capacity than the more epithelial cells. The framework is multiscale in that we start with an individual-based model that tracks the dynamics of single cells, which is based on a branching random walk over a lattice representing both physical and phenotype space. We formally derive the corresponding continuum model, which takes the form of a coupled system comprising a partial integro-differential equation for the local cell population density function, a partial differential equation for the MDE concentration and an infinite-dimensional ordinary differential equation for the ECM density. Despite the intricacy of the model, we show, through formal asymptotic techniques, that for certain parameter regimes it is possible to carry out a detailed travelling wave analysis and obtain invading fronts with spatial structuring of phenotypes. Precisely, the most mesenchymal cells dominate the leading edge of the invasion wave and the most epithelial (and most proliferative) dominate the rear, representing a bulk tumour population. As such, the model recapitulates similar observations into a front to back structuring of invasion waves into leader-type and follower-type cells, witnessed in an increasing number of experimental studies over recent years.

Opinion dynamics is an important and very active area of research that delves into the complex processes through which individuals form and modify their opinions within a social context. The ability to comprehend and unravel the mechanisms that drive opinion formation is of great significance for predicting a wide range of social phenomena such as political polarisation, the diffusion of misinformation, the formation of public consensus and the emergence of collective behaviours. In this paper, we aim to contribute to that field by introducing a novel mathematical model that specifically accounts for the influence of social media networks on opinion dynamics. With the rise of platforms such as Twitter, Facebook, and Instagram and many others, social networks have become significant arenas where opinions are shared, discussed and potentially altered. To this aim after an analytical construction of our new model and through incorporation of real-life data from Twitter, we calibrate the model parameters to accurately reflect the dynamics that unfold in social media, showing in particular the role played by the so-called influencers in driving individual opinions towards predetermined directions.

The connection between Residual Neural Networks (ResNets) and continuous-time control systems (known as NeurODEs) has led to a mathematical analysis of neural networks, which has provided interesting results of both theoretical and practical significance. However, by construction, NeurODEs have been limited to describing constant-width layers, making them unsuitable for modelling deep learning architectures with layers of variable width. In this paper, we propose a continuous-time Autoencoder, which we call AutoencODE, based on a modification of the controlled field that drives the dynamics. This adaptation enables the extension of the mean-field control framework originally devised for conventional NeurODEs. In this setting, we tackle the case of low Tikhonov regularisation, resulting in potentially non-convex cost landscapes. While the global results obtained for high Tikhonov regularisation may not hold globally, we show that many of them can be recovered in regions where the loss function is locally convex. Inspired by our theoretical findings, we develop a training method tailored to this specific type of Autoencoders with residual connections, and we validate our approach through numerical experiments conducted on various examples.

In this paper, we propose a novel approach that employs kinetic equations to describe the collective dynamics emerging from graph-mediated pairwise interactions in multi-agent systems. We formally show that for large graphs and specific classes of interactions a statistical description of the graph topology, given in terms of the degree distribution embedded in a Boltzmann-type kinetic equation, is sufficient to capture the collective trends of networked interacting systems. This proves the validity of a commonly accepted heuristic assumption in statistically structured graph models, namely that the so-called connectivity of the agents is the only relevant parameter to be retained in a statistical description of the graph topology. Then, we validate our results by testing them numerically against real social network data.