Computation of non-Newtonian quadratic convection in electro-magneto-hydrodynamic (EMHD) duct flow with temperature-dependent viscosity

Abstract

Motivated by novel developments in smart non-Newtonian thermal duct systems, a theoretical study has been presented in this article for electro-magneto-hydrodynamic (EMHD) buoyancy-driven flow of a fourth-grade viscoelastic fluid in a vertical duct with quadratic convection. The viscosity of the fourth-grade fluid model is assumed to be temperature-dependent, and the Reynolds exponential model is deployed. Viscous heating and Joule dissipation effects are included. The duct comprises a pair of parallel electrically insulated vertical flat plates located a finite distance apart. Via suitable scaling transformations, a nonlinear boundary value problem is derived for the momentum and heat transport. A homotopy perturbation method (HPM) solution is obtained coded in Mathematica symbolic software. There is a considerable enhancement in wall skin friction with an increment in fourth-grade fluid parameter, Brinkman number, electrical field parameter, thermal buoyancy parameter, and quadratic thermal convection parameter. However, skin friction is strongly reduced with a rise in variable viscosity parameter, Hartmann (magnetic) number, and electromagnetic heat generation to conduction ratio. Nusselt number magnitudes are elevated with increase in variable viscosity parameter, thermal buoyancy parameter, and quadratic thermal convection parameter, whereas they are significantly suppressed with increment in fourth-grade fluid parameter, Brinkman number, and Hartmann magnetic number.

Keywords

Fourth-grade viscoelastic fluid EMHD quadratic convection temperature-dependent viscosity viscous and Joule heating homotopy perturbation solution

Introduction

Magnetohydrodynamics (MHD) has emerged as a major discipline in 21st century engineering sciences, tracing its roots to astrophysics.¹ It involves the interaction between magnetic fields and either viscous or inviscid electrically conducting fluids. In modern energy systems a major application is the heat transfer and flow control in nuclear ducts and other energy systems (solar) where excessive temperatures can arise. MHD ducts flows feature in lithium blankets,² turbulence suppression,³ magnetic sensors,⁴ ferromagnetic solar absorber collectors,⁵ liquid metal flow couplers,⁶ coolant deployment in sodium fast reactors (SFRs),⁷ and pressure drop regulation.⁸ In mathematical models of MHD duct flows frequently dimensionless parameters arise which include the Hartmann number (ratio of magnetic body force to viscous force), Stuart magnetic interaction number (ratio of magnetic force to inertial force), and magnetic Reynolds number (when induction effects are significant). Electrohydrodynamics (EHD)⁹ is a sub-discipline of the more general continuum electromechanics and also known as electro-kinetics. It focuses on the interaction between electrical fields and electrically charged liquids. While in MHD usually transverse magnetic fields are utilized to generate a perpendicular Lorentz electromagnetic body force, in EHD the electrical field is often implemented along the axis of the conduit or duct. Many exciting applications for EHD duct flows (e.g. parallel plate geometries) have emerged in recent decades including atomizers,¹⁰ corona-based duct flows,¹¹ electrode gas pump channels,¹² electrostatic precipitators,¹³ dielectric barriers for electro-fluid injection,¹⁴ EHD heat exchangers,¹⁵ coolants for circuits,¹⁶ thermal enhancement,¹⁷ and bio-inspired microfluidic pumping systems.¹⁸ As with MHD, some popular electrohydrodynamic dimensionless numbers arise in numerical and theoretical studies of EHD duct flows and these include the electrical Hartmann number and electrical Reynolds number. The many studies communicated on either MHD duct or EHD duct flows have demonstrated that generally transverse magnetic field damps flows whereas axial (longitudinal) electrical field accelerates flows. These findings instigated interest in combining both electrical and magnetic fields in the science of electro-magnetohydrodynamics (EMHD) which emerged in the late 1960s.¹⁹ The most comprehensive treatment of the mechanics of EMHD flows was presented by Eringen and Maugin²⁰ who considered many types of fluids including ferromagnetic, dielectric, electro-viscous, ferroelectric, electro-strictive, electro-magneto-rheological (EMR), and hybrid magnetic ionic liquids. Subsequently many researchers have investigated the combined effects of magnetic and electrical fields (both static and alternating) in viscous duct flows. Jian and Chang²¹ developed theoretical solutions for EMHD incompressible viscous flow in a horizontal a slit microchannel duct under a lateral uniform electrical field and a spatially non-uniform vertical magnetic field. They deployed the variation of parameter method (VPM), Gaussian quadrature and Chebyshev spectral collocation techniques, and obtained good correlation with experiments. They found that axial flow deceleration is induced with increasing magnetic field decay factor whereas strong acceleration is generated with larger electrical field strength parameter. Chakraborty and Paul²² computed the EMHD transport in a microchannel, showing that at relatively weak magnetic field intensities, significant enhancement in volumetric flow rate is produced. They also noted that excessive electrical field can ramp up the Joule heating and induce wall degradation. Umavathi and Bég²³ investigated the thermo-solutal EMHD nanofluid convection in a vertical duct for nuclear cooling applications. They showed that the polarity of the electric circuit strongly modifies the velocity distribution and interaction with magnetic field. For negative electrical field, the flow is directed downwards whereas for positive electrical field it is upwards. For the case of a short circuit (vanishing electrical field), a linear velocity distribution is computed. They also found that skin friction is depleted at the left duct wall with increasing Hartmann number and nanoparticle solid volume fraction but enhanced with those these parameters at the right wall. They also examined the relative performance of copper, titanium oxide, and silver nanoparticles and noted that silver nanoparticles generate the highest Nusselt number, whereas titanium oxide nanoparticles produce the lowest. Nusselt numbers were also shown to be suppressed with electrical field but enhanced with magnetic field.

The above studies ignored non-Newtonian effects. However, many studies have identified that in for example nuclear slurry duct flows and solar direct absorber tubes, the presence of particles suspended may lead to rheological behavior. This has further motivated researchers to explore the use of different non-Newtonian formulations to better characterize the nonlinear shear-stress-strain behavior of such fluids. Soundalgekar and Aranake²⁴ used the Stokeian polar (couple stress) rheological model to study heat transfer in vertical duct flows. Muzychka and Edge²⁵ developed an elegant model for power-law liquids based on the Rabinowitsch-Mooney formulation and computed volumetric flow rates in a variety of ducts including circular and elliptic tubes, parallel channels, rectangular ducts, isosceles triangular ducts, circular annular, and polygonal ducts. Other relevant studies include Chaudhuri et al.²⁶ (on forced convection heat transfer in third-grade viscoelastic fluids in ducts) and Rajagopal and Na²⁷ (third grade fluids with thermal buoyancy effects in ducts). All these studies confirmed the significant deviation in heat and momentum characteristics due to non-Newtonian behavior. A more sophisticated viscoelastic model of the differential type is the fourth-grade model. This is superior to the second grade and third grade Reiner-Rivlin models. In addition to predicting normal stress differences, the fourth-grade model can also accommodate the shear thinning or thickening characteristics of actual industrial slurries since it does not assume shear viscosity to be a constant). Several investigations have implemented this model in thermal duct flows. Kazemi et al.²⁸ studied the hydromagnetic flow of a fourth-grade viscoelastic fluid in a duct with forced convection. They showed that there is a decrement in average velocity and the bulk temperature inside the duct with an increase in non-Newtonian parameter. Fatima et al.²⁹ presented perturbation solutions for two-phase ionic fourth-grade flow in a convergent/divergent duct containing spherically homogeneous metallic particles under an axial electrical field. Nadeem and Ali³⁰ conducted a second law thermodynamic analysis of fourth-grade viscoelastic flow with heat transfer in a vertical pipe using a homotopy method. Further investigations include Sajid et al.³¹, Salawu et al.³² (on thermal stability of magnetized fourth-grade fluids), and Khan et al. (on bio-inspired pumping of chemically reacting fourth-grade fluids with species diffusion). These studies generally showed that with increasing viscoelastic material parameter associated with the fourth-grade model, deceleration is induced in the flow.

In nuclear and other energy duct applications, complex thermophysical phenomena arise. These include thermal dispersion, stratification, radiative heat transfer, and convective heating at the boundaries. Another interesting characteristic is quadratic convection (nonlinear convection or non-Boussinesq convection) which is associated with nonlinear density temperature differences in the buoyancy force term in natural convection. This generalizes the conventional Boussinesq approximation and provides more sophistication in thermal modeling of such systems. Temperature gradients can be significantly modified with quadratic effects and therefore considering the nonlinear variation of density with temperature is preferable to the simpler approach of a linear variation of density with temperature. A number of articles have addressed quadratic convection in multi-physical flows. Kunnegowda et al.³³ used homotopy perturbation method to compute the MHD viscoplastic micro-duct transport with an exponential space dependent heat source and quadratic convection effects. Triveni and Mahanthesh³⁴ simulated the nanofluid flow in an annular duct with quadratic convection and quadratic thermal radiation. Kumar and Sood³⁵ studied the non-Boussinesq mixed convection stagnation point flow on a contracting surface to a porous medium permeated with a magnetic field. They showed that increasing quadratic convection parameter strongly modifies the transport and both accelerates the boundary layer and encourages thermal diffusion. Okoya et al.³⁶ computed the natural convection flow from a moving wall with variable thermal conductivity and quadratic non-Boussinesq effects. Al-Kouz et al.³⁷ evaluated quadratic convection and thermal radiation influences on second order viscoelastic boundary layer flow with a bvp5c finite difference algorithm in MATLAB. They showed that flow is accelerated with stronger quadratic convection and a cooling effect is induced. Further studies include Sabu et al.³⁸ (on inclined plate MHD nanofluids), Bhatti et al.³⁹ (on viscolastic Maxwell radiative-convection flow from a stretching wall with thermal relaxation effects), and Mahanthesh et al.⁴⁰ (on dusty nanofluid convection from a vertical wall).

In the present work we theoretically study the fourth-grade EMHD natural convection flow in a vertical parallel plate duct with temperature-dependent viscosity and quadratic convection. Previous studies have considered alternative simpler rheological models including the Williamson model,⁴¹ Phan-Thien-Tanner (PPT) model,⁴² couple stress model,⁴³ and second order Reiner-Rivlin model.^44,45 Furthermore, we incorporate the effects of viscous dissipation and Joule heating (Ohmic dissipation)^46,47 and the Reynolds exponential temperature-dependent viscosity model^48–52 all of which can exert marked modifications in nuclear duct electro-magneto-hydrodynamic flows. The novelty of the present work is therefore the collective consideration of fourth-grade viscoelasticity, viscous heating, exponential viscosity, and Joule dissipation, which have not been considered hitherto anywhere in the scientific literature, despite immediate applications in nuclear duct electromagnetic flow control. Fourth-grade fluids are materials with exceptional flow characteristics that fall between regular liquids and solids. These materials exhibit a complicated rheological behavior known as non-Newtonian flow, in which the viscosity (resistance to flow) varies with shear rate, temperature, and pressure. Fourth-grade fluid flow has applications in numerous scientific disciplines, including physics, chemistry, engineering, and material science. Understanding fourth-grade fluid is essential for comprehending the behavior of complex fluids such as emulsions, foams, and suspensions. Scientists can design sophisticated materials with the desired properties for a variety of applications, including drug delivery systems, cosmetics, and paint formulations, by studying the flow properties of these substances. Viscoelastic fourth-grade rheology is also applicable to geology and the study of geological formations. It explains phenomena such as lava flow, mudslides, and glacier movement. By comprehending non-Newtonian dynamics, scientists can predict and manage such natural occurrences more precisely. Fourth-grade fluid flow is crucial in a variety of other industries, including food processing, polymer manufacturing, and hydrocarbon extraction. Understanding the flow behavior of non-Newtonian substances is crucial for optimizing manufacturing processes, regulating product quality, and increasing productivity. In particular the modern thrust in electromagnetic smart polymer-based conducting materials in high-risk industries including nuclear reactor technology has mobilized great interest in electro-magneto-fluid dynamics of viscoelastic fourth-grade liquids. This provides a strong justification for the current study. The normalized nonlinear ordinary differential conservation equations and boundary conditions are solved with a homotopy perturbation technique. Detailed solutions for axial velocity, temperature, skin friction, and Nusselt number are presented graphically and in tables for the effects of key electromagnetic, thermal, and dissipation parameters. Detailed interpretation is included. The simulations as noted earlier are pertinent to providing a deeper understanding of laminar transport phenomena in nuclear and solar ducts deploying advanced functional rheological working liquids.

Mathematical model for EMHD non-Newtonian quadratic dissipative duct flow

The regime under investigation is visualized in Figure 1. An electromagnetic fourth-grade fluid flows under natural convection between the vertical parallel electrically insulated plates which are located a finite distance 2d apart. A Cartesian coordinate system ( $x'$ , $Y'$ ) is used. The walls of the duct are located at $x' = d$ and $x' = - d$ and are sustained at constant but non-similar temperatures $θ_{2}$ and $θ_{1},$ where $θ_{1} > θ_{2}$ . The fourth-grade fluid is electrically conducting and incompressible. An axial electric field (E) and a transverse magnetic field ( $\overset{⇀}{B}$ ) are imposed which are mutually orthogonal. Maxwell displacement, Hall current, and magnetic induction effects are ignored.

Figure 1.

Physical model for EMHD fourth-grade fluid flow in a vertical parallel duct.

Under these assumptions, the mass conservation (continuity) equation takes the form:

\nabla \cdot \vec{U} = 0,

(1)

The momentum equation using Ohm’s law is defined as:

\begin{matrix} ρ (\vec{U} \cdot \nabla \vec{U} + \frac{\partial \vec{U}}{\partial t'}) = - \nabla \cdot p' + \nabla τ' + \overset{⇀}{J} \times \overset{⇀}{B} + g \\ [{(ρ \bar{β})}_{0} (θ - θ_{m}) + {(ρ \bar{β})}_{1} {(θ - θ_{m})}^{2}], \end{matrix}

(2)

Here the stress tensor is designated by $τ'$ , pressure is denoted by $p',$ g is gravity, linear and nonlinear thermal expansion coefficients are designated by ${\bar{β}}_{0}, {\bar{β}}_{1}$ , $θ$ is the temperature, $θ_{m} [= (θ_{1} + θ_{2}) / 2]$ is the mean temperature, density is designated by $ρ$ , time is $t'$ , and magnetic field is designated by $\overset{⇀}{B}$ . The current density vector is designated by $\overset{⇀}{J}$ and it is expressed as:

\vec{U} \times \overset{⇀}{B} + \overset{⇀}{E}),

(3)

Here electrical conductivity is designated by $σ$ , and electrical field strength is denoted by $\overset{⇀}{E}$ .

The constitutive equation for fourth-grade Reiner-Rivlin differential viscoelastic fluid is expressed following⁵³ as:

\begin{matrix} τ' = \sum_{i = 1}^{4} S_{i}, \\ S_{1} = μ (T) ℏ_{1}, \\ S_{2} = α_{1} ℏ_{2} + α_{2} ℏ_{1}^{2}, \\ S_{3} = β_{1} ℏ_{3} + β_{2} (ℏ_{1} ℏ_{2} + ℏ_{2} ℏ_{1}) + β_{3} (trace ℏ_{2}) ℏ_{1}, \\ \begin{array}{l} S_{4} = γ_{1} ℏ_{4} + γ_{2} (ℏ_{3} ℏ_{1} + ℏ_{1} ℏ_{3}) + γ_{3} ℏ_{2}^{2} + γ_{4} (ℏ_{2} ℏ_{1}^{2} + ℏ_{1}^{2} ℏ_{2}) \\ + γ_{5} (tr ℏ_{2}) ℏ_{2} + γ_{6} (tr ℏ_{2}) ℏ_{1}^{2} + [γ_{7} tr ℏ_{3} + γ_{8} tr ℏ_{2} ℏ_{1}] ℏ_{1}, \end{array} \end{matrix}}

(4)

Here $α_{j}, β_{j}, γ_{j}, j = 1, 2, 3, . . ., 8$ indicate the material constants, and the kinematical tensors $ℏ_{1}, ℏ_{2}, ℏ_{3}$ are expressed as:

\begin{matrix} ℏ_{1} = Ξ + Ξ^{t}, Ξ = grad \vec{U}, \\ ℏ_{n} = \frac{d ℏ_{n - 1}}{d t'} + Ξ^{t} ℏ_{n - 1} + ℏ_{n - 1} Ξ, n = 2, 3, \dots \end{matrix}}

(5)

The energy equation considering the effects of Joule heating and viscous dissipation is specified as:

(ρ c_{p}) \frac{d θ}{d t^{'}} = k g r a d θ + τ^{'} : g r a d \vec{U} + \frac{\overset{⇀}{J} \cdot \overset{⇀}{J}}{σ},

(6)

Here thermal conductivity of the fluid is designated by $k$ , $c_{p}$ designates the specific heat, and the symbols “·” and “:” represents the single and double dot vector products.

For the proposed problem, we seek the velocity field and the temperature field in the following form:

\vec{U} = [0, V (x'), 0], θ = θ (x') .

(7)

The general vector equations may be reduced to the following component form:

\begin{matrix} \frac{\partial p'}{\partial x'} = (2 α_{1} + α_{2}) \frac{d}{d x'} {(\frac{d V}{d x'})}^{2} \\ + 4 (γ_{3} + γ_{4} + γ_{5} + \frac{γ_{6}}{2}) \frac{d}{d x'} {(\frac{d V}{d x'})}^{4}, \end{matrix}

(8)

\begin{matrix} \frac{\partial p ″}{\partial y ″} = \frac{d}{d x'} [μ (θ) \frac{d V}{d x'}] + 6 (β_{2} + β_{3}) \frac{d^{2} V}{d x^{' 2}} (\frac{d V}{d x^{' 2}}) - σ B^{2} V + σ E B \\ + [{(ρ \bar{β})}_{0} (θ - θ_{m}) + {(ρ \bar{β})}_{1} {(θ - θ_{m})}^{2}] g, \end{matrix}}

(9)

\frac{\partial p'}{\partial z'} = 0

(10)

In equation (9), the viscoelastic fluid viscosity depends on temperature and the Reynolds exponential model is used which is defined as:

μ (θ) = μ_{0} \exp [- \bar{h} (θ - θ_{m})],

(11)

Here $μ_{0}$ indicates the dynamic viscosity, and $\bar{h}$ represents the strength dependency between $θ$ and $μ (θ)$ .

The modified pressure is stated as follows:

\begin{matrix} \bar{p} = p' - (2 α_{1} + α_{2}) {(\frac{d V}{d x'})}^{2} \\ - 4 (γ_{3} + γ_{4} + γ_{5} + \frac{γ_{6}}{2}) {(\frac{d V}{d x'})}^{4} . \end{matrix}

(12)

In view of equation (12), equations (8–10) can be reduced to the following form as:

\frac{\partial \bar{p}}{\partial x'} = 0,

(13)

\begin{array}{l} \frac{\partial \bar{p}}{\partial y'} = \frac{d}{d x ″} [μ (θ) \frac{d V}{d x ″}] + 6 (β_{2} + β_{3}) \frac{d^{2} V}{d X^{″ 2}} {(\frac{d V}{d X ″})}^{2} - σ B^{2} V + σ E B \\ + [{(ρ \bar{β})}_{0} (θ - θ_{m}) + {(ρ \bar{β})}_{1} {(θ - θ_{m})}^{2}] g, \end{array}}

(14)

\frac{\partial \bar{p}}{\partial z'} = 0 .

(15)

The preceding equation indicates that the pressure gradient $\partial \hat{p} / \partial y'$ is constant. We obtain:

\begin{array}{l} \frac{d}{d x ″} [μ (θ) \frac{d V}{d x ″}] + 6 (β_{2} + β_{3}) \frac{d^{2} V}{d x^{' 2}} {(\frac{d V}{d x'})}^{2} - σ B^{2} V + σ E B \\ + [{(ρ β)}_{0} (θ - θ_{m}) + {(ρ β)}_{1} {(θ - θ_{m})}^{2}] g = \frac{\partial \bar{p}}{\partial y'}, \end{array}}

(16)

The following are the velocity boundary conditions at the duct walls:

V (\pm d) = 0,

(17)

The energy equation including viscous dissipation and Joule heating is defined as:

\begin{matrix} k \frac{d^{2} θ}{d x'^{2}} + μ (θ) {(\frac{d V}{d x'^{2}})}^{2} + 2 (β_{2} + β_{3}) {(\frac{d V}{d x'})}^{4} \\ + σ (B^{2} V^{2} - 2 BEV + E^{2}) = 0 . \end{matrix}

(18)

The relevant thermal boundary conditions prescribed at the duct walls are as follows:

θ (- d) = θ_{1}, θ (+ d) = θ_{2},

(19)

The next step is to transform the above equations into dimensionless form using the scaling variables listed below:

v = \frac{V}{v_{0}}, x = \frac{x'}{d}, T = \frac{θ - θ_{m}}{θ_{1} - θ_{2}},

(20)

Here $v_{0}$ indicates the reference velocity. Using equation (20), the dimensionless form of equations (16) and (19) emerge as:

\begin{matrix} \frac{d}{d x} (\exp [- ε T] \frac{d v}{d x}) + 6 χ {(\frac{d v}{d x})}^{2} \frac{d^{2} v}{d x^{2}} - γ^{2} v \\ + η + λ_{r} (T + λ_{c} T^{2}) = Ω_{p}, \end{matrix}

(21)

\frac{d^{2} T}{d x^{2}} + δ \exp [- ε T] {(\frac{d v}{d x})}^{2} + 2 δ χ {(\frac{d v}{d x})}^{4} + δ γ^{2} v^{2} - ζ_{1} v + ζ_{2} = 0 .

(22)

The boundary conditions (17) and (19) now assume the dimensionless form:

\begin{matrix} v (- 1) = 0, v (+ 1) = 0, \\ T (- 1) = 0.5, T (+ 1) = - 0.5, \end{matrix}}

(23)

And

\begin{matrix} χ = \frac{v_{0}^{2} (β_{2} + β_{3})}{μ_{0} d^{2}}, Re = \frac{v_{0} d}{υ}, G r = \frac{g d^{3} {(ρ \bar{β})}_{0} (θ_{1} - θ_{2})}{υ μ_{0}}, \\ ε = \bar{h} (θ_{1} - θ_{2}), η = \frac{σ B E d^{2}}{μ_{0} v_{0}}, λ_{c} = \frac{{(ρ \bar{β})}_{0} (θ_{1} - θ_{2})}{{(ρ β)}_{1}}, \\ γ = B d \sqrt{\frac{σ}{μ_{0}}}, λ_{r} = \frac{G r}{Re}, Ω_{p} = \frac{d^{2}}{v_{0} μ_{0}} \frac{\partial \bar{p}}{\partial y'} δ = \frac{μ_{0} v_{0}^{2}}{k (θ_{1} - θ_{2})}, \\ ζ_{1} = \frac{2 B E d^{2} v_{0} σ}{(θ_{1} - θ_{2}) k}, ζ_{2} = \frac{d^{2} E^{2} σ}{(θ_{1} - θ_{2}) k}, υ = \frac{μ_{0}}{ρ} . \end{matrix}}

(24)

Here $γ$ is Hartmann number, Re is Reynolds number, Gr is the thermal Grashof number, $ε$ indicates the variable viscosity parameter, $Ω_{p}$ is the dimensionless pressure gradient, $η$ is electrical field strength parameter, $λ_{c}$ quadratic thermal convection parameter, $χ$ is fourth-grade fluid parameter, $λ_{r}$ is thermal buoyancy parameter, $δ$ is Brinkman number (which indicates the ratio of heat created by viscous dissipation to heat transferred by molecular conduction), $ζ_{1}$ is a parameter relating the effects of heat generation owing to the interaction of magnetic and electric fields on heat conduction, $ζ_{2}$ is ratio of Joule heating to heat conduction, and ν is kinematic viscosity of the fourth-grade fluid.

On the left plate $x' = - d$ , the skin friction coefficient and Nusselt number are expressed in dimensionless format following Rajagopal and Na²⁷:

\begin{matrix} S ~ {\exp [- ε T] \frac{dv}{dx} |}_{x = - 1} + {2 γ {(\frac{dv}{dx})}^{3} |}_{x = - 1}, \\ h ~ {- \frac{dT}{dx} |}_{x = - 1} . \end{matrix}}

(25)

Perturbation solutions using HPM

The transformed ordinary differential equations for momentum and energy, that is, equations (21) and (22) are nonlinear and coupled, rendering exact solutions intractable. As a result, we utilize the homotopy perturbation method (HPM) for determining series solutions. This method was introduced as a special case of the homotopy analysis method (HAM) and is very adaptive and yields excellent accuracy and stability for nonlinear differential equation systems. It has been utilized in a wide spectrum of complex engineering science problems, including structural vibration,⁵⁴ peristaltic pumping,^55,56 and more recently, electromagnetic duct flows containing porous media.⁴⁴ The perturbation formulation for the present problem is as follows:

\begin{array}{l} h (u, ζ) = (1 - ζ) [L_{0} (u) - L_{0} ({\bar{v}}_{0})] + ζ \\ [L_{0} (u) - ε Θ \frac{d^{2} u}{d x^{2}} - ε \frac{d Θ}{d x} \frac{d u}{d x} + 6 χ {(\frac{d u}{d x})}^{2} \frac{d^{2} u}{d x^{2}} - γ^{2} u + η + λ_{r} (Θ + λ_{c} Θ^{2}) - Ω_{p}], \end{array}

(26)

\begin{array}{l} h (Θ, ζ) = (1 - ζ) [L_{0} (Θ) - L_{0} ({\bar{T}}_{0})] + ζ \\ [L_{0} (Θ) + δ (1 - ε Θ) {(\frac{d u}{d x})}^{2} + 2 δ χ {(\frac{d u}{d x})}^{4} + (δ γ^{2} u^{2} - ζ_{1} u + ζ_{2})], \end{array}

(27)

Where $ζ$ is the embedding parameter. The linear operator $L_{0}$ is selected and the following tentative guesses, $u_{0}, Θ_{0}$ are made:

\begin{matrix} L_{0} = \frac{d^{2}}{d x^{2}}}, Linear operator \\ {\bar{v}}_{0} = \frac{1}{4} (x^{2} - 1), Θ_{0} = - \frac{x}{2}}, Initial guess \end{matrix}

(28)

Let us now propose the series expansions for equations (26) and (27), that is,

{\begin{matrix} u = u_{0} + ζ u_{1} + ζ^{2} u_{2} + \dots, \\ θ = θ_{0} + ζ θ_{1} + ζ^{2} θ_{2} + \dots . \end{matrix}

(29)

By inserting equation (29) into equations (26) and (27), we acquire a set of linear differential equations at every level. After the accomplishment of this stage, we will utilize the symbolic software Mathematica to determine the linear differential equation solutions for each order.

Zeroth other system $ζ^{0}$ and solutions

By substituting equation (29) into equations (26) and (27), we get the zeroth order system as follows:

\begin{matrix} L_{0} (u_{0}) - L_{0} ({\bar{v}}_{0}) = 0, \\ u_{0} (\pm 1) = 0, \end{matrix}}

(30)

\begin{matrix} L_{0} (Θ_{0}) - L_{0} ({\bar{T}}_{0}) = 0, \\ {\bar{T}}_{0} (- 1) = 0.5, {\bar{T}}_{0} (+ 1) = - 0.5, \end{matrix}}

(31)

The zeroth order solution is found to be:

\begin{matrix} u_{0} = \frac{x^{2} - 1}{4}, \\ Θ_{0} = - \frac{x}{2}, \end{matrix}}

(32)

First order system $ζ^{1}$ and solutions

The first order system is obtained in the following form:

\begin{matrix} L_{0} (u_{1}) + L_{0} ({\bar{v}}_{0}) - ε Θ_{0} \frac{d^{2} u_{0}}{d x^{2}} - ε \frac{d Θ_{0}}{d x} \frac{d u_{0}}{d x} + 6 χ {(\frac{d u_{0}}{d x})}^{2} \frac{d^{2} u_{0}}{d x^{2}} \\ - γ^{2} u_{0} + η + λ_{r} (Θ_{0} + λ_{c} Θ_{0}^{2}) - Ω_{p}, \\ u_{1} (\pm 1) = 0, \end{matrix}}

(33)

\begin{matrix} L_{0} (Θ_{1}) + L_{0} ({\bar{T}}_{0}) + δ (1 - ε Θ_{0}) {(\frac{d u_{0}}{d x})}^{2} + 2 δ χ {(\frac{d u_{0}}{d x})}^{4} \\ + (δ γ^{2} {u_{0}}^{2} - ζ_{1} u_{0} + ζ_{2}), \\ Θ_{1} (\pm 1) = 0, \end{matrix}}

(34)

The solution of the first order system is obtained as:

\begin{matrix} u_{1} = \frac{x^{2} - 1}{48} [H a^{2} (- 5 + x^{2}) - λ_{c} λ_{r} + 4 x (- ε + λ r) \\ - x^{2} (λ_{c} λ_{r} + 3 χ) - 3 (4 + 8 η + χ - 8 Ω_{p})], \end{matrix}

(35)

\begin{matrix} Θ_{1} = \frac{1 - x^{2}}{480} [γ^{2} (11 - 4 x^{2} + x^{4}) δ - 10 (- 5 + x^{2}) ζ_{1} \\ + 240 ζ_{2} + 2 δ {(1 + x^{2}) (5 + χ) + x^{4} χ}], \end{matrix}

(36)

Second order system $ζ^{2}$ and solutions

The second order system is obtained as:

\begin{matrix} L_{0} (u_{2}) - ε Θ_{1} \frac{d^{2} u_{0}}{d x^{2}} - ε Θ_{0} \frac{d^{2} u_{1}}{d x^{2}} - ε \frac{d Θ_{1}}{d x} \frac{d u_{0}}{d x} - ε \frac{d Θ_{0}}{d x} \frac{d u_{1}}{d x} \\ + 6 χ [2 \frac{d {\bar{v}}_{0}}{d x} \frac{d^{2} {\bar{v}}_{0}}{d x^{2}} \frac{d {\bar{v}}_{1}}{d x} + \frac{d^{2} {\bar{v}}_{1}}{d x^{2}} {(\frac{d {\bar{v}}_{0}}{d x})}^{2}] \\ - γ^{2} u_{1} + λ_{r} (Θ_{1} + 2 λ_{c} Θ_{0} Θ_{1}), \\ u_{2} (\pm 1) = 0, \end{matrix}}

(37)

\begin{matrix} L_{0} (Θ_{1}) + δ [2 \frac{d u_{1}}{d x} + 8 χ \frac{d u_{1}}{d x} {(\frac{d u_{0}}{d x})}^{2} - ε Θ_{1} - 2 ε Θ_{0} \frac{d u_{1}}{d x}] \\ \frac{d u_{0}}{d x} + 2 δ γ^{2} u_{1} u_{0} - ζ_{1} u_{1}, \\ Θ_{1} (\pm 1) = 0, \end{matrix}}

(38)

The solution of the second order system is obtained as:

\begin{matrix} u_{2} = u_{2, 1} + u_{2, 2} x^{1} + u_{2, 3} x^{2} + u_{2, 4} x^{3} + u_{2, 5} x^{4} + u_{2, 6} x^{5} \\ + u_{2, 7} x^{6} + u_{2, 8} x^{7} + u_{2, 9} x^{8} + u_{2, 10} x^{9}, \end{matrix}

(39)

\begin{matrix} Θ_{2} = Θ_{2, 1} + Θ_{2, 2} x^{1} + Θ_{2, 3} x^{2} + Θ_{2, 4} x^{3} + Θ_{2, 5} x^{4} \\ + Θ_{2, 6} x^{5} + Θ_{2, 7} x^{6} + Θ_{2, 8} x^{7} + Θ_{2, 9} x^{8} + Θ_{2, 10} x^{9} . \end{matrix}

(40)

The constants $u_{2, n}, Θ_{2, n}; n = 1, 2, \dots 10,$ derived in the preceding calculations can be found in the Appendix. We may deduce the final form of the outcomes employing the HPM condition as:

v = lim_{ζ \to 1} u = u_{0} + ζ u_{1} + ζ^{2} u_{2} + \dots .

(41)

θ = lim_{ζ \to 1} θ_{0} + ζ θ_{1} + ζ^{2} θ_{2} + \dots .

(42)

Results and discussion

In this section, we present detailed graphical and tabular solutions for velocity, temperature, skin friction, and the Nusselt number in the duct, for the influence of all key control parameters. All computations are executed in Mathematica and shown in Table 1 and Figures 2 to 15. As one particular case of our investigation, Table 2 displays the numerical comparison with previously published data. It is noticed that the current findings are precisely consistent with Abbasi et al.⁵⁷.

Table 1.

Comparison of the velocity profile with previously published results by fixing the remaining parameters as $ε = γ = η = λ_{r} = λ_{c} = 0, χ = 0.5, Ω_{p} = 1$ .

x	Present results	Abbasi et al.⁵⁷
0	−0.395353044732200	−0.395353045
0.1	−0.390278044901173	−0.390377559
0.2	−0.375724400825540	−0.375724401
0.3	−0.352088141541283	−0.352088141
0.4	−0.320513599352345	−0.320336095
0.5	−0.281320880996496	−0.281320881
0.6	−0.235797999318292	−0.235797999
0.7	−0.184408971096291	−0.184408971
0.8	−0.127691231510506	−0.127691231
0.9	−0.066095053134708	−6.61E-02
1	0	0

Figure 2.

Coefficient of skin friction for multiple values of variable viscosity parameter $ε$ , fourth-grade fluid parameter $χ$ , Brinkman number $δ$ , electric field strength parameter $η$ , Hartmann number $γ$ , thermal buoyancy parameter $λ_{r}$ , and quadratic thermal convection parameter $λ_{c}$ .

Figure 3.

Nusselt number for multiple values of variable viscosity parameter $ε$ , fourth-grade fluid parameter $χ$ , Brinkman number $δ$ , electric field strength parameter $η$ , Hartmann number $γ$ , thermal buoyancy parameter $λ_{r}$ , and quadratic thermal convection parameter $λ_{c}$ .

Figure 4.

Influence of fourth-grade material parameter $χ$ on velocity distribution.

Figure 5.

Influence of electric field strength parameter $η$ on velocity distribution.

Figure 6.

Influence of variable viscosity parameter $ε$ on velocity distribution.

Figure 7.

Influence of thermal buoyancy (natural) convection parameter $λ_{r}$ on velocity distribution.

Figure 8.

Effect of quadratic thermal convection parameter $λ_{c}$ on velocity distribution.

Figure 9.

Influence of Hartmann (magnetic) number $γ$ on velocity distribution.

Figure 10.

Effect of fourth-grade material fluid parameter $χ$ on temperature distribution.

Figure 11.

Influence of Brinkman number $δ$ on temperature distribution.

Figure 12.

Influence of electromagnetic heat generation to conduction ratio $ζ_{1}$ on temperature distribution.

Figure 13.

Influence of variable viscosity parameter $ε$ on temperature distribution.

Figure 14.

Effect of thermal buoyancy parameter $λ_{r}$ on temperature distribution.

Figure 15.

Effect of quadratic thermal convection parameter $λ_{c}$ on temperature distribution.

Table 2.

HPM results for the effects of key parameters on skin friction and Nusselt number.

Parameters	Values	Coefficient of skin friction	Nusselt number
$ε$	0	0.363579642	0.058035341
	0.4	0.346855352	0.064830629
	0.6	0.327962242	0.073336828
$χ$	0	0.335917849	0.070981940
	0.5	0.343657214	0.066632321
	1	0.348663956	0.063809635
$δ$	0.3	0.345738770	0.070702503
	0.5	0.346855352	0.064830629
	0.7	0.347993394	0.058866085
$η$	1	0.346855352	0.064830629
	1.5	0.776411238	0.110749135
	2	1.240973934	0.096432223
$γ$	0	0.363702482	0.069451489
	0.5	0.346855352	0.064830629
	1	0.306876751	0.053988768
$λ_{r}$	0.5	0.176197991	0.040934872
	1	0.346855352	0.064830629
	2	0.685748678	0.082449520
$λ_{c}$	0	0.285794633	0.057846555
	0.5	0.346855352	0.064830629
	1	0.407345992	0.070554752
$ζ_{1}$	0	0.362327892	0.025694391
	0.5	0.346855352	0.064830629
	1	0.333674090	0.099221419

Skin friction and Nusselt number

Table 2 shows skin friction coefficient and Nusselt number for all the emerging parameters.

A significant enhancement in skin friction accompanies an elevation in fourth-grade fluid parameter $χ$ , Brinkman number $δ$ , electric field parameter $η$ , thermal buoyancy parameter $λ_{r},$ and quadratic thermal convection parameter $λ_{c}$ . However, skin friction is strongly reduced with a rise in variable viscosity parameter $ε$ , Hartmann (magnetic) number $γ$ , and electromagnetic heat generation to conduction ratio, $ζ_{1}$ . Nusselt number magnitudes are elevated with increase in variable viscosity parameter $ε$ , electromagnetic heat generation to conduction ratio, $ζ_{1}$ thermal buoyancy parameter $λ_{r}$ , and quadratic thermal convection parameter $λ_{c}$ , whereas they are significantly suppressed with increment in fourth-grade fluid parameter $χ$ , Brinkman number $δ$ , and Hartmann magnetic number $γ .$ A more complex response is computed for increment in electric field parameter $η$ − initially the Nusselt number is boosted from $η$ = 0 (vanishing electrical field) to $η$ = 1, but a subsequent increase to $η$ = 1.5 actually decreases heat transfer to the wall of the duct.

The results In Table 2 are further visualized in Figures 2 and 3. Figure 2 shows the variation in the skin friction profile. It is noticed here that increment in the variable viscosity parameter $ε$ and the Hartmann magnetic field produce a consistent decrement in skin friction profile. Clearly a boost in exponential viscosity will enhance overall viscosity of the fourth-grade fluid and this will inhibit shearing along the wall, manifesting in a lower skin friction. A stronger magnetic field will also oppose momentum transfer and will decelerate the flow reducing the skin friction at the wall but simultaneously enhancing the momentum boundary layer thickness. The fourth-grade material parameter has a different influence on skin friction since it is related to stress relaxation and retardation, not merely viscosity. An increase in this parameter causes a gradual increase in skin friction implying a thinner momentum boundary layer. An increment followed by a weak increase is observed in skin friction with Brinkman number due to kinetic energy dissipation (internal friction). Increasing thermal buoyancy parameter strongly elevates the skin friction due to the intensity of natural convection currents which assist momentum development. This assists the flow and manifests in an elevation in velocity which implies faster shearing of the fluid against the duct walls and produces the associated hike in skin friction. Similarly increasing electrical field strength parameter and quadratic convection both enhance skin friction and produce thinner momentum boundary layers. Clearly the stronger electrical field intensity assists axial momentum development and is the primary reason for deploying such an orientation of the electrical field. The nonlinear effect of quadratic convection therefore avoids the under-prediction in skin friction associated with the classical Boussinesq (linear) model.

Figure 3 shows that Nusselt number is elevated with variable viscosity parameter but is suppressed with fourth-grade material parameter. A significant decrement in Nusselt number is also associated with greater Brinkman number since temperatures are elevated in the boundary layer and the net effect is to draw heat away from the wall, leading to a cooling there. A thinner thermal boundary layer thickness will result with higher Brinkman numbers associated with stronger viscous dissipation. Hartmann number elevation also induces a strong downturn in Nusselt number since the electromagnetic fluid is heated with higher magnetic field. Increasing the magnetic field in the fluid flow amplifies Lorentz force that dominates over the viscous force. This results in a net transfer of heat away from the wall into the core region of the duct and decreases Nusselt number. A similar though less prominent increment in Nusselt number accompanies an increase in thermal buoyancy and quadratic convection parameters. Again, the initial increase in Nusselt number is computed with a rise in electrical field strength parameter but subsequent increase in this parameter produces the opposite trend which is associated with the presence of Joule dissipation in the flow.

Velocity distribution

Figures 4 to 9 show the impact of all key parameters on velocity profiles. The velocity profile in Figure 4 exhibits a strong depletion with larger values of the fourth-grade fluid parameter. In practical terms, an increase in $χ$ increases the stress relaxation relative to retardation in a fourth-grade which counteracts momentum diffusion and decelerates the flow. This will increase hydrodynamic boundary layer thickness at the duct walls. Furthermore, when $χ = 0$ , the classical case of a Newtonian viscous fluid is retrieved. Figure 5 indicates the velocity profile is significantly enhanced with greater electric field intensity $(η)$ . When an electric field is introduced into a working fluid containing charged particles, the particles are attracted to the electrodes, causing the fluid to be induced in the same direction as the applied electric field. This technique is most often seen in microfluidic devices when a strong electric field is present in relatively small channels. However, it can also be used in larger nuclear reactor ducts to regulate momentum distribution. Strong acceleration is therefore produced with axial electrical field. Figure 6 depicts the effect of the varying viscosity parameter $ε$ on the velocity profile. It is evident that the velocity increases in the left half space but then decreases in the right half space with increment in viscosity parameter. The peak velocity arises consistently near the centerline of the duct. Asymmetrical parabolic profiles span the duct width. Figure 7 depicts the effect of the thermal buoyancy convection parameter $λ_{r}$ on the velocity profile. With increasing values of this parameter, the flow is accelerated across the channel width, that is, velocity profile is greatly enhanced owing to the enhanced contribution of thermal buoyancy relative to viscous forces. For $λ_{r}$ = 0 forced convection is present, and velocities are minimized. For $λ_{r}$ = 1 both thermal buoyancy and viscous force contribute equally. When $λ_{r}$ > 1, thermal buoyancy dominates the regime and significant momentum flux escalation is produced. It is also noteworthy that as viscosity varies with temperature and shear rate, the fluid becomes less dense in the presence of heat and this induces a boost to momentum resulting in a hike in velocity. Figure 8 demonstrates that acceleration is also produced in the axial flow with increasing quadratic convection parameter $λ_{c}$ . The nonlinear temperature dependence of the non-Boussinesq buoyancy term directly contributes to the elevation in velocity. Figure 9 depicts the impact of the Hartmann number $γ$ on the velocity field. As can be observed, the velocity profile decreases consistently owing to the significant retarding influence of the magnetic field. When a magnetic field is applied to a working fluid, it experiences a Lorentz force, which is a drag force that opposes fluid motion. Strong flow deceleration is therefore induced with higher Hartmann number across the entire span of the duct. Clearly the non-magnetic scenario ( $γ$ = 0), that is, vanishing magnetic field produces the maximum velocity. For $γ$ = both Lorentz magnetic drag force and viscous hydrodynamic force contribute equally. For $γ$ > 1 the Lorentz magnetic drag force dominates the viscous force leading to strong flow deceleration. Symmetric parabolic topologies are computed across the duct. The excellent flow regulation properties of transverse magnetic field are thereby confirmed.

Temperature distribution

Figures 10 to 15 depict the evolution of the temperature profiles with key parameters. Figure 10 demonstrates that the temperature is boosted with elevation in fourth-grade fluid $χ$ parameter rises. However, we can see that the implications are insignificant closer to the channel’s walls, with only a minor increase in the middle of the channel. The heating of the fluid is associated with stress relaxation and retardation effects. Thermal diffusion is encouraged with these phenomena. Figure 11 illustrates that enhancing the Brinkman number $δ$ causes the temperature profile to strongly increase, in particular, near the left wall as dissipation enhances the conversion of kinetic energy to heat via internal friction. As a result, the fourth-grade fluid also decelerates. The case of vanishing Brinkman number, $δ$ = 0 clearly corresponds to the minimal temperature and confirms that non-inclusion of viscous heating leads to an under-prediction in temperature in the bulk fluid. For this case the fourth order derivative term, $2 δ γ {(\frac{dv}{dx})}^{4}$ in equation (22) disappears. Figure 12 demonstrates the fact that the electromagnetic heat generation relative to thermal conduction parameter $ζ_{1}$ strongly suppresses temperature profile. Clearly when this parameter vanishes $ζ_{1}$ → 0, and thermal conduction dominates over the electrical and magnetic field heating contribution and maximum temperatures are computed. Figure 13 shows that the variable viscosity parameter produces a opposite behavior in temperature, closer to the duct centerline. The exponential viscosity variation exerts a much less dramatic effect on temperature than the fourth-grade material parameter. Figure 14 shows that there is a slight decrement in temperatures across the width of the duct with increment in thermal buoyancy (natural convection) parameter, $λ_{r}$ . As with the other temperature profiles, the weak temperature is always computed at the left hotter duct wall and decays monotonically to the right cooler wall. Thermal buoyancy therefore primarily influences the momentum field, and indeed features only in that equation. Via coupling with the energy (heat) equation there is an indirect effect on the temperature distribution. The temperature, however, is not as markedly influenced as the velocity field. Similarly in Figure 15 there is a relatively weak decrease in temperature across the duct width with greater values of quadratic convection parameter $λ_{r}$ . Again, the minor decrement is associated with the coupling effect between the momentum and energy equations (21) and (22), specifically via the term, $λ_{r} (T + λ_{c} T^{2})$ in equation (21). Generally smooth profiles are computed for both temperature and velocity distributions testifying to the excellent numerical stability of HPM.

Concluding remarks

A new theoretical study has been presented for fourth-grade viscoelastic EMHD natural convection flow in a vertical parallel plate duct with non-Boussinesq quadratic convection. The effects of viscous dissipation and Joule heating (Ohmic dissipation) and Reynolds exponential temperature-dependent viscosity model have been included, all of which arise in nuclear and solar hybrid duct transport phenomena. The normalized nonlinear ordinary differential boundary value problem has been solved with a homotopy perturbation method (HPM) up to second order approximation which achieves exceptional accuracy and stability in terms of power series solutions. Detailed solutions for axial velocity, temperature, skin friction, and Nusselt number have been depicted graphically and in tables for the effects of key electromagnetic, thermal, and dissipation parameters. The key findings of the current analysis may be summarized as follows:

(i) There is a considerable enhancement in wall skin friction with an increment in fourth-grade fluid parameter $χ$ , Brinkman number $δ$ , electrical field parameter $η$ , thermal buoyancy parameter $λ_{r}$ , and quadratic thermal convection parameter $λ_{c}$ . However, skin friction is strongly reduced with a rise in variable viscosity parameter $ε$ , Hartmann (magnetic) number $γ$ , and electromagnetic heat generation to conduction ratio, $ζ_{1}$ .

(ii) Nusselt number magnitudes are elevated with increase in variable viscosity parameter $ε$ and electromagnetic heat generation to conduction ratio, $ζ_{1}$ , thermal buoyancy parameter $λ_{r}$ , and quadratic thermal convection parameter $λ_{c},$ whereas they are significantly suppressed with increment in fourth-grade fluid parameter $χ$ , Brinkman number $δ$ , and Hartmann magnetic number $γ$ . The latter parameters therefore achieve better thermal control at the duct walls which is critical in regulating high temperatures (“thermal management”) in real nuclear reactor duct transport.

(iii) A variable response in Nusselt number is computed for increment in electric field parameter $η;$ initially the Nusselt number is boosted from $η$ = 0 (vanishing electrical field) to $η$ = 1, but a subsequent increase to $η$ = 1.5 actually decreases heat transfer to the wall of the duct.

(iv) Increasing fourth-grade viscoelastic parameter $χ$ elevates the viscosity and damps the axial velocity strongly as does higher magnetic Hartmann number.

(v) A minor increase in temperature is induced with greater viscoelastic parameter $χ$ , while a minor decrement is noticed for thermal buoyancy parameter $λ_{r}$ and quadratic thermal convection parameter $λ_{c}$ .

(vi) A marked boost in temperature is generated with increasing Brinkman number $δ$ whereas a very minor reduction is produced with a large increment in variable viscosity parameter.

(vii) Increasing electrical strength parameter $η$ produces significant flow acceleration.

(viii) In the left half space of the duct an increase in variable viscosity parameter $ε$ produces significant axial flow acceleration whereas in the right half space weak deceleration is produced.

(ix) Consistent acceleration, that is, higher velocities are associated with an increment in thermal buoyancy parameter $λ_{r}$ (ratio of thermal Grashof number and Reynolds number).

(x) A weaker but distinct increase in axial velocity is also observed with higher values of quadratic thermal convection parameter $λ_{c}$ .

(xi) The greatest elevation in skin friction number is produced with high electrical field strength parameter whereas the largest Nusselt number at the wall corresponds to low Brinkman number.

The current study has probed deeper into emerging smart electromagnetic liquids for nuclear and other energy systems duct transport phenomena. Attention has, however, been confined to steady state flows. Future work may consider time-dependent effects and additionally the use of other non-Newtonian microstructural models such as the polar and micropolar family of formulations. Furthermore, alternating electrical and magnetic fields may also be addressed to reflect further industrial scenarios.

Footnotes

Appendix

Handling Editor: Chenhui Liang

Author’s Note

Hakan F Öztop is also affiliated to Department of Mechanical and Nuclear Engineering, College of Engineering, University of Sharjah, United Arab Emirates.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Lijun Zhang and Muhammad Mubashir Bhatti are supported by the National Natural Science Foundation of China No. 12172199.

ORCID iD

Muhammad Mubashir Bhatti

References

Lorraine

Lorrain

Houle

. Magneto-fluid dynamics: fundamentals and case studies of natural phenomena. New York, NY: Springer, 2006.

Walker

Moon

, et al. Electromagnetic control of flow distribution between parallel ducts of liquid-metal-cooled blankets. Fus Technol 1994; 25: 398–410.

Liu

Shen

H-B

, et al. Applications of magneto-strictive, magneto-optical, magnetic fluid materials in optical fiber current sensors and optical fiber magnetic field sensors: a review. Opt Fiber Technol 2021; 65: 102634.

Dahmani

Muñoz-Cámara

Laouedj

, et al. Heat transfer enhancement of ferrofuid fow in a solar absorber tube under non-uniform magnetic field created by a periodic current-carrying wire. Sustain Energy Technol Assess 2022; 52: 101996.

Ikeda

Aoyama

Gotou

, et al. Development of analytical method for DC electromagnetic flow coupler including end effects. J Nucl Sci Tech 1987; 24: 988–998.

Murakami

Araseki

. Characteristic evaluation of electromagnetic flow couplers using a liquid metal MHD analysis code. Nucl Eng Des 2005; 235: 1503–1515.

Kim

Cha

Kim

, et al. DC magnetic field effect on a liquid sodium channel flow. Nucl Eng Des 2008; 238: 280–284.

Miyzaki

Inoue

Yamaoka

, et al. Magneto-hydro-dynamic pressure drop of lithium flow in rectangular ducts. Fus Technol 1986; 10: 830–836.

Castellanos

. Electrohydrodynamics. New York, NY: Springer, 1998.

10.

Kourmatzis

Shrimpton

. Electrohydrodynamic inter-electrode flow and liquid jet characteristics in charge injection atomizers. Exp Fluid 2014; 55: 1688.

11.

Dau

Dinh

Biu

, et al. Corona based airflow using parallel discharge electrodes. Exp Therm Fluid Sci 2016; 79: 52–53.

12.

Chang

Tsubone

Chung

, et al. Mechanism of electrohydrodynamically induced flow in a wire-non-parallel plate electrode type gas pump. J Electrost 2009; 67: 335–339.

13.

Zhao

Adamiak

. Numerical simulation of the electrohydrodynamic flow in a single wire-plate electrostatic precipitator. IEEE Trans Ind Appl 2008; 44: 683–691.

14.

Tanski

Kocik

Mizeraczyk

. Electrohydrodynamic gas pump with both insulated electrodes driven by dielectric barrier discharge. IEEE Trans Dielectr Electr Insul 2011; 18: 1429–1432.

15.

Gazaryan

Chirkov

. Numerical estimation of the performance of a flow-type electrohydrodynamic heat exchanger with the streamlined electrode configuration. J Electrostat 2019; 97: 31–36.

16.

Tanski

Kocik

Barbucha

, et al. A system for cooling electronic elements with an EHD coolant flow. J Phys Conf Ser 2014; 494: 12010.

17.

Kasayapanand

Kiatsiriroat

. EHD enhanced heat transfer in wavy channel. Int Commun Heat Mass Transfer 2005; 32: 809–821.

18.

Narla

Tripathi

Bhandari

, et al. Electrokinetic insect-bioinspired membrane pumping in a high aspect ratio bio-microfluidic system. Microfluid Nanofluid 2022; 26: 85.

19.

Carstoiu

. Fundamental equations of electromagnetodynamics of fluids: various consequences. Proc Natl Acad Sci USA 1968; 59: 326–331.

20.

Eringen

Maugin

. Electrodynamics of continua. New York, NY: Springer, 1990.

21.

Jian

Chang

. Electromagnetohydrodynamic (EMHD) micropumps under a spatially non-uniform magnetic field. AIP Adv 2015; 5: 057121.

22.

Chakraborty

Paul

. Microchannel flow control through a combined electromagnetohydrodynamic transport. J Phys D Appl Phys 2006; 39: 5364.

23.

Umavathi

Bég

. Double diffusive convection in a dissipative electrically conducting nanofluid under orthogonal electrical and magnetic field: a numerical study. Nanosci Technol 2021; 12: 59–90.

24.

Soundalgekar

Aranakey

. Effects of couple stresses in MHD channel flow. Nucl Eng Des 1977; 44: 301–308.

25.

Muzychka

Edge

. Laminar non-Newtonian fluid flow in noncircular ducts and microchannels. ASME J Fluids Eng 2008; 130: 111201.

26.

Chaudhuri

Sinha

Chakraborty

, et al. Thermal characteristics of forced convection in combined pressure and shear-driven flow of a non-Newtonian third-grade fluid through parallel plates. Heat Transfer 2021; 50: 6737–6756.

27.

Rajagopal

. Natural convection flow of a non-Newtonian fluid between two vertical flat plates. Acta Mech 1985; 54: 239–246.

28.

Kazemi

Javanmard

Taheri

, et al. Heat transfer investigation of the fourth-grade non-Newtonian MHD fluid flow in a plane duct considering the viscous dissipation, Joule heating and forced convection on the walls. SN Appl Sci 2020; 2: 1752.

29.

Fatima

Nazeer

Lashin

MMA

, et al. Developments of electro-osmotic two-phase flows of fourth-grade fluid through convergent and divergent channels. Mathematics 2023; 11: 1832.

30.

Nadeem

Ali

. Analytical solutions for pipe flow of a fourth-grade fluid with Reynold and Vogel’s models of viscosities. Commun Nonlinear Sci Numer Simul 2009; 14: 2073–2090.

31.

Sajid

Hayat

Asghar

. On the analytic solution of the steady flow of a fourth-grade fluid. Phys Lett A 2006; 355: 18–26.

32.

Salawu

Fatunmbi

Ayanshola

. On the diffusion reaction of fourth-grade hydromagnetic fluid flow and thermal criticality in a plane Couette medium. Results Eng 2020; 8: 100169.

33.

Kunnegowda

Mahanthesh

Lorenzini

, et al. Signifcance of induced magnetic field and exponential space dependent heat source on quadratic convective flow of Casson fluid in a micro-channel via HPM. Math Model Eng Probl 2019; 6: 369–384.

34.

Triveni

Mahanthesh

. Optimization and sensitivity analysis of heat transport of hybrid nanoliquid in an annulus with quadratic Boussinesq approximation and quadratic thermal radiation. Eur Phys J Plus 2020; 135: 459–522.

35.

Kumar

Sood

. Effect of quadratic density variation on mixed convection stagnation point heat transfer and MHD fluid flow in porous medium towards a permeable shrinking sheet. J Porous Media 2016; 19: 1083–1097.

36.

Okoya

Hassan

Salawu

. On free convection flow of a moving vertical permeable plate with quadratic Boussinesq approximation and variable thermal conductivity. Heat Transf Res 2021; 52: 55–66.

37.

Al-Kouz

Mahanthesh

Alqarni

, et al. A study of quadratic thermal radiation and quadratic convection on viscoelastic material flow with two different heat source modulations. Int Comm Heat Mass Transfer 2021; 126: 105364.

38.

Sabu

Mackolil

Mahanthesh

, et al. Nanoparticle aggregation kinematics on the quadratic convective magnetohydrodynamic flow of nanomaterial past an inclined flat plate with sensitivity analysis. Proc Inst Mech Eng Part E J Process Mech Eng 2022; 236: 1056–1066.

39.

Bhatti

Shahid

Sarris

, et al. Spectral relaxation computation of Maxwell fluid flow from a stretching surface with quadratic convection and non-Fourier heat flux using Lie symmetry transformations. Int J Mod Phys B 2023; 37: 2350082.

40.

Mahanthesh

Mackolil

Radhika

, et al. Significance of quadratic thermal radiation and quadratic convection on boundary layer two-phase flow of a dusty nano-liquid past a vertical plate. Int Commun Heat Mass Transf 2021; 120: 105029.

41.

Fasheng

Bhatti

Bég

, et al. Numerical study of natural convection dissipative electro-magnetic non-Newtonian flow through a non-Darcy channel. ZAMM 2022; 102: e202100608.

42.

Escandón

Santiago

Bautista

, et al. Hydrodynamics and thermal analysis of a mixed electromagnetohydrodynamic-pressure driven flow for Phan–Thien–Tanner fluids in a microchannel. Int J Therm Sci 2014; 86: 246–257.

43.

Siva

Jangili

Kumbhakar

, et al. Unsteady electromagnetohydrodynamic flow of couple stress fluid through a microchannel: a theoretical analysis. Eur J Mech B Fluids 2022; 95: 83–93.

44.

Bhatti

Bég

Ellahi

, et al. Natural convection non-Newtonian EMHD dissipative flow through a microchannel containing a non-Darcy porous medium: homotopy perturbation method study. Qual Theory Dyn Syst 2022; 21: 97.

45.

Gangadhar

Rupa Lavanya

Chamkha

. EMHD flow of second-grade fluid through a heated permeable disk with space dependent heat source. Proc Inst Mech Eng Part E J Process Mech Eng 2023: 095440892311592. DOI: 10.1177/09544089231159203.

46.

Chakraborty

Dey

Chakraborty

. Thermal characteristics of electromagneto-hydrodynamic flows in narrow channels with viscous dissipation and Joule heating under constant wall heat flux. Int J Heat Mass Transfer 2020; 67: 1151–1162.

47.

Huang

Lai

. Effect of Joule heating on EHD-enhanced natural convection. AIAA J Thermophys Heat Transfer 2006; 20: 939–945.

48.

Kalpana

Madhura

Kudenatti

. Impact of temperature-dependant viscosity and thermal conductivity on MHD boundary layer flow of two-phase dusty fluid through permeable medium. Eng Sci Technol Int J 2019; 22: 416–427.

49.

Ahmad

. Flow control of non-Newtonian fluid using Riga plate: Reiner-Phillipoff and Powell-Eyring viscosity models. J Appl Fluid Mech 2019; 12: 127–133.

50.

Hussain

Akbar

Sarwar

, et al. Effect of time dependent viscosity and radiation efficacy on a non-Newtonian fluid flow. Heliyon 2019; 5: e01203.

51.

Akbar

Tripathi

Khan

, et al. A numerical study of magnetohydrodynamic transport of nanofluids from a vertical stretching sheet with exponential temperature-dependent viscosity and buoyancy effects. Chem Phys Lett 2016; 661: 20–30.

52.

Prakash

Siva

Tripathi

, et al. Peristaltic pumping of magnetic nanofluids with thermal radiation and temperature-dependent viscosity effects: modelling a solar magneto-biomimetic nanopump. Renew Energy 2019; 133: 1308–1326.

53.

Sobamowo

Akinshilo

. Analysis of flow, heat transfer and entropy generation in a pipe conveying fourth-grade fluid with temperature dependent viscosities and internal heat generation. J Mol Liq 2017; 241: 188–198.

54.

Rashidi

Shooshtari

Bég

. Homotopy perturbation study of nonlinear vibration of von Karman rectangular plates. Comput Struct 2012; 106/107: 46–55.

55.

Alarabi

Elsayed

Bég

. Homotopy perturbation method for heat transfer in peristaltic flow of viscoelastic fluid in an eccentric cylinder with variable viscosity effects. Life Sci J 2014; 11: 197–206.

56.

Tripathi

Bég

. A numerical study of oscillating peristaltic flow of generalized Maxwell viscoelastic fluids through a porous medium. Transp Porous Media 2012; 95: 337–348.

57.

Abbasi

Ahmadian

Ppetroudi

, et al. Analysis of a fourth-grade fluid flow in a channel by application of VIM and HAM. Indian J Sci Res 2014; 1: 389–395.

Computation of non-Newtonian quadratic convection in electro-magneto-hydrodynamic (EMHD) duct flow with temperature-dependent viscosity

Abstract

Keywords

Introduction

Mathematical model for EMHD non-Newtonian quadratic dissipative duct flow

Perturbation solutions using HPM

Zeroth other system ζ 0 and solutions

First order system ζ 1 and solutions

Second order system ζ 2 and solutions

Results and discussion

Skin friction and Nusselt number

Velocity distribution

Temperature distribution

Concluding remarks

Footnotes

Appendix

Author’s Note

Declaration of conflicting interests

Funding

ORCID iD

References

Zeroth other system $ζ^{0}$ and solutions

First order system $ζ^{1}$ and solutions

Second order system $ζ^{2}$ and solutions