# Transformation Magnetohydrodynamics in Presence of a Channel Filled with Porous Medium and Heat Transfer of Non-Newtonian Fluid by Using Lie Group Transformations

Authors: Nazek Alessa
Journal: Journal of Function Spaces
 (2020)

---

## Abstract

In this paper, the numerical results are presented by using Lie group transformations, to be more efficient and sophisticated. To solve various fluid dynamic problems numerically, we present the numerical results in a field of velocity and distribution of temperature for different parameters regarding the problem of radiative heat, a magnetohydrodynamics, and non-Newtonian viscoelasticity for the unstable flow of optically thin fluid inside a channel filled with nonuniform wall temperature and saturated porous medium, including Hartmann number, porous medium and frequency parameter, and radiation parameter, with a comparison of the corresponding flow problems for a Newtonian fluid. Moreover, the effects of the pertinent parameters on the friction coefficient of skin and local Nusselt number were discussed numerically and also illustrate that graphically.


---

## Body

## 1. Introduction

The electromagnetic propulsion is an important application of magnetohydrodynamic flow of non-Newtonian viscoelastic fluids. Also, the magnetohydrodynamics is one of the applications of an electrically conducting fluid, as well as it will help the fluids which have thixotropic behavior in the flow of the liquid as the blood, lubrication with greases and heavy oils, and plastic extension. They have applications in chemical engineering as well.Lie group transformation was developed [1, 2] to use in mechanics of fluid and also in transformations of equations to solve some nonlinear problems.Several years ago, simple flow problems in hydrodynamics were study, and new attention has been given to the magnetohydrodynamics by Moreau [3]. Control transfer unsteadies the heat in the vertical porous channel [4]. There are many applications for the motion of Newtonian fluids in many areas, such as the flow and heat transfer magnetohydrodynamics, the nuclear fuel flow, flow of blood, and also by using the multideck asymptotic technique in Poiseuille flow [5–9]. Considering the electromagnetic propulsion is another important application. This system is consists of plasma and source of power (such as a reactor of nuclear) and also tube in which the electromagnetic forces accelerate the plasma. In studying such a system, it requires understanding the properties of transport and state equation, the relationship of shear rate-shear stress, the conductivity of electrical, thermal conductivity, and radiation. It is worth to mention that this study is closely associated with magnetochemistry.The presence of an external magnetic field will likely affect some of those properties where the plasma sets in motion of hydrodynamic [10, 11].In the last years, researchers have been interested in the issue of the flow of fluid through porous media, because the pores of the reservoir rocks recover the crude oil and at the same time because of the magnetic field effect on the work and quality of the systems. There are several studies on the transfer of heat in porous media and the flow of magnetohydrodynamics [12]. Chamkha [13] has considered the effects of Hall on a vertical plate during the flow of magnetohydrodynamic free convection through a thermally stratified porous medium. For instance, Raptis et al., in the presence of a porous medium through two parallel panels, [14] have analyzed the flow of hydromagnetic free convection [15] through vertical plate and studied the flow of mixed convection in a porous medium under the effect of a magnetic field. Through a porous medium in a channel filled with it, Makinde et al. [16] have studied the possibility of heat transfer to the magnetohydrodynamic flow inside the cannel.The present paper is aimed at modifying the study of Makinde et al. [16] for non-Newtonian fluid and investigating the combined effects of non-Newtonian viscoelasticity and radiative heat transfer through a channel filled with porous medium for optically thin fluid. In this work, we have used Tanhed method to transform a system of partial differential equations to a system of nonlinear ordinary differential equations, and then solve it numerically by using the shooting method. Numerical results were present for skin friction coefficient, velocity profile, temperature distribution, and local Nusselt number, with a comparison of the corresponding flow problems for a Newtonian fluid. The effects of the different parameters of the problem on the flow have been shown and discussed graphically.


## 2. Mathematical Analysis

Assume that the unsteady flow of non-Newtonian fluid of conducting an optically thin fluid in a porous medium inside a channel filled with it, under the influence of radiative heat transfer and field of homogeneous magnetic an externally applied (see Figure1).Figure 1
Described the problem.Suppose that the force of electromagnetic produced is too small and the fluid has a weak electrical conductivity. Leta,b system of Cartesian, such that a lies along the center of the channel, and b is the distance measured. Let a model of Boussinesq incompressible fluid, then give as the equations of governing the motion.
(1)∂z∂i=−1λ∂u∂a+ϑ∂2z∂b2+υ∂3z∂i∂b2−ϑP−ρE02λz+hαr−hαr0,∂r∂i=tλS∂2r∂b2−1λS∂v∂b.The problem’s boundary conditions are(2)z=0,r=rjatb=m,z=0,r=r0atb=0.Let the temperaturesr0,rj of both walls are high enough to transfer heat of radiative. [17] Cogley et al. supposed that the optically thin fluid has a relatively low density with a flux of radiative heat shown as follows:
(3)∂v∂i=4ϕ2r0−rj,such that ϕ is the mean coefficient of radiation absorption. Now, introducing these dimensionless variables,
(4)a¯=am,b¯=bm,z¯=zw,η=r−r0rj−r0,i′=iwm,u¯=muλϑw,such that w is the flow, and in another meaning, it is velocity. The governing equations of dimensionless with ignoring the bars (for clarity) been given by
(5)y∂z¯∂i′=−∂u¯∂a¯+∂2z¯∂b¯2+A∂3z¯∂i′∂b¯2−qz¯−Nz¯+fgη,(6)L∂η∂i′=∂2η∂b¯2+Rη.The transformed boundary conditions may be written as follows:(7)z¯=0,η=1atb=1,z¯=0,η=0atb=0,where the dimensionless parameters are defined as follows:
(8)Grashofnumber:fg=hαrjm2−hαrm02νU,Hartmannnumber:N=m2ρE2ρϑ,Radiationparameter:R=4ϕ2m2t,Pecletnumber:L=wmλSt,Darcynumber:C=Pm2,Porousmediumshapefactorparameter:q=1C,Non−Newtonianviscoelasticityparameter:A=υwmϑ.


### 2.1. Lie Group Transformations

We are using this method to solve the problems of fluid dynamics numerically, with efficiently and high quality [18].Let(9)−∂u¯∂a¯=μtanhji′,z¯b¯,i′=z0b¯tanhji′,ηb¯,i′=η0b¯tanhji′,such that μ is a constant and j is the oscillation frequency.


### 2.2. Solution Mechanism

To solve equations (5) and (6) under the boundary condition (7), replace the expressions above in equation (9) into equations (5) to (7), and we obtain
(10)d2z0db¯2−n1z0=−n2μ+fgη0,(11)d2η0db¯2+n3η0=0,with boundary conditions
(12)z0=0,η0=cothji′atb¯=1,z0=0,η0=0atb¯=0,where n1=q+N+2jycosech2ji′/1+2Ajcosech2ji′, n2=1/1+2Ajcosech2ji′, and n3=R−2jLcosech2ji′.


### 2.3. Skin Friction and Nusselt Number

The upper wall of the channel defines the shear stress (13) more clearly; is δ considered in the case that b=0?
(13)δ=wm∂z¯∂b¯+γ1w2γm2∂2z¯∂i′∂b¯b=0.The local nondimensional skin friction coefficientψ is given by
(14)ψ=δw/m=tanhji′z0′0+εsecζ2ji′z0″0,where ε=γ1wj/γm.The heat transfer rate is defined as follows:(15)vj=−t∂r∂b¯b=0.The coefficient of local heat transfer is defined as follows:(16)ζb¯=mvjb¯rj−r0.The local Nusselt number may be written as(17)Rz=ζb¯t=−tanhji′η0′0.


## 2.1. Lie Group Transformations

We are using this method to solve the problems of fluid dynamics numerically, with efficiently and high quality [18].Let(9)−∂u¯∂a¯=μtanhji′,z¯b¯,i′=z0b¯tanhji′,ηb¯,i′=η0b¯tanhji′,such that μ is a constant and j is the oscillation frequency.


## 2.2. Solution Mechanism

To solve equations (5) and (6) under the boundary condition (7), replace the expressions above in equation (9) into equations (5) to (7), and we obtain
(10)d2z0db¯2−n1z0=−n2μ+fgη0,(11)d2η0db¯2+n3η0=0,with boundary conditions
(12)z0=0,η0=cothji′atb¯=1,z0=0,η0=0atb¯=0,where n1=q+N+2jycosech2ji′/1+2Ajcosech2ji′, n2=1/1+2Ajcosech2ji′, and n3=R−2jLcosech2ji′.


## 2.3. Skin Friction and Nusselt Number

The upper wall of the channel defines the shear stress (13) more clearly; is δ considered in the case that b=0?
(13)δ=wm∂z¯∂b¯+γ1w2γm2∂2z¯∂i′∂b¯b=0.The local nondimensional skin friction coefficientψ is given by
(14)ψ=δw/m=tanhji′z0′0+εsecζ2ji′z0″0,where ε=γ1wj/γm.The heat transfer rate is defined as follows:(15)vj=−t∂r∂b¯b=0.The coefficient of local heat transfer is defined as follows:(16)ζb¯=mvjb¯rj−r0.The local Nusselt number may be written as(17)Rz=ζb¯t=−tanhji′η0′0.


## 3. Results and Discussion

We used fourth-order Runge-Kutta integration and Newton’s shooting method to solve equations (10) and (11) under the boundary conditions (12) numerically, for different values of the parameters, and used the following parameter values L=0.71, fg=1.0, y=1.0, μ=1.0, and j=1.0. Figures 2–6, we noticed the profile of velocity for various values in the parameter of the radiation R, Hartmann number N, non-Newtonian (viscoelasticity) parameter A, porous medium parameter q, and frequency parameter ji′. In Figures 2 and 3, we can notice that the profile of fluid velocity is parabolic where the minimum value concentrated at the walls, and the maximum value concentrated along the centerline of the channel. Nevertheless, we can observe that the fluid velocity profile decreases when the Hartmann number is increased and increases if the radiation parameter is increased, which means a magnetic field intensity. Figures 4 and 5 clearly show that the fluid velocity decreases when there is an increase in a porous medium and increases when there is a decrease in the non-Newtonian (viscoelasticity) parameter. In Figure 6, we can note that the fluid velocity increases when there is an increase in the frequency parameter. In Figures 7 and 8, we can note that the fluid temperature increases when there is an increase in the radiation parameter and increases when there is a decrease in the frequency parameter.Figure 2
Velocity profile for various values ofR with q=1.0, N=1.0, A=1.0, and ji′=π/4.Figure 3
Velocity profile for various values ofN with q=1.0, R=1.0, A=1.0, and ji′=π/4.Figure 4
Velocity profile for various values ofA with q=1.0, N=1.0, R=1.0, and ji′=π/4.Figure 5
Velocity profile for various values ofq with N=1.0, R=1.0, A=1.0, and ji′=π/4.Figure 6
Velocity profile for various values ofji′ with q=1.0, N=1.0, R=1.0, and A=1.0.Figure 7
Fluid temperature for various values ofR with q=1.0, N=1.0, A=1.0, and ji′=π/4.Figure 8
Fluid temperature for various values ofji/ with q=1.0, N=1.0, R=1.0, and A=1.0.The magnetic field can produce a resistive type force that appears to withstand the fluid flow, and its velocity reduces, where radiation parameter valueR increases with fixed magnetic field parameter N values (see Figures 9 and 10).Figure 9
Velocity ofZ when N=0.5, R=1.5.Figure 10
Velocity ofZ when N=0.5, R=2.5.Table1 illustrates the effects of A, N, q, R, and ji′ on the coefficient of skin friction (wall shear stress) and the number of local Nusselt (wall heat transfer rate).Table 1
Numerical values ofψ and Rz for various values of A, N, q, R, and ji′ with j=1.0, y=1.0, L=0.7, μ=1.0, fg=1.0, and ε=0.1.
ANqRji′ψ−Rz1.01.01.01.0π/40.177350.700642.00.126710.700643.00.098560.700641.01.01.01.0π/40.177350.700642.00.154480.700643.00.126800.700641.02.01.01.0π/40.154480.700642.00.136660.700643.00.114310.700641.02.01.00.5π/40.172920.618241.00.177350.700641.50.186600.876741.01.01.01.0π/40.177350.70064π/20.445491.066653π/40.550361.1620π0.576051.18286Coefficient decreases asA, q, and N increase but increases as R and ji′ increase, and the local Nusselt number increases as R and ji′ increase.


## 4. Conclusion

This paper reviewed the combined effects of a transverse magnetic field, non-Newtonian viscoelasticity, and radiative heat transfer to unsteady the flow of optically thin fluid in a porous medium inside a channel filled with it and nonuniform channel wall temperature. This paper presented the results of the temperature and velocity numerically and then used to compute the rate’s transfer of heat at the channel walls and shear stress of the wall. In general, these results show that the increase in non-Newtonian fluid parameter reduces the velocity profile and the fluid velocity is parabolic within maximum value along the centerline of the channel and minimum at the walls. Nevertheless, we observed that the fluid velocity profile increases when there is an increase in the frequency parameter, and the increase of non-Newtonian fluid parameter reduces the skin friction coefficient, while the frequency parameter increases the wall heat transfer rate. Also, it was interesting to compare the results with the corresponding flow problems for a Newtonian fluid.



---

*Source: 8840287-2020-10-22.xml*