|
|
Line 2: |
Line 2: |
| | | |
| | | |
− | == OBJECTIVE: Liquid/Gaz Interface Simulations == | + | == Scope of the code == |
− | | + | == Numerics == |
− | | + | == Post-Processsing == |
− | <div class="infobox floatright" style="width: 500px;">
| + | == Team == |
− | [[File:Imag1.png|right|thumb|300px|]]
| + | == Publications == |
− | </div>
| + | == Partnership == |
− | | + | |
− | | + | |
− | {| border="0" cellspacing="0" cellpadding="2" align="left"
| + | |
− | ! Describe the interface motion precisely
| + | |
− | ! [[File:fleche_vert.png|50 px]]
| + | |
− | ! Level Set Method
| + | |
− | |}
| + | |
− | | + | |
− | | + | |
− | | + | |
− | {| border="0" cellspacing="0" cellpadding="2" align="left"
| + | |
− | ! Handle jump conditions at the interface
| + | |
− | ! without artificial smoothing
| + | |
− | ! [[File:fleche_vert.png|50 px]]
| + | |
− | ! Ghost Fluid Method
| + | |
− | |-
| + | |
− | |}
| + | |
− | | + | |
− | | + | |
− | | + | |
− | {| border="0" cellspacing="0" cellpadding="2" align="left"
| + | |
− | ! Respect mass conservation
| + | |
− | ! [[File:fleche_vert.png|50 px]]
| + | |
− | ! VOF Method
| + | |
− | |-
| + | |
− | |}
| + | |
− | | + | |
− | | + | |
− | | + | |
− | {| border="0" cellspacing="0" cellpadding="2" align="left"
| + | |
− | ! Solve incompressible Navier Stokes equations
| + | |
− | ! [[File:fleche_vert.png|50 px]]
| + | |
− | ! Projection Method
| + | |
− | |-
| + | |
− | |}
| + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | [[File:fleche_roueg.png|50 px]] '''Archer is a 3D incompressible Navier Stokes solver with Level Set-Ghost Fluid-VOF coupling and MPI parallelization'''
| + | |
− | | + | |
− | | + | |
− | | + | |
− | == NUMERICS ==
| + | |
− | | + | |
− | * Cartesian mesh
| + | |
− | * Mac grid (velocities on cell boundaries)
| + | |
− | * WENO 5 scheme for convective terms, Adams Bashforth or RK3 or split algorithm in time
| + | |
− | * Multigrid algorith for preconditionning Conjugate Gradient Method in Poisson equation solver
| + | |
− | * Ghost Fluid method for variable discontinuities at the interface
| + | |
− | * CLSVOF for mass conservation
| + | |
− | | + | |
− | *MPI parallelization
| + | |
− | | + | |
− | | + | |
− | | + | |
− | == UNDER DEVELOPMENT ==
| + | |
− | | + | |
− | * Immersed boundaries
| + | |
− | * Adaptive mesh refinement
| + | |
− | * Coupling with Lagrangian solver
| + | |
− | * '''New convection scheme (modified Rudman's method)'''
| + | |
− | | + | |
− | == GALLERY ==
| + | |
− | | + | |
− | | + | |
− | | + | |
− | === Turbulent jet ===
| + | |
− | <div class="infobox floatright" style="width: 500px;">
| + | |
− | [[File:jet_2048.jpg|right|thumb|300px|]]
| + | |
− | </div>
| + | |
− | | + | |
− | | + | |
− | | + | |
− | Diesel type Jet Atomization :
| + | |
− | Diameter: 100µm, Liquid velocity 100 ms-1, Turbulence 5%, gas velocity: 0 ms-1,
| + | |
− | <math>\rho_{liq}</math>=696 kgm-3, <math>\rho_{gas}</math>=25 kgm-3
| + | |
− | | + | |
− | <math>\mu_{liq}</math> =1.210-3kgm-1s-1, <math>\mu_{gas}</math>=10-5kgm-1s-1, <math>\sigma</math>=0.06Nm-1
| + | |
− | Numerical simulation by coupling Level Set / VOF / Ghost fluid methods
| + | |
− | | + | |
− | Mesh : 256x256x2048 (130 millions points) MPI parallelization 128 procs.
| + | |
− | *[[Media:Jet_Fin.avi|Jet.avi]]
| + | |
− | | + | |
− | | + | |
− | === Triple disk injector ===
| + | |
− | <div class="infobox floatright" style="width: 500px;">
| + | |
− | {| class="wikitable" style="margin: 1em auto 1em auto;"
| + | |
− | |[[File:Mesh_intern_LES.png|center|thumb|Mesh for internal flow LES simulations|150px|]]
| + | |
− | |[[File:Velocity.png|center|thumb|150px|]]
| + | |
− | |[[File:triple.png|center|thumb|150px|]]
| + | |
− | |}
| + | |
− | | + | |
− | | + | |
− | </div>
| + | |
− | *Coupling between internal flow simulations and DNS of primary breakup
| + | |
− | *Radius of the jet exit: 90 microns
| + | |
− | *Mesh 256x1024x1024 (~270 millions) dx~1.44 micron
| + | |
− | *~3000 nodes on the (half) jet exit
| + | |
− | *2048 procs
| + | |
− | *[[Media:toto.avi|Triple.avi]]
| + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | === Liquid/Gas mixing layer ===
| + | |
− | <div class="infobox floatright" style="width: 680px;">
| + | |
− | {| class="wikitable" style="margin: 1em auto 1em auto;"
| + | |
− | |[[File:Couche1.png|center|thumb||100px|]]
| + | |
− | |[[File:couche2.png|center|thumb|150px|]]
| + | |
− | |[[File:nappe.png|center|thumb|150px|]]
| + | |
− | |}
| + | |
− | | + | |
− | </div>
| + | |
− | *Air/water mixing layer
| + | |
− | *30m/s in air, 0.3 in water
| + | |
− | *Grid is 512x512x1024 nodes
| + | |
− | *Mesh size equals 48 µm.
| + | |
− | *2048 procs, 20 hours for 3500 time steps; Simulation on 70000 time steps, dt~=10-6 s
| + | |
− | | + | |
− | *[[Media:shear.avi|Shear.avi]]
| + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | === Air assisted atomization === | + | |
− | <div class="infobox floatright" style="width: 500px;">
| + | |
− | [[File:0193.jpeg|right|thumb|150px|]]
| + | |
− | </div>
| + | |
− | | + | |
− | *Air/water Jet « Marmottant » Modified Rudman's method two grids
| + | |
− | *Ugas=35m/s and Uliq=0.2m/s
| + | |
− | *<math>\rho_{gas}</math>=1,2 kg/m3 <math>\rho_{liq}</math> =1000 kg/m3
| + | |
− | *Lx=45.6mm (4mm+1.7mm)x8 that is 4 the « jet +gas » diameter
| + | |
− | *Grid: 256x256x256,(512x512x512 for level set and VOF, mesh is 90µm)
| + | |
− | *1024 procs on Curie, 12h wall-clock for 5000 time steps.
| + | |
− | *Simulation is 250000 time steps.
| + | |
− | *[[Media:Anim1 294 crop.avi |Marmotte.avi]]
| + | |
− | | + | |
− | === 2D gas/liquid shear layer ===
| + | |
− | <div class="infobox floatright" style="width: 500px;">
| + | |
− | [[File:9_11.jpeg|right|thumb|300px|]]
| + | |
− | </div>
| + | |
− | *M=<math>\rho_{gas}</math><math>U_{gas}</math>^2 / <math>\rho_{liq}</math><math>U_{liq}</math>^2 = 16
| + | |
− | *Ugas=30 m/s Uliq=0.26 m/s
| + | |
− | *<math>\rho_{gas}</math>=1,2 kg/m3 <math>\rho_{liq}</math> =1000 kg/m3
| + | |
− | *Lx=80mm, Ly=40 mm 10mm water + 10mm air + 20mm quiet
| + | |
− | *Grid: 1024x512(that is 2048x1024 level set and VOF, mesh is 40 µm)
| + | |
− | *512 procs on Curie
| + | |
− | *18h wall clock for 100000 time steps.
| + | |
− | *Simulation on 2 000 000 time steps
| + | |
− | *[[Media:Ite1_20_1sur10_25im.avi|2Dliquid_shear_layer.avi]]
| + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | | + | |
− | === Air assisted atomization === | + | |
− | <div class="infobox floatright" style="width: 500px;">
| + | |
− | [[File:ite304.jpg|right|thumb|150px|]]
| + | |
− | </div>
| + | |
− | | + | |
− | *Air/water Jet « Marmottant » RUDMAN-TYPE Technique with one grid
| + | |
− | *Ugas=22.6m/s and Uliq=0.27m/s
| + | |
− | *<math>\rho_{gas}</math>=1,2 kg/m3 <math>\rho_{liq}</math> =1000 kg/m3
| + | |
− | *Lx=34.2mm (4mm+1.7mm)x6 that is 3 x the « jet +gas » diameter
| + | |
− | *Grid: 256x256x256,(mesh is 134µm)
| + | |
− | *512 procs on Curie, 12h wall-clock for 30000 time steps.
| + | |
− | *Grid: 512x512x512
| + | |
− | *Grid: 512x512x1024
| + | |
− | *8192 procs on TURING, 16h wall clock for 6000 time steps.
| + | |
− | *[[Media:TOUT_comp2.avi |Marmotte.avi]]
| + | |
− | | + | |
− | === Air assisted atomization ===
| + | |
− | <div class="infobox floatright" style="width: 500px;">
| + | |
− | [[File:im1000.jpg|right|thumb|150px|]]
| + | |
− | </div>
| + | |
− | *Same configuration
| + | |
− | *Initial boundary layers included
| + | |
− | *Grid:512x512x1024
| + | |
− | *[[Media:With_bound_lay.avi |Oct_15_2015.avi]]
| + | |
− | *Helicoïdal behaviourin simulations (experimentally observed)
| + | |
− | *[[Media:Tourne960.avi |helicoid.avi]]
| + | |