Modeling and Simulation of Mutiphase Flows (CFD Tips) : Pouriya Niknam
Modeling and Simulation of Mutiphase Flows (CFD Tips) : Pouriya Niknam
Modeling and Simulation of Mutiphase Flows (CFD Tips) : Pouriya Niknam
Mutiphase flows
(CFD tips)
Pouriya Niknam
Supervisor: Dr. Daniele Fiaschi
2
Outline
Flow Specific
bubbly
? Process Specific
droplet Separation
particle-laden Filtration
slug Suspension
annular Evaporation
stratified/free surface Reaction
rapid granular flow
multiphase or multicomponent?
– Distinguish multiphase and/or multicomponent
• Liquid/Gas or Gas/Liquid
• Gas/Solid
• Liquid/Liquid
– Technically, two immiscible liquids are “multi-fluid”, but are
often referred to as a “multiphase” flow due to their similarity
in behavior
Single component Multi-component
Water Air
Single phase
Pure nitrogen H2O+oil emulsions
Coal particles in air
Multi-phase Steam bubble in H2O
Sand particle in H2O
– “separated” or “interfacial”
• both fluids are more or less contiguous throughout the
domain
– “dispersed”
• One of the fluids is dispersed as non-contiguous
isolated regions within the other (continuous) phase.
• The former is the “dispersed” phase, while the latter
is the “carrier” phase.
• One can now describe/classify the geometry of the
dispersion:
• Size & geometry
• Volume fraction
Gas-Liquid Flow
Aeration:
-produced by wave action
- used as reactor in chemical processing
- enhanced gas-liquid mass transfer
Gas-Liquid Flow
Energy production – liquid fuel combustion
Biomedical – inhalant drug delivery
Gas-Solid Flow
Increased complexity
Lagrangian
Track individual point particles.
Particles do not interact.
Algebraic slip model
Mixture model
Dispersed phase in a continuous phase.
Solve one momentum equation for the mixture.
Neither particle-wall interaction nor particle-particle are taken into account
Two-fluids theory (multi-fluids)
Eulerian-Eulerian models: two co-existing fluids
Solve as many momentum equations as there are phases.
Particle-wall interaction taken into account, particle-particle usually not.
Eulerian-granular model (EGM)
Both particle-wall and particle-particle interaction are taken into account
dispersed phase model (DPM)
Eulerian/Lagrangian
Solve the trajectories of individual objects and their collisions, inside a continuous phase.
Particle-wall interaction always taken into account, particle-particle usually not
Fully resolved and coupled.
Modeling approach, DPM
Trajectories of particles/droplets are
computed in a Lagrangian frame.
Exchange (couple) heat, mass, and momentum
with Eulerian frame gas phase.
Discrete phase volume fraction should
preferably be less than 10%.
Mass loading can be large (+100%).
No particle-particle interaction or break up.
Turbulent dispersion modeled by:
Stochastic tracking.
Particle cloud model.
Model particle separation, spray drying,
liquid fuel or coal combustion, etc.
Multiphase flow regimes
Solve U-Momentum
• There are two kinds of solvers available
in FLUENT – Pressure based and Density Solve V-Momentum Solve Mass,
based. Solve Mass Momentum,
Solve W-Momentum & Momentum Energy,
Species
• The pressure-based solvers take Solve Mass
momentum and pressure (or pressure Continuity;
Update Velocity
correction) as the primary variables.
– Pressure-velocity coupling algorithms are
derived by reformatting the continuity Solve Energy
equation
• Two algorithms are available with the Solve Species
pressure-based solvers:
– Segregated solver – Solves for pressure
Solve Turbulence Equation(s)
correction and momentum sequentially.
– Coupled Solver (PBCS) – Solves pressure
and momentum simultaneously. Solve Other Transport Equations as required
Choosing a Solver
• The pressure-based solver is applicable for a wide range of flow regimes from low
speed incompressible flow to high-speed compressible flow.
– Requires less memory (storage).
– Allows flexibility in the solution procedure.
• The pressure-based coupled solver (PBCS) is applicable for most single phase
flows, and yields superior performance to the standard pressure-based solver.
– Now available for multiphase (Eulerian)
– Requires 1.5–2 times more memory than the segregated solver.
• The density-based coupled solver (DBCS) is applicable when there is a strong
coupling, or interdependence, between density, energy, momentum, and/or
species.
– Examples: High speed compressible flow with combustion, hypersonic flows, shock
interactions.
Convergence Difficulties with multiphase
• Numerical instabilities can arise with an ill-posed problem, poor-quality mesh and/or
inappropriate solver settings.
– Exhibited as increasing (diverging) or “stuck” residuals.
– Unconverged results are very misleading!
• Troubleshooting
1. Ensure that the problem is well-posed.
2. Compute an initial solution using a
first-order discretization scheme.
3. For the pressure-based solver, decrease Continuity equation convergence
underrelaxation factors for equations trouble affects convergence of
having convergence problems. all equations.
4. For the density-based solver, reduce
the Courant number.
5. Disabling Volume fraction &phase equations
6. Remesh or refine cells which have large
aspect ratio or large skewness.
(Remember that you cannot improve
cell skewness by using mesh adaption!)
Modifying Under-Relaxation Factors
• Under-relaxation factor, α, is included to stabilize the
iterative process for the pressure-based solver
• Use default under-relaxation factors to start a calculation.
• Decreasing under-relaxation
for momentum often aids
convergence.
– Default settings are suitable for a
wide range of problems, you can
reduce the values when necessary.
– Appropriate settings are best learned
from experience!
Modifying the Courant Number
• A transient term is included in the density-based solver even for
steady state problems.
– The Courant number defines the
time step size.
• For density-based explicit solver:
– Stability constraints impose a
maximum limit on the Courant
number.
• Cannot be greater than 2
(default value is 1).
• Reduce the Courant number when
having difficulty converging.
• For density-based implicit solver:
– The Courant number is not limited
by stability constraints.
• Default value is 5.
Starting from a Previous Solution
• A previously calculated solution can be
used as an initial condition when
changes are made to the case setup.
– Use solution interpolation to
initialize a run (especially
useful for starting fine-mesh cases
when coarse-mesh solutions are
available).
– Once the solution is initialized,
additional iterations always use
the current data set as the starting
Actual Problem Initial Condition
point.
Heat Transfer Isothermal
– Some suggestions on how to
Natural convection Low Rayleigh number
provide initial conditions for some
actual problems: Turbulence Inviscid (Euler) solution
Steady or Unsteady
Prel,max=100,001 Pa Prel,max=1 Pa
Prel,min=99,999 Pa Prel,min=-1 Pa
Pref
F F
c 0
t x
n 1 n 1 n 1
Fi Fi n
Fi F
c i 1
0
t x
Propagation Direction
i -1 i
Let’s revisit 1-D Vorticity Transport Equation:
F F
c 0
t x
Fi Fi Fi 1 Fi x Fi 2 Fi Fi
x x 2x 2 x 2
Second order
First order
Fourth order
Third order(MUSCL)
Sixth order
Fifth order (WENO)
Eighth order
MUSCL
A q R q L
1
Fi 1/ 2 Symmetric Part
2
This part is used to control This part is used to
dispersion and truncation errors control dissipation errors
2nd order:
1
Fi 1 Fi
2
4th order:
1
Fi 2 7Fi1 7Fi Fi1
12
6th order:
1
Fi3 8Fi2 37Fi1 37Fi 8Fi1 Fi2
30
MUSCL
A q R q L
1
Fi 1/ 2 Symmetric Part 2
LEFT RIGHT
STENCIL
qL qR STENCIL
3 6
qR qi 1 q i 2 qi 1 qi 1 qi
1 1
6 3
qi qi 1 qi q i qi 1
1 1
qL
qi qi 1 qi q i qi 1
1
3 1
6
1st order MUSCL: qL
31 6 1
qR qi 1 q i 2 qi 1 qi 1 qi
q q q q q
1
6 1
3
qR i 1 i2 i 1 i 1 i
33
Question?