Boundary Layer Calculator

Download as xls, pdf, or txt
Download as xls, pdf, or txt
You are on page 1of 10

EVBSPBLC --- EXCEL VISUAL BASIC STAGNATION POINT

BOUNDARY LAYER CALCULATOR

Hypersonic Cold Wall (Fig. 8-3, 12,000 ft/sec @ 70,000 ft, Hayes & Probstein Hypersonic Flow Theory, 1959)

Solution: Converged
Iterations:
22

Flow Properties:
Minf
12.3597
Pinf (lbf/ft)
93.7255
Tinf (R)
392.246
Twall (R)
540
Wall Type
Temp
RhoVbar
0

Axisymmetric Spherical Stagnation Point Profiles:


Eta
Y (ft)
u/Uedge H/Hedge T/Tedge
0
0
0 0.043637 0.053034
0.1 3.65E-07 0.023521 0.059723 0.072501
0.2 8.46E-07
0.04887 0.077144 0.093292
0.3 1.45E-06 0.076011 0.095916 0.115153
0.4 2.19E-06 0.104854
0.11602 0.137886
0.5 3.06E-06 0.135279
0.13742 0.161399
0.6 4.07E-06 0.167153 0.160078
0.18567
0.7 5.22E-06 0.200336
0.18395 0.210705
0.8 6.52E-06 0.234683 0.208992 0.236513
0.9 7.98E-06 0.270041 0.235153 0.263097
1 9.59E-06 0.306249 0.262378 0.290448
1.1 1.14E-05 0.343136 0.290603 0.318541
1.2 1.33E-05 0.380522 0.319753 0.347336
1.3 1.54E-05 0.418215 0.349741 0.376776
1.4 1.77E-05 0.456014 0.380467 0.406788
1.5 2.01E-05 0.493711 0.411817
0.43728
1.6 2.28E-05
0.53109 0.443663 0.468145
1.7 2.56E-05 0.567932
0.47586 0.499259
1.8 2.86E-05 0.604017
0.50825 0.530483
1.9 3.18E-05 0.639128 0.540664 0.561664
2 3.51E-05 0.673057
0.57292
0.59264
2.1 3.87E-05 0.705609 0.604831 0.623237
2.2 4.24E-05 0.736605 0.636205 0.653282
2.3 4.63E-05
0.76589
0.66685 0.682596
2.4 5.03E-05 0.793334 0.696579 0.711008
2.5 5.45E-05 0.818837 0.725218 0.738356
2.6 5.89E-05 0.842333 0.752603 0.764489
2.7 6.34E-05 0.863788 0.778593 0.789275
2.8 6.81E-05 0.883202 0.803067 0.812604
2.9 7.29E-05 0.900609 0.825933
0.83439
3 7.78E-05
0.91607 0.847124 0.854573
3.1 8.28E-05 0.929675 0.866603
0.87312
3.2
8.8E-05 0.941533 0.884362 0.890024
3.3 9.32E-05
0.95177 0.900419 0.905304

EXCEL Worksheet STAGPTCBL

Gas Properties:
Gamma
1.4
Gas Con
1716
Thetav
5500
Pr
0.71
Suther
2.27E-08
Tconstant
198.6

Gamma
1.399375
1.394898
1.38415
1.369585
1.354888
1.34205
1.331547
1.32318
1.316567
1.311333
1.30717
1.303835
1.301144
1.298955
1.297162
1.295684
1.294458
1.293435
1.292577
1.291855
1.291243
1.290724
1.290282
1.289905
1.289581
1.289304
1.289066
1.288861
1.288685
1.288534
1.288404
1.288292
1.288196
1.288114

Page 1

Solution Properties:
Etainf
6
Delta Eta
0.005
Conv Test 1.00E-10
Max Iter
30
Print Eta
0.1

Stagnation Conditions:
P0' (psf)
T0' (R) H0' (Btu/lbm)
18802.76 10182.22 2968.859
Qdot (Btu/ft-sec)
EVBSPBLC
Fay-Riddell (Sphere)
-2595.19
-2617.75
Profile Integral Parameters:
Delta* (ft) Theta (ft) Delta*/Theta
-8.9E-06 3.14E-05 -0.28452
Reinf(Rnose)
Stinf(H0'-Hwall)
116613.9
0.017004
Stinf * Sqrt( Reinf / C* )
0.676947

Velocity Profile
7
6
5
Eta

Body Properties:
Geometry
Sphere
Rnose (ft) 0.020833

4
3
2
1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1
u / Uedge

10/5/2013 @ 8:30 PM

Temperature Profile
7
6
5
4

Eta

Gamma
1.288043
1.287983
1.287932
1.287889
1.287852
1.287822
1.287796
1.287775
1.287757
1.287742
1.28773
1.287721
1.287713
1.287706
1.287701
1.287697
1.287694
1.287692
1.28769
1.287688
1.287687
1.287686
1.287685
1.287685
1.287684
1.287684
1.287684

3
2
1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1

T / Tedge

Total Enthalpy Profile


7
6
5

Eta

Axisymmetric Spherical Stagnation Point Profiles:


Eta
Y (ft)
u/Uedge H/Hedge T/Tedge
3.4 9.85E-05 0.960526 0.914818 0.919003
3.5 0.000104 0.967944 0.927623 0.931184
3.6 0.000109 0.974171 0.938917 0.941927
3.7 0.000115 0.979349 0.948798 0.951323
3.8
0.00012 0.983616 0.957372 0.959476
3.9 0.000126 0.987101 0.964752 0.966493
4 0.000132 0.989922 0.971054 0.972485
4.1 0.000137 0.992185 0.976394 0.977561
4.2 0.000143 0.993987 0.980882 0.981828
4.3 0.000149 0.995408 0.984626 0.985387
4.4 0.000154
0.99652 0.987725 0.988333
4.5
0.00016 0.997383 0.990271 0.990753
4.6 0.000166 0.998047 0.992347 0.992726
4.7 0.000172 0.998555 0.994027 0.994323
4.8 0.000178 0.998939 0.995377 0.995606
4.9 0.000183 0.999229 0.996453 0.996629
5 0.000189 0.999445 0.997306 0.997439
5.1 0.000195 0.999604 0.997975 0.998076
5.2 0.000201 0.999722 0.998498 0.998572
5.3 0.000207 0.999808 0.998903 0.998957
5.4 0.000212 0.999869 0.999214 0.999253
5.5 0.000218 0.999914 0.999452 0.999479
5.6 0.000224 0.999945 0.999632
0.99965
5.7
0.00023 0.999967 0.999768 0.999779
5.8 0.000236 0.999983 0.999869 0.999876
5.9 0.000242 0.999993 0.999945 0.999947
6 0.000247
1
1
1

4
3
2
1
0
0

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

1.1

H / Hedge

Ratio of Specific Heats


7
6
5
Eta

4
3
2
1
0
1.25

EXCEL Worksheet STAGPTCBL

Page 2

1.3

1.35
Gamma

1.4

10/5/2013 @ 8:30 PM

1.45

NOMENCLATURE

Definition
Body Properties:
Geometry

Rnose (ft)

Flow Properties:
Minf
Pinf (lbf/ft)
Tinf (R)
Twall (R)
Wall Type

RhoVbar

Gas Properties:
Gamma
Gas Con
Thetav

Geometry Type
For Two-Dimensional Cylindrical Enter: Cylinder
For Axisymmetric Spherical Enter: Sphere
Cylinder or Sphere Nose Radius

Free-Stream Mach Number


Free-Stream Static Pressure
Free-Stream Static Temperature
Wall Temperature
For Adiabatic Wall: Twall Will Be Automatically Calculated
Prescribed Wall Temperature or Adiabatic Wall
For Prescribed Wall Temperature Enter: Temp
For Adiabatic Wall Enter: Adi
Wall Mass Transfer Parameter: RhowallVwall / RhoinfVinf
For Wall Mass Addition/Blowing/Injection Enter: (+) RhoVbar
For Wall Mass Removal/Suction Enter: (-) RhoVbar

Pr
Suther
Tconstant

Ratio of Specific Heats (Cp / Cv)


Specific Gas Constant
Characteristic Vibrational Temperature
For Perfect Gas Solution, Enter Thetav = 0.0
Prandtl Number
Sutherland Viscosity Law Constant
Sutherland Viscosity Law Reference Temperature

Solution Properties:
Etainf
Delta Eta
Conv Test
Max Iter
Print Eta

Maximum Eta Value ( 0 <= Eta <= Etainf)


Eta Grid Increment in Solution
Solution Convergence Test
Maximum Number of Allowed Solution Iterations
Eta Grid Increment in Profile Printout

Stagnation Conditons:
P0'
T0'
H0'

Stagnation Point Total Pressure


Stagnation Point Total Temperature
Stagnation Point Total Enthalpy

Qdot:
EVBSPBLC
Fay-Riddell

Stagnation Point Heat Transfer


Computed From EVBSPBLC Boundary-Layer Profile Solution
Computed From Fay-Riddell Correlation for a Sphere or Cylinder

Profile Integral Parameters:


Delta*
Theta
Delta*/Theta

Reinf(Rnose)
Stinf(H0'-Hwall)
Correlation Parameter

Boundary-Layer Displacement Thickness


Boundary-Layer Momentum Thickness
Boundary-Layer Shape Factor (Delta* / Theta)

Free-Stream Reynolds Number Based on Nose Radius


Reinf = Rhoinf * Vinf * Rnose / Muinf
Free-Stream Stanton Number Based on ( H0' - Hwall )
Stinf = -Qdot / [Rhoinf * Vinf * ( H0' - Hwall )]
Stinf * Sqrt( Reinf / C* )
C* = ((Rhoe * Mue) ^ 0.8) * ((Rhowall * Muwall) ^ 0.2) / (Rhoinf * Muinf)
Rhoe = Rho( P0', T0' )
Mue = Mu( T0' )
Rhowall = Rho( P0', Twall )
Muwall = Mu( Twall )
Rhoinf = Rho( Pinf, Tinf )
Muinf = Mu( Tinf )

Units

Applicable Values:
Air
Nitrogen

Not Applicable

ft

Nondimensional
lbf / ft or psf
R
R
Not Applicable

Nondimensional

Nondimensional
ft / sec - R
R
Nondimensional
lbf - sec / ft - R^0.5
R

Nondimensional
Nondimensional
Not Applicable
Not Applicable
Nondimensional

lbf / ft or psf
R
Btu / lbm

Btu / ft - sec
Btu / ft - sec
Btu / ft - sec

1.4
1716
5500

1.4
1774
6102

0.71
2.27E-08
198.6

0.71
2.18E-08
192

ft
ft
Nondimensional

Nondimensional
Nondimensional
Nondimensional

EXAMPLE INPUT

Body Properties:

Flow Properties:

Gas Properties:

Solution Properties:

Hypersonic Cold Wall (Fig. 8-3, 12,000 ft/sec @ 70,000 ft, Hayes & Probstein Hypersonic Flow Theory, 1959)
Geometry Sphere
Minf
12.3597
Gamma
1.4
Rnose (ft) 0.020833
Pinf (lbf/ft) 93.7255
Gas Con
1716
Tinf (R)
392.246
Thetav
5500
Twall (R)
540
Pr
0.71
Wall Type
Temp
Suther
2.27E-08
RhoVbar
0
Tconstant
198.6

Cylinder With Prescribed Hot Wall Under Subsonic Flow Conditions (Note Overshoot in Velocity Profile)
Geometry Cylinder
Minf
0.25
Gamma
1.4
Rnose (ft)
0.5
Pinf (lbf/ft)
2200
Gas Con
1716
Tinf (R)
550
Thetav
0
Twall (R)
900
Pr
0.71
Wall Type
Temp
Suther
2.27E-08
RhoVbar
0
Tconstant
198.6

Cylinder With Adiabatic Wall (DO NOT INPUT Twall for Adiabatic Wall Cases)
Geometry Cylinder
Minf
3
Gamma
1.4
Rnose (ft)
0.5
Pinf (lbf/ft)
1000
Gas Con
1716
Tinf (R)
400
Thetav
0
Twall (R)
Pr
0.71
Wall Type
Adi
Suther
2.27E-08
RhoVbar
0
Tconstant
198.6

Sphere With Massive Wall Injection @ Prescribed Wall Temp (Note Boundary Layer Effectively Blown Off)
Geometry Sphere
Minf
3
Gamma
1.4
Rnose (ft)
0.1
Pinf (lbf/ft)
1000
Gas Con
1716
Tinf (R)
400
Thetav
0
Twall (R)
540
Pr
0.71
Wall Type
Temp
Suther
2.27E-08
RhoVbar
0.01
Tconstant
198.6

Cylinder With Wall Suction Under Subsonic Flow Conditions and Adiabatic Wall
Geometry Cylinder
Minf
0.25
Gamma
1.4
Rnose (ft)
0.5
Pinf (lbf/ft)
2000
Gas Con
1716
Tinf (R)
540
Thetav
0
Twall (R)
Pr
0.71
Wall Type
Adi
Suther
2.27E-08
RhoVbar
-0.001
Tconstant
198.6

Reference Sphere (1 ft Radius) in Nitrogen @ Minf = 13.5 & Reinf/ft = 1.37E+06


Geometry Sphere
Minf
13.5
Gamma

1.4

Rnose (ft)

Pinf (lbf/ft)
Tinf (R)
Twall (R)
Wall Type
RhoVbar

2.4768
91.2
540
Temp
0

Gas Con
Thetav
Pr
Suther
Tconstant

1774
6102
0.71
2.18E-08
192

Solution Properties:

onic Flow Theory, 1959)


Etainf
Delta Eta
Conv Test
Max Iter
Print Eta

6
0.005
1.00E-10
30
0.1

t in Velocity Profile)
Etainf
Delta Eta
Conv Test
Max Iter
Print Eta

5
0.005
1.00E-10
30
0.1

Etainf
Delta Eta
Conv Test
Max Iter
Print Eta

4
0.005
1.00E-10
30
0.1

Effectively Blown Off)


Etainf
Delta Eta
Conv Test
Max Iter
Print Eta

15
0.1
1.00E-10
30
0.1

Etainf
Delta Eta
Conv Test
Max Iter
Print Eta

3
0.005
1.00E-10
30
0.05

Etainf

Delta Eta
Conv Test
Max Iter
Print Eta

0.005
1.00E-10
30
0.1

You might also like