Next Article in Journal
On the Determination of the 3D Velocity Field in Terms of Conserved Variables in a Compressible Ocean
Next Article in Special Issue
Adsorption of Precursors on Substrates in the Presence of scCO2 for the Synthesis of Supported Metallic Nanoparticles: Experiments and Modeling
Previous Article in Journal
Detection and Recognition of the Underwater Object with Designated Features Using the Technical Stereo Vision System
Previous Article in Special Issue
A Meshless Algorithm for Modeling the Gas-Dynamic Interaction between High-Inertia Particles and a Shock Layer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Control of MHD Flow and Heat Transfer of a Micropolar Fluid through Porous Media in a Horizontal Channel

by
Miloš Kocić
*,
Živojin Stamenković
,
Jelena Petrović
and
Jasmina Bogdanović-Jovanović
Faculty of Mechanical Engineering, University of Niš, A. Medvedeva 14, 18000 Niš, Serbia
*
Author to whom correspondence should be addressed.
Fluids 2023, 8(3), 93; https://doi.org/10.3390/fluids8030093
Submission received: 1 February 2023 / Revised: 27 February 2023 / Accepted: 2 March 2023 / Published: 8 March 2023
(This article belongs to the Special Issue Focus on Supercritical Fluids: Control and Extraction)

Abstract

:
The problem considered in this paper is a steady micropolar fluid flow in porous media between two plates. This model can be used to describe the flow of some types of fluids with microstructures, such as human and animal blood, muddy water, colloidal fluids, lubricants and chemical suspensions. Fluid flow is a consequence of the constant pressure gradient along the flow, while two parallel plates are fixed and have different constant temperatures during the fluid flow. Perpendicular to the flow, an external magnetic field is applied. General equations of the problem are reduced to ordinary differential equations and solved in the closed form. Solutions for velocity, microrotation and temperature are used to explain the influence of the external magnetic field (Hartmann number), the characteristics of the micropolar fluid (coupling and spin gradient viscosity parameter) and the characteristics of the porous medium (porous parameter) using graphs. The results obtained in the paper show that the increase in the additional viscosity of micropolar fluids emphasizes the microrotation vector. Moreover, the analysis of the effect of the porosity parameter shows how the permeability of a porous medium can influence the fluid flow and heat transfer of a micropolar fluid. Finally, it is shown that the influence of the external magnetic field reduces the characteristics of micropolar fluids and tends to reduce the velocity field and make it uniform along the cross-section of the channel.

1. Introduction

Magnetohydrodynamics studies the flow of an electrically conductive fluid in the presence of a magnetic field, either external or induced in the fluid during its movement. Such flows have wide application, from the steel casting industry, through heat exchangers in nuclear fusion reactors, to bioplasma, medicine and nanotechnologies. Electrically conductive fluids are used in a large number of industrial processes and technological procedures. Liquid metals, such as aluminum, mercury or crucible steel, are good examples of such fluids. The interaction between a moving fluid and a magnetic field causes a variety of MHD phenomena that can be employed in different ways. Although certain pioneering experiments in magnetohydrodynamics were conducted as early as the 19th century, it was not before the 1960s that the use of magnetic fields in the metallurgical industry was truly considered. In continuous steel casting, the flow in jets that fill molds is first controlled, and then the fluid is steadied while solidifying so as to yield a high-quality product. Another example is the construction of semiconductor crystals, where a high level of purity is achieved through electromagnetic control. Electromagnetic flow control in technological processes is also used to stabilize melts and their free surfaces, and to produce very fine powders, semiconductors, aluminum and superalloys with exquisite characteristics.
The Increased interest in studying MHD phenomena is today related to the development of fusion reactors, where plasma is controlled by a high-intensity magnetic field. Recent research has shown that electromagnetohydrodynamic (EMHD) flows may present a very reliable means to transport low-conductivity fluids in microsystems. A number of fluids can be transported in such microfluid devices with different goals. For example, the velocity of one fluid can be increased in its direct interaction with another, more mobile fluid, or the problem of the flow of two fluids can be applied in heat transfer devices, or one can perform the controlled mixing of fluids. Most of the scientific community associates the beginning of magnetohydrodynamics with Hartmann and Lazarus [1,2], who conducted the first experiments with liquid metals in 1937. Julius Hartmann also worked on different technical processes and tried to develop an electromagnetic conduction pump to transport electrically conductive fluids. After Hartmann and Lazarus, Bo Lehnert [3] studied a case similar to the Couette flow of a conductive viscous fluid between two parallel plates in the presence of a perpendicular magnetic field. He was the first to show that a magnetic field and the resulting induced currents significantly alter the velocity profile, and that boundary layers with high velocity gradients are formed in the vicinity of plates. Vast possibilities for the application and development of magnetohydrodynamics are provided by the diversity of electrically conductive fluids. One group of electrically conductive fluids that drew the attention of researchers in the second half of the 20th century was micropolar fluids.
Physically, micropolar fluids represent fluids that contain solid, randomly oriented (mainly spherical) particles dissolved in a viscous fluid, with the deformation of these particles being ignored. The Navier–Stokes model of classical hydrodynamics is largely limited—it cannot explain fluids that contain microstructures. To precisely describe the behavior of such fluids, one needs a theory that takes into consideration the internal movement of individual material particles. One of the best formulated theories of fluids with microstructures is the theory of micropolar fluids introduced by A.C. Eringen [4].
The theory of micropolar fluids that includes the influence of the presence of microstructures and their movement (translation and rotation) on fluid flow was first defined by Eringen [5]. The first results of his research were published in 1964 and 1965 in scientific journals and conference proceedings, while the theory of micropolar fluids was finally presented to the public in 1966 in the paper titled “Theory of Micropolar Fluids” [4]. Later, Eringen [6] (1972) generalized the theory of micropolar fluids by including heat transfer. Before Eringen, there were other theories that attempted to describe the flow of fluids with microstructures; however, not one was successful. Willson [7] (1970) was the first to introduce the boundary layer problem in micropolar fluids. In the same year, Peddiseon and McNitt [8] (1970) derived the boundary layer equations for a micropolar fluid, which are important for a large number of technical processes, and applied these equations to the steady flow around the front stagnation point, steady flow past an infinite plate and flow initiated in the boundary layer on an infinite plate.
Ariman [9,10] (1973, 1974) studied microcontinuum fluid mechanics, as well as a special case of micropolar fluid flow. On the other hand, Nath [11] (1975, 1976) later developed similar (automodel) and approximative solutions for an unsteady laminar boundary layer of micropolar fluids, while Ahmadi [12] (1976) also considered only similar solutions for an incompressible boundary layer of a micropolar fluid around a semi-infinite plate, by using the Runge–Kutta method with Newton’s iteration. Recently, there have been studies that consider effects of the thermal radiation of a micropolar fluid flow [13,14]. Many studies have used the model of micropolar fluids to describe the behavior of biological fluids (human and animal blood). Since Bayliss’ early research [15] (1952), several theoretical and experimental studies have been performed in order to understand the characteristics of blood flow [16,17,18]. All this research has been performed under the assumption that blood acts as a Newtonian fluid. However, it is now widely accepted that blood behaves as a non-Newtonian fluid under certain circumstances, particularly when subjected to low shear stress. The theory of micropolar fluids was first applied as a model for human blood by Turk [19,20] (1973, 1974), followed by Hogen et al. [21] (1989). Previous experimental research [22,23,24] into blood flow has shown that in given flow conditions, blood flow deviates to a certain extent from the flow of Newtonian fluids. More recently, fluid flow through channels has attracted a lot of attention, bearing in mind its application in the cardiovascular system [25], treating blood as a micropolar fluid. Ariman [20] examined blood flow in a circular rigid tube and determined that the micropolar fluid is a better model for blood flow since it also takes into account the microrotation of particles suspended in blood, which was later confirmed by the experimental results of Bugliarello and Sevilla [22].
A porous medium is a heterogeneous system comprising a solid matrix (skeleton, support) with its holes (pores), of relatively small dimensions compared with the characteristic dimensions of the matrix, filled with a fluid. This medium can be treated as a continuum by adequately observing the role of each phase (solid and fluid) during the transport through this system of phases. Due to the strong presence of porous media in nature, the research into them began early, first experimentally and then theoretically. In 1855 and 1856, Henry Darcy led a group of researchers who experimentally examined the flow of water through sand. On the basis of these experiments, they obtained the result that the velocity of water flow through sand is proportional to the pressure gradient, while the proportionality coefficient is sand permeability, whose value depends on the fluid characteristics (viscosity, temperature, etc.). This result was published in [26] and became known as Darcy’s law. Although Darcy’s law was initially developed to describe flow through sand, it has since been generalized for various situations and is in wide use today.
Research into fluid flow and heat transfer in porous media under the influence of an external magnetic field began in the 1980s. This research was certainly prompted by MHD flows and heat transfer in porous media that occur in numerous engineering processes, such as compact heat exchangers, metallurgy, casting, liquid metal filtration, nuclear reactor cooling, fusion control, etc. The problem of micropolar fluid flow through porous media has many practical applications [27], such as the flow through porous rocks, foam and foamy/spongy bodies, as well as alloys, polymer mixtures and microemulsions. For this reason, many researchers have recently studied unsteady free convective micropolar fluid flow in the presence or absence of a magnetic field. For example, Srinivasacharya [28] (2001) considered the effects of microrotation on an unsteady micropolar fluid flow between two parallel porous plates, with periodic suction. Bhargara [29] (2003) obtained numerical solutions for a free convective MHD micropolar fluid flow between two vertical parallel porous plates, by using the quasi-linearization method. Finally, Zueco et al. [30] (2009) analyzed an unsteady hydromagnetic micropolar fluid flow between two vertical and porous walls, by using the Network Simulation Method. The conclusion of all three studies is that the magnetic field is an important parameter in controlling the velocity of heat transfer in many applications of magnetohydrodynamics.
The research into blood flow in human and animal organs is also relevant. In 2003, Khaled and Vafai [31] examined the flow and heat transfer in biological tissues. Their main concept was the transfer in porous media using mass diffusion and different models of convective flows, such as Darcy’s and Brinkman’s models. They also analyzed energy transfer in tissues. They determined that the theory of porous media in biological tissues is the most efficient one since it contains fewer assumptions in relation to other biothermal models. In 2012, Mehmod et al. [32] examined the characteristics of an unsteady two-dimensional blood flow in a diseased porous arterial segment with flexible walls. The porous medium represented a lumen that contained blood clots, cholesterol and fatty platelets. It was noticed that a decrease in permeability led to a significant deceleration of flow, while the pressure drop and wall stress increased.
All the above studies imply the importance and comprehensiveness of the application of magnetohydrodynamic micropolar fluid flow in porous media. Therefore, this paper deals with the influence of magnetic fields and micropolar fluid characteristics on the nature of the flow of such fluids in porous media.

2. Physical and Mathematical Model

The model of micropolar fluid flow can be used to explain the flow of human or animal blood [33]. Furthermore, for example, there are numerous different approaches where porous media are used to represent lungs [34,35]. Thence, MHD micropolar fluid flow in porous media can be used to describe and take into consideration the numerous flows of some types of fluids with microstructures, such as human and animal blood, muddy water, colloidal fluids, lubricants and chemical suspensions.
The basic equations governing the motion of micropolar fluid flow through porous media are [36]
ρ t + d i v ρ v = 0
ρ D v D t = ρ f p + j × B + μ v + ε + μ λ d i v v + μ + λ Δ v + 2 λ r o t ω
ρ D v D t = ρ f p + j × B + μ v + ε + μ λ d i v v + μ + λ Δ v + 2 λ r o t ω
ρ c p T t + v T = k 2 T + Φ + j 2 σ + ρ g v
where ρ —fluid density, v —vector of fluid velocity, f —vector of volume forces, p —pressure field, j —density of current, B —magnetic field vector, ε —dielectric constant, μ —dynamic viscosity, α , β , γ , λ —additional viscosities for micropolar fluids, ω —vector of microrotation, c p —mass heat capacity, T —temperature field and σ —electrical conductivity.
As can be seen from the basic equations for the flow of micropolar fluids (1–4), exact solutions for the Navier–Stokes equations are very difficult to obtain due to the nonlinearity that occurs in the partial differential equations [37]. On the other hand, most of the real fluid flows are of a non-Newtonian nature, which highlights the importance of this research. In addition, an important characteristic of micropolar fluids should be emphasized here. As we are dealing with fluids that contain microstructures (randomly oriented particles that can rotate around their axes), the Cauchy stress tensor is not symmetric. The presence of coupled stresses is a consequence of taking into account the rotational degrees of freedom for the elementary volume of the micropolar fluid [38]. This once again highlights the importance of considering the flow of micropolar fluids.
The physical model of the considered problem consists of two parallel and fixed plates extending in the x and z directions, with constant and different temperatures T w 1 and T w 2 , while, between them, there is a micropolar fluid flow in a porous medium (Figure 1). Steady fluid flow is a consequence of the constant pressure gradient along the x axis. Perpendicular to this fluid flow, in the direction of the y axis, an external magnetic field, of strength B , is applied. Flow analysis is usually performed assuming the constant fluid electrical conductivity and treating the problem as one-dimensional, i.e., a fully developed flow takes place in the x direction. The problem considered in the paper is given in the inductionless approximation and, due to a small difference in the temperature of the plates, buoyance force is neglected. Taking these approximations into consideration, the vectors of fluid velocity v and external magnetic field B are
v = u i ,
B = B j .
The governing equations for the flow, under the above-stated assumptions, can now be written as
μ + λ d 2 u d y 2 + λ d ω d y σ B 2 u d p d x μ χ u = 0
γ d 2 ω d y 2 λ d u d y 2 λ ω = 0
k d 2 T d y 2 + μ + λ d u d y 2 + σ B 2 u 2 + μ χ u 2 = 0
From the previous system of Equations (7)–(9), it can be noted that Equations (7) and (8) need to be solved together, i.e., as coupled equations.
Appropriate fluid and thermal boundary conditions for this problem are represented by expressions (10):
u = 0 , ω = 0 , T = T w 2   f o r   y = 0 u = 0 , ω = 0 , T = T w 1   f o r   y = h
The no-slip conditions require that the fluid velocities are equal to the plate velocities, that the boundary temperature conditions are isothermal conditions and that there is no microrotation at the plates.
To solve the system of equations for the observed problem, it needs to be transformed into the non-dimensional form. Therefore, the following transformations are introduced:
y * = y h , u * = u U , ω * = ω ω 0 , ω 0 = U h , U = h 2 P μ , P = p x = c o n s t , θ = T T w 2 T w 1 T w 2 , H a = B h σ μ , P r = μ c p k , E c = U 2 c p T w 1 T w 2 , K = λ μ , Γ = γ μ h 2 , R = h 2 χ .
The introduced non-dimensional parameters are well known in the theory of MHD micropolar fluid flow. The given parameters are H a —Hartmann number, P r —Prandtl number, E c —Eckert number, K —coupling parameter, Γ —spin gradient viscosity parameter, R —porosity parameter ( χ —porous medium permeability).
Equations (7)–(9), having been transformed into the non-dimensional form, have the following form:
1 + K d 2 u * d y * 2 + K d ω d y H a 2 + R u + 1 = 0
Γ d 2 ω * d y * 2 K d u * d y 2 K ω * = 0
d 2 θ d y * 2 + 1 + K P r E c d u * d y * 2 + P r E c H a 2 + R u * 2 = 0
The boundary dimensionless conditions for the previous equations are
u * = 0 , ω * = 0 , θ = 0   f o r   y * = 0 u * = 0 , ω * = 0 , θ = 1   f o r   y * = 1
As mentioned before, impulse Equation (12) and equation of angular impulse (13) are coupled. Therefore, after some mathematical transformations and operations, the fourth-order differential equation for non-dimensional velocity is obtained:
d 4 u * d y * 4 a d 2 u * d y * 2 + b u * d = 0
In Equation (16), the newly introduced constants a , b and d , are defined as
a = B * + E A D * , b = E B * , d = E C ,
where
A = K 1 + K , B * = H a 2 + R 1 + K , C = 1 1 + K , D * = K Γ , E = 2 K Γ .
Finding a solution for the fourth-order differential Equation (16) is not an easy task. In this case, there are three possible solutions, which depend on the roots of the characteristic equation.
An appropriate particular solution of Equation (16) is
u p * = d b = C B * ,
Meanwhile, the characteristic equation of the homogeneous part of Equation (15) is given in the form of
δ 4 a δ 2 + b = 0 ,
or
f o r   δ 2 = r r 2 a r + b = 0 .
The roots of characteristic Equation (20) are
r 1 = 1 2 a + a 2 4 b , r 2 = 1 2 a a 2 4 b
The mathematical model applied in the paper provides the possibility of predicting the fluid behavior based on the root of the characteristic equation for velocity (16). The exact analytical solution, obtained in this way, enables the analysis of the influence of all significant parameters without the need to repeat the numerical solution of partial differential equations for each case. Furthermore, the solution can be used to check the numerical solution of some more complex partial differential equations.
The existence and uniqueness of the solution are analyzed when analyzing the roots of the characteristic equation of the partial differential Equation (16). Depending on the value of a 2 4 b (it can be higher, lower or equal to zero), characteristic Equation (20), and therefore Equation (16), has three possible solutions.
For example, as expression a 2 4 b depends on the physical properties of the fluid ( μ , λ , γ , σ ) , the strength of the external magnetic field B , the permeability of the porous media χ and the distance between the plates h , for different values of these variables, expression a 2 4 b can be positive, negative or equal to zero.
Let us consider a problem of certain fluid flow in porous media between plates [36]. The physical properties of the fluid are
μ = 0.002 P a · s ; λ = 2 × 10 4 P a · s ; γ = 8 × 10 9 N · s ; σ = 1 S / m
where the strength of the external magnetic field is B = 10 [ T ] , the porous media permeability is χ = 1 × 10 5 m 2 and the distance between the plates is h = 0.01 m . For the set values, the numerical value of expression a 2 4 b becomes
a 2 4 b = 10.33
However, if we change the strength of the external magnetic field from B = 10 [ T ] to B = 5 [ T ] , the value of expression a 2 4 b becomes
a 2 4 b = 21.43
or, for the value of B = 8.99 [ T ] ,
a 2 4 b = 0 .
Based on the previous discussion, three possible solutions for non-dimensional velocity are given in the following equations:
-
in the case a 2 4 b > 0 ,
u * = C 1 e x p δ 1 y * + C 2 e x p δ 2 y * + C 3 e x p δ 3 y * + C 4 e x p δ 4 y * + d b ,
-
in the case a 2 4 b = 0 ,
u * = C 5 + C 6 y * e x p ξ 1 y * + C 7 + C 8 y * e x p ξ 2 y * + d b ,
-
in the case a 2 4 b < 0 ,
u * = C 9 c o s β 1 y * + C 10 s i n β 1 y * e x p α 1 y * + C 11 c o s β 1 y * + C 12 s i n β 1 y * e x p α 1 y * + d b .
Now, three possible solutions for microrotation are given, respectively:
ω * = C 1 D 1 e x p δ 1 y * + C 2 D 2 e x p δ 2 y * + C 3 D 3 e x p δ 3 y * + C 4 D 4 e x p δ 4 y * ,
ω * = E 13 + E 14 y * e x p ξ 1 y * + E 15 + E 16 y * e x p ξ 2 y * , ω * = P 3 * s i n β 1 y * + P 4 * c o s β 1 y * e x p α 1 y * +
P 5 * s i n β 1 y * + P 6 * c o s β 1 y * e x p α 1 y *
as well as for non-dimensional temperatures:
θ = P r E c [ 1 4 δ 1 2 C 1 e x p 2 δ 1 y * + 1 4 δ 2 2 C 2 e x p 2 δ 2 y * + 1 4 δ 3 2 C 3 e x p 2 δ 3 y * + 1 4 δ 4 2 C 4 e x p 2 δ 4 y * + D 2 δ 1 + δ 3 2 e x p ( ( δ 1 + δ 3 ) y * ) + D 3 δ 1 + δ 4 2 e x p δ 1 + δ 4 y * + E 1 δ 2 + δ 3 2 e x p δ 2 + δ 3 y * + E 2 δ 2 + δ 4 2 e x p δ 2 + δ 4 y * + ( 1 δ 1 2 F e x p δ 1 y * 1 + 1 δ 2 2 F e x p δ 2 y * 2 + 1 δ 3 2 F e x p δ 3 y * 3 + 1 δ 4 2 F e x p δ 4 y * 4 ) + 1 2 D y * 2 1 + 1 2 E y * 2 3 + 1 2 F y * 2 5 + H 1 y * 1 + H 2 1 ] ,
θ = PrEc [ ( ( Ω 28 + Ω 29 y * + Ω 30 y * 2 ) exp 2 ξ 1 y * + ( Ω 31 + Ω 32 y * + Ω 33 y * 2 ) exp 2 ξ 2 y * + Ω 36 * + Ω 37 * y * exp ξ 2 y 2 + Ω 38 * y * 2 + Ω 39 * y * 3 + Ω 40 * y * 4 + H 2 1 y * + H 2 2 ,
θ = PrEc { [ ( 1 2 α 1 Ω 45 + 1 2 χ 1 Ω 47 χ 2 Ω 49 cos 2 β 1 y * + 1 2 χ 2 Ω 47 χ 1 Ω 49 sin 2 α 1 y * exp 2 α 1 y * + 1 2 α 1 Ω 46 1 2 χ 1 Ω 48 + χ 2 Ω 50 cos 2 β 1 y * + 1 2 χ 2 Ω 48 χ 1 Ω 50 sin 2 β 1 y * exp 2 α 1 y * 1 2 β 1 Ω 51 sin 2 β 1 y * 1 2 β 1 Ω 52 cos 2 β 1 y * + Ω 53 χ 1 Ω 55 χ 2 cos β 1 y * + Ω 53 χ 2 + Ω 55 χ 1 sin β 1 y * exp α 1 y * + Ω 54 χ 1 Ω 56 χ 2 cos β 1 y * Ω 54 χ 2 Ω 56 χ 1 sin β 1 y * exp α 1 y * + 1 2 Ω 57 y * 2 + H 1 3 y * + H 2 3 .
Together with the expressions for velocity, microrotation and temperature, the shear stress at the plates and the dimensionless heat transfer coefficient—the Nusselt number on the plates—are given as important characteristics of the flow and heat transfer.
It is very interesting to note that for a fluid in a microcontinuum (couple stress fluids, micropolar fluids, polar fluids, dipolar fluids, etc.), the stress tensor is not symmetric. The shear stress in the considered problem is defined by the following expression:
τ = μ + λ d u * d y * + λ ω
-
in the case a 2 4 b > 0 ,
τ = ( μ + λ ) ( C 1 δ 1 exp ( δ 1 y * ) + C 2 δ 2 exp ( δ 2 y * ) + C 3 δ 3 exp ( δ 3 y * ) + C 4 δ 4 exp ( δ 4 y * ) ) + λ ( C 1 D 1 exp ( δ 1 y * ) + C 2 D 2 exp ( δ 2 y * ) + C 3 D 3 exp ( δ 3 y * ) + C 4 D 4 exp ( δ 4 y * )
-
in the case a 2 4 b = 0 ,
τ = ( μ + λ ) [ C 6 exp ( ξ 1 y * ) + ξ 1 ( C 5 + C 6 y * ) exp ( ξ 1 y * ) + C 8 exp ( ξ 2 y * ) + ( C 7 + C 8 y * ) exp ( ξ 2 y * ) ] + λ [ ( E 13 + E 14 y * ) exp ( ξ 1 y * ) + ( E 15 + E 16 y * ) exp ( ξ 2 y * ) ]
-
in the case a 2 4 b < 0 ,
τ = ( μ + λ ) { [ ( α 1 C 10 β 1 C 9 ) sin β 1 y * + β 1 C 10 + α 1 C 9 cos β 1 y * exp α 1 y * + β 1 C 12 α 1 C 11 cos β 1 y α 1 C 12 + β 1 C 11 sin β 1 y * exp α 1 y + λ P 3 * sin β 1 y * + P 4 * cos β 1 y * exp α 1 y * + P 5 * sin β 1 y * + P 6 * cos β 1 y * exp α 1 y *
As the Nusselt number represents the ratio of convective to conductive heat transfer at a boundary in a fluid, it is very important to understand heat transfer on the plates. Hence, the Nusselt number on the plates in the considered problem is
N u = d θ d y *
-
in the case a 2 4 b > 0 ,
N u = PrEc 1 2 δ 1 Cexp 1 2 δ 1 y * + 1 2 δ 2 Cexp 2 2 δ 2 y * + 1 2 δ 3 Cexp 3 2 δ 3 y * + 1 2 δ 4 Cexp 4 2 δ 4 y * + D 1 y * D 2 δ 1 + δ 3 exp δ 1 + δ 3 y * + D 3 δ 1 + δ 4 exp δ 1 + δ 4 y * + E 1 δ 2 + δ 3 exp δ 2 + δ 3 y * + E 2 δ 2 + δ 4 exp δ 2 + δ 4 y * + + E 3 y * + δ 1 Fexp 1 δ 1 y * + 1 δ 2 Fexp 2 δ 2 y * + 1 δ 3 Fexp 3 δ 3 y * + 1 δ 4 Fexp 4 δ 4 y * + F 5 y * + H 1 1
-
in the case a 2 4 b = 0 ,
N u = PrEc [ Ω 15 exp 2 ξ 1 y * + Ω 16 y * exp 2 ξ 1 y * + Ω 17 y * 2 exp 2 ξ 1 y * + Ω 18 exp 2 ξ 2 y * + Ω 19 y * exp 2 ξ 2 y * + Ω 20 y * 2 exp 2 ξ 2 y * + Ω 21 exp ξ 1 + ξ 2 y * + Ω 22 y * exp ξ 1 + ξ 2 y * + Ω 23 y * 2 exp ξ 1 + ξ 2 y * + Ω 24 exp ξ 1 y * + Ω 25 y * exp ξ 1 y * + Ω 26 exp ξ 2
-
in the case a 2 4 b < 0 ,
N u = PrEc [ Ω 45 exp 2 α 1 y * Ω 46 exp 2 α 1 y * + Ω 47 cos 2 β 1 y * exp 2 α 1 y * + Ω 48 cos 2 β 1 y * exp 2 α 1 y * + Ω 49 sin 2 β 1 y * exp 2 α 1 y * + Ω 50 sin 2 β 1 y * exp 2 α 1 y * Ω 51 cos 2 β 1 y * + Ω 52 sin 2 β 1 y * + Ω 53 cos β 1 y * exp α 1 y * Ω 54 cos β 1 y * exp α 1 y * + Ω 55 sin β 1 y * exp α 1 y * + Ω 56 sin β 1 y * exp α 1 y * + Ω 57 y * + 3 H 1
The constants that occur in Equations (22)–(38) are given in the Supplementary Materials.

3. Results and Discussion

In this section, the results for velocity, microrotation and temperature, as well as for shear stress and the Nusselt number, are presented graphically. The influence of the external magnetic field, the characteristics of the micropolar fluid and the characteristics of the porous medium are investigated and presented using the Hartmann number ( H a ), coupling ( K ) and spin gradient viscosity parameter ( Γ ), as well as the porosity parameter ( R ). Throughout this investigation and based on the obtained results, some basic conclusions on the flow of a micropolar fluid in porous media can be made. The values of the dimensionless parameters were chosen in such a way that the graphs would show the tendency of changes in physical quantities, with the change in the values of the given parameters.
The first three figures show the influence of the additional viscosity gamma ( γ ), through the spin viscosity parameter Γ (expressions 11), on velocity, microrotation and temperature.
From Figure 2, it can be noted that the velocity decreases with an increase in the spin viscosity parameter Γ from the value of Γ = 0.05 to Γ = 5 . However, at the same time, with a further increase in the spin gradient viscosity parameter ( Γ = 50 ), the velocity starts to increase. This means that the velocity function has reached its extreme for a value of the spin viscosity parameter between Γ = 5 and Γ = 50 . Figure 3, which represents the influence of the spin viscosity parameter Γ on the vector of microrotation, shows that a decrease in the absolute value of microrotation occurs in the entire height of the channel with an increase in the spin viscosity parameter Γ . Both figures show that with an increase in the additional viscosity gamma ( γ ), the motion characteristics of the micropolar fluid start to decrease. This is the expected phenomenon of the viscous braking of the fluid flow, which is a consequence of the increase in the total viscosity of the fluid.
An increase in the additional viscosity gamma ( γ ) leads to a decrease in temperature in the middle of the channel (Figure 4). As a consequence of a decrease in velocity and microrotation with an increase in the spin gradient viscosity parameter Γ , Joule heating in the middle of the channel decreases. An increase in temperature near the plates, on the other hand, happens due to the viscous heat.
Together with the additional viscosity gamma ( γ ), one more viscosity is taken into consideration when micropolar fluids are investigated. The influence of the additional viscosity lambda ( λ ) is presented through the coupling parameter K . The next three figures show the velocity, microrotation and temperature dependency on the coupling parameter K .
From Figure 5 and Figure 6, it is obvious that the additional viscosity lambda emphasizes the characteristics of micropolar fluids [39]. The vector of microrotation, defined with angular impulse Equation (13), is the landmark of micropolar fluids. Because of this, from Figure 6, it can be noted that an increase in the coupling parameter K leads to an increase in the absolute value of microrotation. Figure 5 shows that the velocity decreases with an increase in the coupling parameter K . This means that the additional viscosity lambda ( λ ) emphasizes the microrotation vector of small particles embedded in the fluid.
Figure 7 shows the influence of the coupling parameter K on the dimensionless temperature. With an increase in the coupling parameter K , conduction becomes the dominant form of heat transfer. This is also the consequence of the velocity decrease with an increase in the additional viscosity lambda ( λ ), which has already been explained in Figure 5.
The next five graphs (from Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12) present the influence of the external magnetic field through the Hartmann number (Ha) on velocity, microrotation, temperature, shear stress and the Nusselt number.
From Figure 8 and Figure 9, it is obvious that the velocity and microrotation decrease in absolute value as the Hartmann number increases [40,41]. This means that the motion characteristics of the micropolar fluid become lower in the entire width between the plates with an increase in the external magnetic field. This is something that is expected, as, with an increase in the external magnetic field in electrically conducting fluids, the Lorentz force acts in the direction opposite to the fluid flow and tends to retard the flow. It is very interesting to consider this through the roots of characteristic Equation (20), given by expression (21). In the previous section, it is determined that the value of expression a 2 4 b depends on the physical properties of the fluid ( μ , λ , γ , σ ) , the strength of the external magnetic field B , the permeability of the porous media χ and the distance between the plates h . It is shown that with a change in the intensity of the external magnetic field B , there is a change in the value of expression a 2 4 b (it becomes positive, negative or equal to zero), which means that the tendency of the velocity or microrotation vector, for certain fluids, can be determined simply by knowing the change in the value of expression a 2 4 b .
Figure 10 presents the influence of the external magnetic field on the dimensionless temperature. As the Hartmann number increases, conduction becomes the dominant form of heat transfer between the plates. For lower values of the Hartmann number, Joule heating increases the temperature of the fluid in the middle of the channel. The increase in temperature near the plates is the consequence of viscous heating.
The next figure, Figure 11, shows the behavior of the Nusselt number on the lower and upper plates for different values of the Hartmann number (Ha) and coupling parameter (K). For the lower as well as for the upper plate, the Nusselt number decreases in absolute value with an increase in the Hartmann number, which is expected in the case of the Hartmann flow, as heat transfer is mainly in the form of conduction. This is the consequence of the dominant influence of the magnetic field. On the other hand, for larger values of the coupling parameter (K), heat transfer on the plates is lower, and this is due to an increase in the additional viscosity lambda ( λ ), which is characteristic for micropolar fluids.
The influence of the external magnetic field and additional viscosity lambda on shear stress is presented in Figure 12 [42]. The increase in the Hartmann number causes a decrease in shear stress in absolute value, both on the lower and the upper plates. As expected, for higher values of the Hartmann number, the Lorentz force decreases the velocity in the entire width of the channel (Figure 8), and because of this, the shear stress becomes lower on the plates. In the other case, with an increase in the additional viscosity lambda ( λ ), which implies an increase in the coupling parameter (K), based on expressions (11), shear stress increases as the total viscosity of the fluid increases (the interpretation based on Equation (31)).
Finally, the last group of figures display the influence of the porosity parameter R on micropolar fluid flow and heat transfer. It is well known that as the value of the porosity parameter increases, the holes of the porous medium become smaller. The first two figures, Figure 13 and Figure 14, show the influence of the porosity parameter on the motion characteristics of the micropolar fluid, i.e., on the velocity and microrotation vectors.
As already mentioned, with an increase in the porosity parameter R , the fluid flow space becomes smaller so the micropolar fluid flow starts to decelerate; see Figure 13. Velocity decreases with an increase in the porosity parameter R , and the microrotation vector decreases as well in absolute value in the entire height between the plates. This is something that is expected and has been shown by many authors [43,44]. This means that a physically smaller space causes the braking of the fluid flow in a porous medium.
The next figure, Figure 15, presents the influence of the porosity parameter R on the shear stress on the plates.
From Figure 15, it is obvious that shear stress, both on the lower and the upper plates, decreases in absolute value with an increase in the porosity parameter. As already shown in the previous two figures, with an increase in the porosity parameter, the micropolar fluid flow decelerates, which leads to a decrease in shear stress on the plates. On the other hand, with an increase in the coupling parameter K , shear stress on the plates increases in absolute value on both plates. An increase in the coupling parameter K means that the additional viscosity lambda ( λ ) is increasing, and further that the total viscosity of the micropolar fluid becomes higher. Because of this, and according to expression (31), shear stress on the plates increases with an increase in the coupling parameter K . This is due to the friction of the fluid against the plates with an increase in the overall viscosity of the fluid.
The last two figures, Figure 16 and Figure 17, provide answers to the influence of the porosity parameter on heat transfer.
From Figure 16, it can be noted that for smaller values of the porosity parameter R , viscous heating increases the temperature near the plates, while Joule heating increases it in the middle of the channel. However, for larger values of the porosity parameter ( R = 50 ), Joule heating in the middle of the channel becomes dominant, as viscous heating near the plates starts to decrease.
The same can be seen in Figure 17, where, with an increase in the porosity parameter R , heat transfer on the plates decreases in absolute value, both on the upper and the lower plates. Regarding the influence of the coupling parameter K , with an increase in this parameter, heat transfer on the plates becomes smaller as the total viscosity of the micropolar fluid increases.

4. Conclusions

The object of consideration in this paper is the flow of a micropolar fluid in porous media. First, the physical model of the considered problem is defined. It is a steady micropolar fluid flow in porous media between two plates, where the plates have been kept at constant and different temperatures T w 1 and T w 2 , and an external magnetic field, of strength B , is applied perpendicularly to the fluid flow. The general equations for the considered problem are reduced to ordinary differential equations and solved together with appropriate fluid and thermal boundary conditions. Closed-form solutions for velocity, microrotation and dimensionless temperature are obtained in three possible cases, depending on the roots of characteristic equations. Numerical results are then presented in the form of graphs. These figures show the influence of the external magnetic field (Hartmann number), characteristics of the micropolar fluid (coupling and spin gradient viscosity parameter) and characteristics of the porous medium (porous parameter) on the flow and heat transfer characteristics of micropolar fluid flow in porous media. Some general conclusions can be made based on the discussion given in this paper. First, the external magnetic field reduces the characteristics of micropolar fluids and tends to retard the fluid flow as well. Second, an increase in the additional viscosity lambda emphasizes the vector of microrotation, as one of the main characteristics of micropolar fluids. Finally, the porosity parameter is used to show how the permeability of a porous medium can influence the fluid flow and heat transfer of a micropolar fluid. One of the ideas in this research is its application to the flow of biological fluids in porous environments, which are a natural environment for them. It is particularly interesting that such research is applied to blood flow in artificial tissues or membranes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/fluids8030093/s1.

Author Contributions

Methodology, M.K. and Ž.S.; Validation, J.P. and J.B.-J.; Writing – original draft, M.K.; Project administration, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia (Contract No. 451-03-47/2023-01/200109).

Data Availability Statement

All necessary data are in included in Supplementary Materials, as well as in the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hartmann, J.; Hg-Dynamics, I. Theory of the laminar flow of an electrically conductive liquid in a homogeneous magnetic field. Det Kgl. Dan. Vid. Sels. Mat.-Fys. Medd. 1937, 15, 1–27. [Google Scholar]
  2. Hartmann, J.; Lazarus, F. Hg-Dynamics II. Experimental investigations on the flow of mercury in a homogeneous magnetic field. Det Kgl. Dan. Vid. Sels. Mat.-Fys. Medd. 1937, 15, 1–45. [Google Scholar]
  3. Lehnert, B. On the behavior of an electrically conductive liquid in a magnetic field. Ark. Fys. 1952, 5, 69–90. [Google Scholar]
  4. Eringen, A.C. Theory of micropolar fluids. J. Math. Mech. 1966, 16, 1–18. [Google Scholar] [CrossRef]
  5. Eringen, A.C. Simple micro fluids. Int. J. Eng. Sci. 1964, 2, 205–217. [Google Scholar] [CrossRef]
  6. Eringen, A.C. Theory of thermomicrofluids. J. Math. Anal. Appl. 1972, 38, 480–496. [Google Scholar] [CrossRef] [Green Version]
  7. Willson, A.J. Boundary Layers in Micropolar liquids. Math. Proc. Camb. Philos. Soc. 1970, 67, 469–476. [Google Scholar] [CrossRef]
  8. Peddieson, J.; McNitt, R.P. Boundary layer theory for a micropolar fluid. Recent Adv. Eng. Sci. 1970, 5, 405–426. [Google Scholar]
  9. Ariman, T.; Turk, M.A.; Sylvester, N.D. Microcontinuum fluid mechanics—A review. Int. J. Eng. Sci. 1973, 11, 905–930. [Google Scholar] [CrossRef]
  10. Ariman, T.; Turk, M.A.; Sylvester, N.D. Applications of microcontinuum fluid mechanics. Int. J. Eng. Sci. 1974, 12, 273–293. [Google Scholar] [CrossRef]
  11. Nath, G. Similar solutions for the incompressible laminar boundary layer with pressure gradient in Micropolar fluids. Rheol. Acta 1975, 14, 190–198. [Google Scholar] [CrossRef]
  12. Ahmadi, G. Self-similar solution of incompressible micropolar boundary layer flow over a semi-infinite plate. Int. J. Eng. Sci. 1976, 14, 639–646. [Google Scholar] [CrossRef]
  13. Sadia, S.; Naheed, B.; Anwar, H.; Muhammad, A.; Rama, G.; Qasem, A.-M. Effect of thermal ra-diation on conjugate natural convection flow of a micropolar fluid along a vertical surface. Comput. Math. Appl. 2020, 83, 74–83. [Google Scholar] [CrossRef]
  14. Sadia, S.; Faryad, A.; Begum, N.; Hossain, M.A.; Gorla, R.S.R. Periodic magnetohydrodynamic natural convection flow of a micropolar fluid with radiation. Int. J. Therm. Sci. 2017, 111, 215–222. [Google Scholar] [CrossRef]
  15. Bayliss, L.E. Rheology of Blood and Lymph, Deformation and Flow in Biological Systems; Fry-Wyssling, A., Ed.; North-Holland Publishing Company: Amsterdam, The Netherlands, 1952; p. 354. [Google Scholar]
  16. Fung, Y.C. Theoretical considerations of the elasticity of red cells and small blood vessels. Fed. Am. Soc. Exp. Biol. 1966, 25, 1761–1772. [Google Scholar]
  17. Lew, H.S.; Fung, Y.C. Entry flow into blood vessels at arbitrary Reynolds number. J. Biomech. 1970, 3, 23–28. [Google Scholar] [CrossRef]
  18. Whitmore, R.L. Rheology of the Circulation; Pergamon Press Oxford: Oxford, UK, 1968. [Google Scholar]
  19. Turk, M.A.; Ariman, T.; Sylvester, N.D. On Pulsatile Blood Flow. Trans. Soc. Rheol. 1973, 17, 1–21. [Google Scholar] [CrossRef]
  20. Ariman, T.; Turk, M.A.; Sylvester, N.D. On Steady and Pulsatile Flow of Blood. ASME J. Appl. Mech. 1974, 41, 1. [Google Scholar] [CrossRef]
  21. Hogen, H.A.; Hensiksen, M. An Elevation of a Micropolar Model for Blood Flow Through an Idealized Stenosis. J. Biomech. 1989, 22, 211–218. [Google Scholar] [CrossRef] [PubMed]
  22. Bugliarello, G.; Sevilla, J. Velocity distribution and other characteristics of steady and pulsatile blood flow in fine glass tubes. Biorheology 1970, 7, 85–107. [Google Scholar] [CrossRef] [PubMed]
  23. Goldsmith, H.L.; Skalak, R. Hemodynamics. Annu. Rev. Fluid Mech. 1975, 7, 231–247. [Google Scholar] [CrossRef]
  24. Cocklet, G.R. The Rheology of Human Blood. In Biomechanics: Its Foundation and Objectives; Prentice-Hall: Hoboken, NJ, USA, 1972; p. 63. [Google Scholar]
  25. Radhakrishnamacharya, G. Flow of micropolar fluid through a constricted channel. Int. J. Eng. Sci. 1977, 15, 719. [Google Scholar]
  26. Darcy, H.; de la ville de Dijan, F. Librairie des Corps Imperianx des Pontset Chausses et des Mines. 1856. Available online: https://gallica.bnf.fr/ark:/12148/bpt6k624312/f1n657.pdf (accessed on 1 February 2023).
  27. Ouni, M.; Selimefendigil, F.; Hatem, B.; Kolsi, L.; Omri, M. Utilization of wavy porous layer, magnetic field and hybrid nanofluid with slot jet impingement on the cooling performance of conductive panel. Int. J. Numer. Methods Heat Fluid Flow 2022, 33, 360–384. [Google Scholar] [CrossRef]
  28. Srinivasacharya, D.; Ramana Murthy, J.V.; Venugopalan, D. Unsteady stokes flow of micropolar fluid between two parallel porous plates. Int. J. Eng. Sci. 2001, 39, 1557–1563. [Google Scholar] [CrossRef]
  29. Bhargara, R.; Kumar, L.; Takhar, H.S. Numerical solution of free convection MHD micropolar fluid flow between two parallel porous vertical plates. Int. J. Eng. Sci. 2003, 41, 123–136. [Google Scholar] [CrossRef]
  30. Zueco, J.; Eguia, P.; Lopez-Ochoa, L.M.; Collazo, J.; Patino, D. Unsteady MHD free convection of a micropolar fluid between two parallel porous vertical walls with convection from the ambient. Int. Commun. Heat Mass Transf. 2009, 36, 3–9. [Google Scholar] [CrossRef]
  31. Khaled, A.R.A.; Vafai, K. The role of porous media in modeling flow and heat transfer in biological tissues. Int. J. Heat Mass Transf. 2003, 46, 4989–5003. [Google Scholar] [CrossRef]
  32. Mehmod, O.U.; Mustapha, N.; Snafie, S. Unsteady Two–Dimensional Blood Flow in Porous Artery with Multiirregular Stenoses. Transp. Porous Med. 2012, 92, 259–275. [Google Scholar] [CrossRef]
  33. Abdullah, E.Y. Study of pressure applied to blood vessels using a mathematical model. Int. J. Nonlinear Anal. Appl. 2022, 13, 1341–1350. [Google Scholar] [CrossRef]
  34. Bhatti, M.M.; Abbas, M.A. Simultaneous effects of slip and MHD on peristaltic blood flow of Jeffrey fluid model through a porous medium. Alex. Eng. J. 2016, 55, 1017–1023. [Google Scholar] [CrossRef] [Green Version]
  35. Misra, J.C.; Adhikary, S.D. MHD oscillatory channel flow, heat and mass transfer in a physiological fluid in presence of chemical reaction. Alex. Eng. J. 2016, 55, 287–297. [Google Scholar] [CrossRef] [Green Version]
  36. Wang, L.; Chu, X.; Wan, J.; Xiu, C. Implementation of micropolar fluids model and hydrodynamic behavior analysis using user-defined function in FLUENT. Adv. Mech. Eng. 2020, 12, 1687814020943052. [Google Scholar] [CrossRef]
  37. Joseph, S.P. Some exact solutions for incompressible couple stress fluid flows. Malaya J. Mat. 2020, S, 648–652. [Google Scholar] [CrossRef]
  38. Baranovskii, E.S.; Burmasheva, N.V.; Prosviryakov, E.Y. Exact Solutions to the Navier–Stokes Equations with Couple Stresses. Symmetry 2021, 13, 1355. [Google Scholar] [CrossRef]
  39. Kocić, M.M.; Stamenković, Ž.M.; Petrović, J.D.; Bogdanović-Jovanović, J.B.; Nikodijević, M.D. Heat transfer in micropolar fluid flow under the influence of magnetic field. Therm. Sci. 2016, 20, S1391–S1404. [Google Scholar] [CrossRef]
  40. Ahmad, S.; Ali, K.; Bashir, H. Interaction of micropolar fluid structure with the porous media in the flow due to a rotating cone. Alex. Eng. J. 2021, 60, 1249–1257. [Google Scholar] [CrossRef]
  41. Zohra, F.T.; Uddin, M.J.; Ismail, A.I.M. Magnetohydrodynamic bio-nanoconvective Naiver slip flow of micropolar fluid in a stretchable horizontal channel. Heat Transf.—Asian Res. 2019, 48, 3636–3656. [Google Scholar] [CrossRef]
  42. Mekheimer, K.S.; El Kot, M.A. The micropolar fluid model for blood flow through a tapered artery with a stenosis. Acta Mech. Sin. 2008, 24, 637–644. [Google Scholar] [CrossRef]
  43. Nisar, K.S.; Faridi, A.A.; Ahmad, S.; Khan, N.; Ali, K.; Jamshed, W.; Abdel-Aty, A.-H.; Yahia, I.S. Cumulative Impact of Micropolar Fluid and Porosity on MHD Channel Flow: A Numerical Study. Coatings 2022, 12, 93. [Google Scholar] [CrossRef]
  44. Petrović, J.D.; Stamenković, Ž.M.; Kocić, M.M.; Nikodijević, M.D. Porous medium magnetohydrodynamic flow and heat transfer of two immiscible fluids. Therm. Sci. 2016, 20, S1405–S1417. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Physical model.
Figure 1. Physical model.
Fluids 08 00093 g001
Figure 2. Velocity profiles for different values of spin viscosity parameter Γ .
Figure 2. Velocity profiles for different values of spin viscosity parameter Γ .
Fluids 08 00093 g002
Figure 3. Microrotation for different values of spin viscosity parameter Γ .
Figure 3. Microrotation for different values of spin viscosity parameter Γ .
Fluids 08 00093 g003
Figure 4. Influence of spin viscosity parameter Γ on dimensionless temperature θ .
Figure 4. Influence of spin viscosity parameter Γ on dimensionless temperature θ .
Fluids 08 00093 g004
Figure 5. Influence of coupling parameter K on velocity u * .
Figure 5. Influence of coupling parameter K on velocity u * .
Fluids 08 00093 g005
Figure 6. Influence of coupling parameter K on microrotation ω * .
Figure 6. Influence of coupling parameter K on microrotation ω * .
Fluids 08 00093 g006
Figure 7. Dimensionless temperature θ in function of coupling parameter K .
Figure 7. Dimensionless temperature θ in function of coupling parameter K .
Fluids 08 00093 g007
Figure 8. Influence of Hartmann number on velocity u * .
Figure 8. Influence of Hartmann number on velocity u * .
Fluids 08 00093 g008
Figure 9. Influence of Hartmann number on microrotation ω * .
Figure 9. Influence of Hartmann number on microrotation ω * .
Fluids 08 00093 g009
Figure 10. Dimensionless temperature θ as a function of Hartmann number.
Figure 10. Dimensionless temperature θ as a function of Hartmann number.
Fluids 08 00093 g010
Figure 11. Nusselt number at upper and lower plates as a function of Hartmann number and coupling parameter K .
Figure 11. Nusselt number at upper and lower plates as a function of Hartmann number and coupling parameter K .
Fluids 08 00093 g011
Figure 12. Shear stress at upper and lower plates as a function of Hartmann number and coupling parameter K .
Figure 12. Shear stress at upper and lower plates as a function of Hartmann number and coupling parameter K .
Fluids 08 00093 g012
Figure 13. Velocity profiles u * for different values of porosity parameter R .
Figure 13. Velocity profiles u * for different values of porosity parameter R .
Fluids 08 00093 g013
Figure 14. Microrotation vector ω * as a function of porosity parameter R .
Figure 14. Microrotation vector ω * as a function of porosity parameter R .
Fluids 08 00093 g014
Figure 15. Shear stress on upper and lower plate for different values of porosity parameter R and coupling parameter K .
Figure 15. Shear stress on upper and lower plate for different values of porosity parameter R and coupling parameter K .
Fluids 08 00093 g015
Figure 16. Dimensionless temperature θ for different values of porosity parameter R .
Figure 16. Dimensionless temperature θ for different values of porosity parameter R .
Fluids 08 00093 g016
Figure 17. Nusselt number on upper and lower plate for different values of porosity parameter R and coupling parameter K .
Figure 17. Nusselt number on upper and lower plate for different values of porosity parameter R and coupling parameter K .
Fluids 08 00093 g017
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kocić, M.; Stamenković, Ž.; Petrović, J.; Bogdanović-Jovanović, J. Control of MHD Flow and Heat Transfer of a Micropolar Fluid through Porous Media in a Horizontal Channel. Fluids 2023, 8, 93. https://doi.org/10.3390/fluids8030093

AMA Style

Kocić M, Stamenković Ž, Petrović J, Bogdanović-Jovanović J. Control of MHD Flow and Heat Transfer of a Micropolar Fluid through Porous Media in a Horizontal Channel. Fluids. 2023; 8(3):93. https://doi.org/10.3390/fluids8030093

Chicago/Turabian Style

Kocić, Miloš, Živojin Stamenković, Jelena Petrović, and Jasmina Bogdanović-Jovanović. 2023. "Control of MHD Flow and Heat Transfer of a Micropolar Fluid through Porous Media in a Horizontal Channel" Fluids 8, no. 3: 93. https://doi.org/10.3390/fluids8030093

APA Style

Kocić, M., Stamenković, Ž., Petrović, J., & Bogdanović-Jovanović, J. (2023). Control of MHD Flow and Heat Transfer of a Micropolar Fluid through Porous Media in a Horizontal Channel. Fluids, 8(3), 93. https://doi.org/10.3390/fluids8030093

Article Metrics

Back to TopTop