Physics:Vortex sheet

From HandWiki

A vortex sheet is a term used in fluid mechanics for a surface across which there is a discontinuity in fluid velocity, such as in slippage of one layer of fluid over another.[1] While the tangential components of the flow velocity are discontinuous across the vortex sheet, the normal component of the flow velocity is continuous. The discontinuity in the tangential velocity means the flow has infinite vorticity on a vortex sheet. At high Reynolds numbers, vortex sheets tend to be unstable. In particular, they may exhibit Kelvin–Helmholtz instability.

The formulation of the vortex sheet equation of motion is given in terms of a complex coordinate [math]\displaystyle{ z = x + iy }[/math]. The sheet is described parametrically by [math]\displaystyle{ z(s,t) }[/math] where [math]\displaystyle{ s }[/math] is the arclength between coordinate [math]\displaystyle{ z }[/math] and a reference point, and [math]\displaystyle{ t }[/math] is time. Let [math]\displaystyle{ \gamma(s,t) }[/math] denote the strength of the sheet, that is, the jump in the tangential discontinuity. Then the velocity field induced by the sheet is

[math]\displaystyle{ \frac{\partial z^*}{\partial t} = -\frac{\imath}{2\pi} \int\limits_{-\infty}^{\infty} \frac{\gamma(s',t)\mathrm{d}s'}{z(s, t)-z(s', t)} }[/math]

The integral in the above equation is a Cauchy principal value integral. We now define [math]\displaystyle{ \Gamma }[/math] as the integrated sheet strength or circulation between a point with arc length [math]\displaystyle{ s }[/math] and the reference material point [math]\displaystyle{ s=0 }[/math] in the sheet.

[math]\displaystyle{ \Gamma(s,t) = \int\limits_{0}^{s}\gamma(s',t)\mathrm{d}s' \qquad\mathrm{and}\qquad \frac{\mathrm{d}\Gamma}{\mathrm{d}s}=\gamma(s,t) }[/math]

As a consequence of Kelvin's circulation theorem, in the absence of external forces on the sheet, the circulation between any two material points in the sheet remains conserved, so [math]\displaystyle{ \mathrm{d}\Gamma/\mathrm{d}t=0 }[/math]. The equation of motion of the sheet can be rewritten in terms of [math]\displaystyle{ \Gamma }[/math] and [math]\displaystyle{ t }[/math] by a change of variable. The parameter [math]\displaystyle{ s }[/math] is replaced by [math]\displaystyle{ \Gamma }[/math]. That is,

[math]\displaystyle{ \frac{\partial z^*}{\partial t}=-\frac{\imath}{2\pi} \int\limits_{-\infty}^{\infty}\frac{d\Gamma'}{z(\Gamma, t)-z(\Gamma', t)} }[/math]

This nonlinear integro-differential equation is called the Birkoff-Rott equation. It describes the evolution of the vortex sheet given initial conditions. Greater details on vortex sheets can be found in the textbook by Saffman (1977).

Diffusion of a vortex sheet

Once a vortex sheet, it will diffuse due to viscous action. Consider a planar unidirectional flow at [math]\displaystyle{ t=0 }[/math],

[math]\displaystyle{ u= \begin{cases} +U, & \text{for } y\gt 0\\ -U, & \text{for } y\lt 0 \end{cases} }[/math]

impling the presence of a vortex sheet at [math]\displaystyle{ y=0 }[/math]. The velocity discontinuity smooths out according to[2]

[math]\displaystyle{ u(y,t)=\frac{U}{2\sqrt{\pi\nu t}}\left[\int_0^\infty e^{-(s-y)^2/4\nu t}\mathrm{d}s - \int_0^\infty e^{-(s+y)^2/4\nu t}\mathrm{d}s\right]. }[/math]

where [math]\displaystyle{ \nu }[/math] is the kinematic viscosity. The only non-zero vorticity component is in the [math]\displaystyle{ z }[/math] direction, given by

[math]\displaystyle{ \omega_z = - \frac{U}{\sqrt{\pi \nu t}} e^{-y^2/4\nu t} }[/math].

Vortex sheet with periodic boundaries

A flat vortex sheet with periodic boundaries in the streamwise direction can be used to model a temporal free shear layer at high Reynolds number. Let us assume that the interval between the periodic boundaries is of length [math]\displaystyle{ 1 }[/math]. Then the equation of motion of the vortex sheet reduces to

[math]\displaystyle{ \frac{\partial z^*}{\partial t}=-\frac{\imath}{2} \int\limits_{0}^{1}\cot \pi(z(\Gamma, t)-z(\Gamma', t))\;d\Gamma' }[/math]

Continuous vortex sheet approximation by panel method. Roll-up of a vortex sheet due to an initial sinusoidal perturbation.

Note that the integral in the above equation is a Cauchy principal value integral. The initial condition for a flat vortex sheet with constant strength is [math]\displaystyle{ z(\Gamma, 0) = \Gamma }[/math]. The flat vortex sheet is an equilibrium solution. However, it is unstable to infinitesimal periodic disturbances of the form [math]\displaystyle{ \sum^{\infty}_{k=-\infty}A_k\mathrm{e}^{\imath2\pi k\Gamma} }[/math]. Linear theory shows that the Fourier coefficient [math]\displaystyle{ A_k }[/math] grows exponentially at a rate proportional to [math]\displaystyle{ k }[/math]. That is, higher the wavenumber of a Fourier mode, the faster it grows. However, a linear theory cannot be extended much beyond the initial state. If nonlinear interactions are taken into account, asymptotic analysis suggests that for large [math]\displaystyle{ k }[/math] and finite [math]\displaystyle{ t\lt t_c }[/math], where [math]\displaystyle{ t_c }[/math] is a critical value, the Fourier coefficient [math]\displaystyle{ A_k }[/math] decays exponentially. The vortex sheet solution is expected to lose analyticity at the critical time. See Moore (1979), and Meiron, Baker and Orszag (1983).

The vortex sheet solution as given by the Birkoff-Rott equation cannot go beyond the critical time. The spontaneous loss of analyticity in a vortex sheet is a consequence of mathematical modeling since a real fluid with viscosity, however small, will never develop singularity. Viscosity acts a smoothing or regularization parameter in a real fluid. There have been extensive studies on a vortex sheet, most of them by discrete or point vortex approximation, with or without desingularization. Using a point vortex approximation and delta-regularization Krasny (1986) obtained a smooth roll-up of a vortex sheet into a double branched spiral. Since point vortices are inherently chaotic, a Fourier filter is necessary to control the growth of round-off errors. Continuous approximation of a vortex sheet by vortex panels with arc wise diffusion of circulation density also shows that the sheet rolls-up into a double branched spiral.

In many engineering and physical applications the growth of a temporal free shear layer is of interest. The thickness of a free shear layer is usually measured by momentum thickness, which is defined as

[math]\displaystyle{ \theta = \int\limits_{y=-\infty}^{\infty}\left(\frac{1}{4} - \left(\frac{\left\langle u\right\rangle}{2U}\right)^2\right)\mathrm{d}y }[/math]

where [math]\displaystyle{ \left\langle u\right\rangle = \frac{1}{L}\int_{0}^{L}\ u(x,y,t) dx }[/math] and [math]\displaystyle{ U }[/math] is the freestream velocity. Momentum thickness has the dimension of length and the non-dimensional momentum thickness is given by [math]\displaystyle{ \theta_{ND} = \theta/L }[/math]. Momentum thickness can be used to measure the thickness of a vortex layer.

See also

References

  1. McGraw-Hill Dictionary of Scientific and Technical Terms Retrieved July 2012
  2. Drazin, P. G., & Riley, N. (2006). The Navier-Stokes equations: a classification of flows and exact solutions (No. 334). Cambridge University Press.