Menu
Home
People
Places
Arts
History
Plants & Animals
Science
Life & Culture
Technology
Reference.org
Blasius boundary layer
open-in-new
<h2 id="prandtls-boundary-layer-equations">Prandtl's boundary layer equations</h2> <p>Using scaling arguments, <a href="/facts/Ludwig_Prandtl/I4b4E7sq">Ludwig Prandtl</a><a class="footnote-ref" id="fnref:1" href="#fn:1"><sup>1</sup></a> argued that about half of the terms in the <a href="/facts/Navier-Stokes_equations/48pQugEa">Navier-Stokes equations</a> are negligible in boundary layer flows (except in a small region near the leading edge of the plate). This leads to a reduced set of equations known as the <a href="/facts/Boundary_layer/go8gqWM1">boundary layer equations</a>. For steady incompressible flow with constant viscosity and density, these read: </p> <ul><li>Mass Continuity: ∂ u ∂ x + ∂ v ∂ y = 0 {\displaystyle {\dfrac {\partial u}{\partial x}}+{\dfrac {\partial v}{\partial y}}=0} </li> <li> x {\displaystyle x} -Momentum: u ∂ u ∂ x + v ∂ u ∂ y = − 1 ρ ∂ p ∂ x + ν ∂ 2 u ∂ y 2 {\displaystyle u{\dfrac {\partial u}{\partial x}}+v{\dfrac {\partial u}{\partial y}}=-{\dfrac {1}{\rho }}{\dfrac {\partial p}{\partial x}}+{\nu }{\dfrac {\partial ^{2}u}{\partial y^{2}}}} </li> <li> y {\displaystyle y} -Momentum: 0 = − ∂ p ∂ y {\displaystyle 0=-{\dfrac {\partial p}{\partial y}}} </li></ul> <p>Here the coordinate system is chosen with x {\displaystyle x} pointing parallel to the plate in the direction of the flow and the y {\displaystyle y} coordinate pointing normal to the plate, u {\displaystyle u} and v {\displaystyle v} are the x {\displaystyle x} and y {\displaystyle y} velocity components, p {\displaystyle p} is the <a href="/facts/Pressure/eyc37GRG">pressure</a>, ρ {\displaystyle \rho } is the <a href="/facts/Density/92p6hh9I">density</a> and ν {\displaystyle \nu } is the <a href="/facts/Viscosity/QLVfrdCz">kinematic viscosity</a>. </p><p>A number of similarity solutions to this set of equations have been found for various types of flow, including flow on a thin flat-plate. The term <i>similarity</i> refers to the property that the velocity profiles at different positions in the flow are the same apart from scaling factors. Similarity scaling factors reduce the set of partial differential equations to a relatively easily solved set of non-linear ordinary differential equations. <a href="/facts/Paul_Richard_Heinrich_Blasius/ENWkX3sK">Paul Richard Heinrich Blasius</a>,<a class="footnote-ref" id="fnref:2" href="#fn:2"><sup>2</sup></a> one of Prandtl's students, developed the similarity model corresponding to the flow for the case where the pressure gradient, ∂ p {\displaystyle {\partial p}} / ∂ x {\displaystyle {\partial x}} , along a thin flat-plate is negligible compared to any pressure gradient in the boundary layer region. </p> <h2 id="blasius-equation---first-order-boundary-layer">Blasius equation - first-order boundary layer</h2> <p>Blasius showed that for the case where ∂ p / ∂ x = 0 {\displaystyle {\partial p}/{\partial x}=0} , the Prandtl x {\displaystyle x} -momentum equation has a self-similar solution. The self-similar solution exists because the equations and the boundary conditions are invariant under the transformation </p><p> x → c 2 x , y → c y , u → u , v → v c {\displaystyle x\rightarrow c^{2}x,\quad y\rightarrow cy,\quad u\rightarrow u,\quad v\rightarrow {\frac {v}{c}}} </p><p>where c {\displaystyle c} is any positive constant. He introduced the self-similar variables </p> <p> η = y δ ( x ) = y U ν x , ψ = ν U x f ( η ) {\displaystyle \eta ={\dfrac {y}{\delta (x)}}=y{\sqrt {\dfrac {U}{\nu x}}},\quad \psi ={\sqrt {\nu Ux}}f(\eta )} </p><p>where δ ( x ) ∝ ν x / U {\textstyle \delta (x)\propto {\sqrt {\nu x/U}}} is the <a href="/facts/Boundary_layer_thickness/TVC5Szbm">boundary layer thickness</a>, U {\displaystyle U} is the free stream velocity, and ψ {\displaystyle \psi } is the <a href="/facts/Stream_function/nEoYTGAl">stream function</a>. The stream function is directly proportional to the normalized function, f ( η ) {\displaystyle f(\eta )} , which is only a function of the similarity thickness variable. This leads directly to the velocity components:<a class="footnote-ref" id="fnref:3" href="#fn:3"><sup>3</sup></a>: 136 </p><p> u ( x , y ) = ∂ ψ ∂ y = U f ′ ( η ) , v ( x , y ) = − ∂ ψ ∂ x = 1 2 ν U x [ η f ′ ( η ) − f ( η ) ] {\displaystyle u(x,y)={\dfrac {\partial \psi }{\partial y}}=Uf'(\eta ),\quad v(x,y)=-{\dfrac {\partial \psi }{\partial x}}={\frac {1}{2}}{\sqrt {\dfrac {\nu U}{x}}}[\eta f'(\eta )-f(\eta )]} </p><p>Where the prime denotes derivation with respect to η {\displaystyle \eta } . Substitution into the x {\displaystyle x} -momentum equation gives the Blasius equation </p><p> 2 f ( 3 ) + f ″ f = 0 {\displaystyle 2f^{(3)}+f''f=0} </p><p>The boundary conditions are the <a href="/facts/No-slip_condition/jg1xo685">no-slip condition</a>, the impermeability of the wall and the free stream velocity outside the boundary layer </p><p> u ( x , 0 ) = 0 → f ′ ( 0 ) = 0 v ( x , 0 ) = 0 → f ( 0 ) = 0 u ( x , ∞ ) = U → f ′ ( ∞ ) = 1 {\displaystyle {\begin{aligned}u(x,0)&=0&\rightarrow &&f'(0)&=0\\v(x,0)&=0&\rightarrow &&f(0)&=0\\u(x,\infty )&=U&\rightarrow &&f'(\infty )&=1\end{aligned}}} </p><p>This is a third-order non-linear <a href="/facts/Ordinary_differential_equation/ygKGQ8kD">ordinary differential equation</a> which can be solved numerically, e.g. with the <a href="/facts/Shooting_method/vlcGhGVB">shooting method</a>. The boundary condition at infinity is converted to f ″ ( 0 ) = 0.332043934904293 {\displaystyle f''(0)=0.332043934904293} . </p><p>With the solution for f {\displaystyle f} and its derivatives in hand, the Prandtl y {\displaystyle y} -momentum equation can be non-dimensionalized and rearranged to obtain the y {\displaystyle y} -pressure gradient, ∂ p {\displaystyle {\partial p}} / ∂ y {\displaystyle {\partial y}} , as<a class="footnote-ref" id="fnref:4" href="#fn:4"><sup>4</sup></a>: 46 </p><p> x 2 U 2 δ ∗ 1 ρ ∂ P ∂ y = 1 2 η f ( 3 ) + 1 2 f ″ − 1 4 f f ′ + 1 4 η f ′ 2 + 1 4 η f f ″ , {\displaystyle {\frac {x^{2}}{U^{2}\delta ^{*}}}{\frac {1}{\rho }}{\frac {\partial P}{\partial y}}\quad =\quad {\frac {1}{2}}\eta f^{(3)}\;+\;{\frac {1}{2}}f''\;-\;{\frac {1}{4}}ff'\;+\;{\frac {1}{4}}\eta f'^{2}\;+\;{\frac {1}{4}}\eta ff''\quad ,} where δ ∗ {\displaystyle \delta ^{*}} is the Blasius displacement thickness. </p><p>The Blasius normal velocity v ( x , y ) {\displaystyle v(x,y)} and the y {\displaystyle y} -pressure gradient asymptotes to a value of 0.86 and 0.43, respectively, at large η {\displaystyle \eta } -values whereas u ( x , y ) {\displaystyle u(x,y)} asymptotes to the free stream velocity U {\displaystyle U} . As η {\displaystyle \eta } goes to zero, the scaled y {\displaystyle y} -pressure gradient goes to 0.16603. </p> <p>The limiting form for small η ≪ 1 {\displaystyle \eta \ll 1} is </p><p> f ( η ) = 1 2 α η 2 + O ( η 5 ) , α = 0.332057336215196 {\displaystyle f(\eta )={\frac {1}{2}}\alpha \eta ^{2}+O(\eta ^{5}),\qquad \alpha =0.332057336215196} </p><p>and the limiting form for large η ≫ 1 {\displaystyle \eta \gg 1} is<a class="footnote-ref" id="fnref:5" href="#fn:5"><sup>5</sup></a> </p><p> f ( η ) = η − β + O ( ( η − β ) − 2 e − 1 2 ( η − β ) 2 ) , β = 1.7207876575205 {\displaystyle f(\eta )=\eta -\beta +O\left((\eta -\beta )^{-2}e^{-{\frac {1}{2}}(\eta -\beta )^{2}}\right),\qquad \beta =1.7207876575205} </p><p>The characteristic parameters for boundary layers are the two sigma viscous boundary layer thickness,<a class="footnote-ref" id="fnref:6" href="#fn:6"><sup>6</sup></a> δ v {\displaystyle \delta _{v}} , the displacement thickness δ ∗ {\displaystyle \delta ^{*}} , the momentum thickness θ {\displaystyle \theta } , the wall <a href="/facts/Shear_stress/aWWnzIhE">shear stress</a> τ w {\displaystyle \tau _{w}} and the drag force F {\displaystyle F} acting on a length l {\displaystyle l} of the plate. For the Blasius solution, they are given by δ 99 ≈ δ v = 5.29 ν x U δ ∗ = δ 1 = ∫ 0 ∞ ( 1 − u U ) d y = 1.72 ν x U θ = δ 2 = ∫ 0 ∞ u U ( 1 − u U ) d y = 0.665 ν x U τ w = μ ∂ u ∂ y | y = 0 = 0.332 ρ μ U 3 x F = 2 ∫ 0 l τ w d x = 1.328 ρ μ l U 3 {\displaystyle {\begin{aligned}\delta _{99}&\approx \delta _{v}=5.29{\sqrt {\frac {\nu x}{U}}}\\[1ex]\delta ^{*}&=\delta _{1}=\int _{0}^{\infty }\left(1-{\frac {u}{U}}\right)dy=1.72{\sqrt {\frac {\nu x}{U}}}\\[1ex]\theta &=\delta _{2}=\int _{0}^{\infty }{\frac {u}{U}}\left(1-{\frac {u}{U}}\right)dy=0.665{\sqrt {\frac {\nu x}{U}}}\\[1ex]\tau _{w}&=\mu \left.{\frac {\partial u}{\partial y}}\right|_{y=0}=0.332{\sqrt {\frac {\rho \mu U^{3}}{x}}}\\[1ex]F&=2\int _{0}^{l}\tau _{w}dx=1.328{\sqrt {\rho \mu lU^{3}}}\end{aligned}}} </p><p>The factor 2 {\displaystyle 2} in the drag force formula is to account both sides of the plate. </p><p>The <a href="/facts/Boundary_layer/go8gqWM1">Von Kármán Momentum integral</a> and the energy integral for Blasius profile reduce to </p><p> τ w ρ U 2 = ∂ δ 2 ∂ x + v w U 2 ε ρ U 3 = ∂ δ 3 ∂ x + v w U {\displaystyle {\begin{aligned}{\frac {\tau _{w}}{\rho U^{2}}}&={\frac {\partial \delta _{2}}{\partial x}}+{\frac {v_{w}}{U}}\\{\frac {2\varepsilon }{\rho U^{3}}}&={\frac {\partial \delta _{3}}{\partial x}}+{\frac {v_{w}}{U}}\end{aligned}}} </p><p>where τ w {\displaystyle \tau _{w}} is the wall shear stress, v w {\displaystyle v_{w}} is the wall injection/suction velocity, ε {\displaystyle \varepsilon } is the energy dissipation rate, δ 2 {\displaystyle \delta _{2}} is the momentum thickness and δ 3 {\displaystyle \delta _{3}} is the energy thickness. </p> <h3>Uniqueness of Blasius solution</h3> <p>The Blasius solution is not unique from a mathematical perspective,<a class="footnote-ref" id="fnref:7" href="#fn:7"><sup>7</sup></a>: 131 as <a href="/facts/Ludwig_Prandtl/I4b4E7sq">Ludwig Prandtl</a> himself noted it in his <a href="/facts/Boundary_layer/go8gqWM1">transposition theorem</a> and analyzed by series of researchers such as <a href="/facts/Keith_Stewartson/YWlhK9lK">Keith Stewartson</a>, <a href="/facts/Paul_A._Libby/Cu8QMFet">Paul A. Libby</a>.<a class="footnote-ref" id="fnref:8" href="#fn:8"><sup>8</sup></a> To this solution, any one of the infinite discrete set of eigenfunctions can be added, each of which satisfies the linearly perturbed equation with homogeneous conditions and exponential decay at infinity. The first of these eigenfunctions turns out to be the x {\displaystyle x} derivative of the first order Blasius solution, which represents the uncertainty in the effective location of the origin. </p> <h3>Second-order boundary layer</h3> <p>This boundary layer approximation predicts a non-zero vertical velocity far away from the wall, which needs to be accounted in next order outer inviscid layer and the corresponding inner boundary layer solution, which in turn will predict a new vertical velocity and so on. The vertical velocity at infinity for the first order boundary layer problem from the Blasius equation is </p><p> v = 0.86 ν U x {\displaystyle v=0.86{\sqrt {\frac {\nu U}{x}}}} </p><p>The solution for second order boundary layer is zero. The solution for outer inviscid and inner boundary layer are<a class="footnote-ref" id="fnref:9" href="#fn:9"><sup>9</sup></a>: 134 </p><p> ψ ( x , y ) ∼ { y − ν U x β ℜ x + i y , outer ν U x f ( η ) + 0 , inner {\displaystyle \psi (x,y)\sim {\begin{cases}y-{\sqrt {\frac {\nu }{Ux}}}\beta \ \Re {\sqrt {x+iy}},&{\text{outer }}\\{\sqrt {\nu Ux}}f(\eta )+0,&{\text{inner}}\end{cases}}} </p><p>Again as in the first order boundary problem, any one of the infinite set of eigensolution can be added to this solution. In all the solutions R e = U x / ν {\displaystyle Re=Ux/\nu } can be considered as a <a href="/facts/Reynolds_number/4sdungsf">Reynolds number</a>. </p> <h3>Third-order boundary layer</h3> <p>Since the second order inner problem is zero, the corresponding corrections to third order problem is null i.e., the third order outer problem is same as second order outer problem.<a class="footnote-ref" id="fnref:10" href="#fn:10"><sup>10</sup></a>: 139 The solution for third-order correction does not have an exact expression, but the inner boundary layer expansion is of the form, </p><p> ψ ( x , y ) ∼ 2 ν U x f ( η ) + 0 + ( ν U x ) 3 / 2 [ log ( U x ν ) x 2 f 32 ( η ) + 1 2 x f 31 ( η ) ] + ⋅ ⋅ ⋅ {\displaystyle \psi (x,y)\sim {\sqrt {2\nu Ux}}f(\eta )+0+\left({\frac {\nu }{Ux}}\right)^{3/2}\left[\log \left({\frac {Ux}{\nu }}\right){\sqrt {\frac {x}{2}}}f_{32}(\eta )+{\frac {1}{\sqrt {2x}}}f_{31}(\eta )\right]+\cdot \cdot \cdot } </p><p>where f 32 {\displaystyle f_{32}} is the first eigensolution of the first order boundary layer solution (which is x {\displaystyle x} derivative of the first order Blasius solution) and solution for f 31 {\displaystyle f_{31}} is nonunique and the problem is left with an undetermined constant. </p> <h2 id="blasius-boundary-layer-with-suction">Blasius boundary layer with suction</h2> <p>Suction is one of the common methods to postpone the boundary layer separation.<a class="footnote-ref" id="fnref:11" href="#fn:11"><sup>11</sup></a> Consider a uniform suction velocity at the wall v ( 0 ) = − V {\displaystyle v(0)=-V} . <a href="/facts/Bryan_Thwaites/hMLeEXgN">Bryan Thwaites</a><a class="footnote-ref" id="fnref:12" href="#fn:12"><sup>12</sup></a> showed that the solution for this problem is same as the Blasius solution without suction for distances very close to the leading edge. Introducing the transformation </p><p> ψ = 2 U ν x f ( ξ , η ) , ξ = V x 2 U ν , η = U 2 ν x y {\displaystyle \psi ={\sqrt {2U\nu x}}f(\xi ,\eta ),\quad \xi =V{\sqrt {\frac {x}{2U\nu }}},\quad \eta ={\sqrt {\frac {U}{2\nu x}}}y} </p><p>into the boundary layer equations leads to </p><p> u = U ∂ f ∂ η , v = − U ν 2 x ( f + ξ ∂ f ∂ ξ − η ∂ f ∂ η ) , {\displaystyle u=U{\frac {\partial f}{\partial \eta }},\quad v=-{\sqrt {\frac {U\nu }{2x}}}\left(f+\xi {\frac {\partial f}{\partial \xi }}-\eta {\frac {\partial f}{\partial \eta }}\right),} ∂ 3 f ∂ η 3 + f ∂ 2 f ∂ η 2 + ξ ( ∂ f ∂ ξ ∂ 2 f ∂ η 2 − ∂ 2 f ∂ ξ ∂ η ∂ f ∂ η ) = 0 {\displaystyle {\frac {\partial ^{3}f}{\partial \eta ^{3}}}+f{\frac {\partial ^{2}f}{\partial \eta ^{2}}}+\xi \left({\frac {\partial f}{\partial \xi }}{\frac {\partial ^{2}f}{\partial \eta ^{2}}}-{\frac {\partial ^{2}f}{\partial \xi \partial \eta }}{\frac {\partial f}{\partial \eta }}\right)=0} </p><p>with boundary conditions, </p><p> f ( ξ , 0 ) = ξ , ∂ f ∂ η ( ξ , 0 ) = 0 , ∂ f ∂ η ( ξ , ∞ ) = 0. {\displaystyle f(\xi ,0)=\xi ,\quad {\frac {\partial f}{\partial \eta }}(\xi ,0)=0,\quad {\frac {\partial f}{\partial \eta }}(\xi ,\infty )=0.} </p> <h3>Von Mises transformation</h3> <p>Iglisch obtained the complete numerical solution in 1944.<a class="footnote-ref" id="fnref:13" href="#fn:13"><sup>13</sup></a> If further <a href="/facts/Richard_von_Mises/LFVGPRad">von Mises</a> transformation<a class="footnote-ref" id="fnref:14" href="#fn:14"><sup>14</sup></a> is introduced </p><p> σ = 2 ξ , ψ − V x = U ν 2 V σ τ 2 , ϕ = 4 u 2 U 2 , χ = U 2 − u 2 = U 2 ( 1 − V 4 ) , {\displaystyle \sigma =2\xi ,\quad \psi -Vx={\frac {U\nu }{2V}}\sigma \tau ^{2},\quad \phi ={\frac {4u^{2}}{U^{2}}},\quad \chi =U^{2}-u^{2}=U^{2}\left(1-{\frac {V}{4}}\right),} </p><p>then the equations become </p><p> ϕ ∂ 2 ϕ ∂ τ 2 + ( 2 σ τ + τ 3 − ϕ τ ) ∂ ϕ ∂ τ = 2 σ τ 2 ∂ ϕ ∂ σ {\displaystyle {\sqrt {\phi }}{\frac {\partial ^{2}\phi }{\partial \tau ^{2}}}+\left(2\sigma \tau +\tau ^{3}-{\frac {\sqrt {\phi }}{\tau }}\right){\frac {\partial \phi }{\partial \tau }}=2\sigma \tau ^{2}{\frac {\partial \phi }{\partial \sigma }}} </p><p>with boundary conditions, </p><p> ϕ ( 0 , τ ) = 4 , ϕ ( σ , 0 ) = 0 , ϕ ( σ , ∞ ) = 4. {\displaystyle \phi (0,\tau )=4,\quad \phi (\sigma ,0)=0,\quad \phi (\sigma ,\infty )=4.} </p><p>This <a href="/facts/Parabolic_partial_differential_equation/gJQfRUlr">parabolic partial differential equation</a> can be marched starting from σ = 0 {\displaystyle \sigma =0} numerically. </p> <h3>Asymptotic suction profile</h3> <p>Since the convection due to suction and the diffusion due to the solid wall are acting in the opposite direction, the profile will reach steady solution at large distance, unlike the Blasius profile where boundary layer grows indefinitely. The solution was first obtained by <a href="/facts/Alan_Arnold_Griffith/6Hp56TQ4">Griffith</a> and F.W. Meredith.<a class="footnote-ref" id="fnref:15" href="#fn:15"><sup>15</sup></a> For distances from the leading edge of the plate x ≫ ν U / V 2 {\displaystyle x\gg \nu U/V^{2}} , both the boundary layer thickness and the solution are independent of x {\displaystyle x} given by </p><p> δ = ν V , u = U ( 1 − e − y V / ν ) , v = − V . {\displaystyle \delta ={\frac {\nu }{V}},\quad u=U(1-e^{-yV/\nu }),\quad v=-V.} </p><p><a href="/facts/Keith_Stewartson/YWlhK9lK">Stewartson</a><a class="footnote-ref" id="fnref:16" href="#fn:16"><sup>16</sup></a> studied matching of full solution to the asymptotic suction profile. </p> <h2 id="compressible-blasius-boundary-layer">Compressible Blasius boundary layer</h2> <p>Here Blasius boundary layer with a specified <a href="/facts/Specific_enthalpy/bBou8lCE">specific enthalpy</a> h {\displaystyle h} at the wall is studied. The <a href="/facts/Density/92p6hh9I">density</a> ρ {\displaystyle \rho } , <a href="/facts/Viscosity/QLVfrdCz">viscosity</a> μ {\displaystyle \mu } and thermal conductivity κ {\displaystyle \kappa } are no longer constant here. The equation for conservation of mass, momentum and energy become </p><p> ∂ ( ρ u ) ∂ x + ∂ ( ρ v ) ∂ y = 0 , ρ ( u ∂ u ∂ x + v ∂ u ∂ y ) = ∂ ∂ y ( μ ∂ u ∂ y ) , ρ ( u ∂ h ∂ x + v ∂ h ∂ y ) = ∂ ∂ y ( μ P r ∂ h ∂ y ) + μ ( ∂ u ∂ y ) 2 {\displaystyle {\begin{aligned}{\frac {\partial (\rho u)}{\partial x}}+{\frac {\partial (\rho v)}{\partial y}}&=0,\\\rho \left(u{\frac {\partial u}{\partial x}}+v{\frac {\partial u}{\partial y}}\right)&={\frac {\partial }{\partial y}}\left(\mu {\frac {\partial u}{\partial y}}\right),\\\rho \left(u{\frac {\partial h}{\partial x}}+v{\frac {\partial h}{\partial y}}\right)&={\frac {\partial }{\partial y}}\left({\frac {\mu }{Pr}}{\frac {\partial h}{\partial y}}\right)+\mu \left({\frac {\partial u}{\partial y}}\right)^{2}\end{aligned}}} </p><p>where P r = c p ∞ μ ∞ / κ ∞ {\displaystyle Pr=c_{p_{\infty }}\mu _{\infty }/\kappa _{\infty }} is the <a href="/facts/Prandtl_number/UUPDVoeB">Prandtl number</a> with suffix ∞ {\displaystyle \infty } representing properties evaluated at infinity. The boundary conditions become </p><p> u = v = h − h w ( x ) = 0 for y = 0 , {\displaystyle u=v=h-h_{w}(x)=0\ {\text{for}}\ y=0,} u − U = h − h ∞ = 0 for y = ∞ or x = 0. {\displaystyle u-U=h-h_{\infty }=0\ {\text{for}}\ y=\infty \ {\text{or}}\ x=0.} </p><p>Unlike the incompressible boundary layer, similarity solution exists only if the transformation </p><p> x → c 2 x , y → c y , u → u , v → v c , h → h , ρ → ρ , μ → μ {\displaystyle x\rightarrow c^{2}x,\quad y\rightarrow cy,\quad u\rightarrow u,\quad v\rightarrow {\frac {v}{c}},\quad h\rightarrow h,\quad \rho \rightarrow \rho ,\quad \mu \rightarrow \mu } </p><p>holds and this is possible only if h w = constant {\displaystyle h_{w}={\text{constant}}} . </p> <h3>Howarth transformation</h3> <p>Introducing the self-similar variables using <a href="/facts/Howarth%E2%80%93Dorodnitsyn_transformation/aw1cCqPG">Howarth–Dorodnitsyn transformation</a> </p><p> η = U 2 ρ ∞ μ ∞ U x ∫ 0 y ρ d y , f ( η ) = ψ 2 ρ ∞ μ ∞ U x , h ~ ( η ) = h h ∞ , h ~ w = h w h ∞ , ρ ~ = ρ ρ ∞ , μ ~ = μ μ ∞ {\displaystyle \eta ={\frac {U}{\sqrt {2\rho _{\infty }\mu _{\infty }Ux}}}\int _{0}^{y}\rho \,dy,\quad f(\eta )={\frac {\psi }{\sqrt {2\rho _{\infty }\mu _{\infty }Ux}}},\quad {\tilde {h}}(\eta )={\frac {h}{h_{\infty }}},\quad {\tilde {h}}_{w}={\frac {h_{w}}{h_{\infty }}},\quad {\tilde {\rho }}={\frac {\rho }{\rho _{\infty }}},\quad {\tilde {\mu }}={\frac {\mu }{\mu _{\infty }}}} </p><p>the equations reduce to ( ρ ~ μ ~ f ″ ) ′ + f f ″ = 0 , ( ρ ~ μ ~ h ~ ′ ) ′ + P r f h ~ ′ + P r ( γ − 1 ) M 2 ρ ~ μ ~ f ″ 2 = 0 {\displaystyle {\begin{aligned}({\tilde {\rho }}{\tilde {\mu }}f'')'+ff''=0,\\({\tilde {\rho }}{\tilde {\mu }}{\tilde {h}}')'+Prf{\tilde {h}}'+Pr(\gamma -1)M^{2}{\tilde {\rho }}{\tilde {\mu }}f''^{2}=0\end{aligned}}} where γ {\displaystyle \gamma } is the <a href="/facts/Specific_heat_ratio/jFcMaukb">specific heat ratio</a> and M = U / c ∞ {\displaystyle M=U/c_{\infty }} is the <a href="/facts/Mach_number/Ee01MIgw">Mach number</a>, where c ∞ {\displaystyle c_{\infty }} is the <a href="/facts/Speed_of_sound/dYtbt1nf">speed of sound</a>. The equation can be solved once ρ ~ = ρ ~ ( h ~ ) , μ ~ = μ ~ ( h ~ ) {\displaystyle {\tilde {\rho }}={\tilde {\rho }}({\tilde {h}}),\ {\tilde {\mu }}={\tilde {\mu }}({\tilde {h}})} are specified. The boundary conditions are </p><p> f ( 0 ) = f ′ ( 0 ) = θ ( 0 ) − h ~ w = f ′ ( ∞ ) − 1 = h ~ ( ∞ ) − 1 = 0. {\displaystyle f(0)=f'(0)=\theta (0)-{\tilde {h}}_{w}=f'(\infty )-1={\tilde {h}}(\infty )-1=0.} </p><p>The commonly used expressions for air are γ = 1.4 , P r = 0.7 , ρ ~ = h ~ − 1 , μ ~ = h ~ 2 / 3 {\displaystyle \gamma =1.4,\ Pr=0.7,\ {\tilde {\rho }}={\tilde {h}}^{-1},\ {\tilde {\mu }}={\tilde {h}}^{2/3}} . If c p {\displaystyle c_{p}} is constant, then h ~ = θ ~ = T / T ∞ {\displaystyle {\tilde {h}}={\tilde {\theta }}=T/T_{\infty }} . The temperature inside the boundary layer will increase even though the plate temperature is maintained at the same temperature as ambient, due to dissipative heating and of course, these dissipation effects are only pronounced when the <a href="/facts/Mach_number/Ee01MIgw">Mach number</a> M {\displaystyle M} is large. </p> <h2 id="first-order-blasius-boundary-layer-in-parabolic-coordinates">First-order Blasius boundary layer in parabolic coordinates</h2> <p>Since the boundary layer equations are <a href="/facts/Parabolic_partial_differential_equation/gJQfRUlr">Parabolic partial differential equation</a>, the natural coordinates for the problem is <a href="/facts/Parabolic_coordinates/MLledsNu">parabolic coordinates</a>.<a class="footnote-ref" id="fnref:17" href="#fn:17"><sup>17</sup></a>: 142 The transformation from <a href="/facts/Cartesian_coordinates/lkTqvMmB">Cartesian coordinates</a> ( x , y ) {\displaystyle (x,y)} to <a href="/facts/Parabolic_coordinates/MLledsNu">parabolic coordinates</a> ( ξ , η ) {\displaystyle (\xi ,\eta )} is given by </p><p> x + i y = 1 2 ( ξ + i η ) 2 , x = 1 2 ( ξ 2 − η 2 ) , y = ξ η . {\displaystyle x+iy={\tfrac {1}{2}}(\xi +i\eta )^{2},\quad x={\tfrac {1}{2}}(\xi ^{2}-\eta ^{2}),\quad y=\xi \eta .} </p> <h2 id="see-also">See also</h2> <ul><li><a href="/facts/Falkner%E2%80%93Skan_boundary_layer/a0pXb55H">Falkner–Skan boundary layer</a></li> <li><a href="/facts/Emmons_problem/Xg4dDYv3">Emmons problem</a></li></ul> <h2 id="external-links">External links</h2> <ul><li><a href="https://ntrs.nasa.gov/citations/20050028493">[1]</a> - English translation of Blasius' original paper - NACA Technical Memorandum 1256.</li></ul> <h2 id="footnotes">Footnotes</h2> <ul><li>Parlange, J. Y.; Braddock, R. D.; Sander, G. (1981). "Analytical approximations to the solution of the Blasius equation". <i>Acta Mech</i>. 38 (1–2): 119–125. <a href="/facts/Bibcode_(identifier)/9HtdQSGB">Bibcode</a>:<a href="https://ui.adsabs.harvard.edu/abs/1981AcMec..38..119P">1981AcMec..38..119P</a>. <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1007%2FBF01351467">10.1007/BF01351467</a>. <a href="/facts/S2CID_(identifier)/ldJsHa2Y">S2CID</a> <a href="https://api.semanticscholar.org/CorpusID:122114667">122114667</a>.</li> <li>Pozrikidis, C. (1998). <i>Introduction to Theoretical and Computational Fluid Dynamics</i>. Oxford. <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0-19-509320-9.</li> <li>Schlichting, H. (2004). <i>Boundary-Layer Theory</i>. Springer. <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-3-540-66270-9.</li> <li>Wilcox, David C. <i>Basic Fluid Mechanics</i> DCW Industries Inc. 2007</li> <li>Boyd, John P. (1999), <a href="http://projecteuclid.org/getRecord?id=euclid.em/1047262359">"The Blasius function in the complex plane"</a>, <i>Experimental Mathematics</i>, 8 (4): 381–394, <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1080%2F10586458.1999.10504626">10.1080/10586458.1999.10504626</a>, <a href="/facts/ISSN_(identifier)/DPAflDvU">ISSN</a> <a href="https://search.worldcat.org/issn/1058-6458">1058-6458</a>, <a href="/facts/MR_(identifier)/uP137L11">MR</a> <a href="https://mathscinet.ams.org/mathscinet-getitem?mr=1737233">1737233</a></li></ul> <h2 id="references">References</h2> <ol> <li id="fn:1"><p>Prandtl, L. (1904). "Über Flüssigkeitsbewegung bei sehr kleiner Reibung". Verhandlinger 3. Int. Math. Kongr. Heidelberg: 484–491. <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li> <li id="fn:2"><p>Blasius, H. (1908). "Grenzschichten in Flüssigkeiten mit kleiner Reibung". Z. Angew. Math. Phys. 56: 1–37. <a href="#fnref:2" class="footnote-back-ref">↩</a></p></li> <li id="fn:3"><p>Schlichting, H., (1979). Boundary-Layer Theory, 7th ed., McGraw-Hill, New York. <a href="#fnref:3" class="footnote-back-ref">↩</a></p></li> <li id="fn:4"><p>Weyburne, David (2022). Aspects of Boundary Layer Theory. ISBN 978-0-578-98334-9. <a href="978-0-578-98334-9" target="_blank">978-0-578-98334-9</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></p></li> <li id="fn:5"><p>Boyd, J. (2008). "The Blasius function: computations before computers, the value of tricks, undergraduate projects, and open research problems". SIAM Rev. 50 (4): 791–804. Bibcode:2008SIAMR..50..791B. doi:10.1137/070681594. <a href="/wiki/Bibcode_(identifier)" target="_blank">/wiki/Bibcode_(identifier)</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></p></li> <li id="fn:6"><p>Weyburne, D. (2014). "New thickness and shape parameters for the boundary layer velocity profile". Experimental Thermal and Fluid Science. 54: 22–28. doi:10.1016/j.expthermflusci.2014.01.008. <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></p></li> <li id="fn:7"><p>Van Dyke, Milton (1975). Perturbation methods in fluid mechanics. Parabolic Press. ISBN 9780915760015. <a href="9780915760015" target="_blank">9780915760015</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></p></li> <li id="fn:8"><p>Libby, Paul A., and Herbert Fox. "Some perturbation solutions in laminar boundary-layer theory." Journal of Fluid Mechanics 17.3 (1963): 433-449. <a href="#fnref:8" class="footnote-back-ref">↩</a></p></li> <li id="fn:9"><p>Van Dyke, Milton (1975). Perturbation methods in fluid mechanics. Parabolic Press. ISBN 9780915760015. <a href="9780915760015" target="_blank">9780915760015</a> <a href="#fnref:9" class="footnote-back-ref">↩</a></p></li> <li id="fn:10"><p>Van Dyke, Milton (1975). Perturbation methods in fluid mechanics. Parabolic Press. ISBN 9780915760015. <a href="9780915760015" target="_blank">9780915760015</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></p></li> <li id="fn:11"><p>Rosenhead, Louis, ed. Laminar boundary layers. Clarendon Press, 1963. <a href="#fnref:11" class="footnote-back-ref">↩</a></p></li> <li id="fn:12"><p>Thwaites, Bryan. On certain types of boundary-layer flow with continuous surface suction. HM Stationery Office, 1946. <a href="#fnref:12" class="footnote-back-ref">↩</a></p></li> <li id="fn:13"><p>Iglisch, Rudolf. Exakte Berechnung der laminaren Grenzschicht an der längsangeströmten ebenen Platte mit homogener Absaugung. Oldenbourg, 1944. <a href="#fnref:13" class="footnote-back-ref">↩</a></p></li> <li id="fn:14"><p>Von Mises, Richard. "Bemerkungen zur hydrodynamik." Z. Angew. Math. Mech 7 (1927): 425-429. <a href="#fnref:14" class="footnote-back-ref">↩</a></p></li> <li id="fn:15"><p>Griffith, A. A., and F. W. Meredith. "The possible improvement in aircraft performance due to the use of boundary layer suction." Royal Aircraft Establishment Report No. E 3501 (1936): 12. <a href="#fnref:15" class="footnote-back-ref">↩</a></p></li> <li id="fn:16"><p>Stewartson, K. "On asymptotic expansions in the theory of boundary layers." Studies in Applied Mathematics 36.1-4 (1957): 173-191. <a href="#fnref:16" class="footnote-back-ref">↩</a></p></li> <li id="fn:17"><p>Van Dyke, Milton (1975). Perturbation methods in fluid mechanics. Parabolic Press. ISBN 9780915760015. <a href="9780915760015" target="_blank">9780915760015</a> <a href="#fnref:17" class="footnote-back-ref">↩</a></p></li> </ol>