99 research outputs found
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 -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
Utilisation des équations Euler-PML en milieu hétérogène borné pour la résolution d'un problème inverse en géophysique.
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 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
-stability result (with ) is presented assuming the
aforementioned condition.Comment: 19 page
- …