1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
{
"metadata": {
"name": "",
"signature": "sha256:8adfc21cf4552e14b563ca7c0f2e56ece56cd7f1691cf9486f2cfbd789df7be1"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "heading",
"level": 1,
"metadata": {},
"source": [
"Chapter 11: Boundary Layer"
]
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 11.1, Page 383"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from __future__ import division\n",
"import math\n",
"\n",
"\n",
" #Initializing the variables\n",
"rho = 860;\n",
"v = 10**-5;\n",
"Us = 3;\n",
"b = 1.25;\n",
"l = 2;\n",
"\n",
" #Calculations\n",
"x = 1; # At x =1 \n",
"Rex = Us*x/v;\n",
"ReL = Us*l/v ;\n",
"mu = rho*v;\n",
"T0 = 0.332*mu*Us/x*Rex**0.5;\n",
"Cf = 1.33*ReL**-0.5;\n",
"F = rho*Us**2*l*b*Cf ;\n",
"print \"shear stress (N/m^2) :\" ,round(T0,1)\n",
"print \"Total, double-sided resistance of the plate (N) :\",round(F,3)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"shear stress (N/m^2) : 4.7\n",
"Total, double-sided resistance of the plate (N) : 33.224\n"
]
}
],
"prompt_number": 1
},
{
"cell_type": "heading",
"level": 2,
"metadata": {},
"source": [
"Example 11.2, Page 387"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from __future__ import division\n",
"import math\n",
"\n",
"\n",
" #Initializing the variables\n",
"Us = 6;\n",
"b = 3;\n",
"l = 30;\n",
"rho = 1000;\n",
"mu = 10**-3;\n",
"T = 20+273; # Temperature in Kelvin\n",
"\n",
" #Calculations\n",
" \n",
"ReL = rho*Us*l/mu; \n",
"Cf = 0.455*math.log10(ReL)**-2.58 ;\n",
"\n",
"F = rho*Us**2*l*b*Cf ;\n",
"Lt = 10**5*mu/(rho*Us); # Assuming transition at Rel = 10**5\n",
"\n",
"print \"Drag Force (kN) :\", round(F/1000,2)\n",
"print \"Transition length (m) :\",round(Lt,4)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Drag Force (kN) : 6.36\n",
"Transition length (m) : 0.0167\n"
]
}
],
"prompt_number": 2
}
],
"metadata": {}
}
]
}
|
