Melting and viscous dissipation effects on MHD flow over a moving surface with constant heat source (Q2317126)

The authors consider a 2D magnetohydrodynamics (MHD) boundary layer flow of an incompressible and electrically conducting fluid towards a moving surface which melts at a steady state into a warm liquid. They write the system as \(\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0\), \(u \frac{\partial u}{\partial x}+v\frac{\partial v}{\partial y}=\vartheta \frac{ \partial ^{2}u}{\partial y^{2}}-\frac{\sigma B^{2}(u)}{\rho }(u-U_{\infty })\), \(u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{\partial ^{2}T}{\partial y^{2}}-\frac{1}{\rho C_{p}}\frac{\partial q_{r}}{\partial y}-\frac{Q_{H}}{\rho C_{p}}(T-T_{\infty })+\frac{\vartheta }{ \rho C_{p}}(\frac{\partial u}{\partial y})^{2}\), where \((u,v)\) is the velocity vector, \(B\) is the applied magnetic field, \(\sigma \) is the electrical conductivity of the fluid, \(\rho \) its density, \(T\) its temperature, \(\alpha \) its thermal conductivity and \(q_{r}\) the radiative heat flux. The boundary conditions are \(u=U_{w}\), \(T=T_{m}\) at \(y=0\), \( u\rightarrow U_{\infty }\), \(T\rightarrow T_{\infty }\) as \(y\rightarrow \infty \) and \(k(\frac{\partial T}{\partial y})_{y=0}=\rho \{\lambda +c_{s}(T_{m}-T_{\infty })\}v(x,0)\). They introduce the Rosseland approximation and the approximation \(T^{4}=4TT_{\infty }^{3}-3T_{\infty }^{4} \) to simplify the problem. They finally introduce a change of variables and of unknowns and the stream function \(\psi \) so that \(u=\frac{\partial \psi }{ \partial y}\) and \(v=-\frac{\partial \psi }{\partial x}\) and they end with the problem \(f^{\prime \prime \prime }+\frac{1}{2}ff^{\prime \prime }-M(f^{\prime }-1)=0\), \(\frac{1}{Pr}(1+R)\theta ^{\prime \prime }+\frac{1}{ 2}f\theta ^{\prime }-Q\theta +E_{c}f^{\prime \prime 2}=0\), with the boundary conditions \(f^{\prime }=\varepsilon \), \(\theta =0\), \(Pr f+H\theta ^{\prime }=0\) at \(\eta =0\) and \(f^{\prime }\rightarrow 1\), \(\theta \rightarrow 1\) as \( \eta \rightarrow \infty \). The authors solve this last problem using a fourth-order Runge-Kutta method with a shooting method. The purpose of the paper is to analyze the influence of the viscous dissipation and of the thermal radiation, through different parameters, on the solution of the final problem. The authors illustrate their results with figures derived from the numerical simulations.