Fluent Tutorial
Fluent Tutorial
Gabriel Wȩcel
1
0.1 Air heater
This example shows how to simplify 3D geometry of the real object and perform
calculation in 2D space. Air heater geometry is given in Figure 1. The air stream
flows from left to the right. There are only two opening in the heater, inlet (left)
and outlet (right). The other sides (front, back, top, bottom) are insulated walls.
Inside the heater 3 cylindrical pipes are positioned perpendicular to the flow.
Data of the operation of the heater are given in Table 1.
The flow patern in such configuration of the heater in any cross section (aligned
with flow direction) is almost the same. The width of the heater is high enough
to neglect influence of the side walls (front and back) on the flow in the middle
part of the heater. Hence we can simulate flow in the heater with good accuracy
assuming 2D geometry (see the heater cross section in the Figure 2, we do not
utilize symmetry of the cross section geometry since we want later to analyze
different cylinder alignments). Even for simplified geometry flow over cylinders
emerge to be complex, with stagnation zones and revers flow close to the cylin-
ders. That feature require proper treatment of the mesh. First of all it needs
to be symmetric as the flow is symmetric. The best is to try generate fully
structured mesh and if possible with cells edges aligned with the direction of
the flow. This is not possible in all area of the flow, but at least at the cylinder
boundaries cells alignment should follow flow direction.
heaters
inlet outflow
2
m
04
0.
heaters
0.2 m
insulated walls
inlet outflow
3
Enter Radius = 0.02, Start Angle = 225, End Angle = 315
Press Apply
4
Face 1
Face 4
Face 2
Face 3
mesh faces around tube, first set distribution of the nodes on the edges
Mesh → Edge → Mesh Edges
Select with mouse (holding Shift button) edges creating tube
Deselect Grading
Select Spacing and enter Interval size = 0.001
Select Option Mesh
Press Apply
set distribution of the nodes on the edges radially connected with tube
Select with mouse (holding Shift button) radial edges
Select Grading
Select Type First Length and enter Length = 0.001
Select Spacing and enter Interval count = 20
Select Option Mesh
Press Apply
The mesh generated should have similar form of that shown in Figure 6, however
number of elements is different.
make 2 copies of mesh around the tubes
Geometry → Face → Move/Copy/Align
Select with mouse all faces
Select Copy and enter 2 (this is number of copies)
Select Operation → Translate
Enter X = 0.1, Y = 0, Z = 0
5
Figure 6: Mesh around the tube.
Press Apply
Coping of the faces in the way presented above result in double Edges lying at
the same position between copied meshes. If we did not set them as interfaces
Fluent will treat them as walls (no flow between these part of mesh). Simple
solution to this problem is connecting this edges.
As the result the double edges will be connected and one of them be deleted.
Remaining part of the mesh is generated by simply creating rectangular faces.
Since procedure is very simple only picture showing consequent steps is given
in Figure 8.
6
Figure 8: Steps of creating air heater geometry.
7
Enter, Name: heater 01
Select Type → WALL
Pick Entity : Edges, edge creating first pipe of the heater
press Apply
Read mesh file (mesh files have extension ṁsh) created in previous section.
File → Read → Case...
Define solver settings as default.
Define → Models → Solver...
Set turbulence modell
Define → Models → Viscous...
Select k − ² Standard turbulence model with option Standard Wall Function
Define material properties Define → Materials...
Check material properties for air
Enter Density (kg/m3) equal to 1.225
Enter Cp (j/kgK) equal to 1006.43
8
Enter Thermal Conductivity (W/mK) equal to 0.0242
Enter Viscosity (kg/ms) equal to 1.7894e-5
Confirm changes pressing Change Create
Close Material panel pressing Close
Define boundary condition Define → Boundary Conditions...
Select Zone → inlet, press Set
Enter Velocity Magnitude (m/s) equal to 5
Enter Temperature (K) equal to 300
Select Turbulence Specification Method → Intensity and Hydraulic Diameter
Enter Turbulence Intensity (%) equal to 10
Enter Hydraulic Diameter (m) equal to 0.4
Accept settings pressing OK
Select Zone → heater 01, press Set
From Thermal tab for Thermal Conditions select Heat Flux
Enter Heat Flux (W/m2) equal to 16000
Accept settings pressing OK
Press Copy from Boundary Condition panel
Under From Zone select heater 01
Under To Zone select heater 02, heater 03
Press Copy
Accept selection pressing OK
Close Copy BCs panel pressing Close
Close Boundary Conditions panel pressing Close
set up solver parameters
Solve → Controls → Solution...
assume all default settings
Initialize solution Solve → Initialize → Initialize...
Press Init and close Solution Initialization panel
Set solution monitoring option Solve → Monitors → Residual...
Under Option select Plot
For Residual → continuity → Convergence Criterion enter value equal to 10e-9
Accept settings pressing OK
save Fluent settings parameter in case file (case files have extension .cas) File
→ Write → Case..., enter file name and accept settings pressing OK
Type in Fluent command window it 100, (this command executes 100 iterations)
Observe in Fluent result window residuals of the solved equations
9
Post processing
All calculated variable can be displayed in the Fluent result window in the form
of colored field. Additionally profiles of these variables can be created at arbi-
tral position inside computational domain. The result window can be zoomed
in and out in order to observe particular regions of the flow. Displayed variable
scale can be adjusted to arbitral ranges.
Surface → Line/Rake...
From Type select Rake, advantage of using Rake instead of Line is that we prede-
fine number of points at which values are plotted or printed, in case of using Line
number of points for dense meshes can be large, additionally Line extract values
from the closes volume cell center, when Rake interpolate values from closes vol-
ume cells and calculate it for position at which point of a rake is placed.
For Number of Point enter required value, we suggest at least 20, but that strongly
depend on the problem solved, variable printed and position of the Rake
From Points enter,
x0(m) = 0.1, x1(m) = 0.5
y0(m) = 0, y1(m) = 0.2
For New Surface Name enter desired name of the rake and press Create
Close Line/Rake Surface panel pressing Close
Created Rake can be used from Plot → XY Plot... panel in order to plot, print
or write to a file solution variables placed on the defined by rake positions.
10
discussion of the near wall treatment and simply standard wall approach is
utilized. However main purpose of this tutorial is to show how to relatively easy
create structured mesh for the geometry which automatic mesh generator are
not able to handle such meshes.
One can try to build geometry with different heater alignments and observe how
that influence, the flow, pressure drop, and temperature profile at the outlet of
the heater.
0.2 Cyclone.
In many industrial processes emerge a need of cleaning gases from dispersed
inert particles suspended within gas (eg. removal of flying ash from flue gases
in industrial coal fired boilers). Easiest and most commonly used method of
separation takes advantage of gravitation forces. Device which work on this
basis is a cyclone. With this tutorial we build simple cyclone geometry, mesh it,
and run Fluent simulations. Instead of flue gases we will use air stream polluted
with ash. Data for boundary conditions are given in Table 2.
11
0.6 m
0.2 m
0.2 m
0.5 m
0.8 m
0.2 m
m
0.7
0.2 m
1.0 m
0.2 m
0.05 m
0.2 m
0.4 m
12
0.6 m
0.2 m
0.5 m
m
0.2 m 0.7
0.6 m
0.2 m
0.8 m
1.0 m
0.2 m
0.05 m
0.2 m
0.2 m
0.4 m
press Apply
13
Check Move, Translate
Enter X = 0, Y = 0, Z = 1.55
press Apply
14
Figure 11: Faces to be selected for Face Connect oparation.
15
Zones → Specify Boundary Types
Check (Add)
Enter, Name: ash
Select Type WALL
Pick Entity : Faces, faces creating ash hopper, see Figure 12
press Apply
outlet
velocity inlet
ash hopper
16
procedure for meshing cyclone geometry.
Mesh → Face
Pick Faces, select all the faces
Select Elements: Tri
Select Type: Pave
Check Spacing: Apply
Enter Interval size 0.05
Press Apply
Mesh → Volume
Pick Volumes, select all the volumes
Select Elements: Tet/Hybrid
Select Type: Tgrid
Uncheck Spacing: Apply
Press Apply
Final mesh should contain around 20 000 cells. The last task to perform in
Gambit is to export generated mesh to the file.
File → Export → Mesh
Press Browse to select destination folder.
Enter name of the file, extension will be given by default.
Press Accept
Read mesh file (mesh files have extension ṁsh) created in previous section.
File → Read → Case...
Define solver settings as default.
Define → Models → Solver...
Set turbulence modell
Define → Models → Viscous...
Select k − ² RNG turbulence model with option Swirl Dominated Flow
In the Discrete Phase Model panel change Maximum Number of Steps to 10000
Set Injections
Select Injection Type → surface
Select Release From Surfaces → in (inis an inlet face)
Select Material → ash
Select Diameter Distribution → rosin-rammler-logarithmic
Select tab Point Properties
Enter Total Flow Rate (kg/s) equal to 0.001
Enter Min. Diameter (m) equal to 1e-6
Enter Max. Diameter (m) equal to 300e-6
17
Enter Mean Diameter (m) equal to 150e-6
Enter Spread Parameter equal to 2.8
Enter Number of Diameters equal to 15
Select tab Turbulent Dispersion
From Stochastic Tracking select Discrete Random Walk Model
Enter Number of Tries equal to 5
Accept settings pressing OK
Define material properties Define → Materials...
Change density for air
Enter Density (kg/m3) equal to 1.094
Confirm changes pressing Change Create
Change density for inert-particle ash
Enter Density (kg/m3) equal to 2100
Confirm changes pressing Change Create
Define operating condition Define → Operating Conditions...
Select Gravity
Enter gravitation acceleration Z (m/s2) equal to 9.81
Accept settings pressing OK
Define boundary condition Define → Boundary Conditions...
Select Zone → in, press Set
Enter Velocity Magnitude (m/s) equal to 7.98
Select Turbulence Specification Method → Intensity and Hydraulic Diameter
Enter Hydraulic Diameter (m) equal to 0.2
Accept settings pressing OK
Select Zone → ash, press Set
Select tab DPM
Under Discrete Phase Model Condition select Boundary Cond. Type → trap
Accept settings pressing OK
Select Zone → wall, press Set
Select tab DPM
Under Discrete Phase Reflection Condition select Normal → constant
Enter value equal to 0.8
Under Discrete Phase Reflection Condition select Tangent → constant
Enter value equal to 0.8
Accept settings pressing OK
Close Boundary Condition panel
set up solver parameters Solve → Controls → Solution...
From Discretization select Momentum → Second Order Upwind
From Discretization select Turbulent Kinetic Energy → Second Order Upwind
From Discretization select Turbulent Dissipation Rate → Second Order Upwind
Accept settings pressing OK
Initialize solution Solve → Initialize→ Initialize...
Press Init and close Solution Initialization panel
Set solution monitoring option Solve → Monitors → Residual...
Under Option select Plot
For Residual → continuity → Convergence Criterion enter value equal to 10e-9
18
Accept settings pressing OK
save Fluent settings parameter in case file (case files have extension .cas) File
→ Write → Case..., enter file name and accept settings pressing OK
Type in Fluent command window it 100, (this command executes 100 iterations)
Observe in Fluent result window residuals of the solved equations
The cyclone axis is aligned with Z axis and crossing X = 0 and Y = 0 cartesian
coordinates. Now we create plane crossing cyclone for Y = 0.
19
For New Surface Name enter desired name of the surface and press Create
Close Plane Surface panel pressing Close
You can repeat procedure above to create more planes in the arbitral positions
inside analyzed domain. You can see created planes by displaying them in the
Fluent result window.
20
Lagrangian model. See below for executing tracing procedure.
The most interesting is Mass Transfer Summary which shows mass fluxes of
Incomplete, Trapped and Escaped particle streams. General report shows only
number of incomplete, trapped and escaped particle stream. Sometimes even
21
large number of escaped particle streams not impose low cyclone efficiency,
because these streams could be low diameter particle streams. Hence in order
to asses cyclone efficiency mass fluxes of trapped and escaped particle streams
needs to be compared. Escaped particle stream number indicate how many trace
of the particle streams has not been completed. There are neither trapped nor
escaped and traced has been finished inside domain. If number and mass flux
of incomplete stream is large we need to increase Max. Number of Steps under
Discrete Phase Model panel opened from Define menu.
22