New traveling wave solutions of MHD micropolar fluid in porous medium

*Correspondence: jqrza@neduet.edu.pk; jqrza26@yahoo.com Department of Mathematics, NED University of Engineering & Technology, Karachi 75270, Pakistan Abstract This paper aims to determine the exact solutions of non-Newtonian fluid namely micropolar fluid with MHD in a porous medium by traveling wave method. The governing equations of incompressible micropolar fluid with MHD in a porous medium are non-linear PDEs reduced to ODEs through wave parameter ξ = mx + ny + Ut. The set of new exact solutions are determined for five different cases. In special cases, the solution for micropolar fluid with and without MHD and porous effects can also be obtained from general solutions. Furthermore, these solution reduces to a Newtonian solution if we put vortex viscosity κ → 0. Finally, the influence of the material and other parameters of interest on the fluid motion, as well as a comparison among micropolar and Newtonian fluids is also analyzed by 2D and 3D graphical illustrations.


Introduction
In the present, many researchers are working on the non-Newtonian fluid from both essential and sensible point of view [1]. These fluids have immediate effects on the processing of polymer, animal blood, liquid crystal, and geological flows in the earth mantle. The general conditions of non-Newtonian fluid are exceptionally non-linear and higherarrangement than Navier-Stokes equations. Therefore, many analytical and numerical solutions are accessible to non-Newtonian fluid on the topic.
The electrically conducting fluid and magnetic properties are sufficiently studied in magnetohydrodynamics (MHD). MHD has specific applications like engineering science, metallurgical industry, electromagnetic pump, power generation, and flow meter. Hydromagnetic movements have a major part of study in the fields of the aerospace, astronomical, and planetary magnetosphere. The cleansing of liquid metals from non-metallic presence through the use of the attractive field is another basic component of MHD [2,3]. Khalid et al. [4] have evaluated the exact solution of wall couple stress in MHD by Laplace transform and convolution. Fatunmbi et al. [5] have numerically studied the MHD stagnation point flow of a micropolar fluid by applying RK (Runge-Kutta) method. Hammouch [6] investigated the numerous solution of steady magnetohydrodynamics flow of dilatant (2020) 28:23 Page 2 of 22 fluids by various method techniques. Interesting recent studies on non-newtonian fluids with magnetohydrodynamics have been given by different researchers [7][8][9][10][11]. Micropolar fluid theory was proposed by Eriggen [12]. Mekheimer et al. [13] showed that the micropolar fluid model is a kind of non-Newtonian fluid which depends upon a microstructure and belongs to non-symmetrical stress tensor. Physically, micropolar fluid may be rigid particles, at random oriented (spherical) elements suspended in a viscous medium where the change of fluid elements is disregarded. Micropolar fluid can perform a better model for animal blood. Raza et al. [14] studied the different branches of the solution of micropolar fluid in a channel with permeable walls. Hussanan et al. [15] examined the heat and mass transfer in a micropolar fluid. Sheikholeslami et al. [16] evaluated the micropolar fluid flow in a channel by homotopy perturbation method (HPM). Lukaszewicz [17] explored the nature of the micropolar fluid. Ravi and Prasad [18] studied the relations of pulsatile and peristaltic stream of couple stress fluid through a porous medium in an elastic route.
There are many methods for solving non-linear partial differential equations such as bilinear transformation [19,20], homotopy perturbation method [21,22], (G \ G) expansion method [23], Exp-function [24,25], modified simple equation (MSE) method [26], and fingero-imbibition method [27]. The traveling wave method is useable for solving non-linear partial differential equations because it provides a simple exact solution of non-linear partial differential equations. Due to the many importance of this method in theory of non-linear partial differential equations, this method is applicable in the field of fluid mechanics, chemical kinematics, electromagnetic, and nonlinear mechanics. Some important work-related traveling wave solutions recently appeared in [28][29][30][31][32][33][34] and references presented there.
In the present communication, the traveling wave solution, and transformation method are used to find the exact solutions of unsteady incompressible two-dimensional laminar flow of MHD micropolar fluid in a porous medium. The mathematical model that considered here are carefully entrenched and judge with many general cases. It is noticed that the present model is capable of covering such benchmark cases and can give new information on the parameter determination of the dynamical system which governed micropolar fluid equations. To the best of our knowledge, only one study, which deals with the flows of unsteady incompressible two-dimensional laminar flow of micropolar fluid via traveling wave transformation appeared in [28]. Therefore, this communication aims to extend the results of [28] for unsteady two-dimensional incompressible laminar flow of MHD micropolar fluid in a porous medium. The family of new exact solutions is found for five distinct cases. Particularly, the solutions for micropolar fluid with and without MHD and porous effects can also be established from present general solutions. The solutions for MHD Newtonian fluid in the porous medium may for all cases easily be determined if we substitute vortex viscosity κ → 0 in respective equations. Finally, the impact of parameters defines micropolar fluid and other parameters are analyzed through 2D and 3D graphical illustrations and differences among micropolar and Newtonian fluids are also discussed. For the present paper, we adopted the method [28][29][30][31][32][33] in the following pattern. In the "Basic governing equations" section, we provide the basic governing equations and the theoretical development of the traveling wave method. In the "Traveling wave solutions" section, we present the solution of our governing equations by the traveling wave methods along with five different cases for MHD micropolar fluid in a porous medium. Some initial and boundary conditions are discussed in the "Making initial and boundary conditions" section. In the "Numerical results and discussions" section, we discuss the numerical results with the help of graphical illustrations.

Basic governing equations
The equation of continuity and the conservation equations of linear momentum and angular momentum for an incompressible unsteady micropolar fluid, in the presence of magnetohydrodynamic through a porous medium, by neglecting the body force and body couple are as follows: In the above equations, − → V is a velocity; − → is a microrotation; p be the pressure; ρ is the fluid density; j be the gyration parameters of the fluid; μ and κ are the dynamic and vortex viscosities; α, β, and γ are the respective coefficients of coupled viscosities [29]; − → B is the total magnetic force vector; − → J is the current density; and κ * is the porosity parameter. The velocity, microrotation, and magnetic component porous plate are where σ is the electrically conductively of fluid and then Eqs. (1-3) provide the following governing equations (see Fig. 1

Procedure of traveling wave method
Suppose the non-linear partial differential equations include four dependent variables u , v, N, and p three variables t, x, and y.
Here, M i is a polynomial function of u, v, N, and p which contain the non-linear terms and higher order of derivative; where i = 1, 2, 3, 4., we display traveling wave solution as follows: where ξ = mx + ny + Ut. The system Eqs. (8)(9)(10)(11) can be converted into the system of ordinary differential equations.
This system of ordinary differential equations are may or may not be solvable. In our case, these system are solvable.

Traveling wave solutions
In this section, we will determine the traveling wave solutions of governing Eqs. (4-7). The traveling wave parameter type can be taken, and then, the solution has the following representation.
Here, U is the constant phase velocity and m and n are constants. On using the traveling wave parameter ξ into Eqs. (4-7), we get representation of the form where the prime denotes the differentiation with respect to ξ . Integration to Eq. (17), w.r.t ξ on both side yields where c 0 is an arbitrary constant. Using Eq. (21) in Eqs. (18)- (20) reduced to To find u, v, N, and p from the above three equations, the following five cases will be discussed.

Making initial and boundary conditions
Some exact traveling solutions are determined for Eqs. (4)-(7) in the "Traveling wave solutions" section, without initial and boundary conditions. Definitely, the initial and boundary conditions can be made from these solutions. For instance, if governing Eqs.
Similarly, we can made the initial and boundary conditions to the other many cases.

Numerical results and discussions
In this paper, we obtain exact traveling wave solutions of MHD flow of micropolar fluid in a porous medium. The methodology in the recent work is easily reduced to the non-linear partial differential equations of MHD micropolar fluid to linearized ordinary differential equations. The method has been used directly without any restrictive assumption and laborious calculation. These exact traveling solutions are obtained for five different cases. From these solutions, we can find special solutions for with and without MHD and porous effects. It is also noted that when κ, N −→ 0, we obtain the solutions for viscous fluid flow. Furthermore, in all cases, we write general solutions in the exponential form; however, we can also write general solutions in the oscillatory form by assuming some suitable restriction on the roots of the equations, but we avoid it due to the length of the paper. Now, in order to reveal some relevant physical aspect of the determine results the diagram of the velocity fluid components u(x, t), v(x, t), and the microrotation N(x, t) are present against x for different values of t and of the important parameters of the micropolar fluid. For the sake of convenience, we present a graph only for case I, U + c 0 = 0, and a similar prediction can be done by other cases. Figures 2, 3, and 4 showed the influence of time, space variable x, and three-dimensional representation of u, v, and N. From Figs. 2 and 3, it is clear that u is increasing function of time t and space variable x; however, v and N are also increasing functions (in absolute value) of these variables. The combined effects of these variables are shown in 3D graphs of Fig. 4. It is more clear from 3D pictures of u, v, and N that they are becoming strengthen with the strengthening values of t and x. The influence of the vortex viscosity κ is shown in Fig. 5; from these figures, we can see that velocity components u, v, and microrotation N (in absolute sense) increase with the increasing values of κ. The influence of kinematic viscosity ν is shown in Fig. 6. The effect of kinematic viscosity ν on a fluid motion is quite opposite to that of vortex viscosity κ. The effect of magnetic parameter B 0 is shown in Fig. 7. It is clear that the effect of the magnetic parameter B 0 is similar to u and v. For instance, the velocity component u is increasing while v is also increasing function (in absolute value) of B 0 because the particle of micropolar fluid always travel rotationally, which may cause increment in velocity. However, the influence of magnetic parameter B 0 on the microrotation profile N is decreasing in the two-third domain and increasing in the remaining domain. In Fig. 8, the effect of porosity parameter k * is depicted on the velocity components and microrotation profiles. It is noted that the porosity parameter k has opposite effects on the profiles of u and v in comparison to the magnetic parameter B 0 . The microrotation N is decreasing the function of the porosity parameter. The influence of the parameter m appears in the traveling wave parameter ξ is shown in Fig. 9. It comes to the notice that increasing values of the parameter m reduce the values of u and N; however, the values of v are increased very shortly. Finally for comparison, the profile of the velocity components u(x, t) and v(x, t) corresponding to the motion of the three types of fluids (MHD micropolar fluid in porous medium, micropolar, and MHD Newtonian in porous medium fluids) are together presented in Fig. 10 for similar values of material constants. It is clear from these figures that MHD Newtonian in a porous medium fluid is the fastest, and the simple micropolar fluid is slowest. It is also brought to the knowledge that magnetic and porosity parameters fasten the fluid motion in this case that we have considered here. SI units are used in making all these graphs and Mathcad software is used for making these graphs.