Fluent 13.0 workshop03-Multiple-Species PDF

Customer Training Material
W k h 3
Workshop
Modeling
g Multi-Species
p
Flow
Introduction to
ANSYS FLUENT
ANSYS, Inc. Proprietary
2010 ANSYS, Inc. All rights reserved.
WS3-1
Release 13.0
December 2010
WS3: Modelling Multi Species Flow
Introduction (1)
In this workshop you will analyze the release of heat and combustion
gases from a single car with an engine fire in a ventilated parking garage.
The simulation will be run as a steady state analysis assuming the fire
has reached a stable developed stage.

WS3-2
Release 13.0
December 2010
Introduction (2)
Simulation Physics & Boundary Conditions
Mixture of N2, O2, CO2 and H20

6 m/s exhaust
80 N/m3 momentum source in jets
0.1 kg/s flue gases (H2O and CO2) at 1200 K
k
k-epsilon
il turbulence
b l
model
d l
DO radiation model with WSGGM for absorption coefficient
300 K inlet temperature
Fresh Air Inlet
pressure_inlet
Symmetry Plane
symmetry
Jet Fan
fluid_jet
Heat and gas release from
engine bay
mass-flow-inlet-flue
mass flow inlet flue_gas
gas
Air Outlet
inlet_exhaust
WS3-3
Release 13.0
December 2010
Goals
This Workshop will demonstrate how to:

Set up a case with multiple species
Define absorption coefficient for flue gases
Define gravity and operating conditions for buoyancy driven flows
Define volumetric sources
Ch
Change
solver
l
settings
i
and
d boundary
b
d
conditions
di i
d i the
during
h calculation
l l i
by using calculation activities
Judge solution convergence using monitors
Postprocess a case with multiple species

WS3-4
Release 13.0
December 2010
Preparation
Files
Copy the following files to your working directory.
garage.msh
garage msh
calc_activities.jou
multi-species-solved.wbpz

WS3-5
Release 13.0
December 2010
Starting Workbench
Open Workbench
Start > Programs > ANSYS 13.0 > ANSYS Workbench
Drag a FLUENT Component System into the project schematic
Rename the system to Garage (RMB on cell A1 to rename the
system)

WS3-6
Release 13.0
December 2010
User Files Directory
Manage External Files

This project requires the import of existing files. Workbench creates a
directory named User
User Files
Files where files created externally can be
stored and associated with the project.
This enables Workbench to synchronise projects with external files.
From the Workbench project page save the project to your working
directory using the name multi-species.
The user_files directory will be found in the project files directory.
Copy the two files garage.msh and calc_activities.jou to this

directory. All files associated with the project are listed when the
Files option is checked from the Workbench View menu.

WS3-7
Release 13.0
December 2010
Starting FLUENT in Workbench
Start FLUENT
Double click on the Setup cell to open FLUENT
Choose 3D and Double Precision under Options and retain the
other default settings (if your computer has 2 or more nodes and
parallel licenses are available, you also could start FLUENT
parallel)

WS3-8
Release 13.0
December 2010
Import Mesh and General setup
Import the existing mesh file

Under the Fluent File menu select Import >
Mesh.
Select the file garage.msh and click OK to
import the mesh.
To check for any problems in the mesh click
Check. There should be no problems
reported in the TUI window.
Reorder the mesh using Mesh > Reorder >
Domain (from the menu).
Reordering the domain can improve the
computational performance of the solver by
rearranging the nodes, faces and cells in memory.
Retain defaults for the solver

Enable gravity and set z = -9.81 m/s2
WS3-9
Release 13.0
December 2010
Setting Physics (Turbulence)
Specify turbulence model

Select Models in the navigation pane.
Double click Viscous in the model
selection pane. The Viscous Model
panel will open.
Select k-epsilon (2 eqn) under Model
Model,
Standard under k-epsilon Model and
Enhanced Wall Treatment under NearWall Treatment.
Turbulence modelling, as with all physics
modelling, is a complex area. There are
many application specific options. The
standard k-epsilon model is a simple but
robust model. The documentation provides
further guidance on which models to use for
specific applications.
Click
Cli k OK
WS3-10
Release 13.0
December 2010
Setting Physics (Radiation)
Specify radiation model

Double click Radiation in the model selection pane. The Radiation Model panel
will open.
Select Discrete Ordinates (DO) under Model.
Set Flow Iterations per Radiation Iteration to 1 and click OK.
ANSYS FLUENT provides five radiation models. The DO model spans the entire range of
optical
ti l thicknesses,
thi k
and
d allows
ll
you tto solve
l problems
bl
ranging
i ffrom surface-to-surface
f
t
f
radiation
di ti
to participating radiation in combustion problems.
You will be informed, that material properties or
methods have changed. The Text User Interface
(TUI) will list new material properties required for
radiation modelling. Click OK. Material
properties will be defined later
later.
Note also that the energy equation is switched on

automatically when any heat transfer related
physics are requested.

WS3-11
Release 13.0
December 2010
Setting Physics (Species / Mixture)
Specify species model

Double click Species in the model
selection pane. The Species Model
panel will open.
In this workshop the products of
combustion (heat & gases) will be
modelled rather than the reaction itself.
Select Species Transport and click

OK.
Switching on the species model will, as
with the radiation model, introduce new
material properties. An information
box will appear
appear. Click OK to accept
this.
This setup will enable the tracking of

non-reacting
non
reacting chemical species
species.
WS3-12
The species model requires the definition

of a mixture representing the chemical
species of interest. The default mixture
contains nitrogen, oxygen and water
vapour.
Release 13.0
December 2010
Defining Materials (1)
Specify mixture
Select Materials in the navigation
pane.
Note that under Mixture in the Materials
pane the default mixture is listed as
containing nitrogen, oxygen and watervapour.
vapour
Double click mixture-template, this will

open the Create/Edit Materials panel
with the mixture preselected
preselected.

WS3-13
Release 13.0
December 2010
Add a mixture species

In the Create/Edit Materials panel click
on FLUENT Database.
The FLUENT Database Materials panel will
appear.
Select Fluid as Material Type.

All predefined fluids materials will be listed
under FLUENT Fluid Materials.
Select carbon-dioxide (co2).

Click Copy.
The TUI will inform you that carbon-dioxide
has been added.
Cl
Close the
th FLUENT D
Database
t b
M
Materials
t i l
panel.

WS3-14
Release 13.0
December 2010
Add a mixture species

Copying carbon-dioxide(co2) from the fluid database has made the
species available to the simulation, now add it to the mixture.
In the Create/Edit Materials panel ensure Material Type is set to Mixture and click
the Edit button for the Mixture Species.
In the Species panel select co2 from the Available Materials list and use the Add
b tton to add it to the Selected Species list
button
list. The Selected Species defines the
component species of the mixture.
The order of the species listed under Selected Species is important. The most
abundant species should be listed last. Use the Remove button to remove n2
followed by the Add button to replace n2 as the last species. Click OK.

WS3-15
Release 13.0
December 2010
Specify mixture
Modify the existing settings for Thermal Conductivity, Viscosity and
Absorption Coefficient as shown below. Click Change/Create to apply
the changes.

WS3-16
Release 13.0
December 2010
Specify solid material

Select Solid under Material Type (The default solid material is aluminum(al)).
Modify aluminum (Name,
(Name Chemical Formula & Properties) as shown
below.
Click Change/Create and choose No for overwriting. Close the
Create/Edit Materials panel
panel.
Selecting No preserves the original material (aluminum) and adds the new modified
material.

WS3-17
Release 13.0
December 2010
Cell Zone Conditions (1)
Set cell zone conditions
Select Cell Zone Conditions from the Navigation Pane.

Double click fluid_jet
fluid jet from the Zone list
list.
Activate Source Terms.
Select the Source Terms tab and click on Edit for Y Momentum.
Add 1 source, select
l
C
Constant and
d enter a value
l off -80
80 ((n/m
/ 3)
Click OK in both panels.

WS3-18
Release 13.0
December 2010
Boundary Conditions (exhaust 1)
Set boundary conditions
Select Boundary Conditions from the Navigation Pane.

Double click inlet_exhaust
inlet exhaust from the Zone list
list.
Apply Momentum settings as shown.
Click OK.

WS3-19
Release 13.0
December 2010
Boundary Conditions (inlet 1)

Double click pressure_inlet from the Zone list.
Apply Momentum and Thermal settings as shown
shown.
Continued on next page..

WS3-20
Release 13.0
December 2010
Boundary Conditions (inlet 2)

Apply Radiation and Species settings as shown.
Click OK
OK.

WS3-21
Release 13.0
December 2010
Boundary Conditions (smoke release 1)

Double click mass-flow-inlet-flue_gas from the Zone list.
Apply Momentum and Thermal settings as shown
shown.
Continued on next page..

WS3-22
Release 13.0
December 2010
Boundary Conditions (smoke release 2)

Apply Species settings as shown.
Click OK
OK.

WS3-23
Release 13.0
December 2010
Boundary Conditions (walls 1)

Double click wall from the Zone list.
Apply Thermal settings as shown
shown.
Click OK.

WS3-24
The Shell Conduction option enables thin

walls to solve for heat transfer in both the
normall and
d planar
l
di
directions
ti
without
ith t th
the
need to volume mesh them.
Release 13.0
December 2010
Boundary Conditions (walls 2)

Select Copy from the Boundary Conditions Pane.
Select wall in the From Boundary Zone list and wall_ceiling
wall ceiling and
wall_floor in the To Boundary Zones list.
Click Copy, click OK in the question dialog box then Close.
This will copy all boundary settings from the boundary zone wall to both wall_ceiling
wall ceiling
and wall_floor.

WS3-25
Release 13.0
December 2010
Operating Conditions
Set operating conditions

Select Operating Conditions.
Apply Specified Operating Density settings as shown
shown.
Click OK.

WS3-26
Release 13.0
December 2010
Operating Conditions
Notes
ANSYS FLUENT avoids the problem of roundoff error by subtracting the operating
pressure (generally a large pressure roughly equal to the average absolute pressure in
th flow)
the
fl ) from
f
the
th absolute
b l t pressure, and
d using
i th
the resultlt (t
(termed
d th
the gauge pressure).
)
The absolute pressure is simply the sum of the operating pressure and the gauge
pressure.
Operating temperature is only used when using the Boussinesq density model, so in this
case, it has no meaning.
Operating density is also a value for avoiding roundoff errors.
errors For simulations where
pressure boundary conditions are present it is important to set the value correctly
otherwise the pressure at the boundary will be incorrect and may lead to unphysical flow
conditions. Here you have to set it to the density for the conditions at the pressure-inlet a gas at 293
293.15
15 K with 23% O2 and 67% N2. You can initialize your flow field with these
conditions to get the value for the operating density from the postprocessor (e. g.
Reports -> Volume Integral). See the Users Guide Natural Convection and BuoyancyDriven Flows for more details.

WS3-27
Release 13.0
December 2010
Solution Methods
Set solution methods

Select Solution Methods from the Navigation Pane.
Under Spatial Discretization set Pressure to Body Force Weighted and
Energy to Second Order Upwind.
The Body Force Weighted scheme is recommended for problems involving large
body forces.
The Second Order Upwind scheme is generally recommended for final solutions
(first order is less accurate but numerically more stable and is often used for initial
approximate solutions).

WS3-28
Release 13.0
December 2010
Solution Controls
Set solution controls

Select Solution Controls from the Navigation
Pane.
Set Under-Relaxation Factors as listed below.
Lower Under-Relaxation Factors will reduce the
solution change between iterations leading to more
stable convergence though requiring more iterations
to reach convergence.
Pressure = 0.3
Density = 1
Body Forces = 1
Momentum = 0.7
Turbulent Kinetic Energy = 0.5
T b l t Di
Turbulent
Dissipation
i ti R
Rate
t =0
0.5
5
Turbulent Viscosity = 0.7
h2o = 1
o2 = 1
co2 = 1
Energy = 1
Discrete Ordinates = 1

WS3-29
Release 13.0
December 2010
Monitors (1)
Set solution monitors

Select Monitors from the Navigation Pane.
Click Edit and set the Residual Monitors as shown below
below.
By default ANSYS FLUENT will plot residuals to the window and print to the console. The default setting for the
convergence criterion is Absolute which means that the solver will continue until all residuals fall below the
Absolute Criteria values specified in the Equations box. Switching the Convergence Criterion to none will cause
the solver to continue until a maximum number of iterations is reached. Click OK.

WS3-30
Release 13.0
December 2010
Monitors (2)
Set surface monitors

Under Surface Monitors, click Create.
Residuals monitor the numerical accuracy of the solution. However, it is equally important to ensure that
solution variables have converged to sensible stable values. Creating Surface Monitors enables solution values
of interest to be monitored on specific surfaces within the domain. Here, the mass flow rate will be monitored on
the surface interior_jet_in. The monitor will be plotted to window 2 and printed to the console.
Set the Surface Monitor as shown below and click OK.

WS3-31
Release 13.0
December 2010
Monitors (3)
Set surface monitors

Create another Surface Monitor.
There are many different types of calculations that can be performed over surfaces listed under Report Type.
This time we will monitor the integral of the Total Surface Heat Flux on the surface wall_floor plotting to
window 3 and printing to the console.
Set the Surface Monitor as shown below and click OK.

Created Surface Monitors are listed after creation.

WS3-32
Release 13.0
December 2010
Monitors (3)
Notes on monitors
Buoyancy driven flows often show transient behavior. For this reason,
the residuals often will oscillate. Because of this, convergence should
always be judged by solution variable monitors and flux reports. The
residuals will however show overall for convergence behavior and
stability.
y
In cases of an oscillating steady state solution, a common approach is
to continue the simulation in transient mode
mode. In many cases the
oscillations will reduce significantly after a few time steps.
The use of surface/volume monitors combined with residuals will
provide the best overall judge of solution convergence.

WS3-33
Release 13.0
December 2010
Solution Initialization
Set initialization values

Select Solution Initialization from the Navigation Pane.
The iterative process requires initial values for all solution variables. The default (standard) initialization method
applies single values for each variable globally throughout the domain. These values should be within a
reasonable range of what is expected in the domain. Specifying unrealistic values could cause convergence
problems or require more iterations. Values can be specified manually or can be calculated from any or all
boundaries using the Compute from option.
Specify the values manually as follows.
Initial Values
Gauge Pressure = 0
X Velocity = 0
Y Velocity = 0
Z Velocity = 0
Turbulent Kinetic Energy = 0.01
Turbulent Dissipation Rate = 0.01
h2o = 0
o2 = 0.23
co2 = 0
Temperature = 300

WS3-34
Release 13.0
December 2010
Calculation Activities (1)
Changing settings during the calculation
Although simulations can often be run to convergence without

adjusting the solver controls it is often beneficial to change these
settings during the calculation. For example, some equations may not
need as many iterations to converge. Also, in many thermal cases, it
is advantageous
g
in terms of solution stability
y to solve the momentum
equations first before switching on the energy equation. This step
illustrates how to use the Case Modification panel to adjust solver
settings during a calculation.
In this case we will first solve the full equation set for a small number
of iterations,, then onlyy solve for mass/momentum before finallyy
solving the full set again.

WS3-35
Release 13.0
December 2010
Select Calculation Activities from the Navigation Pane, check

the Automatically Initialize and Modify Case box and click Edit.
Select the Case Modification tab.
Using the defaults would result in the simulation being run

using the existing settings (Original Settings) for 1 iteration.
We will
W
ill be
b running
i using
i th
the original
i i l settings
tti
th
then apply
l 2
modifications. Increase Defined Modifications to 2.

WS3-36
Release 13.0
December 2010
Commands are entered using the Text User Interface (TUI) syntax. Options selected using the
Graphical User Interface (GUI) have equivalent TUI commands.
Rather than type these out a journal file is provided for this workshop that will enter these for you.
From the main menu, File>Read>Journal, Select the journal file calc_activities.jou from the user_files
directory. Answer No to the question dialog box asking if there are any read/write commands in the
journal.
Modification 1 will switch off all equations except the flow equations and run for 100 iterations.
Modification 2 will switch all equations back on and run for 20 iterations.
Set the Number of Iterations for Original Settings to 20.
Click OK and answer No to the question dialog box

box.

WS3-37
Release 13.0
December 2010
The modifications now contain the TUI commands for switching various
equations
q
on and off.
For example, the TUI command for switching off the energy equation is of the
form; (substituting n for y to switch on)
/solve/set/equations
q
temperature
p
n
The list of available TUI commands is extensive and goes beyond options
available in the GUI.
All aspects of the solver can be controlled using

g TUI commands.
The journal used here contained the TUI commands to setup up the 2 case
modifications.
Consult the documentation for a more comprehensive description of TUI

command usage.

WS3-38
Release 13.0
December 2010
Run Calculation
Start the calculation

Select Run Calculation from the Navigation Pane.
The Number of Iterations is set to Automatic since this has been defined in the Case
Modification panel.
Before starting the calculation it is good practice to save the

project.
File>Save Project.
Project
The calculation will take approximately 40 minutes in serial.
However if you are short of time, we have supplied a converged

model that you can use instead (see next slide).

WS3-39
Release 13.0
December 2010
Reading in an Solution computed elsewhere
There are times when you may want to bring a solution (perhaps made on a remote
cluster) into your workbench session.
1. Stop the iterations. Close the existing FLUENT session, and return to Workbench
2. Drag a new FLUENT session into the Workbench page. Rename as
Imported_Result.
3. In the Setup Cell (B2), right click
and
d select
l t I
Importt FLUENT case
Pick the file workshop3converged.cas.gz
4 FLUENT will
4.
ill start
t t up be
b sure to
t
select double precision. Once
open, close this session again.
5 In
5.
I the
th Solution
S l ti cellll (B3)
(B3), right
i ht click
li k
and Import Final Data, selecting the
file workshop3-converged.dat.gz.
6 D
6.
Double
bl click
li k on th
the solution
l ti cellll (B3)
and continue with this workshop.
WS3-40
Release 13.0
December 2010
Check Residuals and Monitors (1)
Check Monitors, Residuals and Fluxes
If progressing from the pre solved project select Monitors from the Navigation Pane, click
Edit under Residuals, Statistics and Force Monitors and click Plot in the Residual Monitors
Panel that appears.
The residuals are still quite high after 140 iterations and the surface monitors
are still changing.

WS3-41
Release 13.0
December 2010
Check Residuals and Monitors (2)
Check Monitors, Residuals and Fluxes
Select Reports from the Navigation Pane and select

Fluxes under Reports.
p
Click Setup.
p
In the Flux Reports panel check the Mass balance

over all boundaries as shown below. Click compute.
Note that the Net Results indicate a relatively large imbalance.

Thi iis an iindication
This
di ti th
thatt th
the solution
l ti iis nott ffully
ll converged.
d N
Nett
heat flux can be calculated in the same way. Close the panel.

WS3-42
Release 13.0
December 2010
Exit
Close Fluent
Close ANSYS FLUENT by selecting File>Close FLUENT
ANSYS FLUENT contains basic built in post processing capabilities which can
be used to quickly assess results graphically and numerically during and after
the solution.
solution
CFD POST is a powerful separate post processing application containing many

more advanced features. The remainder of this workshop will introduce some
features of CFD POST.

WS3-43
Release 13.0
December 2010
Starting CFD POST
Add a Results System
Drag a Results System onto the Solution cell of the FLUENT system.
This will link a CFD POST (Results) system to the FLUENT system solution
data as shown.
Start CFD POST by double clicking on cell B2.
This will start CFD POST and load the FLUENT solution data.

WS3-44
Release 13.0
December 2010
Post Processing Wall Temperature (1)
Add a Location representing a group of surfaces
For displaying the temperature on all walls the most effective way is to use a
surface group.
g p
Select Locations > Surface Group, enter a name as shown below and click OK.

WS3-45
Release 13.0
December 2010
Post Processing Wall Temperature (2)
Define Location details
The details of the new location

will be displayed
p y in the bottom left
pane.
Select Locations, click on and

select all walls (CTRL click to
multiple
lti l select).
l t)
Click OK in the Location Selector.

WS3-46
Release 13.0
December 2010
Postprocessing Wall Temperature (3)
Apply contour to Location
Select the Color tab, click on Mode and select Variable Temperature and click
on Apply.
pp y

WS3-47
Release 13.0
December 2010
Postprocessing Wall Temperature (4)
Make changes to the legend
If you want to change the appearance of the legend, select it in the tree (double
click)) and change
g the details accordingly
g y for example,
p set Precision to Fixed
instead of Scientific.

WS3-48
Release 13.0
December 2010
Postprocessing Section Plane (1)
Add a Location representing a section plane
Select Locations > Isosurface, enter a name as shown below and click OK.

WS3-49
Release 13.0
December 2010
Postprocessing Section Plane (2)
Define Location details
Under Geometry select Y as Variable and set the value to 9 m and set Color
and Render as shown below.
Click Apply. Switch off the previous Location (Walls) by unchecking it in the
Outline.

WS3-50
Release 13.0
December 2010
Postprocessing 3d Isosurface
Add a 3d Isosurface representing gases
Use the same procedure - select Locations > Isosurface, enter a name (Gas)
and click OK.

WS3-51
Release 13.0
December 2010
Postprocessing Streamlines (1)
Add 3d streamlines to visualize flow
Click the Streamline button
Apply the Details as shown and click Apply.

and enter a name (Streamlines).
WS3-52
Release 13.0
December 2010
Postprocessing Volume Rendering
Add Volume Rendering
Click the Volume Rendering button

Volume).
)
Apply the settings shown below under the Geometry tab, switch off the other
Locations and click Apply.

WS3-53
and enter a name (Gas
Release 13.0
December 2010
Save and Close
Other Features
There are many other features available in CFD POST
Animations
Expressions & Functions
Charts & Tables
Automatic Report Generation
Multiple Case Comparison
Close CFD POST (File > Quit) and save the project from the Workbench Project
Page.

WS3-54
Release 13.0
December 2010

Fluent 13.0 workshop03-Multiple-Species PDF

Uploaded by

Copyright:

Available Formats

Fluent 13.0 workshop03-Multiple-Species PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Fluent 13.0 workshop03-Multiple-Species PDF

Uploaded by

Copyright:

Available Formats

Customer Training Material

WS3: Modelling Multi Species Flow

Customer Training Material

ANSYS, Inc. Proprietary

WS3: Modelling Multi Species Flow

Customer Training Material

Simulation Physics & Boundary Conditions

Mixture of N2, O2, CO2 and H20

WS3: Modelling Multi Species Flow

Customer Training Material

This Workshop will demonstrate how to:

ANSYS, Inc. Proprietary

WS3: Modelling Multi Species Flow

Customer Training Material

ANSYS, Inc. Proprietary

WS3: Modelling Multi Species Flow

Customer Training Material

ANSYS, Inc. Proprietary

WS3: Modelling Multi Species Flow

User Files Directory

Customer Training Material

Manage External Files

Copy the two files garage.msh and calc_activities.jou to this

ANSYS, Inc. Proprietary

WS3: Modelling Multi Species Flow

Starting FLUENT in Workbench

Customer Training Material

ANSYS, Inc. Proprietary

WS3: Modelling Multi Species Flow

Import Mesh and General setup

Customer Training Material

Import the existing mesh file

Retain defaults for the solver

WS3: Modelling Multi Species Flow

Setting Physics (Turbulence)

Customer Training Material

Specify turbulence model

WS3: Modelling Multi Species Flow

Setting Physics (Radiation)

Customer Training Material

Specify radiation model

Note also that the energy equation is switched on

ANSYS, Inc. Proprietary

WS3: Modelling Multi Species Flow

Setting Physics (Species / Mixture)

Customer Training Material

Specify species model

Select Species Transport and click

This setup will enable the tracking of

The species model requires the definition

WS3: Modelling Multi Species Flow

Defining Materials (1)

Customer Training Material

Double click mixture-template, this will

ANSYS, Inc. Proprietary

WS3: Modelling Multi Species Flow

Defining Materials (2)

Customer Training Material

Add a mixture species