Tutorials Report Computational Fluid Dynamics - CH5101: Tutorial 1: - Problem Statement - Fluid Flow in A 2D Pipe
Tutorials Report Computational Fluid Dynamics - CH5101: Tutorial 1: - Problem Statement - Fluid Flow in A 2D Pipe
Assumptions –
Laminar flow
Constant Viscosity
Neglect backflow
Incompressible fluid
Neglecting energy equations
4m
Inlet Velocity 0.005 m/s
40 m
The mesh created in this case is refined for finer mesh elements, thus acquiring more accurate
data.
Figure 1 – 2D Pipe Mesh
For initialisation, we use standard initialisation for this simulation. However, hybrid
initialisation can also be used if the convergence is too sensitive to initial values.
With a go-to guess for the start of the inlet velocity, we set the initial x and y component of the
velocities to zero (0 m/s). Similarly, the initial gauge pressure is also set to zero (0 Pa).
For the boundary conditions, apart from the zero velocity at the wall (no slip condition at the
walls), the other parameters are defined as follows –
1.) Inlet velocity at the left entrance = 0.0005 m/sec
2.) Exit gauge pressure at the exit on the right boundary = 0 Pa
3.) Zero velocity or no slip condition at the walls.
The solver for this simulation relies on pressure based solving with the aim to converge solution
at steady state condition. Also, the following conditions and parameters were set for this
simulation.
1) Study type – Steady state
2) Viscous model – Laminar
3) Domain material – Water (liquid)
4) Inlet velocity – 0.0005 m/sec
5) Pressure – velocity coupling scheme: SIMPLE
6) spatial discretisation
a. Gradient – least square cell based
b. Pressure – Second order
c. Momentum – Second order upwind
Results and Analysis –
The simulation is run with a maximum limit of 1000 iterations with reporting and update
interval on 1. The solution was converged in just 107 iterations, with the standard error
tolerance fulfilling the criterion.
Figure 2
Figure 3
Remarks –
The results depict the expected velocity profile of a laminar flow of a fluid in a 2D pipe. The
parabolic profile is clearly visible in figure 3. The results are in agreement with the theoretical
profile defined as follows.