Magnetohydrodynamic (MHD) Flow of Micropolar Fluid with Effects of Viscous Dissipation and Joule Heating Over an Exponential Shrinking Sheet: Triple Solutions and Stability Analysis

Liaquat Ali Lund,Zurni Omar,Ilyas Khan,Jawad Raza,El-Sayed M. Sherif,Asiful H. Seikh

DOI: https://doi.org/10.3390/sym12010142

2020-01-10

Symmetry

Abstract:A numerical study was carried out to examine the magnetohydrodynamic (MHD) flow of micropolar fluid on a shrinking surface in the presence of both Joule heating and viscous dissipation effects. The governing system of non-linear ordinary differential equations (ODEs) was obtained from the system of partial differential equations (PDEs) by employing exponential transformations. The resultant equations were transformed into initial value problems (IVPs) by shooting technique and then solved by the Runge–Kutta (RK) method. The effects of different parameters on velocity, angular velocity, temperature profiles, skin friction coefficient, and Nusselt number were obtained and demonstrated graphically. We observed that multiple solutions occurred in certain assortments of the parameters for suction on a surface. The stability analysis of solutions was performed, and we noted that the first solution was stable while the remaining two solutions were not. The results also showed that the velocity of the fluid increased as the non-Newtonian parameter rose in all solutions. Furthermore, it was detected that the temperature of fluid rose at higher values of the Eckert number in all solutions.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to consider the influence of Joule heating and viscous dissipation effects when a micropolar fluid flows on an exponentially contracting surface in a magnetohydrodynamics (MHD) environment. Specifically, the researchers are concerned with:

1. **Multiple - solution phenomenon**: The paper explores whether multiple solutions will occur under certain parameter combinations and whether these solutions are stable.
2. **Stability analysis**: For the multiple solutions found, a stability analysis was carried out to determine which solutions are stable.
3. **Parameter influence**: The influence of different parameters (such as magnetic parameter \(M\), micropolar parameter \(K\), suction/injection parameter \(S\), Prandtl number \(Pr\) and Eckert number \(Ec\)) on velocity, angular velocity, temperature distribution, wall - friction coefficient and Nusselt number was studied.

### Main objectives
1. **Expand previous work**: Expand the work of Aurangzaib et al. by considering the influence of magnetic field, viscous dissipation and Joule heating.
2. **Find multiple solutions**: Find the maximum number of multiple solutions.
3. **Stability analysis**: Conduct a stability analysis on multiple solutions to determine the stable solutions.

### Mathematical model
The following mathematical models are used in the paper to describe the problem:
- **Continuity equation**:
\[
\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y} = 0
\]
- **Momentum equation**:
\[
u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=(\nu+\kappa)\frac{\partial^{2}u}{\partial y^{2}}+\kappa\frac{\partial N}{\partial y}-\sigma B^{2}\frac{u}{\rho}
\]
- **Angular momentum equation**:
\[
u\frac{\partial N}{\partial x}+v\frac{\partial N}{\partial y}=\gamma\frac{\partial^{2}N}{\partial y^{2}}-\kappa(2N+\frac{\partial u}{\partial y})
\]
- **Energy equation**:
\[
u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha\frac{\partial^{2}T}{\partial y^{2}}+(\mu+\kappa)\frac{1}{\rho C_{p}}(\frac{\partial u}{\partial y})^{2}+\frac{\sigma B^{2}}{\rho}u^{2}
\]

### Boundary conditions
\[
\begin{cases}
v = -v_{w}(x),&u = u_{w}(x),&N=-m\frac{\partial u}{\partial y},&T = T_{w}(x),&\text{at }y = 0\\
u\rightarrow0,&N\rightarrow0,&T\rightarrow T_{\infty},&\text{as }y\rightarrow\infty
\end{cases}
\]

### Similarity transformation
Use the following similarity transformation:
\[
\begin{cases}
u = U_{w}e^{x/l}f'(\eta)\\
v = -\sqrt{\frac{\nu U_{