Source author record

Mohammed Lemou

Mohammed Lemou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2020arXiv

Uniformly accurate numerical schemes for a class of dissipative systems

We consider a class of relaxation problems mixing slow and fast variations which can describe population dynamics models or hyperbolic systems, with varying stiffness (from non-stiff to strongly dissipative), and develop a multi-scale method by decomposing this problem into a micro-macro system where the original stiffness is broken. We show that this new problem can therefore be simulated with a uniform order of accuracy using standard explicit numerical schemes. In other words, it is possible to solve the micro-macro problem with a cost independent of the stiffness (a.k.a. uniform cost), such that the error is also uniform. This method is successfully applied to two hyperbolic systems with and without non-linearities, and is shown to circumvent the phenomenon of order reduction.

preprint2019arXiv

Nonlinear instability of inhomogeneous steady states solutions to the HMF Model

In this work we prove the nonlinear instability of inhomogeneous steady states solutions to the Hamiltonian Mean Field (HMF) model. We first study the linear instability of this model under a simple criterion by adapting the techniques developed in [19]. In a second part, we extend to the inhomogeneous case some techniques developed in [14, 17, 18] and prove a nonlinear instability result under the same criterion.

preprint2018arXiv

Uniformly accurate methods for Vlasov equations with non-homogeneous strong magnetic field

In this paper, we consider the numerical solution of highly-oscillatory Vlasov and Vlasov-Poisson equations with non-homogeneous magnetic field. Designed in the spirit of recent uniformly accurate methods, our schemes remain insensitive to the stiffness of the problem, in terms of both accuracy and computational cost. The specific difficulty (and the resulting novelty of our approach) stems from the presence of a non-periodic oscillation, which necessitates a careful ad-hoc reformulation of the equations. Our results are illustrated numerically on several examples.

preprint2016arXiv

Averaging of highly-oscillatory transport equations

In this paper, we develop a new strategy aimed at obtaining high-order asymptotic models for transport equations with highly-oscillatory solutions. The technique relies upon recent developments averaging theory for ordinary differential equations, in particular normal form expansions in the vanishing parameter. Noteworthy, the result we state here also allows for the complete recovery of the exact solution from the asymptotic model. This is done by solving a companion transport equation that stems naturally from the change of variables underlying high-order averaging. Eventually, we apply our technique to the Vlasov equation with external electric and magnetic fields. Both constant and non-constant magnetic fields are envisaged, and asymptotic models already documented in the literature and re-derived using our methodology. In addition, it is shown how to obtain new high-order asymptotic models.

preprint2016arXiv

Nonlinear Geometric Optics method based multi-scale numerical schemes for highly-oscillatory transport equations

We introduce a new numerical strategy to solve a class of oscillatory transport PDE models which is able to captureaccurately the solutions without numerically resolving the high frequency oscillations {\em in both space and time}.Such PDE models arise in semiclassical modeling of quantum dynamics with band-crossings, and otherhighly oscillatory waves. Our first main idea is to use the nonlinear geometric optics ansatz, which builds theoscillatory phase into an independent variable. We then choose suitable initial data, based on the Chapman-Enskog expansion, for the new model. For a scalar model, we prove that so constructed model will have certain smoothness, and consequently, for a first order approximation scheme we prove uniform error estimates independent of the (possibly small) wave length. The method is extended to systems arising from a semiclassical model for surface hopping, a non-adiabatic quantum dynamic phenomenon. Numerous numerical examples demonstrate that the method has the desired properties.

preprint2016arXiv

Uniformly accurate numerical schemes for the nonlinear Dirac equation in the nonrelativistic limit regime

We apply the two-scale formulation approach to propose uniformly accurate (UA) schemes for solving the nonlinear Dirac equation in the nonrelativistic limit regime. The nonlinear Dirac equation involves two small scales $\varepsilon$ and $\varepsilon^2$ with $\varepsilon\to0$ in the nonrelativistic limit regime. The small parameter causes high oscillations in time which brings severe numerical burden for classical numerical methods. We transform our original problem as a two-scale formulation and present a general strategy to tackle a class of highly oscillatory problems involving the two small scales $\varepsilon$ and $\varepsilon^2$. Suitable initial data for the two-scale formulation is derived to bound the time derivatives of the augmented solution. Numerical schemes with uniform (with respect to $\varepsilon\in (0,1]$) spectral accuracy in space and uniform first order or second order accuracy in time are proposed. Numerical experiments are done to confirm the UA property.

preprint2015arXiv

Asymptotic Preserving numerical schemes for multiscale parabolic problems

We consider a class of multiscale parabolic problems with diffusion coefficients oscillating in space at a possibly small scale $\varepsilon$. Numerical homogenization methods are popular for such problems, because they capture efficiently the asymptotic behaviour as $\varepsilon \rightarrow 0$, without using a dramatically fine spatial discretization at the scale of the fast oscillations. However, known such homogenization schemes are in general not accurate for both the highly oscillatory regime $\varepsilon \rightarrow 0$ and the non oscillatory regime $\varepsilon \sim 1$. In this paper, we introduce an Asymptotic Preserving method based on an exact micro-macro decomposition of the solution which remains consistent for both regimes.

preprint2015arXiv

Extended Rearrangement inequalities and applications to some quantitative stability results

In this paper, we prove a new functional inequality of Hardy-Littlewood type for generalized rearrangements of functions. We then show how this inequality provides {\em quantitative} stability results of steady states to evolution systems that essentially preserve the rearrangements and some suitable energy functional, under minimal regularity assumptions on the perturbations. In particular, this inequality yields a {\em quantitative} stability result of a large class of steady state solutions to the Vlasov-Poisson systems, and more precisely we derive a quantitative control of the $L^1$ norm of the perturbation by the relative Hamiltonian (the energy functional) and rearrangements. A general non linear stability result has been obtained in \cite{LMR} in the gravitational context, however the proof relied in a crucial way on compactness arguments which by construction provides no quantitative control of the perturbation. Our functional inequality is also applied to the context of 2D-Euler system and also provides quantitative stability results of a large class of steady-states to this system in a natural energy space.

preprint2015arXiv

Multiscale numerical schemes for kinetic equations in the anomalous diffusion limit

We construct numerical schemes to solve kinetic equations with anomalous diffusion scaling. When the equilibrium is heavy-tailed or when the collision frequency degenerates for small velocities, an appropriate scaling should be made and the limit model is the so-called anomalous or fractional diffusion model. Our first scheme is based on a suitable micro-macro decomposition of the distribution function whereas our second scheme relies on a Duhamel formulation of the kinetic equation. Both are \emph{Asymptotic Preserving} (AP): they are consistent with the kinetic equation for all fixed value of the scaling parameter $\varepsilon >0$ and degenerate into a consistent scheme solving the asymptotic model when $\varepsilon$ tends to $0$. The second scheme enjoys the stronger property of being uniformly accurate (UA) with respect to $\varepsilon$. The usual AP schemes known for the classical diffusion limit cannot be directly applied to the context of anomalous diffusion scaling, since they are not able to capture the important effects of large and small velocities. We present numerical tests to highlight the efficiency of our schemes.

preprint2015arXiv

Nonlinear stability criteria for the HMF Model

We study the nonlinear stability of a large class of inhomogeneous steady state solutions to the Hamiltonian Mean Field (HMF) model. Under a simple criterion, we prove the nonlinear stability of steady states which are decreasing functions of the microscopic energy. To achieve this task, we extend to this context the strategy based on generalized rearrangement techniques which was developed recently for the gravitational Vlasov-Poisson equation. Explicit stability inequalities are established and our analysis is able to treat non compactly supported steady states to HMF, which are physically relevant in this context but induces additional difficulties, compared to the Vlasov-Poisson system.

preprint2015arXiv

Numerical schemes for kinetic equations in the diffusion and anomalous diffusion limits. Part I: the case of heavy-tailed equilibrium

In this work, we propose some numerical schemes for linear kinetic equations in the diffusion and anomalous diffusion limit. When the equilibrium distribution function is a Maxwellian distribution, it is well known that for an appropriate time scale, the small mean free path limit gives rise to a diffusion type equation. However, when a heavy-tailed distribution is considered, another time scale is required and the small mean free path limit leads to a fractional anomalous diffusion equation. Our aim is to develop numerical schemes for the original kinetic model which works for the different regimes, without being restricted by stability conditions of standard explicit time integrators. First, we propose some numerical schemes for the diffusion asymptotics; then, their extension to the anomalous diffusion limit is studied. In this case, it is crucial to capture the effect of the large velocities of the heavy-tailed equilibrium, so that some important transformations of the schemes derived for the diffusion asymptotics are needed. As a result, we obtain numerical schemes which enjoy the Asymptotic Preserving property in the anomalous diffusion limit, that is: they do not suffer from the restriction on the time step and they degenerate towards the fractional diffusion limit when the mean free path goes to zero. We also numerically investigate the uniform accuracy and construct a class of numerical schemes satisfying this property. Finally, the efficiency of the different numerical schemes is shown through numerical experiments.

preprint2014arXiv

Models of dark matter halos based on statistical mechanics: I. The classical King model

We consider the possibility that dark matter halos are described by the Fermi-Dirac distribution at finite temperature. This is the case if dark matter is a self-gravitating quantum gas made of massive neutrinos at statistical equilibrium. This is also the case if dark matter can be treated as a self-gravitating collisionless gas experiencing Lynden-Bell's type of violent relaxation. In order to avoid the infinite mass problem and carry out a rigorous stability analysis, we consider the fermionic King model. In this paper, we study the non-degenerate limit leading to the classical King model. This model was initially introduced to describe globular clusters. We propose to apply it also to large dark matter halos where quantum effects are negligible. We determine the caloric curve and study the thermodynamical stability of the different configurations. Equilibrium states exist only above a critical energy $E_c$ in the microcanonical ensemble and only above a critical temperature $T_c$ in the canonical ensemble. For $E<E_c$, the system undergoes a gravothermal catastrophe and, for $T<T_c$, it undergoes an isothermal collapse. We compute the profiles of density, circular velocity, and velocity dispersion. We compare the prediction of the classical King model to the observations of large dark matter halos. Because of collisions and evaporation, the central density increases while the slope of the halo density profile decreases until an instability takes place. We show that large dark matter halos are relatively well-described by the King model at, or close to, the point of marginal microcanonical stability. At that point, the King model generates a density profile that can be approximated by the modified Hubble profile. This profile has a flat core and decreases as $r^{-3}$ at large distances, like the observational Burkert profile. Less steep halos are unstable.

preprint2014arXiv

Models of dark matter halos based on statistical mechanics: II. The fermionic King model

We discuss the nature of phase transitions in the fermionic King model which describes tidally truncated quantum self-gravitating systems. This distribution function takes into account the escape of high energy particles and has a finite mass. On the other hand, the Pauli exclusion principle puts an upper bound on the phase space density of the system and stabilizes it against gravitational collapse. As a result, there exists a statistical equilibrium state for any accessible values of energy and temperature. We plot the caloric curves and investigate the nature of phase transitions as a function of the degeneracy parameter in both microcanonical and canonical ensembles. We consider stable and metastable states and emphasize the importance of the latter for systems with long-range interactions. Phase transitions can take place between a "gaseous" phase unaffected by quantum mechanics and a "condensed" phase dominated by quantum mechanics. The phase diagram exhibits two critical points, one in each ensemble, beyond which the phase transitions disappear. There also exist a region of negative specific heats and a situation of ensemble inequivalence for sufficiently large systems. We apply the fermionic King model to the case of dark matter halos made of massive neutrinos. The gaseous phase describes large halos and the condensed phase describes dwarf halos. Partially degenerate configurations describe intermediate size halos. We argue that large dark matter halos cannot harbor a fermion ball because these nucleus-halo configurations are thermodynamically unstable (saddle points of entropy). Large dark matter halos may rather contain a central black hole resulting from a dynamical instability of relativistic origin occurring during the gravothermal catastrophe.

preprint2014arXiv

The fermionic King model

We study the fermionic King model which may provide a relevant model of dark matter halos. The exclusion constraint can be due to quantum mechanics (for fermions such as massive neutrinos) or to Lynden-Bell's statistics (for collisionless systems undergoing violent relaxation). This model has a finite mass. Furthermore, a statistical equilibrium state exists for all accessible values of energy. Dwarf and intermediate size halos are degenerate quantum objects stabilized against gravitational collapse by the Pauli exclusion principle. Large halos at sufficiently high energies are in a gaseous phase where quantum effects are negligible. They are stabilized by thermal motion. Below a critical energy they undergo gravitational collapse (gravothermal catastrophe). This may lead to the formation of a central black hole that does not affect the structure of the halo. This may also lead to the formation of a compact degenerate object surrounded by a hot massive atmosphere extending at large distances. We argue that large dark matter halos should not contain a degenerate nucleus (fermion ball) because these nucleus-halo structures are thermodynamically unstable. We compare the rotation curves of the classical King model to observations of large dark matter halos (Burkert profile). Because of collisions and evaporation, the central density increases while the slope of the halo density profile decreases until an instability takes place. We find that the observations are compatible with a King profile at, or close to, the point of marginal stability in the microcanonical ensemble. At that point, the King profile can be fitted by the modified Hubble profile that has a flat core and a halo in which the density decreases as $r^{-3}$. This is qualitatively similar to the Burkert profile.

preprint2013arXiv

Uniformly accurate numerical schemes for highly oscillatory Klein-Gordon and nonlinear Schrödinger equations

This work is devoted to the numerical simulation of nonlinear Schrödinger and Klein-Gordon equations. We present a general strategy to construct numerical schemes which are uniformly accurate with respect to the oscillation frequency. This is a stronger feature than the usual so called "Asymptotic preserving" property, the last being also satisfied by our scheme in the highly oscillatory limit. Our strategy enables to simulate the oscillatory problem without using any mesh or time step refinement, and the orders of our schemes are preserved uniformly in all regimes. In other words, since our numerical method is not based on the derivation and the simulation of asymptotic models, it works in the regime where the solution does not oscillate rapidly, in the highly oscillatory limit regime, and in the intermediate regime with the same order of accuracy. In the same spirit as in [5], the method is based on two main ingredients. First, we embed our problem in a suitable "two-scal" reformulation with the introduction of an additional variable. Then a link is made with classical strategies based on Chapman-Enskog expansions in kinetic theory despite the dispersive context of the targeted equations, allowing to separate the fast time scale from the slow one. Uniformly accurate (UA) schemes are eventually derived from this new formulation and their properties and performances are assessed both theoretically and numerically.

preprint2012arXiv

Asymptotic preserving schemes for highly oscillatory kinetic equation

This work is devoted to the numerical simulation of a Vlasov-Poisson model describing a charged particle beam under the action of a rapidly oscillating external electric field. We construct an Asymptotic Preserving numerical scheme for this kinetic equation in the highly oscillatory limit. This scheme enables to simulate the problem without using any time step refinement technique. Moreover, since our numerical method is not based on the derivation of the simulation of asymptotic models, it works in the regime where the solution does not oscillate rapidly, and in the highly oscillatory regime as well. Our method is based on a "double-scale" reformulation of the initial equation, with the introduction of an additional periodic variable.

preprint2012arXiv

Micro-macro schemes for kinetic equations including boundary layers

We introduce a new micro-macro decomposition of collisional kinetic equations in the specific case of the diffusion limit, which naturally incorporates the incoming boundary conditions. The idea is to write the distribution function $f$ in all its domain as the sum of an equilibrium adapted to the boundary (which is not the usual equilibrium associated with $f$) and a remaining kinetic part. This equilibrium is defined such that its incoming velocity moments coincide with the incoming velocity moments of the distribution function. A consequence of this strategy is that no artificial boundary condition is needed in the micro-macro models and the exact boundary condition on $f$ is naturally transposed to the macro part of the model. This method provides an 'Asymptotic preserving' numerical scheme which generates a very good approximation of the space boundary values at the diffusive limit, without any mesh refinement in the boundary layers. Our numerical results are in very good agreement with the exact so-called Chandrasekhar value, which is explicitely known in some simple cases.

preprint2012arXiv

Stable ground states and self-similar blow-up solutions for the gravitational Vlasov-Manev system

In this work, we study the orbital stability of steady states and the existence of blow-up self-similar solutions to the so-called Vlasov-Manev (VM) system. This system is a kinetic model which has a similar Vlasov structure as the classical Vlasov-Poisson system, but is coupled to a potential in $-1/r- 1/r^2$ (Manev potential) instead of the usual gravitational potential in $-1/r$, and in particular the potential field does not satisfy a Poisson equation but a fractional-Laplacian equation. We first prove the orbital stability of the ground states type solutions which are constructed as minimizers of the Hamiltonian, following the classical strategy: compactness of the minimizing sequences and the rigidity of the flow. However, in driving this analysis, there are two mathematical obstacles: the first one is related to the possible blow-up of solutions to the VM system, which we overcome by imposing a sub-critical condition on the constraints of the variational problem. The second difficulty (and the most important) is related to the nature of the Euler-Lagrange equations (fractional-Laplacian equations) to which classical results for the Poisson equation do not extend. We overcome this difficulty by proving the uniqueness of the minimizer under equimeasurabilty constraints, using only the regularity of the potential and not the fractional-Laplacian Euler-Lagrange equations itself. In the second part of this work, we prove the existence of exact self-similar blow-up solutions to the Vlasov-Manev equation, with initial data arbitrarily close to ground states. This construction is based on a suitable variational problem with equimeasurability constraint.

preprint2011arXiv

A boundary matching micro/macro decomposition for kinetic equations

We introduce a new micro/macro decomposition of collisional kinetic equations which naturally incorporates the exact space boundary conditions. The idea is to write the distribution fonction $f$ in all its domain as the sum of a Maxwellian adapted to the boundary (which is not the usual Maxwellian associated with $f$) and a reminder kinetic part. This Maxwellian is defined such that its 'incoming' velocity moments coincide with the 'incoming' velocity moments of the distribution function. Important consequences of this strategy are the following. i) No artificial boundary condition is needed in the micro/macro models and the exact boundary condition on $f$ is naturally transposed to the macro part of the model. ii) It provides a new class of the so-called 'Asymptotic preserving' (AP) numerical schemes: such schemes are consistent with the original kinetic equation for all fixed positive value of the Knudsen number $\eps$, and if $\eps \to 0 $ with fixed numerical parameters then these schemes degenerate into consistent numerical schemes for the various corresponding asymptotic fluid or diffusive models. Here, the strategy provides AP schemes not only inside the physical domain but also in the space boundary layers. We provide a numerical test in the case of a diffusion limit of the one-group transport equation, and show that our AP scheme recovers the boundary layer and a good approximation of the theoretical boundary value, which is usually computed from to the so-called Chandrasekhar function.

preprint2010arXiv

A new variational approach to the stability of gravitational systems

We consider the three dimensional gravitational Vlasov Poisson system which describes the mechanical state of a stellar system subject to its own gravity. A well-known conjecture in astrophysics is that the steady state solutions which are nonincreasing functions of their microscopic energy are nonlinearly stable by the flow. This was proved at the linear level by several authors based on the pioneering work by Antonov in 1961. Since then, standard variational techniques based on concentration compactness methods as introduced by P.-L. Lions in 1983 have led to the nonlinear stability of subclasses of stationary solutions of ground state type. In this paper, inspired by pioneering works from the physics litterature (Lynden-Bell 94, Wiechen-Ziegler-Schindler MNRAS 88, Aly MNRAS 89), we use the monotonicity of the Hamiltonian under generalized symmetric rearrangement transformations to prove that non increasing steady solutions are local minimizer of the Hamiltonian under equimeasurable constraints, and extract compactness from suitable minimizing sequences. This implies the nonlinear stability of nonincreasing anisotropic steady states under radially symmetric perturbations.

preprint2010arXiv

Escape of stars from gravitational clusters in the Chandrasekhar model

We study the evaporation of stars from globular clusters using the simplified Chandrasekhar model. This is an analytically tractable model giving reasonable agreement with more sophisticated models that require complicated numerical integrations. In the Chandrasekhar model: (i) the stellar system is assumed to be infinite and homogeneous (ii) the evolution of the velocity distribution of stars f(v,t) is governed by a Fokker-Planck equation, the so-called Kramers-Chandrasekhar equation (iii) the velocities |v| that are above a threshold value R>0 (escape velocity) are not counted in the statistical distribution of the system. In fact, high velocity stars leave the system, due to free evaporation or to the attraction of a neighboring galaxy (tidal effects). Accordingly, the total mass and energy of the system decrease in time. If the star dynamics is described by the Kramers-Chandrasekhar equation, the mass decreases to zero exponentially rapidly. Our goal is to obtain non-perturbative analytical results that complement the seminal studies of Chandrasekhar, Michie and King valid for large times $t\to+\infty$ and large escape velocities $R\to +\infty$. In particular, we obtain an exact semi-explicit solution of the Kramers-Chandrasekhar equation with the absorbing boundary condition f(R,t)=0. We use it to obtain an explicit expression of the mass loss at any time t when $R\to +\infty$. We also derive an exact integral equation giving the exponential evaporation rate $λ(R)$, and the corresponding eigenfunction $f_λ(v)$, when $t\to +\infty$ for any sufficiently large value of the escape velocity R. For $R\to +\infty$, we obtain an explicit expression of the evaporation rate that refines the Chandrasekhar results.

preprint2010arXiv

Orbital stability of spherical galactic models

We consider the three dimensional gravitational Vlasov Poisson system which is a canonical model in astrophysics to describe the dynamics of galactic clusters. A well known conjecture is the stability of spherical models which are nonincreasing radially symmetric steady states solutions. This conjecture was proved at the linear level by several authors in the continuation of the breakthrough work by Antonov in 1961. In a previous work (arXiv:0904.2443), we derived the stability of anisotropic models under {\it spherically symmetric perturbations} using fundamental monotonicity properties of the Hamiltonian under suitable generalized symmetric rearrangements first observed in the physics litterature. In this work, we show how this approach combined with a {\it new generalized} Antonov type coercivity property implies the orbital stability of spherical models under general perturbations.

preprint2009arXiv

Stable ground states for the relativistic gravitational Vlasov-Poisson system

We consider the three dimensional gravitational Vlasov-Poisson (GVP) system in both classical and relativistic cases. The classical problem is subcritical in the natural energy space and the stability of a large class of ground states has been derived by various authors. The relativistic problem is critical and displays finite time blow up solutions. Using standard concentration compactness techniques, we however show that the breaking of the scaling symmetry allows the existence of stable relativistic ground states. A new feature in our analysis which applies both to the classical and relativistic problem is that the orbital stability of the ground states does not rely as usual on an argument of uniqueness of suitable minimizers --which is mostly unknown-- but on strong rigidity properties of the transport flow, and this extends the class of minimizers for which orbital stability is now proved.