Researcher profile

Marc Massot

Marc Massot contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
8topics
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

12 published item(s)

preprint2022arXiv

Does the multiresolution lattice Boltzmann method allow to deal with waves passing through mesh jumps?

We consider an adaptive multiresolution-based lattice Boltzmann scheme, which we have recently introduced and studied from the perspective of the error control and the theory of the equivalent equations. This numerical strategy leads to high compression rates, error control and its high accuracy has been explained on uniform and dynamically adaptive grids. However, one key issue with non-uniform meshes within the framework of lattice Boltzmann schemes is to properly handle acoustic waves passing through a level jump of the grid. It usually yields spurious effects, in particular reflected waves. In this paper, we propose a simple mono-dimensional test-case for the linear wave equation with a fixed adapted mesh characterized by a potentially large level jump. We investigate this configuration with our original strategy and prove that we can handle and control the amplitude of the reflected wave, which is of fourth order in the space step of the finest mesh. Numerical illustrations show that the proposed strategy outperforms the existing methods in the literature and allow to assess the ability of the method to handle the mesh jump properly.

preprint2022arXiv

Multiresolution-based mesh adaptation and error control for lattice Boltzmann methods with applications to hyperbolic conservation laws

Lattice Boltzmann Methods (LBM) stand out for their simplicity and computational efficiency while offering the possibility of simulating complex phenomena. While they are optimal for Cartesian meshes, adapted meshes have traditionally been a stumbling block since it is difficult to predict the right physics through various levels of meshes. In this work, we design a class of fully adaptive LBM methods with dynamic mesh adaptation and error control relying on multiresolution analysis. This wavelet-based approach allows to adapt the mesh based on the regularity of the solution and leads to a very efficient compression of the solution without loosing its quality and with the preservation of the properties of the original LBM method on the finest grid. This yields a general approach for a large spectrum of schemes and allows precise error bounds, without the need for deep modifications on the reference scheme. An error analysis is proposed. For the purpose of assessing the approach, we conduct a series of test-cases for various schemes and scalar and systems of conservation laws, where solutions with shocks are to be found and local mesh adaptation is especially relevant. Theoretical estimates are retrieved while a reduced memory footprint is observed. It paves the way to an implementation in a multi-dimensional framework and high computational efficiency of the method for both parabolic and hyperbolic equations, which is the subject of a companion paper.

preprint2022arXiv

Numerical challenges in the simulation of 1D bounded low-temperature plasmas with charge separation in various collisional regimes

We study a 1D geometry of a plasma confined between two conducting floating walls with applications to laboratory plasmas. These plasmas are characterized by a quasi-neutral bulk that is joined to the wall by a thin boundary layer called sheath that is positively charged. Although analytical solutions are available in the sheath and the pre-sheath, joining the two areas by one analytical solution is still an open problem which requires the numerical resolution of the fluid equations coupled to Poisson equation. Current numerical schemes use high-order discretizations to correctly capture the electron current in the sheath, presenting unsatisfactory results in the boundary layer and they are not adapted to all the possible collisional regimes. In this work, we identify the main numerical challenges that arise when attempting the simulations of such configuration and we propose explanations for the observed phenomena via numerical analysis. We propose a numerical scheme with controlled diffusion as well as new discrete boundary conditions that address the identified issues.

preprint2022arXiv

Second-order uniformly asymptotic-preserving space-time-ImEx schemes for hyperbolic balance laws with stiff relaxation

We consider hyperbolic systems of conservation laws with relaxation source terms leading to a diffusive asymptotic limit under a parabolic scaling. We introduce a new class of secondorder in time and space numerical schemes, which are uniformly asymptotic preserving schemes. The proposed Implicit-Explicit (ImEx) approach, does not follow the usual path relying on the method of lines, either with multi-step methods or Runge-Kutta methods, or semi-discretized in time equations, but is inspired from the Lax-Wendroff approach with the proper level of implicit treatment of the source term. As a result, it yields a very compact stencil in space and time and we are able to rigorously show that both the second-order accuracy and the stability conditions are independent of the fast scales in the asymptotic regime, including the study of boundary conditions. We provide an original derivation of l 2 and l $\infty$ stability conditions of the scheme that do not deteriorate the second order accuracy without relying on a limiter of any type in the linear case, in particular for shock solutions, and extend such results to the nonlinear case, showing the novelty of the method. The prototype system for the linear case is the hyperbolic heat equation, whereas barotropic Euler equations of gas dynamics with friction are the one for the nonlinear case. The method is also able to yield very accurate steady solutions in the nonlinear case when the relaxation coefficient in the source term depends on space. A thorough numerical assessment of the proposed strategy is provided by investigating smooth solutions, solutions with shocks and solutions leading to a steady state with space dependent relaxation coefficient.

preprint2020arXiv

Travelling wave mathematical analysis and efficient numerical resolution for a one-dimensional model of solid propellant combustion

We investigate a model of solid propellant combustion involving surface pyrolysis coupled to finite activation energy gas phase combustion. Existence and uniqueness of a travelling wave solution are established by extending dynamical system tools classically used for premixed flames, dealing with the additional difficulty arising from the surface regression and pyrolysis. An efficient shooting method allows to solve the problem in phase space without resorting to space discretisation nor fixed-point Newton iterations. The results are compared to solutions from a CFD code developed at ONERA, assessing the efficiency and potential of the method, and the impact of the modelling assumptions is evaluated through parametric studies.

preprint2019arXiv

An asymptotic preserving well-balanced scheme for the isothermal fluid equations in low-temperature plasma applications

We present a novel numerical scheme for the efficient and accurate solution of the isothermal two-fluid (electron and ion) equations coupled to Poisson's equation for low-temperature plasmas. The model considers electrons and ions as separate fluids, comprising the electron inertia and charge separation. The discretization of this system with standard explicit schemes is constrained by very restrictive time steps and cell sizes related to the resolution of the Debye length, electron plasma frequency, and electron sound waves. Both sheath and electron inertia are fundamental to fully explain the physics in low-pressure and low-temperature plasmas. However, most of the phenomena of interest for fluid models occur at speeds much slower than the electron thermal speed and are quasi-neutral, except in small charged regions. In this work, we present a scheme based on the Lagrange-projection operator splitting that preserves the asymptotic regime where the plasma is quasi-neutral with massless electrons. As a result, the quasi-neutral regime is treated without the need of an implicit solver nor the resolution of the Debye length and electron plasma frequency. Additionally, the scheme proves to accurately represent the dynamics of the electrons both at low speeds and when the electron speed is comparable to the thermal speed. In addition, a well-balanced treatment of the ion source terms is proposed in order to tackle problems where the ion temperature is very low compared to the electron temperature. The scheme significantly improves the accuracy both in the quasi-neutral limit and in the presence of plasma sheaths when the Debye length is resolved. The scheme is assessed in two specifically designed test-cases: a quasi-neutral two-stream periodic perturbation with analytical solution and a low-temperature discharge that includes sheaths.

preprint2019arXiv

Entropy supplementary conservation law for non-linear systems of PDEs with non-conservative terms: application to the modelling and analysis of complex fluid flows using computer algebra

In the present contribution, we investigate first-order nonlinear systems of partial differential equations which are constituted of two parts: a system of conservation laws and non-conservative first order terms. Whereas the theory of first-order systems of conservation laws is well established and the conditions for the existence of supplementary conservation laws, and more specifically of an entropy supplementary conservation law for smooth solutions, well known, there exists so far no general extension to obtain such supplementary conservation laws when non-conservative terms are present. We propose a framework in order to extend the existing theory and show that the presence of non-conservative terms somewhat complexifies the problem since numerous combinations of the conservative and non-conservative terms can lead to a supplementary conservation law. We then identify a restricted framework in order to design and analyze physical models of complex fluid flows by means of computer algebra and thus obtain the entire ensemble of possible combination of conservative and non-conservative terms with the objective of obtaining specifically an entropy supplementary conservation law. The theory as well as developed computer algebra tool are then applied to a Baer-Nunziato two-phase flow model and to a multicomponent plasma fluid model. The first one is a first-order fluid model, with non-conservative terms impacting on the linearly degenerate field and requires a closure since there is no way to derive interfacial quantities from averaging principles and we need guidance in order to close the pressure and velocity of the interface and the thermodynamics of the mixture. The second one involves first order terms for the heavy species coupled to second order terms for the electrons, the non-conservative terms impact the genuinely nonlinear fields and the model can be rigorously derived from kinetic theory. We show how the theory allows to recover the whole spectrum of closures obtained so far in the literature for the two-phase flow system as well as conditions when one aims at extending the thermodynamics and also applies to the plasma case, where we recover the usual entropy supplementary equation, thus assessing the effectiveness and scope of the proposed theory.

preprint2018arXiv

Consistent transport properties in multicomponent two-temperature magnetized plasmas:

A fluid model is developed for multicomponent two-temperature magnetized plasmas in chemical non-equilibrium from the partially- to fully-ionized collisional regimes. We focus on transport phenomena aiming at representing the chromosphere of the Sun. Graille et al. [M3AS 19(04):527-599, 2009] have derived an asymptotic fluid model for multicomponent plamas from kinetic theory, yielding a rigorous description of the dissipative effects. The governing equations and consistent transport properties are obtained using a multiscale Chapman-Enskog perturbative solution to the Boltzmann equation based on a non-dimensional analysis. The mass disparity between the electrons and heavy particles is accounted for, as well as the influence of the electromagnetic field. We couple this model to the Maxwell equations for the electromagnetic field and derive the generalized Ohm's law for multicomponent plasmas. The model inherits a well-identified mathematical structure leading to an extended range of validity for the Sun chromosphere conditions. We compute consistent transport properties by means of a spectral Galerkin method using the Laguerre-Sonine polynomial approximation. Two non-vanishing polynomial terms are used when deriving the transport systems for electrons, whereas only one term is retained for heavy particles. In a simplified framework where the plasma is fully ionized, we compare the transport properties for the Sun chromosphere to conventional expressions for magnetized plasmas due to Braginskii, showing a good agreement between both results. For more general partially ionized conditions, representative of the Sun chromosphere, we compute the muticomponent transport properties corresponding to the species diffusion velocities, heavy-particle and electron heat fluxes, and viscous stress tensor of the model, for a Helium-Hydrogen mixture in local thermodynamic equilibrium. The model is assessed for the 3D radiative magnetohydrodynamic simulation of a pore, in the highly turbulent upper layer of the solar convective zone. The resistive term is found to dominate mainly the dynamics of the electric field at the pore location. The battery term for heavy particles appears to be higher at the pore location and at some intergranulation boundaries.

preprint2011arXiv

A new numerical strategy with space-time adaptivity and error control for multi-scale streamer discharge simulations

This paper presents a new resolution strategy for multi-scale streamer discharge simulations based on a second order time adaptive integration and space adaptive multiresolution. A classical fluid model is used to describe plasma discharges, considering drift-diffusion equations and the computation of electric field. The proposed numerical method provides a time-space accuracy control of the solution, and thus, an effective accurate resolution independent of the fastest physical time scale. An important improvement of the computational efficiency is achieved whenever the required time steps go beyond standard stability constraints associated with mesh size or source time scales for the resolution of the drift-diffusion equations, whereas the stability constraint related to the dielectric relaxation time scale is respected but with a second order precision. Numerical illustrations show that the strategy can be efficiently applied to simulate the propagation of highly nonlinear ionizing waves as streamer discharges, as well as highly multi-scale nanosecond repetitively pulsed discharges, describing consistently a broad spectrum of space and time scales as well as different physical scenarios for consecutive discharge/post-discharge phases, out of reach of standard non-adaptive methods.

preprint2011arXiv

Adaptive time splitting method for multi-scale evolutionary partial differential equations

This paper introduces an adaptive time splitting technique for the solution of stiff evolutionary PDEs that guarantees an effective error control of the simulation, independent of the fastest physical time scale for highly unsteady problems. The strategy considers a second order Strang method and another lower order embedded splitting scheme that takes into account potential loss of order due to the stiffness featured by time-space multi-scale phenomena. The scheme is then built upon a precise numerical analysis of the method and a complementary numerical procedure, conceived to overcome classical restrictions of adaptive time stepping schemes based on lower order embedded methods, whenever asymptotic estimates fail to predict the dynamics of the problem. The performance of the method in terms of control of integration errors is evaluated by numerical simulations of stiff propagating waves coming from nonlinear chemical dynamics models as well as highly multi-scale nanosecond repetitively pulsed gas discharges, which allow to illustrate the method capabilities to consistently describe a broad spectrum of time scales and different physical scenarios for consecutive discharge/post-discharge phases.

preprint2007arXiv

Kinetic Theory of Plasmas: Translational Energy

In the present contribution, we derive from kinetic theory a unified fluid model for multicomponent plasmas by accounting for the electromagnetic field influence. We deal with a possible thermal nonequilibrium of the translational energy of the particles, neglecting their internal energy and the reactive collisions. Given the strong disparity of mass between the electrons and heavy particles, such as molecules, atoms, and ions, we conduct a dimensional analysis of the Boltzmann equation. We then generalize the Chapman-Enskog method, emphasizing the role of a multiscale perturbation parameter on the collisional operator, the streaming operator, and the collisional invariants of the Boltzmann equation. The system is examined at successive orders of approximation, each of which corresponding to a physical time scale. The multicomponent Navier-Stokes regime is reached for the heavy particles, which follow a hyperbolic scaling, and is coupled to first order drift-diffusion equations for the electrons, which follow a parabolic scaling. The transport coefficients exhibit an anisotropic behavior when the magnetic field is strong enough. We also give a complete description of the Kolesnikov effect, i.e., the crossed contributions to the mass and energy transport fluxes coupling the electrons and heavy particles. Finally, the first and second principles of thermodynamics are proved to be satisfied by deriving a total energy equation and an entropy equation. Moreover, the system of equations is shown to be conservative and the purely convective system hyperbolic, thus leading to a well-defined structure.

preprint2007arXiv

Numerical simulation of spray coalescence in an eulerian framework : direct quadrature method of moments and multi-fluid method

The scope of the present study is Eulerian modeling and simulation of polydisperse liquid sprays undergoing droplet coalescence and evaporation. The fundamental mathematical description is the Williams spray equation governing the joint number density function f(v, u; x, t) of droplet volume and velocity. Eulerian multi-fluid models have already been rigorously derived from this equation in Laurent et al. (2004). The first key feature of the paper is the application of direct quadrature method of moments (DQMOM) introduced by Marchisio and Fox (2005) to the Williams spray equation. Both the multi-fluid method and DQMOM yield systems of Eulerian conservation equations with complicated interaction terms representing coalescence. In order to validate and compare these approaches, the chosen configuration is a self-similar 2D axisymmetrical decelerating nozzle with sprays having various size distributions, ranging from smooth ones up to Dirac delta functions. The second key feature of the paper is a thorough comparison of the two approaches for various test-cases to a reference solution obtained through a classical stochastic Lagrangian solver. Both Eulerian models prove to describe adequately spray coalescence and yield a very interesting alternative to the Lagrangian solver.