Together with the diffusion of scientific computing, there has been a recent and impressive increase of the demand for numerical methods. The problems to be addressed are everyday more complex and require specific numerical algorithms. The quality of the results has to be accurately assessed, so that in-silico experiments results can be trusted. Nowadays, producing such reliable numerical results goes way beyond the abilities of isolated researchers, and must be carried out by structured teams.

The topics addressed by the RAPSODI project-team belong to the broad theme of numerical methods for the approximation of solutions of systems of partial differential equations (PDEs). Besides standard convergence properties, a good numerical method for approximating a physical problem has to satisfy at least the following three criteria:

We aim to develop methods fulfilling the above quality criteria for physical models which all display a dissipative behavior, and that are motivated by industrial collaborations or multidisciplinary projects. In particular, we have identified a couple of specific situations we plan to investigate: models from corrosion science (in the framework of nuclear waste repository) , low-frequency electromagnetism , and mechanics of complex inhomogeneous fluids arising in avalanches or in porous media .

Ideally, we should allow ourselves to design entirely new numerical methods. For some applications however (often in the context of industrial collaborations), the members of the team have to compose with existing codes. The numerical algorithms have thus to be optimized under this constraint.

Some technological bottlenecks related to points (a)–(c) mentioned above are well identified. In particular, it appears that a good numerical method should handle general meshes, so that dynamic mesh adaptation strategies can be used in order to achieve (b). But it should also be of the highest possible order while remaining stable in the sense of (a), and robust in the sense of (c). There have been numerous research contributions on each point of (a)–(c) in the last decades, in particular for solving each difficulty separately, but combining them still leads to unsolved problems of crucial interest.

Let us mention for example the review paper by Jérôme Droniou , where it is highlighted that all the linear methods for solving diffusion equations on general meshes suffer from the same lack of monotonicity and preserve neither the positivity of the solutions nor the decay of the entropy. Moreover, there is no complete convergence proof for the nonlinear methods exposed in . The first convergence proof for a positivity preserving and entropy diminishing method designed to approximate transient dissipative equations on general meshes was proposed recently in . The idea and the techniques introduced in should be extended to practical applications.

In systems of PDEs, the values of the physical parameters often change the qualitative behavior of the solutions.
Then, one challenge in the numerical approximation of such systems is the design of methods which can be applied for a large range of parameters,
as in particular in the regime of singular perturbations.
Such schemes, called asymptotic-preserving (AP) schemes , are powerful tools as they allow the use of the same scheme for a given problem and for its limit with fixed discretization parameters. In many cases, the AP property of numerical schemes is just empirically established, without any rigorous proof. We aim to extend the techniques recently introduced in for the drift-diffusion system, and based on the control of the numerical dissipation of entropy, to other dissipative systems in order to prove the AP property of numerical schemes.

The question of the robustness of the numerical methods with respect to the physical parameters is also fundamental for fluid mixture models. The team already developed such schemes for the variable density Navier–Stokes system , . We aim to propose new ones for more complex models with the same philosophy in mind. On the one hand, we will be interested in high-order schemes, which must be as simple as possible in view of 3D practical implementations. Let us stress that combining high order accuracy and stability is very challenging. On the other hand, the optimization of the computations will have to be considered, in particular with the development of some a posteriori error estimators. Impressive progresses have been achieved in this field , allowing important computational savings without compromising the accuracy of the results. Recently, we successfully applied this strategy to the Reissner–Mindlin model arising in solid mechanics , the dead-oil model for porous media flows , or the Maxwell equations in their quasi-static approximation for some eddy current problems , . We aim to develop new a posteriori estimators for other dissipative systems, like fluid mixture models.

In a nutshell, our goal is to take advantage of and extend the most recent breakthroughs of the mathematical community to tackle in an efficient way some application-guided problems coming either from academics or from industrial partners. To this end, we shall focus on the following objectives, which are necessary for the applications we have in mind:

Up to now, almost all numerical methods dedicated to degenerate parabolic problems that the mathematicians are able to analyze rely on the use of mathematical transformations (like, e.g., the Kirchhoff's transform). It forbids the extension of the analysis to complex realistic models. The methods used in the industrial codes for solving such complex problems rely on the use of what we call NNM, i.e., on methods that preserve all the nonlinearities of the problem without reducing them thanks to artificial mathematical transforms. Our aim is to take advantage of the recent breakthrough proposed by C. Cancès and Cindy Guichard in and to develop efficient new numerical methods with a full numerical analysis (stability, convergence, error estimates, robustness with respect to physical parameters, etc).

There has been an extensive effort in the recent years to develop numerical methods for diffusion equations that are robust with respect to heterogeneities, anisotropy, and the mesh (see for instance for a comprehensive discussion on such methods). On the other hand, the understanding of the role of nonlinear stability properties in the asymptotic behaviors of dissipative systems increased significantly in the last decades (see for instance , ).

Recently, C. Chainais-Hillairet and co-authors , , developed
a strategy based on the control of the numerical counterpart
of the physical entropy to develop and analyze AP numerical methods.
In particular, these
methods show great promises for capturing accurately the behavior of the solutions to dissipative problems
when some physical parameter is small with respect to the discretization characteristic parameters, or in the long-time asymptotics.
Since it requires the use of nonlinear test functions in the analysis,
strong restrictions on the physics (isotropic problems) and on the mesh (Cartesian grids, Voronoï boxes, etc)
are required in , , . The schemes proposed in and allow
to handle nonlinear test functions in the analysis without restrictions on the mesh and on the anisotropy of the problem.
Combining the nonlinear schemes à la
with the methodology of , , would provide schemes that are robust both
with respect to the meshes and to the parameters. Therefore, they would also be robust under adaptive mesh refinement.

We aim at extending the range of the NS2DDV software by introducing new physical models, like for instance the low-Mach model, which gives intermediate solutions between the compressible Navier–Stokes model and the incompressible Navier–Stokes one.
This model was introduced in as a limit system which describes combustion processes at low Mach number in a confined region.
Within this scope, we will propose a theoretical study for proving the existence of weak solutions for a particular class of models for which the dynamic viscosity of the fluid is a specific function of the density. We will also propose the extension of a combined Finite Volume-Finite Element method, initially developed for the simulation of incompressible and variable density flows, to this class of models.

The numerical experiments carried out in show that in case of very strong anisotropy, the convergence of the proposed NNM becomes too slow (less than first order). Indeed, the method appears to strongly overestimate the dissipation. In order to make the method more competitive, it is necessary to estimate the dissipation in a more accurate way. Preliminary numerical results show that second-order accuracy in space can be achieved in this way. One also aims to obtain (at least) second-order accuracy in time without jeopardizing the stability. For many problems, this can be done by using so-called two-step backward differentiation formulas (BDF2) .

Concerning the inhomogeneous fluid models, we aim to investigate new methods for the solution of the mass equation. Indeed, we aim at increasing the accuracy while maintaining some positivity-like properties and the efficiency for a wide range of physical parameters. To this end, we will consider Residual Distribution schemes, that appear as an alternative to Finite Volume methods. Residual Distribution schemes enjoy very compact stencils. Therefore, their extension from 2D to 3D entails reasonable difficulties. These methods appeared twenty years ago, but recent extensions to unsteady problems , , with high-order accuracy , , or for parabolic problems , make them very competitive. Relying on these breakthroughs, we aim at designing new Residual Distribution schemes for fluid mixture models with high-order accuracy while preserving the positivity of the solutions.

The question of the a posteriori error estimators will also have to be addressed in this optimization context. Since the pioneering papers of Babuška and Rheinboldt more than thirty years ago , a posteriori error estimators have been widely studied. We will take advantage of the huge corresponding bibliographical database in order to optimize our numerical results.

For example, we would like to generalize the results we derived for the harmonic magnetodynamic case (e.g., , ) to the temporal magnetodynamic one, for which space-time a posteriori error estimators have to be developed. A space-time refinement algorithm should consequently be proposed and tested on academic as well as industrial benchmarks.

We also want to develop a posteriori estimators for the variable density Navier–Stokes model or some of its variants. To do so, several difficulties have to be tackled: the problem is nonlinear, unsteady, and the numerical method , we developed combines features from Finite Elements and Finite Volumes. Fortunately, there exists a significant literature on the subject. Some recent references are devoted to the unsteady Navier–Stokes model in the Finite Element context , . In the Finite Volume context, recent references deal with unsteady convection-diffusion equations , , , . We want to adapt some of these results to the variable density Navier–Stokes system, and to be able to design an efficient space-time remeshing algorithm.

Many systems are modeled as a large number (

Porous media flows are of great interest in many contexts, like, e.g., oil engineering, water resources management, nuclear waste repository management, or carbon dioxide sequestration. We refer to , for an extensive discussion on porous media flow models.

From a mathematical point of view, the transport of complex fluids in porous media often leads to possibly degenerate parabolic conservation laws. The porous rocks can be highly heterogeneous and anisotropic. Moreover, the grids on which one intends to solve numerically the problems are prescribed by the geological data, and might be nonconformal with cells of various shapes. Therefore, the schemes used for simulating such complex flows must be particularly robust.

The team is interested in the theoretical and numerical analysis of mathematical models describing the degradation of materials, as concrete carbonation and corrosion. The study of such models is an important environmental and industrial issue. Atmospheric carbonation degrades reinforced concretes and limits the lifetime of civil engineering structures. Corrosion phenomena issues occur for instance in the reliability of nuclear power plants and the nuclear waste repository. The study of the long time evolution of these phenomena is of course fundamental in order to predict the lifetime of the structures.

From a mathematical point of view, the modeling of concrete carbonation (see ) as the modeling of corrosion in an underground repository (DPCM model developed by Bataillon et al. ) lead to systems of PDEs posed on moving domains. The coupling between convection-diffusion-reaction equations and moving boundary equations leads to challenging mathematical questions.

The team is interested in numerical methods for the simulation of systems of PDEs describing complex flows, like for instance mixture flows, granular gases, rarefied gases, or quantum fluids.

Variable-density, low-Mach flows have been widely studied in the recent literature because of their applicability in various phenomena such as flows in high-temperature gas reactors, meteorological flows, flows with convective and/or conductive heat transfer or combustion processes. In such cases, the resolution of the full compressible Navier–Stokes system is not adapted, because of the sound waves' speed. The Boussinesq incompressible model is not a better alternative for such low-speed phenomena, because the compressibility effects cannot be totally cancelled due to large variations of temperature and density. Consequently, some models have been formally derived, leading to the filtering of the acoustic waves by the use of some formal asymptotic expansions and two families of methods have been developed in the literature in order to compute these flows. We are interested in particular in the so-called pressure-based methods, which are more robust than density-based solvers, although their range of validity is in general more limited.

The kinetic theory of molecular gases models a gas as a system of elastically colliding spheres, conserving mechanical energy during impact. Once initialized, it takes to a molecular gas no more than a few collisions per particle to relax to its equilibrium state, characterized by a Maxwellian velocity distribution and a certain homogeneous density (in the absence of external forces). A granular gas is a system of dissipatively colliding, macroscopic particles (grains). This slight change in the microscopic dynamics (converting energy into heat) causes drastic changes in the behavior of the gas: granular gases are open systems, which exhibit self-organized spatio-temporal cluster formations, and have no equilibrium distribution. They can be used to model silos, avalanches, pollen or planetary rings.

Quantum models can be used to describe superfluids, quantum semiconductors, weakly interacting Bose gases, or quantum trajectories of Bohmian mechanics. They have attracted considerable attention in the last decades, due in particular to the development of nanotechnology applications. To describe quantum phenomena, there exists a large variety of models. In particular, there exist three different levels of description: microscopic, mesoscopic, and macroscopic. The quantum Navier–Stokes equations deal with a macroscopic description in which the quantum effects are taken into account through a third-order term called the quantum Bohm potential. This Bohm potential arises from the fluid dynamical formulation of the single-state Schrödinger equation. The nonlocality of quantum mechanics is approximated by the fact that the equations of state do not only depend on the particle density but also on its gradient. These equations were employed to model field emissions from metals and steady-state tunneling in metal-insulator-metal structures, and to simulate ultra-small semiconductor devices.

The knowledge of the geology is a prerequisite before simulating flows within the subsoil. Numerical simulations of the geological history thanks to stratigraphy numerical codes allow to complete the knowledge of the geology where experimental data are lacking. Stratigraphic models consist in a description of the erosion and sedimentation phenomena at geological scales.

The characteristic time scales for the sediments are much larger than the characteristic time scales for the water in the river. However, the (time-averaged) water flux plays a crucial role in the evolution of the stratigraphy. Therefore, defining appropriate models that take the coupling between the rivers and the sediments into account is fundamental and challenging. Once the models are at hand, efficient numerical methods must be developed.

Numerical simulation is nowadays an essential tool in order to design electromagnetic systems, by estimating the electromagnetic fields generated in a wide variety of devices. An important challenge for many applications is to quantify the intensity of the electric field induced in a conductor by a current generated in its neighborhood. In the low-frequency regime, we can for example cite the study of the impact on the human body of a high-tension line or, for higher frequencies, the one of a smartphone. But the ability to simulate accurately some electromagnetic fields is also very useful for nondestructive control, in the context of the maintenance of nuclear power stations for example. The development of efficient numerical tools, among which a posteriori error estimators, is consequently necessary to reach a high precision of calculation in order to provide estimations that are as reliable as possible.

C. Cancès was promoted senior Inria researcher (DR2) as of October 2021.

I. Lacroix-Violet was promoted from associate professor at Université de Lille to full professor at Université de Lorraine (Nancy, France). She left the project-team at the end of August 2021.

On November 23, 2021 C. Cancès was awarded the from the Académie des Sciences: see . The Blaise Pascal Prize is awarded every year by the Académie des Sciences after consultation of the SMAI-GAMNI group. It aims at rewarding outstanding achievements on the devising and mathematical analysis of numerical methods for the solution of PDEs realized by an under-40 researcher in France.

From November 16 to November 19, 2021 was held in Lille (at Polytech Lille engineering school) the 4th edition of the conference "Asymptotic Behaviors of systems of PDEs arising in Physics and Biology" (), organized by C. Cancès, C. Chainais-Hillairet, M. Herda, I. Lacroix-Violet, Alexandre Mouton (Université de Lille), and T. Rey. This event gathered around 70 participants. The ABPDE conference has now become a well-established recurrent event (every 3 years); see also Section .

For the self-assessment of our platforms and codes, we adopt the framework defined by (Software family, Audience, Evolution and maintenance, Duration of the development, Contribution of the team, Web page, Description).

Family={research, vehicle}; Audience={partners}; Evolution={lts}; Duration={2}; Contribution={leader}; URL={}

ParaSkel++ is a C++ platform (freely distributed under LGPL), developed by L. Beaude and S. Lemaire, for the high-performance, arbitrary-order, 2/3D numerical approximation of PDEs on general polytopal meshes using skeletal Galerkin methods.
Skeletal Galerkin methods are a vast family of numerical methods for the approximation of PDE-based models that satisfy the following two building principles (see ):

The very structure underpinning skeletal methods grants them the property of being amenable to static condensation, i.e., locally to each cell, bulk DOF can be eliminated in terms of the local skeleton DOF by means of a Schur complement. The final global system to solve thus writes in terms of the skeleton DOF only.
The skeletal family encompasses in particular standard FE methods and virtual-like Galerkin methods (VEM, HHO, HDG...). It does not contain (plain vanilla) DG methods.
ParaSkel++ offers a high-performance factorized C++ architecture for the implementation of arbitrary-order skeletal methods on general 2/3D polytopal meshes.
A first version (v1, August 2021) of the platform is operational, featuring a sequential implementation of all the main skeletal methods. The next crucial development steps will be the parallelization on shared memory (before considering distributed memory), and the implementation of efficient quadrature formulas on polytopal cells.
Eventually, the ParaSkel++ platform is expected to possess five main assets with respect to other codes of the same nature from the community: (i) a unified 2/3D implementation, (ii) the native support of any type of DOF (vertex-, edge-, face-, and cell-based), (iii) an ultra-factorized architecture (with common-to-all-methods local elimination and global assembly steps), (iv) the use of efficient quadrature formulas on general polytopal cells (without the need for subtessellation), and (v) the embedding of parallel computation facilities.

Family={vehicle}; Audience={partners}; Evolution={lts}; Duration={2}; Contribution={leader};

URL={}

Family={vehicle}; Audience={team}; Evolution={nofuture}; Duration={1}; Contribution={leader};

URL={}

The Matlab code FPfrac, developed by M. Herda, was used to produce the numerical simulations of the article "On a structure-preserving numerical method for fractional Fokker–Planck equations" .

Family={vehicle}; Audience={team}; Evolution={nofuture}; Duration={0.5}; Contribution={instigator}; URL={}

The Julia code FV4SM , developed by C. Cancès, Virginie Ehrlacher (ENPC & Inria Paris), and Laurent Monasse (Inria Nice), was used to produce the numerical simulations of the article "Finite Volumes for the Stefan–Maxwell cross-diffusion system" . It simulates the Stefan–Maxwell equations on Cartesian grids thanks to an entropy-diminishing TPFA Finite Volume scheme.

Family={vehicle}; Audience={partners}; Evolution={basic}; Duration={1}; Contribution={instigator};

URL={}

The Matlab code DPCM-proof, developed by Maxime Breden (École Polytechnique), C. Chainais-Hillairet, and Antoine Zurek (Université de Technologie de Compiègne), was used to perform the computer-assisted proofs of the article "Existence of traveling wave solutions for the Diffusion Poisson Coupled Model: a computer-assisted proof" .

In , C. Cancès and his co-author prove the existence of weak solutions to a system of two diffusion equations that are coupled by a pointwise volume constraint. The time evolution is given by gradient dynamics for a free energy functional. Their primary example is a model for the demixing of polymers, the corresponding energy is the one of Flory, Huggins and de Gennes. Due to the nonlocality in the equations, the dynamics considered here is qualitatively different from the one found in the formally related Cahn–Hilliard equations. Their angle of attack stems from the theory of optimal mass transport, that is, they consider the evolution equations for the two components as two gradient flows in the Wasserstein distance with one joint energy functional that has the volume constraint built in. The main difference with their previous work is the nonlinearity of the energy density in the gradient part, which becomes singular at the interface between pure and mixed phases.

In , C. Chainais-Hillairet et al. develop their work . They establish a priori estimates which lead to the existence of a solution to the scheme, and they prove the exponential decay of the discrete relative entropy towards the thermal equilibrium. Moreover, numerical results assess the good behavior of the numerical schemes.

In , C. Cancès and his co-author study a time-implicit Finite Volume scheme for the degenerate Cahn–Hilliard model proposed in and studied mathematically in . The scheme is shown to preserve the key properties of the continuous model, namely mass conservation, positivity of the concentrations, the decay of the energy, and the control of the entropy dissipation rate. This allows to establish the existence of a solution to the nonlinear algebraic system corresponding to the scheme. Furthermore, thanks to compactness arguments, the approximate solution is shown to converge towards a weak solution to the continuous problem as the discretization parameters tend to 0. Numerical results illustrate the behavior of the numerical scheme.

In , I. Lacroix-Violet et al. focus on the numerical integration in time of nonlinear Schrödinger equations using different methods preserving the energy or a discrete analog of it. In particular, they give a rigorous proof of the order of the relaxation method (presented in for cubic nonlinearities) and they propose a generalized version that allows to deal with general power law nonlinearities. Numerical simulations for different physical models show the efficiency of these methods.

In , S. Lemaire presents a unifying viewpoint on Hybrid High-Order (HHO) and Virtual Element (VE) methods on general polytopal meshes in dimension 2 or 3, in terms of both formulation and analysis. The focus is on a model Poisson problem. To bridge the two paradigms, (i) he transcribes the (conforming) VE method into the HHO framework and (ii) proves

In , S. Lemaire et al. establish the equivalence between the Multiscale Hybrid-Mixed (MHM) and the Multiscale Hybrid High-Order (MsHHO) methods for a variable diffusion problem with piecewise polynomial source term. Under the idealized assumption that the local problems defining the multiscale basis functions are exactly solved, they prove that the equivalence holds for general polytopal (coarse) meshes and arbitrary approximation orders. They finally leverage the interchange of properties to perform a unified convergence analysis, as well as to improve on both methods.

In , I. Lacroix-Violet et al. consider global weak solutions to compressible Navier–Stokes–Korteweg equations with density dependent viscosities, in a periodic domain

In , B. Merlet et al. establish new results on the approximation of

The PhD thesis of S. Bassetto was funded by IFPEn. The contract followed the lines of the bilateral contract between Inria and IFPEn. S. Bassetto defended on December 16, 2021.

The PhD thesis of M. Jonval, that started in October 2021, is co-funded by Inria (salaries) and IFPEn (overhead costs). The contract follows the lines of the bilateral contract between Inria and IFPEn.

CEA (Christian Bataillon) and ANDRA (Laurent Trenty) are involved in the EURAD project on corrosion modeling together with the RAPSODI project-team (C. Cancès, C. Chainais-Hillairet, B. Merlet, and F. Raimondi). More details on the project can be found in Section .

Between October 19 and November 19, C. Chainais-Hillairet and M. Herda invited Patricio Farrell, research group leader at WIAS (Berlin, Germany), and Dilara Abdel, PhD student within P. Farrell's group, to work on the numerical simulation of perovskite semiconductors. P. Farrell was invited as a visiting professor funded by the LabEx CEMPI (see Section ), whereas D. Abdel received funding from the French Embassy in Berlin.

The team was also visited by national researchers.

S. Lemaire invited Théophile Chaumont-Frelet (Inria Nice) to visit him in Lille on November 8-10.

C. Cancès invited Flore Nabet (École Polytechnique) to visit him in Lille on June 15-17, as well as on December 5-6.

C. Cancès, C. Chainais-Hillairet and B. Merlet are involved in the H2020 project (EUropean Joint Programme on RADioactive Waste Management). Inside EURAD, the work-package is concerned with the development and improvement of numerical methods and tools for modelling coupled processes. The task of the RAPSODI project-team inside EURAD/DONUT is to establish an energetic formulation of the Diffusion Poisson Coupled Model leading to new long-time robust numerical methods for the simulation of the corrosion processes in an underground repository. The project started in 2019, and the RAPSODI project-team received a grant of 138 750 euros. The first technical report (2020) is available . The post-doc position of F. Raimondi is funded by EURAD.

M. Herda is the French P. I. of a bilateral French-Austrian 2021 program. The two-year project is entitled "Design and analysis of structure-preserving numerical schemes for cross-diffusion systems" and has been submitted in collaboration with an Austrian research team at the Institute for Analysis and Scientific Computing, T. U. Vienna. The project involves other members of the RAPSODI project-team (C. Cancès, C. Chainais-Hillairet, B. Gaudeul, and T. Rey). The grant of 4 400 euros is dedicated to cover travel expenses.

C. Cancès and M. Herda are members of the ANR JCJC project . This project aims at gathering PDE analysts, probability theorists, and theoretical physicists to work on the derivation of macroscopic properties of physical systems from their microscopic description. The rigorous microscopic description of moving interfaces, the understanding of macroscopic nonlocal effects, and the mathematical apprehension of the underlying atomic mechanisms, are particularly important matters of this project.

C. Cancès is a member of the ANR JCJC project . This project focuses on the mathematical and numerical study of cross-diffusion systems in moving domains. The targeted application is the simulation of the production of photovoltaic devices by a vapor deposition process.

C. Chainais-Hillairet and T. Rey are members of the ANR JCJC project . This project is related to the analysis and simulation of multiscale models of semiconductors. As almost all current electronic technology involves the use of semiconductors, there is a strong interest for modeling and simulating the behavior of such devices, which was recently reinforced by the development of organic semiconductors used for example in solar panels or in mobile phones and television screens (among others).

Through their affiliation to the Laboratoire Paul Painlevé of Université de Lille, RAPSODI team members benefit from the support of the LabEx .

The "Laboratoire d'Excellence" CEMPI (Centre Européen pour les Mathématiques, la Physique et leurs Interactions), a project of the Laboratoire de mathématiques Paul Painlevé (LPP) and the laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM), was created in the context of the "Programme d'Investissements d'Avenir" in February 2012.
The association Painlevé-PhLAM creates in Lille a research unit for fundamental and applied research and for training and technological development that covers a wide spectrum of knowledge stretching from pure and applied mathematics to experimental and applied physics.
The CEMPI research is at the interface between mathematics and physics.
It is concerned with key problems coming from the study of complex behaviors in cold atoms physics and nonlinear optics, in particular fiber optics. It deals with fields of mathematics such as algebraic geometry, modular forms, operator algebras, harmonic analysis, and quantum groups, that have promising interactions with several branches of theoretical physics.

The post-doc position of M. Pegon is funded by the LabEx CEMPI. The research stay in Lille of Patricio Farrell (WIAS, Berlin, Germany) between October and November 2021 was also supported by the LabEx CEMPI.

C. Chainais-Hillairet has been a member of the CNRS (Nuclear power, Energy, Environment, Waste, Society) project POCO (Preuves assistées par Ordinateur pour un modèle de COrrosion) from 2020 to 2021. The project was coordinated by Maxime Breden (École Polytechnique). It focused on computer-assisted proofs for a corrosion model (see ).

J. Candau-Tilh and M. Pegon received funding (

The MANAKINEQO (MAthematical and Numerical Advances in KINetic EQuatiOns) project (referenced as R-ERCGEN-19-007-REY) is a proposal funded (116 545 euros) within the ERC Generator program from I-SITE ULNE. Between February 2020 and January 2022, T. Rey, P. I. of the project, aims at investigating mathematical properties, as well as developing efficient numerical schemes, for multiscale collisional kinetic equations of the Boltzmann type. The 18-month post-doc of R. Bailo has been funded by this grant, as well as most of the ABPDE IV conference (see Section ).

S. Lemaire is the P. I. of the ADT project ParaSkel++ (funded by Inria Lille), that started in February 2020. The aim of the project is to develop an optimized parallel C++ platform for the arbitrary-order numerical approximation of PDEs by skeletal methods on general 2/3D polytopal meshes (more details in Section ). L. Beaude was hired as a temporary engineer for this project. She left the team at the end of August 2021 for a permanent engineer position at BRGM (Orléans, France).

C. Calgaro, E. Creusé and T. Rey are members of the SQUAW (Super QUantum fluids and shAllow Water equations) project, headed by Olivier Goubet (Université de Lille) and funded by the 2021 regional call.
The SQUAW project brings together the skills of applied mathematicians from four research units of the FMHF (Fédération de recherche Mathématique des Hauts-de-France). One scientific direction of the projet is the coupling of quantum models with classical fluid models (Navier–Stokes type). The challenge is to successfully concatenate two completely different types of models, the numerical simulation of which currently requires completely different approaches. The fields of application of these models range from superfluids to semiconductors. Two Master 2 internships (one in Lille and one in Valenciennes) will be funded by the SQUAW project.

C. Cancès, C. Chainais-Hillairet, M. Herda, I. Lacroix-Violet, and T. Rey co-organized the 4th edition of the conference "Asymptotic Behaviors of systems of PDEs arising in Physics and Biology" (), that was held in Lille (at Polytech Lille engineering school) from November 16 to November 19, 2021. The event gathered around 70 participants and featured 10 plenary talks, 19 contributed talks, and 12 posters.

E. Creusé co-organized the event , that was held in Valenciennes on November 2-5, 2021.

E. Creusé also co-organized the 3rd edition of the event "Analyse Appliquée en Hauts-de-France", that was held in Valenciennes on July 6, 2021.

C. Cancès co-organized the DONUT scientific days within the EURAD project (see Section ), held online on January 27-28, 2021.

A. Natale co-organized the , that was held in La Piscine Museum (Roubaix) on November 25, 2021.

R. Bailo co-organized, together with Sergio Perez (Imperial College London, UK), a mini-symposium on the "Challenges in structure-preserving numerical methods for PDEs" at the , that was held remotely on April 6-9, 2021.

C. Cancès co-organized, together with Jakub W. Both (Bergen University, Norway), a mini-symposium on "Dissipation-driven nonlinear and coupled processes in porous media" at the SIAM Conference on Mathematical and Computational Issues in the Geosciences (), that was held remotely on June 21-24, 2021.

S. Lemaire co-organized, together with Andrea Borio (Politecnico di Torino, Italy), Ilario Mazzieri (Politecnico di Milano, Italy), and Giuseppe Vacca (Università di Milano Bicocca, Italy), a mini-symposium on the "Advances in polygonal and polyhedral methods" at the conference, that was held remotely on January 11-15, 2021.

S. Lemaire co-organized, together with Alexandre Ern (ENPC & Inria Paris) and Théophile Chaumont-Frelet (Inria Nice), a mini-symposium on "High-order face-based discretization methods" at the conference, that was held remotely on July 12-16, 2021.

C. Chainais-Hillairet was a member of the scientific committees of the tenth held in June 2021 in La Grande-Motte, and of the to be held in Paris in January 2022.

C. Chainais-Hillairet is a member of the editorial board of the .

RAPSODI permanent team members are regular reviewers for all the main international journals in PDEs, numerical analysis, and scientific computing.

R. Bailo gave an online talk at the kinetic mini-symposium (co-organized by himself) of the held remotely on April 6-9. He also gave an online talk in the ANEDP seminar of the Laboratoire Paul Painlevé.

C. Cancès was an invited speaker at the held in Madrid between August 30 and September 3. He also gave a contributed talk at the conference held online on September 6-9. C. Cancès was finally invited to take part in the Oberwolfach workshop held online on February 22-26.

J. Candau-Tilh presented a poster at the event held on December 8-10 in Nancy.

C. Chainais-Hillairet was an invited keynote speaker at the conference held online on September 6-9. She also gave a seminar at the Université de Nantes.

B. Gaudeul gave online talks at the in WIAS (Berlin, Germany), as well as at the PhD Students Seminar of the Laboratoire Jacques-Louis Lions (Sorbonne Université).

M. Herda gave a talk in the Summer School , within the mini-symposium "Kinetic approaches in biological systems". The event was originally planned at the Mittag-Leffler Institute (Stockholm, Sweden) but was eventually held online on July 19-23. M. Herda also gave an online talk at the (Mathematical MOdelling in Ecology) held remotely on April 2.

I. Lacroix-Violet gave an online seminar for the PDEs Webinar at Université de Lorraine.

S. Lemaire gave online talks at the conference held remotely on January 11-15, as well as in the conference (in the mini-symposium co-organized by himself) held remotely on July 12-16, 2021.

B. Merlet gave talks at the workshop held between August 30 and September 3 in Cortona (Italy), and at the workshop held on September 6-10 at the Centro di Ricerca Matematica Ennio De Giorgi in Pisa (Italy). He also gave an online talk for the PDE team seminar at Université de Poitiers.

J. Moatti presented posters at the conference held online on September 6-9, at the held in École Polytechnique (Palaiseau) on October 27-29, and for the conference held in Lille on November 16-19.

A. Natale gave a talk at the Oberwolfach workshop held online on February 22-26. He also gave a talk for the held in the Laboratoire Paul Painlevé on July 1st, as well as for the ANEDP seminar of the Laboratoire Paul Painlevé.

M. Pegon gave talks for the held online on February 4, 2021, for the held online on March 24-26, at the tenth held on June 21-25 in La Grande-Motte, at the workshop held on September 6-10 at the Centro di Ricerca Matematica Ennio De Giorgi in Pisa (Italy), at the held on November 25 in La Piscine Museum (Roubaix), and at the event held on December 8-10 in Nancy. He also gave an online seminar at the Virginia Commonwealth University (see ).

F. Raimondi gave a talk for the held online on February 4, 2021. She also presented posters at the held in École Polytechnique (Palaiseau) on October 27-29, as well as for the conference held in Lille on November 16-19.

T. Rey gave online talks at the kinetic mini-symposium (co-organized by R. Bailo) of the held remotely on April 6-9, and at the conference held remotely on June 21-25. He also gave an online seminar for the KinetiCAM work group of Cambridge University, as well as a talk in the Applied Mathematics Seminar at Université de Strasbourg.

C. Cancès is the leader of the task "Numerical methods for high-performance computing of coupled processes" in the DONUT work-package on the development and improvement of numerical methods and tools for modelling coupled processes within the H2020 project EURAD on the management of nuclear waste at the European level (see Section ).

C. Calgaro was part of the selection committee for a tenured assistant professor (MCF) position at Université du Littoral Côte d'Opale (Calais), as well as of the selection committee for 4 tenured junior research scientist (CRCN and ISFP) positions at Inria Lille.

C. Chainais-Hillairet was part of the selection committee for a tenured assistant professor (MCF) position at Université de Lille, as well as of the selection committee for a full professor (PR) position at Université de Versailles Saint-Quentin-en-Yvelines.

E. Creusé was the president of the selection committee for a temporary assistant professor (MCF article 19) position at Université Polytechnique Hauts-de-France (Valenciennes).

E. Creusé was also an HCERES expert for the evaluation of the Laboratoire de Mathématiques Jean Leray in Nantes.

I. Lacroix-Violet was part of the selection committee for a tenured assistant professor (MCF) position at Université de Montpellier, as well as of the selection committee (together with C. Calgaro) for a tenured assistant professor (MCF) position at Université du Littoral Côte d'Opale (Calais).

C. Cancès is a member of the scientific advisory board (BSC) of the Inria Lille research center.

C. Chainais-Hillairet is vice-director of the Laboratoire Paul Painlevé, in charge of human resources for researchers, professors and assistant professors. She was also an elected member of the Conseil de la Faculté des Sciences et Technologies at Université de Lille until the end of April 2021.

E. Creusé is the director of the Département de Mathématiques et Applications de Valenciennes (DEMAV), within the Laboratoire de Matériaux Céramiques et de Mathématiques (CERAMATHS) at Université Polytechnique Hauts-de-France (Valenciennes). He is also responsible for the MAS (Modélisation, Simulation, Aléa) team of the DEMAV.

Until the end of August 2021, B. Gaudeul was the delegate of the PhD students at the Commission Mixte, and was a member of the Commission Égalité, whose aim is to fight against gender-based discriminations.

Since 2019, M. Herda is the co-organizer of the weekly of the Laboratoire Paul Painlevé. He is also an elected member of the Conseil de Laboratoire and of the Commission Mixte which, every few months, provide consultative advice and votes on matters (budget, promotions...) related to the math laboratory and the math department. M. Herda is also substitute member of the Inria Lille Comité de Centre.

Until the end of August 2021, I. Lacroix-Violet was a member of the Conseil de la Fédération de Recherche des Hauts-de-France, and of the Commission Emploi Recherche (CER) of the Inria Lille research center.

S. Lemaire is a member of the Commission de Développement Technologique (CDT) of the Inria Lille research center.

B. Merlet is in charge of the Numerical Analysis and PDEs (ANEDP) team of the Laboratoire Paul Painlevé. He is also a member of the Commission Mixte.

T. Rey is a member of the , the local correspondent of the biomath-oriented research group , a member of the council of the , and the point of contact concerning sustainable development at Laboratoire Paul Painlevé.

RAPSODI team members are strongly involved in teaching at Université de Lille (and Université Polytechnique Hauts-de-France, Valenciennes).

Faculty members of the project-team ensure their teaching duties (

Inria members of the project-team also take an important part in teaching activities. In 2021, C. Cancès taught the course "Fundamental notions in Mathematics" (32h) in the framework of the Master 1 "Data Science" of Université de Lille and École Centrale Lille, as well as a course on differential equations and their approximation (28h) for first-year students at École Centrale Lille. M. Herda gave tutorials (16h) in the introductory course to scientific computing for first-year students at École Centrale Lille. He also gave tutorials (36h) on multivariate calculus to second-year undergraduate students at Université de Lille. M. Jonval taugh a refresher course in maths (30h) at SKEMA Business School for first-year students. He also taught a course on Fourier and Laplace transforms (20h) at ISEN Lille for third-year students. S. Lemaire taught the course "Mathematical Tools for Simulation" (44h) in the framework of the Master 2 "Scientific Computing" at Université de Lille. J. Moatti gave tutorials on linear algebra (36h) for L1 MIASHS students at Université de Lille. He also gave tutorials (18h) for a refresher course in L2 LAS at Université de Lille. He finally gave tutorials (10h) in L1 SESI at Université de Lille. A. Natale taught the course "Refresher in Mathematics" (16h) in the framework of the Master 2 "Data Science" of Université de Lille and École Centrale Lille. He also taugh the course "Numerical approximation of nonlinear problems" (22h) in the framework of the Master 1 "Scientific Computing" at Université de Lille.

Post-doc of E. Facca (Inria): "Locally conservative methods for the approximation of dynamical transport on unstructured meshes", advised by A. Natale, since October 2021, funded by Inria.

Post-doc of F. Raimondi (CNRS): "Variational modeling of corrosion", co-advised by C. Cancès, C. Chainais-Hillairet and B. Merlet, since October 2020, funded by the H2020 project EURAD.

Post-doc of M. Pegon (Université de Lille): "Theoretical shape optimization problems", advised by B. Merlet, since September 2020, funded by the LabEx CEMPI.

PhD of M. Jonval (Inria/IFPEn): "Advanced numerical methods for stiff problems in the context of reactive transport", co-supervised by Ibtihel Ben Gharbia (IFPEn), C. Cancès, Thibault Faney (IFPEn) and Quang-Huy Tran (IFPEn), since October 2021, co-funded by Inria and IFPEn in the framework of the bilateral contract. Prior to his recruitment as a PhD student, M. Jonval joined the team as a research assistant between May 1st and September 30.

PhD of T. Laidin (Université de Lille): "Hybrid kinetic/fluid numerical methods and discrete hypocoercivity for the Boltzmann equation for semi-conductors", co-supervised by Marianne Bessemoulin-Chatard (CNRS & Université de Nantes), C. Chainais-Hillairet and T. Rey, since October 2021, co-funded by the LabEx CEMPI and the Hauts-de-France region.

PhD of J. Candau-Tilh (Université de Lille): "Isoperimetric problems with Wasserstein interactions", co-supervised by Michael Goldman (Université Paris Diderot) and B. Merlet, since September 2021, funded by an ENS fellowship.

PhD of J. Moatti (Inria): "Design and analysis of high-order methods for convection-diffusion models, study of the long-time behavior", co-supervised by C. Chainais-Hillairet, M. Herda and S. Lemaire, since October 2020, funded by Inria.

Post-doc of R. Bailo (Université de Lille): "Projective integration of the multiple-species Boltzmann
equation", advised by T. Rey, from June 2020 until October 2021, funded by the ERC Generator project MANAKINEQO.

PhD of S. Bassetto (IFPEn): "Towards a more robust and accurate treatment of capillary effects in multiphase flow simulations in porous media", co-supervised by C. Cancès, Guillaume Enchéry (IFPEn) and Quang-Huy Tran (IFPEn), defended on December 16 , co-funded by IFPEn and Inria in the framework of the bilateral contract.

PhD of B. Gaudeul (Université de Lille): "Numerical approximation of cross-diffusion systems arising in physics and biology", co-supervised by C. Cancès and C. Chainais-Hillairet, defended on August 30 , funded by an ENS fellowship.

Engineer position of L. Beaude (Inria) on the development of the ParaSkel++ platform, supervised
by S. Lemaire, from February 2020 until August 2021, funded by Inria (ADT ParaSkel++).

M2 internship of T. Laidin (Université de Nantes): "Hybrid kinetic/fluid numerical methods for the linear BGK equation", co-supervised by Marianne Bessemoulin-Chatard (CNRS & Université de Nantes) and T. Rey, April-September 2021.

M2 internship of A. Salheddine (Université Polytechnique Hauts-de-France, Valenciennes): "Numerical simulation of incompressible fluid flows by a BGK method", supervised by E. Creusé, July-December 2021.

M2 internship of A. Spadotto (Politecnico di Milano, Italy): "Hybrid High-Order methods for magnetostatics", co-supervised by Daniele A. Di Pietro (Université de Montpellier) and S. Lemaire, April-September 2021.

M2 project of Fabian Polvin, Chaimae Elbaraka and Ismail Bouhmala (Université Polytechnique Hauts-de-France, Valenciennes): "Implementation of the 2D MAC-FD scheme for the numerical simulation of the boundary layer equations", supervised by E. Creusé, January-May 2021.

M1 internship of L. Cleenewerck (Université de Lille): "Implementation of numerical schemes for simulations in population dynamics", supervised by M. Herda, June-July 2021.

M1 internship of S. Despierres (Université de Lille): "Projective integration for kinetic chemotaxis", supervised by T. Rey, June-July 2021.

M1 internship of J. Drappier (Université de Lille): "Finite element approximation of mean field games", supervised by A. Natale, June-July 2021.

M1 internship of K.-H. N'Gakosso (Université de Lille): "Implementation of some numerical schemes for a corrosion model", supervised by C. Chainais-Hillairet, June-July 2021.

M1 internship of F. Pigot (Université de Lille): "Micro-macro correspondence of models for pedestrian dynamics", supervised by R. Bailo, June-July 2021.

M1 internship of J. Rouzé (Université de Lille): "Inverse models for sachet salads", supervised by C. Calgaro, in collaboration with Bonduelle, June-July 2021.

M1 project of S. Despierres and F. Pigot (Université de Lille): "Asymptotic-preserving schemes for linear kinetic equations", supervised by T Rey, February-May 2021.

M1 project of Charbel Ghosn (Université de Lille): "Projective integration for nonlinear ODEs", supervised by T. Rey, January-May 2021.

L3 internship of C. Delbergue (Université de Lille): "Integration schemes for mass action kinetics", supervised by M. Jonval, October-December 2021.

C. Cancès was a jury member for the PhD defense of Aya Oussaily (Université de Technologie de Compiègne) on October 11. He was also part of the jury, as co-superviser, for the PhD defenses of B. Gaudeul (Université de Lille) on August 30, of Gabriele Todeschi (Université Paris Dauphine) on December 13, and of S. Bassetto (Université de Lille) on December 16.

C. Chainais-Hillairet reported on 3 PhD theses, respectively written by Gopikrishnan Chirappurathu Remesan (IIT Bombay, India & Monash University, Australia), Guissel Dongmo (Université Paris-Saclay), and Virgile Dubos (Sorbonne Université). She was also a jury member for the PhD defenses of Guissel Dongmo on December 13 and Virgile Dubos on December 14, and a jury member as well as the president of the jury for the PhD defenses of Hélène Bloch (Université Paris-Saclay) on September 23 and S. Bassetto (Université de Lille) on December 16. She was finally part of the jury, as co-superviser, for the PhD defense of B. Gaudeul (Université de Lille) on August 30.

E. Creusé reported on the PhD thesis of Jérémy Alloul (Université d'Orléans), defended on June 21. He was also a jury member for the PhD defense of Joanna Bisch (Université de Lille) on October 22.

M. Herda was a jury member for the PhD defense of Mohamad Rachid (Université de Nantes) on December 7.

B. Merlet reported on the PhD thesis of Chih-Kang Huang (Université Lyon 1), defended on October 15. He also reported on the Habilitation of Matthieu Bonnivard (Université Paris Diderot), to be defended on January 7, 2022.

C. Calgaro is in charge of the communication of the Laboratoire Paul Painlevé. She regularly gives/organizes conferences in high schools in the framework of the .

T. Laidin welcomed a schoolboy (3ème) for half a day to introduce him to the research world.

In October 2021, C. Chainais-Hillairet took part in the organized in Lille, where she gave a talk to an audience of high school female students. C. Chainais-Hillairet and F. Raimondi also participated to speed meeting exchanges with the schoolgirls.

In October 2021, J. Venel gave a popularization conference for an audience of high school female students during the event organized in Arras.