Customer Training Material

L t
Lecture
9
Advanced Physics
y

Introduction to
ANSYS FLUENT
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-1

Release 13.0
December 2010

Advanced Modelling Options

Outline

Customer Training Material

Multiphase Flow Modelling

Discrete phase model

Eulerian model
Mixture model
Volume-of-Fluid (VOF) model

Reacting Flow Modelling

Eddy dissipation model

Non-premixed, premixed and partially premixed combustion models
Detailed chemistry models
Pollutant formation
Surface reactions

Modelling Moving Parts

Single
Si
l and
d multiple
lti l reference
f
fframes
Mixing planes
Sliding meshes
Dynamic
y
meshes
Six-degree-of-freedom solver

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-2

Release 13.0
December 2010

Advanced Modelling Options

Customer Training Material

Multiphase Flows

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-3

Release 13.0
December 2010

Advanced Modelling Options

Multiphase Flows - Introduction

Customer Training Material

In many flows, there is more than one fluid present in the domain
Different substances (eg oil & water, or methane & air)
Different phases of same substance (water & steam)

The key
y issue is how these two fluids are mixed
If they are mixed at a molecular level, the problem is a multi-species flow.
A common example is where two gases are present (methane and air)
A diffusivityy (material
(
property)
p p y ) is set for the mixture,, and one extra transport
p
equation is solved for the mass fraction of primary component.
If the mixing is more macroscopic, then it is a multiphase flow.
In such cases there is an identifiable boundary between the two phases
The user must therefore indicate to the solver how this boundary performs
maybe a free surface (VOF model), or a typical droplet size (mixture model)

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-4

Release 13.0
December 2010

Advanced Modelling Options

Multiphase Flows - Introduction

Customer Training Material

The fluid system is defined by a primary and multiple secondary phases.

One of the phases is considered continuous (primary)
The others (secondary) are considered to be dispersed within the continuous
phase.
(Note that for free-surface flows, using the Volume of Fluid model, a distinct
interface is defined between the phases and both could be considered
continuous)
Secondary phase(s)

Primary Phase

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-5

Release 13.0
December 2010

Advanced Modelling Options

Multiphase Flow Regimes

Gas/Liquid
Liquid/Liquid

Gas / Solid

Liquid / Solid

Customer Training Material

Bubbly flow Discrete gaseous bubbles

in a continuous fluid, e.g. absorbers,
evaporators, sparging devices.
Droplet flow Discrete fluid droplets in a
continuous gas, e.g. atomizers,
combustors
Slug flow Large bubbles in a continuous
liquid
Stratified / free-surface flow
Immiscible fluids separated by a clearly
defined interface, e.g. free-surface flow
Particle-laden flow Discrete solid
particles in a continuous fluid, e.g. cyclone
separators, air classifiers, dust collectors,
dust-laden environmental flows
Fluidized beds Fluidized bed reactors
Slurry flow Particle flow in liquids
liquids, solids
suspension, sedimentation, and hydrotransport

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-6

Slug Flow

Bubbly, Droplet, or
Particle-Laden Flow

Stratified / FreePneumatic Transport,

Surface Flow Hydrotransport, or Slurry Flow

Sedimentation

Fluidized Bed
Release 13.0
December 2010

Advanced Modelling Options

Multiphase Models Available in FLUENT

Customer Training Material

FLUENT contains four distinct multiphase modeling approaches:

Discrete Phase Model (DPM)

Volume of Fluid Model (VOF)
Eulerian Model
Mixture Model

It is important to select the most appropriate modelling approach when

attempting to model a multiphase flow.
Depends
p
on whether the flow is stratified or disperse
p
length
g scale of the
interface between the phases dictates this.
Also the Stokes number (the ratio of the particle relaxation time to the
characteristic time scale of the flow) should be considered.
where

and

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-7

Release 13.0
December 2010

Advanced Modelling Options

DPM Example Spray Dryer

Customer Training Material

Spray drying involves the transformation of a liquid spray into dry powder
in a heated chamber. The flow, heat, and mass transfer are simulated
using the DPM model in FLUENT.

Initial particle
Diameter: 2 mm

1.1 mm

0.2 mm

Contours of
Evaporated
Water
Stochastic Particle Trajectories for Different Initial Diameters

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-8

Release 13.0
December 2010

Advanced Modelling Options

Eulerian Model Example 3D Bubble Column

Customer Training Material

z = 20 cm

z = 15 cm

z = 10 cm

z = 5 cm

Isosurface of Gas
Volume Fraction = 0.175
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

Liquid Velocity Vectors

L9-9

Release 13.0
December 2010

Advanced Modelling Options

The Granular Option in the Eulerian Model

Customer Training Material

Granular flows occur when a high concentration of

solid particles is present. This leads to high
frequency of interparticle collisions.
Particles are assumed to behave similarly to a
dense cloud of colliding molecules. Molecular
cloud theory is applied to the particle phase
phase.

Gravity

Application of this theory leads to the appearance

of additional stresses in the momentum equations
q
for continuous and particle phases
These stresses (granular viscosity, pressure
etc.) are determined by intensity of particle velocity
fluctuations
Kinetic energy associated with particle velocity
fluctuations is represented by a pseudo-thermal or
granular temperature
Inelasticity of the granular phase is taken into
account
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-10

Gas / Sand
Gas

Contours of Solids Volume

Fraction for High Velocity
Gas/Sand Production
Release 13.0
December 2010

Advanced Modelling Options

Mixture Model Example Gas Sparging

Customer Training Material

The sparging of nitrogen

gas into a stirred tank is
simulated by the mixture
multiphase model
model. The
rotating impeller is
simulated using the
multiple reference frame
(MRF) approach.
FLUENT simulation
provided
id d a good
d prediction
di ti
on the gas-holdup of the
agitation system.

Animation of Gas
Volume Fraction
Contours
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-11

Water Velocity Vectors on

a Central Plane at
t = 15 sec.
Release 13.0
December 2010

Advanced Modelling Options

VOF Example Automobile Fuel Tank Sloshing

Sloshing (free surface
movement) of liquid in an
automotive fuel tank
under various
accelerating conditions is
simulated by the VOF
model in FLUENT.
Simulation shows that the
tank with internal baffles
will
ill kkeep th
the ffuell iintake
t k
orifice fully submerged at
all times compared to the
tank with no baffles.

Customer Training Material

Fuel Tank Without Baffles

t = 1.05 sec

t = 2.05 sec

Fuel Tank With Baffles

ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-12

Release 13.0
December 2010

Advanced Modelling Options

Learning More

Customer Training Material

Many more workshops are available.

Look at www.fluentusers.com (same password as ANSYS portal)
Follow link to ANSYS FLUENT > Advanced / Intermediate Tutorials

Heat and Mass Transfer with the Mixture Model

Hydrodynamics of Bubble Column Reactors
Horizontal Film Boiling
Dam-Break Simulation Using FLUENT's Volume of Fluid Model
Using FLUENT's Erosion Model to Investigate Erosion in a 90 degree Elbow Bend
Using the Eulerian Multiphase Model with Species Transport
Modeling Flow and Heat Transfer in Packed Bed Reactor
Modeling Rapid Condensation of Steam in a 2D Laval Nozzle
S l i a 2D B
Solving
Box F
Falling
lli iinto W
Water
Modeling Uniform Fluidization in 2D Fluidized Bed
Modeling Bubble Breakup and Coalescence in a Bubble Column Reactor
Spin Coating of a Rotating Circular Glass Substrate
F l Tank
Fuel
T k Sloshing
Sl hi
Continuous Steel Casting of a Round Billet
Modeling the Effect of Sedimentation Concentration using a UDF
Modeling Nucleate Boiling using FLUENT

We also strongly recommend attending our Advanced Multiphase Training

ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-13

Release 13.0
December 2010

Advanced Modelling Options

Customer Training Material

Reacting Flows

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-14

Release 13.0
December 2010

Advanced Modelling Options

Reacting Flows Introduction

Customer Training Material

So far we have assumed that whatever materials are entering the domain
are the same as those leaving the domain.
However in some cases the materials entering will react with each other to
form new p
products (CO2 , H2O , NO
Ox etc))
By defining the reaction chemistry and kinetics, FLUENT can compute the
chemical reaction
reaction.
Within the flow domain, we have already seen how the solver can compute the
species concentration and temperature.
This can then be combined with knowledge of the reaction to form new species
in the model, with a corresponding transfer of energy.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-15

Release 13.0
December 2010

Advanced Modelling Options

Applications of Reacting Flow Systems

Customer Training Material

FLUENT contains models which are applicable to a wide range of

homogeneous and heterogeneous reacting flows

Furnaces
Boilers
Process heaters
Gas turbines
Rocket engines
IC engine
CVD, catalytic reactions

Temperature in a Gas Furnace

CO2 Mass Fraction

Predictions of
Flow field and mixing
characteristics
Temperature field
Species
p
concentrations
Particulates and pollutants
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

Stream Function

L9-16

Release 13.0
December 2010

Advanced Modelling Options

Background

Customer Training Material

Modeling Chemical Kinetics in Combustion

Fast Chemistry
Global chemical reaction mechanisms (Finite Rate / Eddy Dissipation)
Equilibrium/flamelet model (Mixture fraction)
Finite rate chemistry
Fuel

Reactor

Flow configuration
Non-premixed reaction systems
Can be simplified to a mixing problem

Oxidizer

Fuel
+
Oxidizer

Premixed reaction systems

Outlet

Reactor

Outlet

Cold reactants propagate into hot products.

Secondary
Fuel or Oxidizer

Partially premixed systems

Reacting system with both non-premixed
and premixed inlet streams.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-17

Fuel
+
Oxidizer

Reactor

Outlet

Release 13.0
December 2010

Advanced Modelling Options

Overview of Reacting Flow Models in FLUENT

Customer Training Material

FLOW CONFIGURATION
Premixed

Non-Premixed

Partially Premixed

CHEMISTRY
Y

Eddy Dissipation Model

(Species Transport)

Fast
Chemistry

Premixed
Combustion
Model
Reaction Progress
Variable*

Non-Premixed
Equilibrium
Model

Partially Premixed
Model

Mixture Fraction

Reaction Progress
Variable
+
Mixture Fraction

Laminar Flamelet Model

Finite-Rate
Ch i t
Chemistry

Laminar Finite-Rate Model

Eddy-Dissipation Concept (EDC) Model
Composition PDF Transport Model

*Rate classification not truly applicable since species mass fraction is not determined.
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-18

Release 13.0
December 2010

Advanced Modelling Options

Pollutant Formation Models

Customer Training Material

NOx formation models (predict qualitative trends of NOx formation).

FLUENT contains three mechanisms for calculating NOx production.
Thermal NOx
Prompt NOx
Fuel NOx

NOx reburning model

Selective Non-Catalytic
y Reduction ((SNCR)) model
Ammonia and urea injection

Soot formation models

Moos-Brookes model
One step and two steps model
Soot affects the radiation absorption (Enable the Soot-Radiation option in the
Soot panel)

SOx formation models

Additional equations for SO2
SO2, H2S
H2S, and
and, optionally
optionally, SO3 are solved
solved.
In general, SOx prediction is performed as a post-process.
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-19

Release 13.0
December 2010

Advanced Modelling Options

Discrete Phase Model (DPM)

Customer Training Material

The reactions need not be restricted to the fluids in the model, the Discrete Phase
model can also be used to simulate reacting particles for example the
combustion of pulverised coal particulates.

Numerous submodels are available.

Heating/cooling of the discrete phase

Vaporization and boiling of liquid droplets
Volatile evolution and char combustion for combusting particles
Droplet breakup and coalescence using spray models
Erosion/Accretion

Numerous applications
Particle separation
p
and classification,, spray
p y drying,
y g, aerosol dispersion,
p
, bubble sparging
p g g
of liquids, liquid fuel and coal combustion.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-20

Release 13.0
December 2010

Advanced Modelling Options

Surface Reactions

Customer Training Material

Some reactions occur just on the wall surface (for example, in a vehicles
catalytic converter).
Chemical species deposited onto surfaces are treated as distinct from the
same chemical species in the gas.
Site balance equation is solved for every surface-adsorbed (or site)
species.
species
Detailed surface reaction mechanisms can be considered (any number of
reaction steps and any number of gas-phases or/and site species).
Surface chemistry mechanism in Surface CHEMKIN format can be imported into
FLUENT.
Surface reaction can occur at a wall or in porous media.
Different surface reaction mechanisms can be specified on different surfaces.

Application examples
Catalytic reactions
CVD

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-21

Release 13.0
December 2010

Advanced Modelling Options

Learning More

Customer Training Material

Many more workshops are available.

Look at www.fluentusers.com (same password as ANSYS portal)
Follow link to ANSYS FLUENT > Advanced / Intermediate Tutorials

Coal Combustion with Eddy Break Up (EBU) Model

Modeling Liquid Reactions in a CIJR using the Unsteady Laminar Flamelet Model
Simulation of a Piloted Jet Flame using Unsteady Laminar Flamelet Model
2D Simulation of a 300 KW BERL Combustor Using the Magnussen Model
Premixed Flow in a Conical Chamber using the Finite-Rate Chemistry Model
PDF Transport Simulation of a Piloted Jet Diffusion Flame
Liquid Fuel Combustion
3D Simulation of 300 kW BERL Combustor Using the Laminar Flamelet Model
NO Modeling
NOx
M d li with
i h the
h SNCR M
Model
d lU
Using
i U
Urea IInjection
j i
Modeling Liquid Reactions in a CIJR using the Unsteady Laminar Flamelet Model
Simulation of a Piloted Jet Flame using Unsteady Laminar Flamelet Model
EDC Simulation of a Piloted Jet Diffusion Flame
P
Premixed
i dC
Combustion
b ti iin a C
Conical
i l Ch
Chamber
b using
i th
the Zimont
Zi
t Model
M d l
Modeling Surface Reaction in a Single Circular Channel
Multiple Char Reactions
Partially Premixed Combustion in a Coaxial Combustor
M d li E
Modeling
Evaporation
ti off Liquid
Li id Droplets
D l t in
i a Circular
Ci l Ch
Channell

We also strongly recommend attending our Advanced Combustion Course

ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-22

Release 13.0
December 2010

Advanced Modelling Options

Customer Training Material

Modelling Moving Parts

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-23

Release 13.0
December 2010

Advanced Modelling Options

Introduction

Customer Training Material

Many flow problems involve domains which

exhibit forms of motion.
Two types off motion are possible translational
and rotational.
Th
There are two
t
modeling
d li approaches
h ffor moving
i
domains:
Moving Reference Frames
Frame of reference is attached to the moving domain
domain.
Governing equations are modified to account for
moving frame.

Moving / Deforming Domains

Domain position and shape are tracked with respect
to a stationary reference frame.
Solutions are inherently transient.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-24

Release 13.0
December 2010

Advanced Modelling Options

Single Reference Frame (SRF) Modeling

Customer Training Material

SRF attaches a reference frame to a single

moving domain.
All fluid motion is defined with respect to the
g frame.
moving
Rotating frames introduce additional
accelerations to the equations of fluid
mechanics, which are added by Fluent
when you activate a moving reference
frame.

Why use a moving reference frame?

Centrifugal
Compressor
(single blade passage)

Flow field which is transient when viewed in

a stationary frame can become steady
when viewed in a moving frame.
Advantages

Steady state solution*

so ut o
Simpler BCs
Faster turn-around time
Easier to post-process and analyze

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

* NOTE: You may still have unsteadiness

in the rotating frame due to turbulence,
circumferentially non-uniform variations in
flow, separation, etc. example: vortex
shedding from fan blade trailing edge
L9-25

Release 13.0
December 2010

Advanced Modelling Options

Multiple Reference Frame (MRF) Modeling

Many moving zone problems involve
stationary components which cannot be
described by surfaces of revolution
(SRF not valid).
valid)
Systems like these can be solved by
dividing the domain into multiple fluid
zones some zones will be rotating,
g
others stationary.
The multiple zones communicate
across one or more interfaces (these
may or may nott be
b non-conformal)
f
l)
The way in which the interface is
treated leads to one of following
approaches for multiple zone models:
Multiple Reference Frame Model (MRF)
Mixing Plane Model (MPM)
Sliding mesh model (SMM)

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-26

Customer Training Material

interface

Multiple Component
(blower wheel + casing)
Steady State (Approximate)
Transient (Best Accuracy)

Release 13.0
December 2010

Advanced Modelling Options

The Mixing Plane Model (MPM)

Customer Training Material

The MPM is a technique which permits

steady-state solutions for multistage axial
and centrifugal turbomachines.
Domain is comprised of multiple zones
Each zone is self contained with a inlet,
outlet, wall, periodic BCs
Steady-state SRF solutions are obtained in
each domain, with the domains linked by
passing boundary conditions.
The BC links between the domains are
called mixing planes.
Mixing plane
(P
(Pressure
outlet
tl t linked
li k d with
ith
a mass flow inlet)

BCs are passed as circumferentially

averaged profiles of flow variables, which
are updated at each iteration.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

ADVANTAGE of MPM: Requires only a

single blade passage per blade row
regardless of the number of blades.
L9-27

Release 13.0
December 2010

Advanced Modelling Options

The Sliding Mesh Model (SMM)

Customer Training Material

Although MRF and MPM are good for

indicating the overall behaviour arising from
the motion, they will not capture the transient
j
((and its wake)) p
passes
detail as one object
another.

Shock
interaction

potential
t ti l
interaction

So, for example, the image shows where

there may be a transient issue with shock
interaction between rotor and stator.
If transient interaction can not be neglected,
we can employ the Sliding Mesh model
(SMM) to account for the relative motion
between the stationary and rotating
components.
p

Stator

Rotor

wake interaction

Sliding Mesh cases are always run in a

transient manner, and one region of grid cells
is moved (rotated or translated) each time
step.
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-28

Release 13.0
December 2010

Advanced Modelling Options

How the Sliding Mesh Model Works

Customer Training Material

Like the MRF model, the domain is divided into moving and stationary zones,
separated by non-conformal interfaces.
moving mesh zone

cells at time t

cells at time t + t

Governing equations have a new moving mesh form, and are solved in the
stationary reference frame for absolute quantities (see Appendix for more
details).
Moving reference frame formulation is NOT used here (i.e. no additional accelerations
acting as sources terms in the momentum equations).
Equations are a special case of the general moving/deforming mesh formulation.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-29

Release 13.0
December 2010

Advanced Modelling Options

The Dynamic Mesh (DMM) Model

Customer Training Material

FLUENT can take mesh motion one step further, and modify the mesh in
the solver at every time step so as to resolve larger motions (that would be
impossible to account for with sliding meshes or MRF
Examples:
Automotive p
piston moving
g
inside a cylinder
Positive displacement
pumps
A flap moving on an
airplane wing
A valve opening and
closing
An artery expanding
and contracting

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-30

Release 13.0
December 2010

Advanced Modelling Options

Dynamic Mesh (DM) Methods

Customer Training Material

Internal node positions are automatically calculated based on user

specified boundary/object motion, cell type, and meshing schemes
Basic Schemes
Smoothing (Spring analogy)
Local remeshing
g
Layering

Other Methods

2.5 D
User defined mesh motion
In cylinder motion (RPM
In-cylinder
(RPM, stroke length
length, crank angle
angle, ))
Prescribed motion via profiles or UDF
Coupled motion based on hydrodynamic forces from the flow solution, via
FLUENTss six-degree-of-freedom
FLUENT
six degree of freedom (6DOF) solver.
solver

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-31

Release 13.0
December 2010

Advanced Modelling Options

Dynamic Mesh Methods

Layering
Layers of cells are generated
and collapsed as they are
overrun by the moving
boundary. Layering is
appropriate for quad/hex/prism
meshes with linear or rotational
motion and can tolerate small or
large boundary deflections.
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

Customer Training Material

Local Remeshing

Spring Analogy

In local remeshing,
g as cells become Spring
p g analogy
gy is useful when there
skewed due to moving boundaries,
are small boundary deformations.
cells are collapsed and the skewed
The connectivity and cell count is
region is remeshed. Local remeshing unchanged during motion. Spring
is appropriate for tri/tet meshes with
analogy is appropriate for tri/tet
large range of boundary motion
motion.
meshes with small deformations
deformations.

L9-32

Release 13.0
December 2010

Advanced Modelling Options

Summary

Customer Training Material

Five different approaches may be used to model flows over moving parts.

Single (Rotating) Reference Frame Model

Multiple Reference Frame Model
Mixing Plane Model
Sliding Mesh Model
Dynamic Mesh Model

First three methods are primarily steady-state approaches while sliding

mesh and dynamic
y
mesh are inherently
y transient.
Enabling these models, involves in part, changing the stationary fluid
zones to either Moving Reference Frame or Moving Mesh
Mesh.
Most physical models are compatible with moving reference frames or
moving meshes (e.g.
(e g multiphase,
multiphase combustion
combustion, heat transfer,
transfer etc
etc.))
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-33

Release 13.0
December 2010

Advanced Modelling Options

Learning More

Customer Training Material

Many more workshops are available.

Look at www.fluentusers.com (same password as ANSYS portal)
Follow link to ANSYS FLUENT > Advanced / Intermediate Tutorials

2D Adiabatic Compression (Remeshing and Spring Smoothing)

2D Adiabatic Compression (Layering)
3D Adiabatic Compression (Layering, Remeshing, and Spring Smoothing)
Solving a 2D Box Falling into Water
Simulating a 2D Check Valve using FLUENT's DMM and Spring Smoothing
Submarine Docking Simulation Using MDM Model
Using
Us
gaU
UDF to
o Co
Control
o the
e Dynamic
y a c Mesh
es o
of a Flexible
e b e Osc
Oscillating
a g Membrane
e ba e
Solving a 2D Vibromixer Problem Using the Dynamic Mesh Model
Store Separation from a 3D Delta Wing
Projectile Moving Inside a Barrel

We also strongly recommend attending our Dynamic Mesh Training

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-34

Release 13.0
December 2010

Appendix 1

Multiphase Flow
Modeling
ANSYS, Inc. Proprietary
2009 ANSYS, Inc. All rights reserved.

9-35

April 28, 2009

Inventory #002600

Advanced Modelling Options

Volume and Particulate Loading

Customer Training Material

Volume loading dilute vs. dense

Refers to the volume fraction of secondary phase(s)

For dilute loading (less than around 10%), the

average inter-particle
i t
ti l di
distance
t
iis around
d ttwice
i
the particle diameter. Thus, interactions among
particles can be neglected.

Particulate loading ratio of dispersed and continuous phase inertia

inertia.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-36

Release 13.0
December 2010

Advanced Modelling Options

Turbulence Modeling in Multiphase Flows

Customer Training Material

Turbulence modeling with multiphase flows is challenging.

single phase turbulence models (such as k
k or RSM) are used
Presently, single-phase
to model turbulence in the primary phase only.
Turbulence equations may contain additional terms to account for
turbulence modification by secondary phase(s).
If phases are separated and the density ratio is of order 1 or if the particle
volume fraction is low (< 10%), then a single-phase model can be used to
represent the mixture.
In other cases, either single phase models are still used or particlepresence-modified models are used.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-37

Release 13.0
December 2010

Advanced Modelling Options

Phases as Mixtures of Species

Customer Training Material

In all multiphase models within FLUENT, any phase can be composed of

either a single material or a mixture of species.
Material definition of phase mixtures is the same as in single phase flows.
It is possible to model heterogeneous reactions (reactions where the
reactants and products belong to different phases).
This means that heterogeneous reactions will lead to interfacial mass transfer.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-38

Release 13.0
December 2010

Advanced Modelling Options

Discrete Phase Model (DPM) Overview

Customer Training Material

Trajectories of particles, droplets or bubbles are computed in a Lagrangian frame.

Particles can exchange heat, mass, and momentum with the continuous gas phase.
Each trajectory represents a group of particles, all with the same initial conditions.
DPM neglects collisions and other inter-particle interactions.
Turbulent dispersion of particles can be modeled using either stochastic tracking (the most common
method) or a particle cloud model.

Many submodels are available Heat transfer, vaporization/boiling, combustion,

breakup/coalescence, erosion/accretion.
Applicability of DPM

Flow regime:
Volume loading:
P ti l t L
Particulate
Loading:
di
Stokes Number:

Bubbly flow, droplet flow, particle-laden flow

Must be dilute (volume fraction < 12%), otherwise use Dense DPM Model
L
Low
tto moderate
d t
All ranges of Stokes number

Application examples

Cyclones
C
Spray dryers
Particle separation and classification
Aerosol dispersion
Liquid fuel
Coal combustion

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-39

Release 13.0
December 2010

Advanced Modelling Options

Discrete Phase Model (DPM) Setup

Define

Define

Models

Customer Training Material

Discrete Phase

Injections

Display
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

Particle Tracks
L9-40

Release 13.0
December 2010

Advanced Modelling Options

DPM Boundary Conditions

Customer Training Material

Escape

Trap

Reflect
Wall-jet

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-41

Release 13.0
December 2010

Advanced Modelling Options

The Eulerian Multiphase Model

Customer Training Material

The Eulerian multiphase model is a multi-fluid model. This means that all phases
are assumed to exist simultaneously.
Conservation equations for each phase contain single-phase terms (pressure gradient,
thermal conduction etc.))
Conservation equations also contain interfacial terms (drag, lift, mass transfer, etc.).

Interfacial terms are generally nonlinear and therefore, convergence can

sometimes be difficult.
Eulerian Model applicability
Flow regime

Bubbly flow, droplet flow, slurry flow,

fluidized bed, particle-laden flow
Volume loading
Dilute to dense
Particulate
P ti l t lloading
di
L
Low
tto hi
high
h
Stokes number
All ranges

Application examples

High particle loading flows

Slurry flows
Sedimentation
Fluidized beds
Risers
Packed bed reactors

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-42

Release 13.0
December 2010

Advanced Modelling Options

Eulerian Multiphase Model Equations

Continuity:
( q q )
t

Customer Training Material

Volume fraction for the qth phase

+ ( q q u q ) = m& pq
n

p =1

Momentum for qth phase:

( q q u q )
t
transient

+ ( q q u q u q ) = q p + q q g + q +

(R

pq

+ m& pq u q ) + q q (Fq + Flift ,q + Fvm,q )

p =1

convection

pressure

body

shear

Solids pressure term is included

for granular model
model.

interphase
interphase mass
forces exchange

external, lift, and

virtual mass forces

exchange

The inter-phase exchange forces are expressed as:

In g
general: Fpq = Fqp
R pq = K pq (u p u q )
Energy equation for the qth phase can be similarly formulated.

Exchange coefficient

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-43

Release 13.0
December 2010

Advanced Modelling Options

Eulerian Multiphase Model Equations

Customer Training Material

Multiphase species transport for species i belonging to mixture of qth

phase
Mass fraction of species i in qth phase

q q q
Yi + q q u qYi q = q J iq + q Riq + q Siq +
t

ttransient
a se t

convective

diffusion

(m&
n

pi q j

m& q j p i

p =1

homogeneous
reaction
heterogeneous
homogeneous
reaction
production

H
Homogeneous and
dh
heterogeneous
t
reactions
ti
are setup
t th
the same as iin
single phase
The same species may belong to different phases without any relation
between themselves

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-44

Release 13.0
December 2010

Advanced Modelling Options

Eulerian Model Setup

Define

Define

Customer Training Material

Phases

Models

Viscous

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-45

Release 13.0
December 2010

Advanced Modelling Options

Mixture Model Overview

Customer Training Material

The mixture model is a simplified Eulerian approach, based on the

assumption of small Stokes number.
Solves the mixture momentum equation (for mass-averaged mixture velocity)
Solves a volume fraction transport equation for each secondary phase.

Mixture model applicability

Flow regime:
Volume loading:
Particulate Loading:
Stokes Number:

Bubbly, droplet, and slurry flows

Dilute to moderately dense
Low to moderate
St << 1

Application examples

Hydrocyclones
Bubble column reactors
Solid suspensions
Gas sparging

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-46

Release 13.0
December 2010

Advanced Modelling Options

Mixture Model Equations

Customer Training Material

Solves one equation for continuity of the mixture

Solves for the transport of volume fraction of each secondary phase

Drift velocity
Solves one equation for the momentum of the mixture

The mixture properties are defined as:

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-47

Release 13.0
December 2010

Advanced Modelling Options

Mixture Model Setup

Define

Define

Models

Customer Training Material

Multiphase

Phases

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-48

Release 13.0
December 2010

Advanced Modelling Options

Mixture Model Setup

Customer Training Material

Boundary Conditions

Volume fraction defined for each

secondary phase.
phase location,
To define initial p
patch volume fractions after
solution initialization.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-49

Release 13.0
December 2010

Advanced Modelling Options

Cavitation Submodel

Customer Training Material

The Cavitation model models the

formation of bubbles when the
local liquid pressure is below the
vapor pressure.
The effect of non-condensable
gases is included.
Mass conservation equation for
the vapor phase includes vapor
generation and condensation
t
terms
which
hi h d
depend
d on th
the sign
i off th
the diff
difference between
b t
local
l
l pressure
and vapor saturation pressure (corrected for on-condensable gas
presence).
Generally
G
ll used
d with
ith th
the mixture
i t
model,
d l iincompatible
tibl with
ith VOF
VOF.
Tutorial is available for learning the in-depth setup procedure.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-50

Release 13.0
December 2010

Advanced Modelling Options

Eulerian-Granular Model Setup

Customer Training Material

Granular option must be enabled when defining

the secondary phases.
Granular properties require definition.
Phase interaction models appropriate for
granular flows must be selected.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-51

Release 13.0
December 2010

Advanced Modelling Options

The Volume of Fluid (VOF) Model Overview

Customer Training Material

The VOF model is designed to track the location and motion of a free surface
between two or more immiscible fluids.
VOF model can account for:
Turbulence, energy and species transport
Surface tension and wall adhesion effects.
Compressibility of phase(s)

VOF model applicability:

Flow regime
Volume loading
Particulate loading
Turbulence modeling
Stokes number

Slug flow, stratified/free-surface flow

Dilute to dense
Low to high
Weak to moderate coupling between phases
All ranges

Application examples

Large slug flows

Tank filling
Offshore separator sloshing
Coating

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-52

Release 13.0
December 2010

Advanced Modelling Options

VOF Model Setup

Define

Models

Define

Phases

Customer Training Material

Multiphase

Define

Operating Conditions

Operating Density should be set to that off

O
lightest phase with body forces enabled.
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-53

Release 13.0
December 2010

Advanced Modelling Options

UDFs for Multiphase Applications

When a multiphase model is
enabled, storage for properties and
variables is set aside for mixture as
phases.
well as for individual p

Domain ID = 1

Mixture Domain

Additional thread and domain data

structures required.
In general the type of DEFINE

macro determines which thread or

domain (mixture or phase) gets
passed to your UDF.
C_R(cell,thread) will return the
mixture density if thread is the
mixture thread or the phase
densities if it is the phase thread.
Numerous macros exist for data
retrieval.
ti
l

Customer Training Material

Phase 1
Domain

Phase 2
Domain

Interaction Domain

Mixture Thread

Phase
ase 3
Domain

Phase
Thread

Domain ID

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-54

Release 13.0
December 2010

Advanced Modelling Options

Heterogeneous Reaction Setup

Define

Customer Training Material

Phases

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-55

Release 13.0
December 2010

Appendix 2

Reacting Flow
Modeling
ANSYS, Inc. Proprietary
2009 ANSYS, Inc. All rights reserved.

9-56

April 28, 2009

Inventory #002600

Advanced Modelling Options

Eddy Dissipation Model (EDM)

Customer Training Material

Applicability
Flow Regime:
Chemistry:
Configuration:
C fi
ti

Turbulent flow (high Re)

Fast chemistry
P
Premixed
i d/N
Non-Premixed
P
i d/P
Partially
ti ll P
Premixed
i d

Application examples
Gas reactions
Coal combustion

Limitations
Unreliable when mixing and kinetic time scales are of similar order of magnitude
Does not predict kinetically-controlled intermediate species and dissociation effects.
Cannot realisticallyy model p
phenomena which depend
p
on detailed kinetics such as
ignition, extinction.

Solves species transport equations. Reaction rate is controlled by turbulent mixing.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-57

Release 13.0
December 2010

Advanced Modelling Options

Non-Premixed Model

Customer Training Material

Applicability
Flow Regime:
Chemistry:
Configuration:

Turbulent flow (high Re)

Equilibrium or moderately
q
((flamelet))
non-equilibrium
Non-Premixed only

Fuel
Reactor

Outlet

Oxidizer

Application examples
G
Gas reaction
ti (furnaces,
(f
burners).
b
) This
Thi iis usually
ll th
the model
d l off choice
h i if assumptions
ti
are
valid for gas phase combustion problems. Accurate tracking of intermediate species
concentration and dissociation effects without requiring knowledge of detailed reaction
rates (equilibrium).

Limitations
Unreliable when mixing and kinetic time scales are comparable
Cannot realisticallyy model p
phenomena which depend
p
on detailed kinetics ((such as
ignition, extinction).

Solves transport equations for mixture fraction and mixture fraction variance
(instead of the individual species equations)
equations).

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-58

Release 13.0
December 2010

Advanced Modelling Options

Premixed Combustion Model

Customer Training Material

Applicability
Flow Regime:
Chemistry:
Configuration:
C fi
ti

Turbulent flow (high Re)

Fast chemistry
P
Premixed
i d only
l

Fuel
+
Oxidizer

Reactor

Outlet

Application examples
Premixed reacting flow systems
Lean premixed gas turbine combustion chamber

Limitations
Cannot realistically model phenomena which depend on detailed kinetics (such as
ignition, extinction).

Uses a reaction progress variable which tracks the position of the flame front
(Zimont model).

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-59

Release 13.0
December 2010

Advanced Modelling Options

Partially Premixed Combustion Model

Customer Training Material

Applicability
Flow Regime:
Chemistry:
Configuration:
C fi
ti

Turbulent flow (high Re)

Equilibrium or moderately non-equilibrium (flamelet)
P ti ll premixed
Partially
i d only
l
Secondary
Fuel or Oxidizer

Application examples
Gas turbine combustor with dilution cooling holes
holes.
Systems with both premixed and non-premixed streams

Fuel
+
Oxidizer

Reactor

Outlet

Limitations
Unreliable when mixing and kinetic time scales are comparable.
Cannot realistically model phenomena which depend on detailed kinetics (such as
ignition, extinction).

In the partially premixed model, reaction progress variable and mixture fraction
approach are combined. Transport equations are solved for reaction progress
variable,
i bl mixture
i t
ffraction,
ti
and
d mixture
i t
ffraction
ti variance.
i
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-60

Release 13.0
December 2010

Advanced Modelling Options

Detailed Chemistry Models

Customer Training Material

The governing equations for detailed chemistry are generally stiff and difficult to
solve.
Tens of species
Hundreds of reactions
Large spread in reaction time scales.

Detailed kinetics are used to model:

Flame ignition and extinction

P ll t t (NOx,
Pollutants
(NO CO,
CO UHC
UHCs))
Slow (non-equilibrium) chemistry
Liquid/liquid reactions

Available Models:

Laminar finite rate

Eddy Dissipation Concept (EDC) Model
PDF transport
KINetics model (requires additional license feature)

CHEMKIN-format reaction mechanisms and thermal properties can be imported

directly.
g
in order to
FLUENT uses the In-Situ Adaptive Tabulation ((ISAT)) algorithm
accelerate calculations (applicable to laminar, EDC, PDF transport models).
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-61

Release 13.0
December 2010

Appendix 3

Moving and
Deforming Mesh
ANSYS, Inc. Proprietary
2009 ANSYS, Inc. All rights reserved.

9-62

April 28, 2009

Inventory #002600

Advanced Modelling Options

Absolute and Relative Velocities The Velocity Triangle

Customer Training Material

Absolute Velocity velocity measured w.r.t. the stationary frame

Relative Velocity velocity measured w.r.t. the moving frame.
The relationship between the absolute and relative velocities is given by
the Velocity Triangle rule:

In turbomachinery, this relationship can be illustrated using the laws of

vector addition
addition.
Blade motion

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-63

Release 13.0
December 2010

Advanced Modelling Options

Geometry Constraints for SRF

Customer Training Material

Single Fluid Domain

Walls and flow boundaries
Walls and flow boundaries (inlets and
outlets) which move with the fluid
domain may assume any shape.
Walls and flow boundaries which are
stationary (with respect to the fixed
frame) must be surfaces of revolution
about the rotational axis.
You can also impose a tangential
component of velocity on a wall
provided the wall is a surface of
revolution.

Shroud wall is stationary

(surface of revolution)

You can employ rotationally-periodic

boundaries if geometry and flow permit

Axis of
rotation

Hub, blade walls rotate

with moving frame

Advantage - reduced domain size

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-64

Release 13.0
December 2010

Advanced Modelling Options

SRF Set-up: Cell Zones

Customer Training Material

Use fluid BC panel to define rotational

axis origin and direction vector for
rotating reference frame
Direction vectors should be unit vectors
but Fluent will normalize them if they
arent

Select Moving Reference Frame as the

Motion Type for SRF
Enter Moving Frame Velocities
Rotational and Translational velocities
Rotation direction defined by right-hand
rule
Negative speed implies rotation in
opposite direction

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-65

Release 13.0
December 2010

Advanced Modelling Options

SRF Set-up: Boundary Conditions

Customer Training Material

Inlets:
Choice of specification of velocity vector or flow
direction in absolute or relative frames.
NOTE: Total pressure and temperature definitions
depend on velocity formulation!

Outlets
Static pressure or outflow.
Radial
R di l equilibrium
ilib i
option.
ti

Other Flow BCs

Periodics
Non-reflecting BCs
Target mass flow outlet

Walls
Specify walls to be
Moving with the domain
Stationary

NOTE: Stationary wall for Wall

M ti means stationary
Motion
t ti
w.r.t.
t the
th
cell zone!
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-66

Release 13.0
December 2010

Advanced Modelling Options

Interfaces

Customer Training Material

Fluid zones in multiple zone models

communicate across interface
boundaries.
Conformal interfaces
An interior mesh surface separates
cells from adjacent fluid zones.
Face mesh must be identical on either
side of the interface.

Non-conformal (NC) interfaces

Cells zones are physically disconnected
from each other.
Interface consists of two overlapping
surfaces (type = interface)
Fluent NC interface algorithm passes
fluxes from on surface to the other in a
conservative fashion (i.e. mass,
momentum, energy fluxes are
conserved).
d)
User creates interfaces using
Periodic repeat interface
DefineGrid Interfaces

Conformal interface

Interfaces may be periodic

C
Called
ll d periodic
i di repeatt iinterface.
t f
Require identical translational or
rotational offset.
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

Non-conformal interface

L9-67

Release 13.0
December 2010

Advanced Modelling Options

The MRF Model

Customer Training Material

The computational domain is divided into stationary and rotating fluid zones.
Interfaces separate zones from each other.
Interfaces can be Conformal or Non-Conformal.

Flow
Fl
equations
ti
are solved
l d in
i each
h fluid
fl id zone.
Flow is assumed to be steady in each zone (clearly an approximation).
SRF equations used in rotating zones.
At the interfaces between the rotating and stationary zones
zones, appropriate transformations
of the velocity vector and
velocity gradients are performed to compute
fluxes of mass, momentum, energy, and
other scalars.
scalars

MRF ignores the relative motions of the zones

with respect to each other.
Does not account for fluid dynamic
y
interaction
between stationary and rotating components.
For this reason MRF is often referred to as the
frozen rotor approach.

Ideally
Ideally, the flow at the MRF interfaces
should be relatively uniform or mixed out.
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-68

Mesh interface between

pump rotor and housing

Release 13.0
December 2010

Advanced Modelling Options

Geometric Constraints for MRF

Customer Training Material

Walls and flow boundaries which are contained within the rotating fluid zone
interfaces are assumed to be moving with the fluid zones and may assume any
shape.
Stationary walls and flow boundaries are allowed if they are surfaces of revolution
revolution.

The interface between two zones must be a surface of revolution with respect to
the axis of rotation of the rotating zone.
Periodic repeat interfaces are permitted but the periodic angles (or offsets) must be
identical for all zones.
stationary zone
rotating zone

Wrong!

Correct
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-69

Interface is not a
surface
or revolution
Release 13.0
December 2010

Advanced Modelling Options

Equations of Fluid Dynamics for Moving Frames

Customer Training Material

Equations of fluid dynamics can be transformed to a moving reference

frame with a choice of the velocities which are solved.
Relative Velocity Formulation (RVF)
Uses the relative velocity and relative total internal energy as the dependent
variables.

Absolute Velocity Formulation (AVF)

Uses the relative velocity and relative total internal energy as the dependent
variables.

Source terms appear in the momentum equations

for rotating frames.
Refer to Appendix for detailed listing of equations.
Relative formulation of x momentum equation:

Absolute formulation of x momentum equation:

Momentum
source terms

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-70

Release 13.0
December 2010

Advanced Modelling Options

MRF Set-Up

Customer Training Material

Generate mesh with appropriate

stationary and rotating fluid zones
Can choose conformal or nonconformal interfaces between cell
zones

For each rotating fluid zone

(Fluid BC), select Moving
Reference Frame as the Motion
Type and enter the rotational axis and
moving frame speed.
Identical to SRF except for multiple
zones
Stationary zones remain with
Stationary option enabled

Set up for BCs and solver settings

same as SRF.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-71

Release 13.0
December 2010

Advanced Modelling Options

MPM Set-up

Customer Training Material

Assign motion types and speeds to fluid zones and appropriate BCs for
each zone (like SRF).
Select upstream and downstream
zones which will comprise
mixing plane pair.
Upstream will always be
Pressure Outlet.
Downstream can be any inlet
BC type.

Set the number of Interpolation

P i t ffor profile
Points
fil resolution.
l ti
Should be about the same axial/radial
resolution as the mesh.

Mixing Plane Geometry determines

method of profile averaging.
Mixing plane controls
Under-relaxation Profile changes
g are under-relaxed from one iteration to the
next using factor between 0 and 1.
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.

L9-72

Release 13.0
December 2010

Advanced Modelling Options

SMM Set-up

Customer Training Material

Enable transient solver.

For moving zones, select Moving Mesh
as Motion Type in Fluid BC panel.
Define sliding zones as non-conformal
interfaces.
Enable Periodic Repeat option if
sliding/rotating motion is periodic.

Oth
Other BCs
BC and
d solver
l
settings
tti
are
same as the SRF, MRF models.
Run calculation until solution becomes
time-periodic

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-73

Release 13.0
December 2010

Advanced Modelling Options

Dynamic Mesh Setup

Customer Training Material

Enable transient solver.

Enable Dynamic Mesh model in
DefineDynamic Mesh.
Activate desired Mesh Methods and set
parameters as appropriate.
Define boundary motion in the Dynamic
Mesh Zones GUI.
UDF may be required.

Other models, BCs, and solver settings

are same as SMM models.
Mesh motion can be previewed using
SolveMesh Motion utility.

ANSYS, Inc. Proprietary

2010 ANSYS, Inc. All rights reserved.

L9-74

Release 13.0
December 2010