## Section:
Scientific Foundations2>
### Numerical schemes for fluid mechanics3>

Participants : Rémi Abgrall, Mario Ricchiuto, Pietro Marco Congedo, Cécile Dobrzynzki, Héloïse Beaugendre, Pierre-Henri Maire, Luc Mieussens.

A large number of engineering problems involve fluid
mechanics. They may involve the coupling of one or more physical
models. An example is provided by aeroelastic problems, which have
been studied in details by other Inria teams. Another example is given by
flows in pipelines where the fluid (a mixture of air–water–gas) does
not have
well-known physical properties, and there are even more exotic situations
that will be discussed later.
Another application is the influence of fluid flow on noise production.
Problems in aeroacoustics are indeed becoming more and more important in
everyday life. In some occasions, one needs specific numerical tools
to take into account *e.g.* a fluids' exotic equation of state, or
because the
amount of required computational resources becomes huge, as in unsteady flows. Another situation
where specific tools are needed is when one is interested in very specific physical quantities, such as
*e.g.* the lift and drag of an airfoil, a situation where commercial tools can only provide a very crude answer.

It is a fact that there are many commercial codes. They allow users to
simulate a lot of different flow types. The quality of the results is however
far from optimal in many cases. Moreover, the numerical technology implemented in these codes
is often not the most recent. To give a few examples, consider the noise generated
by wake vortices in supersonic flows (external aerodynamics/aeroacoustics),
or the direct simulation of a 3D compressible mixing layer in a complex geometry (as in combustion chambers).
Up to our knowledge, due to the very different temporal
and physical scales that need to be captured,
a direct simulation of these phenomena is
not in the reach of the most recent technologies because the numerical
resources required are currently unavailable.
*We need to invent specific algorithms for this purpose.*

In order to simulate efficiently these complex physical
problems, we are working on some fundamental aspects of the numerical
analysis of non linear hyperbolic problems. *Our goal is to develop
more accurate and more efficient schemes that can adapt to modern computer architectures*.

More precisely, *we are working on a class of numerical schemes*, known in literature as
Residual Distribution schemes, *specifically
tailored to unstructured and hybrid meshes*. They have the most possible compact
stencil that is compatible with the expected order of accuracy.
This *accuracy is at least of second order, and it can go up to any order
of accuracy, even though fourth order is considered for practical applications.*
Since the stencil is compact, the implementation on parallel machines becomes
simple. These schemes are very flexible in nature, which is so far one of the most important advantage
over other techniques. This feature has allowed us to adapt the schemes to the requirements of different
physical situations (*e.g.* different formulations allow either en efficient explicit
time advancement for problems involving small time-scales, or a fully implicit space-time
variant which is unconditionally stable and allows to handle stiff problems
where only the large time scales are relevant). This flexibility has also enabled
to devise a variant using the same data structure of the popular Discontinuous Galerkin
schemes, which are also part of our scientific focus.

The compactness of the second order version of the schemes enables us to use efficiently the high performance parallel linear algebra tools developed by the team. However, the high order versions of these schemes, which are under development, require modifications to these tools taking into account the nature of the data structure used to reach higher orders of accuracy. This leads to new scientific problems at the border between numerical analysis and computer science. In parallel to these fundamental aspects, we also work on adapting more classical numerical tools to complex physical problems such as those encountered in interface flows, turbulent or multiphase flows, geophysical flows, and material science. An effort for developing a more predictive tool for multiphase compressible flows is also underway. Within this project, several advancements have been performed, i.e. considering a more complete systems of equations including viscosity, working on the thermodynamic modeling of complex fluids, and developing stochastic methods for uncertainty quantification in compressible flows.

We expect within a few years to be able to demonstrate the potential of our developments on applications ranging from the the reproduction of the complex multidimensional interactions between tidal waves and estuaries, unsteady aerodynamics and aeroacoustics associated to both external and internal compressible flows, compressible ideal and non-ideal MHD (in relation with the ITER project), and the behavior of complex materials. This will be achieved by means of a multi-disciplinary effort involving our research on residual discretizations schemes, the parallel advances in algebraic solvers and partitioners, and the strong interactions with specialists in computer science, scientific computing, physics, mechanics, and mathematical modeling.

Our research in numerical algorithms has led to the development of the
`RealfluiDS` platform which is described in
section
5.3 .
New software developments are under way in the field of free surface flows and complex materials modeling.
These developments are performed in the code `SLOWS` (Shallow-water fLOWS) for free surface flows,
and in the solver `COCA` (CodeOxydationCompositesAutocicatrisants)
for the simulation of the self-healing process in composite materials.
These developments will be described in sections
5.10 and
5.2 .

This work is supported by the EU-Strep IDIHOM, various research contracts and in
part by the `ANEMOS` project and the ANR-Emergence `RealFluids`
grant. A large part of the team also beneficiates of the `ADDECCO` ERC grant.