Numerical Study of Lorentz Force Interaction with Micro Structure in Channel Flow

: The heat transfer Magnetohydrodynamics ﬂows have been potentially used to enhance the thermal characteristics of several systems such as heat exchangers, electromagnetic casting, adjusting blood ﬂow, X-rays, magnetic drug treatment, cooling of nuclear reactors, and magnetic devices for cell separation. Our concern in this article is to numerically investigate the ﬂow of an incompressible Magnetohydrodynamics micropolar ﬂuid with heat transportation through a channel having porous walls. By employing the suitable dimensionless coordinates, the ﬂow model equations are converted into a nonlinear system of dimensionless ordinary differential equations, which are then numerically treated for different preeminent parameters with the help of quasi-linearization. The system of complex nonlinear differential equations can efﬁciently be solved using this technique. Impact of the problem parameters for microrotation, temperature, and velocity are interpreted and discussed through tables and graphs. The present numerical results are compared with those presented in previous literature and examined to be in good contact with them. It has been noted that the imposed magnetic ﬁeld acts as a frictional force which not only increases the shear stresses and heat transfer rates at the channel walls, but also tends to rotate the micro particles in the ﬂuid more rapidly. Furthermore, viscous dissipation may raise ﬂuid temperature to such a level that the possibility of thermal reversal exists, at the geometric boundaries of the domain. It is therefore recommended that external magnetic ﬁelds and viscous dissipation effects may be considered with caution in applications where thermal control is required.


Introduction
The fluids involving nano-sized polymeric additives, mixed in a non-symmetric manner, are recognized as micropolar fluid. This fluid involves spinning micro fragments that provide a better mechanism for the successful employment of micropolar fluid in biotechnology and engineering. At first and almost at the same time, a hypothesis was proposed by Vogel and Patterson [1], Fabula and Hoyt [2] on the liquids encasing polymeric added substances in a small amount. Their examinations portrayed that the liquids having polymeric added substances show a reduction in the shear stresses in closeness to a rigid surface. The Newtonian fluids could not describe this characteristic. Eringen [3,4] initiated the framework of micropolar fluids. This work was further developed by Ariman et al. [5,6] who also provided examples and applications of micropolar fluids. Micropolar liquids can communicate the flow conduct of exotic lubricants, paints, polymeric materials, flow in capillaries, animal blood, and ferro liquids.
Exploration exercises expecting to investigate fluid material science at micro and nano scales have been expanding in recent years. The fluid flows in micromachined fluid systems (e.g., valves and pumps) and channels utilizing Navier-Stokes equations have been deeply analyzed in the existing literature [7]. In certain circumstances, the Navier-Stokes conditions (obtained from the classical continuum) cannot explain the fluid flow attributes at a large scale [8]. It may be due to the reason that the spinning molecules substantially affect the flow field when the molecular size is comparable to the channel size. The Navier-Stokes equation does not involve the effect of molecular spin. The theory of microcontinuum comprised of micromorphic, microstretch, and micropolar theories was established by Lee et al. [9] and Eringen [10][11][12]. This theory offered a mathematical structure to investigate such motions. The degree of freedom called gyration is used, in the micropolar theory, to find the rotation of micro-structured molecules. Therefore, the system of micropolar flow model equations consists of an additional transport equation to solve gyration. The molecular spin can be determined from this equation. The theory of micropolar fluid is an alternative persuasive methodology to numerically explore the microscale fluid dynamics, and is computationally effective rather than the molecular dynamics (MD) simulation. Papautsky et al. [13] adopted the model of micropolar fluid for the first time, and experimentally observed that the volume flow rate decreased for the flow in a rectangular microchannel.
Micropolar fluid flows have been explored analytically and numerically by several scholars. A persuasive numerical scheme "RKF 45-method" was applied by Souayeh et al. [14] to solve the flow model equations of micropolar nanofluid with allowance for warm nonlinear radiation. A framework of biothermal characteristics was established for all the preeminent parameters whose effects were elaborated through several graphs. An analytical and computational study of hybrid nano-micropolar flow was premeditated by Tassaddiq [15]. The findings of this study depicted that the macro-velocity field continuously reduced against increasing values of the material parameters. A blood flow of a two-phase model, in which blood is assumed to be a micropolar fluid, was examined by Jaiswal and Yadav [16]. A twodimensional microscopic seepage channel flow model having a real core structure was established by Wang et al. [17]. The finite element method was used to simulate the flow of solid particles in the fluid and the fact of clogging reservoir was examined. The impact of uniform magnetic field was also taken in their analysis. The authors graphically analyzed the flow model parameters for several Darcy and Brinkman regions. A micropolar flow, with allowance for thermal radiation, in a channel, was carried out by Ahmad et al. [18]. They adopted finite difference discretization along with quasi linearization to solve the dimensionless differential equations and spotted that thermal radiation caused an increase in the mass transmission rate on the lower wall of the channel. Lund et al. [19] considered the viscous dissipation and thermal radiation effects with first-order slip in a micropolar fluid flow over a linear shrinking sheet. It was noticed that the fluid velocity increased in the presence of a strong magnetic field. The fourth-order PC4-FDM (predictor-corrector finite difference method) was employed by Khader et al. [20] to solve the non-dimensional model dynamical equations of micropolar fluid. The flow was considered over a contracting or expanding surface. Mass and heat transmission rates, as well as the friction factor, were enumerated under some specific conditions. The bioconvective flow of micropolar nanofluid and gyrotactic microorganisms with the coexistence of solutal and thermal stratifications was developed by Tlili et al. [21]. The flow model problem, in this study, was innovated by using the partial slip boundary conditions. The coupled and highly nonlinear differential equations were treated numerically via the Homotopy Analysis Method (HAM). It was profound that the effect of bio-convection Lewis number was to deteriorate the local density number. Gangadhar et al. [22] developed a model problem for the flow of micropolar fluid over an extending surface. The slip conditions were also taken into account in the flow. A numerical technique, bvp4c MATLAB solvers, was used to attain numerical solutions. A micropolar fluid flow model to solve the dynamical boundary layer problems was established by Sui et al. [23]. They introduced the power-law function in this model. In this study, the numerical solutions were found using the homotopy analysis method. The interaction of applied Lorentz force with the micropolar fluid inside a magnetohydrodynamic micropump was examined by Alizadeh et al. [24]. Their outcomes were correlated with that of experimental ones and were found to be almost similar to each other for special cases. Maleki et al. [25] presented a comprehensive study describing the magnetic nanoparticles effects on microfluidics. They proposed that the shear as well as magnetic force provided a better mechanism of vortices' manipulation inside the droplet.
Flows across the channels play a pivotal role in biotechnology and industry because of their practical employments in binary gas diffusion, ablation cooling, air circulation in the respiratory system, and the combustion process in rocket motors (see [26]). The micropolar flows through channels have been deliberated by various authors. An incompressible, steady, and laminar flow of micropolar within a resistive permeable medium between walls of a channel with mass and heat deportation, by considering the effect of heat generation, was examined numerically by Ahmad et al. [27]. A micropolar flow between parallel plates of an inclined channel was analyzed by Srinivasacharya et al. [28]. The upper plate was placed at a fixed heat flux and the lower plate was maintained at a constant temperature. They used the spectral quasi-linearization method to solve the governing equations numerically. The flow of an unsteady micropolar fluid, by taking the Dirichlet boundary conditions over a domain Ωε, was scrutinized by Boukrouche et al. [29]. This flow problem was thus interpreted by the angular micro-rotation field, the pressure, and the fluid velocity. The MHD micropolar flow generated by peristaltic waves with heat and mass transfer attributes through a curved channel was determined by Ahmed et al. [30]. The lubrication approximation was utilized to reduce the relevant governing equations. The impacts of various parameters such as curvature parameter, Brinkman number, and micropolar parameter on the flow were analyzed in detail. Ding et al. [31] presented an interesting novel study of micropolar fluid flow in nanofluidic channels, in which it was argued that the microstructure nature of fluid particles had a profound effect on the electro kinetic phenomena. The effect of EDLs (electrical double layers) and velocity slip were also taken in this channel flow problem to investigate their effects on microrotation. Singh et al. [32] proposed a numerical study regarding the flow of micropolar fluid to report the effects of entropy generation on the flow within an inclined channel in the presence of thermal conductivity, changeable dynamic viscosity, and steady vortex viscosity. A channel flow involving micropolar fluid was considered by Ahmad et al. [33]. An increase in the micropolar material parameters increased the flow velocity, microrotation, and temperature while an increase in the porosity parameter produced an opposite effect.
The magneto-hydrodynamic phenomenon is widely used in the fields of astrophysics, engineering, geophysics, and aerospace engineering. Examples include photochemical reactors, plasma confinement, fiber coating, transportation, magnetic drug targeting, heat exchangers, electromagnetic casting, X-rays, cooling of nuclear reactors, sensors, and so forth. The flow of electrically conducting fluids immersed in external magnetic fields was numerically explored by Fonseca et al. [34] using the finite volume method. The influence of the magnetic field on the flow velocity, along with other parameters, was observed. The nature of the non-Newtonian pulsatile flow of micropolar-Casson fluid, influenced by Lorentz force, in a restrained channel subject to Darcy's law was surveyed by Ali et al. [35]. In their study, the maximum velocity of the flow was attained at the constriction throat of a channel. The effect of Lorentz force on a flow-through cavity was determined by Sheikholeslami and Rokni [36] by using the finite element method which was based on control volume methodology. Their results portrayed that the Rayleigh number caused an increase in convective flow. The peristaltic flow of MHD micropolar fluid inside a tapered channel, to examine the impacts of thermal radiation and entropy generation, was explored by Asha and Deepa [37]. The relevant equations of motion and heat transfer were solved by Adomian Decomposition Method (ADM). A numerical investigation of two different types of immiscible fluids (one of which was micropolar fluid) flow through a vertical channel in the coexistence of magnetohydrodynamic was presented by Tetbirt et al. [38]. They discussed the comparative results of magnetic and non-magnetic terms. Umar et al. [39] analyzed the heat transfer of micropolar fluid flow in a constricted channel influenced by thermal radiation and the Lorentz force. A finite difference-based flow solver, on a Cartesian grid, was used to find the numerical solution.
Research is carried out by Ahmad et al. [40] to comprehend how much magnetic field affects the flow through a cavity in the presence of a dipole. Additionally, it was determined that the effect of the magnetic field shifted the temperature field towards the zone of higher temperature around the dipole location. Taking the flow of micropolar fluid along a two-dimensional channel, Mirzaaghaian and Ganji [41] found that the Reynolds number significantly affected the streamwise velocity and negligibly affected the concentration and temperature. They verified the validity and precision of the Differential Transformation Method (DTM) and declared that this method is prominent to find the approximate solutions of the fluid problems. The same flow was presumed in a channel having porous walls [42]. In this problem, the concentration equation contained the chemical reaction term. The numerical solution was determined via the homotopy perturbation method (HPM). It was concluded that the Sherwood and the Nusselt number are directly proportional to the Peclet number and suction/injection parameter. An effort was taken out by Mahian et al. [43] to investigate the magnetohydrodynamic effects on the distributions of temperature, entropy generation, and velocity. The flow was taken within the gap of two concentric rotating cylinders. The governing equations were solved analytically to determine the average and local entropy generation rate. A novel study on the fractal-based approaches and fractal models to elaborate the convective heat transfer, critical heat flux, and effective thermal conductivity was presented by Cai et al. [44]. Combined power and heat systems were discussed by Mahian et al. [45]. A technique (Grey Wolf Optimization) was introduced to model several structures of the system. The downward flow under the influence of a strong magnetic field through a vertical duct was explored numerically by Zhang and Zikanov [46] using the two-dimensional approximation. Ahmad et al. [47,48] offered novel studies involving the flow of hybrid nanoparticles under the magnetohydrodynamic environment. They determined that the effect of Lorentz force suppressed the velocity and enhanced the temperature.
The current study regarding the flow of micropolar liquids within permeable media (channel in the present case) has a scope in several fields of technology and industry like microemulsions, lubrication, micro machines, petrochemical, foams and aerogels, polymer blends, and alloys. A survey of the existing literature denotes that no work has been performed so far to numerically investigate the micropolar flow through porous channels taking into account the simultaneous effects of viscous dissipation and magnetohydrodynamics.
An important characteristic, viscous dissipation, is taken into account in the recent analysis. An efficient computational algorithm is presented to obtain a meaningful numerical solution. The novel results for the fluid flow through porous channels have been investigated and physically interpreted through graphs and tables.

Problem Formulation
Consider an incompressible and steady micropolar flow through the gap between the walls of a channel. Both channel walls are porous with different permeability. The fixed pressure gradient causes the motion of the fluid between channel walls and the magnetic field is applied along the y-axis. The impacts of the body, as well as couple forces, are considered to be negligible. Taking into account these flow assumptions, we intend to investigate the two-dimensional micropolar fluid flow problem involving the magnetic and viscous dissipation effects at any channel cross-section. The components u and v are appropriately taken as velocity components are taken along and across the channel, as appeared in Figure 1. The upper wall of the channel is situated at y = −c and the lower wall is placed at y = c, so that 2c is the entire channel width. The injection velocity at the lower wall is U 1 and that taken at the upper wall is U 2 . The permeability parameter, to analyze the effects of different permeability at both walls, is defined as The injection velocities through the porous walls can be determined by this parameter. The relevant flow model equations of micropolar fluid, as suggested by Eringen [3], can be written as: The terms involved in the above equations such as V, π, l, f , ρ, v, and j represent the fluid velocity vector, pressure, body couple per unit mass, body force per unit mass, density, fluid velocity vector, and the microinertia, respectively. The material constants are denoted by λ, µ, α, β, γ, κ.
For the 2D problem under observation, the normal velocity v and the streamwise velocity u are two non-zero components of velocity. Hence, we may precise the microrotation and the velocity vectors as: Putting (5) in the Equations (2)-(4), we obtain: ∂v ∂y where κ 0 is the thermal conductivity, σ is the electrical conductivity, B 0 is the strength of the magnetic field, T is the temperature, and c p is the specific heat at constant pressure.
The boundary conditions are given as: where the injection/suction velocity at the lower and upper wall is U 1 and U 2 , respectively. It may be important to mention that one of these velocities will be positive, and the other will be negative. The microrotation and velocity field are acquired by solving Equations (6)-(10) under adequate boundary conditions (11). However, the following similarity variables are suggested, as proposed by Berman [49], to alter the relevant governing PDEs (6)-(10) into ordinary ones: where η = y c is the similarity variable, and Uc is a constant reference velocity. Furthermore, g(η) is the dimensionless angular velocity, and χ is the spatial variable along the horizontal direction. The components u and v in correlation with the stream function ψ are: The continuity Equation (6) is satisfied identically by the above-proposed velocity field and, thus, these components describe possible fluid motion. The equation of continuity is not considered further because of the elimination of the pressure term from the governing equations. Using (12)- (14) in (7)-(10), we find The notation of prime expresses the differentiation w.r.t.η. Now, boundary conditions (11) in view of (12) and (13) become: Excluding the terms of pressure from Equations (15) and (16), we get Equation (17) can be rearranged as: Equation (10) now takes the form: In the above dimensionless equations, M expresses the magnetic field parameter, Re which represents the parameter that associates the injection velocity at both walls. Furthermore, the parameters such as C 3 , C 2 , and C 1 are the material parameters describing the microinertia density, spin gradient viscosity, and vortex viscosity, respectively. These parameters are written as As U 2 is included in Re, thus the value of Re will be less than zero at the upper wall. It is noticed that, if we eliminate the microrotation, then Equation (19) will become the equation of Newtonian fluid, whereas Equation (20) reduces to zero (identically) for j = 0, as mentioned in [4][5][6]. This leads us towards the validation of our micropolar fluid flow model.

Numerical Approach
A typical approach to numerically solving the governing differential equations is as follows: Setting t = θ , s = g , p = f , q = f , r = f in Equations (19) and (20), we have With respect to the following BCs: Here, the initial conditions to be determined are g (−1), f (−1), θ (−1) and f (−1). Thus, it may be better to incorporate a shooting technique to solve the above system. This technique may associate an algorithm of fourth-dimensional zero (to determine the missing conditions) and the Runge-Kutta method (to determine the first-order differential equations). In order to satisfy the boundary conditions (e.g., p(1) = 0, f (1) = 0, g(1) = 0, θ(1) = 0), the unknown initial conditions of the original boundary value problem are computed. A troubling situation appears when the shooting method diverges the solution before the completion of the iterative process. Usually, it happens even for exact guesses of the initial conditions. This singularity may be attributed due to the reason that the solution of the problem explicitly depends on the initial conditions or the differential equations are not stable. On the other side, a finite difference approach may not suffer from these inadequacies and can easily provide a better solution as compared to the shooting and RK methods. Along with these features, we employed a finite difference scheme in our previous work (see Refs. [50][51][52]), which we describe as follows: We integrate Equation (19), which takes the following form: where β 0 denotes the integration constant to find out. This equation is further reduced into the following equation after substituting f = p: Equation (25) is solved with respect to the following BCs: To attain the numerical solution, first, Equations (20), (21) and (25) are discretized at a specific grid point by utilizing central differences instead of derivatives. Consequently, the resultant system is solved numerically by the SOR technique under the BCs (26). In the previous work, we employed a trial and error scheme to find the integration constant β 0 . It took an enormous amount of time to find β 0 by this method. Some manual calculation was also performed in each simulation. This is why we search for another persuasive approach that is based on the FDM and does not need to determine any unknown. In the present work, the Quasi-linearization method is intended to apply which is a simple numerical solution to the problem.

Numerical Solution Using Quasi-Linearization
The sequences of the vectors f (k) , g (k) , θ (k) , and φ (k) are assembled to incorporate the quasi-linearization for the system of Equations (19)- (21). Equation (19) will be linearized by defining the function Q such as: Using Taylor's series expansion, we obtain: Replacing the derivatives with the central differences in the ordinary differential Equation (28), we will obtain the equation that will determine the sequence s (k) : The matrix is linked as where the number of grid points is denoted by n. The matrices A n×n and B n×1 are initialized as follows: and for 2 < i < n − 1 and On the contrary, Equations (20) and (21) are readily linear. Hence, these equations (to find θ (k) and g (k) ) can be rearranged as: The above equations s (k+1) are assumed to be known. The computational analysis is described below:

•
The suggested guesses for s (0) , g (0) , and θ (0) are provided, so that the conditions given in Equation (26) are satisfied. • The linear system (28) is solved to gain s (1) . • The linear system (32) is obtained from the FD discretization, and then the value of known s (1) is used to get θ (1) and g (1) .

•
The iteration process is continued until and unless the numerical solutions of the Equations (24) and (25) are attained.

•
The three sequences are produced until max It may be essential to mention here that matrix A (in Equation (31)) is not diagonally dominant and is a pentadiagonal matrix. Hence, the SOR technique may not work well or fail. Thus, some other direct methods like Gaussian elimination or the LU factorization method can be used.

Results and Discussion
This section is committed to exploring the numerical outcomes with the support of physical interpretations. The flow and heat transport attributes through a permeable channel are characterized using micropolar fluid. A struggle is made to check out the impacts of pertinent parameters on the flow velocities f (η) & f (η), microrotation g(η), and temperature θ(η) as well as on shear stresses, couple stresses, and rate of heat transportation. The step sizes as well as the edges (for boundary layers) are balanced such that the flow, microrotation, and temperature profiles demonstrate an asymptotic behavior. We have taken fixed values of the parameter in the simulation (computational) analysis such as Re = 1, C 1 = 0.1, C 2 = 0.2, C 3 = 0.3, M = 1, A = 0.5, Pr = 6.2, and Ec = 0.2 otherwise identified. An excellent comparison of our numerical results (please see Table 1) with the literature for the classical Newtonian fluid is a source of validation for our computational technique. In order to validate our mathematical model and the computational scheme for non-Newtonian fluid, our numerical results for f (η) and g(η) against the fixed values C 1 = 0.1, C 2 = 1, C 3 = 0.1, M = 0, A = 0 and Re = 0.1 of the governing parameters, for the liming case, have been plotted (please see Figure 2.) against the ones reported in literature (Mirzaaghaian and Ganji [41]), by employing the Runge-Kutta method and the Differential Transform method. Excellent comparison confirms the correctness of the micropolar model and the accuracy of our numerical approach. The simulation values of the material constants are given in Table 2. The numerical consequences for several steps of η are enumerated in Table 3. This table specifies that our numerical results converge in the best way with the values of step-size η, and it confirms the accuracy of our numerical procedure. Table 4 reveals that the Reynolds number generates an opposite effect on the lower as well as an upper wall for shear and couple stresses. Its effect substantially reduces the couple and shear stresses on the lower wall (y = −c) and increases on the upper wall (y = c), whereas the heat transfer rate increases on both walls with its effect. The values of shear stresses, heat transport rates, and couple stresses on both channel walls for distinct values of the magnetic field are portrayed in Table 5. All values of f (±1) and θ (±1) seem to be increasing on both walls with the impact of M. The couple stress enhances for g (−1) and reduces for g (1). In view of mechanics, the magnetic field applies a frictional force, known as the Lorentz force. Due to this force, fluid is dragged towards the walls. This anomaly not only increases shear stresses at the walls but also tends to rotate microfluid particles rapidly. Moreover, the Lorentz force produces the temperature difference (between the temperature at channel walls and temperature of fluid) and hence increases the heat transport at both walls. Table 2. Non-dimensional parameters C 1 , C 2 , and C 3 for different cases.    Table 6 reveals the influence of permeability parameter A on heat transfer rate, couple, and shear stress for fixed values of M, Pr, Ec, Re, and material constants C 1 , C 2 and C 3 . All of the physical quantities enhance on upper as well as lower channel walls with the effect of A. The analysis of skin frictions, couple stresses, and rates of heat transfer for material parameters on both the channel walls is provided in Table 7. The diverse values of micropolar material constants are given in Table 2. The first case (see Table 2) where all values of material parameters are taken as zero represents the Newtonian case, whereas the remaining values are randomly chosen to determine their impacts on the flow [54,55]. The combined impact of C 1 , C 2 and C 3 is to devaluate the skin friction and enhance the couple stresses, whereas the joint effect of these parameters causes a reduction in the heat transport rate on the lower wall and escalation on the upper permeable wall.  The rates of heat transfer for diverse Ec and Pr values recorded in Tables 8 and 9 portray that the Eckert number upsurges the heat transport rate on both walls; however, the Prandtl number causes a decrease in θ (−1) and an increase in θ (1). This study leads to the results that the heat transportation rates, shear, and couple stresses may be accommodated by selecting suitable parametric values and the permeability on both walls to achieve the desired consequences.  The streamwise velocity f , the angular velocity g, and the normal velocity f are depicted in Figure 3 for the assigned values of M and Re. The velocity f initially decreases and then increases gradually for Re when concavity changes, but, contrarily, the velocity f initially increases and then decreases gradually for M when concavity changes (see Figure 3a,b). It can be seen from Figure 3c,d that the effect of the Reynolds number is to downturn the normal velocity, but the magnetic field parameter enhances the normal velocity. In the same way, the angular velocity is increased by the Reynolds number and decreased by the magnetic field parameter (as depicted in Figure 3e,f). It is worth mentioning here that the value of a magnetic parameter is fixed (e.g., M = 10) in Figure 3a,c,e. Similarly, the value of Reynolds number is fixed (e.g., Re = 10) in Figure 3b,d,f. It is observed here that both the magnetic parameter and the Reynolds number yield opposite effects on the profiles of streamwise, angular, and normal velocities. The retarding force elevates with amplifying values of magnetic field parameters. Consequently, the magnetic field suppresses the streamwise velocity (after a change in concavity) and improves the heat transport rate. The effect of A (permeability parameter) on non-dimensional velocities is elucidated in Figure 4a,c,e. The streamwise velocity and the microrotation increase, whereas the normal velocity reduces with an increase in the permeability parameter. The same velocities ( f , f and g) are displayed in Figure 4b,d,f for various guesses of micropolar material parameters C 1 , C 2 &C 3 (see values in Table 2). The numerical results designate that the normal velocity is depressed while streamwise velocity (after concavity) and microrotation increase by escalating the material parameters. Velocity f tends to decrease at the start, but, after concavity, it attains maximum value near the upper wall of the channel (Figure 4b). It may be attributed because the permeability drives streamwise moving fluid beyond the wall and causes an increase in the velocity near the upper wall. The microrotation profile for micropolar parameters is concave upward in the upper half of the channel and concave downward in the lower half of the channel. When C 1 = C 2 = C 3 = 0 at both permeable walls, the effect of skin frictions equally spreads from the walls such that the angular velocity (microrotation) is zero at the center of the channel (see Figure 4f). The skin friction tends to rotate the micropolar fluid in opposite routes on the walls. This is why the microrotation has opposite signs near the lower channel wall.
Temperature θ(η) is depicted in Figure 5a,b for multiple values of Ec (Eckert number) and Pr (Prandtl number) as a function of the coordinate η. For both of these parameters, temperature enhances near the upper wall of the channel. The results of Figure 5c expresses that the temperature profile θ accelerates with an increment in the material parameters, while Figure 6 shows the velocity distribution across the channel for Newtonian and micropolar fluids.

Conclusions
An inclusive computational analysis of a micropolar flow with an allowance for heat transfer through a channel is presented. The effects of magnetohydrodynamics (MHD) and viscous dissipation have also been taken in the flow. The quasi-linearization is utilized to attain the solution against microrotation, velocity, and thermal energy. The key points of the present work may be mentioned as:

•
The permeability parameter and the material parameters tend to enhance the microrotation, but these parameters depreciate the normal velocities.

•
The Reynolds number inserts a low effect on couple stresses while it yields a significant effect on skin friction and heat transport rates.

•
The permeability parameter A substantially enhances shear stresses, couple stresses, and the rates of heat transfer on both the channel walls.

•
The impact of the Eckert and the Prandtl number is to uplift the temperature curves.

•
The micropolar constants C 1 , C 2 &C 3 intensively affect the microrotation rather than the streamwise and the normal velocities.

•
The micropolar fluid causes an escalation in couple stresses and a reduction in the shear stresses. Data Availability Statement: Data will be provided by the authors on request.