# Analysis of the Physical Behavior of the Periodic Mixed-Convection Flow around a Nonconducting Horizontal Circular Cylinder Embedded in a Porous Medium

Authors: Muhammad Ashraf; Zia Ullah; Saqib Zia; Sayer O. Alharbi; Dumitru Baleanu; Ilyas Khan
Journal: Journal of Mathematics
 (2021)

---

## Abstract

An oscillatory mixed-convection fluid flow mechanism across a nonconducting horizontal circular cylinder embedded in a porous medium has been computed. For this purpose, a model in the form of partial differential equations is formulated, and then, the governing equations of the dimensionless model are transformed into the primitive form for integration by using primitive variable formulation. The impact of emerging parameters such as porous medium parameterΩ, Richardson number λ, magnetic force parameter ξ, and Prandtl number Pr on skin friction, heat transfer, and current density is interpreted graphically. It is demonstrated that accurate numerical results can be obtained by the present method by treating nonoscillating and oscillating parts of coupled partial differential equations simultaneously. In this study, it is well established that the transient convective heat transfer, skin friction, and current density depend on amplitude and phase angle. One of the objects of the present study is to predict the mechanism of heat and fluid flow around different angles of a nonconducting horizontal circular cylinder embedded in a porous medium.


---

## Body

## 1. Introduction

The mixed convection mechanism and thermally driven flow in a porous medium have important applications in utilizing geothermal energy, industrial process, and mechanical, civil, and chemical engineering. Moreover, in insulation materials, the flow of helium in pebble-bed nuclear reactors, underground disposal of nuclear or nonnuclear waste, food processing and storage, metallurgy, crude oil extraction, flow in the eye of glaucoma patients, and flow through filtering media are important examples of fluid flows in the porous medium. Porous medium effects cause resistance to the fluid flow due to the difference between the viscosity of a porous medium and the fluid. In the view of these applications, Joshi and Gebhart [1] presented a technical note on the mixed-convection mechanism through heat flux surface with a porous medium by using Hamming’s predictor-corrector method. The heat and fluid flow mechanism over a flat plate which is placed in a porous medium has been performed in Hossain et al.’s study [2]. They used Keller box technique to obtain numerical results with free-stream velocity and small amplitude variations in the wall temperature. The effect of transverse magnetic field on mixed-convection phenomena over a semi-infinite permeable vertical plate in a porous medium at constant heat flux has been studied by Chamkha [3]. Later, Chamkha et al. [4] proposed a nonsimilar boundary-layer problem for the natural convection flow around an isothermal sphere placed in a saturated porous medium.The mixed-convection boundary-layer flow considering Cu,Al2O3, and TiO2 as a nanoparticles in the base fluid along the heat and cooled cylinder has been investigated in Nazar et al.’s study [5]. Gorla et al. [6] discussed a boundary layer phenomenon of the mixed convection flow of the nanofluid past a vertical wedge that is placed in a porous medium theoretically. Cimpean and Pop [7] studied the effect of a porous medium on the fully developed mixed-convection flow along a channel for the saturated nanofluid. Imran et al. [8] presented an analysis of mixed-convection unsteady flow over a stretching vertical sheet placed in a porous medium with a heat source. Sumaily and Thomson [9] studied a steady and pulsating forced-convection flow around a circular cylinder placed in an empty or porous medium and concluded that heat transfer enhances due to the porous medium. Hayat et al. [10] analyzed a mixed-convection flow mechanism over a porous sheet in the presence of convective boundary conditions with variable thermal conductivity effects. They examined similar effects of velocity and temperature profiles in the qualitative sense. Darcy and porous particle diameter-affected heat transfer rate in the cavity has been computed by Behzadi et al. in [11].Ullah et al. [12] considered a magnetohydrodynamic mixed-convection flow of Casson fluid past a wedge by assuming that the wedge is moving inside a porous medium with the effects of convective boundary conditions and chemical reaction numerically. Yadav and Singh [13] considered a transient flow mechanism around a circular cylinder that was filled with porous material. They observed that the velocity gradient decreases slowly from the boundary to the central region of the cylinder for large values of the viscosity ratio. Ashraf and Fatima [14] discussed the effect of viscous dissipation on the transient flow to obtain numerical solutions around different stations of the sphere. Garg et al. [15] examined the convective boundary-condition outcomes on the unsteady MHD flow over a porous plate in the slip-flow regime. Ashraf et al. [16, 17] performed the natural convection mechanism on different shapes numerically. Rashad et al. [18] investigated the heat transfer problem by means of micropolar nanofluids past a cylinder in the presence of a saturated porous medium. They obtained oscillatory behavior of microrotation velocity for large values of the material parameter Δ.Rashad et al. [19] illustrated the heat transfer behavior of the mixed convection flow of a water-copper nanofluid around different positions of the trapezoidal enclosure that is consigned in a porous medium with variable properties. The variable density effect on an oscillatory fluid flow mechanism around a nonconducting shape has been investigated by Ashraf and Ullah in [20], and then they performed magnetothermo analysis in [21]. Recently, Ullah et al. [22] computed an oscillatory flow model across a nonconducting shape embedded in a thermally stratified medium. They depicted the considerable fluctuating response in heat transport for large value of Prandtl number.Taking idea from the current literature review, it is illustrated that the oscillatory mixed-convection fluid flow around a nonconducting horizontal circular cylinder embedded in a porous medium has not yet been computed in the literature. By following the current literature [2, 18], we establish a mathematical model in terms of coupled nonlinear partial differential equation to highlight the physical behavior of the oscillatory mixed-convection flow across a circular cylinder embedded in a porous medium. Later, we will solve this model by using the very accurate numerical scheme finite difference method. We explore the effects of different parameters included in the flow model on skin friction τw, heat transfer qw, and current density jw.


## 2. Governing Equations and Problem Formulation

Considering a two-dimensional, unsteady, viscous, and incompressible boundary-layer fluid flow phenomenon across a nonconducting horizontal circular cylinder embedded in a porous medium, as shown in Figure1, here, the coordinate system (x, y) represents the distance along the surface and normal to the surface, respectively. We assume the magnetic field coordinate system (Hx, Hy) along the surface and normal to the surface, respectively, where T is the temperature variable throughout the assumed domain. We establish the concept of oscillatory by assuming external fluid velocity as u=Uτ, far from the surface of the domain. Moreover, to justify our claim of given phenomena, we assume the effect of the magnetic field from outreach to the surface of the cylinder. By following the current literature [18], we extend this work to highlight the physical behavior of the oscillatory mixed-convection flow across a circular cylinder embedded in a porous medium. In keeping view of the constraint of the length of the article, we assume the set of the coupled nonlinear dimensionless system of partial differential equations as given as follows:(1)∂u¯∂x¯+∂v¯∂y¯=0,(2)∂u¯∂τ+u¯∂u¯∂x¯+v¯∂u¯∂y¯=∂2u¯∂y¯2+ξh¯x∂h¯x∂y¯+h¯y∂h¯x∂y¯+Ω1−u¯+λθ¯sinα,(3)∂h¯x∂x¯+∂h¯y∂y¯=0,(4)∂h¯∂τ+u¯∂h¯∂x¯+v¯∂h¯∂y¯−h¯x∂u¯∂x¯−h¯y∂u¯∂y¯=1γ∂2h¯∂y¯2,(5)∂θ¯∂τ+u¯∂θ¯∂x¯+v¯∂θ¯∂y¯=1Pr∂2θ¯∂y¯2.Figure 1
Geometry of the flow process.The dimensionalized selected boundary conditions are given as follows:(6)u¯=v¯=0,h¯y=h¯x=0,θ¯=1aty¯=0,u¯⟶U¯τ,θ¯⟶0,h¯y⟶1asy¯⟶∞.In equations (1)–(6), ξ is the magnetic force parameter, the Richardson number is λ, the magnetic Prandtl number is γ, the Prandtl parameter is Pr, Ω is the permeability parameter (the porous medium parameter), DaL is the Darcy number, ReL is Reynold’s number, θ is the dimensionless fluid temperature, and ε+ represents the porosity of the medium. We take the stream velocity equation as Uτ=1+εeiωτandε≪1, in which ε is the oscillating condition of small magnitude and the parameter ω is the frequency parameter. The velocity, magnetic field, and temperature components are u,v,hx,hy, and θ in the sum form of nonoscillating and oscillating parts given as follows:(7)u¯=us+εuteiωτ,v¯=vs+εvteiωτ,h¯x=hxs+εhxteiωτ,h¯y=hys+εhyteiωτ,θ¯=θs+εθteiωτ.We substitute (7) into equations (1)–(5) along with boundary conditions (6). Then, by following Ashraf and Ullah [20], we separate nonoscillating and oscillating components of the order analysis subject to Oε0 and Oεeiωτ. Later, by following Ullah et al.’s work [21], we use primitive variable formulation to obtain nonoscillating and oscillating parts of the primitive coupled nonlinear PDEs.


## 3. Computational Technique

As discussed in the preceding paragraph, the obtained form of the primitive system of the nonoscillating and oscillating differential equation is solved numerically by adopting the finite difference method. First, we solve the nonoscillating part, and then, obtained results are used in the oscillating part to overcome the nonlinearity involved in the system of equations. From the difference form of each equation, we obtained the system of algebraic equations. This system of algebraic equation has coefficient tridiagonal matrix and unknown variablesU,V,θ, and φ. In keeping view of the structure of the coefficient matrix, we use the Gaussian elimination technique to find the values of unknown variables. Once we secure values of these variables, we use these values to determine the oscillatory form of skin friction τw, heat transfer qw, and current density jw with the help of the following expression:(8)τw=∂U∂Yy=0+εAscosωt+αs,qw=∂θ∂Yy=0+εAtcosωt+αt,jw=∂φ∂Yy=0+εAmcosωt+αm,As=u12+u221/2,At=θ12+θ221/2,Am=φx12+φx221/2,αs=tan−1u2u1,αt=tan−1θ2θ1,αm=tan−1φx2φx1,where u1andu2 are real and imaginary parts of velocity, θ1 and θ2are the real and imaginary parts of temperature, and φx1 and φx2 are the real and imaginary parts of the magnetic field, respectively. Moreover, As, At, and Am are amplitudes, while αs, αt, and αm are phase angles in the abovementioned mechanism.


## 4. Results and Discussion

The current work is addressed on an oscillatory mixed-convection fluid flow mechanism for the viscous and incompressible fluid across a nonconducting horizontal circular cylinder embedded in a porous medium. Here, we precisely thrashed out the substantial effect of emerging parameters graphically within the prescribed boundary conditions. We examine the effect of all parameters on skin frictionτw, heat transfer qw, and current density jw across three different stations of the cylinder. Here, we accord the numerical solution for the effects of all parameters engaged in the flow model gradually.Figures2(a)–2(c) are presented to illustrate the oscillating behavior of τw, qw, and jw for three selected values of Ω=0.1,0.5, and 1.0 at α=0.5,1.0, and π position for Pr = 7.0. In Figure 2(a), the small amplitude of skin friction with distinct variations at each station is obtained, respectively. Figure 2(b) depicts that an enhancement in oscillation of heat transfer with a similar trend is obtained for lower Ω=0.1 at each angle. The fluctuating behavior in current density is uniformly distributed with high amplitude at each station α=0.5,1.0,andπ in the presence of Ω for Pr = 7.0 in Figure 2(c). Increasing Ω means the medium is more porous and the fluid permeability in the porous layer is increased and thus yields the resistance in the fluid flow. However, due to the strong Richardson number λ which acts like a pressure gradient and dominated over the resistance, skin friction is increased and slight changes in qw and jw are noted. Figure 3(a) shows the effects of the magnetic force parameter ξ on all abovementioned physical quantities. In this plot, it is noted that an increment in ξ causes reduction in skin friction. Figure 3(b) describes the higher oscillation in heat transfer with the same variation in ξ as observed in Figure 3(a). With increasing values of ξ that enhanced the Lorentz forces which opposes the flow and generates the current in fluid layers. Thus, Lorentz forces increase the friction between the layers of the fluid which yields an increment in qw and jw in Figure 3(c). In Figure 4(a), we claim that the increase in Richardson number λ has same oscillating effects on τw at each position. The maximum oscillations in qw and jw at the position with good variations are examined in Figures 4(b) and 4(c). The large value of Richardson number λ acts like a pressure gradient to increase the motion of the fluid that causes an increase in qw and jw. The impact of different values of Prandtl number Pr is highlighted for oscillating skin friction in Figure 5(a). The fluctuations in τw are deliberated in the prominent form, and good variations are obtained at α=1.0 and π station. The minimum amplitude in heat transfer for lower Pr =1.0 is noted, but good amplitude with certain height is plotted for large Pr =7.0 at each station in Figure 5(b). The current density with prominent amplitude at each position is depicted in Figure 5(c). Also, it is highlighted that the fluctuation of qw increases as Pr is increased, but mild variation in current density is observed. The most favorable position is α=1.0 for all oscillating results in the abovesaid mechanism.Figure 2
The physical profiles for (a)τw, (b) qw, and (c) jw at positions α=0.5,1.0, and π with three values of the porous medium parameter Ω=0.1,0.5, and 1.0 where others are γ=0.1, Pr =7.0, ξ=1.0, and λ=3.5.
(a)(b)(c)Figure 3
The physical profiles for (a)τw, (b) qw, and (c) jw at positions α=0.5,1.0, and π with three values of ξ=1.0,5.0, and 10.0 where others are γ=0.01, Pr=7.0, Ω=0.1, and λ=6.0.
(a)(b)(c)Figure 4
The physical profiles for (a)τw, (b) qw, and (c) jw at positions α=0.5,1.0, and π with three values of λ=0.5,1.5, and 2.5 where others are γ=0.01, Pr=0.1, ξ=0.2, and Ω=0.1.
(a)(b)(c)Figure 5
The physical profiles for (a)τw, (b) qw, and (c) jw at positions α=0.5,1.0, and π with three values of Pr =1.0,4.0, and 7.0 where others are γ=0.01, λ=2.0, ξ=0.1, and Ω=0.1.
(a)(b)(c)


## 5. Conclusions

The current work is based on an oscillatory mixed-convection heat transfer mechanism for the viscous and incompressible electrically conducting fluid across a nonconducting horizontal circular cylinder embedded in a porous medium. The dimensionless model is converted into the primitive form by applying primitive variable formulation. The final primitive formed equations are solved with the help of the implicit finite difference method. The oscillatory skin frictionτw, heat transfer qw, and current density jw are displayed around different α=0.5,1.0, and π positions of the cylinder for various physical parameters graphically.(1
)The small amplitude of oscillation in skin friction is noted at each station for each physical parameter, but prominent fluctuations are examined in heat transfer and current density(2)
The oscillating behavior in current density is uniformly distributed at each stationα=0.5,1.0,andπ for each value of the porous medium parameter and Prandtl number(3)
The prominent oscillating response in heat transfer and current density is noted with good variations for each value of the mixed convection parameterλ at each station(4)
The large value of the parameterλ acts like a pressure gradient to increase the motion of the fluid which caused an increase in qw and jw and dominance over porosity(5)
The maximum fluctuation in heat transfer is obtained for large values of Prandtl number Pr=4.0 and 7.0 at each position of the nonconducting cylinder(6)
The most favorable position for prominent amplitude and maximum value of heat transfer and current density isα=1.0 radian



---

*Source: 8839146-2021-01-25.xml*