PRECCINSTA Burner (Vincent Moureau)
Direct Numerical Simulation of an aeronautical burner . The mesh features 2.6 billion tetrahedrons and a resolution of 100 microns.
KIAI burner (Vincent Moureau)
Large-Eddy Simulations of a swirl burner designed and operated at CORIA (J.P. Frenillot, G. Cabot, B. Renou, M. Boukhalfa).
Stratified combustion (Catherine Gruselle, Vincent Moureau, Ghislain Lartigue & Yves D'Angelo)
Large-Eddy Simulation and Direct Numerical Simulation of flame kernel development in a stratified propane/air mixture. The turbulent simulation (left movie) reproduces the experimental measurements of Balusamy S., Lecordier B. and Cessou A. from CORIA.
Two phase flow tabulated chemistry
2D Large-Eddy Simulation, injection of a premixed kerosene/air mixture on the left with a high level of turbulence. Some kerosene droplets are added to this premixing.
Two phase flow in the MERCATO burner (Benjamin Farcy)
3D simulation of the MERCATO burner under reactive conditions. Particles are two-way coupled with the gaseous phase.
Reactive flow in the MESOCORIA burner (Pierre Benard, Vincent Moureau, Ghislain Lartigue & Yves D'Angelo)
3D simulation of the MESOCORIA burner under reactive conditions: H2/CH4/air.
Formula One (David Taieb, Guillaume Ribert & Vincent Moureau)
Computation of a Formula 1 meeting with the 2010 regulations.
The design is based on the 2008 car which was simulated with the Fluent software with less than one million cells. The new car has the main features observed during the early part of F1 season, like the coca bottle shaped sidepods, the double-deck diffuser, the outer mirror disposition (forbidden by the FIA in the second part of the season), the three elements front wing.
The body of the car is discretized with 6.5mm element leading to 36 M cells in the computational domain.
Interaction between two Le Mans Series prototypes (David Taieb, Guillaume Ribert & Vincent Moureau)
NTNU wind tunnel (Félix Houtin-Mongrolle, Pierre Bénard, Ghislain Lartigue & Vincent Moureau)
This simulation represent results led in NTNU wind tunnel  by the YALES2 flow solver. This configuration present two aligned wind turbines to investigate wakes interaction. In this case the first wind turbine is yawed, i.e not aligned with the mean streamwise velocity. Moreover, the wind tunnel inflow is triggered by a upstream turbulence grid generating a sheared turbulent inflow. Such grid is modelled using an interesting method based on the actuator line method and can be found in .
Similar study for various inflows and yaw angles can be found in this paper: https://iopscience.iop.org/article/10.1088/1742-6596/1618/6/062064
NREL5MW under yaw with tower and nacelle (Félix Houtin-Mongrolle, Pierre Bénard, Ghislain Lartigue & Vincent Moureau)
Initialization of the simulation of an increasingly yawed NREL5MW (Horizontal axis Wind turbine, Diameter=126m). The inlet velocity is 8m/s. Wind turbine blades, tower and nacelle are repsresented using the actuator line method.
Volumic rendering of vorticity magnitude allow to observe the main vortices generated in the wake. Between 10s and 30s, the wind turbine is slowly reaching a yaw angle of 30°. The generated tip vortices are interacting with the tower and nacelle wake, triggering the vortices destabilization. The bottom boundary condition use a wall law model, which slowly generate a boundary layer type flow.
This case had been run on TGCC Irene Joliot-Curie supercomputer, on the Skylake partition. The cost of the simulation time (60s) is 3.7khCPU with a wall clock time of 6h on 616 cores.
DTU10MW under yaw and turbulence (Félix Houtin-Mongrolle, Pierre Bénard, Ghislain Lartigue & Vincent Moureau)
High fidelity Large Eddy Simulations of the DTU10MW (Horizontal Axis Wind Turbine, Diameter=178,3m) under two differents inflows and two yaw angles. The inlet velocity is 10m/s, for one of the inflow Mann synthetic turbulence is added to the mean. Wind turbine blades are repsresented using the actuator line method.
Q-criterion contour colored by streamwise velocity allow to observe the rotor tip vortices destabilization in the wake. For the cases with synthetic turbulence, the background vortices are in transparency. With meshes gathering up to 1.7 Billions elements, the resolution allow to capture the fine interaction of tip vortices and the wake evolution up to 12 diameter (~2.2km). The results are discussed in the following paper: https://iopscience.iop.org/article/10.1088/1742-6596/1934/1/012011
This work is part of the PRACE project WIMPY - Wind turbine Multi Physics. These cases have been run on TGCC Irene Joliot-Curie supercomputer, on the AMD Rome partition on 8448 cores.
Vertical Axis Turbine Simulation (Félix Houtin-Mongrolle, Pierre Bénard, Ghislain Lartigue & Vincent Moureau)
Initialization of the simulation of the interaction between the vortices shed by a bluff body and a vertical axis hydro turbine (VAHT, Diameter=3m). The inlet flow velocity is 5m/s.
The Strouhal number of the bluff body is synchronized with the VAHT rotation speed (TSR=3.14). The generated vortices are then transported and impact the VAHT blades profiles (DU40) when they are facing downstream.
This case had been run on TGCC Irene Joliot-Curie supercomputer, on the AMD Rome partition. The cost of one flow though time (5s) is 5.4khCPU with a wall clock time of 6.1h on 896 cores.
T7.2 Blade (Nicolas Maheu)
Large-Eddy Simulation of heat exchanges on a turbine blade.
Triple disk injector (Vincent Moureau)
Computation of a Triple Disk injector (Grout et al 2007). The densities and viscosities are those of water and air at atmospheric pressure and temperature. The video on the left was performed with 203 million tets and the one on the right with 1.6 billion tets with a resolution of 2.5 microns.
Pouring flow (Vincent Moureau and Olivier Desjardins)
Sample computation of a 2D two-phase flow with realistic properties for air and water to highlight the robustness of the method developed by Desjardins and Moureau at the 2010 CTR Summer Program.
Splashing (Vincent Moureau)
2D computation with YALES2 of a Lagrangian spray splashing on a wall and forming a film modeled with a level set and the Ghost Fluid Method. The grey particles and the grey film have the properties of water and the color represents the velocity magnitude in the gas. The Lagrangian particle are one-way coupled to the gas through drag for sake of simplicity.
Liquid jet in cross flow (Leparoux et al., ICLASS 2018)
Computation of a liquid kerosene jet in cross flow [Ragucci et al., 2007] using conservative level-set solver (SPS) coupled with parallel dynamic mesh adaptation with revisited numerical methods by [Janodet el al., JCP, 2022].
Oil churning (Cailler et al., ICMF 2019)
Simulation performed with the spray solver (SPS) of YALES2 with an accurate conservative level set (ACLS) method. The lowest part of the cylinder is immersed in an oil bath: as the cylinder rotates, the oil is entrained and then some droplets are ejected because of inertia. Moreover, some air engulfment in the oil bath can also be observed at the front part of the cylinder. Experiments from [Changenet et al., 2011] and [Leprince et al., 2012].
Lagrangian simulation of the MERCATO burner (Lola Guedot)
3D simulation of the MERCATO burner under isothermal conditions. Particles are two-way coupled with the gaseous phase. The mesh consists of 326 million tetrahedra. Velocity magnitude (top) and evaporated fuel mass fraction (bottom) are displayed in the mid-plane.
Bio-mechanics from I3M lab in Montpellier
Simulation of a cardiac cycle (Christophe Chnafa, Simon Mendez, Franck Nicoud)
3D computation of a cardiac cycle with the Arbitrary-Lagrangian Eulerian solver of YALES2. This solver and the calculations were done in the I3M lab of the University of Montpellier by C. Chnafa, S. Mendez and F. Nicoud. The color in the movie represents the vorticity.
The grid on which the fluid problem is computed is extracted from 4D (3D + time) medical images from a patient. Ten 3D images are taken from different times during the heart cycle. A grid is extracted from one medical image using a segmentation protocol. Then, grid deformations are computed from the combination of an image registration algorithm and of interpolations process. Hence, boundary movements are extracted from medical images and applied as boundary conditions for the fluid problem, resulting in a patient-specific computation. The spatial resolution is imposed to be close to 0.8 mm in all three spatial directions along the cycle, which yields grids of approximately three-million tetrahedral elements. Valves are modelled by immersed boundaries, and the heart is handled by a conformal mesh.
Settling of spherical particles (Yann Dufresne)
These results are obtained with the granular flow solver of YALES2 developed during the PhD thesis of Y. Dufresne funded by the ANR project MORE4LESS coordinated by IFP-EN. The flow solver is highly scalable and enables to perform simulations of the settling of 10 million soft spheres on 512 cores of the Curie machine (GENCI, CEA).
LES simulation of a chordwise-flexible blade (Thomas Fabbri)
Simulation of an experiment of a pitching foil in a water channel. The chord based Reynolds is approximately 200k, and the simulation results were in good agreement with experimental data. The simulation is performed using the Fluid-Structure Interaction solver developed during the PhD of Thomas Fabbri (LEGI). It uses a partitioned coupling between an Arbitrary Lagrangian Eulerian (ALE) solver and a Structural Mechanics Solver (SMS), both included inside the YALES2 library. The fluid solver uses LES approach and an unstructured grid of 27 millions tetrahedrons while 1.9k 10-nodes tetrahedrons were used by the solid solver. The coupling is performed via the CWIPI library. The mesh movement is computed using a pseudo solid method, combined with the dynamic mesh adaptation module based on MMG3D.
Turbulent flow around a flexible thin rubber plate (Thomas Fabbri)
Simulation of a Fluid-Structure Interaction case of a turbulent flow around a fixed cylinder with a flexible thin rubber plate and a rear mass (Re=30 470). The fluid mesh is composed by about 40 millions tetrehedrons while the solid mesh is made of 900 27-nodes hexahedrons.
Immersed boundaries on unstructured grids (Vincent Moureau)
On the left, 2D computation with YALES2 of the flow around two moving cylinders with an immersed boundary technique implemented for unstructured grids. The color represents the velocity magnitude. On the right, simulation of a stirred-tank reactor with YALES2. The mesh consists of 31 million tetrahedra. Simulation performed by V. Moureau from CORIA and N. Perret from Rhodia-Solvay.
Dynamic mesh adaptation (Vincent Moureau)
Demonstration of 2D and 3D dynamic mesh adaptation with YALES2. 2D remeshing is based on in-house Delaunay triangulation and 3D remeshing is based on the MMG3D library developed by C. Dobrzynski at INRIA.
LES of a Hydropower Pump Turbine (Eliott Alloin and Guillaume Balarac)
Simulations on a full PT geometry are done for several operating points. An Arbitrary-Lagrangian-Eulerian (ALE) approach coupled with a Dynamic Mesh Adaptation (DMA) methodology is used to take into account the runner rotation. Moreover, a mesh convergence strategy is used to guarantee the accuracy of the mean field discretization, and to ensure the resolution of a sufficient part of the turbulent scales.