Physics:Stuart–Landau equation

The Stuart–Landau equation describes the behavior of a nonlinear oscillating system near the Hopf bifurcation, named after John Trevor Stuart and Lev Landau. In 1944, Landau proposed an equation for the evolution of the magnitude of the disturbance, which is now called as the Landau equation, to explain the transition to turbulence based on a phenomenological argument and an attempt to derive this equation from hydrodynamic equations was done by Stuart for plane Poiseuille flow in 1958. The formal derivation to derive the Landau equation was given by Stuart, Watson and Palm in 1960. The perturbation in the vicinity of bifurcation is governed by the following equation

[math]\displaystyle{ \frac{dA}{dt} = \sigma A - \frac{l}{2} A |A|^2. }[/math]

where

• [math]\displaystyle{ A = |A| e^{i\phi} }[/math] is a complex quantity describing the disturbance,
• [math]\displaystyle{ \sigma = \sigma_r + i\sigma_i }[/math] is the complex growth rate,
• [math]\displaystyle{ l = l_r + i l_i }[/math] is a complex number and [math]\displaystyle{ l_r }[/math] is the Landau constant.

The evolution of the actual disturbance is given by the real part of [math]\displaystyle{ A(t) }[/math] i.e., by [math]\displaystyle{ |A|\cos\phi }[/math]. Here the real part of the growth rate is taken to be positive, i.e., [math]\displaystyle{ \sigma_r\gt 0 }[/math] because otherwise the system is stable in the linear sense, that is to say, for infinitesimal disturbances ([math]\displaystyle{ |A| }[/math] is a small number), the nonlinear term in the above equation is negligible in comparison to the other two terms in which case the amplitude grows in time only if [math]\displaystyle{ \sigma_r\gt 0 }[/math]. The Landau constant is also taken to be positive, [math]\displaystyle{ l_r\gt 0 }[/math] because otherwise the amplitude will grow indefinitely (see below equations and the general solution in the next section). The Landau equation is the equation for the magnitude of the disturbance,

[math]\displaystyle{ \frac{d|A|^2}{dt} = 2\sigma_r |A|^2 - l_r |A|^4, }[/math]

which can also be re-written as

[math]\displaystyle{ \frac{d|A|}{dt} = \sigma_r |A| -\frac{l_r}{2} |A|^3. }[/math]

Similarly, the equation for the phase is given by

[math]\displaystyle{ \frac{d\phi}{dt}= \sigma_i-\frac{l_i}{2} |A|^2. }[/math]

Due to the universality of the equation, the equation finds its application in many fields such as hydrodynamic stability, Belousov–Zhabotinsky reaction, etc.

General solution

The Landau equation is linear when it is written for the dependent variable [math]\displaystyle{ |A|^{-2} }[/math],

[math]\displaystyle{ \frac{d|A|^{-2}}{dt} + 2\sigma_r |A|^{-2} = l_r. }[/math]

The general solution for [math]\displaystyle{ \sigma_r\neq 0 }[/math] of the above equation is

[math]\displaystyle{ |A(t)|^{-2} = \frac{l_r}{2\sigma_r} + \left(|A(0)|^{-2} - \frac{l_r}{2\sigma_r}\right)e^{-2\sigma_r t}. }[/math]

As [math]\displaystyle{ t\rightarrow\infty }[/math], the magnitude of the disturbance [math]\displaystyle{ |A| }[/math] approaches a constant value that is independent of its initial value, i.e., [math]\displaystyle{ |A|_{\mathrm{max}}\rightarrow(2\sigma_r/l_r)^{1/2} }[/math] when [math]\displaystyle{ t\gg 1/\sigma_r }[/math]. The above solution implies that [math]\displaystyle{ |A| }[/math] does not have a real solution if [math]\displaystyle{ l_r\lt 0 }[/math] and [math]\displaystyle{ \sigma_r\gt 0 }[/math]. The associated solution for the phase function [math]\displaystyle{ \phi(t) }[/math] is given by

[math]\displaystyle{ \phi(t)-\phi(0) = \sigma_i t - \frac{l_i}{2l_r} \ln \left[1+ \frac{|A(0)|^2l_r}{2\sigma_r}(e^{2\sigma_r t}-1)\right]. }[/math]

As [math]\displaystyle{ t\gg 1/\sigma_r }[/math], the phase varies linearly with time, [math]\displaystyle{ \phi \sim (\sigma_i/\sigma_r-l_i/l_r)\sigma_rt. }[/math]

It is instructive to consider a hydrodynamic stability case where it is found that, according to the linear stability analysis, the flow is stable when [math]\displaystyle{ Re\leq Re_{\mathrm{cr}} }[/math] and unstable otherwise, where [math]\displaystyle{ Re }[/math] is the Reynolds number and the [math]\displaystyle{ Re_{\mathrm{cr}} }[/math] is the critical Reynolds number; a familiar example that is applicable here is the critical Reynolds number, [math]\displaystyle{ Re_{\mathrm{cr}}\approx 50 }[/math], corresponding to the transition to Kármán vortex street in the problem of flow past a cylinder. The growth rate [math]\displaystyle{ \sigma_r }[/math] is negative when [math]\displaystyle{ Re\lt Re_{\mathrm{cr}} }[/math] and is positive when [math]\displaystyle{ Re\gt Re_{\mathrm{cr}} }[/math] and therefore in the neighbourhood [math]\displaystyle{ Re\rightarrow Re_{\mathrm{cr}} }[/math], it may written as [math]\displaystyle{ \sigma_r=\text{const}.\times (Re-Re_{\mathrm{cr}}) }[/math] wherein the constant is positive. Thus, the limiting amplitude is given by

[math]\displaystyle{ |A|_{\mathrm{max}} \propto \sqrt{Re-Re_{\mathrm{cr}}}. }[/math]

Negative Landau constant

When the Landau constant is negative, [math]\displaystyle{ l_r\lt 0 }[/math], we must include a negative term of higher order to arrest the unbounded increase of the perturbation. In this case, the Landau equation becomes

[math]\displaystyle{ \frac{d|A|^2}{dt} = 2\sigma_r |A|^2 - l_r |A|^4 - \beta_r|A|^6, \quad \beta_r\gt 0. }[/math]

The limiting amplitude then becomes

[math]\displaystyle{ |A|_{\mathrm{max}}\rightarrow \frac{|l_r|}{2\beta_r} \pm \sqrt{\frac{l_r^2}{4\beta_r^2}+\frac{2|l_r|\sigma_r}{\beta_r}}, \quad \text{as} \quad t\gg 1/\sigma_r }[/math]

where the plus sign corresponds to the stable branch and the minus sign to the unstable branch. There exists a value of a critical value [math]\displaystyle{ Re_{\mathrm{cr}}' }[/math] where the above two roots are equal ([math]\displaystyle{ \sigma_r = -|l_r|/8\beta_r }[/math]) such that [math]\displaystyle{ Re_{\mathrm{cr}}'\lt Re_{\mathrm{cr}} }[/math], indicating that the flow in the region [math]\displaystyle{ Re_{\mathrm{cr}}'\lt Re\lt Re_{\mathrm{cr}} }[/math] is metastable, that is to say, in the metastable region, the flow is stable to infinitesimal perturbations, but not to finite amplitude perturbations.

See also

• Landau's phase transition theory
• Ginzburg–Landau theory