Intermediate Lab 1 Manual
Intermediate Lab 1 Manual
Intermediate Lab 1 Manual
1. Purpose
The Purpose of CFD Lab 1 is to simulate steady laminar and turbulent pipe flow following the
CFD Process by an interactive step-by-step approach. Students will have hands-on experiences
using ANSYS to compute axial velocity profile, centerline velocity, centerline pressure, and
friction factor. Students will conduct verification studies for friction factor and axial velocity
profile of laminar pipe flows, including iterative error and grid uncertainties and effect of
refinement ratio on verification. Students will validate turbulent pipe flow simulation using
EFD data, analyze the differences between laminar and turbulent flows, and present results in
CFD Lab report.
Geometry
Pipe (ANSYS
Design Modeler)
Laminar
Turbulent
Physics
Mesh/Grid
Solution
Results
General (ANSYS
Fluent - Setup)
Structure
(ANSYS Mesh)
Plots (ANSYS
Fluent- Results)
Model (ANSYS
Fluent - Setup)
Uniform
(ANSYS Mesh)
Solution
Methods
(ANSYS Fluent
- Solution)
Boundary
Conditions
(ANSYS Fluent Setup)
(ANSYS Mesh)
Non-uniform
Monitors
(ANSYS Fluent Solution)
Reference Values
(ANSYS Fluent Setup)
Solution
Initialization
(ANSYS Fluent Solution)
Graphics and
Animations
(ANSYS FluentResults)
2. Simulation Design
In CFD Lab 1, simulation will be conducted for laminar and turbulent pipe flows. Reynolds
number is 655 for laminar flow and 111,569 for turbulent pipe flow, respectively. The schematic
of the problem and the parameters for the simulation are shown below.
Parameter
Radius of Pipe
Diameter of Pipe
Length of the Pipe
Unit
m
m
m
Value
0.02619
0.05238
7.62
Uniform Grid
Non-uniform Grid
Symmetry Axis
Inlet
Outlet
Velocity Profile
Pipe Wall
Since the flow is axisymmetric we only need to solve the flow in a single plane from the
centerline to the pipe wall. Boundary conditions need to be specified include inlet, outlet, wall,
and axis, as will be described details later. Uniform flow was specified at inlet, the flow will
reach the fully developed regions after a certain distance downstream. No-slip boundary
condition will be used on the wall and constant pressure for outlet. Symmetric boundary
condition will be applied on the pipe axis. Uniform grids will be used for the laminar flow
whereas non-uniform grid will be used for the turbulent flow.
Table 2 - Grids
Grid/Mesh
Grid/Mesh
Type
# of Divisions
453
45
319
32
226
23
113
11
80
56
28
564
15
Uniform
0
T
Nonuniform
Experimental, analytical results and simulation results will be compared. Additionally, detailed
verification and validation study will be conducted. All the studies are detailed in the Table 3. In
this manual, detailed instructions are given for the turbulent flow simulation and laminar flow
simulations using non-uniform grid and uniform grid 8 respectively. The grid and simulation
setups have been provided for rest of the simulation with workbench on the class website
(user.engineering.uiowa.edu/~me_160/Lab/Intermediate-2013/Lab%201-Pipe/Intermediate%20Lab%201%20Student.wbpz).
Please refer to exercise at the end of the manual to determine data and figures that need to be
saved.
Table 3 - Simulation matrix
Study
V&V of friction factor and axial velocity profile
V&V of friction factor
V&V of friction factor
V&V of friction factor
Axial velocity, centerline velocity
Axial velocity, centerline pressure, centerline velocity
Grid
2,3,4
6,7,8
0,2,4
4,6,8
8
T
Model
Laminar
Turbulent
All analytical data (AFD) for Laminar Pipe Flow and EFD data for turbulent pipe flow can be
downloaded from the class website (http://css.engineering.uiowa.edu/~me_160).
3.2. Toolbox > Component Systems. Drag and drop Geometry, Mesh and Fluent
components to Project Schematic as per below.
3.3. Right click on the upper corner of the components on the drop down arrow then select
rename. Change the names as per below.
3.4. Create connections between component as per below. You can select components part
and drop it onto the target component part to create connections.
3.5. File > Save As. Save the workbench file to H drive. The H drive is shared between the
computers in engineering labs.
4. Geometry Creation
4.1. Right click Geometry and select New Geometry. (Since all the geometries are linked
together, only one geometry creation is required)
4.3. Select the XYPlane under the Tree Outline and click New Sketch button.
4.5. Enable the auto constraints option to pick the exact point as below
4.6. Select Sketching > Rectangle. Create a rectangle geometry as per below.
4.7. Select Dimensions > General. Click on top edge then click anywhere. Repeat the same
thing for one of the vertical edges. You should have a similar figure as per below.
10
4.8. Click on H1 under Details View and change it to 7.62 m. Click on V2 and change it t
0.02619 m.
4.9. Concept > Surface From Sketches and select the sketch and hit Apply.
4.10.
11
4.11.
4.12.
If you see the lightning sign next to Geometry in the workbench then right click
on Geometry and click Update as shown below. If you dont see the check mark after
you update then you may have made a mistake when you created the geometry.
12
5. Mesh Generation
5.1. Right click Mesh and select Edit.
5.2. Right click on Mesh then select Insert > Mapped Face Meshing.
13
5.4. Click on the Edge Button. This will allow you to select edges of your geometry.
14
5.6. Hold Ctrl and select the top and bottom edge then click Apply. Specify details of sizing
as per below.
Laminar
Turbulent
15
5.7. Repeat step 5.5. Select the left and right edge and click Apply for uniform grid flow and
change sizing parameters as per below. Change the sizing parameters separately for nonuniform grid as per below. Make sure to select edges individually when changing sizing
parameters for non-uniform grid.
Uniform Grid 8
16
5.8. Click on Generate Mesh button and select Mesh under Outline.
Uniform Grid 8
Non-uniform Grid
17
5.9. Change the edge names by right clicking and selecting Create Named Selection. Name
left, right, bottom and top edges as inlet, outlet, axis and wall respectively. Your outline
should look same as the figure below.
18
Uniform Grid 8
Non-uniform Grid
5.10.
File > Save Project. Save the project and close the window. Update Mesh on
Workbench if necessary.
19
6. Solve
6.1. Right click Setup and select Edit.
20
6.3. Solution Setup > General > Check. You may ignore the warning messages. (Note: If
you get and error message you may have made a mistake while creating you mesh)
21
6.4. Solution Setup > General > Solver. Choose options shown below.
Axis Boundary Condition
Model
Laminar
Turbulent
Variable
u (m/s)
v (m/s)
P (Pa)
u (m/s)
v (m/s)
P (Pa)
k (m^2/s^2)
e (m^2/s^3)
Magnitude
Zero Gradient
22
6.5. Solution Setup > Models >Viscous-Laminar > Edit. Select parameters as per below.
Laminar flow
Turbulent flow
23
6.6. Solution Setup > Materials > air > Create/Edit Change the Density and Viscosity
as per below and click Change/Create. Close the dialog box when finished.
6.7. Solution Setup > Cell Zone Conditions > Zone > surface_body. Change type to fluid
and click OK. Select Material Name as air and click OK.
24
6.8. Solution Setup > Boundary Conditions > inlet > Edit Change parameters as per
below and click OK.
P (Pa)
Y
u (m/s)
34.08
N
v (m/s)
0
N
Laminar flow
S l
Model
Variable
Magnitude
Zero Gradient
Turbulent
P (Pa)
Intensity
0.01
Y
N
Length Scale
0.000294
N
25
Turbulent flow
S l
6.9. Solution Setup > Boundary Conditions > outlet > Edit. Change parameters as per
below and click OK.
Laminar
u (m/s) v (m/s)
Y
Y
P (Pa)
0
N
u (m/s)
Y
v (m/s)
Y
Laminar flow
Turbulent
P (Pa)
k (m^2/s^2)
400
1
N
Y
e (m^2/s^3)
1
Y
26
Turbulent flow
6.10.
Solution Setup > Boundary Conditions > wall > Edit Change parameters as
per below and click OK. No need to change for laminar cases.
Laminar
u (m/s) v (m/s)
0
0
N
N
P (Pa)
Y
u (m/s)
0
N
v (m/s)
0
N
P (Pa)
Y
Turbulent
k (m^2/s^2)
N
e (m^2/s^3)
N
Roughness
2.50E-05
-
27
Turbulent flow
6.11.
Solution Setup > Boundary Conditions > Operating Condition. Change
parameters as per below and click OK.
28
6.12.
Turbulent flow
29
6.13.
Turbulent flow
30
6.14.
Solution > Monitors > Residuals > Edit. Change convergence criterion to 1e-6
for all five equations as per below and click OK. (Note: for iterative error study you will
need to use 1e-5)
Laminar flow
Turbulent flow
31
6.15.
Solution > Solution Initialization. Change parameters as per below and click
Initialize.
Laminar flow
32
Turbulent flow
6.16.
Solution > Run calculation. Change number of iterations to 1000 and click
Calculate.
6.17.
File > save project. Make sure to save the project for later use.
33
7. Results
This section shows how to analyze your results. You do not need to do all of the analysis for every case.
Please refer to exercises to determine what analysis you need to for each simulation.
7.1.
Displaying Mesh
Select all the surface you want to display, lines and points you create can be displayed
here as well.
34
7.2.
35
File > Save Picture. Using option as per below save the residuals.
36
7.3.
Surface > Point. Change x and y values as per below click Create. Repeat this for other
lines shown in the table below.
Point
Name
point-1
point-2
point-3
point-4
point-5
point-6
point-7
point-8
point-9
point-10
x0
y0
7.62
7.62
7.62
7.62
7.62
7.62
7.62
7.62
7.62
7.62
0
0.005
0.010
0.015
0.020
0.021
0.022
0.023
0.024
0.025
Surface > Line/Rake. Change x and y values as per below click Create. Repeat this for
other lines shown in the table below.
37
Surface
Name
x=10d
x=20d
x=40d
x=60d
x=100d
X0
Y0
X1
Y1
0.5238
1.0476
2.0952
3.1428
5.238
0
0
0
0
0
0.5238
1.0476
2.0952
3.1428
5.238
0.02619
0.02619
0.02619
0.02619
0.02619
38
7.4.
Plotting Results
Results > XY Plot > Setup > Curves. For Curve # 0 select the Pattern as per below and
click Apply. Repeat this for all the curves 0 through 7.
Download the experimental data for the simulation from the class website.
(http://user.engineering.uiowa.edu/~me_160/documents/axialvelocityAFD-laminar-pipe.xy)
(http://user.engineering.uiowa.edu/~me_160/documents/axialvelocityEFD-turbulent-pipe.xy)
(http://user.engineering.uiowa.edu/~me_160/documents/pressure-EFD-turbulent-pipe.xy)
39
Results > XY Plot > Setup > Load File. Select axialvelocityEFD-turbulent-pipe.xy and
click OK.
40
41
7.5.
Exporting Data
Select Plots > XY Plot. Then change parameter as per below and click Write. This will
export the shear stress along the wall of the pipe. You will need this data to compute the
shear stress coefficient at the developed region.
42
7.6.
Results > Graphics and Animations > Vectors > Set Up Change the vector
parameters as per below and click Display.
43
Results > Graphics and Animations > Contours > Set Up Change the vector
parameters as per below and click Display.
44
7.7.
V&V Instructions
45
Change parameters as per below and click Write Make sure to select points 1 through.
46
Download V&V excel sheet for CFD Lab 1 from class website
(http://user.engineering.uiowa.edu/~me_160/Lab/Intermediate-2013/Lab%201-Pipe/Intermediate%20Lab%201%20V&V.xlsx)
Open file using Wordpad, copy points to input into V&V Excel file.
47
Paste value into V&V Excel file according to its y position and its grid number. Use the
Keep Text Only paste function by right clicking in the cell and selecting it from the paste
options.
Repeat this process for the remaining y location points and then the two remaining grid
solutions. All yellow cells should be filled.
V&V Instructions for Friction
Right click Solution > Select Edit
48
49
Name the file according to grid number and save to project folder.
50
Open file with a text editor such as TexPad and copy wall shear stress at the x location of
approximately 7m.
51
Make sure when pasting you select Keep Text Only and you select the proper cell
corresponding to the grid number.
Repeat this process for the remaining six grids. Each yellow cell should be filled.
52
8. Exercises
You need complete the following assignments and present results in your lab reports following
the lab report instructions.
* 8.1.-8.4. and 8.6. are for laminar flows, 8.5. is for turbulent flows
8.1. Iterative error studies: Use grid 4 and 8 with laminar flow conditions. Use two
different convergent limits 10-5 and 10-6 and fill in the following table for the values on friction
factors (grid 4 is given on workbench file which can be found on the class website). Find the
relative error between AFD friction factor (0.097747231) and friction factor computed by CFD,
which is computed by:
FactorCFD FactorAFD
100%
FactorAFD
To get the value of FactorCFD , you need export wall shear stress data. Then use the wall shear
stress at the developed region to calculate the friction factor. The equation for the friction factor
is C=8*/(r*U^2), where C is the friction factor, is wall shear stress, r is density and U is the
inlet velocity. Discuss the effect of convergent limit on results for these two meshes
Mesh No.
4
(
%)
( %)
8
(
%)
( %)
Figure need to be reported: residuals history for mesh 8 for two convergent limits.
Data need to be reported: the above table with values.
8.2.
Verification study for friction factor of laminar pipe with refinement ratio
1.41: Use the simulations with the meshes for grid 0, 2, 3, 4, 6, 7, and 8 with convergence limit
1e-6 (Except for mesh 8 other meshes and their setup is provided on the workbench file in the
class website). Export friction factor and insert the values into V&V excel sheet (Refer to V&V
instructions). For each parameter, refer to Nomenclature sheet in V&V excel sheet.
Which set of meshes is closer to the asymptotic range? Which set has a lower grid uncertainty
(Ug)? Which set is closer to the theoretical value of order of accuracy (2nd order)? For the fine
mesh 8, also compare its relative error of the friction factor (the one using convergent limit 10-6
in the table in exercise 8.1) with the grid uncertainty for 6,7,8, which is higher and what does that
mean?
53
8.3.
Verification study for friction factor of laminar pipe with refinement ratio 2:
Use the simulation for the meshes 0, 2, 4, 6 and 8 with convergence limit 1e-6. Export friction
factor and insert the values into V&V excel sheet (Refer to V&V instructions).
Compared to results in 8.2, which set of meshes is sensitive to grid refinement ratio? Why?
8.4.
Verification study of axial velocity profile: Use mesh 4 as the fine mesh, use
grid refinement ratio 1.414 and convergence limit 10-6. Follow the V&V in the results section.
Save the figures and discuss if the simulation has been verified. Discuss which mesh solution is
closest to the AFD data, why?
Figures need to be reported: Figures and tables in the V&V excel sheet.
Figures need to be reported: Axial velocity profile with EFD data, normalized axial
velocity profile at x=100D, centerline pressure distribution with EFD data, centerline velocity
distribution, contour of axial velocity, velocity vectors showing the developing region and
developed regions.
Data need to be reported: Developing length and compared it with that using formula in
textbook.
Comparison between laminar and turbulent pipe flow
8.6.
Compare the results of laminar pipe flow using mesh 8 in exercise 8.1 (convergent limit 10-6)
with results of turbulent pipe flow in exercise 8.5. Analyze the difference in normalized axial
velocity profile and developing length for laminar and turbulent pipe flows.
Figures need to be reported: Axial velocity profile with AFD data, normalized axial
velocity profile at x=100D, centerline velocity distribution for laminar flows.
Data need to be reported: Developing length for laminar pipe flow and compared it with
that using formula in textbook.
8.7.
Questions need to be answered in CFD Lab report
8.7.1. Answer all the questions in exercises 8.1 to 8.6
8.7.2. Analyze the difference between CFD/AFD and CFD/EFD and possible error sources.