HPC enabling of OpenFOAM R for CFD applications HPC simulation of volcanic ash plumes and application of OpenFOAM to CFD volcanological problems 06-08 April 2016, Casalecchio di Reno, BOLOGNA. Matteo Cerminara – [email protected]Tomaso Esposti Ongaro – [email protected]Mattia de’ Michieli Vitturi – [email protected]Istituto Nazionale di Geofisica e Vulcanologia Istituto Nazionale di Oceanografia e Geofisica Sperimentale
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Istituto Nazionale di Geofisica e VulcanologiaIstituto Nazionale di Oceanografia e Geofisica Sperimentale
Effusive and explosive eruptions:
The example of Etna volcano
effusive (lava flows)
(2006)
(lava fountains)
(2014)
explosive (ash plumes)
(2015)
Explosivity is mostly controlled by multiphase processes in the magma(melt+gas+crystals):
• Gas phase transitions (gas exsolution, bubble nucleation andexpansion).
• Non-linear magma rheology (brittle transition = fragmentation).
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 2 / 44
Challenges in
magma ascent modeling with OpenFOAM
Effusive regime
• Introduce multicomponent physics and phase transitions
• Bubble nucleation
• Manage phase separation and degassing
Explosive regime
• Wave propagation
• Fragmentation conditions
• Manage different domains (below/above fragmentation)
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 3 / 44
TwoPhaseChangeEulerFoam
A new multiphase multicomponent model with phase change has beendeveloped. Each phase is the mixture of several components andexsolution/evaporation laws can be defined for each component.
Test 1 (decompression experiment):
• validation through comparisonswith decompression experimentsperformed at LMU, Munich.
• silicon oil chosen as analogue formagmatic melt, and saturatedwith Argon at 10MPA
• slow decompression toatmospheric conditions
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 4 / 44
TwoPhaseChangeEulerFoam
A new multiphase multicomponent model with phase change has beendeveloped. Each phase is the mixture of several components andexsolution/evaporation laws can be defined for each component.
Test 1 (decompression experiment):
• validation through comparisonswith decompression experimentsperformed at LMU, Munich.
• silicon oil chosen as analogue formagmatic melt, and saturatedwith Argon at 10MPA
• slow decompression toatmospheric conditions
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 5 / 44
TwoPhaseChangeEulerFoam
A new multiphase multicomponent model with phase change has beendeveloped. Each phase is the mixture of several components andexsolution/evaporation laws can be defined for each component.
• Number of cells is necessarily large! Effective parallelism.
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 9 / 44
Model-related uncertainty
Weak Plume Strong Plume
Does uncertainty blur accuracy?
• A ”blind” intercomparison test (4 models)
• Model-related uncertainty is still much larger than• Error associated with grid size• Error associated with SGS model
• Measurement uncertainty is within the model error
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 10 / 44
Table of contents
1 The ASHEE model
2 Verification and validation study
3 Volcanologic application
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 11 / 44
Eulerian fields
A new model ASHEE (Ash Equilibrium Eulerian) has been developedtaking advantage of the OpenFOAM infrastructure
• it models a mixture of I gas species and J solid particle species, eachwith diameter dj
• the ith gas species is characterized by the following fields in eachpoint (x, t)
• bulk density ρi• velocity ug
• temperature Tg
• the jth solid species is characterized by the following fields in eachpoint (x, t)
• bulk density ρj• velocity uj
• temperature Tj
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 12 / 44
Mass averaged Eulerian fields
The conservation equations for mass, momentum and energy are writtenfocusing on the mass averaged mixture properties:
Mass averaged field ψm
Firstly we define the mixture density as ρm =∑I
i=1 ρi +∑J
j=1 ρj , so thatthe mass fractions of each phase are defined:
• yi = ρi/ρm
• yj = ρj/ρm.
Thus, given a generic field ψ(x, t) for the ith gas species (ψi ) or for thejth solid species (ψj), we define
ψm =I∑
i=1
yiψi +J∑
j=1
yjψj
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 13 / 44
Physical configuration
The Eulerian model in “mixture” formulation
∂tρm +∇ · (ρmum) =∑j∈J
Sj ;
∂t(ρmyi ) +∇ · (ρmugyi ) = 0, i ∈ I ;
∂t(ρmyj) +∇ · (ρmujyj) = Sj , j ∈ J ;
∂t(ρmum) +∇ · (ρmum ⊗ um + ρmTr) =
= −∇p +∇ · T + ρmg +∑j∈J
Sjuj ;
∂t(ρmhm) +∇ · [ρmhm(um + vh)] +
+ ∂t(ρmKm) +−∇ · [ρmKm(um + vK )] =
= ∂tp +∇ · (T · ug − q) + ρm(g · um) +∑j∈J
Sj(hj + Kj) .
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 14 / 44
Physical configuration
• The “mixture” formulation is useful in two-way coupled multiphasesystems, where the mass of the dispersed phase has a non-negligibleeffect on the dynamics.
• The problem is formulated in the dispersed regime εs . 10−3,collisions between particles are disregarded
• It focuses the mathematical problem on the mass averaged fields:ρm , um , hm (improving stability).
• All the effects due to the kinematic decoupling are confined in theterms Tr(yj , vj) , vh(yj , vj , hj , hm) , vK (yj , vj ,Kj ,Km) , keeping intoaccount the effect of the relative velocity vj = uj − ug
• No explicit dependence on the drag functional expression isnecessary at this level
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 15 / 44
Equilibrium Eulerian model
• if particles are perfectly coupled both thermally (Tj = Tg = T ) andkinematically (uj = ug = u) we recover the dusty gas model, whereall the decoupling terms are zero
• in ASHEE we adopt the equilibrium-Eulerian model, where thesystem is thermally perfectly coupled while the kinematic decouplingis approximated via an asymptotic expansion relative to the Stokestime τj :
uj = ug + wj − τj(ag + wj · ∇ug) + O(τj)
where wj = τjg is the settling velocity and ag is the gas phaseacceleration
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 16 / 44
Equilibrium Eulerian model
• The contribution of the particle inertia can be accurately taken intoaccount by using standard Navier-Stokes numerical algorithm,without needing to implicitly solve the drag term
• Particle decoupling and preferential concentration well modeled upto St . 0.2, keeping the advantages of the dusty gas model
• Total number of equation highly reduced for a polydispersed mixture(4J PDEs less):
• Eulerian: I + 3 + 5J• Equilibrium-Eulerian: I + 3 + J.
• Allows to solve efficiently the multiphase dynamics at geophysicalscale for particles of size up to ' 1 mm.
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 17 / 44
Paper on the ASHEE model
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 18 / 44
Numerical configuration
We implemented the following subgrid-scale LES models for the subgridterms of the compressible equilibrium-Eulerian model:
• Compressible Smagorinsky (static and dynamic)
• Turbulent Kinetic Energy model (static and dynamic)
• WALE model (static and dynamic)
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 19 / 44
Numerical configuration
We modified the compressible monophase PISO-PIMPLE algorithm:
• predictor for the mixture density ρm• PIMPLE loop
• solve for the mass fractions yi,j• predictor for the mixture velocity um
• solve for the enthalpy hm, fixing the compressibility• PISO loop
• solve the decoupling• solve for pressure p• correct velocity and fluxes
• solve LES models
• correct mixture density
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 20 / 44
Table of contents
1 The ASHEE model
2 Verification and validation study
3 Volcanologic application
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 21 / 44
DHITDecaying Homogeneous and Isotropic Turbulence (2563 cells)The solver is able to simulate accurately the turbulence.
k
E(k)
100 101 10210-10
10-9
10-8
10-7
10-6
10-5
10-4
10-3
10-2
DNS Bernardini & Pirozzoli
OpenFoam
Figure: Test case validation: comparison with an eight order DNS after one large-eddy turnovertime (10000 time steps).
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 22 / 44
Figure: Scalability test on PLX and FERMI environments (with little outputwork). Collaboration between us and Paride Dagna.
In order to fix ideas, the solver reach a velocity in the range 1÷10 Mcells/s on 1024 cores
(multiphase÷monophase) .
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 23 / 44
SGS LES models
1e-06
1e-05
0.0001
0.001
1 10
ener
gy s
pect
rum
E(k
)
wavenumber k
[noM] 2563
[noM] 323
[sma] 323
[oneEqEddy] 323
[wale] 323
Figure: DHIT using static SGS LES models in a box with 323 cells
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 24 / 44
SGS LES models
1e-06
1e-05
0.0001
0.001
1 10
ener
gy s
pect
rum
E(k
)
wavenumber k
[noM] 2563
[dynSma] 643
[moin] 323
[dynSma] 323
[dynOneEqEddy] 323
[dynWale] 323
[cubic] [dynOneEqEddy] 323
Figure: DHIT using dynamic SGS LES models in a box with 323 cells
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 25 / 44
Kinematic decoupling
Slice of the turbulent box at t/τe ' 2.2. The two panels representrespectively a logarithmic color map of yj=2 (Stmax = 0.5) and of |ag|
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 26 / 44
Kinematic decoupling
1e-05
0.0001
0.001
0.01
0.1
1
0.01 0.1 1
τ j τξ ⟨P⟩ j
(LE
S);
τ
j τη ⟨
P⟩ j
(D
NS
)
Stξ (LES); Stη (DNS)
[noM], 2563
[dynWale], 323
1.52*St2
Rani 2003
fit of Rani 2003
Figure: Evolution of the degree of preferential concentration with Stξ (LES) orStη (DNS). We obtain a good agreement between equilibriumEulerianLES/DNS and Lagrangian DNS simulations.
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 27 / 44
Advection schemes
Figure: Advection test case from Holzmann-cfd solved with ASHEE
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 28 / 44
Advection schemes
-0.2
0
0.2
0.4
0.6
0.8
1
1.2
-0.02 -0.015 -0.01 -0.005 0 0.005 0.01 0.015 0.02
trac
er
time
upwindlinearUpwind
linearlimitedLinear
limitedLimitedLinearvanLeer
limitedVanLeercubic
limitedCubicUMISTSFCD
QUICKfilteredLinear
Figure: Tracer mass fraction along the cavity section.
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 29 / 44
Advection schemes
10-6
10-5
10-4
10-3
1 10
E(k
)
k
[noM], 2563
[dynWale], 323
[dynWale], [MUSCL], 323
[dynWale], [UMIST], 323
[dynWale], [vanLeer], 323
[dynWale], [rk4blended], 323
[dynWale], [filtlin0200], 323
Figure: Performances of OpenFOAM schemes in DHIT LES
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 30 / 44
Other benchmarks
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
ρ
x/ct
analytic solutionrhoCentralFoam
ASHEE
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1
-2 -1.5 -1 -0.5 0 0.5 1 1.5 2
ρ
x/ct
Sod’s shock tube
1
10
10 100 1000
Nu
t [s]
Ra = 1e6, Nu = 8.800
Ra = 1e5, Nu = 4.519
Ra = 1e4, Nu = 2.243
Ra = 1e3, Nu = 1.118
Natural convection
Experimental plumeMixing
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 31 / 44
Table of contents
1 The ASHEE model
2 Verification and validation study
3 Volcanologic application
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 32 / 44
Plinian eruption
We have used ASHEE to simulate volcanic eruptions (scale ≈ 100 km)from the vent (mass eruption rate up to Mton/s) to the atmosphere
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 33 / 44
Volcanic plumes
We are participating to an international benchmark initiative involvingtwo numerical simulations:
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 34 / 44
Effects of resolution and SGS modelTime and horizontal average of the velocity field, with 8, 16 and 32 cellsin a vent diameter; with and without decoupling model (weakPlume)
Acknowledgments: we acknowledge CINECA for the availability ofhigh-performance computing resources and technical support on portingOpenFOAM on HPC architectures in the framework of ISCRA projects:IsB06 VolcFOAM, IsC26 VolcAshP and IsC07 GEOFOAM.
Matteo Cerminara et al. / HPC enabling of OpenFOAM for CFD applications / Modeling 44 / 44