Maxime Herda

Francis Filbet, Maxime Herda

We propose a finite volume scheme for a class of nonlinear parabolic equations endowed with non-homogeneous Dirichlet boundary conditions and which admit relative en-tropy functionals. For this kind of models including porous media equations, Fokker-Planck equations for plasma physics or dumbbell models for polymer flows, it has been proved that the transient solution converges to a steady-state when time goes to infinity. The present scheme is built from a discretization of the steady equation and preserves steady-states and natural Lyapunov functionals which provide a satisfying long-time behavior. After proving well-posedness, stability, exponential return to equilibrium and convergence, we present several numerical results which confirm the accuracy and underline the efficiency to preserve large-time asymptotic.

Martin Campos Pinto, Frédérique Charles, Bruno Després et al.

The motivation of this work stems from the numerical approximation of bounded functions by polynomials satisfying the same bounds. The present contribution makes use of the recent algebraic characterization found in [B. Després, Numer. Algorithms, 76(3), (2017)] and [B. Després and M. Herda, Numer. Algorithms, 77(1), (2018)] where an interpretation of monovariate polynomials with two bounds is provided in terms of a quaternion algebra and the Euler four-squares formulas. Thanks to this structure, we generate a new nonlinear projection algorithm onto the set of polynomials with two bounds. The numerical analysis of the method provides theoretical error estimates showing stability and continuity of the projection. Some numerical tests illustrate this novel algorithm for constrained polynomial approximation.

Claire Chainais-Hillairet, Maxime Herda

We are interested in the large-time behavior of solutions to finite volume discretizations of convection-diffusion equations or systems endowed with non-homogeneous Dirichlet and Neumann type boundary conditions. Our results concern various linear and nonlinear models such as Fokker-Planck equations, porous media equations or drift-diffusion systems for semiconductors. For all of these models, some relative entropy principle is satisfied and implies exponential decay to the stationary state. In this paper we show that in the framework of finite volume schemes on orthogonal meshes, a large class of two-point monotone fluxes preserve this exponential decay of the discrete solution to the discrete steady state of the scheme. This includes for instance upwind and centered convections or Scharfetter-Gummel discretizations. We illustrate our theoretical results on several numerical test cases.