Misplaced Pages

Elementary flow

Article snapshot taken from[REDACTED] with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
This article includes a list of general references, but it lacks sufficient corresponding inline citations. Please help to improve this article by introducing more precise citations. (February 2018) (Learn how and when to remove this message)

In the larger context of the Navier-Stokes equations (and especially in the context of potential theory), elementary flows are basic flows that can be combined, using various techniques, to construct more complex flows. In this article the term "flow" is used interchangeably with the term "solution" due to historical reasons.

The techniques involved to create more complex solutions can be for example by superposition, by techniques such as topology or considering them as local solutions on a certain neighborhood, subdomain or boundary layer and to be patched together. Elementary flows can be considered the basic building blocks (fundamental solutions, local solutions and solitons) of the different types of equations derived from the Navier-Stokes equations. Some of the flows reflect specific constraints such as incompressible or irrotational flows, or both, as in the case of potential flow, and some of the flows may be limited to the case of two dimensions.

Due to the relationship between fluid dynamics and field theory, elementary flows are relevant not only to aerodynamics but to all field theory in general. To put it in perspective boundary layers can be interpreted as topological defects on generic manifolds, and considering fluid dynamics analogies and limit cases in electromagnetism, quantum mechanics and general relativity one can see how all these solutions are at the core of recent developments in theoretical physics such as the ads/cft duality, the SYK model, the physics of nematic liquids, strongly correlated systems and even to quark gluon plasmas.

Two-dimensional uniform flow

Uniform
Potential flow streamlines for an ideal uniform flow

For steady-state, spatially uniform flow of a fluid in the xy plane, the velocity vector is

v = v 0 cos ( θ 0 ) e x + v 0 sin ( θ 0 ) e y {\displaystyle \mathbf {v} =v_{0}\cos(\theta _{0})\,\mathbf {e} _{x}+v_{0}\sin(\theta _{0})\,\mathbf {e} _{y}}

where

v 0 {\displaystyle v_{0}} is the absolute magnitude of the velocity (i.e., v 0 = | v | {\displaystyle v_{0}=|\mathbf {v} |} );
θ 0 {\displaystyle \theta _{0}} is the angle the velocity vector makes with the positive x axis ( θ 0 {\displaystyle \theta _{0}} is positive for angles measured in a counterclockwise sense from the positive x axis); and
e x {\displaystyle \mathbf {e} _{x}} and e y {\displaystyle \mathbf {e} _{y}} are the unit basis vectors of the xy coordinate system.

Because this flow is incompressible (i.e., v = 0 {\displaystyle \nabla \cdot \mathbf {v} =0} ) and two-dimensional, its velocity can be expressed in terms of a stream function, ψ {\displaystyle \psi } :

v x = ψ y {\displaystyle v_{x}={\frac {\partial \psi }{\partial y}}}
v y = ψ x {\displaystyle v_{y}=-{\frac {\partial \psi }{\partial x}}}

where

ψ = ψ 0 v 0 sin ( θ 0 ) x + v 0 cos ( θ 0 ) y {\displaystyle \psi =\psi _{0}-v_{0}\sin(\theta _{0})\,x+v_{0}\cos(\theta _{0})\,y}

and ψ 0 {\displaystyle \psi _{0}} is a constant.

In cylindrical coordinates:

v r = 1 r ψ θ {\displaystyle v_{r}=-{\frac {1}{r}}{\frac {\partial \psi }{\partial \theta }}}
v θ = ψ r {\displaystyle v_{\theta }={\frac {\partial \psi }{\partial r}}}

and

ψ = ψ 0 + v 0 r sin ( θ θ 0 ) {\displaystyle \psi =\psi _{0}+v_{0}\,r\sin(\theta -\theta _{0})}

This flow is irrotational (i.e., × v = 0 {\displaystyle \nabla \times \mathbf {v} =\mathbf {0} } ) so its velocity can be expressed in terms of a potential function, ϕ {\displaystyle \phi } :

v x = ϕ x {\displaystyle v_{x}=-{\frac {\partial \phi }{\partial x}}}
v y = ϕ y {\displaystyle v_{y}=-{\frac {\partial \phi }{\partial y}}}

where

ϕ = ϕ 0 v 0 cos ( θ 0 ) x v 0 sin ( θ 0 ) y {\displaystyle \phi =\phi _{0}-v_{0}\cos(\theta _{0})\,x-v_{0}\sin(\theta _{0})\,y}

and ϕ 0 {\displaystyle \phi _{0}} is a constant.

In cylindrical coordinates

v r = ϕ r {\displaystyle v_{r}={\frac {\partial \phi }{\partial r}}}
v θ = 1 r ϕ θ {\displaystyle v_{\theta }={\frac {1}{r}}{\frac {\partial \phi }{\partial \theta }}}
ϕ = ϕ 0 v 0 r cos ( θ θ 0 ) {\displaystyle \phi =\phi _{0}-v_{0}\,r\cos(\theta -\theta _{0})}

Two-dimensional line source

Point-source
Potential flow streamlines for an ideal line source

The case of a vertical line emitting at a fixed rate a constant quantity of fluid Q per unit length is a line source. The problem has a cylindrical symmetry and can be treated in two dimensions on the orthogonal plane.

Line sources and line sinks (below) are important elementary flows because they play the role of monopole for incompressible fluids (which can also be considered examples of solenoidal fields i.e. divergence free fields). Generic flow patterns can be also de-composed in terms of multipole expansions, in the same manner as for electric and magnetic fields where the monopole is essentially the first non-trivial (e.g. constant) term of the expansion.

This flow pattern is also both irrotational and incompressible.

This is characterized by a cylindrical symmetry:

v = v r ( r ) e r {\displaystyle \mathbf {v} =v_{r}(r)\mathbf {e} _{r}}

Where the total outgoing flux is constant

S v d S = 0 2 π ( v r ( r ) e r ) ( e r r d θ ) = 2 π r v r ( r ) = Q {\displaystyle \int _{S}\mathbf {v} \cdot d\mathbf {S} =\int _{0}^{2\pi }(v_{r}(r)\,\mathbf {e} _{r})\cdot (\mathbf {e} _{r}\,r\,d\theta )=\!2\pi \,r\,v_{r}(r)=Q}

Therefore,

v r = Q 2 π r {\displaystyle v_{r}={\frac {Q}{2\pi r}}}

This is derived from a stream function

ψ ( r , θ ) = Q 2 π θ {\displaystyle \psi (r,\theta )=-{\frac {Q}{2\pi }}\theta }

or from a potential function

ϕ ( r , θ ) = Q 2 π ln r {\displaystyle \phi (r,\theta )=-{\frac {Q}{2\pi }}\ln r}

Two-dimensional line sink

The case of a vertical line absorbing at a fixed rate a constant quantity of fluid Q per unit length is a line sink. Everything is the same as the case of a line source a part from the negative sign.

v r = Q 2 π r {\displaystyle v_{r}=-{\frac {Q}{2\pi r}}}

This is derived from a stream function

ψ ( r , θ ) = Q 2 π θ {\displaystyle \psi (r,\theta )={\frac {Q}{2\pi }}\theta }

or from a potential function

ϕ ( r , θ ) = Q 2 π ln r {\displaystyle \phi (r,\theta )={\frac {Q}{2\pi }}\ln r}

Given that the two results are the same a part from a minus sign we can treat transparently both line sources and line sinks with the same stream and potential functions permitting Q to assume both positive and negative values and absorbing the minus sign into the definition of Q.

Two-dimensional doublet or dipole line source

Potential flow streamlines for an ideal doublet, or dipole, line

If we consider a line source and a line sink at a distance d we can reuse the results above and the stream function will be

ψ ( r ) = ψ Q ( r d / 2 ) ψ Q ( r + d / 2 )   d ψ Q ( r ) {\displaystyle \psi (\mathbf {r} )=\psi _{Q}(\mathbf {r} -\mathbf {d} /2)-\psi _{Q}(\mathbf {r} +\mathbf {d} /2)\ \simeq \mathbf {d} \cdot \nabla \psi _{Q}(\mathbf {r} )}

The last approximation is to the first order in d.

Given

d = d [ cos ( θ 0 ) e x + sin ( θ 0 ) e y ] = d [ cos ( θ θ 0 ) e r + sin ( θ θ 0 ) e θ ] {\displaystyle \mathbf {d} =d=d}

It remains

ψ ( r , θ ) = Q d 2 π sin ( θ θ 0 ) r {\displaystyle \psi (r,\theta )=-{\frac {Qd}{2\pi }}{\frac {\sin(\theta -\theta _{0})}{r}}}

The velocity is then

v r ( r , θ ) = Q d 2 π cos ( θ θ 0 ) r 2 {\displaystyle v_{r}(r,\theta )={\frac {Qd}{2\pi }}{\frac {\cos(\theta -\theta _{0})}{r^{2}}}}
v θ ( r , θ ) = Q d 2 π sin ( θ θ 0 ) r 2 {\displaystyle v_{\theta }(r,\theta )={\frac {Qd}{2\pi }}{\frac {\sin(\theta -\theta _{0})}{r^{2}}}}

And the potential instead

ϕ ( r , θ ) = Q d 2 π cos ( θ θ 0 ) r {\displaystyle \phi (r,\theta )={\frac {Qd}{2\pi }}{\frac {\cos(\theta -\theta _{0})}{r}}}

Two-dimensional vortex line

Potential flow streamlines for an ideal vortex line

This is the case of a vortex filament rotating at constant speed, there is a cylindrical symmetry and the problem can be solved in the orthogonal plane.

Dual to the case above of line sources, vortex lines play the role of monopoles for irrotational flows.

Also in this case the flow is also both irrotational and incompressible and therefore a case of potential flow.

This is characterized by a cylindrical symmetry:

v = v θ ( r ) e θ {\displaystyle \mathbf {v} =v_{\theta }(r)\,\mathbf {e} _{\theta }}

Where the total circulation is constant for every closed line around the central vortex

v d s = 0 2 π ( v θ ( r ) e θ ) ( e θ r d θ ) = 2 π r v θ ( r ) = Γ {\displaystyle \oint \mathbf {v} \cdot d\mathbf {s} =\int _{0}^{2\pi }(v_{\theta }(r)\,\mathbf {e} _{\theta })\cdot (\mathbf {e} _{\theta }\,r\,d\theta )=\!2\pi \,r\,v_{\theta }(r)=\Gamma }

and is zero for any line not including the vortex.

Therefore,

v θ = Γ 2 π r {\displaystyle v_{\theta }={\frac {\Gamma }{2\pi r}}}

This is derived from a stream function

ψ ( r , θ ) = Γ 2 π ln r {\displaystyle \psi (r,\theta )={\frac {\Gamma }{2\pi }}\ln r}

or from a potential function

ϕ ( r , θ ) = Γ 2 π θ {\displaystyle \phi (r,\theta )=-{\frac {\Gamma }{2\pi }}\theta }

Which is dual to the previous case of a line source

Generic two-dimensional potential flow

Given an incompressible two-dimensional flow which is also irrotational we have:

2 ψ = 0 {\displaystyle \nabla ^{2}\psi =0}

Which is in cylindrical coordinates

1 r r ( r ψ r ) + 1 r 2 2 ψ θ 2 = 0 {\displaystyle {\frac {1}{r}}{\frac {\partial }{\partial r}}\left(r{\frac {\partial \psi }{\partial r}}\right)+{\frac {1}{r^{2}}}{\frac {\partial ^{2}\psi }{\partial \theta ^{2}}}=0}

We look for a solution with separated variables:

ψ ( r , θ ) = R ( r ) Θ ( θ ) {\displaystyle \psi (r,\theta )=R(r)\Theta (\theta )}

which gives

r R ( r ) d d r ( r d R ( r ) d r ) = 1 Θ ( θ ) d 2 Θ ( θ ) d θ 2 {\displaystyle {\frac {r}{R(r)}}{\frac {d}{dr}}\left(r{\frac {dR(r)}{dr}}\right)=-{\frac {1}{\Theta (\theta )}}{\frac {d^{2}\Theta (\theta )}{d\theta ^{2}}}}

Given the left part depends only on r and the right parts depends only on θ {\displaystyle \theta } , the two parts must be equal to a constant independent from r and θ {\displaystyle \theta } . The constant shall be positive. Therefore,

r d d r ( r d d r R ( r ) ) = m 2 R ( r ) {\displaystyle r{\frac {d}{dr}}\left(r{\frac {d}{dr}}R(r)\right)=m^{2}R(r)}
d 2 Θ ( θ ) d θ 2 = m 2 Θ ( θ ) {\displaystyle {\frac {d^{2}\Theta (\theta )}{d\theta ^{2}}}=-m^{2}\Theta (\theta )}

The solution to the second equation is a linear combination of e i m θ {\displaystyle e^{im\theta }} and e i m θ {\displaystyle e^{-im\theta }} In order to have a single-valued velocity (and also a single-valued stream function) m shall be a positive integer.

therefore the most generic solution is given by

ψ = α 0 + β 0 ln r + m > 0 ( α m r m + β m r m ) sin [ m ( θ θ m ) ] {\displaystyle \psi =\alpha _{0}+\beta _{0}\ln r+\sum _{m>0}{\left(\alpha _{m}r^{m}+\beta _{m}r^{-m}\right)\sin {}}}

The potential is instead given by

ϕ = α 0 β 0 θ + m > 0 ( α m r m β m r m ) cos [ m ( θ θ m ) ] {\displaystyle \phi =\alpha _{0}-\beta _{0}\theta +\sum _{m\mathop {>} 0}{(\alpha _{m}r^{m}-\beta _{m}r^{-m})\cos {}}}

References

Specific
  1. Oliver, David (2013-03-14). The Shaggy Steed of Physics: Mathematical Beauty in the Physical World. Springer Science & Business Media. ISBN 978-1-4757-4347-0.
  2. Laplace operator

Further reading

External links

Category:
Elementary flow Add topic