Among all aspects of geosciences, we mainly focus on gravity driven flows arising in many situations such as

There exists a strong demand from scientists and engineers in fluid mechanics for models and numerical tools able to simulate not only the water depth and the velocity field but also the distribution and evolution of external quantities such as pollutants or biological species and the interaction between flows and structures (seashores, erosion processes...). The key point of the researches carried out within ANGE is to answer this demand by the development of efficient, robust and validated models and numerical tools.

Due to the variety of applications with a wide range of spatial scales, reduced-size models like the shallow water equations are generally required. From the modelling point of view, the main issue is to describe the behaviour of the flow with a reduced-size model taking into account several physical processes such as non-hydrostatic terms, biological species evolution, topography and structure interactions within the flow. The mathematical analysis of the resulting model do not enter the field of hyperbolic equations anymore and new strategies have to be proposed. Moreover, efficient numerical resolutions of reduced-size models require particular attention due to the different time scales of the processes and in order to recover physical properties such as positivity, conservativity, entropy dissipation and equilibria.

The models can remain subject to uncertainties that originate from incomplete description of the physical processes and from uncertain parameters. Further development of the models may rely on the assimilation of observational data and the uncertainty quantification of the resulting analyses or forecasts.

The research activities carried out within the ANGE team strongly couple the development of methodological tools with applications to real–life problems and the transfer of numerical codes. The main purpose is to obtain new models adapted to the physical phenomena at stake, identify the main properties that reflect the physical meaning of the models (uniqueness, conservativity, entropy dissipation, ...), propose effective numerical methods to approximate their solution in complex configurations (multi-dimensional, unstructured meshes, well-balanced, ...) and to assess the results with data in the purpose of potentially correcting the models.

The difficulties arising in gravity driven flow studies are threefold.

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 Navier-Stokes equations, but relaxing the shallow water assumptions. The main purpose is then to obtain models well-adapted 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 discretisation 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 Saint-Venant system from the Navier-Stokes equations is based on two approximations (the so-called shallow water assumptions), namely

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 Navier-Stokes equations.

On the other hand, many applications require the coupling with non-hydrodynamic equations, as in the case of micro-algae production or erosion processes. These new equations comprise non-hyperbolic features and a special analysis is needed.

As for the first shallow water assumption, multi-layer systems
were proposed to describe the flow as a superposition of Saint-Venant type systems 28, 31, 32.
Even if this approach has provided interesting results, layers are considered separate and non-miscible fluids, which implies strong limitations.
That is why we proposed a slightly different approach 29, 30 based on a Galerkin type decomposition along the vertical axis of all variables and leading, both for the model and its discretisation, to more accurate results.

A kinetic representation of our multilayer model allows to derive robust numerical schemes endowed with crucial 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 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.

Processing an asymptotic expansion (w.r.t. the aspect ratio for shallow water flows) into the Navier-Stokes equations, we obtain at the leading order the Saint-Venant system. Going one step further leads to a vertically averaged version of the Euler/Navier-Stokes equations involving some non-hydrostatic terms. This model has several advantages:

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 multi-physics systems are investigated. More details about the industrial impact of these studies are presented in the following section.

In environmental applications, the most accurate numerical models remain subject to uncertainties that originate from their parameters and shortcomings in their physical formulations. It is often desirable to quantify the resulting uncertainties in a model forecast. The propagation of the uncertainties may require the generation of ensembles of simulations that ideally sample from the probability density function of the forecast variables. Classical approaches rely on multiple models and on Monte Carlo simulations. The applied perturbations need to be calibrated for the ensemble of simulations to properly sample the uncertainties. Calibrations involve ensemble scores that compare the consistency between the ensemble simulations and the observational data. The computational requirements are so high that designing fast surrogate models or metamodels is often required.

In order to reduce the uncertainties, the fixed or mobile observations of various origins and accuracies can be merged with the simulation results. The uncertainties in the observations and their representativeness also need to be quantified in the process. The assimilation strategy can be formulated in terms of state estimation or parameter estimation (also called inverse modelling). Different algorithms are employed for static and dynamic models, for analyses and forecasts. A challenging question lies in the optimization of the observational network for the assimilation to be the most efficient at a given observational cost.

The main challenge in the study of the non-hydrostatic model is to design a robust and efficient numerical scheme endowed with properties such as: positivity, wet/dry interfaces treatment, consistency. It must be noticed that even if the non-hydrostatic model looks like an extension of the Saint-Venant system, most of the known techniques used in the hydrostatic case are not efficient as we recover strong difficulties encountered in incompressible fluid mechanics due to the extra pressure term. These difficulties are reinforced by the absence of viscous/dissipative terms.

In the quest for a better balance between accuracy and efficiency,
a strategy consists in the adaptation of models.
Indeed, the systems of partial differential equations we consider result from a hierarchy of simplifying
assumptions. However, some of these hypotheses may turn out to be irrelevant locally. The adaptation of
models thus consists in determining areas where a simplified model
(e.g. shallow water type) is valid and
where it is not. In the latter case, we may go back to the “parent” model
(e.g. Euler) in the corresponding area.
This implies to know how to handle the coupling between the aforementioned models from both theoretical
and numerical points of view. In particular, the numerical treatment of transmission conditions is a key point. It requires the estimation of characteristic values (Riemann invariant) which have to be determined according to the regime (torrential or fluvial).

Hydrodynamic models comprise advection and sources terms. The conservation of the balance between source terms, typically viscosity and friction, has a significant impact since the overall flow is generally a perturbation around an equilibrium. The design of numerical schemes able to preserve such balances is a challenge from both theoretical and industrial points of view. The concept of Asymptotic-Preserving (AP) methods is of great interest in order to overcome these issues.

Another difficulty occurs when a term, typically related to the pressure,
becomes very large compared to the order of magnitude of the velocity.
At this regime, namely the so-called low Froude (shallow water)
or low Mach (Euler) regimes, the difference between the speed of the gravity waves
and the physical velocity makes classical numerical schemes inefficient:
firstly because of the error of truncation which is inversely proportional to the small parameters,
secondly because of the time step governed by the largest speed of the gravity wave.
AP methods made a breakthrough in the numerical resolution
of asymptotic perturbations of partial-differential equations concerning the first point.
The second one can be fixed using partially implicit scheme.

Coupling problems also arise within the fluid when it contains pollutants, density variations or biological species. For most situations, the interactions are small enough to use a splitting strategy and the classical numerical scheme for each sub-model, whether it be hydrodynamic or non-hydrodynamic.

The sediment transport raises interesting issues from a numerical aspect. This is an example of coupling between the flow and another phenomenon, namely the deformation of the bottom of the basin that can be carried out either by bed load where the sediment has its own velocity or suspended load in which the particles are mostly driven by the flow. This phenomenon involves different time scales and nonlinear retroactions; hence the need for accurate mechanical models and very robust numerical methods. In collaboration with industrial partners (EDF–LNHE), the team already works on the improvement of numerical methods for existing (mostly empirical) models but our aim is also to propose new (quite) simple models that contain important features and satisfy some basic mechanical requirements. The extension of our 3D models to the transport of weighted particles can also be here of great interest.

Numerical simulations are a very useful tool for the design of new processes, for instance in renewable energy or water decontamination. The optimisation of the process according to a well-defined objective such as the production of energy or the evaluation of a pollutant concentration is the logical upcoming challenge in order to propose competitive solutions in industrial context. First of all, the set of parameters that have a significant impact on the result and on which we can act in practice is identified. Then the optimal parameters can be obtained using the numerical codes produced by the team to estimate the performance for a given set of parameters with an additional loop such as gradient descent or Monte Carlo method. The optimisation is used in practice to determine the best profile for turbine pales, the best location for water turbine implantation, in particular for a farm.

Sustainable development and environment preservation have a growing importance and scientists have to address difficult issues such as: management of water resources, renewable energy production, bio/geo-chemistry of oceans, resilience of society w.r.t. hazardous flows, urban pollutions, ...

As mentioned above, the main issue is to propose models of reduced complexity, suitable for scientific computing and endowed with stability properties (continuous and/or discrete). In addition, models and their numerical approximations have to be confronted with experimental data, as analytical solutions are hardly accessible for these problems/models. A. Mangeney (IPGP) and N. Goutal (EDF) may provide useful data.

Reduced models like the shallow water equations are particularly well-adapted to the modelling of geophysical flows since there are characterized by large time or/and space scales. For long time simulations, the preservation of equilibria is essential as global solutions are a perturbation around them. The analysis and the numerical preservation of non-trivial equilibria, more precisely when the velocity does not vanish, are still a challenge. In the fields of oceanography and meteorology, the numerical preservation of the so-called geostrophic state, which is the balance between the gravity field and the Coriolis force, can significantly improve the forecasts. In addition, data assimilation is required to improve the simulations and correct the dissipative effect of the numerical scheme.

The sediment transport modelling is of major interest in terms of applications, in particular to estimate the sustainability of facilities with silt or scour, such as canals and bridges. Dredging or filling-up operations are expensive and generally not efficient in the long term. The objective is to determine a configuration almost stable for the facilities. In addition, it is also important to determine the impact of major events like emptying dam which is aimed at evacuating the sediments in the dam reservoir and requires a large discharge. However, the downstream impact should be measured in terms of turbidity, river morphology and flood.

It is a violent, sudden and destructive flow. Between 1996 and 2005, nearly 80% of natural disasters in the world have meteorological or hydrological origines. The main interest of their study is to predict the areas in which they may occur most probably and to prevent damages by means of suitable amenities. In France, floods are the most recurring natural disasters and produce the worst damages. For example, it can be a cause or a consequence of a dam break. The large surface they cover and the long period they can last require the use of reduced models like the shallow water equations. In urban areas, the flow can be largely impacted by the debris, in particular cars, and this requires fluid/structure interactions be well understood. Moreover, underground flows, in particular in sewers, can accelerate and amplify the flow. To take them into account, the model and the numerical resolution should be able to treat the transition between free surface and underground flows.

Tsunamis are another hydrological disaster largely studied. Even if the propagation of the wave is globally well described by the shallow water model in oceans, it is no longer the case close to the epicenter and in the coastal zone where the bathymetry leads to vertical accretions and produces substantial dispersive effects. The non-hydrostatic terms have to be considered and an efficient numerical resolution should be induced.

While viscous effects can often be neglected in water flows, they have to be taken into account in
situations such as avalanches, debris flows, pyroclastic flows, erosion processes, ...i.e. when the fluid rheology
becomes more complex. Gravity driven granular flows consist of solid particles commonly mixed with
an interstitial lighter fluid (liquid or gas) that may interact with the grains and decrease the intensity of
their contacts, thus reducing energy dissipation and favoring propagation. Examples include subaerial or
subaqueous rock avalanches (e.g. landslides).

Nowadays, simulations of the hydrodynamic regime of a river, a lake or an estuary, are not restricted to the determination of the water depth and the fluid velocity. They have to predict the distribution and evolution of external quantities such as pollutants, biological species or sediment concentration.

The potential of micro-algae as a source of biofuel and as a technological solution for CO2 fixation is the subject of intense academic and industrial research. Large-scale production of micro-algae has potential for biofuel applications owing to the high productivity that can be attained in high-rate raceway ponds. One of the key challenges in the production of micro-algae is to maximize algae growth with respect to the exogenous energy that must be used (paddlewheel, pumps, ...). There is a large number of parameters that need to be optimized (characteristics of the biological species, raceway shape, stirring provided by the paddlewheel). Consequently our strategy is to develop efficient models and numerical tools to reproduce the flow induced by the paddlewheel and the evolution of the biological species within this flow. Here, mathematical models can greatly help us reduce experimental costs. Owing to the high heterogeneity of raceways due to gradients of temperature, light intensity and nutrient availability through water height, we cannot use depth-averaged models. We adopt instead more accurate multilayer models that have recently been proposed. However, it is clear that many complex physical phenomena have to be added to our model, such as the effect of sunlight on water temperature and density, evaporation and external forcing.

Many problems previously mentioned also arise in larger scale systems like lakes. Hydrodynamics of lakes is mainly governed by geophysical forcing terms: wind, temperature variations, ...

One of the booming lines of business is the field of renewable and decarbonated energies. In particular in the marine realm, several processes have been proposed in order to produce electricity thanks to the recovering of wave, tidal and current energies. We may mention water-turbines, buoys turning variations of the water height into electricity or turbines motioned by currents. Although these processes produce an amount of energy which is less substantial than in thermal or nuclear power plants, they have smaller dimensions and can be set up more easily.

The fluid energy has kinetic and potential parts. The buoys use the potential energy whereas the water-turbines are activated by currents. To become economically relevant, these systems need to be optimized in order to improve their productivity. While for the construction of a harbour, the goal is to minimize swell, in our framework we intend to maximize the wave energy.

This is a complex and original issue which requires a fine model of energy exchanges and efficient numerical tools. In a second step, the optimisation of parameters that can be changed in real-life, such as bottom bathymetry and buoy shape, must be studied. Eventually, physical experiments will be necessary for the validation.

The urban environment is essentially studied for air and noise pollutions. Air pollution levels and noise pollution levels vary a lot from one street to next. The simulations are therefore carried out at street resolution and take into account the city geometry. The associated numerical models are subject to large uncertainties. Their input parameters, e.g. pollution emissions from road traffic, are also uncertain. Quantifying the simulation uncertainties is challenging because of the high computational costs of the numerical models. An appealing approach in this context is the use of metamodels, from which ensembles of simulations can be generated for uncertainty quantification.

The simulation uncertainties can be reduced by the assimilation of fixed and mobile sensors. High-quality fixed monitoring sensors are deployed in cities, and an increasing number of mobile sensors are added to the observational networks. Even smartphones can be used as noise sensors and dramatically increase the spatial coverage of the observations. The processing and assimilation of the observations raises many questions regarding the quality of the measurements and the design of the network of sensors.

There is a growing interest for environmental problems at city scale, where a large part of the population is
concentrated and where major pollutions can occur. Numerical simulation is well established to study the
urban environment, e.g. for road traffic modelling. As part of the smartcity movement, an increasing number
of sensors collect measurements, at traditional fixed observation stations, but also on mobile devices, like
smartphones. They must properly be taken into account given their number but also their potential low quality.

Pratical applications include air pollution and noise pollution. These directly relate to road traffic. Data assimilation and uncertainty propagation are key topics in these applications.

Only few travels were done last year (including one flight) as a consequence of a will of the team to avoid this type of transportation.

Part of ANGE activity is devoted to research on renewable energy. In this way, the team took part to the organization of the EMRSim 22 conference, which devoted to Marine Energy Techniques and Simulation.

2023 was a year of evaluation for our project-team. This evaluation, coordinated by Inria's Evaluation Committee, is an important step in the team's life cycle. In particular, it enabled us to elaborate collectively on the future of all teams involved in the environmental sciences theme. We have produced a summary of this foresight, which raises more questions than it answers, but which reflects the topic's complexity and the debate that has opened up within the institute itself.

This report can be found here.

Polyphemus is a modeling system for air quality. As such, it is designed to yield up-to-date simulations in a reliable framework: data assimilation, ensemble forecast and daily forecasts. Its completeness makes it suitable for use in many applications: photochemistry, aerosols, radionuclides, etc. It is able to handle simulations from local to continental scales, with several physical models. It is divided into three main parts:

libraries that gather data processing tools (SeldonData), physical parameterizations (AtmoData) and post-processing abilities (AtmoPy),

programs for physical pre-processing and chemistry-transport models (Polair3D, Castor, two Gaussian models, a Lagrangian model),

model drivers and observation modules for model coupling, ensemble forecasting and data assimilation.

Coll. with Olivier Bernard
In 8, we focus on a periodic resource allocation problem applied to a dynamical system which comes from a biological system. More precisely, we consider a system with

other participants: Elise Basquin, Alain Rabaute, Maud Thomas, Sara Lafuerza, Abdelmounim El M'Rini, Denis Mercier, Elia D'acremont, Marie-Odile Bristeau, Axel Creach

The Alboran Basin may be subject to tsunami hazards. If such an event were to occur, it is expected that the urbanised and densely populated areas of northern Moroccan coastline would be affected. Precise inundation hazard maps are needed for tsunami risk management in this region. In 6, we argue that the diversity of hazard mapping methods ensures the robustness of the scientific knowledge about the exposure of a territory. Hence, the main objective of this study is to analyse the exposure of the plain of Martil (north of Morocco), by using four hazard mapping methods to create inundation maps for two scenarios of tsunamis generated by extreme submarine mass failure (SMF) in the Alboran Sea, of 0.9 km3 and 3.8 km3 respectively. A digital terrain model of the plain was used to explore four methods of inundation mapping. The static method identified 4.32 km2 and 19.83 km2 of flooded areas for each scenario using water height values as inundation thresholds. The hybrid and the volumetric methods use the volume of water to determine the inundation extent. For the first scenario, 3.51 km2 of the plain were inundated using the hybrid method, and 20.11 km2 for the second scenario. The results of the volumetric methods are 2.32 km2 and 7.82 km2 respectively for the first and second scenario. Finally, the fourth method relies on numerical hydrodynamic modelling of tsunami inundation (Freshkiss3d® code). With this method, 4.55 km2 of the plain were flooded in the first scenario, and 24.12 km2 for the second. The comparison of the results highlights that the most sensitive areas to tsunami inundation are the lowest topographic ones, being the beaches and the wadis floodplains. This result raises questions on the current coastal development and the preparedness of its population, thus calling for more attention to engage on tsunami risk management related questions.

other participant: Etienne Mémin

In 10, we characterize numerically through two criteria the Lévy area related to unresolved fluctuation velocities associated to a stochastic coarse-scale representation of geophysical fluid flow dynamics. We study in particular whether or not the process associated to the random unresolved velocity components exhibits a Lévy area corresponding to a Wiener process, and if the law of this process can reasonably be approached by a centered Dirac measure. This exploration enables us to answer positively to a conjecture made for the constitution of high-order discrete time evolution schemes for stochastic representation defined from stochastic transport.

other participant: Patrick Gallinari

Forecasting physical signals in long time range is among the most challenging tasks in Partial Differential Equations (PDEs) research. To circumvent limitations of traditional solvers, many different Deep Learning methods have been proposed. They are all based on auto-regressive methods and exhibit stability issues. Drawing inspiration from the stability property of implicit numerical schemes, we introduce in 14 a stable auto-regressive implicit neural network. We develop a theory based on the stability definition of schemes to ensure the stability in forecasting of this network. It leads us to introduce hard constraints on its weights and propagate the dynamics in the latent space. Our experimental results validate our stability property, and show improved results at long-term forecasting for two transports PDEs.

other participant: Louis Serrano, Yuan Yin, Jocelyn Ahmed Mazari, Patrick Gallinari

For numerical design, the development of efficient and accurate surrogate models is paramount. They allow us to approximate complex physical phenomena, thereby reducing the computational burden of direct numerical simulations. In 15, we propose INFINITY, a deep learning model that utilizes implicit neural representations (INRs) to address this challenge. Our framework encodes geometric information and physical fields into compact representations and learns a mapping between them to infer the physical fields. We use an airfoil design optimization problem as an example task and we evaluate our approach on the challenging AirfRANS dataset, which closely resembles real-world industrial use-cases. The experimental results demonstrate that our framework achieves state-of-the-art performance by accurately inferring physical fields throughout the volume and surface. Additionally we demonstrate its applicability in contexts such as design exploration and shape optimization: our model can correctly predict drag and lift coefficients while adhering to the equations.

Coll. with François Bouchut
In 9, we are interested in free surface flows where density variations coming e.g. from temperature or salinity differences play a significant role in the hydrody-namic regime. In water, acoustic waves travel much faster than gravity and internal waves, hence the study of models arising from compressible fluid mechanics often requires a decoupling between these waves. Starting from the compressible Navier-Stokes system, we derive the so-called Navier-Stokes-Fourier system in an "incompressible" regime using the low-Mach scaling, hence filtering the acoustic waves, neglecting the density dependency on the fluid pressure but keeping its variations in terms of temperature and salinity. A slightly modified low-Mach asymptotics is proposed to obtain a model with thermo-mechanical compatibility. The case when the density depends only on the temperature is studied first. Then the variations of the fluid density with respect to temperature and salinity are considered, and it seems to be the first time that salinity dependency is considered in this low Mach limit. We give a layer-averaged formulation of the obtained models in an hydrostatic context, allowing to derive numerical schemes endowed with strong stability properties that are presented in a companion paper. Several stability properties of the layer-averaged Navier-Stokes-Fourier system are proved.

In 11, we construct an accurate linear model describing the propagation of both acoustic and gravity waves in water. This original model is obtained by the linearization of the compressible Euler equations, written in Lagrangian coordinates. The system is studied in the isentropic case, with a free surface, an arbitrary bathymetry, and vertical variations of the background temperature and density. We show that our model is an extension of some models from the literature to the case of a non-barotropic fluid with a variable sound speed. Other models from the literature are recovered from our model through two asymptotic analyses, one for the incompressible regime and one for the acoustic regime. We also propose a method to write the model in Eulerian coordinates. Our model includes many physical properties, such as the existence of internal gravity waves or the variation of the sound speed with depth.

Other participants: Vladimír Šverák

In 13, we consider the Cauchy problem for the incompressible Navier–Stokes equations in for a one-parameter family of explicit scale-invariant axi-symmetric initial data, which is smooth away from the origin and invariant under the reflection with respect to the xy-plane. Working in the class of axi-symmetric fields, we calculate numerically scale-invariant solutions of the Cauchy problem in terms of their profile functions, which are smooth. The solutions are necessarily unique for small data, but for large data we observe a breaking of the reflection symmetry of the initial data through a pitchfork-type bifurcation. By a variation of previous results by Jia and Šverák (Invent Math 196(1):233-265, 2013) it is known rigorously that if the behavior seen here numerically can be proved, optimal non-uniqueness examples for the Cauchy problem can be established, and two different solutions can exists for datum which is divergence-free, smooth away from the origin, compactly supported, and locally -homogeneous near the origin. In particular, assuming our (finite-dimensional) numerics represents faithfully the behavior of the full (infinite-dimensional) system, the problem of uniqueness of the Leray–Hopf solutions (with non-smooth initial data) has a negative answer and, in addition, the perturbative arguments such those by Kato (Math Z 187(4):471-480, 1984) and Koch and Tataru (Adv Math 157(1):22-35, 2001), or the weak-strong uniqueness results by Leray, Prodi, Serrin, Ladyzhenskaya and others, already give essentially optimal results. There are no singularities involved in the numerics, as we work only with smooth profile functions. It is conceivable that our calculations could be upgraded to a computer-assisted proof, although this would involve a substantial amount of additional work and calculations, including a much more detailed analysis of the asymptotic expansions of the solutions at large distances.

Other participants: Antoine Belgodere, Jean-Baptiste Filippi, Florian Guéniot

Avoiding catastrophic wildfires is a natural rationale for fighting fires in their early stage. Beside this benefit, may a marginal decrease in the duration of smaller wildfires be worthwhile? The present article addresses this topic by estimating the marginal damage of the duration of forest fires. In 7, we perform two sets of wildfire simulations in Corsica, and estimate the damage based on the type of land use in burned areas. Results suggest that the marginal cost of the duration of fires rises by a factor of 4 during the first 400 minutes. The two reasons appear to be the increase in the marginal burned area (a physical mechanism) and the increase in the value of the marginal burned area, due to the ignition points being located in low-value places (a human mechanism). Using a conservative calibration, our results corroborate the principle of early initial attack already in use in countries with sufficient fire fighting forces, but subject to debate because of its cost.

Other participants: Laura Grigori, Sever A Hirstoaga

In 12, we introduce a new strategy for solving highly oscillatory two-dimensional Vlasov-Poisson systems by means of a specific version of the parareal algorithm. The novelty consists in using reduced models, obtained from the two-scale convergence theory, for the coarse solving. The reduced models are useful to approximate the original Vlasov-Poisson model at a low computational cost since they are free of high oscillations. Both models are numerically solved in a particle-in-cell framework. We illustrate this strategy with numerical experiments based on long time simulations of a charged beam in a focusing channel and under the influence of a rapidly oscillating external electric field. On the basis of computing times, we provide an analysis of the efficiency of the parareal algorithm in terms of speedup.

Yohan Penel, who left the team in 2022, supervised the PhD thesis of Giuseppe Parasiliti about the Physical, mathematical and numerical modelling of a gas flow for the transportation of liquified natural gas. This work is the result of a close collaboration with the corporation GTT, which has already collaborated with ANGE in the last years, through the Carnot institute SMILE. The defens

Jacques Sainte-Marie has a contract with Eaux de Paris about Hydraulic modeling, calibration and diagnosis. (2020-2023, with S. Labbé, Laboratoire D'Alembert and LPSM)

The goal of the ALLOWAPP project is the design of space-time parallel algorithms for large-scale optimization problems associated with wave propagation phenomena. Such problems appear in seismology, geophysics, but also in various applications from data assimilation. The large amount of data and the volume of computations required for the accurate numerical solution of wave propagation problems, within an optimization loop, requires the use of massively parallel computers. Time-parallel methods have experienced a great development in the last ten years, and for parabolic problems an almost perfect efficiency for a large number of processors has been achieved (scalability). It is quite different for wave propagation problems. In this project, we propose to develop robust, efficient and scalable methods for space-time parallelization of these optimization problems.

The GeoFun project aims to improve the modeling and simulation of geophysical flows involving at least two different processes. Numerical simulation of watersheds and estimation of water resources is the main application of the project's achievements. In this context, a free surface flow (rivers, lakes) is the upper part of a groundwater flow (water table). Our vision of river transport is often naive, because we think first of rivers, lakes and floods, but in reality, 80 % of the water of the continents is underground. Sometimes, the porous substratum is covered by an impermeable rock stratum, which confines the flow as in pipes, except for some points where springs and resurgences appear.

The objective of SingFlows is to develop mathematical and numerical tools for the analysis of three problems in fluid dynamics: the behaviour of anisotropic flows (boundary layers, shallow water flows), the dynamics of vortical structures, and the evolution of fixed or floating structures in water waves. Our will to unify these different problems is natural, because they share many mathematical features. The underlying keypoint is that they are described by singular solutions of Euler or Navier-Stokes equations. The word singular refers here: - either to a lack of smoothness: it applies for instance to vortex filaments, which are Dirac masses along curves, or to the contact line between water and the floating structure, - or to a singular dependence of the solution with respect to a parameter, typically the Reynolds number (like in boundary layers). The connection between the two points of view is usually made by viscous regularization of the non-smooth structure, or conversely by taking the vanishing limit of the parameter. More generally, the three problems considered in SingFlows involve flows with very small scales. A relevant description then requires the derivation of reduced models.

The project aims at developing the interplay between Deep Neural Networks (DNNs) and Differential Equations (DEs), with the goal of modeling complex dynamical systems arising from the observation of natural phenomena. Two application domains are targeted, environment and healthcare. We address three fundamental questions: how to adapt and apply numerical analysis theory to DNNs for analyzing them, providing theoretical guaranties and improving their robustness, how to combine simulation and data based models into hybrid systems, how could DNNs help solving DEs and complement numerical solvers. In addition, we evaluate our methods on simulation and real world data in the environment and health domains. DeepNuM gathers partners with complementary skills: DEs and Environment (INRIA-ANGE), Machine Learning and DNNs (Sorbonne – MLIA), DEs and Biophysics (INRIA-EPIONE).

Giant submarine landslides (10-2000 km3) are found in the thick Quaternary sediment succession of passive continental margins. Their ages coincide with periods of sea-level fall and rise, but it is unclear how such vast failures can be triggered on low seafloor slopes (<2?) in the absence of a triggering factor such as seismicity. Key hypotheses involve excess pore pressures linked to reductions in gas-hydrate stability, driven by changes either in climate or in subsurface fluid flow. The MEGA project wants to explore such hypotheses through the first modelling of linked changes in gas hydrate and slope stability in response to ocean pressure and temperature changes, using an innovative comparison of the Nile and Amazon deep-sea fans that experience different forms of climate forcing over glacial-interglacial timescales. As such megaslides have never triggered in historical times, MEGA will provide input for the first modelling of their tsunamogenic consequences on coastal zones.

The MathGeoPhy interdisciplinary research group was created in January 2022, for five years. It is funded by the French National Center for Scientific Research (CNRS), with the mission of animating the French scientific community around the theme of mathematics in interaction with the geophysics of fluid and solid envelopes. The members of the GdR are interested in mathematical modeling, scientific computing and the development of new numerical methods applied in particular to :

The objectives of this project are the following:

The GDR is an entry and orientation point. It provides a forum for the exchange of information concerning industrial needs and the kills of the academic community; and enables the bringing together of players. Valuing the national scientific community The GDR EMR gives visibility to the community, in particular through the development of a mapping of the actors and themes available on the web platform

Member of the organizing committees

The team members served as referees for the following journals:

JS is editor in chief of MATAPLI (national journal of the applied maths community). EA is responsible for its section "Du côté des écoles d’ingénieurs".

EG is the president of commission française pour l'enseignement des mathématiques (CFEM)

Teaching activities of ANGE are summarized in the following.

The format is here supervisior, grade, student, institution, period, title.