Michelson–Sivashinsky equation

In combustion, Michelson–Sivashinsky equation describes the evolution of a premixed flame front, subjected to the Darrieus–Landau instability, in the small heat release approximation. The equation was derived by Gregory Sivashinsky in 1977, who along the Daniel M. Michelson, presented the numerical solutions of the equation in the same year. Let the planar flame front, in a uitable frame of reference be on the $xy$-plane, then the evolution of this planar front is described by the amplitude function $u(\mathbf {x} ,t)$ (where $\mathbf {x} =(x,y)$) describing the deviation from the planar shape. The Michelson–Sivashinsky equation, reads as

${\frac {\partial u}{\partial t}}+{\frac {1}{2}}(\nabla u)^{2}-\nu \nabla ^{2}u-{\frac {1}{8\pi ^{2}}}\int |\mathbf {k} |e^{i\mathbf {k} \cdot (\mathbf {x} -\mathbf {x} ')}u(\mathbf {x} ,t)d\mathbf {k} d\mathbf {x} '=0,$

where $\nu$ is a constant. Incorporating also the Rayleigh–Taylor instability of the flame, one obtains the Rakib–Sivashinsky equation (named after Z. Rakib and Gregory Sivashinsky),

${\frac {\partial u}{\partial t}}+{\frac {1}{2}}(\nabla u)^{2}-\nu \nabla ^{2}u-{\frac {1}{8\pi ^{2}}}\int |\mathbf {k} |e^{i\mathbf {k} \cdot (\mathbf {x} -\mathbf {x} ')}u(\mathbf {x} ,t)d\mathbf {k} d\mathbf {x} '+\gamma \left(u-\langle u\rangle \right)=0,\quad$

where $\langle u\rangle (t)$ denotes the spatial average of $u$, which is a time-dependent function and $\gamma$ is another constant.

N-pole solution

The equations, in the absence of gravity, admits an explicit solution, which is called as the N-pole solution since the equation admits a pole decomposition,as shown by Olivier Thual, Uriel Frisch and Michel Hénon in 1988. Consider the 1d equation

$u_{t}+uu_{x}-\nu u_{xx}=\int _{-\infty }^{+\infty }e^{ikx}{\hat {u}}(k,t)dk,$

where ${\hat {u}}$ is the Fourier transform of $u$. This has a solution of the form

${\begin{aligned}u(x,t)&=-2\nu \sum _{n=1}^{2N}{\frac {1}{x-z_{n}(t)}},\\{\frac {dz_{n}}{dt}}&=-2\nu \sum _{l=1,l\neq n}^{2N}{\frac {1}{z_{n}-z_{l}}}-i\mathrm {sgn} (\mathrm {Im} z_{n}),\end{aligned}}$

where $z_{n}(t)$ (which appear in complex conjugate pairs) are poles in the complex plane. In the case periodic solution with periodicity $2\pi$, the it is sufficient to consider poles whose real parts lie between the interval $0$ and $2\pi$. In this case, we have

${\begin{aligned}u(x,t)&=-\nu \sum _{n=1}^{2\pi }\cot {\frac {x-z_{n}(t)}{2}},\\{\frac {dz_{n}}{dt}}&=-\nu \sum _{l\neq n}\cot {\frac {z_{n}-z_{l}}{2}}-i\mathrm {sgn} (\mathrm {Im} z_{n})\end{aligned}}$

These poles are interesting because in physical space, they correspond to locations of the cusps forming in the flame front.

Dold–Joulin equation

In 1995, John W. Dold and Guy Joulin generalised the Michelson–Sivashinsky equation by introducing the second-order time derivative, which is consistent with the quadratic nature of the dispersion relation for the Darrieus–Landau instability. The Dold–Joulin equation is given by

${\frac {\partial ^{2}\varphi }{\partial t^{2}}}+{\mathcal {I}}\left({\frac {\partial \varphi }{\partial t}}-{\frac {1}{2}}(\nabla \varphi )^{2}-\nu \nabla ^{2}\varphi -\nu {\mathcal {I}}(\varphi )\right)=0,$

where ${\mathcal {I}}(e^{i\mathbf {k} \cdot \mathbf {x} })=|\mathbf {k} |e^{i\mathbf {k} \cdot \mathbf {x} }$ corresponds to the non-local integral operator.

See also

• Kuramoto–Sivashinsky equation