Weakly nonlinear stability analysis of magnetohydrodynamic channel flow using an efficient numerical approach

We analyze weakly nonlinear stability of a flow of viscous conducting liquid driven by pressure gradient in the channel between two parallel walls subject to a transverse magnetic field. Using a non-standard numerical approach, we compute the linear growth rate correction and the first Landau coefficient, which in a sufficiently strong magnetic field vary with the Hartmann number as \documentclass[12pt]{minimal}\begin{document}$\mu _{1}\sim (0.814-\mathrm{i}19.8)\times 10^{-3}\textit {Ha}$\end{document}μ1∼(0.814−i19.8)×10−3Ha and \documentclass[12pt]{minimal}\begin{document}$\mu _{2}\sim (2.73-\mathrm{i}1.50)\times 10^{-5}\textit {Ha}^{-4}.$\end{document}μ2∼(2.73−i1.50)×10−5Ha−4. These coefficients describe a subcritical transverse velocity perturbation with the equilibrium amplitude \documentclass[12pt]{minimal}\begin{document}$|A|^{2}=\Re [\mu _{1}]/\Re [\mu _{2}](\textit {Re}_{c}-\textit {Re})\sim 29.8\textit {Ha}^{5}(\textit {Re}_{c}-\textit {Re}),$\end{document}|A|2=ℜ[μ1]/ℜ[μ2](Rec−Re)∼29.8Ha5(Rec−Re), which exists at Reynolds numbers below the linear stability threshold \documentclass[12pt]{minimal}\begin{document}$\textit {Re}_{c}\sim 4.83\times 10^{4}\textit {Ha}.$\end{document}Rec∼4.83×104Ha. We find that the flow remains subcritically unstable regardless of the magnetic field strength. Our method for computing Landau coefficients differs from the standard one by the application of the solvability condition to the discretized rather than continuous problem. This allows us to bypass both the solution of the adjoint problem and the subsequent evaluation of the integrals defining the inner products, which results in a significant simplification of the method.