Asymptotic Behavior of systems of PDE arising in physics and biology: theoretical and numerical points of view (ABPDE II)

15 juin 2016, 09:00 → 17 juin 2016, 16:00 Europe/Paris

Alexandre MOUTON (Laboratoire Paul Painlevé), Claire CHAINAIS-HILLAIRET (Laboratoire Paul Painlevé), Clément CANCES (INRIA Lille), Felix Ali Mehmeti (LAMAV), Ingrid LACROIX-VIOLET (Laboratoire Paul Painlevé), Thomas REY (Laboratoire Paul Painlevé)

Aim and scope

The main goals of this workshop are the theoretical study of asymptotic behaviors (in large time or with respect to some parameters) of problems arising in physics and biology and the development of asymptotic preserving numerical methods.

The second edition of this workshop features ten plenary speakers. In addition, several contributed talks and a poster session will complete the program.

 

 

Program

10 plenary speakers are confirmed :

• Marianne Bessemoulin-Chatard (Nantes, France)
• Rémi Carles (Montpellier, France)
• Frédérique Charles (Paris, France)
• Giacomo Dimarco (Ferrara, Italy)
• Klemens Fellner (Graz, Austria)
• Mohammed Lemou (Rennes, France)
• Serge Nicaise (Valenciennes, France)
• Magali Ribot (Orléans, France)
• Jan ten Thije Boonkkamp (Eindhoven, Netherlands)
• Alexis Vasseur (Austin, Texas)

The complete program is available here. A printable version of the program is downloadable here.

 

Registration

Registrations are now closed.

Funding

Limited funding for the local expenses of students and young researchers is available. If you wish to apply for such support, please register and send a CV and a publication list by email to Thomas Rey.

 

A participant list is available here.

Transparents complete_programme.pdf

mer. 15 juin

09:30 → 10:10 Welcome

10:10 → 10:55 Global Existence of Solutions to the 3D Navier-Stokes Equations with Degenerate Viscosities

We prove the existence of global weak solutions for 3D compressible Navier-Stokes equations with degenerate viscosities. The method is based on the Bresch and Desjardins entropy. The solutions are obtained as limits of the quantic Navier-Stokes system. The main contribution is to derive the Mellet-Vasseur type inequality for the weak solutions, even if it is not verified by the first level of approximation. This provides existence of global solutions in time, for the compressible Navier-Stokes equations, for any gamma bigger than one, in three dimensional space, with large initial data, possibly vanishing on the vacuum. This is a joint work with Cheng Yu. The paper will appear in Inventiones.

Orateur: Alexis Vasseur (Department of Mathematics)

10:55 → 11:30 Coffee break

11:30 → 12:15 A (mainly numerical) study of a hyperbolic model for chemotaxis

The aim of this talk is to give some first results on the behaviour of the solutions of a 1D hyperbolic type chemotaxis system, based on incompressible Euler equation. More precisely, I will completely describe the stationary solutions with vacuum for this system and I will study numerically the stability of these steady states after the presentation of an adapted numerical scheme. A comparison with a limit parabolic system will also be performed.

Orateur: Magali Ribot (Laboratoire de Mathématiques - Analyse, Probabilités, Modélisation - Orléans)

12:15 → 12:50 Modelling and numerical approximation for the nonconservative bi-temperature MHD model

In order to achieve inertial confinement fusion, one has to improve the knowledge of the laser-plasma interaction. There exists two main ways of describing this phenomenon, the microscopic (kinetic) approach and the macroscopic (hydrodynamic) approach. The kinetic approach is not competitive since it is too expensive in computational time. This is why we investigate an intermediate model in thermal nonequilibrium, which is between the kinetic model and the hydrodynamic model. In the first stage of the confinement the magnetic field is negligible, the relevant intermediate model is than the nonconservative bitemperature Euler model. Recently in [1], an entropic approximation of this system has been derived thanks to numerical schemes based on an underlying conservative kinetic model. However in the last stage of the confinement the target is penetrated by relativistic electrons, which induces a strongly variable magnetic field. This is why we want to deal with an intermediate model which takes into account the magnetic field. In this work we propose to study a bitemperature MHD model. This system consists in four conservation equations for mass, impulsion and magnetic field and two nonconservation equations, that is to say, one for each energy. Physically, the model describes the interaction of a mixture of one species of ions and one species of electrons in thermal nonequilibrium subjected to a transverse variable magnetic field. A first result is to have been able to established the hydrodynamic model from an underlying kinetic model. More precisely, using an out of equilibrium Chapman-Enskop procedure, the bitemperature MHD model is constructed from a BGK model coupled with Maxwell equations with full Lorentz force, which includes the magnetic field. Finally, we approximate the weak solutions of the bitemperature MHD model using a kinetic scheme, based on the underlying kinetic model. References [1] D. Aregba-Driollet, S. Brull, J.Breil, B. Dubroca and E. Estibal, Modelling and numerical ap- proximation for the nonconservative bitemperature Euler model, preprint, 2015.

Orateur: Xavier Lhebrard (Centre Lasers Intenses et Applications)

12:50 → 14:35 Lunch

14:35 → 15:20 A multiscale numerical approach for a class of time-space oscillatory problems

High oscillations may arise in many physical problems: Schrödinger equations, kinetic equations, or more generally high frequency waves. In this talk, we will present a general strategy that allows the construction of uniformly (with respect to the oscillation frequency) accurate numerical schemes in the following situations: i) time oscillations with applications to kinetic and Schrödinger equations. ii) time-space oscillations with applications to some high frequency waves and semi-classical quantum models. Some numerical tests will be presented to illustrate the efficiency of the strategy.

Orateur: Mohammed Lemou (Institut de Recherche Mathématiques de Rennes)

15:20 → 15:55 Global existence for small data of the viscous Green-Naghdi equations

We consider the Cauchy problem for the Green-Naghdi equations with viscosity, for small initial data. It is well-known that adding a second order dissipative term to a hyperbolic system leads to the existence of global smooth solutions, once the hyperbolic system is symmetrizable and the so-called Kawashima-Shizuta condition is satisfied. We first show that the Green-Naghdi equations can be written in a symmetric form, using the associated Hamiltonian. This system being dispersive, in the sense that it involves third order derivatives, the symmetric form is based on symmetric differential operators. Then, we use this structure for an appropriate change of variable to prove that adding viscosity effects through a second order term leads to global existence of smooth solutions, for small data. We also deduce that constant solutions are asymptotically stable.

Orateur: Dena Kazerani (Laboratoire Jacques-Louis Lions)

15:55 → 16:30 Coffee break

16:30 → 17:05 Convergence to equilibrium for gradient-like systems with analytic features

A celebrated result of S. Lojasiewicz states that every bounded solution of a gradient flow associated to an analytic function converges to a steady state as time goes to infinity. Convergence rates can also be obtained. These convergence results have been generalized to a large variety of finite or infinite dimensional gradient-like flows. The fundamental example in infinite dimension is the semilinear heat equation with an analytic nonlinearity. In this talk, we show how some of these results can be adapted to time discretizations of gradient-like flows, in view of applications to PDEs such as the Allen-Cahn equation, the sine-Gordon equation, the Cahn-Hilliard equation, the Swift-Hohenberg equation, or the phase-field crystal equation.

Orateur: Morgan Pierre (Laboratoire de Mathématiques et Applications)

17:05 → 17:40 Numerical convergence rate for the diffusive limit of the p-system with damping

We are interested in the study of the diffusive limit of the $p$-system with damping and its approximation by an Asymptotic Preserving (AP) Finite Volume scheme. Provided the system is endowed with an entropy-entropy flux pair, we give the convergence rate of classical solutions of the p-system with damping towards the smooth solutions of the porous media equation using a relative entropy method. Adopting a semi-discrete scheme, we establish that the convergence rate is preserved by the approximated solutions. Several numerical experiments illustrate the relevance of this result.

Orateur: Hélène Mathis (Laboratoire de Mathématiques Jean Leray)

jeu. 16 juin

08:45 → 09:30 Coffee

09:30 → 10:15 Implicit-explicit linear multistep methods for stiff kinetic equations

We consider the development of high order asymptotic-preserving linear multistep methods for kinetic equations and related problems. The methods are first developed for BGK-like kinetic models and then extended to the case of the full Boltzmann equation. The behavior of the schemes in the Navier-Stokes regime is also studied and compatibility conditions derived. We show that, compared to IMEX Runge-Kutta methods, the IMEX multistep schemes have several advantages due to the absence of coupling conditions and to the greater computational efficiency. The latter is of paramount importance when dealing with the time discretization of multidimensional kinetic equations.

Orateur: Giacomo Dimarco (Department of Mathematics and Computer Science)

10:15 → 11:00 From particle methods to hybrid semi-Lagrangian schemes

Particle methods for transport equations consist in pushing forward particles along the characteristic lines of the flow, and to describe then the transported density as a sum of weighted and smoothed particles. Conceptually simple, standard particle methods have the main drawback to produce noisy solutions or to require frequent remapping. In this talk we present two classes of particle methods which aim at improving the accuracy of the numerical approximations with a minimal amount of smoothing. The idea of the Linearly Transformed Particle method is to transform the shape functions of particles in order to follow the local variation of the flow. This method has been adapted and analyzed for the Vlasov- Poisson system and for a compressible aggregation equation. In both cases the error estimate is improved compared to classical particle methods, with the gain of a strong convergence of the numerical solution. However, for long remapping periods, shapes of particles could become to much stretched out. The second method solve this problem of locality by combining a backward semi-Lagrangian approach and local linearizations of the flow. The convergence properties are improved and validated by numerical experiments. This is a joint work with Martin Campos-Pinto (LJLL, UPMC).

Orateur: Frédérique Charles (Laboratoire Jacques-Louis Lions)

11:00 → 11:30 Coffee break

11:30 → 12:15 Time splitting methods and the semi-classical limit for nonlinear Schrödinger equations

We consider the time discretization based on Lie-Trotter splitting, for the nonlinear Schrödinger equation, in the semi-classical limit, with initial data under the form of WKB states. Both the exact and the numerical solutions keep a WKB structure, on a time interval independent of the Planck constant. We prove error estimates, which show that the quadratic observables can be computed with a time step independent of the Planck constant. We give a flavor of the functional framework, based on time-dependent analytic spaces.

Orateur: Rémi Carles (Institut Montpelliérain Alexander Grothendieck)

12:15 → 12:50 Dimensional reduction of a multiscale model based on long time asymptotics

Depending on their velocity field, some models lead to moment equations that enable one to compute monokinetic solutions economically. We detail the example of a multiscale structured cell population model, consisting of a system of 2D transport equations. The reduced model, a system of 1D transport equations, is obtained by computing the moments of the 2D model with respect to one variable. The 1D solution is defined from the solution of the 2D model starting from an initial condition that is a Dirac mass in the direction removed by reduction. Long time properties of the 1D model solution are obtained in connection with properties of the support of the 2D solution for general case initial conditions. Finite volume numerical approximations of the 1D reduced model can be used to compute the moments of the 2D solution with proper accuracy. The numerical robustness is studied in the scalar case, and a full scale vector case is presented.

Orateur: Marie Postel (Laboratoire Jacques-Louis Lions)

12:50 → 14:40 Lunch

14:40 → 15:25 Global existence and large-time behaviour for reaction-diffusion models

Systems of nonlinear reaction-diffusion equations are encountered frequently as models in chemistry, physics, populations dynamics and biology. However, due to the lack of comparison principles for general reaction-diffusion systems, already the existence of global weak/classical solutions poses many open problems, in particular in 3D. In the absence of comparison principles, so called duality methods have recently proven to be one of the most powerful tools in obtaining global solutions for nonlinear reaction-diffusion systems. The first part of this talk will present recent advances and results concerning the existence of global solutions via duality methods. The second part of the talk will then consider reaction-diffusion systems, which feature an entropy functional and discuss the convergence to equilibrium states with computable rates for large classes of such reaction-diffusion models.

Orateur: Klemens Fellner (Institut für Mathematik und Wissenschaftliches Rechnen)

15:25 → 16:00 Asymptotic analysis for a simplified model of model of dynamical perfect plasticity

In this talk, we will present an initial boundary value problem for a hyperbolic system under constraints, coming from mechanics. To study the solutions of such a system, we will use a viscous approach that relaxes the constraints. We will explain the asymptotic analysis, when the viscous parameter tends to zero, which leads to an interaction between the boundary condition and the constraints for the constrained system. If time permits, we will show some numerical results.

Orateur: Clément Mifsud (Laboratoire Jacques Louis Lions)

16:00 → 16:30 Coffee break

16:00 → 18:00 Poster session

Orateurs: Ahmed Ait Hammou Oulhaj (Laboratoire Paul Painlevé), An Zhang (CEREMADE), Andrea Bondesan (Laboratoire MAP5), Judith Berendsen (Institute for Computational and Applied Mathematics), Maxime Herda (Institut Camille Jordan), Samia Zermani (Institut aux Etudes d'Ingénieurs el Manar)

19:30 → 22:30 Gala dinner

Google Maps url : http://tinyurl.com/hoyqyp5

ven. 17 juin

08:45 → 09:30 Coffee

09:30 → 10:15 Stability results of dissipative systems via the frequency domain approach

The frequency domain approach goes back to J. Prüss [Trans. Amer. Math. Soc. 284 (1984), 847-857] and F. L. Huang [Ann. Differential Equations 1 (1985), 43-56] that show that a $C_0$ semigroup $(e^{tA})_{t\geq 0}$ of contractions in a Hilbert space $H$ is exponentially stable if and only if the resolvent of $A$ is uniformly bounded on the imaginary axis. Afterwards Z. Liu and B. Rao [Z. Angew. Math. Phys. 56 (2005), 630-644], C. J. K. Batty and T. Duyckaerts [J. Evol. Equ. 8 (2008), 765-780], and A. Bátkai, K.-J. Engel, J. Prüss and R. Schnaubelt [Math. Nachr. 279 (2006), 1425-1440] have given some sufficient conditions on the behavior of the resolvent of $A$ on the imaginary axis that guarantee an almost polynomial decay of the semigroup. Finally an optimal result about the polynomial decay was found by A. A. Borichev and Yu. V. Tomilov [Math. Ann. 347 (2010), 455-478]. This approach is a powerful tool for the study of the decay of the semigroup associated with concrete dissipative systems since it reduces to the study of the resolvent on the imaginary axis. In our talk, we will first recall these two results and then illustrate them on two particular dissipative systems, namely a generalized telegraph equation [Z. Angew. Math. Phys. 66 (2015), 3221-3247] and a dispersive medium model (joint work with C. Scheid (Univ. Nice)).

Orateur: Serge Nicaise (Laboratoire de Mathématiques et de leurs Applications de Valenciennes)

10:15 → 10:50 Uniform asymptotic preserving scheme for hyperbolic systems in 2D

In this work, we are interested by the discretization of hyperbolic system with stiff source term. Firstly we consider a simple linear case : the damped wave equation which can be approximative by a diffusion equation at the limit. For this equation we propose a asymptotic preserving scheme which converge uniformly on general and unstructured 2D meshes contrary to the classical extension of the AP which are not consistent in the limit regime on unstructured meshes. After that we propose to extend this method to a nonlinear problem: the Euler equations with friction. At the end the link with the well-balanced scheme (for Euler-Poisson) will be introduced.

Orateur: Emmanuel Franck (Inria Nancy Grand-est)

10:50 → 11:25 Coffee break

11:25 → 12:10 Complete flux schemes for conservation laws of advection-diffusion-reaction type

Complete flux schemes are recently developed numerical flux approximation schemes for conservation laws of advection-diffusion-reaction type; see e.g. [1, 2]. The basic complete flux scheme is derived from a local one-dimensional boundary value problem for the entire equation, including the source term. Consequently, the integral representation of the flux contains a homogeneous and an inhomogeneous part, corresponding to the advection-diffusion operator and the source term, respectively. Suitable quadrature rules give the numerical flux. For time-dependent problems, the time derivative is considered a source term and is included in the inhomogeneous flux, resulting in an implicit semi-discretisation. The implicit system proves to have much smaller dissipation and dispersion errors than the standard semidiscrete system, especially for dominant advection. Just as for scalar equations, for coupled systems of conservation laws, the complete flux approximation is derived from a local system boundary value problem, this way incorporatin the coupling between the constituent equations in the discretization. Also in the system case, the numerical flux (vector) is the superpostion of a homogeneous and an inhomogeneous component, corresponding to the advection-diffusion operator and the source term vector, respectively. The scheme is applied to multi-species diffusion and satisfies the mass constraint exactly. References [1] J.H.M. ten Thije Boonkkamp and M.J.H. Anthonissen, The finite volume-complete flux scheme for advection-diffusion-reaction equations, J. Sci. Comput. 46, pp. 47-70 (2011). [2] J.H.M. ten Thije Boonkkamp, J. van Dijk, L. Liu and K.S.C. Peerenboom, Extension of the complete flux scheme to systems of comservation laws, J. Sci. Comput. 53, pp. 552-568 (2012).

Orateur: Jan ten Thije Boonkkamp (Department of Mathematics and Computer Science)

12:10 → 12:45 Non linear stability of Minkowski space-time with massive scalar field

In this talk we will present some recent work about the system of Einstein equation coupled with a massive scalar field and the system of $f(R)$ field equation (partially published in [2]). More precisely, on the nonlinear global stability of the Minkowski space-time within these two similar contexts. In a PDE point of view, they are equivalent to the global existence of a special class of quasi-linear wave-Klein-Gordon system with small initial data. To the author’s knowledge there is not so much choice to deal with this kind of system (for a detailed explication of the major difficulty, see for example in [1] page 2), and we apply the “hyperboloidal foliation method” introduced by the author in [1] combined with some newly developed tools such as $L^∞$ estimates on Klein-Gordon equations in curved space-time and $L^∞$ estimates on wave equations based on the expression of spherical means. We also adapt some tools developed in classical framework for the analysis of Einstein equation into our hyperboloidal foliation framework, such as the estimates based on wave gauge conditions and the L$^∞$ estimates on wave equations based on integration along characteristics. References [1] P. LeFloch and Y. Ma, The hyperboloidal foliation method, World Scientific, 2015 [2] P. LeFloch and Y. Ma, The nonlinear stability of Minkowski space for self-gravitating massive field. The wave-Klein-Gordon model, Comm. Math. Phys., published online.

Orateur: Yue Ma (School of Mathematics and Statistics)

12:45 → 14:40 Lunch

14:40 → 15:25 Exponential decay of a finite volume scheme to the thermal equilibrium for drift-diffusion systems

We are interested in the large-time behavior of a numerical scheme discretizing drift-diffusion systems for semiconductors. The considered scheme is finite volume in space, and the numerical fluxes are a generalization of the classical Scharfetter-Gummel scheme, which allows to consider both linear or nonlinear pressure laws. We study the convergence of approximate solutions towards an approximation of the thermal equilibrium state as time tends to infinity, and obtain a decay rate by controlling the discrete relative entropy with the entropy production. This result is proved under assumptions of existence and uniform-in-time $L^\infty$ estimates for numerical solutions, which will be discussed. This is a joined work with Claire Chainais-Hillairet.

Orateur: Marianne Bessemoulin (Laboratoire de Mathématiques Jean Leray)

15:25 → 16:00 Entropy methods for degenerate diffusions and weighted functional inequalities

We will present results on large time asymptotics for some fast diffusion equations with power law weights. We will show that, for such diffusions, new phenomena appear : the asymptotic rates of convergence, obtained by linearization, are not global, the underlying functional inequalities may experience symmetry breaking and the Barenblatt self-similar profiles is not optimal.

Orateur: Bruno Nazaret (SAMM - Statistique, Analyse, Modélisation multidisciplinaire)