Английская Википедия:Bickley jet
In fluid dynamics, Bickley jet is a steady two-dimensional laminar plane jet with large jet Reynolds number emerging into the fluid at rest, named after W. G. Bickley, who gave the analytical solution in 1937,[1] to the problem derived by Schlichting in 1933[2] and the corresponding problem in axisymmetric coordinates is called as Schlichting jet. The solution is valid only for distances far away from the jet origin.
Flow description
Consider a steady plane emerging into the same fluid, a type of submerged jets from a narrow slit, which is supposed to be very small (such that the fluid loses memory of the shape and size of the slit far away from the origin, it remembers only the net momentum flux). Let the velocity be <math>(u,v)</math> in Cartesian coordinate and the axis of the jet be <math>x</math> axis with origin at the orifice. The flow is self-similar for large Reynolds number (the jet is so thin that <math>u(x,y)</math> varies much more rapidly in the transverse <math>y</math> direction than the streamwise <math>x</math> direction) and can be approximated with boundary layer equations.
- <math>
\begin{align} \frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} &= 0,\\ u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} &= \nu \frac{\partial^2 u}{\partial y^2}, \end{align} </math>
where <math>\nu</math> is the kinematic viscosity and the pressure is everywhere equal to the outside fluid pressure. Since the fluid is at rest far away from the center of the jet
- <math> u \rightarrow 0</math> as <math>y\rightarrow\pm\infty</math>,
and because the flow is symmetric about <math>x</math> axis
- <math>v=0</math> at <math>y=0</math>,
and also since there is no solid boundary and the pressure is constant, the momentum flux <math>M</math> across any plane normal to the <math>x</math> axis must be the same
- <math>M = 2\rho \int_0^\infty u^2 \, dy</math>
is a constant, where <math>\rho</math> which also constant for incompressible flow.
Proof of constant axial momentum flux
The constant momentum flux condition can be obtained by integrating the momentum equation across the jet.
<math display="block"> \begin{align} & \int_{-\infty}^\infty u \frac{\partial u}{\partial x} \, dy + \int_{-\infty}^\infty v \frac{\partial u}{\partial y} \, dy = \left[\nu \frac{\partial u}{\partial y}\right]_{-\infty}^\infty, \\[10pt] & \frac{d}{dx}\int_{-\infty}^\infty u^2 \, dy = 0, \quad \Rightarrow \quad \int_0^\infty u^2 \, dy = \text{constant}. \end{align} </math>
where <math display="block">v \frac{\partial u }{\partial y} = \frac{\partial (uv)}{\partial y} - u \frac{\partial v}{\partial y}= \frac{\partial (uv)}{\partial y} + u \frac{\partial u}{\partial x}</math> is used to simplify the above equation. The mass flux <math>Q</math> across any cross section normal to the <math>x</math> axis is not constant, because there is a slow entrainment of outer fluid into the jet, and it's a part of the boundary layer solution. This can be easily verified by integrating the continuity equation across the boundary layer.
<math display="block"> \begin{align} \int_{-\infty}^\infty \frac{\partial u}{\partial x} \, dy + \int_{-\infty}^\infty \frac{\partial v}{\partial y} \, dy &= 0,\\[8pt]
\frac{d}{dx} \int_{-\infty}^\infty u \, dy = - \Big[ v \Big]_{-\infty}^\infty &= - 2 v(x,\infty).
\end{align} </math>
where symmetry condition <math>v(x,-\infty) = -v(x,\infty)</math> is used.[3][4]
Self-similar solution
The self-similar solution is obtained by introducing the transformation <math display="block">\eta = \frac{y}{x^{2/3}}, \quad u = \frac{6\nu}{x^{1/3}} F'(\eta), \quad v = \frac{2\nu}{x^{2/3}} (2\eta F'(\eta) - F(\eta))</math> the equation reduces to <math display="block">F' + 2FF + 2F'^2 = 0,</math> while the boundary conditions become <math display="block">F'(\pm\infty) = 0, \quad F(0) = 0, \quad M = 72 \nu^2 \rho \int_0^\infty F'^2 \, d\eta.</math>
The exact solution is given by <math display="block">F(\eta) = \alpha \tanh \alpha\eta</math> where <math>\alpha</math> is solved from the following equation <math display="block">M = 72 \nu^2 \rho \int_0^\infty \operatorname{sech}^4\eta \, d\eta = 48 \nu^2\rho\alpha^3, \quad \Rightarrow \quad \alpha = \left(\frac M {48\nu^2\rho} \right)^{1/3}. </math>
Letting <math display="block">\xi = \alpha \eta = 0.2752\left(\frac M {\nu^2\rho}\right)^{1/3} \frac y {x^{2/3}}, </math>
the velocity is given by <math display="block"> \begin{align} u &= 0.4543 \left(\frac{M^2}{\nu\rho^2 x}\right)^{1/3} \operatorname{sech}^2\xi,\\ v &= 0.5503 \left(\frac{M\nu}{\rho x^2}\right)^{1/3} (2\xi \operatorname{sech}^2 \xi - \tanh \xi). \end{align} </math>
The mass flow rate <math>Q</math> across a plane at a distance <math>x</math> from the orifice normal to the jet is[5][6][7] <math display="block">Q = 2\rho \int_0^\infty u \, dy = 3.3019 \left(M\nu \rho^2 x\right)^{1/3}.</math>
See also
References
- ↑ Bickley, W. G. "LXXIII. The plane jet." The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 23.156 (1937): 727-731.(Original paper:http://www.tandfonline.com/doi/abs/10.1080/14786443708561847?journalCode=tphm18)
- ↑ Schlichting, Hermann. "Laminare strahlausbreitung." ZAMM‐Journal of Applied Mathematics and Mechanics/Zeitschrift für Angewandte Mathematik und Mechanik 13.4 (1933): 260-263.
- ↑ Kundu, P. K., and L. M. Cohen. "Fluid mechanics, 638 pp." Academic, Calif (1990).
- ↑ Pozrikidis, Costas, and Joel H. Ferziger. "Introduction to theoretical and computational fluid dynamics." (1997): 72–74.
- ↑ Rosenhead, Louis, ed. Laminar boundary layers. Clarendon Press, 1963.
- ↑ Acheson, David J. Elementary fluid dynamics. Oxford University Press, 1990.
- ↑ Drazin, Philip G., and Norman Riley. The Navier–Stokes equations: a classification of flows and exact solutions. No. 334. Cambridge University Press, 2006.