Efficient numerical calculation of drift and diffusion coefficients in the diffusion approximation of kinetic equations

In this paper we study the diffusion approximation of a swarming model given by a system of interacting Langevin equations with nonlinear friction. The diffusion approximation requires the calculation of the drift and diffusion coefficients that are given as averages of solutions to appropriate Poisson equations. We present a new numerical method for computing these coefficients that is based on the calculation of the eigenvalues and eigenfunctions of a Schr\"odinger operator. These theoretical results are supported by numerical simulations showcasing the efficiency of the method

Stochastic and deterministic models for age-structured populations with genetically variable traits

Understanding how stochastic and non-linear deterministic processes interact is a major challenge in population dynamics theory. After a short review, we introduce a stochastic individual-centered particle model to describe the evolution in continuous time of a population with (continuous) age and trait structures. The individuals reproduce asexually, age, interact and die. The 'trait' is an individual heritable property (d-dimensional vector) that may influence birth and death rates and interactions between individuals, and vary by mutation. In a large population limit, the random process converges to the solution of a Gurtin-McCamy type PDE. We show that the random model has a long time behavior that differs from its deterministic limit. However, the results on the limiting PDE and large deviation techniques \textit{\`a la} Freidlin-Wentzell provide estimates of the extinction time and a better understanding of the long time behavior of the stochastic process. This has applications to the theory of adaptive dynamics used in evolutionary biology. We present simulations for two biological problems involving life-history trait evolution when body size is plastic and individual growth is taken into account.Comment: This work is a proceeding of the CANUM 2008 conferenc

Two-dimensional pseudo-gravity model: particles motion in a non-potential singular force field

We analyze a simple macroscopic model describing the evolution of a cloud of particles confined in a magneto-optical trap. The behavior of the particles is mainly driven by self-consistent attractive forces. In contrast to the standard model of gravitational forces, the force field does not result from a potential; moreover, the nonlinear coupling is more singular than the coupling based on the Poisson equation. We establish the existence of solutions under a suitable smallness condition on the total mass or, equivalently, for a sufficiently large diffusion coefficient. When a symmetry assumption is fulfilled, the solutions satisfy strengthened estimates (exponential moments). We also investigate the convergence of the $N$-particles description towards the PDE system in the mean field regime

A numerical approach for modelling thin cracked plates with XFEM

The modelization of bending plates with through the thickness cracks is investigated. We consider the Kirchhoff-Love plate model which is valid for very thin plates. We apply the eXtended Finite Element Method (XFEM) strategy: enrichment of the finite element space with the asymptotic bending and with the discontinuity across the crack. We present two variants and their numerical validations and also a numerical computation of the stress intensity factors

Efficient numerical calculation of drift and diffusion coefficients in the diffusion approximation of kinetic equations

In this paper we study the diffusion approximation of a swarming model given by a system of interacting Langevin equations with nonlinear friction. The diffusion approximation requires the calculation of the drift and diffusion coefficients that are given as averages of solutions to appropriate Poisson equations. We present a new numerical method for computing these coefficients that is based on the calculation of the eigenvalues and eigenfunctions of a Schr¨odinger operator. These theoretical results are supported by numerical simulations showcasing the efficiency of the method

Quantitative lower bounds for the full Boltzmann equation, Part I: Periodic boundary conditions

We prove the appearance of an explicit lower bound on the solution to the full Boltzmann equation in the torus for a broad family of collision kernels including in particular long-range interaction models, under the assumption of some uniform bounds on some hydrodynamic quantities. This lower bound is independent of time and space. When the collision kernel satisfies Grad's cutoff assumption, the lower bound is a global Maxwellian and its asymptotic behavior in velocity is optimal, whereas for non-cutoff collision kernels the lower bound we obtain decreases exponentially but faster than the Maxwellian. Our results cover solutions constructed in a spatially homogeneous setting, as well as small-time or close-to-equilibrium solutions to the full Boltzmann equation in the torus. The constants are explicit and depend on the a priori bounds on the solution.Comment: 37 page

Global existence for a system of non-linear and non-local transport equations describing the dynamics of dislocation densities

In this paper, we study the global in time existence problem for the Groma-Balogh model describing the dynamics of dislocation densities. This model is a two-dimensional model where the dislocation densities satisfy a system of transport equations such that the velocity vector field is the shear stress in the material, solving the equations of elasticity. This shear stress can be expressed as some Riesz transform of the dislocation densities. The main tool in the proof of this result is the existence of an entropy for this syste

Distributional and classical solutions to the Cauchy Boltzmann problem for soft potentials with integrable angular cross section

This paper focuses on the study of existence and uniqueness of distributional and classical solutions to the Cauchy Boltzmann problem for the soft potential case assuming $S^{n-1}$ integrability of the angular part of the collision kernel (Grad cut-off assumption). For this purpose we revisit the Kaniel--Shinbrot iteration technique to present an elementary proof of existence and uniqueness results that includes large data near a local Maxwellian regime with possibly infinite initial mass. We study the propagation of regularity using a recent estimate for the positive collision operator given in [3], by E. Carneiro and the authors, that permits to study such propagation without additional conditions on the collision kernel. Finally, an $L^{p}$-stability result (with $1\leq p\leq\infty$) is presented assuming the aforementioned condition.Comment: 19 page