BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Rémi Abgrall

Rémi Abgrall contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phphysics.flu-dyn

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)

preprint2023arXiv

"A Posteriori" Limited High Order and Robust Residual Distribution Schemes for Transient Simulations of Fluid Flows in Gas Dynamics

In this paper, we propose a novel approximation strategy for time-dependent hyperbolic systems of conservation laws for the Euler system of gas dynamics that aims to represent the dynamics of strong interacting discontinuities. The goal of our method is to allow an approximation with a high-order of accuracy in smooth regions of the flow, while ensuring robustness and a non-oscillatory behaviour in the regions of steep gradients, in particular across shocks. Following the Multidimensional Optimal Order Detection (MOOD) (Clain et al., Journal of Computational Physics 2011; Diot et al., Computer & Fluids 2012) approach, a candidate solution is computed at a next time level via a high-order accurate explicit scheme (Abgrall, Journal of Scientific Computing 2017; Abgrall et al., Computers & Mathematics with Applications 2018). A so-called detector determines if the candidate solution reveals any spurious oscillations or numerical issue and, if so, only the troubled cells are locally recomputed via a more dissipative scheme. This allows to design a family of "a posteriori" limited, robust and positivity preserving, as well as high accurate, non-oscillatory and effective scheme. Among the detecting criteria of the novel MOOD strategy, two different approaches from literature, based on the work of Clain et al. (Journal of Computational Physics 2011), Diot et al. (Computer & Fluids 2012) and of Vilar (Journal of Computational Physics 2019), are investigated. Numerical examples in 1D and 2D, on structured and unstructured meshes, are proposed to assess the effective order of accuracy for smooth flows, the non-oscillatory behaviour on shocked flows, the robustness and positivity preservation on more extreme flows.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

A hybrid finite element - finite volume method for conservation laws

We propose an arbitrarily high-order accurate numerical method for conservation laws that is based on a continuous approximation of the solution. The degrees of freedom are point values at cell interfaces and moments of the solution inside the cell. To lowest ($3^\text{rd}$) order this method reduces to the Active Flux method. The update of the moments is achieved immediately by integrating the conservation law over the cell, integrating by parts and employing the continuity across cell interfaces. We propose two ways how the point values can be updated in time: either by first deriving a semi-discrete method that uses a finite-difference-type formula to approximate the spatial derivative, and integrating this method e.g. with a Runge-Kutta scheme, or by using a characteristics-based update, which is inspired by the original (fully discrete) Active Flux method. We analyze stability and accuracy of the resulting methods.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Extensions of Active Flux to arbitrary order of accuracy

Active Flux is a recently developed numerical method for hyperbolic conservation laws. Its classical degrees of freedom are cell averages and point values at cell interfaces. These latter are shared between adjacent cells, leading to a globally continuous reconstruction. The update of the point values includes upwinding, but without solving a Riemann Problem. The update of the cell average requires a flux quadrature at the cell interface, which can be immediately performed using the point values. This paper explores different extensions of Active Flux to arbitrarily high order of accuracy, while maintaining the idea of global continuity. We propose to either increase the stencil while keeping the same degrees of freedom, or to increase the number of point values, or to include higher moments as new degrees of freedom. These extensions have different properties, and reflect different views upon the relation of Active Flux to the families of Finite Volume, Finite Difference and Finite Element methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A pressure-based method for weakly compressible two-phase flows under a Baer-Nunziato type model with generic equations of state and pressure and velocity disequilibrium

Within the framework of diffuse interface methods, we derive a pressure-based Baer-Nunziato type model well-suited to weakly compressible multiphase flows. The model can easily deal with different equation of states and it includes relaxation terms characterized by user-defined finite parameters, which drive the pressure and velocity of each phase toward the equilibrium. There is no clear notion of speed of sound, and thus, most of the classical low Mach approximation cannot easily be cast in this context. The proposed solution strategy consists of two operators: a semi-implicit finite-volume solver for the hyperbolic part and an ODE integrator for the relaxation processes. Being the acoustic terms in the hyperbolic part integrated implicitly, the stability condition on the time step is lessened. The discretization of non-conservative terms involving the gradient of the volume fraction fulfills by construction the non-disturbance condition on pressure and velocity to avoid oscillations across the multimaterial interfaces. The developed simulation tool is validated through one-dimensional simulations of shock-tube and Riemann-problems, involving water-aluminum and water-air mixtures, vapor-liquid mixture of water and of carbon dioxide, and almost pure flows. The numerical results match analytical and reference ones, except some expected discrepancies across shocks, which however remain acceptable (errors within some percentage points). All tests were performed with acoustic CFL numbers greater than one, and no stability issues arose, even for CFL greater than 10. The effects of different values of relaxation parameters and of different amount equations of state -- stiffened gas and Peng-Robinson -- were investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reinterpretation and Extension of Entropy Correction Terms for Residual Distribution and Discontinuous Galerkin Schemes: Application to Structure Preserving Discretization

For the general class of residual distribution (RD) schemes, including many finite element (such as continuous/discontinuous Galerkin) and flux reconstruction methods, an approach to construct entropy conservative/ dissipative semidiscretizations by adding suitable correction terms has been proposed by Abgrall (J.~Comp.~Phys. 372: pp. 640--666, 2018). In this work, the correction terms are characterized as solutions of certain optimization problems and are adapted to the SBP-SAT framework, focusing on discontinuous Galerkin methods. Novel generalizations to entropy inequalities, multiple constraints, and kinetic energy preservation for the Euler equations are developed and tested in numerical experiments. For all of these optimization problems, explicit solutions are provided. Additionally, the correction approach is applied for the first time to obtain a fully discrete entropy conservative/dissipative RD scheme. Here, the application of the deferred correction (DeC) method for the time integration is essential. This paper can be seen as describing a systematic method to construct structure preserving discretization, at least for the considered example.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Assessment of a non-conservative four-equation multiphase system with phase transition

This work focuses on the formulation of a four-equation model for simulating unsteady two-phase mixtures with phase transition and strong discontinuities. The main assumption consists in a homogeneous temperature, pressure and velocity fields between the two phases. Specifically, we present the extension of a residual distribution scheme to solve a four-equation two-phase system with phase transition written in a non-conservative form, i.e. in terms of internal energy instead of the classical total energy approach. This non-conservative formulation allows avoiding the classical oscillations obtained by many approaches, that might appear for the pressure profile across contact discontinuities. The proposed method relies on a Finite Element based Residual Distribution scheme which is designed for an explicit second-order time stepping. We test the non-conservative Residual Distribution scheme on several benchmark problems and assess the results via a cross-validation with the approximated solution obtained via a conservative approach, based on a HLLC scheme. Furthermore, we check both methods for mesh convergence and show the effective robustness on very severe test cases, that involve both problems with and without phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High Order Residual Distribution Conservative Finite Difference HWENO Scheme for Steady State Problems

In this paper, we develop a high order residual distribution (RD) method for solving steady state conservation laws in a novel Hermite weighted essentially non-oscillatory (HWENO) framework recently developed in [24]. In particular, we design a high order HWENO integration for the integrals of source term and fluxes based on the point value of the solution and its spatial derivatives, and the principles of residual distribution schemes are adapted to obtain steady state solutions. Two advantages of the novel HWENO framework have been shown in [24]: first, compared with the traditional HWENO framework, the proposed method does not need to introduce additional auxiliary equations to update the derivatives of the unknown variable, and just compute them from the current point value of the solution and its old spatial derivatives, which saves the computational storage and CPU time, and thereby improve the computational efficiency of the traditional HWENO framework. Second, compared with the traditional WENO method, reconstruction stencil of the HWENO methods becomes more compact, their boundary treatment is simpler, and the numerical errors are smaller at the same grid. Thus, it is also a compact scheme when we design the higher order accuracy, compared with that in [11] Chou and Shu proposed. Extensive numerical experiments for one- and two-dimensional scalar and systems problems confirm the high order accuracy and good quality of our scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

General polytopal H(div) conformal finite elements and their discretisation spaces

We present a class of discretisation spaces and H(div)-conformal elements that can be built on any polytope. Bridging the flexibility of the Virtual Element spaces towards the element's shape with the divergence properties of the Raviart-Thomas elements on the boundaries, the designed frameworks offer a wide range of H(div)-conformal discretisations. As those elements are set up through degrees of freedom, their definitions are easily amenable to the properties the approximated quantities are wished to fulfill. Furthermore, we show that one straightforward restriction of this general setting share its properties with the classical Raviart-Thomas elements at each interface, for any order and any polytopial shape. Then, to close the introduction of those new elements by an example, we investigate the shape of the basis functions corresponding to particular elements in the two dimensional case.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Asymptotic preserving Deferred Correction Residual Distribution schemes

This work aims to extend the residual distribution (RD) framework to stiff relaxation problems. The RD is a class of schemes which is used to solve hyperbolic system of partial differential equations. Up to our knowledge, it was used only for systems with mild source terms, such as gravitation problems or shallow water equations. What we propose is an IMEX (implicit--explicit) version of the residual distribution schemes, that can resolve stiff source terms, without refining the discretization up to the stiffness scale. This can be particularly useful in various models, where the stiffness is given by topological or physical quantities, e.g. multiphase flows, kinetic models, viscoelasticity problems. Moreover, the provided scheme is able to catch different relaxation scales automatically, without losing accuracy. The scheme is asymptotic preserving and this guarantees that in the relaxation limit, we recast the expected macroscopic behaviour. To get a high order accuracy, we use an IMEX time discretization combined with a Deferred Correction (DeC) procedure, while naturally RD provides high order space discretization. Finally, we show some numerical tests in 1D and 2D for stiff systems of equations.

0 citations0 reviewsNo code linkedOpen access