Gold:UFR3-05 instruct: Difference between revisions
No edit summary |
No edit summary |
||
Line 5: | Line 5: | ||
<div class="Section1"> | |||
Instructions for transonic axisymmetric bump flow calculation | |||
<nowiki>=============================================================</nowiki> | |||
Grid | |||
<nowiki>====</nowiki> | |||
A Fortran program for generating a single-block two-dimensional grid, | |||
together with sufficient documentation, can be found in files | |||
[../C/gridaxibump.f.txt gridaxibump.f] and [../C/gridaxibump.htm gridaxibump] | |||
For an axisymmetric calculation the 2D plane should be rotated as | |||
mentioned at the top of gridaxibump.f | |||
Boundary Conditions | |||
<nowiki>===================</nowiki> | |||
Boundary conditions for the variables are as follows: | |||
<span style="mso-tab-count: 1"> </span>X=Xmin: inflow - Uniform inlet Mach Number of 0.875 for axial | |||
<span style="mso-tab-count: 2"> </span>component and zero value for others | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>X=Xmax: outflow - zero longitudinal gradient | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>Y=Ymin: no-slip wall | |||
<span style="mso-tab-count: 1"> </span>Y=YMAX: Euler wall<span style="mso-tab-count: 1"> </span> | |||
Notes about the dimension of the computational domain: | |||
<span style="mso-tab-count: 1"> </span>1. XMAX set at x/c=3.5 from bump trailing edge. This is sufficiently far <span style="mso-tab-count: 2"> </span> | |||
<span style="mso-tab-count: 2"> </span> from the zone of interest; here c is the bump chord length. | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>2. YMAX set at 4.5*c ensures that there is no shock reflection. | |||
<span style="mso-tab-count: 2"> </span>There is, however, a fluctuation of ~1% of free-stream Mach No. | |||
<span style="mso-tab-count: 2"> </span>on the top boundary. This is found to have negligible effect | |||
<span style="mso-tab-count: 2"> </span>on the critical flow features such as CP. | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>3. XMIN set at 4.0*c upstream from the bump leading edge. | |||
<span style="mso-tab-count: 2"> </span>After several trials, we found that if we specify | |||
<span style="mso-tab-count: 2"> </span>a plug velocity profile at this location, the corresponding | |||
<span style="mso-tab-count: 2"> </span>profile at x/c=-0.25 matches with experiment reasonably well. | |||
<span style="mso-tab-count: 2"> </span>However, other inlet profiles with different XMIN location | |||
<span style="mso-tab-count: 2"> </span>may be possible.<span style="mso-spacerun: yes"> </span> | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>4. All of above observations are based on high-Re k-e calculations. | |||
<span style="mso-tab-count: 2"> </span> | |||
Experimental Data | |||
<nowiki>=================</nowiki> | |||
The experimental data at different axial locations are given in files | |||
Experiment-CP.txt and Experiment-UV.txt | |||
<span style="mso-tab-count: 1"> </span>Wall static Pressure (CP) is calculated as | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>CP = (p-p0)/0.5*rho0*u0**2 | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>p0, rho0 and u0 are the free-stream quantities | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>X: Normalised distance along the flow. (=x/c, where c is the | |||
<span style="mso-tab-count: 2"> </span>bump chord length. X=1.0 corresponds to the bump trailing edge. | |||
<span style="mso-tab-count: 1"> </span>Y: Vertical distance from the bottom solid wall (=y/c) | |||
<span style="mso-tab-count: 1"> </span>U: Normalised streamwise velocity (=u/u0) | |||
<span style="mso-tab-count: 1"> </span>V: Normalised transverse velocity (=v/u0) | |||
<span style="mso-tab-count: 1"> </span>UU: Normalised streamwise component of normal stress (= u'u'/u0**2) | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>VV: Normalised transverse component of normal stress (= v'v’/u0**2) | |||
<span style="mso-tab-count: 1"> </span>MUV: Normalised Reynolds shear stress (= -u'v'/uo**2) | |||
<span style="mso-tab-count: 1"> </span>TKE: Turbulent Kinetic Energy (=k/u0**2) | |||
<span style="mso-tab-count: 1"> </span>where, k = 0.5*(u'**2+v'**2+w'**2) | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>Since only u' and v' were measured, the third component was | |||
<span style="mso-tab-count: 1"> </span>calculated from : | |||
<span style="mso-tab-count: 1"> </span> | |||
<span style="mso-tab-count: 1"> </span>w'**2=0.5*(u'**2+v'**2) | |||
NOTE: Please note that at some locations, data for all the above | |||
variables were not always available. This may be recognized in the | |||
data sets below by the appearance of a '999' which does not represent | |||
a real value. | |||
CFD Calculations | |||
<nowiki>=================</nowiki> | |||
The data derived from CFD calculations using a number of different turbulence models can be found in: | |||
[../U3-05des.htm#CFD_Data CFD Files] | |||
The interpretation of the tabulated data is the same as that above for the experimental data with the following additions. | |||
CF =<span style="mso-spacerun: yes"> </span>(wall shear stress)/(0.5*rho0*u0**2) | |||
NUT = <span style="mso-spacerun: yes"> </span>(Turbulent Viscosity)/(rho*u0*c) | |||
The final column of data is the normalised second scale determining variable (e.g e or w etc.) | |||
</div> | |||
{{UFR|front=UFR 3-05|description=UFR 3-05 Description|references=UFR 3-05 References|testcase=UFR 3-05 Test Case|evaluation=UFR 3-05 Evaluation|qualityreview=UFR 3-05 Quality Review|bestpractice=UFR 3-05 Best Practice Advice|relatedACs=UFR 3-05 Related ACs}} | {{UFR|front=UFR 3-05|description=UFR 3-05 Description|references=UFR 3-05 References|testcase=UFR 3-05 Test Case|evaluation=UFR 3-05 Evaluation|qualityreview=UFR 3-05 Quality Review|bestpractice=UFR 3-05 Best Practice Advice|relatedACs=UFR 3-05 Related ACs}} |
Revision as of 13:49, 11 April 2010
Instructions for transonic axisymmetric bump flow calculation
=============================================================
Grid
====
A Fortran program for generating a single-block two-dimensional grid,
together with sufficient documentation, can be found in files
[../C/gridaxibump.f.txt gridaxibump.f] and [../C/gridaxibump.htm gridaxibump]
For an axisymmetric calculation the 2D plane should be rotated as
mentioned at the top of gridaxibump.f
Boundary Conditions
===================
Boundary conditions for the variables are as follows:
X=Xmin: inflow - Uniform inlet Mach Number of 0.875 for axial
component and zero value for others
X=Xmax: outflow - zero longitudinal gradient
Y=Ymin: no-slip wall
Y=YMAX: Euler wall
Notes about the dimension of the computational domain:
1. XMAX set at x/c=3.5 from bump trailing edge. This is sufficiently far
from the zone of interest; here c is the bump chord length.
2. YMAX set at 4.5*c ensures that there is no shock reflection.
There is, however, a fluctuation of ~1% of free-stream Mach No.
on the top boundary. This is found to have negligible effect
on the critical flow features such as CP.
3. XMIN set at 4.0*c upstream from the bump leading edge.
After several trials, we found that if we specify
a plug velocity profile at this location, the corresponding
profile at x/c=-0.25 matches with experiment reasonably well.
However, other inlet profiles with different XMIN location
may be possible.
4. All of above observations are based on high-Re k-e calculations.
Experimental Data
=================
The experimental data at different axial locations are given in files
Experiment-CP.txt and Experiment-UV.txt
Wall static Pressure (CP) is calculated as
CP = (p-p0)/0.5*rho0*u0**2
p0, rho0 and u0 are the free-stream quantities
X: Normalised distance along the flow. (=x/c, where c is the
bump chord length. X=1.0 corresponds to the bump trailing edge.
Y: Vertical distance from the bottom solid wall (=y/c)
U: Normalised streamwise velocity (=u/u0)
V: Normalised transverse velocity (=v/u0)
UU: Normalised streamwise component of normal stress (= u'u'/u0**2)
VV: Normalised transverse component of normal stress (= v'v’/u0**2)
MUV: Normalised Reynolds shear stress (= -u'v'/uo**2)
TKE: Turbulent Kinetic Energy (=k/u0**2)
where, k = 0.5*(u'**2+v'**2+w'**2)
Since only u' and v' were measured, the third component was
calculated from :
w'**2=0.5*(u'**2+v'**2)
NOTE: Please note that at some locations, data for all the above
variables were not always available. This may be recognized in the
data sets below by the appearance of a '999' which does not represent
a real value.
CFD Calculations
=================
The data derived from CFD calculations using a number of different turbulence models can be found in:
[../U3-05des.htm#CFD_Data CFD Files]
The interpretation of the tabulated data is the same as that above for the experimental data with the following additions.
CF = (wall shear stress)/(0.5*rho0*u0**2)
NUT = (Turbulent Viscosity)/(rho*u0*c)
The final column of data is the normalised second scale determining variable (e.g e or w etc.)