Difference between revisions of "Archer"
From www.coria-cfd.fr
Aniszewski (Talk | contribs) (creation of a page) |
|||
(95 intermediate revisions by 2 users not shown) | |||
Line 1: | Line 1: | ||
− | '''Archer''' is | + | [[File:Archer.png]] |
+ | |||
+ | |||
+ | == OBJECTIVE: Liquid/Gaz Interface Simulations == | ||
+ | |||
+ | |||
+ | <div class="infobox floatright" style="width: 500px;"> | ||
+ | [[File:Imag1.png|right|thumb|300px|]] | ||
+ | </div> | ||
+ | |||
+ | |||
+ | {| 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]] |
Revision as of 18:23, 26 October 2015
Contents
OBJECTIVE: Liquid/Gaz Interface Simulations
Describe the interface motion precisely | Level Set Method |
---|
Handle jump conditions at the interface | without artificial smoothing | Ghost Fluid Method |
---|
Respect mass conservation | VOF Method |
---|
Solve incompressible Navier Stokes equations | Projection Method |
---|
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
Diesel type Jet Atomization : Diameter: 100µm, Liquid velocity 100 ms-1, Turbulence 5%, gas velocity: 0 ms-1, =696 kgm-3, =25 kgm-3
=1.210-3kgm-1s-1, =10-5kgm-1s-1, =0.06Nm-1 Numerical simulation by coupling Level Set / VOF / Ghost fluid methods
Mesh : 256x256x2048 (130 millions points) MPI parallelization 128 procs.
Triple disk injector
- 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
- Triple.avi
Liquid/Gas mixing layer
- 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
Air assisted atomization
- Air/water Jet « Marmottant » Modified Rudman's method two grids
- Ugas=35m/s and Uliq=0.2m/s
- =1,2 kg/m3 =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.
- Marmotte.avi
2D gas/liquid shear layer
- M=^2 / ^2 = 16
- Ugas=30 m/s Uliq=0.26 m/s
- =1,2 kg/m3 =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
- 2Dliquid_shear_layer.avi
Air assisted atomization
- Air/water Jet « Marmottant » RUDMAN-TYPE Technique with one grid
- Ugas=22.6m/s and Uliq=0.27m/s
- =1,2 kg/m3 =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.
- Marmotte.avi
Air assisted atomization
- Same configuration
- Initial boundary layers included
- Grid:512x512x1024
- Oct_15_2015.avi
- Helicoïdal behaviourin simulations (experimentally observed)
- helicoid.avi