Researcher profile

Maxime Herda

Maxime Herda contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2022arXiv

Long-time behaviour of hybrid finite volume schemes for advection-diffusion equations: linear and nonlinear approaches

We are interested in the long-time behaviour of approximate solutions to heterogeneous and anisotropic linear advection-diffusion equations in the framework of hybrid finite volume (HFV) methods on general polygonal/polyhedral meshes. We consider two linear methods, as well as a new, nonlinear scheme, for which we prove the existence and the positivity of discrete solutions. We show that the discrete solutions to the three schemes converge exponentially fast in time towards the associated discrete steady-states. To illustrate our theoretical findings, we present some numerical simulations assessing long-time behaviour and positivity. We also compare the accuracy of the schemes on some numerical tests in the stationary case.

preprint2020arXiv

A Fokker-Planck approach to the study of robustness in gene expression

We study several Fokker-Planck equations arising from a stochastic chemical kinetic system modeling a gene regulatory network in biology. The densities solving the Fokker-Planck equations describe the joint distribution of the messenger RNA and micro RNA content in a cell. We provide theoretical and numerical evidences that the robustness of the gene expression is increased in the presence of micro RNA. At the mathematical level, increased robustness shows in a smaller coefficient of variation of the marginal density of the messenger RNA in the presence of micro RNA. These results follow from explicit formulas for solutions. Moreover, thanks to dimensional analyses and numerical simulations we provide qualitative insight into the role of each parameter in the model. As the increase of gene expression level comes from the underlying stochasticity in the models, we eventually discuss the choice of noise in our models and its influence on our results.

preprint2020arXiv

A note on hypocoercivity for kinetic equations with heavy-tailed equilibrium

In this paper we are interested in the large time behavior of linear kinetic equations with heavy-tailed local equilibria. Our main contribution concerns the kinetic Lévy-Fokker-Planck equation, for which we adapt hypocoercivity techniques in order to show that solutions converge exponentially fast to the global equilibrium. Compared to the classical kinetic Fokker-Planck equation, the issues here concern the lack of symmetry of the non-local Lévy-Fokker-Planck operator and the understanding of its regularization properties. As a complementary related result, we also treat the case of the heavy-tailed BGK equation.

preprint2020arXiv

Computation of sum of squares polynomials from data points

We propose an iterative algorithm for the numerical computation of sums of squares of polynomials approximating given data at prescribed interpolation points. The method is based on the definition of a convex functional $G$ arising from the dualization of a quadratic regression over the Cholesky factors of the sum of squares decomposition. In order to justify the construction, the domain of $G$, the boundary of the domain and the behavior at infinity are analyzed in details. When the data interpolate a positive univariate polynomial, we show that in the context of the Lukacs sum of squares representation, $G$ is coercive and strictly convex which yields a unique critical point and a corresponding decomposition in sum of squares. For multivariate polynomials which admit a decomposition in sum of squares and up to a small perturbation of size $\varepsilon$, $G^\varepsilon$ is always coercive and so it minimum yields an approximate decomposition in sum of squares. Various unconstrained descent algorithms are proposed to minimize $G$. Numerical examples are provided, for univariate and bivariate polynomials.

preprint2020arXiv

Large time behavior of nonlinear finite volume schemes for convection-diffusion equations

In this contribution we analyze the large time behavior of a family of nonlinear finite volume schemes for anisotropic convection-diffusion equations set in a bounded bidimensional domain and endowed with either Dirichlet and / or no-flux boundary conditions. We show that solutions to the two-point flux approximation (TPFA) and discrete duality finite volume (DDFV) schemes under consideration converge exponentially fast toward their steady state. The analysis relies on discrete entropy estimates and discrete functional inequalities. As a biproduct of our analysis, we establish new discrete Poincaré-Wirtinger, Beckner and logarithmic Sobolev inequalities. Our theoretical results are illustrated by numerical simulations.

preprint2019arXiv

A projection algorithm on the set of polynomials with two bounds

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.

preprint2018arXiv

Anisotropic Boltzmann-Gibbs dynamics of strongly magnetized Vlasov-Fokker-Planck equations

We consider various sets of Vlasov-Fokker-Planck equations modeling the dynamics of charged particles in a plasma under the effect of a strong magnetic field. For each of them in a regime where the strength of the magnetic field is effectively stronger than that of collisions we first formally derive asymptotically reduced models. In this regime, strong anisotropic phenomena occur; while equilibrium along magnetic field lines is asymptotically reached our asymptotic models capture a non trivial dynamics in the perpendicular directions. We do check that in any case the obtained asymptotic model defines a well-posed dynamical system and when self consistent electric fields are neglected we provide a rigorous mathematical justification of the formally derived systems. In this last step we provide a complete control on solutions by developing anisotropic hypocoercive estimates.

preprint2017arXiv

Large-time behavior of solutions to Vlasov-Poisson-Fokker-Planck equations: from evanescent collisions to diffusive limit

The present contribution investigates the dynamics generated by the two-dimensional Vlasov-Poisson-Fokker-Planck equation for charged particles in a steady inhomogeneous background of opposite charges. We provide global in time estimates that are uniform with respect to initial data taken in a bounded set of a weighted $L^2$ space, and where dependencies on the mean-free path $τ$ and the Debye length $δ$ are made explicit. In our analysis the mean free path covers the full range of possible values: from the regime of evanescent collisions $τ\to\infty$ to the strongly collisional regime $τ\to0$. As a counterpart, the largeness of the Debye length, that enforces a weakly nonlinear regime, is used to close our nonlinear estimates. Accordingly we pay a special attention to relax as much as possible the $τ$-dependent constraint on $δ$ ensuring exponential decay with explicit $τ$-dependent rates towards the stationary solution. In the strongly collisional limit $τ\to0$, we also examine all possible asymptotic regimes selected by a choice of observation time scale. Here also, our emphasis is on strong convergence, uniformity with respect to time and to initial data in bounded sets of a $L^2$ space. Our proofs rely on a detailed study of the nonlinear elliptic equation defining stationary solutions and a careful tracking and optimization of parameter dependencies of hypocoercive/hypoelliptic estimates.

preprint2016arXiv

On massless electron limit for a multispecies kinetic system with external magnetic field

We consider a three-dimensional kinetic model for a two species plasma consisting of electrons and ions confined by an external nonconstant magnetic field. Then we derive a kinetic-fluid model when the mass ratio $m_e/m_i$ tends to zero. Each species initially obeys a Vlasov-type equation and the electrostatic coupling follows from a Poisson equation. In our modeling, ions are assumed non-collisional while a Fokker-Planck collision operator is taken into account in the electron equation. As the mass ratio tends to zero we show convergence to a new system where the macroscopic electron density satisfies an anisotropic drift-diffusion equation. To achieve this task, we overcome some specific technical issues of our model such as the strong effect of the magnetic field on electrons and the lack of regularity at the limit. With methods usually adapted to diffusion limit of collisional kinetic equations and including renormalized solutions, relative entropy dissipation and velocity averages, we establish the rigorous derivation of the limit model.