We consider a cantilever beam which possesses a possibly non-uniform permanent magnetization, and whose shape is controlled by an applied magnetic field. We model the beam as a plane elastic curve and we suppose that the magnetic field acts upon the beam by means of a distributed couple that pulls the magnetization towards its direction. Given a list of target shapes, we look for a design of the magnetization

We present a posteriori error estimates for inconsistent and non-hierarchical Galerkin methods for linear parabolic problems, allowing them to be used in conjunction with very general mesh modification for the first time. We treat schemes which are non-hierarchical in the sense that the spatial Galerkin spaces between time-steps may be completely unrelated from one another. The practical interest of

In this paper, we establish a new hydrodynamically coupled phase-field model for three immiscible fluid components system. The model consists of the Navier–Stokes equations and three coupled nonlinear Allen–Cahn type equations, to which we add nonlocal type Lagrange multipliers to conserve the volume of each phase accurately. To solve the model, a linear and energy stable time-marching method is constructed

In Lipschitz two- and three-dimensional domains, we study the existence for the so-called Boussinesq model of thermally driven convection under singular forcing. By singular we mean that the heat source is allowed to belong to H−1(ϖ,Ω), where ϖ is a weight in the Muckenhoupt class A2 that is regular near the boundary. We propose a finite element scheme and, under the assumption that the domain is convex

We derive geometrically linearized theories for incompressible materials from nonlinear elasticity theory in the small displacement regime. Our nonlinear stored energy densities may vary on the same (small) length scale as the typical displacements. This allows for applications to multiwell energies as, e.g. encountered in martensitic phases of shape memory alloys and models for nematic elastomers

We develop a method to compute effectively the Young measures associated to sequences of numerical solutions of the compressible Euler system. Our approach is based on the concept of 𝒦-convergence adapted to sequences of parameterized measures. The convergence is strong in space and time (a.e. pointwise or in certain Lq spaces) whereas the measures converge narrowly or in the Wasserstein distance

Coronavirus Disease 2019 (COVID-19) is a zoonotic illness which has spread rapidly and widely since December, 2019, and is identified as a global pandemic by the World Health Organization. The pandemic to date has been characterized by ongoing cluster community transmission. Quarantine intervention to prevent and control the transmission are expected to have a substantial impact on delaying the growth

In this paper, we propose a novel space-dependent multiscale model for the spread of infectious diseases in a two-dimensional spatial context on realistic geographical scenarios. The model couples a system of kinetic transport equations describing a population of commuters moving on a large scale (extra-urban) with a system of diffusion equations characterizing the non-commuting population acting over

We present a uniform-in-time (and in particle numbers as well) error estimate for the random batch method (RBM) [S. Jin, L. Li and J.-G. Liu, Random batch methods (RBM) for interacting particle systems, J. Comput. Phys.400 (2020) 108877] to the Cucker–Smale (CS) model. The uniform-in-time error estimates of the RBM have been obtained for various interacting particle systems, when corresponding flow

In this paper, we present the hydrodynamic limit of a multiscale system describing the dynamics of two populations of agents with alignment interactions and the effect of an internal variable. It consists of a kinetic equation coupled with an Euler-type equation inspired by the thermomechanical Cucker–Smale (TCS) model. We propose a novel drag force for the fluid-particle interaction reminiscent of

This work introduces and analyzes new primal and dual-mixed finite element methods for deformable image registration, in which the regularizer has a nontrivial kernel, and constructed under minimal assumptions of the registration model: Lipschitz continuity of the similarity measure and ellipticity of the regularizer on the orthogonal complement of its kernel. The aforementioned singularity of the

Bounded-confidence models in social dynamics describe multi-agent systems, where each individual interacts only locally with others. Several models are written as systems of ordinary differential equations (ODEs) with discontinuous right-hand side: this is a direct consequence of restricting interactions to a bounded region with non-vanishing strength at the boundary. Various works in the literature

In this paper, we present a computational modeling approach for the dynamics of human crowds, where the spreading of an emotion (specifically fear) has an influence on the pedestrians’ behavior. Our approach is based on the methods of the kinetic theory of active particles. The model allows us to weight between two competing behaviors depending on fear level: the search for less congested areas and

A Markov chain individual-based model for virus diffusion is investigated. Both the virus growth within an individual and the complexity of the contagion within a population are taken into account. A careful work of parameter choice is performed. The model captures very well the statistical variability of quantities like the incubation period, the serial interval and the time series of infected people

In this work, we describe a hyperbolic model with cell–cell repulsion with a dynamics in the population of cells. More precisely, we consider a population of cells producing a field (which we call “pressure”) which induces a motion of the cells following the opposite of the gradient. The field indicates the local density of population and we assume that cells try to avoid crowded areas and prefer locally

Modeling and simulation of disease spreading in pedestrian crowds have recently become a topic of increasing relevance. In this paper, we consider the influence of the crowd motion in a complex dynamical environment on the course of infection of the pedestrians. To model the pedestrian dynamics, we consider a kinetic equation for multi-group pedestrian flow based on a social force model coupled with

In this work, we consider the two-species chemotaxis system with Lotka–Volterra competitive kinetics in a bounded domain with smooth boundary. We construct weak solutions and prove that they become smooth after some waiting time. In addition, the asymptotic behavior of the solutions is studied. Our results generalize some well-known results in the literature.

We formulate a general SEIR epidemic model in a heterogeneous population characterized by some trait in a discrete or continuous subset of a space ℝd. The incubation and recovery rates governing the evolution of each homogeneous subpopulation depend upon this trait, and no restriction is assumed on the contact matrix that defines the probability for an individual of a given trait to be infected by

We propose and study the framework of dissipative statistical solutions for the incompressible Euler equations. Statistical solutions are time-parameterized probability measures on the space of square-integrable functions, whose time-evolution is determined from the underlying Euler equations. We prove partial well-posedness results for dissipative statistical solutions and propose a Monte Carlo type

The goal of this paper is to investigate the optimality of the d-dimensional rock-salt structure, i.e. the cubic lattice V1/dℤd of volume V with an alternation of charges ±1 at lattice points, among periodic distributions of charges and lattice structures. We assume that the charges are interacting through two types of radially symmetric interaction potentials, according to their signs. We first restrict

In this paper, we quantify the asymptotic limit of collective behavior kinetic equations arising in mathematical biology modeled by Vlasov-type equations with nonlocal interaction forces and alignment. More precisely, we investigate the hydrodynamic limit of a kinetic Cucker–Smale flocking model with confinement, nonlocal interaction, and local alignment forces, linear damping and diffusion in velocity

The z-pinch is a classical steady state for the MHD model, where a confined plasma fluid is separated by vacuum, in the presence of a magnetic field which is generated by a prescribed current along the z direction. We develop a variational framework to study its stability in the absence of viscosity effect, and demonstrate for the first time that such a z-pinch is always unstable. Moreover, we discover

This paper is concerned with a spatially two-dimensional version of a chemotaxis system with logistic cell proliferation and death, for a singular tactic response of standard logarithmic type, and with interaction with a surrounding incompressible fluid through transport and buoyancy. Systems of this form are of significant relevance to the understanding of chemotaxis-fluid interaction, but the rigorous

This paper should be considered as an addendum to [A. Buffa and C. Giannelli, Adaptive isogeometric methods with hierarchical splines: Error estimator and convergence, Math. Models Methods Appl. Sci.26 (2016) 1–25] and [A. Buffa and C. Giannelli, Adaptive isogeometric methods with hierarchical splines: Optimality and convergence rates, Math. Models Methods Appl. Sci.27 (2017) 2781–2802] where Poincaré

The Boussinesq problem describing indentation of a rigid punch of arbitrary shape into a deformable solid body is studied within the context of a linear viscoelastic model. Due to the presence of a non-local integral constraint prescribing the total contact force, the unilateral indentation problem is formulated in the general form as a quasi-variational inequality with unknown indentation depth, and

We consider the incompressible flow of two immiscible fluids in the presence of a solid phase that undergoes changes in time due to precipitation and dissolution effects. Based on a seminal sharp interface model a phase-field approach is suggested that couples the Navier–Stokes equations and the solid’s ion concentration transport equation with the Cahn–Hilliard evolution for the phase fields. The

In this paper, we derive a new 2D brittle fracture model for thin shells via dimension reduction, where the admissible displacements are only normal to the shell surface. The main steps include to endow the shell with a small thickness, to express the three-dimensional energy in terms of the variational model of brittle fracture in linear elasticity, and to study the Γ-limit of the functional as the

We study the existence and properties of solutions and free boundaries of the one-phase Stefan problem with fractional diffusion posed in ℝN. In terms of the enthalpy h(x,t), the evolution equation reads ∂th+(−Δ)sΦ(h)=0, while the temperature is defined as u:=Φ(h):=max{h−L,0} for some constant L>0 called the latent heat, and (−Δ)s stands for the fractional Laplacian with exponent s∈(0,1). We prove

We propose a new discontinuous Galerkin (dG) method for a geometrically nonlinear Kirchhoff plate model for large isometric bending deformations. The minimization problem is nonconvex due to the isometry constraint. We present a practical discrete gradient flow that decreases the energy and computes discrete minimizers that satisfy a prescribed discrete isometry defect. We prove Γ-convergence of the

We deduce a model for glioma invasion that accounts for the dynamics of brain tissue being actively degraded by tumor cells via excessive acidity production, but also according to the local orientation of tissue fibers. Our approach has a multiscale character: we start with a microscopic description of single cell dynamics including biochemical and/or biophysical effects of the tumor microenvironment

We study one-dimensional Eulerian dynamics with nonlocal alignment interactions, featuring strong short-range alignment, and long-range misalignment. Compared with the well-studied Euler-alignment system, the presence of the misalignment brings different behaviors of the solutions, including the possible creation of vacuum at infinite time, which destabilizes the solutions. We show that with a strongly

Understanding how knowledge emerges and propagates within groups is crucial to explain the evolution of human populations. In this work, we introduce a mathematically oriented model that draws on individual-based approaches, inhomogeneous Markov chains and learning algorithms, such as those introduced in [F. Cucker and S. Smale, On the mathematical foundations of learning, Bull. Amer. Math. Soc. 39

We introduce a new stochastic differential model for global optimization of nonconvex functions on compact hypersurfaces. The model is inspired by the stochastic Kuramoto–Vicsek system and belongs to the class of Consensus-Based Optimization methods. In fact, particles move on the hypersurface driven by a drift towards an instantaneous consensus point, computed as a convex combination of the particle

We consider atomistic geometry relaxation in the context of linear tight binding models for point defects. A limiting model as Fermi-temperature is sent to zero is formulated, and an exponential rate of convergence for the nuclei configuration is established. We also formulate the thermodynamic limit model at zero Fermi-temperature, extending the results of [H. Chen, J. Lu and C. Ortner, Thermodynamic

An adaptive edge element method is designed to approximate a quasilinear H(curl)-elliptic problem in magnetism, based on a residual-type a posteriori error estimator and general marking strategies. The error estimator is shown to be both reliable and efficient, and its resulting sequence of adaptively generated solutions converges strongly to the exact solution of the original quasilinear system. Numerical

The z-pinch is a classical steady state for the MHD model, where a confined plasma fluid is separated by vacuum, in the presence of a magnetic field which is generated by a prescribed current along the z-direction. We develop a scaled variational framework to study its stability in the presence of viscosity effect, and demonstrate that any such z-pinch is always unstable. We also establish the existence

We study an implicit finite-volume scheme for nonlinear, non-local aggregation-diffusion equations which exhibit a gradient-flow structure, recently introduced in [R. Bailo, J. A. Carrillo and J. Hu, Fully discrete positivity-preserving and energy-dissipating schemes for aggregation-diffusion equations with a gradient flow structure, arXiv:1811.11502]. Crucially, this scheme keeps the dissipation property

This work is concerned with the iterative solution of systems of quasi-static multiple-network poroelasticity equations describing flow in elastic porous media that is permeated by single or multiple fluid networks. Here, the focus is on a three-field formulation of the problem in which the displacement field of the elastic matrix and, additionally, one velocity field and one pressure field for each

This paper focuses on the connections between four stochastic and deterministic models for the motion of straight screw dislocations. Starting from a description of screw dislocation motion as interacting random walks on a lattice, we prove explicit estimates of the distance between solutions of this model, an SDE system for the dislocation positions, and two deterministic mean-field models describing

In this work, we consider the quasilinear attraction–repulsion chemotaxis model with nonlinear signal production and logistic-type source. We present the global existence of classical solutions under appropriate regularity assumptions on the initial data. In addition, the asymptotic behavior of the solutions is studied, and our results generalize and improve some well-known results in the literature

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

Based on a bootstrap instability method, we prove the existence of unstable strong solutions in the sense of L1-norm to an abstract Rayleigh–Taylor (RT) problem arising from stratified viscous fluids in Lagrangian coordinates. In the proof we develop a method to modify the initial data of the linearized abstract RT problem by exploiting the existence theory of a unique solution to the stratified (steady)

The paper introduces a model of collective behavior where agents receive information only from sufficiently dense crowds in their immediate vicinity. The system is an asymmetric, density-induced version of the Cucker–Smale model with short-range interactions. We prove the basic mathematical properties of the system and concentrate on the presentation of interesting behaviors of the solutions. The results

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

A system of partial differential equations for a diffusion interface model is considered for the stationary motion of two macroscopically immiscible, viscous Newtonian fluids in a three-dimensional bounded domain. The governing equations consist of the stationary Navier–Stokes equations for compressible fluids and a stationary Cahn–Hilliard type equation for the mass concentration difference. Approximate

We study the convergence behavior of Dual-Primal Isogeometric Tearing and Interconnecting (IETI-DP) methods for solving large-scale algebraic systems arising from multi-patch Isogeometric Analysis. We focus on the Poisson problem on two-dimensional computational domains. We provide a convergence analysis that covers several choices of the primal degrees of freedom: the vertex values, the edge averages

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

Variational multiscale methods and their precursors, stabilized methods, which are sometimes supplemented with discontinuity-capturing (DC) methods, have been playing their core-method role in flow computations increasingly with isogeometric discretization. The stabilization and DC parameters embedded in most of these methods play a significant role. The parameters almost always involve some local-length-scale

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

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

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 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