Source author record

Florent Renac

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

ResearcherUnclaimed source record

math.NA Numerical Analysis physics.flu-dyn

Catalog footprint

What is connected

5works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Invariant domain preserving high-order spectral discontinuous approximations of hyperbolic systems

We propose a limiting procedure to preserve invariant domains with time explicit discrete high-order spectral discontinuous approximate solutions to hyperbolic systems of conservation laws. Provided the scheme is discretely conservative and satisfy geometric conservation laws at the discrete level, we derive a condition on the time step to guaranty that the cell-averaged approximate solution is a convex combination of states in the invariant domain. These states are then used to define local bounds which are then imposed to the full high-order approximate solution within the cell via an a posteriori scaling limiter. Numerical experiments are then presented with modal and nodal discontinuous Galerkin schemes confirm the robustness and stability enhancement of the present approach.

preprint2021arXiv

Energy relaxation approximation for the compressible multicomponent flows in thermal nonequilibrium

This work concerns the numerical approximation with a finite volume method of inviscid, nonequilibrium, high-temperature flows in multiple space dimensions. It is devoted to the analysis of the numerical scheme for the approximation of the hyperbolic system in homogeneous form. We derive a general framework for the design of numerical schemes for this model from numerical schemes for the monocomponent compressible Euler equations for a polytropic gas. Under a very simple condition on the adiabatic exponent of the polytropic gas, the scheme for the multicomponent system enjoys the same properties as the one for the monocomponent system: discrete entropy inequality, positivity of the partial densities and internal energies, discrete maximum principle on the mass fractions, and discrete minimum principle on the entropy. Our approach extends the relaxation of energy [Coquel and Perthame, \textit{SIAM J. Numer. Anal.}, 35 (1998), 2223--2249] to the multicomponent Euler system. In the limit of instantaneous relaxation we show that the solution formally converges to a unique and stable equilibrium solution to the multicomponent Euler equations. We then use this framework to design numerical schemes from three schemes for the polytropic Euler system: the Godunov exact Riemann solver [Godunov, Math. Sbornik, 47 (1959), 271--306] and the HLL [Harten et al., SIAM Rev., 25 (1983), 35--61] and pressure relaxation based [Bouchut, Nonlinear stability of finite volume methods for hyperbolic conservation laws and well-balanced schemes for sources, Frontiers in Mathematics, Birkhäuser, 2004] approximate Riemann solvers. Numerical experiments in one and two space dimensions on flows with discontinuous solutions support the conclusions of our analysis and highlight stability, robustness and convergence of the scheme.

preprint2016arXiv

A robust high-order Lagrange-projection like scheme with large time steps for the isentropic Euler equations

We present an extension to high-order of a first-order Lagrange-projection like method for the approximation of the Euler equations introduced in Coquel {\it et al.} (Math. Comput., 79 (2010), pp.~1493--1533). The method is based on a decomposition between acoustic and transport operators associated to an implicit-explicit time integration, thus relaxing the constraint of acoustic waves on the time step. We propose here to use a discontinuous Galerkin method for the space approximation. Considering the isentropic Euler equations, we derive conditions to keep positivity of the mean value of density and satisfy an entropy inequality for the numerical solution in each element of the mesh at any approximation order in space. These results allow to design limiting procedures to restore these properties at nodal values within elements. Numerical experiments support the conclusions of the analysis and highlight stability and robustness of the present method, though it allows the use of large time steps.

preprint2015arXiv

Stationary discrete shock profiles for scalar conservation laws with a discontinuous Galerkin method

We present an analysis of stationary discrete shock profiles for a discontinuous Galerkin method approximating scalar nonlinear hyperbolic conservation laws with a convex flux. Using the Godunov method for the numerical flux, we characterize the steady state solutions for arbitrary approximation orders and show that they are oscillatory only in one mesh cell and are parametrized by the shock strength and its relative position in the cell. In the particular case of the inviscid Burgers equation, we derive analytical solutions of the numerical scheme and predict their oscillations up to fourth-order of accuracy. Moreover, a linear stability analysis shows that these profiles may become unstable at points where the Godunov flux is not differentiable. Theoretical and numerical investigations show that these results can be extended to other numerical fluxes. In particular, shock profiles are found to vanish exponentially fast from the shock position for some class of monotone numerical fluxes and the oscillatory and unstable characters of their solutions present strong similarities with that of the Godunov method.

preprint2014arXiv

Computation of eigenvalue sensitivity to base flow modifications in a discrete framework: Application to open-loop control

A fully discrete formalism is introduced to perform stability analysis of a turbulent compressible flow whom dynamics is modeled with the Reynolds-Averaged Navier-Stokes (RANS) equations. The discrete equations are linearized using finite differences and the Jacobian is computed using repeated evaluation of the residuals. Stability of the flow is assessed solving an eigenvalue problem. The sensitivity gradients which indicate regions of the flow where a passive control device could stabilize the unstable eigenvalues are defined within this fully discrete framework. Second order finite differences are applied to the discrete residual to compute the gradients. In particular, the sensitivity gradients are shown to be linked to the Hessian of the RANS equations. The introduced formalism and linearization method are generic: the code used to evaluate the residual of the RANS equations can be used in a black box manner, and the complex linearization of the Hessian is avoided. The method is tested on a two dimensional deep cavity case, the flow is turbulent with a Reynolds number equal to 860 000 and compressible with a Mach number of 0.8. Several turbulence models and numerical schemes are used to validate the method. Physical features of the flow are recovered, such as the fundamental frequency of the natural flow as well as acoustic mechanisms, suggesting the validity of the method. The sensitivity gradients are then computed and validated, the error in predicting the eigenvalue variation being found less than 3\%. Control maps using a small steady control device are finally obtained, indicating that the control area should be chosen in the vicinity of the leading edge of the cavity.