Section: Research Program
Modelling and analysis
Hazardous flows are complex physical phenomena that can hardly be represented by shallow water type systems of partial differential equations (PDEs). In this domain, the research program is devoted to the derivation and analysis of reduced complexity models compared to the NavierStokes equations, but relaxing the shallow water assumptions. The main purpose is then to obtain models welladapted to the physical phenomena at stake.
Even if the resulting models do not strictly belong to the family of hyperbolic systems, they exhibit hyperbolic features: the analysis and discretization techniques we intend to develop have connections with those used for hyperbolic conservation laws. It is worth noticing that the need for robust and efficient numerical procedures is reinforced by the smallness of dissipative effects in geophysical models which therefore generate singular solutions and instabilities.
On the one hand, the derivation of the SaintVenant system from the NavierStokes equations is based on two approximations, socalled shallow water assumptions, namely

the horizontal fluid velocity is well approximated by its mean value along the vertical direction,

the pressure is hydrostatic or equivalently the vertical acceleration of the fluid can be neglected compared to the gravitational effects.
As a consequence the objective is to get rid of these two assumptions, one after the other, in order to obtain models accurately approximating the incompressible Euler or NavierStokes equations.
On the other hand, many applications require the coupling with nonhydrodynamic equations, as in the case of microalgae production or erosion processes. These new equations comprise nonhyperbolic features and must rely on a special analysis.
Multilayer approach
As for the first shallow water assumption, multilayer systems were proposed describing the flow as a superposition of SaintVenant type systems [21] , [24] , [25] . Even if this approach has provided interesting results, layers are considered separate and nonmiscible fluids, which imply strong limitation. That is why we proposed a slightly different approach [22] , [23] based on Galerkin type decomposition along the vertical axis of all variables and leading, both for the model and its discretization, to more accurate results.
A kinetic representation of our multilayer model allows to derive robust numerical schemes endowed with properties such as: consistency, conservativity, positivity, preservation of equilibria,... It is one of the major achievements of the team but it needs to be analyzed and extended in several directions namely:

The convergence of the multilayer system towards the hydrostatic Euler system as the number of layers goes to infinity is a critical point. It is not fully satisfactory to have only formal estimates of the convergence and sharp estimates would enable to guess the optimal number of layers.

The introduction of several source terms due for instance to Coriolis forces or extra terms from changes of coordinates seems necessary. Their inclusion should lead to substantial modifications of the numerical scheme.

Its hyperbolicity has not yet been proved and conversely the possible loss of hyperbolicity cannot be characterized. Similarly, the hyperbolic feature is essential in the propagation and generation of waves.
Nonhydrostatic models
The hydrostatic assumption consists in neglecting the vertical acceleration of the fluid. It is considered valid for a large class of geophysical flows but is restrictive in various situations where the dispersive effects (like wave propagation) cannot be neglected. For instance, when a wave reaches the coast, bathymetry variations give a vertical acceleration to the fluid that strongly modifies the wave characteristics and especially its height.
When processing an asymptotic expansion (w.r.t. the aspect ratio for shallow water flows) into the NavierStokes equations, we obtain at the leading order the SaintVenant system. Going one step further leads to a vertically averaged version of the Euler/NavierStokes equations integrating the nonhydrostatic terms. This model has several advantages:

it admits an energy balance law (that is not the case for most dispersive models available in the literature),

it reduces to the SaintVenant system when the nonhydrostatic pressure term vanishes,

it consists in a set of conservation laws with source terms,
Multiphysics modelling
The coupling of hydrodynamic equations with other equations in order to model interactions between complex systems represents an important part of the team research. More precisely, three multiphysic systems are investigated. More details about the industrial impact of these studies are presented in the following section.

To estimate the risk for infrastructures in coastal zone or close to a river, the resolution of the shallow water equations with moving bathymetry is necessary. The first step consisted in the study of an equation largely used in engineering science: The Exner equation. The analysis enabled to exhibit drawbacks of the coupled model such as the lack of energy conservation or the strong variations of the solution from small perturbations. A new formulation is proposed to avoid these drawbacks. The new model consists in a coupling between conservation laws and an elliptic equation, like the system Euler/Poisson, suggesting to use wellknown strategies for the analysis and the numerical resolution. In addition, the new formulation is derived from classical complex rheology models and allowed physical phenomena such as threshold laws.

Interaction between flows and floating structures is the challenge at the scale of the shallow water equations. This study needs a better understanding of the energy exchanges between the flow and the structure. The mathematical model of floating structures is very hard to solve numerically due to the nonpenetration condition at the interface between the flow and the structure. It leads to infinite potential wave speeds that could not be solved with classical free surface numerical scheme. A relaxation model was derived to overcome this difficulty. It represents the interaction with the floating structure with a free surface modeltype.

If the interactions between hydrodynamics and biology phenomena are known through laboratory experiments, it is more difficult to predict the evolution, especially for the biological quantities, in a real and heterogeneous system. The objective is to model and reproduce the hydrodynamics modifications due to forcing term variations (in time and space). We are typically interested in phenomena such as eutrophication, development of harmful bacteria (cyanobacteria) and upwelling phenomena.