The CORIA lab shared ressources are:

Metrological ressources from RENADIAG


RENADIAG TOFU (tracking of ultrafast phenomena) is an innovative experimental platform dedicated to the use of ultrafast laser sources for the detection and tracking of ultrafast or rare events. This optical metrology platform, funded by the Normandy region and the European Union, relies on CORIA's know-how in the design and use of laser sources, notably based on optical fibres, at different wavelengths in the near (NIR) and medium infrared (MIR). Through this platform, the laboratory has unique means for the preparation of specialty optical fibres and for the realisation of original ultrafast imaging techniques based on these fibres as well as for their characterisation.

RENADIAG MDA (Aerosol Metrology)

The laboratory has, in the optics and lasers department, an advanced experimental platform for the generation and characterization of soot and other aerosols The size distribution of nanoscale aerosols can be accessed by ex-situ (SMPS) or in-situ (angular scattering) measurements. Mass or volume concentrations are also accessible after sampling (PPS Pegasor) or in-situ (Laser Induced Incandescence). It is also possible to take samples of particles by controlling the dilution and perform morphological analyses of them on the basis of electron microscopy (TEM) images. Classification devices by mass (CPMA), mobility diameter (DMA), aerodynamic diameter (AAC) are also available.

Numerical facilities


The YALES2 code is a massively parallel numerical toolbox on which many solvers have been developed since 2009 to address various physical problems: incompressible flows, combustion, two-phase flows, heat transfer, radiative transfer and much more. All these solvers benefit from the underlying high-performance numerical library: highly-optimized linear solvers, automatic mesh refinement, high-order Finite-Volume integration, parallel I/O, etc.

YALES2 relies on unstructured meshes which enable to describe easily and efficiently complex geometries and a fully parallel dynamic mesh adaptation technique to improve the resolution in physically-relevant zones to mitigate the computational cost. As a result it can easily handle meshes composed of billions of tetrahedra, thus enabling the Direct Numerical Simulation of laboratory and semi-industrial configurations.

YALES2 is currently used by more than 300 persons both in academic and industrial sectors. The library is composed of approximately 500'000 lines of Object-Oriented Fortran and its parallelism is currently ensured by a pure MPI paradigm, although a hybrid OpenMP/MPI as well as a GPU version are under development. From a numerical point-of-view, YALES2 relies on a Node-Centered Finite-Volume approach which naturally guarantees the conservation of the transported quantities. It features various numerical methods, and more particularly a 4th-order convection scheme which is particularly well suited for the Large-Eddy Simulation (LES) approach, as it can transport the eddies present in the flow over long distances without damping them artificially.

The central part of the code is the Low-Mach number Navier-Stokes solver for constant density flows which solves the following set of equations:

To do so, it relies on the classical projection method proposed by Chorin (1968) to impose the divergence-free constraint on the velocity, augmented with a recycling of the pressure at the previous timestep to reduce the splitting errors. This method is composed of a first-step known as the prediction:

followed by a second step known as the correction:

In this last step, the pressure variation is obtained by solving a Poisson equation to enforce the divergence-free constraint on the velocity:

The discretization of this Poisson equation on a mesh will then provide a linear system in which the unknowns are the pressure at each node of the mesh. This linear system can thus be extremely large for real-life applications (the standard mesh size for 3D problems is now between 10M and 500M cells) and can only be solved with dedicated methods, fitted for massively parallel environment. Even with tens of years of R&D on linear solvers, solving the Poisson equation in turbulent flows still accounts for most of the computation cost. Moreover, solving this linear system becomes harder and harder when the mesh size and the number of CPU increases. The first reason is linked to the condition number of the linear system which is actually a measure of the ratio between the largest and the smallest scale in the problem.

The other issue is related to the need to perform parallel operation to solve the linear system. In YALES2, this system is solved with an Conjugate Gradient iterative solver coupled to a preconditioning method based on a deflation algorithm, which first solves the system on a coarser grid and then perform a back-projection of this coarse residual on the fine mesh before performing a classical CG step. This coarse mesh is obtained via a double domain decomposition technique: the full domain is first shared among all the available cores and then each core performs another domain decomposition on its own mesh.


ARCHER is a High-Performance-Computing code developed at the CORIA laboratory. It stands for Academic Research Code for Hydrodynamic Equations Resolution. It also takes its name from the fish, the archerfish (Toxotidae), known for preying insects by spitting a jet of water.

ARCHER is aimed at carrying Direct Numerical Simulations of two-phase flows may they be turbulent, incompressible or compressible, with phase change or in presence of solid boundaries. Data from Archer are used for probing the physical properties (either geometrical, morphological, topological, or dynamical) of different phenomena such as atomization, spray formation, dispersion, evaporation, phase separation, capture of solid aerosols. The close connection of the Archer developers with experimentalist is further noticeable. This wide knowledge serves for building or reinforcing physics-informed models, notably the Eulerian Lagragian Spray Atomization model (ELSA).

It was one of the first code worldwide, undertaking the simulation of liquid-jet atomization under a realistic injection configuration.

Responsable: Alexandre Poux
Scientific project: Thibaut Ménard:
The ARCHER code


The Monte Carlo Aggregation Code (MCAC) is a Discrete Element Method tool developed to simulate the aggregation of soot particles in flames.
These aggregates are initially modeled as spherical primary particles. When these particles grow in flames by coagulation and surface reactions, they experience a simultaneous change in flow and aggregation regimes. This is taken into account in MCAC simulations. This code also considers the mass transfer with the gas phase and also the soot aggregates breakage by oxidation. This is accurately done by coupling MCAC with a Computational Fluid Dynamics (CFD) code. In this context, a population of fractal-like aggregates consisting of overlapped spherical primary particles are simulated.
These particles are in remarkable resemblance with those observed experimentally for example by Transmission Electron Microscopy (TEM) image analysis. Finally, this code also allows both primary/aggregate size distribution, and kinetics of aggregation to be studied.