Abaqus/CFD - Sample Problems
Abaqus/CFD - Sample Problems
Abaqus/CFD - Sample Problems
Abaqus 6.10
Contents
1.
2.
3.
4.
5.
6.
1.
Problem description
A rectangular 2-dimensional channel of width = 1m and length = 2m is considered. An
oscillatory pressure gradient (with zero mean) is imposed at the inlet. The analysis is carried
out in two steps. In the first analysis step, a constant pressure gradient is prescribed for the
first 5 seconds of the simulation to initialize the velocity field to match that of the analytical
steady-state solution. In the second analysis step, the flow is subjected to an oscillatory
pressure gradient. A 40x20 uniform mesh is used for this problem. Two dimensional
geometry is modeled as three dimensional with one element in thickness direction.
Wall boundary condition
1m
z
x
2m
Outlet
dP
P0 cos( (t to ))
dx
Boundary conditions
Pressure inlet
t < to : p = 7.024
t > to : p = 10*Cos((t-to)) ; to = 5, = p/5
Pressure outlet (p = 0)
No-slip wall boundary condition on top and bottom (V = 0)
Analytical solution
dP
P0 e it
dx
Po 2s it cos(z h)
u ( y, t ) Re
e 1
cos(h)
2i
References
Fluid Mechanics, Second Edition: Volume 6 (Course of Theoretical Physics), Authors: L. D. Landau,
E.M. Lifshitz
T= 5 sec
T= 10 sec
Files
ex1_oscillatory_planeflow.py
ex1_oscillatory_planeflow_mesh.inp
Problem description
The following cavity configurations and Reynolds numbers are presented
1m
1m
1m
1m
1m
1m
1m
1m
1m
Boundary conditions
Specified velocity at top plane to match flow Reynolds number
No-slip at all other planes (V = 0)
Hydrostatic mode is eliminated by setting reference pressure to zero at a single node
References
1.
High-Re Solutions for Incompressible Flow Using the Navier-Stokes Equations and a Multigrid Method
U. Ghia, K. N. Ghia, and C. T. Shin
Journal of Computational Physics 48, 387-411 (1982)
2.
3.
Results
Velocities along horizontal and vertical centerlines
Re = 100
Re = 3200
10
Results
Velocities along horizontal and vertical centerlines
450
150
11
Results
Velocity vectors
12
13
Problem description
The material properties are chosen to match the desired Rayleigh number, Ra
Ra
q = 0.0
g L3 T
, Kinematic viscosity
, Thermal diffusivity
, Thermal expansion coefficient
g, Acceleration due to gravity
T = 1.0
T = 0.0
x = 0, T = 1.0
q = 0.0
Square Cavity: Rayleigh number = 1e3, 1e6
Cubical Cavity: Rayleigh number = 1e4
x = 1.0, T = 1.0
Boundary conditions
No-slip velocity boundary condition on all the planes (V = 0)
Specified temperatures
Hydrostatic mode is eliminated by setting reference pressure to zero at a single node
References
1.
2.
"Benchmark solutions for natural convection in a cubic cavity using the high-order timespace method"
Shinichiro Wakashima, Takeo S. Saitoh
International Journal of Heat and Mass Transfer, 47, 853864, (2004)
16
Results
Temperature
Velocity
Present
Benchmark
(de Vahl Davis)
3.695
3.697
0.1781
0.178
Present
Benchmark
(de Vahl Davis)
3.654
3.649
0.8129
0.813
Present
Benchmark
(de Vahl Davis)
q = 0.0
T = 1.0
T = 0.0
q = 0.0
Ra = 1000, Pr = 0.71
Temperature
Velocity
Present
Benchmark
(de Vahl Davis)
219.747
219.36
65.9
64.63
0.0375
0.0379
0.85
0.85
Ra = 1.0e6 , Pr = 0.71
17
Results
Ra = 1.0e4, Pr = 0.71
max
Benchmark
(Wakashima &
Saitoh (2004))
1.1018
0.1984
(z = 0.8250)
0.2216
(x = 0.1177)
Abaqus/CFD
1.1017
0.1986
0.2211
Error
-0.009%
0.1%
0.2%
Present
Benchmark
Present
Benchmark
18
19
4. Turbulent Flow in a
Rectangular Channel
20
Problem description
A rectangular 2-dimensional channel of width = 1 unit and length = 10 units is considered. A
pressure gradient is imposed along the length of the channel by means of specified pressure
at the inlet and zero pressure at the outlet. The pressure gradient is chosen to impose the
desired friction Reynolds number for the flow.
Wall boundary condition
Inlet
Outlet
H=1
Wall boundary condition
L = 10
ut
tw
r
Features
Turbulent flow
Spalart-Allmaras turbulence model
Boundary conditions
No-slip velocity boundary condition on channel walls(V = 0)
Set through thickness velocity components to zero
Distance function = 0, Kinematic turbulent viscosity = 0 at No-slip velocity boundary condition
Mesh
Mesh 50 (Streamwise) X 91 (Normal) X 1 (Through thickness)
y+ at first grid point ~ 0.046
22
2.
23
U ut y
,y
ut
24
Velocity profile
25
Note
Abaqus/CFD results at a friction Reynolds number ~ 395 and Reynolds number
(based on mean velocity) ~ 13750 can be obtained by setting a pressure
gradient of 0.0573
26
Problem description
The computational domain for the vortex shedding calculations consists of an interior square
region (-4 < x < 4; -4 < y < 4) surrounding a cylinder of unit diameter. The domain is
extended in the wake of cylinder up to x = 20 units
Inlet
D=1
Outlet
H=8
Tow-tank condition
H = 20
Fluid Properties
Density = 1 unit
Viscosity = 0.01 units
Boundary conditions
Tow tank condition at top and bottom walls:
U = 1, V = 0
29
30
Strouhal Number*
Coarse
1760
0.1749
0.03
Fine
28160
0.1735
0.0075
The results are consistent with the Strouhal numbers reported in the reference
(0.172-0.173)
The frequency of the vortex shedding is obtained as the half of the peak frequency
obtained by a FFT of kinetic energy plot (from t = 200 sec to 300 sec)
31
Kinetic energy
32
Problem description
The computational domain for the flow calculations consists of an rectangular region (0 < x <
30; -0.5 < y < 0.5). The flow enters the solution domain from 0< y < 0.5 while -0.5 < y < 0.0
represents the step. A parabolic velocity profile is specified at the inlet.
No-slip
Vx(y) = 24y(0.5 y)
Vav = 1.0
Outlet
P=0
H=1
No-slip
No-slip
L = 30
Re
rVav H
m
Fluid Properties
Density = 1 unit
Viscosity = 0.00125 units (the viscosity is chosen so as to set the flow Reynolds
number to 800)
Boundary conditions
Set through thickness velocity components to zero
No-slip velocity boundary condition at top and bottom walls
Vx = 0, Vy = 0
Outlet: P = 0
No-slip velocity boundary condition at the step boundary
Vx = 0, Vy = 0
36
37
L3
L2
L1
Zone 1 (0 < x < 15)
Uniform mesh
x, y
Mesh
(Across the channel x along
the channel length)
L1
Length from the step face to the
lower re-attachment point
L2
Length from the step
face to upper
separation point
L3
Length of the upper
separation bubble
Gartling (1990)
40x800
6.1
4.85
10.48
Abaqus/CFD
Fine
80x1200x1 (Zone 1)
80x832x1 (Zone 2)
5.9919
4.9113
10.334
Abaqus/CFD
Medium
40x600x1 (Zone 1)
40x416x1 (Zone 2)
5.7471
4.8379
10.101
Abaqus/CFD
Coarse
20x300x1 (Zone 1)
20x208x1 (Zone 2)
4.5018
3.9659
8.7748
Elements used by Gartling (1990) were biquadratic in velocity and linear discontinuous pressure elements. In
contrast, the fluid elements in Abaqus/CFD use linear discontinuous in velocity and linear continuous in
pressure.
38
40
41
coarse_parabolic_inlet_velocity.inp
medium_parabolic_inlet_velocity.inp
fine_parabolic_inlet_velocity.inp
Note
The models require a parabolic velocity profile at the inlet. This needs to be
manually included as boundary condition in the generated input file.
The parabolic velocity profile required is provided in files
coarse_parabolic_inlet_velocity.inp, medium_parabolic_inlet_velocity.inp and
fine_parabolic_inlet_velocity.inp for coarse, medium and fine meshes,
respectively.
42