Difference between revisions of "YALES2 Gallery"

From www.coria-cfd.fr
Jump to: navigation, search
(Oil churning (Cailler et al., ICMF 2019))
(33 intermediate revisions by 6 users not shown)
Line 23: Line 23:
 
*[[#Two-phase flow tabulated chemistry|Two-phase flow tabulated chemistry]]
 
*[[#Two-phase flow tabulated chemistry|Two-phase flow tabulated chemistry]]
 
*[[#MERCATO burner|MERCATO burner]]
 
*[[#MERCATO burner|MERCATO burner]]
 +
*[[#MESOCORIA burner|MESOCORIA burner]]
 
}}
 
}}
  
Line 35: Line 36:
 
*[[#Splashing|Splashing]]
 
*[[#Splashing|Splashing]]
 
*[[#Isothermal flow in the MERCATO burner|Isothermal flow in the MERCATO burner]]
 
*[[#Isothermal flow in the MERCATO burner|Isothermal flow in the MERCATO burner]]
 +
}}
 +
 +
{{Main Page/Frame
 +
| color      = DEB887
 +
| title      = Granular flows
 +
| content    =
 +
DEM (Discrete Element Method) simulations of granular flows
 +
 +
*[[#Settling of spherical particles|Settling of spherical particles]]
 
}}
 
}}
  
Line 43: Line 53:
  
 
{{Main Page/Frame
 
{{Main Page/Frame
  | color      = 990000
+
  | color      = 87CEFA
 
  | title      = Aerodynamics
 
  | title      = Aerodynamics
 
  | content    =
 
  | content    =
Line 50: Line 60:
 
*[[#Formula One|Formula One]]
 
*[[#Formula One|Formula One]]
 
*[[#Le Mans Series prototypes|Le Mans Series prototypes]]
 
*[[#Le Mans Series prototypes|Le Mans Series prototypes]]
 +
 +
Large-Eddy Simulation of wind turbines wakes
 +
*[[#NTNU wind tunnel|NTNU wind tunnel]]
 +
*[[#NREL5MW under yaw with tower and nacelle|NREL5MW under yaw with tower and nacelle]]
 +
*[[#DTU10MW under yaw and turbulence|DTU10MW under yaw and turbulence]]
 +
*[[#Vertical Axis Turbine Simulation|Vertical Axis Turbine Simulation]]
 
}}
 
}}
  
 
{{Main Page/Frame
 
{{Main Page/Frame
  | color      = 009999
+
  | color      = FFD700
 
  | title      = Heat transfers
 
  | title      = Heat transfers
 
  | content    =
 
  | content    =
Line 77: Line 93:
  
 
*[[#Immersed boundaries on unstructured grids|Immersed boundaries on unstructured grids]]
 
*[[#Immersed boundaries on unstructured grids|Immersed boundaries on unstructured grids]]
*[[#Mesh deformation|Mesh deformation]]
+
*[[#Dynamic mesh adaptation|Dynamic mesh adaptation]]
 
}}
 
}}
  
Line 106: Line 122:
  
 
<span id="KIAI burner"></span>
 
<span id="KIAI burner"></span>
 +
 
=== '''KIAI burner''' ([[User:Moureauv|Vincent Moureau]])===
 
=== '''KIAI burner''' ([[User:Moureauv|Vincent Moureau]])===
 
Large-Eddy Simulations of a swirl burner designed and operated at CORIA (J.P. Frenillot, G. Cabot, B. Renou, M. Boukhalfa).
 
Large-Eddy Simulations of a swirl burner designed and operated at CORIA (J.P. Frenillot, G. Cabot, B. Renou, M. Boukhalfa).
Line 113: Line 130:
 
| [[File:KIAI_382M_U.png|center|thumb|Velocity field for the cold flow - 382M tetrahedrons|350px]]
 
| [[File:KIAI_382M_U.png|center|thumb|Velocity field for the cold flow - 382M tetrahedrons|350px]]
 
| [[File:KIAI_382M_Q.png|center|thumb|Q-criterion for the cold flow - 382M tetrahedrons|350px]]
 
| [[File:KIAI_382M_Q.png|center|thumb|Q-criterion for the cold flow - 382M tetrahedrons|350px]]
 +
| colspan="3" |
 +
{| style="margin: 8px|center;"
 +
| {{#widget:YouTube|id=5SHPYYSow6U|width=500|height=350}}
 +
|}
 
|}
 
|}
  
  
 +
<span id="Stratified combustion"></span>
  
=== '''Stratified combustion''' ([[User:Gruselle|Catherine Gruselle]], [[User:Moureauv|Vincent Moureau]] and [[User:Lartigue|Ghislain Lartigue]])===
+
=== '''Stratified combustion''' ([[User:Gruselle|Catherine Gruselle]], [[User:Moureauv|Vincent Moureau]], [[User:Lartigue|Ghislain Lartigue]] & [[User:Dangelo| Yves D'Angelo]])===
 
Large-Eddy Simulation and Direct Numerical Simulation of flame kernel development in a stratified propane/air mixture.
 
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.
 
The turbulent simulation (left movie) reproduces the experimental measurements of Balusamy S., Lecordier B. and Cessou A. from CORIA.
Line 127: Line 149:
 
| {{#widget:YouTube|id=LdKXaX4d5Uw|width=400|height=300}}
 
| {{#widget:YouTube|id=LdKXaX4d5Uw|width=400|height=300}}
 
|}
 
|}
 +
 +
 +
<span id="Two-phase flow tabulated chemistry"></span>
  
 
=== '''Two phase flow tabulated chemistry''' ===
 
=== '''Two phase flow tabulated chemistry''' ===
Line 138: Line 163:
 
|}
 
|}
  
=== '''Two phase flow in the MERCATO burner''' ([[User:Farcyb|Benjamin Farcy]]) ===
+
 
 +
<span id="MERCATO burner"></span>
 +
=== '''Two phase flow in the MERCATO burner''' ([[User:Farcyb|Benjamin Farcy]])===
  
 
3D simulation of the MERCATO burner under reactive conditions. Particles are two-way coupled with the gaseous phase.  
 
3D simulation of the MERCATO burner under reactive conditions. Particles are two-way coupled with the gaseous phase.  
Line 146: Line 173:
 
|-
 
|-
 
| [[File:blue_flame.png|800px]]
 
| [[File:blue_flame.png|800px]]
 +
|}
 +
 +
 +
<span id="MESOCORIA burner"></span>
 +
=== '''Reactive flow in the MESOCORIA burner''' ([[User:Benard|Pierre Benard]], [[User:Moureauv|Vincent Moureau]], [[User:Lartigue|Ghislain Lartigue]] & [[User:Dangelo| Yves D'Angelo]]) ===
 +
 +
3D simulation of the MESOCORIA burner under reactive conditions: H2/CH4/air.   
 +
 +
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
|+ MESO-CORIA burner with YALES2
 +
|-
 +
| {{#widget:YouTube|id=KiNwKE2t7v0|width=400|height=300}}
 +
| {{#widget:YouTube|id=gey2Dv-WLg4|width=400|height=300}}
 
|}
 
|}
  
 
== Aerodynamics ==
 
== Aerodynamics ==
  
 +
<span id="Formula One"></span>
 
=== '''Formula One''' ([[User:Taieb|David Taieb]], [[User:Ribert|Guillaume Ribert]] & [[User:Moureauv|Vincent Moureau]]) ===
 
=== '''Formula One''' ([[User:Taieb|David Taieb]], [[User:Ribert|Guillaume Ribert]] & [[User:Moureauv|Vincent Moureau]]) ===
  
Line 169: Line 211:
 
|}
 
|}
  
 +
 +
<span id="Le Mans Series prototypes"></span>
 
=== '''Interaction between two Le Mans Series prototypes''' ([[User:Taieb|David Taieb]], [[User:Ribert|Guillaume Ribert]] & [[User:Moureauv|Vincent Moureau]]) ===
 
=== '''Interaction between two Le Mans Series prototypes''' ([[User:Taieb|David Taieb]], [[User:Ribert|Guillaume Ribert]] & [[User:Moureauv|Vincent Moureau]]) ===
 
{| class="wikitable" style="margin: 1em auto 1em auto;"
 
{| class="wikitable" style="margin: 1em auto 1em auto;"
Line 178: Line 222:
 
| [[File:LMS_stream_Umean.jpg|center|Streamlines of averaged velocity colored by velocity RMS.|400px]]
 
| [[File:LMS_stream_Umean.jpg|center|Streamlines of averaged velocity colored by velocity RMS.|400px]]
 
| [[File:LMS_wake_DF.jpg|center|Longitudinal slice of instantaneous velocity and downforce on bodies.|400px]]
 
| [[File:LMS_wake_DF.jpg|center|Longitudinal slice of instantaneous velocity and downforce on bodies.|400px]]
 +
|}
 +
 +
<span id="NTNU wind tunnel"></span>
 +
=== '''NTNU wind tunnel''' ([[User:Houtin|Félix Houtin-Mongrolle]], [[User:Benard|Pierre Bénard]], [[User:Lartigue|Ghislain Lartigue]] & [[User:Moureauv|Vincent Moureau]]) ===
 +
 +
This simulation represent results led in NTNU wind tunnel [1] 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 [2].
 +
 +
[1] https://doi.org/10.5194/wes-3-883-2018
 +
[2] https://doi.org/10.1080/14685248.2020.1803495
 +
 +
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
 +
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
|+  LES of two aligned wind-turbine in a wind tunnel
 +
|-
 +
|
 +
{| style="margin: 10px;"
 +
|{{#widget:YouTube|id=CWpm8Y1bK2c|width=700|height=400}}
 +
|}
 +
|}
 +
 +
<span id="NREL5MW under yaw with tower and nacelle"></span>
 +
=== '''NREL5MW under yaw with tower and nacelle''' ([[User:Houtin|Félix Houtin-Mongrolle]], [[User:Benard|Pierre Bénard]], [[User:Lartigue|Ghislain Lartigue]] & [[User:Moureauv|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.
 +
 +
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
|+  LES of yawed NREL5MW wind turbine with tower and nacelle
 +
|-
 +
|
 +
{| style="margin: 10px;"
 +
|{{#widget:YouTube|id=FScq40rfBp0|width=630|height=280}}
 +
|}
 +
|}
 +
 +
<span id="DTU10MW under yaw and turbulence"></span>
 +
=== '''DTU10MW under yaw and turbulence''' ([[User:Houtin|Félix Houtin-Mongrolle]], [[User:Benard|Pierre Bénard]], [[User:Lartigue|Ghislain Lartigue]] & [[User:Moureauv|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.
 +
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
|+  High fidelity LES of yawed DTU10MW wind turbine under turbulent/non-turbulent inflows
 +
|-
 +
|
 +
{| style="margin: 10px;"
 +
|{{#widget:YouTube|id=ShVmN-ZpyVE|width=580|height=300}}
 +
|}
 +
|}
 +
 +
<span id="Vertical Axis Turbine Simulation"></span>
 +
=== '''Vertical Axis Turbine Simulation''' ([[User:Houtin|Félix Houtin-Mongrolle]], [[User:Benard|Pierre Bénard]], [[User:Lartigue|Ghislain Lartigue]] & [[User:Moureauv|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.
 +
 +
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
|+  Vertical Axis Turbine Simulation
 +
|-
 +
|
 +
{| style="margin: 10px;"
 +
|{{#widget:YouTube|id=vKTHUtx2sxc|width=580|height=280}}
 +
|}
 
|}
 
|}
  
 
== Heat transfers ==
 
== Heat transfers ==
  
 +
<span id="T7.2 blade"></span>
 
=== '''T7.2 Blade''' ([[User:Maheu|Nicolas Maheu]])===
 
=== '''T7.2 Blade''' ([[User:Maheu|Nicolas Maheu]])===
 
Large-Eddy Simulation of heat exchanges on a turbine blade.
 
Large-Eddy Simulation of heat exchanges on a turbine blade.
Line 193: Line 323:
 
| {{#widget:YouTube|id=iZWYfN4vDrQ|width=400|height=300}}
 
| {{#widget:YouTube|id=iZWYfN4vDrQ|width=400|height=300}}
 
|}
 
|}
 +
  
 
== Two-phase flows ==
 
== Two-phase flows ==
  
 +
<span id="Triple Disk Injector"></span>
 
=== '''Triple disk injector''' ([[User:Moureauv|Vincent Moureau]]) ===
 
=== '''Triple disk injector''' ([[User:Moureauv|Vincent Moureau]]) ===
  
Line 210: Line 342:
 
|}
 
|}
  
 +
 +
<span id="Pouring flow"></span>
 
=== '''Pouring flow''' ([[User:Moureauv|Vincent Moureau]] and [http://cmes.colorado.edu/ Olivier Desjardins]) ===
 
=== '''Pouring flow''' ([[User:Moureauv|Vincent Moureau]] and [http://cmes.colorado.edu/ Olivier Desjardins]) ===
  
Line 223: Line 357:
 
|}
 
|}
  
 +
 +
<span id="Splashing"></span>
 
=== '''Splashing''' ([[User:Moureauv|Vincent Moureau]]) ===
 
=== '''Splashing''' ([[User:Moureauv|Vincent Moureau]]) ===
  
Line 235: Line 371:
 
|}
 
|}
 
|}
 
|}
 +
 +
<span id="JICF"></span>
 +
=== '''Liquid jet in cross flow''' ([https://www.researchgate.net/publication/326631733_Primary_atomization_simulation_applied_to_a_jet_in_crossflow_aeronautical_injector_with_dynamic_mesh_adaptation 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].
 +
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
|+ Liquid jet in cross flow with YALES2
 +
|-
 +
|
 +
{| style="margin: 10px;"
 +
|{{#widget:Video|url=https://www.coria-cfd.fr/images/1/10/Ragucci.mp4|width=1024|height=600}}
 +
|}
 +
|}
 +
 +
<span id="Churning"></span>
 +
=== '''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].
 +
 +
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
|+ Oil churning with YALES2 (left: 2D slice colored by cell size; Right: 3D liquid-gaz interface)
 +
|-
 +
|
 +
{| style="margin: 10px;"
 +
|{{#widget:Video|url=https://www.coria-cfd.fr/images/d/d8/Churning_Metric.mp4|width=512|height=300}}{{#widget:Video|url=https://www.coria-cfd.fr/images/c/c7/Churning_3D.mp4|width=512|height=300}}
 +
|}
 +
|}
 +
 +
 +
<span id="Isothermal flow in the MERCATO burner"></span>
  
 
=== '''Lagrangian simulation of the MERCATO burner''' ([[User:Guedot|Lola Guedot]]) ===
 
=== '''Lagrangian simulation of the MERCATO burner''' ([[User:Guedot|Lola Guedot]]) ===
Line 246: Line 415:
 
== Bio-mechanics from  [http://ens.math.univ-montp2.fr/ I3M lab in Montpellier] ==
 
== Bio-mechanics from  [http://ens.math.univ-montp2.fr/ I3M lab in Montpellier] ==
  
 +
<span id="Simulation of a cardiac cycle"></span>
 
=== '''Simulation of a cardiac cycle''' ([[User:Chnafa|Christophe Chnafa]], [[User:Mendez|Simon Mendez]], [[User:Nicoud|Franck Nicoud]]) ===
 
=== '''Simulation of a cardiac cycle''' ([[User:Chnafa|Christophe Chnafa]], [[User:Mendez|Simon Mendez]], [[User:Nicoud|Franck Nicoud]]) ===
 
 
  
 
{| class="wikitable" style="margin: 1em auto 1em auto;"
 
{| class="wikitable" style="margin: 1em auto 1em auto;"
Line 262: Line 430:
 
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 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.
 
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.
 +
 +
 +
== Granular flows ==
 +
 +
<span id="Settling of spherical particles"></span>
 +
=== '''Settling of spherical particles''' ([[User:Ydufresne|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).
 +
 +
{| class="wikitable" style="margin: 1em auto 1em auto;"
 +
|+Granular flow solver of YALES2
 +
|-
 +
|
 +
{| style="margin: 10px;"
 +
|{{#widget:YouTube|id=RddU7d-0Hyw|width=400|height=300}}
 +
|{{#widget:YouTube|id=3XMatY-lM6c|width=400|height=300}}
 +
|}
 +
|}
 +
  
 
== Advanced numerics ==
 
== Advanced numerics ==
  
 +
<span id="Immersed boundaries on unstructured grids"></span>
 
=== '''Immersed boundaries on unstructured grids''' ([[User:Moureauv|Vincent Moureau]]) ===
 
=== '''Immersed boundaries on unstructured grids''' ([[User:Moureauv|Vincent Moureau]]) ===
  
Line 279: Line 467:
 
|}
 
|}
  
=== '''Mesh deformation''' ([[User:Moureauv|Vincent Moureau]]) ===
 
  
Demonstration of 2D mesh deformation with YALES2. Only the velocity of boundaries is prescribed and the movement of the nodes is found by inverting an elliptic system. Edge swapping is also activated in this example.
+
<span id="Dynamic mesh adaptation"></span>
 +
=== '''Dynamic mesh adaptation''' ([[User:Moureauv|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.
  
 
{| class="wikitable" style="margin: 1em auto 1em auto;"
 
{| class="wikitable" style="margin: 1em auto 1em auto;"
|+ Mesh deformation with YALES2
+
|+ Dynamic mesh adaptation with YALES2
 
|-
 
|-
 
|
 
|
 
{| style="margin: 10px;"
 
{| style="margin: 10px;"
 
|{{#widget:YouTube|id=riJM_NOeA_M|width=400|height=300}}
 
|{{#widget:YouTube|id=riJM_NOeA_M|width=400|height=300}}
 +
|{{#widget:YouTube|id=5elSG_CxF6M|width=400|height=300}}
 +
|{{#widget:YouTube|id=Eaw3g-l2HbY|width=400|height=300}}
 
|}
 
|}
 
|}
 
|}

Revision as of 08:19, 1 June 2022


Y2header bis.png

Welcome to the YALES2 gallery

Selected images and videos of high-fidelity simulations

  Combustion
  Two-phase flows
Two-phase flow simulations with the spray solver (Conservative Level Set + Ghost-Fluid Method) and with the Lagrangian spray solver
  Granular flows
DEM (Discrete Element Method) simulations of granular flows
  Aerodynamics
Large-Eddy Simulation of aerodynamics of complex bodies

Large-Eddy Simulation of wind turbines wakes

  Heat transfers
Large-Eddy Simulation of aerodynamics of complex bodies
  Biomechanics
Large-Eddy Simulation of aerodynamics of complex bodies
  Advanced numerics
Large-Eddy Simulation of aerodynamics of complex bodies

Combustion

PRECCINSTA Burner (Vincent Moureau)

Direct Numerical Simulation of an aeronautical burner [1]. The mesh features 2.6 billion tetrahedrons and a resolution of 100 microns.

PRECCINSTA burner with YALES2
Iso-surface of the Q criterion for the isothermal case
Temperature field for the fully premixed reacting case
OH radical field for the fully premixed reacting case
Couverture du Numéro Spécial Calcul Intensif des Comptes Rendus de Mécanique de l'académie des sciences


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).

KIAI burner with YALES2
Velocity field for the cold flow - 382M tetrahedrons
Q-criterion for the cold flow - 382M tetrahedrons


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.

Stratified combustion with YALES2


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 combustion with YALES2


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.

MERCATO burner with YALES2
Blue flame.png


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.


MESO-CORIA burner with YALES2

Aerodynamics

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.

Formula One with YALES2
Formula 1 with 36 Million cells - Streamlines
Formula 1 with 36 Million cells - Iso-Q criterion


Interaction between two Le Mans Series prototypes (David Taieb, Guillaume Ribert & Vincent Moureau)

Interaction between two Le Mans Series prototypes with YALES2
Instantaneous streamlines colored by velocity RMS.
centerContour of pressure on the upper bodywork.
Streamlines of averaged velocity colored by velocity RMS.
Longitudinal slice of instantaneous velocity and downforce on bodies.

NTNU wind tunnel (Félix Houtin-Mongrolle, Pierre Bénard, Ghislain Lartigue & Vincent Moureau)

This simulation represent results led in NTNU wind tunnel [1] 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 [2].

[1] https://doi.org/10.5194/wes-3-883-2018 [2] https://doi.org/10.1080/14685248.2020.1803495

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

LES of two aligned wind-turbine in a wind tunnel

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.


LES of yawed NREL5MW wind turbine with tower and nacelle

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.

High fidelity LES of yawed DTU10MW wind turbine under turbulent/non-turbulent inflows

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.


Vertical Axis Turbine Simulation

Heat transfers

T7.2 Blade (Nicolas Maheu)

Large-Eddy Simulation of heat exchanges on a turbine blade.

T7.2 blade with YALES2
T7.2 Blade - Iso-Q criterion - 240M tetrahedrons
T7.2 Blade - Iso-T 325K - 240M tetrahedrons


Two-phase flows

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.

Primary atomization with YALES2


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.

Primary atomization with YALES2


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.

Wall splashing with YALES2

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].

Liquid jet in cross flow with YALES2

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].


Oil churning with YALES2 (left: 2D slice colored by cell size; Right: 3D liquid-gaz interface)


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.

MERCATO burner with YALES2
Belle image 1.png

Bio-mechanics from I3M lab in Montpellier

Simulation of a cardiac cycle (Christophe Chnafa, Simon Mendez, Franck Nicoud)

Cardiac cycle with YALES2

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.


Granular flows

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).

Granular flow solver of YALES2


Advanced numerics

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.

Immersed boundaries with YALES2


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.

Dynamic mesh adaptation with YALES2