Physics:Rayleigh–Kuo criterion

Short description: Stability condition for fluids

Figure 1: The induced horizontal and vertical circulation patterns induced by a horizontal shear flow at the surface of the fluid. The horizontal shear is indicated by the length and colour density of the arrows at the surface of the fluid. This shear in the velocity at the surface leads to both vertical and horizontal circulation patterns. Due to Ekman transport from the right side of the flow towards the left, there is divergence (convergence) on the right (left) side of the flow which leads to upwelling (downwelling) and together with the horizontal return flow, this is a full circulation pattern. In the horizontal there can be formation of eddies on either side of the flow which have a circular rotation pattern clockwise on the right side of the flow (blue) and anti-clockwise on the left side of the flow (red). The eddies form on alternating sides of the flow and move with the flow.

The Rayleigh–Kuo criterion (sometimes called the Kuo criterion) is a stability condition for a fluid. This criterion determines whether or not a barotropic instability can occur, leading to the presence of vortices (like eddies and storms). The Kuo criterion states that for barotropic instability to occur, the gradient of the absolute vorticity must change its sign at some point within the boundaries of the current. Note that this criterion is a necessary condition, so if it does not hold it is not possible for a barotropic instability to form. But it is not a sufficient condition, meaning that if the criterion is met, this does not automatically mean that the fluid is unstable. If the criterion is not met, it is certain that the flow is stable.

This criterion was formulated by Hsiao-Lan Kuo and is based on Rayleigh's equation named after the Lord Rayleigh who first introduced this equation in fluid dynamics.

Barotropic instability

Figure 2: The shear stress in the flow of the fluid induces eddies. The upper left panel (a), shows the differences in speed with the length of the black arrows and the density of the colours (darker colours mean larger velocity). When this flow becomes unstable (b), there direction of the flow begins to change. This process is further enhanced in panel (c) until finally full eddies occur (d). These eddies form on alternating sides of the flow field with a clockwise motion on the right (blue circles) and an anti-clockwise circulation on the left side of the flow (red circles).

Vortices like eddies are created by instabilities in a flow. When there are instabilities within the mean flow, energy can be transferred from the mean flow to the small perturbations which can then grow. In a barotropic fluid the density is a function of only the pressure and not the temperature (in contrast to a baroclinic fluid, where the density is a function of both the pressure and temperature). This means that surfaces of constant density (isopycnals) are also surfaces of constant pressure (isobars). Barotropic instability can form in different ways. Two examples are; when there is an interaction between the fluid flow and the bathymetry or topography of the domain; when there are frontal instabilities (may also lead to baroclinic instabilities). These instabilities are not dependent on the density and might even occur when the density of the fluid is constant. Instead, most of the instabilities are caused by a shear on the flow as can be seen in Figure 1. This shear in the velocity field induces a vertical and horizontal vorticity within the flow. As a result, there is upwelling on the right of the flow and downwelling on the left. This situation might lead to a barotropic unstable flow. The eddies that form alternatingly on both sides of the flow are part of this instability.

Another way to achieve this instability is to displace the Rossby waves in the horizontal direction (see Figure 2). This leads to a transfer of kinetic energy (not potential energy) from the mean flow towards the small perturbations (the eddies). The Rayleigh–Kuo criterion states that the gradient of the absolute vorticity should change sign within the domain. In the example of the shear induced eddies on the right, this means that the second derivative of the flow in the cross-flow direction, should be zero somewhere. This happens in the centre of the eddies, where the acceleration of the flow perpendicular to the flow changes direction.

Examples

The presence of these instabilities in a rotating fluid have been observed in laboratory experiments. The settings of the experiment were based on the conditions in the Gulf Stream and showed that within the ocean currents such as the Gulf Stream, it is possible for barotropic instabilities to occur. But barotropic instabilities were also observed in other Western Boundary Currents (WBC). In the Agulhas current, the barotropic instability leads to ring shedding. The Agulhas current retroflects (turns back) near the coast of South Africa. At this same location, some anti-cyclonic rings of warm water escape from the mean current and travel along the coast of Africa. The formation of these rings is a manifestation of a barotropic instability.

Derivation

The derivation of the Rayleigh–Kuo criterion was first written down by Hsiao-Lan Kuo in his paper called 'dynamic instability of two-dimensional nondivergent flow in a barotropic atmosphere' from 1949. This derivation is repeated and simplified below.

First, the assumptions made by Hsiao-Lan Kuo are discussed. Second, the Rayleigh equation is derived in order continue to derive the Rayleigh–Kuo criterion. By integrating this equation and filling in the boundary conditions, the Kuo criterion can be obtained.

Assumptions

In order to derive the Rayleigh–Kuo criterion, some assumptions are made on the fluids properties. We consider a nondivergent, two-dimensional barotropic fluid. The fluid has a mean zonal flow direction which can vary in the meridional direction. On this mean flow, some small perturbations are imposed in both the zonal and meridional direction: [math]\displaystyle{ u(y,t) = U(y) + u^*(y,t) }[/math] and [math]\displaystyle{ v = v^* }[/math]. The perturbations need to be small in order to linearize the vorticity equation. Vertical motion and divergence and convergence of the fluid are neglected. When taking into account these factors, a similar result would have been obtained with only a small shift in the position of the criterion within the velocity profile.

The derivation of the Kuo criterion will be done within the domain [math]\displaystyle{ L=[0,y ] }[/math]. On the northern and southern boundary of this domain, the meridional fluid is zero.

Rayleigh Equation

Barotropic vorticity equation

To derive the Rayleigh equation for a barotropic fluid, the barotropic vorticity equation is used. This equation assumes that the absolute vorticity is conserved: [math]\displaystyle{ \frac{d \zeta_a}{dt} = 0 }[/math] here, [math]\displaystyle{ \frac{d}{dt} }[/math] is the material derivative. The absolute vorticity is the relative vorticity plus the planetary vorticity: [math]\displaystyle{ \zeta_a = \zeta + f }[/math]. The relative vorticity, [math]\displaystyle{ \zeta }[/math], is the rotation of the fluid with respect to the Earth. The planetary vorticity (also called Coriolis frequency),[math]\displaystyle{ f }[/math], is the vorticity of a parcel induced by the rotation of the Earth. When applying the beta-plane approximation for the planetary vorticity, the conservation of absolute vorticity looks like:

[math]\displaystyle{ \frac{d\zeta_a}{dt} = \frac{d}{dt}\left(\zeta + \beta y \right) = 0 }[/math] The relative vorticity is defined as [math]\displaystyle{ \zeta = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}. }[/math] Since the flow field consist of a mean flow with small perturbations, it can be written as [math]\displaystyle{ \zeta = \overline{\zeta} + \zeta^* }[/math] with [math]\displaystyle{ \overline{\zeta} = -\frac{\partial U}{\partial y} }[/math] and [math]\displaystyle{ \zeta^* = \frac{\partial v^*}{\partial x} - \frac{\partial u^*}{\partial y}. }[/math] This formulation is used in the vorticity equation:

[math]\displaystyle{ \begin{align} 0 &= \frac{d}{dt}\left(\zeta + \beta y\right) \\ 0 &= \frac{d}{dt}\left(\zeta ' + \overline{\zeta} + \beta y\right) \\ 0 &= \left( \frac{\partial}{\partial t} + u\frac{\partial}{\partial x} + v\frac{\partial}{\partial y} \right) \left(\zeta^* - \frac{\partial U}{\partial y} + \beta y \right) \end{align} }[/math]Here, [math]\displaystyle{ u }[/math] and [math]\displaystyle{ v }[/math] are the zonal and meridional components of the flow and [math]\displaystyle{ \zeta' }[/math] is the relative vorticity induced by the perturbations on the flow ([math]\displaystyle{ u' }[/math] and [math]\displaystyle{ v' }[/math]). [math]\displaystyle{ U }[/math] is the mean zonal flow and [math]\displaystyle{ \beta }[/math] is derivative of the planetary vorticity [math]\displaystyle{ f }[/math] with respect to [math]\displaystyle{ y }[/math] .

Linearization

A zonal mean flow with small perturbations was assumed, [math]\displaystyle{ u = U+ u^* }[/math], and a meridional flow with a zero mean, [math]\displaystyle{ v = v^* }[/math]. Since it was assumed that the perturbations are small, a linearization can be performed on the barotropic vorticity equation above, ignoring all the non-linear terms (terms where two or more small variables, i.e. [math]\displaystyle{ u^*, v^*, \zeta^* }[/math], are multiplied with one another). Also the derivative of [math]\displaystyle{ u }[/math] in the zonal direction, the time derivative of the mean flow [math]\displaystyle{ U }[/math] and the time derivative of [math]\displaystyle{ \beta y }[/math] are zero. This results in a simplified equation:

[math]\displaystyle{ \begin{align} 0 &= \left( \frac{\partial }{\partial t} + U\frac{\partial}{\partial x}\right)\zeta' - v'\frac{\partial}{\partial y}\frac{\partial U}{\partial y} + v'\frac{\partial}{\partial y} \left(\beta y\right)\\ 0 &= \left( \frac{\partial }{\partial t} + U\frac{\partial}{\partial x} \right)\zeta' + v' \left( \beta - \frac{\partial^2 U}{\partial y^2}\right). \end{align} }[/math]

With [math]\displaystyle{ \zeta' }[/math] as defined above ([math]\displaystyle{ \zeta^* = \frac{\partial v^*}{\partial x} - \frac{\partial u^*}{\partial y} }[/math]) and [math]\displaystyle{ u^* }[/math] and [math]\displaystyle{ v^* }[/math] the small perturbations in the zonal and meridional components of the flow.

Stream function

To find the solution to the linearized equation, a stream function was introduced by Lord Rayleigh for the perturbations of the flow velocity:

[math]\displaystyle{ u^* = \frac{\partial \psi}{\partial y}, \;\;\;\; v^* = \frac{-\partial \psi}{\partial x}. }[/math]These new definitions of the stream function are used to rewrite the linearized barotropic vorticity equation. [math]\displaystyle{ \begin{align} 0 &= \left(\frac{\partial}{\partial t} + U\frac{\partial}{\partial x}\right) \left(\frac{\partial v^*}{\partial x} - \frac{\partial u^*}{\partial y}\right) + v^*\left(\beta - \frac{\partial^2 U}{\partial y^2}\right)\\ 0 &= \left(\frac{\partial}{\partial t} + U\frac{\partial}{\partial x}\right) \left(-\frac{\partial^2 \psi}{\partial x^2} - \frac{\partial^2 \psi}{\partial y ^2}\right) - \frac{\partial \psi}{\partial x}\left(\beta - U''\right)\\ 0 &= \left(\frac{\partial}{\partial t} + U\frac{\partial}{\partial x}\right) \nabla^2 \psi + \frac{\partial \psi}{\partial x}\left(\beta - U''\right) \end{align} }[/math]Here, [math]\displaystyle{ U'' }[/math] is the second derivative of [math]\displaystyle{ U }[/math] with respect to [math]\displaystyle{ y }[/math] [math]\displaystyle{ (U'' = \frac{\partial^2 U}{\partial y^2}) }[/math]. To solve this equation for the stream function, a wave-like solution was proposed by Rayleigh which reads [math]\displaystyle{ \psi(x,y,t) = \Psi(y)e^{i\alpha\left(x-ct\right)} }[/math]. The amplitude [math]\displaystyle{ \Psi(y) }[/math] may be complex number, [math]\displaystyle{ \alpha }[/math] is the wave number which is a real number and [math]\displaystyle{ c }[/math] is the phase velocity which may be complex as well. Inserting this proposed solution leads us to the equation which is known as Rayleigh's equation.

[math]\displaystyle{ \begin{align} 0 &= \left(\frac{\partial}{\partial t} + U\frac{\partial}{\partial x}\right) \nabla^2 \psi + \frac{\partial \psi}{\partial x}\left(\beta - U''\right)\\[14pt] 0 &= \left(\frac{\partial}{\partial t} + U\frac{\partial}{\partial x}\right)\left((i\alpha)^2\Psi e^{i\alpha\left(x-ct \right)} + \Psi''e^{i\alpha\left(x-ct\right)}\right) + i\alpha\Psi e^{i\alpha(x-ct)}(\beta - U'')\\[14pt] 0 &= (i\alpha)^2\Psi e^{i\alpha(x-ct)}(-i\alpha c + i\alpha U ) + \Psi''e^{i\alpha(x-ct)} (-i\alpha c+i\alpha U) + i\alpha\Psi e^{i\alpha(x-ct)}(\beta - U'')\\[14pt] 0 &= i\alpha e^{i\alpha(x-ct)}[(\alpha^2c - \alpha^2 U)\Psi + \Psi''(-c + U) + \Psi(\beta - U'')]\\[14pt] 0 &= (U-c)(\Psi'' - \alpha^2\Psi) + (\beta - U'')\Psi)\\[14pt] \end{align} }[/math]To get to this equation, in the last step it was used that [math]\displaystyle{ \alpha }[/math] can't be zero and neither can the exponential. This means that the terms in the square brackets needs to be zero. The symbol [math]\displaystyle{ \Psi'' }[/math] denotes the second derivative of the amplitude of the stream function, [math]\displaystyle{ \Psi }[/math] with respect to [math]\displaystyle{ y }[/math] [math]\displaystyle{ (\Psi'' = \frac{\partial ^2 \Psi}{\partial y^2}) }[/math]. This last equation that was derived, is known as Rayleigh's equation which is a linear ordinary differential equation. It is very difficult to explicitly solve this equation. It is therefore that Hsiao-Lan Kuo came up with a stability criterion for this problem without actually solving it.

Kuo Criterion

Instead of solving Rayleigh's equation, Hsiao-Lan Kuo came up with a necessary stability condition which had to be met in order for the fluid to be able to get unstable. To get to this criterion, Rayleigh's equation was rewritten and the boundary conditions of the flow field are used.

The first step is to divide Rayleigh's equation by [math]\displaystyle{ (U-c) }[/math] and multiplying the equation by the complex conjugate of [math]\displaystyle{ \Psi\;\; (\Psi^*= \Psi_r - i\Psi_i) }[/math].

[math]\displaystyle{ \begin{align} 0 &= (\Psi'' - \alpha^2\Psi) + \left( \frac{\beta - U''}{U-c} \right)\Psi\\ 0 &= \Psi^*(\Psi_r'' + i\Psi_i'' - \alpha^2\Psi_r - \alpha^2 i\Psi_i) + \Psi^*\left(\frac{\beta - U''}{U-c}\right)(\Psi_r - i\Psi_i)\\ 0 &= \Psi_r\Psi_r'' + \Psi_i\Psi_i'' - \alpha^2(\Psi_r^2 - \Psi_i^2) + \left(\frac{\beta - U''}{U-c}\right) (\Psi_r^2 - \Psi_i^2) + i(-\Psi_i\Psi_r'' + \Psi_r\Psi_i'')\\ 0 &= \Psi_r''\Psi_r + \Psi_i''\Psi_i + \left(-\alpha^2 + \frac{U - c_r}{|U-c|^2}(\beta - U'')\right)|\Psi|^2 + i\left(\frac{c_i}{|U-c|^2}(\beta - U'')|\Psi|^2 - \Psi_r''\Psi_i + \Psi_i''\Psi_r\right) \end{align} }[/math]In the last step, [math]\displaystyle{ (U-c) }[/math] is multiplied with its complex conjugate leading to the following equality is used: [math]\displaystyle{ \frac{1}{U-c} =\frac{1}{U-c_r - ic_i} = \frac{U-c_r + ic_i}{|U-c|^2} }[/math]. For the solution of Rayleigh's equation to exist, both the real and imaginary part of the equation above need to be equal to zero.

Boundary conditions

To get to the Kuo criterion, the imaginary part is integrated over the domain ([math]\displaystyle{ y=[0,L] }[/math]) . The stream function at the boundaries of the domain is zero, [math]\displaystyle{ \Psi(0) = \Psi(L) = 0 }[/math], as already stated in the assumptions. The zonal flow must vanish at the boundaries of the domain. This leads to a constant stream function which is set to zero for convenience.

[math]\displaystyle{ \int_0^L(\Psi_r\Psi_i'' - \Psi_i\Psi_r'')dy + \int_0^L\left(c_i \frac{|\Psi|^2}{|U-c|}(\beta - U'')\right) =0 }[/math]

The first integral can be solved:

[math]\displaystyle{ \begin{align} \int_0^L (\Psi_r\Psi_i'' - \Psi_i\Psi_r'')dy &= \int_0^L \frac{\partial }{\partial y}(\Psi_r\Psi_i' - \Psi_i\Psi_r')dy\\[8pt] &= (\Psi_r\Psi_i' - \Psi_i\Psi_r')|_0^L \\[8pt] &= 0 \\[12pt] \end{align} }[/math]

So the first integral is equal to zero. This means that the second integral should also be zero, making it possible to solve this integral numerically.

[math]\displaystyle{ \begin{align} \int_0^L\left(c_i \frac{|\Psi|^2}{|U-c|}(\beta - U'')\right)dy &=0\\ \end{align} }[/math]

When [math]\displaystyle{ c_i }[/math] is zero, we are dealing with a stable amplitude of the solution, this means that the solution is stable. We are looking for un unstable situation, so then [math]\displaystyle{ \frac{|\Psi|^2}{|U-c|}(\beta - U'') }[/math] should be zero. Since the fraction in front of [math]\displaystyle{ (\beta - U'') }[/math] is non-zero and positive, this leads to the conclusion that [math]\displaystyle{ (\beta - U'') }[/math] should be zero. This leads to the final formulation, the Kuo criterion:

[math]\displaystyle{ \begin{align} \beta - U'' &= 0\\[10pt] \beta &= U''\\[10pt] \frac{\partial(\beta y)}{\partial y} &= \frac{\partial^2 U}{\partial y^2} \end{align} }[/math]Here, [math]\displaystyle{ U }[/math] is the mean zonal flow and [math]\displaystyle{ \beta }[/math] is the derivative of the planetary vorticity [math]\displaystyle{ f }[/math] with respect to [math]\displaystyle{ y }[/math].