α = 1.
55o
M∞ = 0.7
Dimension – 2D, double precision, display all on, serial. Working directory G:\Lab4.
Step -1: Mesh
Read the mesh file
File read mesh [Link]
Check the mesh
General scale (-20m, 25m, -30 m 30 m)
General check (no negative volume)
If reorder is present use it (mesh reorder domain)
Step – 2: General Setting
Solver setting: density based, absolute, steady and planar
Step – 3: Models
Models Energy edit (on)
Models Viscous k-omega (2-eqn) and SST
(1, 0.52,.09,1.5)
Step -4: Materials
Maerials air create/edit
Ideal gas, cp=1006.4, thermal conductivity 0.0242, Sutherlandlaw 1.715e-5,273.11, 110.56
Change/create to save these settings and close the materials dialog box
Step -5: operating condition
cell zone conditions, operating pressure = 0 Pa
Step – 6: Boundary condition
Wall bottom – boundary zone type wall, Edit: Momentum stationary wall, no slip, roughness =0 h,
roughness constant =0.5.
Adiabatic wall, heat flux set to zero.
Wall top check noslip, adiabatic
pressure-far-field boundary is applicable only when using the ideal-gas law.
For example, in lifting airfoil calculations, it is not uncommon for the far-field
boundary to be a circle with a radius of 20 chord lengths.
On the ‘Momentum’ tab
set the gauge static pressure to 73048 Pa
set the Mach Number to 0.7
The angle of attack (α) in this case is 1.55 deg. The x-component of the flow
is cos α and the y-component is sin α.
X – component 0.99963
y- component 0.02705
Select ‘Intensity and Viscosity Ratio’
• Set the intensity to 1%
• Set the viscosity ratio to 1
Select the thermal tab
Temperature 283.24 K
Set the reference values: These are not used to compute the flow solution, but they are used to
report coefficients, such as Cp. Use the free-stream as a reference condition
‘compute from farfield’
Set the following to represent a chord length of 1m
with unit depth:
reference length = 1m
reference depth = 1m
reference area = 1m2
These are not updated from the farfield conditions, but the defaults happen to be 1 already.
Monitor
Monitors > Residuals > Edit
– Make sure ‘Plot’ is on
Set up a monitor for the drag coefficient on the airfoil
Select both wall zones and toggle on ‘Print’, ‘Plot’ and ‘Write’.The α is 1.55° so we need to use the force
vector as shown.
• Lift and drag are defined (perpendicular and parallel respectively) relative to the freestream
flow direction, not the airfoil.
0.99963, 0.02705
Set up a monitor for the lift coefficient on the airfoil in the same way
-0.02705, 0.99963
Save – case and data
Run
Graphics and Animations > Contours
– Turn off ‘Filled’ to just display the
contour lines.
– Turn on ‘Filled’, display again
====================================================================================
