We construct a numerical scheme based on the scalar auxiliary variable (SAV) approach in time and the MAC discretization in space for the Cahn–Hilliard–Navier–Stokes phase- field model, prove its energy stability, and carry out error analysis for the corresponding Cahn–Hilliard–Stokes model only. The scheme is linear, second-order, unconditionally energy stable and can be implemented very efficiently

This editorial paper is devoted to present the papers published in a special issue focused on modeling, qualitative analysis and simulation of the collective dynamics of living, self-propelled particles. A critical analysis of the overall contents of the issue is proposed, thus leading to a forward look to research perspectives.

In this paper, we develop a continuum model for the movement of agents on a lattice, taking into account location desirability, local and far-range migration, and localized entry and exit rates. Specifically, our motivation is to qualitatively describe the homeless population in Los Angeles. The model takes the form of a fully nonlinear, nonlocal, non-degenerate parabolic partial differential equation

The goal of this work is to study an infectious disease spreading in a medium size population occupying a confined environment. For this purpose, we consider a kinetic theory approach to model crowd dynamics in bounded domains and couple it to a kinetic equation to model contagion. The interactions of a person with other pedestrians and the environment are modeled by using tools of game theory. The

Collective behaviours in animal communities are the result of inter-individual communication. However, communication signals are not fixed; they evolve to ensure more effective interactions between the emitter and receiver of these signals. In this study, we use a mathematical approach and investigate the effect of changes in communication signals (at both receiver and emitter levels) on the aggregation

Starting from a particle model describing self-propelled particles interacting through nematic alignment, we derive a macroscopic model for the particle density and mean direction of motion. We first propose a mean-field kinetic model of the particle dynamics. After diffusive rescaling of the kinetic equation, we formally show that the distribution function converges to an equilibrium distribution

We introduce a class of one-dimensional linear kinetic equations of Boltzmann and Fokker–Planck type, describing the dynamics of individuals of a multi-agent society questing for high status in the social hierarchy. At the Boltzmann level, the microscopic variation of the status of agents around a universal desired target, is built up introducing as main criterion for the change of status a suitable

Global optimization of a non-convex objective function often appears in large-scale machine learning and artificial intelligence applications. Recently, consensus-based optimization (CBO) methods have been introduced as one of the gradient-free optimization methods. In this paper, we provide a convergence analysis for the first-order CBO method in [J. A. Carrillo, S. Jin, L. Li and Y. Zhu, A consensus-based

We study a critical exponent of the flocking behavior to the one-dimensional 1D Cucker–Smale (C–S) model with a regular inverse power law communication on a general network with a spanning tree. For this, we propose a new nonlinear functional which can control the velocity diameter and decays exponentially fast as time goes on. As an application of the time-evolution of the nonlinear functional, we

We consider the electric field integral equation (EFIE) modeling the scattering of time-harmonic electromagnetic waves at a perfectly conducting screen. When discretizing the EFIE by means of low-order Galerkin boundary methods (BEM), one obtains linear systems that are ill-conditioned on fine meshes and for low wave numbers k. This makes iterative solvers perform poorly and entails the use of preconditioning

Urban crime such as residential burglary is a social problem in every major urban area. As such, many mathematical models have been proposed to study the collective behavior of these crimes. In [V. B. Pasour, G. E. Tita, P. J. Brantingham, A. L. Bertozzi, M. B. Short, M. R. D’Orsogna and L. B. Chayes, A statistical model of crime behavior, Math. Methods Appl. Sci107 (2008) 1249–1267; M. B. Short, A

We consider linear reaction systems with slow and fast reactions, which can be interpreted as master equations or Kolmogorov forward equations for Markov processes on a finite state space. We investigate their limit behavior if the fast reaction rates tend to infinity, which leads to a coarse-grained model where the fast reactions create microscopically equilibrated clusters, while the exchange mass

In this work, merging ideas from compatible discretisations and polyhedral methods, we construct novel fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra. The spaces and operators that appear in these sequences are directly amenable to computer implementation. Besides proving the exactness, we show that the usual three-dimensional sequence of trimmed Finite Element

We consider a class of macroscopic models for the spatio-temporal evolution of urban crime, as originally going back to Ref. 29 [M. B. Short, M. R. D’Orsogna, V. B. Pasour, G. E. Tita, P. J. Brantingham, A. L. Bertozzi and L. B. Chayes, A statistical model of criminal behavior, Math. Models Methods Appl. Sci.18 (2008) 1249–1267]. The focus here is on the question of how far a certain porous medium

We present a new hydrodynamic model for synchronization phenomena which is a type of pressureless Euler system with nonlocal interaction forces. This system can be formally derived from the Kuramoto model with inertia, which is a classical model of interacting phase oscillators widely used to investigate synchronization phenomena, through a kinetic description under the mono-kinetic closure assumption

In this paper, a phase-field approach for structural topology optimization for a 3D-printing process which includes stress constraints and potentially multiple materials or multiscales is analyzed. First-order necessary optimality conditions are rigorously derived and a numerical algorithm which implements the method is presented. A sensitivity study with respect to some parameters is conducted for

We study the existence of global weak solutions in a three-dimensional time-dependent bounded domain for the incompressible Vlasov–Navier–Stokes system which is coupled with two convection–diffusion equations describing the air temperature and its water vapor mass fraction. This newly introduced model describes respiratory aerosols in the human aiways when one takes into account the hygroscopic effects

In this paper, we consider a compressible dissipative Baer–Nunziato-type system for a mixture of two compressible heat conducting gases. We prove that the set of weak solutions is stable, meaning that any sequence of weak solutions contains a (weakly) convergent subsequence whose limit is again a weak solution to the original system. Such type of results is usually considered as the most essential

In this paper, we tackle the problem of constructing conforming Virtual Element spaces on polygons with curved edges. Unlike previous VEM approaches for curvilinear elements, the present construction ensures that the local VEM spaces contain all the polynomials of a given degree, thus providing the full satisfaction of the patch test. Moreover, unlike standard isoparametric FEM, this approach allows

This paper is devoted to the multidisciplinary modelling of a pandemic initiated by an aggressive virus, specifically the so-called SARS–CoV–2 Severe Acute Respiratory Syndrome, corona virus n.2. The study is developed within a multiscale framework accounting for the interaction of different spatial scales, from the small scale of the virus itself and cells, to the large scale of individuals and further

We concentrate on kinetic models for swarming with individuals interacting through self-propelling and friction forces, alignment and noise. We assume that the velocity of each individual relaxes to the mean velocity. In our present case, the equilibria depend on the density and the orientation of the mean velocity, whereas the mean speed is not anymore a free parameter and a phase transition occurs

In this paper, we study the macroscopic behavior of the inertial spin (IS) model. This model has been recently proposed to describe the collective dynamics of flocks of birds, and its main feature is the presence of an auxiliary dynamical variable, a sort of internal spin, which conveys the interaction among the birds with the effect of better describing the turning of flocks. After discussing the

Chemotherapy is a common treatment for advanced prostate cancer. The standard approach relies on cytotoxic drugs, which aim at inhibiting proliferation and promoting cell death. Advanced prostatic tumors are known to rely on angiogenesis, i.e. the growth of local microvasculature via chemical signaling produced by the tumor. Thus, several clinical studies have been investigating antiangiogenic therapy

This paper concerns with the global existence and boundedness of classical solution of the higher-dimensional forager–exploiter model with homogeneous Neumann boundary condition and nonnegative initial data. For cases where there are no forager and exploiter growth sources, it will be shown that if either the initial data and the production rate of nutrient are small or the taxis effects are small

The diffusion of cells in a viscous incompressible fluid (e.g. water) may be viewed like movement in a porous medium and there is a bidirectorial oxygen exchange between water and their surrounding air in thin fluid layers near the air–water contact surface. This leads to the following chemotaxis-Navier–Stokes system with nonlinear diffusion: nt+u⋅∇n=Δnm−∇⋅(n∇c),x∈Ω, t>0,ct+u⋅∇c=Δc−nc,x∈Ω, t>0,ut+(u⋅∇)u=Δu+∇P+n∇ϕ

We study the effect of small random advection in two models in turbulent combustion. Assuming that the velocity field decorrelates sufficiently fast, we (i) identify the order of the fluctuations of the front with respect to the size of the advection; and (ii) characterize them by the solution of a Hamilton–Jacobi equation forced by white noise. In the simplest case, the result yields, for both models

Starting with this paper, we intend to develop a program aiming at construction of boundary conditions (BCs) for hyperbolic relaxation systems. Physically, such BCs are not always available. The construction is based on the assumption that the relaxation systems and well-posed BCs for the corresponding equilibrium systems are given. This paper focuses on the linearized Suliciu model. We obtain strictly

This paper constitutes the first attempt to bridge the evolutionary theory in economics and the theory of active particles in mathematics. It seeks to present a kinetic model for an evolutionary formalization of economic dynamics. The new derived mathematical representation intends to formalize the processes of learning and selection as the two fundamental drivers of evolutionary environments [G. Dosi

This paper is concerned with the parabolic system ut=Δu−∇⋅uw∇w,(x,t)∈Ω×(0,T),vt=Δv−∇⋅vu∇u,(x,t)∈Ω×(0,T),wt=Δw−(u+v)w−λw+g(x,t),(x,t)∈Ω×(0,T),(0.1) in a bounded ball Ω=BR(0)⊂ℝn (n≥2) with R>0. Where λ≥0 and 0≤g∈C1(Ω¯×(0,∞)). It is shown that for arbitrarily radially symmetric initial data (u0,v0,w0), which are nonnegative and suitably regular, the corresponding Neumann initial-boundary problem admits

The main purpose of this work is to explore a general cell–fluid model which is based on a mixture theory formulation that accounts for the interplay between oxytactically (chemotaxis toward gradient in oxygen) moving bacteria cells in water and the buoyance forces caused by the difference in density between cells and fluid. The model involves two mass balance and two general momentum balance equations

The run and tumble process is well established in order to describe the movement of bacteria in response to a chemical stimulus. However, the relation between the tumbling rate and the internal state of bacteria is poorly understood. This study aims at deriving macroscopic models as limits of the mesoscopic kinetic equation in different regimes. In particular, we are interested in the roles of the

We consider the coupled chemotaxis–Navier–Stokes system with logistic source term nt+u⋅∇n=Δn−∇⋅(n∇c)+rn−μnα,ct+u⋅∇c=Δc−nc,ut+(u⋅∇)u=Δu+∇P+n∇Φ,∇⋅u=0 in a bounded, smooth domain Ω⊂ℝ3, where Φ∈W2,∞(Ω) and where r≥0, μ>0 and 1<α<2 are given parameters. Although the degradation here is weaker than the usual quadratic case, it is proved that for any sufficiently regular initial data, the initial-value problem

In this paper, we consider a cascaded taxis model for two proliferating and degrading species which thrive on the same nutrient but orient their movement according to different schemes. In particular, we assume the first group, the foragers, to orient their movement directly along an increasing gradient of the food density, while the second group, the exploiters, instead track higher densities of the

Nonlinear diffusion problems featuring stochastic effects may be described by stochastic partial differential equations of the form dα(u)−div(β1(∇u))dt+β0(u)dt∋f(u)dt+G(u)dW. We present an existence theory for such equations under general monotonicity assumptions on the nonlinearities. In particular, α, β0, and β1 are allowed to be multivalued, as required by the modelization of solid–liquid phase

We consider a spatially extended mean-field model of a FitzHugh–Nagumo neural network, with a rescaled interaction kernel. Our main purpose is to prove that its asymptotic limit in the regime of strong local interactions converges toward a system of reaction–diffusion equations taking account for the average quantities of the network. Our approach is based on a modulated energy argument, to compare

The heterogeneity and plasticity of macrophages have become a topic of great interest, due to their role in various diseases ranging from cancer to bacterial infections. While initial experimental studies assumed an extreme polarisation situation, with the (anti-tumour) M1 and (pro-tumour) M2 macrophages representing the two extreme cell phenotypes, more recent studies showed a continuum of macrophages

We consider a linearised model of incompressible inviscid flow. Using a regularisation based on the Hodge Laplacian we prove existence and uniqueness of weak solutions for smooth domains. The model problem is then discretised using H(div)-conforming finite element methods, for which we prove error estimates for the velocity approximation in the L2-norm of order O(hk+12). We also prove error estimates

This paper analyzes spectral properties of linear Schrödinger operators under oscillatory high-amplitude potentials on bounded domains. Depending on the degree of disorder, we prove the existence of spectral gaps among the lowermost eigenvalues and the emergence of exponentially localized states. We quantify the rate of decay in terms of geometric parameters that characterize the potential. The proofs

This paper presents a novel approach to quantitatively analyzing pedestrian congestion in evacuation management based on the Hughes and social force models. An accurate analysis of crowds plays an important role in illustrating their dynamics. However, the majority of the existing approaches to analyzing pedestrian congestion are qualitative. Few methods focus on the quantification of the interactions

We investigate the effect of linear diffusion and interactions with the domain boundary on swarm equilibria by analyzing critical points of the associated energy functional. Through this process we uncover two properties of energy minimization that depend explicitly on the spatial domain: (i) unboundedness from below of the energy due to an imbalance between diffusive and aggregative forces depends

This paper is concerned with the critical threshold phenomenon for one-dimensional damped, pressureless Euler–Poisson equations with electric force induced by a constant background, originally studied in [S. Engelberg and H. Liu and E. Tadmor, Indiana Univ. Math. J.50 (2001) 109–157]. A simple transformation is used to linearize the characteristic system of equations, which allows us to study the geometrical

We present a stochastic Justh–Krishnaprasad flocking model describing interactions among individuals in a planar domain with their positions and heading angles. The deterministic counterpart of the proposed model describes the formation of nematic alignment in an ensemble of planar particles moving with a unit speed. When the noise is turned off, we show that the nematic alignment state, in which all

This paper has two objectives. On one side, we develop and test numerically divergence-free Virtual Elements in three dimensions, for variable “polynomial” order. These are the natural extension of the two-dimensional divergence-free VEM elements, with some modification that allows for a better computational efficiency. We test the element’s performance both for the Stokes and (diffusion dominated)

We model and analyze a guiding problem, where the drivers try to steer the evaders’ positions toward a target region while the evaders always try to escape from drivers. This problem is motivated by the guidance-by-repulsion model [R. Escobedo, A. Ibañez and E. Zuazua, Optimal strategies for driving a mobile agent in a “guidance by repulsion” model, Commun. Nonlinear Sci. Numer. Simul.39 (2016) 58–72]

In this paper, we discuss higher Sobolev regularity of convex integration solutions for the geometrically nonlinear two-well problem. More precisely, we construct solutions to the differential inclusion ∇u∈K subject to suitable affine boundary conditions for u with K:=SO(2)1δ01∪SO(2)1−δ01 such that the associated deformation gradients ∇u enjoy higher Sobolev regularity. This provides the first result

Statistical solutions are time-parameterized probability measures on spaces of integrable functions, which have been proposed recently as a framework for global solutions and uncertainty quantification for multi-dimensional hyperbolic system of conservation laws. By combining high-resolution finite volume methods with a Monte Carlo sampling procedure, we present a numerical algorithm to approximate

A two-dimensional mathematical model for cells migrating without adhesion capabilities is presented and analyzed. Cells are represented by their cortex, which is modeled as an elastic curve, subject to an internal pressure force. Net polymerization or depolymerization in the cortex is modeled via local addition or removal of material, driving a cortical flow. The model takes the form of a fully nonlinear

We consider a finite-difference semi-discrete scheme for the approximation of internal controls of a one-dimensional evolution problem of hyperbolic type involving the spectral fractional Laplacian. The continuous problem is controllable in arbitrary small time. However, the high frequency numerical spurious oscillations lead to a loss of the uniform (with respect to the mesh size) controllability

This paper first provides an introduction to the mathematical approach to the modeling, qualitative analysis, and simulation of large systems of living entities, specifically self-propelled particles. Subsequently, a presentation of the papers published in this special issue follows. Finally, a critical analysis of the overall contents of the issue is proposed, thus leading to define some challenging

We introduce a class of new one-dimensional linear Fokker–Planck-type equations describing the dynamics of the distribution of wealth in a multi-agent society. The equations are obtained, via a standard limiting procedure, by introducing an economically relevant variant to the kinetic model introduced in 2005 by Cordier, Pareschi and Toscani according to previous studies by Bouchaud and Mézard. The

In this work, we propose a method of opinion pooling based on pairwise interactions. We assume that each agent has a probability measure on the possible outcomes of some situation, and they try to find a single measure aggregating their estimates. This is a classical problem in Decision Theory, where expert opinions contain some degree of uncertainty, and a Decision Taker needs to pool these estimates

In the frame of isogeometric analysis, we consider a Galerkin boundary element discretization of the hyper-singular integral equation associated with the 2D Laplacian. We propose and analyze an adaptive algorithm which locally refines the boundary partition and, moreover, steers the smoothness of the NURBS ansatz functions across elements. In particular and unlike prior work, the algorithm can increase

We consider a hyperbolic–parabolic system arising from a chemotaxis model in tumor angiogenesis, which is described by a Keller–Segel equation with singular sensitivity. It is known to allow viscous shocks (so-called traveling waves). We introduce a relative entropy of the system, which can capture how close a solution at a given time is to a given shock wave in almost L2-sense. When the shock strength

We consider a Landau–de Gennes model for a suspension of small colloidal inclusions in a nematic host. We impose suitable anchoring conditions at the boundary of the inclusions, and we work in the dilute regime — i.e. the size of the inclusions is much smaller than the typical separation distance between them, so that the total volume occupied by the inclusions is small. By studying the homogenised

We consider an Oldroyd-B model which is derived in Ref. 4 [J. W. Barrett, Y. Lu and E. Süli, Existence of large-data finite-energy global weak solutions to a compressible Oldroyd-B model, Commun. Math. Sci.15 (2017) 1265–1323] via micro–macro-analysis of the compressible Navier–Stokes–Fokker–Planck system. The global well posedness of strong solutions as well as the associated time-decay estimates

We propose and analyze several multilevel algorithms for the fast simulation of possibly nonstationary Gaussian random fields (GRFs) indexed, for example, by the closure of a bounded domain 𝒟⊂ℝn or, more generally, by a compact metric space 𝒳 such as a compact n-manifold ℳ. A colored GRF 𝒵, admissible for our algorithms, solves the stochastic fractional-order equation 𝒜β𝒵=𝒲 for some β>n/4, where

This paper is devoted to the study of the nonlinear stability of the composite wave consisting of two rarefaction waves and a viscous contact wave for the Cauchy problem to a one-dimensional compressible non-isentropic Navier–Stokes/Allen–Cahn system which is a combination of the classical Navier–Stokes system with an Allen–Cahn phase field description. We first construct the composite wave through

We consider the compressible Navier–Stokes–Cahn–Hilliard system describing the behavior of a binary mixture of compressible, viscous and macroscopically immiscible fluids. The equations are endowed with dynamic boundary conditions which allows taking into account the interaction between the fluid components and the rigid walls of the physical domain. We establish the existence of global-in-time weak

After the pioneering work of Garrett and Munk, the statistics of oceanic internal gravity waves has become a central subject of research in oceanography. The time evolution of the spectral energy of internal waves in the ocean can be described by a near-resonance wave turbulence equation, of quantum Boltzmann type. In this work, we provide the first rigorous mathematical study for the equation by showing

We obtain regularity results in weighted Sobolev spaces for the solution of the obstacle problem for the integral fractional Laplacian (−Δ)s in a Lipschitz bounded domain Ω⊂ℝn satisfying the exterior ball condition. The weight is a power of the distance to the boundary ∂Ω of Ω that accounts for the singular boundary behavior of the solution for any 0