Thesis Main Library
Thesis Main Library
TRANSFER PROBLEMS
INVOLVING STRATIFICATION
AND CHEMICAL REACTION
Digbash Sahu
Enrolment No. : Math/06/2021
Registration No. : 17083756 of 2021-2022
Date of Admission : 02/02/2022
Department of Mathematics
Gauhati University
I hereby declare that the thesis titled “A Study on Certain Heat Transfer
Problems Involving Stratification and Chemical Reaction” is the result of my
own research work which has been carried out under the guidance of Prof. Rudra
Kanta Deka, Professor, Department of Mathematics, Gauhati University. I further
declare that this thesis as a whole or any part thereof has not been submitted to any
university (or institute) for the award of any degree or diploma.
This thesis contains less than 90,000 (ninety thousand) words excluding bibliography
and captions.
Digbash Sahu
August 2025
Certificate
This is to certify that the thesis titled “A Study on Certain Heat Transfer
Problems Involving Stratification and Chemical Reaction” is the result of
research work of Digbash Sahu, carried under my supervision, submitted to Gauhati
University for the award of the degree of Doctor of Philosophy in Mathematics.
This thesis conforms to the standard of PhD Thesis under Gauhati University
including the standard related to plagiarism and has a similarity index not more than
10% (ten percent), excluding the bibliography.
I extend my heartfelt thanks to Dr. Rupam Shankar Nath and my fellow research
scholars—Pappu Das, Rakesh Rabha, and Nitul Kalita—for their collaborative ideas
and intellectual contributions, which greatly enriched my work.
Digbash Sahu
Gauhati University, Jalukbari
August 2025
Abstract
In Chapter 8, future research directions are outlined to advance heat and mass
transfer in stratified fluids under MHD and porous media effects. Key areas include
generalized analytical solutions, non-linear stratification, variable magnetic fields,
heterogeneous porous media, advanced nanofluids, higher-order chemical reaction,
experimental validation, environmental applications, CFD simulations, multi-physics
coupling, optimization strategies, and data-driven modeling. These efforts aim to
deepen understanding and optimize engineering and environmental systems.
List of Publications
1. Digbash Sahu and Rudra Kanta Deka, “Thermal and Mass Stratification Effects
on MHD Flow Past an Accelerated Vertical Plate with Variable Temperature and
Exponential Mass Diffusion Embedded in a Porous Medium”, East European
Journal of Physics, 2 (2024), 161-171, E-ISSN: 2312-4539, ISSN: 2312-4334,
[Link] (SCOPUS, ESCI)
3. Digbash Sahu, Rudra Kanta Deka, Pappu Das and Rakesh Rabha, “The Com-
prehensive Analysis of Stratification and Chemical Reaction on MHD Flow along
an Impulsively Started Infinite Vertical Plate in Presence of Radiation and Heat
Source Through a Porous Medium”, Heat Transfer, (2025), E-ISSN: 2688-4542,
ISSN: 2688-4534, [Link] (SCOPUS, ESCI).
4. Digbash Sahu and Rudra Kanta Deka, “Influences of Thermal Stratification and
Chemical Reaction on MHD free Convective Flow along an Accelerated Vertical
Plate with Variable Temperature and Exponential Mass Diffusion in a Porous
Medium”, Heat Transfer, 53(7) (2024), 3643-3666, E-ISSN: 2688-4542, ISSN:
2688-4534, [Link] (SCOPUS, ESCI).
5. Digbash Sahu and Rudra Kanta Deka, “Numerical study of MHD Ternary
Hybrid Nanofluid (Ag − CoF e2 O4 − ZnO/C2 H6 O2 + H2 O): Effects of Thermal
Stratification on a Vertically Stretching Cylinder in a Porous Environment”,
Biointerface Research in Applied Chemistry, 14, 6 (2024), ISSN: 2069-5837,
[Link] (SCOPUS)
xvi
6. Digbash Sahu and Rudra Kanta Deka, “MHD Natural Convection of Ternary
Nanofluid in a Porous Medium Incorporating Stratification and Radiation”,
ZAMM - Journal of Applied Mathematics and Mechanics, 105(5), e70036, (2025),
E-ISSN: 1521-4001, ISSN: 0044-2267, [Link] (SCIE,
SCOPUS).
List of Conferences Presented/Attended
1 Introduction 1
1.1 Basic Overview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Classification of Fluids . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2.1 Compressible and Incompressible fluids . . . . . . . . . . . . . . 2
1.2.2 Ideal and Real fluids . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2.3 Newtonian and Non-Newtonian fluids . . . . . . . . . . . . . . . 3
1.3 Types of Flows of Fluid . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3.1 Laminar and Turbulent flows . . . . . . . . . . . . . . . . . . . 3
1.3.2 Steady and Unsteady Flows . . . . . . . . . . . . . . . . . . . . 4
1.3.3 Uniform and Non-Uniform Flows . . . . . . . . . . . . . . . . . 4
1.3.4 Rotational and Irrotational Flows . . . . . . . . . . . . . . . . . 5
1.4 Basic Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4.1 Fluid pressure . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4.2 Fluid density . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.4.3 Fluid Temperature . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.4.4 Thermal Conductivity . . . . . . . . . . . . . . . . . . . . . . . 6
1.4.5 Specific Heat . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.4.6 Kinematic Viscosity . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.5 Stratification Effect . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.5.1 Thermal Stratification . . . . . . . . . . . . . . . . . . . . . . . 8
1.5.2 Mass Stratification . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.6 Magnetohydrodynamics(MHD) . . . . . . . . . . . . . . . . . . . . . . 9
1.7 Porous Medium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
xx Table of contents
References 159
List of figures
6.1 Thermal and Physical Properties of Ethylene glycol+ water and nanopar-
ticles (Paul et al. [1, 2]) . . . . . . . . . . . . . . . . . . . . . . . . . . 124
6.2 Comparable values of −θ′ (0) with δ = γ = M = K = Q = λ = ϕ1 =
ϕ2 = ϕ3 = 0 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125
6.3 Skin friction Coefficient and Local Nusselt number for different values
of δ, M, K, λ, P r, Q, γ . . . . . . . . . . . . . . . . . . . . . . . . . . . . 129
6.4 Comparison of −Cf Re1/2 x . . . . . . . . . . . . . . . . . . . . . . . . . 129
6.5 Comparison of Local Nusselt Number . . . . . . . . . . . . . . . . . . 129
7.1 Thermal and Physical Properties of Ethylene glycol+ water and nanopar-
ticles (Paul et al. [1, 2], and Ghadikolaei et al. [3, 4]) . . . . . . . . . . 141
7.2 Comparable values of −θ′ (0) with δ = γ = M = K = R = Q = λ =
ϕ1 = ϕ2 = ϕ3 = 0 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 142
7.3 Cf Re1/2
x and N ux Re−1/2
x for different values of δ, M, K, λ, P r, R, Q, γ . . 150
7.4 Comparison of −Cf Rex 1/2
. . . . . . . . . . . . . . . . . . . . . . . . . 151
7.5 Comparison of N ux Rex −1/2
. . . . . . . . . . . . . . . . . . . . . . . . . 151
Introduction
“Research is seeing what everybody else has seen and thinking what
1
nobody else has thought” -Albert Szent-Györgyi
and mixing within the fluid. Understanding these two types of flow is crucial in
fluid dynamics, as they play a significant role in applications ranging from industrial
processes to environmental studies.
Understanding these flow types is critical for designing and analyzing systems in fluid
dynamics and engineering applications.
where, δS is an elementary area surrounding the given point and δF is the normal force
exerted by the fluid on δS. The S.I. unit of pressure is Pascal (Pa) or Newton/m2 .
where, δv represents the volume element around a point and δm is the fluid mass
within δv. The S.I. unit of density is Kg/m3 .
behavior. Density stratification significantly influences fluid motion and the deposition
of sediments. For instance, the flow and settling behavior of particles, droplets, bubbles,
and even microscopic organisms are governed by the density gradients present in
the fluid. This phenomenon is critical in numerous fields. Meteorologists rely on an
understanding of stratified flows to predict weather patterns and analyze atmospheric
phenomena, such as the formation of clouds, storms, and temperature inversions.
Similarly, engineers consider density stratification when designing and constructing
structures like dams and coastal defenses to ensure stability and performance under
varying conditions. Stratification has far-reaching implications for environmental,
ecological, and geophysical phenomena. A prime example of natural stratification is
illustrated in Fig. 1.1. In natural water bodies, thermal stratification occurs during the
summer, with warmer, less dense water forming a distinct layer above cooler, denser
water. Conversely, during winter, inverse stratification arises when surface cooling
results in a denser upper layer. These seasonal patterns profoundly affect aquatic
ecosystems and sediment transport dynamics. Consequently, the study of stratification
in fluid dynamics is essential for understanding its pivotal role in influencing natural
processes and the design of human-engineered environments.
γ ′ν
γ=
u0 (Tw′ − T∞
′ )
1.6 Magnetohydrodynamics(MHD)
Magnetohydrodynamics (MHD) is the study of electrically conducting fluid motion
in the presence of a magnetic field, encompassing examples like plasma, electrolytes,
liquid metals, and salt water. This concept highlights the ability of magnetic fields
to induce electric currents within moving fluids, which in turn generate forces that
oppose the fluid’s motion and modify the magnetic field. Pioneered by Hannes Alfvén,
who received the Nobel Prize in 1970 for his contributions, MHD has gained significant
attention for its wide range of applications across scientific disciplines. The field is
grounded in two fundamental principles:
(i) electric currents produce an induced magnetic field that alters the original
magnetic field. (ii) an electromagnetic force arises from the interaction between the
current and the magnetic field perturbs the fluid’s motion.
MHD fluid flow studies are invaluable in various domains, including engineering,
mechanics, and applied sciences such as astrophysics and geophysics.
and heat exchange systems. The study of porous media is essential for understanding
and optimizing their role in various natural and engineered processes.
In Fig 1.2 porous media play a critical role in dam functionality and safety by
controlling seepage through the foundation and surrounding areas, using materials
like sand or gravel to filter and channel water safely. They are integral to drainage
systems and embankment dams, where porous layers reduce pressure and prevent
erosion. Additionally, they influence water infiltration and groundwater recharge in the
reservoir and floodplain, helping manage water resources effectively. Engineers analyze
the permeability of porous materials to ensure structural stability and prevent issues
like waterlogging or failure.
1.11 Nanofluid
A nanofluid is a suspension of nanoparticles that demonstrates significant potential
as a heat transfer fluid, primarily due to its exceptional thermal conductivity and
superior rheological properties, making it highly effective across various applications.
The term "nanofluid" was first introduced by Stephen U. S. Choi in his pioneering work
presented at the Winter Annual Conference of the American Society of Mechanical
Engineers (ASME) in 1995 [6]. Nanofluids have been explored for numerous applications,
including heat exchangers, transportation cooling, refrigeration, electronic equipment
cooling, transformer oil, industrial cooling, nuclear systems, machining operations,
solar energy utilization, desalination, and defense technologies. Their enhanced heat
transfer characteristics and stability make them a promising solution for improving
thermal performance in both industrial and scientific domains.
1.12.1 Conduction
Conduction is the process by which heat is transferred from the hot end to the cool
end of a body through the movement of particles, without any physical displacement
of the particles. It takes place in solid.
14 Introduction
1.12.2 Convection
It is that mode of transmission of heat in which heat flows from the hot region to the
cold region of the medium from particle to particle due to the actual displacement of
the particles of the medium. The phenomenon occurs within a fluid medium. It can
be categorized as follows:
Forced Convection
Forced convection is an interesting process of heat transmission in which there is
an external force, such as a fan, pump, or suction device, that causes fluid motion.
In everyday life and numerous engineering scenarios, forced convection is frequently
encountered. Examples of such scenarios are central heating, air conditioning, steam
turbines, heat exchangers, pipe flow, and flow over a plate.
Natural Convection
Natural convection, sometimes called free convection, is an important process of heat
and mass transport that happens in fluids (gases, liquids, etc.) as a result of density
variations caused by temperature gradients. Natural convection occurs naturally, unlike
forced convection which depends on external sources such as pumps or fans. This
type of convection is commonly observed in various everyday situations, including the
upward movement of hot air near a radiator or heater, the circulation of air within a
room resulting from temperature differences, and the upward flow of heated water in a
pot.
Mixed Convection
In fluid dynamics, mixed convection, also called combined convection, is an interesting
process where heat movement is affected by both forced convection and natural
convection at the same time. It can be observed in many different engineering equipment,
particularly when heating or cooling channel walls at low fluid velocities.
1.12.3 Radiation
Radiation is a mode of heat transfer in which thermal energy moves from a hot body to
a cold body without requiring an intervening medium. A classic example is the transfer
of heat from the Sun to the Earth, where energy travels through the vacuum of space
as electromagnetic waves. Thermal radiation occurs when a body’s thermal energy
is converted into electromagnetic energy, exhibiting wave-particle duality, meaning it
possesses both wave-like and particle-like characteristics. The significance of radiative
heat transfer increases at higher temperatures, often surpassing the effects of conduction
and convection in extreme thermal environments. This phenomenon is essential in
numerous combustion applications, including fires, furnaces, engines, and rocket nozzles.
Additionally, thermal radiation plays a critical role in nuclear reactions, such as those
1.13 Mass Transfer in Fluids 15
in the Sun, fusion reactors, and nuclear weapons, where immense energy release occurs
through radiative processes.
1.13.2 Convection
Mass transfer through convection occurs as a result of the overall motion of fluid.
The subject matter pertains to the exchange of mass between two fluids that are
substantially immersible in motion, or between a fluid in motion and a surface. The
process of alcohol evaporation involves the transfer of mass through free convection.
The process by which water vapour evaporates from an ocean or big reservoir as a result
of airflow over it can be described as a mass transfer facilitated by forced convection.
The compelled movement of the fluid mass results in the transfer of mass through
forced convection. Convective mass transfer and convective heat transfer work in very
similar ways.
gβ ∗ L3 (Cw − C∞ )
Gm =
ν2
where, g, β ∗ , ν, L, Cw , C∞ is the acceleration due to gravity, volumetric coefficient
of species concentration, kinematic viscosity, characteristic length, concentration at
wall and far away from the wall respectively.
1.15 Non-Dimensional Coefficient 17
K ∗ u20
Da =
ν2
where, K ∗ is the permeability of the medium.
νσB02
M=
ρ∞ u20
The Hartmann number provides a quantitative assessment of the relative significance
of drag forces arising from magnetic induction and viscous forces in Hartmann flow,
hence establishing the velocity profile associated with this particular flow.
1.16 Methodology
1.16.1 Laplace Transform Method
The Laplace transform was named after the French mathematician Pierre Simon De
Laplace (1749-1827). The process involves transforming a function representing a
real variable, typically in the time domain, into a function representing a complex
variable, typically in the complex-valued frequency domain, also referred to as the
s-domain or s-plane. The application of the Laplace transform facilitates the resolution
of differential equations by converting them into algebraic mathematical problems. It
has an important effect on the field of control system engineering and the study of
dynamic control systems.
Given, a function F (t) of t for t > 0. Then, the Laplace transformation of F (t),
denoted by L{F (t)} or f (s) is defined as follows
ˆ ∞
L{F (t)} = f (s) = e−st F (t) dt
0
Here, L and L−1 are Laplace transform operator and inverse Laplace transform operator,
respectively.
bvp4c represents “boundary value problem solver using the fourth-order method.” This
method is used to solve systems of differential equations of the form y ′ = f (x, y).
The solver incorporates the boundary conditions as specified by a function (bcfun)
and generates an initial solution estimate (solinit). The integration settings can be
customised by utilising an options argument that is generated through the bvpset
function. The bvp4c algorithm utilises a finite difference technique that depends on
the three-stage Lobatto IIIa formula. The collocation formula presented above yields a
C 1 continuous solution that exhibits uniform fourth-order accuracy within the interval
[a, b].
The study of transient and steady-state free convection flows over vertical surfaces
has remained a key area of research due to its significance in fluid dynamics and heat
transfer applications. Early studies laid the foundation for understanding thermal and
buoyancy-driven flows under different conditions. Siegel [7] provided an analytical
solution for transient free convection from a vertical flat plate, establishing fundamental
principles governing thermal boundary layers. Later, Menold and Yang [8] developed
asymptotic solutions for unsteady laminar free convection over vertical plates, refining
insights into transient boundary layer behavior, while Schetz and Eichhorn [9] examined
unsteady natural convection near an infinitely long vertical plate, offering valuable
findings on large-scale heat transfer. Hossain and Takhar [10] investigated thermal
radiation effects on mixed convection along a vertical plate with uniform surface
temperature, highlighting its role in high-temperature applications. Das et al. [11]
analyzed transient free convection flow past an infinite vertical plate with periodic
temperature variations, emphasizing the influence of thermal oscillations, while Hossain
et al. [12] extended this to porous vertical plates, further elucidating thermal energy
exchange in porous media. Park and Hyun [13] examined the transient behavior
of buoyancy layers in stratified fluids, which Park [14] later expanded to buoyant
flows in vertical channels. Das et al. [15] formulated an analytical solution for
chemically affected fluid flow past an infinite vertical plate, providing a basis for
chemically reactive flow studies. Subsequent studies by Shapiro and Fedorovich [5],
20 Introduction
Magyari et al. [16], and Bhattacharya and Deka [17] explored unsteady flows in
stably stratified fluids over infinite plates, enhancing the understanding of stratification
effects. Rajesh and coworkers [18, 19] and Ramakrishna et al. [20] investigated fluid
motion over exponentially accelerated vertical plates, analyzing different influencing
factors. Research on Magnetohydrodynamics (MHD) parabolic flow was conducted by
Muthucumaraswamy and Sivakumar [21] and Selvaraj et al. [22], focusing on infinite
and accelerating vertical plates under varying conditions. Vemula et al. [23] performed
an analytical study on MHD nanofluids flowing past an exponentially accelerating
vertical plate with variable temperature, while Jeevandhar et al. [24] extended this
analysis to MHD boundary layer nanofluid flow, incorporating the effects of thermal
radiation.
The study of heat transfer through vertical cylinders was first explored by Sparrow
and Gregg [25], laying the foundation for further investigations into convective heat
transfer in cylindrical geometries. Later, Yang [26] examined the potential similarity
solution for laminar free convection on vertical cylinders, providing deeper insights
into the governing equations. Goldstein and Briggs [27] extended this research by
analyzing transient free convective flow past a vertical flat plate and circular cylinder
using the Laplace transform technique, offering a mathematical framework for unsteady
flow problems. Subsequent numerical studies on unsteady natural convection past
a semi-infinite vertical cylinder, incorporating heat and mass transfer under various
physical conditions, were conducted by Ganesan and Loganathan [28–30] and Takhar
et al. [31]. Further, Rani [32] carried out a numerical investigation on transient
natural convection along a vertical cylinder, considering the effects of variable surface
temperature and mass diffusion. The impact of thermal and mass stratification on
infinite vertical cylinders was later examined by Deka and Paul [33] and Paul and Deka
[34], providing crucial insights into stratification effects in convective heat transfer.
The study of stratification effects in fluid mechanics and heat transfer is crucial due
to the prevalence of stratified fluids in both environmental and engineering applications.
Early experimental studies by Hyun [35] and Lin and Armfield [36] demonstrated that
self-stratification can arise from initially unstratified convective motion, highlighting
the complex nature of stratified flows. Further investigations by Kandasamy et al. [37]
and Mansour et al. [38] analyzed the combined effects of chemical reaction and thermal
stratification on MHD flow over vertically stretching surfaces. In another study, Cheng
[39, 40] examined thermal and mass stratification over a vertical wavy truncated cone
and a wavy surface, respectively, providing insights into the influence of geometric
variations on convective heat and mass transfer. Extending this research, Deka and
Paul [33] and Paul and Deka [34] investigated the thermal and mass stratification
effects on infinite vertical cylinders, further contributing to the understanding of
stratified convection in cylindrical geometries. Srinivasacharya and Surender [41–43]
explored mixed convection heat and mass transfer along a vertical plate in a doubly
stratified, non-Darcy porous medium, considering various flow conditions. Additionally,
1.17 Review of the Relevant Literature 21
Srinivasacharya et al. [44, 45] derived numerical solutions for free convection boundary
layer flow past a semi-infinite vertical plate in a stratified porous medium. Their work
was extended to nanofluid-saturated porous media, where Srinivasacharya and Surender
[46, 47] analyzed thermal and mass stratification effects on mixed and natural convection
boundary layer flows over a vertical flat plate. Moreover, Srinivasacharya and Surender
[48] obtained non-similarity solutions for natural convection heat and mass transport
along a vertical plate with uniform wall temperature and concentration, further refining
the theoretical framework for stratified porous media flows. Recent advancements in
hybrid nanofluid flows have introduced new dimensions to stratified convection research.
Paul et al. [49] investigated the thermal stratification effects of a hybrid nanofluid
(Cu − Al2 O3 /H2 O) in a porous medium, focusing on a vertically stretched cylinder
with heat source/sink effects. Sarma and Paul [50] extended this analysis by examining
the thermophoresis and Brownian motion effects on bioconvective cylindrical-shaped
Ag − CuO/H2 O Ellis hybrid nanofluid flow along a radiative stretched tube under
an inclined magnetic field, shedding light on bioconvective nanofluid interactions.
Expanding on hybrid nanofluid flows, Paul et al. [1] analyzed the Darcy–Forchheimer
flow of Ag − ZnO − CoF e2 O4 /H2 O Casson ternary hybrid nanofluid, induced by a
rotating disk under electro-magnetohydrodynamics (EMHD), revealing the impact of
porosity and external forces in advanced nanofluid applications. Similarly, Paul et al.
[51] explored Casson ternary hybrid nanofluid flow in transformer oil, considering the
effects of a porous rotating disk and Hall current, addressing thermal and electrical
performance improvements in energy systems. Additionally, Paul et al. [52] examined
the mixed convective flow of an engine oil-based non-Newtonian tri-hybrid nanofluid
across a porous rotating disk, expanding knowledge in lubrication, industrial cooling,
and thermal management applications.
radiation and chemical reaction on the unsteady MHD flow of a viscous, incompressible
fluid, considering a parabolic starting velocity over an infinite, isothermal, vertical
porous plate.
Choi and Eastman [6] were pioneers in introducing the concept of nanofluids,
proposing that suspending metallic nanoparticles in conventional heat transfer fluids
could lead to a revolutionary class of enhanced heat transfer fluids. Subsequent studies
explored the impact of various nanoparticles on heat transfer and fluid flow character-
istics. Oztop [59] conducted a numerical analysis of different types of nanoparticles
in a partially heated rectangular enclosure, investigating their influence on thermal
behavior. Das and Jana [60, 61] examined natural convective flow of nanofluids with
radiation effects, considering moving vertical plates and vertical channels, respectively.
Abolbashari et al. [62] employed the Homotopy Analysis Method (HAM) to evaluate
the entropy generation of nanofluids composed of water as the base fluid and four types
of nanoparticles (T iO2 , Al2 O3 , Cu, and CuO) flowing over a stretchable permeable
surface. Further investigations into heat and mass transfer in nanofluids were conducted
by Reddy et al. [63] and Mahanthesh et al. [64], who explored the effects of heat
sources and chemical reaction on nanofluids passing a moving or fixed vertical plate.
Rashid et al. [65] analyzed the magnetohydrodynamic (MHD) behavior of Al2 O3 -Cu
hybrid nanofluids in porous media, emphasizing the impact of nanoparticle shape on
flow dynamics. Wang et al. [66] studied the thermal behavior and natural convection
of MWCNT-water nanofluid flow inside a hexagonal cavity under a magnetic field,
furthering the understanding of MHD-driven nanofluid convection. Similarly, Kumar
Rawat et al. [67] investigated the flow of a T iO2 /EG nanofluid over a wedge, offering
insights into nanoparticle-fluid interactions in boundary layer flows. Yaseen et al. [68]
examined the combined effects of buoyancy forces and Hall currents on the stagnation
point flow of Ag − T iO2 /water hybrid nanofluids over a rotating disc, highlighting
the time-dependent flow characteristics. Additionally, Yaseen et al. [69] obtained
numerical solutions for a wall jet flow problem involving MHD in an Eyring-Powell
hybrid nanofluid composed of water, CoF e2 O4 , and T iO2 nanoparticles. Finally, Riaz
et al. [70] discovered new analytical solutions for the fractional Maxwell fluid contain-
ing ternary-hybrid nanoparticles flowing over an infinitely tall vertical heated wall,
contributing to the growing body of research on advanced hybrid nanofluid applications.
tion. This contrasts with many earlier studies that relied on approximate or purely
numerical schemes. For instance, the findings align with the velocity trends reported
by Muralidharn et al. [71] in the absence of stratification, but the introduction of
stratification parameters in this thesis reveals a stabilizing effect that leads to lower
peak velocities and a more stable flow regime.
Furthermore, unlike studies that focused solely on thermal effects, this research
integrates mass stratification and chemical reactions, offering a more comprehensive
model. The results indicate that the interplay between these factors significantly alters
heat and mass transfer coefficients, leading to higher Nusselt and Sherwood numbers
than predicted in studies considering these effects in isolation. This highlights the
novelty of the current approach in providing a more accurate and holistic understanding
of stratified flows in complex environments.
Thermal and Mass Stratification
Effects on MHD flow past an
2
Accelerated Vertical Plate with
Variable Temperature and
Exponential Mass Diffusion
Embedded in a Porous Medium
Published in East European Journal of Physics, 2, 2024, pp. 161-171,
ISSN: 2312-4334. (SCOPUS, ESCI)
2.1 Introduction
The investigation of heat transfer holds importance across various engineering fields,
such as in the cooling of electronic devices, nuclear reactors, and gas turbines. Similarly,
a grasp of mass transfer is crucial for several chemical engineering applications, like
separation processes, distillation, and absorption. This study aims to explore the
combined effects of thermal and mass stratification on the behavior of MHD unsteady
flow along an accelerated vertical plate, considering variations in mass diffusion and a
variable temperature. The influence of thermal and mass stratifications on the behavior
of MHD unsteady flows in fluids is significant, and this research seeks to analyze the
interaction of these stratifications.
26 Thermal and Mass Stratification Effects on MHD flow past an Accelerated...
Early work by Siegel [7] and Menold and Yang [8] laid the groundwork for un-
derstanding transient free convection from vertical flat plates. Schetz and Eichhorn
[9] extended this by investigating unsteady natural convection near doubly infinite
vertical plates. Later studies explored diverse aspects of this phenomenon. Das et al.
[11] and Park et al. [13, 14] investigated transient buoyant flows in stratified fluids,
while Shapiro and Fedorovich [5] explored convectively driven flows in stably stratified
fluids. Porous media were explored in greater detail. Magyari et al. [16] investigated
unsteady free convection in a fluid-saturated porous medium, and Chaudhary and
Jain [72] considered heat and mass diffusion flow by natural convection in a porous
medium. The influence of radiation and magnetic fields was also examined. The
surveyed literature investigates various aspects of magnetohydrodynamic (MHD) and
porous medium flow. Reddy et al. [73] analyzed MHD boundary layer flow along
vertical plates with ramped temperature. Selvaraj and Jothi [74] study heat sources
in MHD flow past an accelerated plate with variable temperature and mass diffusion.
Sarma and Nazibuddin [75] explore thermal diffusion in unsteady MHD convective flow.
Krishna et al. [76–79] and [80] examine effects like Hall and ion slip, elastico-viscous
fluid behavior, and nanofluid flow dynamics in different MHD contexts, emphasizing
the influence of porous mediums and thermal conditions on flow properties and heat
transfer. More recent research has delved into the effects of chemical reaction in these
flows. Kalita et al. [81] and Nath et al. [82, 83] studied unsteady flow past vertical
plates with chemical reaction in the presence of thermal stratification. Lastly, Kumar
and Deka [84] extended this to porous media, considering mass diffusion, showcasing
the breadth and depth of research in this field.
In this study, we present novel contributions by deriving the exact solution through
the Laplace transform technique, achieving this with perfect accuracy, which proves
to be an extremely effective strategy for obtaining precise solutions. Prior to this
work, there has been no exploration of the combined influences of thermal and mass
stratification on the behavior of MHD unsteady flow past a vertically accelerating
plate embedded within a porous medium, where both temperature and mass diffusion
vary. The aim was to analytically explore the dynamics of MHD unsteady flow
past such a plate, taking into account the effects of thermal and mass stratification.
Subsequently, we compare the results concerning fluid stratification in both thermal
and mass aspects with those of the specific instance lacking any stratification. This
study explores and presents the effects of various physical parameters, including
γ, ξ, Gr, Gc, the Darcy number (Da), and the Magnetic parameter (M ), on the observed
profiles through graphical representations. The insights gained from this research hold
potential applications across a wide range of industries. For example, the outcomes
can be leveraged to enhance thermal system efficiency by minimizing the temperature
discrepancy between the walls and the fluid, as well as by accelerating the rate of mass
diffusion.
2.2 Mathematical Analysis 27
Energy Equation
∂T ′
" #
ρCp ′
+ (⃗q.∇)T
⃗ ′ = α∇2 T ′ (2.2)
∂t
Concentration Equation
∂C ′
+ (⃗q.∇)C
⃗ ′ = D∇2 C ′ (2.3)
∂t′
Consider the MHD unsteady flow of a viscous, incompressible, and stratified fluid
over an accelerating vertical plate within a porous medium. The analysis adopts a
rectangular Cartesian coordinate system (x′ , y ′ , t′ ), where the y ′ axis is perpendicular to
the plate and the x′ axis extends vertically upward along the plate. The fluid velocity
at any point (x′ , y ′ , t′ ) is given by q = (u′ , 0). Initially, at t′ = 0, the temperature
and concentration at the plate are T∞ ′
and C∞ ′
, respectively. For t′ > 0, the plate
28 Thermal and Mass Stratification Effects on MHD flow past an Accelerated...
u′ = 0 T ′ = T∞
′
C ′ = C∞
′
∀y ′ , t′ ≤ 0
′ ′
u′ = u0 t′ T ′ = T∞
′
+ (Tw′ − T∞
′
)At′ C ′ = C∞
′
+ (Cw′ − C∞
′
)ea t at y ′ = 0, t′ > 0
u′ = 0 T ′ → T∞ ′
C ′ → C∞ ′
as y ′ → ∞, t′ > 0
!1/3
gβ(Tw′ − T∞
′
) gβ ∗ (Cw′ − C∞
′
) ν ν ν
Gr = , Gc = , Pr = , Sc = , a = a′ ,
u0 u0 α D u20
1 !1/3
σB02 ν 3 u20 γ ′ ν 2/3 ξ ′ ν 2/3
M= 2 , Da = k ′ , γ= 1/3
, ξ= 1/3
ρu03 ν4 u0 (Tw′ − T∞
′ ) u0 (Cw′ − C∞
′ )
2.3 Method of Solution 29
u2 1/3
A = ν0 is the constant. Then non-dimensional forms of the equations (2.4)-(2.6)
are given by
∂U ∂ 2U 1
= Grθ + GcC + − M + U (2.7)
∂t ∂y 2 Da
∂θ 1 ∂ 2θ
= − γU (2.8)
∂t P r ∂y 2
∂C 1 ∂ 2C
= − ξU (2.9)
∂t Sc ∂y 2
Where,
s
1 1
2
y
η= √ , N =A+B =M + , A−B = M+ − 4(γGr + ξGc)
2 t Da Da
A − Gr Gr − B Gc
H1 = , H2 = , H3 =
A−B A−B A−B
γ(A − Gr) γ(Gr − B) γGc γGc
F1 = , F2 = , F3 = , F4 = ,
A(A − B) B(A − B) A(A − B) B(A − B)
ξ(A − Gr) ξ(Gr − B) ξGc ξGc
F5 = , F6 = , F7 = , F8 =
A(A − B) B(A − B) A(A − B) B(A − B)
Also, fi ’s are inverse Laplace’s transforms given by
√ √
e−y s+p e−y s+p
( ) ( )
f1 (p) = L −1
, f2 (p, a) = L −1
s2 s−a
We separate the complex arguments of the error function contained in the previous
expressions into real and imaginary parts using the formulas provided by Abramowitz
and Stegun [86].
2η
( )
2
θc = t 1 + 2η erf c(η) − √ e−η
2
(2.14)
π
eat h −2η√at √ √ √ i
Cc = e erf c η − at + e2η at erf c η + at (2.15)
2
Gr
Gc Gc eat n −2η√at √
Uc = 1− f1 (N ) − f2 (N, a) + e erf c η − at
N N 2N
√ 2η
( )
√ o tGr 2
+e2η at erf c η + at + 1 + 2η 2 erf c(η) − √ e−η (2.16)
N π
2.3.2 Skin-Friction
The non-dimensional Skin-Friction, which is determined as shear stress on the surface,
is obtained by
2.3 Method of Solution 31
dU
τ=−
dy y=0
The solution for the Skin-Friction is calculated from the solution of Velocity profile U ,
represented by (2.11), as follows:
√
√ √ h √ √
s
t erf ( At)
τ = H1 t A erf ( At) + e−At + √ − H2 t B erf ( Bt)
π 2 A
s √
erf ( Bt)
√
t
q
+e−Bt + √ − H3 eat a + A erf ( (a + A)t)
π 2 B
√ e−At − e−Bt
q #
− a + B erf ( (a + B)t) + √
πt
The solution for the Skin-Friction for the classical case is given from the expression
(2.16), which is represented by
q
t −N t erf ( (N )t) 2Gr
Gr √
s s
t
q
τc = 1− t N erf ( (N )t) + e + √ + +
N π 2 Nt N π
Gc at √ √ √ 1
" q #
e a erf ( at) − a + N erf ( (a + N )t) + √ 1 − e−N t
N πt
dθ
Nu = −
dy y=0
The solution for the Nusselt number is calculated from the solution of Temperature
profile θ, represented by (2.12), as follows:
32 Thermal and Mass Stratification Effects on MHD flow past an Accelerated...
√
√ √
s s
t t erf ( At)
N u = 2(1 − F1 − F2 ) + F1 t A erf ( At) + e−At + √
π π 2 A
√
√ √
s
t erf ( Bt)
+F2 t B erf ( Bt) + e−Bt + √ + (F3 − F4 )
π 2 B
√ √ 1 √ e−At
" # " q #
eat
a erf ( at) + √ − F3 e at
a + A erf ( (a + A)t) + √
πt πt
√ −Bt
" #
q e
+F4 eat a + B erf ( (a + B)t) + √
πt
The solution for the Nusselt number for the classical case is given from the expression
(2.14), which is represented by s
t
N uc = 2
π
dC
Sh = −
dy y=0
The solution for the Sherwood number is calculated from the solution of Concentration
profile C, represented by (2.13), as follows:
√ √
s
√ √ 1
" #
t
Sh = (1 + F7 − F8 ) e at
a erf ( at) + √ + F5 t A erf ( At) + e−At
πt π
√
√ √
s s
erf ( At)
#
t t
+ √ − 2(F5 + F6 ) + F6 t B erf ( Bt) + e−Bt
2 A π π
√
erf ( Bt) √ e−At
# " q #
+ √ − F7 eat a + A erf ( (a + A)t) + √
2 B πt
√ −Bt
" #
q e
+F8 eat a + B erf ( (a + B)t) + √
πt
The solution for the Sherwood number for the classical case is given from the
expression (2.15), which is represented by
√ √ 1
Shc = eat a erf ( at) + √
πt
2.4 Results and Discussion 33
Figure 2.2 illustrates the impact of thermal and mass stratification on velocity
profiles, showing that both types of stratification lead to reduced velocities. When one
form of stratification is held constant, an increase in the other type further decreases
velocity. In the same way that Shapiro and Fedorovich [5] shows that fluid velocity
drops for thermal stratification γ > 0, we find that this is also the case for mass
stratification ξ > 0. Enhancing the parameter of thermal stratification (γ) diminishes
the convective potential across the hot plate and the adjacent fluid, reducing the
buoyancy force and, subsequently, the flow velocity. Similarly, an increase in the mass
stratification (ξ) value leads to a lower concentration gradient between the surface and
its environment, diminishing the buoyancy’s upward force and thus slowing down the
fluid flow. Therefore, the presence of both thermal and mass stratification results in
a reduced fluid velocity compared to conditions without stratification. Hence, both
thermal and mass stratification play a crucial role in establishing a stable stratified
flow. As seen in Figure 2.3, Raising the value of Gc leads to a higher velocity, while a
rise in Gr causes a decrease in velocity. We can observe from Figure 2.4 that velocity
34 Thermal and Mass Stratification Effects on MHD flow past an Accelerated...
decreases as (M ) grows. Figure 2.5 shows that as Da and the fluid’s porosity increase,
2.4 Results and Discussion 35
there is an elevation in the velocity profile, enabling the fluid particles to move more
smoothly.
Figure 2.6 demonstrates the impact of thermal and mass stratification on the
temperature profile. While mass stratification leads to an increase in temperature, an
36 Thermal and Mass Stratification Effects on MHD flow past an Accelerated...
Figures 2.14, 2.15, and 2.16 illustrate the effects of thermal and mass stratification
on fluid velocity, temperature, and concentration over time. In the classical case, fluid
velocity increases indefinitely with time, whereas the introduction of both stratifications
causes it to stabilize at a steady state. The temperature profile is influenced by the
presence of stratifications, with thermal stratification exhibiting a more significant
2.4 Results and Discussion 39
In Figure 2.17, skin friction decreases over time in the presence of both thermal and
mass stratifications compared to the case without stratification. As shown in Figure
2.18, skin friction further decreases with time, with the lowest values observed at the
highest Gr while maintaining a constant Gc. Figures 2.19, 2.20, 2.21, and 2.22 depict
2.4 Results and Discussion 41
the temporal behavior of the Nusselt and Sherwood numbers, both of which increase
over time. Higher thermal stratification moderates the rate of increase for the Nusselt
number and, to a greater extent, for the Sherwood number. In contrast, varying mass
42 Thermal and Mass Stratification Effects on MHD flow past an Accelerated...
2.5 Conclusion
The findings from this research reveal that thermal and mass stratification plays a
critical role in reducing velocity over time, allowing the system to attain equilibrium.
Enhanced mass stratification and diminished thermal stratification contribute to the
44 Thermal and Mass Stratification Effects on MHD flow past an Accelerated...
3.1 Introduction
In the intricate world of fluid dynamics, Thermal Stratification and Mass Stratification
play pivotal roles. Thermal Stratification emerges from heat transfer resistance that
creates distinct layers within a fluid, whereas Mass Stratification develops from density
differences due to variable solute concentrations. These stratified layers are dynamically
modulated by periodic temperature fluctuations, echoing the cyclic thermal conditions
observed in both natural and industrial environments. Furthermore, the introduction of
Variable Mass Diffusion in porous media adds complexity, as the medium’s structural
characteristics significantly influence fluid flow and diffusion properties. Grasping
these nuanced interactions is critical for optimizing system performance and effectively
managing fluid behavior in a variety of applications.
Initial studies by Siegel [7] addressed transient free convection from a vertical flat
plate, laying foundational concepts of fluid mechanics related to heat transfer. Building
46 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
on this, Menold and Yang [8] explored asymptotic solutions for unsteady laminar
free convection, which further elucidated the physics of heat transfer along a vertical
plate. Schetz and Eichhorn [9] expanded the understanding by examining unsteady
natural convection in the vicinity of a doubly infinite vertical plate, contributing to
a broader scope of fluid dynamics. Park and Hyun [13], and later Park [14], shifted
focus towards the transient behavior and buoyant flows of stratified fluids, which are
crucial in understanding environmental and industrial processes. This trend continued
with Shapiro and Fedorovich [5] who investigated the convective flows along a vertical
plate immersed in a stably stratified fluid, providing insights relevant to atmospheric
sciences and engineering applications. Magyari, Pop, and Keller [16] addressed the
complex scenario of a vertical plate embedded in a porous medium, enhancing the
understanding of porous media’s impact on convection processes. Chaudhary and
Jain [72], and subsequent works by Reddy et al. [88], Chamkha et al. [19], and
others, progressively incorporated the effects of magnetic fields, chemical reaction,
and thermal radiation, reflecting the evolving complexity of investigations and their
applications in fields like chemical engineering and geophysical flows. Muralidharn
and Muthucumaraswamy [71] initiated the analysis with an investigation of parabolic
flow past a vertical plate with uniform heat flux and variable mass diffusion, laying
the groundwork for understanding such complex fluid flows. Subsequent studies by
Muthucumaraswamy and Sivakumar [21, 57] extended these findings to include the
effects of MHD, thermal radiation, and chemical reaction, illustrating the intricate
interactions within fluid dynamics. Kumaran and Sandeep [89] explored the impact of
thermophoresis and Brownian motion on MHD Casson and Williamson fluids, indicating
significant effects due to these non-Newtonian fluids. Krishna et al. [90] further analyzed
MHD Carreau fluid flows, highlighting how buoyancy and exponential heat sources alter
fluid behavior. Recent works by Gowri and Selvaraj [91], and Selvaraj et al. [22] focused
on rotational effects and the dynamics of accelerated flows in porous media, providing
deeper insights into how rotational forces influence MHD parabolic flows. Endalew and
Sarkar [92], along with Reddy et al. [93], continued to explore these themes within
different configurations of MHD flows in porous environments, emphasizing the role of
aligned magnetic fields and geometric orientations of plates. Sarma and Ahmed [85]
highlighted the Dufour effect in MHD flow with ramped temperature, demonstrating its
significance in modifying the transport phenomena. This theme of environmental and
operational variability was further echoed in studies by Sharma and Chaudhary [94],
Goto et al. [95], and Li et al. [96], each exploring different aspects of fluid responses
to periodic and environmental changes. The work by Jia et al. [97] investigates the
ecological impacts, linking fluid dynamics with environmental science by studying how
temperature variations influence aquatic carbon cycles. Krishna et al. [98] explored the
influence of Hall effects on the peristaltic flow of Jeffrey fluid through a porous medium
within a vertical stratum, setting the stage for understanding MHD flow dynamics
in complex fluid systems. Building on this, Krishna and Chamkha [99] examined the
Hall effects on MHD squeezing flow of a water-based nanofluid between two parallel
disks, further extending the understanding of fluid behavior in confined geometries and
3.1 Introduction 47
their applications in heat transfer. Krishna et al. [100] later addressed the impact of
heat and mass transfer in micropolar fluids, considering the inclined magnetic field and
Hall effects, contributing valuable insights into the complexities of fluid flow through
porous surfaces under magnetic influences. This exploration continued with Krishna
et al. [101], who investigated the combined effects of Soret, Joule, and Hall forces
on MHD mixed convective flow past an infinite vertical porous plate, revealing the
intricate balance of these forces in ocean engineering and environmental applications.
Krishna and Chamkha [102] expanded the scope by considering the Hall and ion slip
effects on unsteady MHD convective rotating flows of nanofluids, a study particularly
relevant for biomedical engineering applications due to its focus on fluid behavior under
rotational and electromagnetic forces. In another significant contribution, Krishna
et al. [77] investigated the heat and mass transfer characteristics of MHD second-
grade fluids through a porous medium over a semi-infinite vertical stretching sheet,
contributing to the growing body of knowledge on non-Newtonian fluid dynamics in
porous media. Further advancing the field, Krishna [103] examined Hall and ion slip
impacts on unsteady MHD free convective rotating flows of Jeffrey fluid with ramped
wall temperature, providing a detailed analysis of how ramped conditions influence
flow stability and convective heat transfer. This research was followed by Krishna
[104], who explored the radiative MHD rotating flow of Jeffrey fluid past an infinite
vertical porous surface, highlighting the interplay between radiation, ion slip, and
Hall effects in fluid mechanics. Krishna [78] further investigated unsteady MHD free
convective flows through a saturated porous medium, focusing on how exponential
acceleration of the plate influences flow dynamics and thermal characteristics. Krishna
and Chamkha [80] addressed the Hall and ion slip effects on the rotating boundary
layer flow of nanofluids past an infinite vertical plate embedded in a porous medium,
enhancing understanding of nanofluid behaviors in porous media. Finally, Krishna and
Chamkha [76] examined the MHD rotating flow of elastico-viscous fluids through a
porous medium, demonstrating how Hall and ion slip effects can alter the behavior
of viscoelastic fluids, further contributing to the study of complex fluid dynamics in
porous environments. Khan et al. [105] explored viscoelastic nanofluid motion in
the Homann stagnation-region, incorporating thermal radiation characteristics, thus
advancing the understanding of heat transfer in viscoelastic flows. Sarfraz and Khan
[106] analyzed entropy generation in CNT-based nanofluid flows induced by a moving
plate, contributing to the optimization of thermal systems using nanofluids. Sarfraz
and Khan [107] further examined the flow of ethylene glycol-based CNT Homann
nanofluids over a biaxially stretching surface, shedding light on the fluid dynamics
of nanofluids in complex configurations. The work of Sarfraz et al. [108] provided
insights into the thermodynamic efficiency of Homann-Agrawal hybrid nanofluids,
broadening the scope of hybrid nanofluid applications. Ahmed et al. [109] investigated
heat transfer in Jeffrey fluid flow over a power law lubricated surface, emphasizing
the influence of solar radiations and magnetic flux, while Khan et al. [110] examined
energy transport near the Homann stagnation point flow over a spiraling disk using the
Cattaneo–Christov theory, offering new perspectives on thermal transport mechanisms
48 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
in nanofluid systems. The analysis in this study aligns with the work by Muhammad
et al. [111], who explored the behavior of electrically conducting hybrid nanomaterials
near the stagnation region, focusing on how varying parameters affect fluid stability and
heat transfer. Similarly, this research investigates how thermal and mass stratification
influence MHD flow in porous media, using exact solutions derived via the Laplace
transform technique. Additionally, the study’s focus on stability parallels Sarfraz and
Khan’s [112] stability analysis of gyrotactic microorganisms in Jeffrey fluid flow, further
extending the understanding of fluid behavior in complex environments influenced by
stratification and magnetic fields. Both research streams contribute to improving the
efficiency of heat and mass transfer processes in engineered systems. Puneeth et al.
[113] conducted a theoretical analysis on the thermal characteristics of Ree–Eyring
nanofluid flow due to bioconvection, advancing the understanding of heat transfer in
bio-inspired systems. Irfan et al. [114] explored the performance of new mass flux
theory on Carreau nanofluid, offering insights into convective heat transport efficiency.
Hussain et al. [115] investigated hybrid nanoliquid flow over a curved surface, focusing
on chemical reaction and melting heat, contributing to the field of non-linear heat
transfer. Anwar et al. [116] studied MHD nanofluid flow through a Darcy medium with
thermal radiation, furthering knowledge in electromagnetic fluid flow. Anwar [117]
developed numerical simulations for fractional nonlinear viscoelastic flow problems,
laying the groundwork for future fractional fluid dynamics research. Anwar et al.
[118] analyzed convection and diffusion in viscoelastic fluid flow, emphasizing variable
thermal conductivity. Hussain et al. [119] provided insights into MHD Williamson
fluid flow over a convective wedge, considering Ohmic dissipation. Hussain et al. [120]
explored Casson fluid flow over a permeable wedge with chemical reaction, expanding
applications in mass transfer. Haider et al. [121] examined fractional Carreau fluid
with nanoparticles, enhancing thermal performance under magnetic fields. Irfan et al.
[122] modeled MHD bio-convection with gyrotactic microbes, highlighting the complex
interactions in nanofluids.
The study utilizes the Laplace transform technique to address the dynamics of
unsteady MHD parabolic flow, revealing limitations when dealing with arbitrary Prandtl
or Schmidt numbers. The technique leads to equations that become intractable unless
specific conditions are imposed, such as setting both the Prandtl number (P r) and
Schmidt number (Sc) to unity. This restricts the generality of the solutions, limiting
their applicability to only those scenarios where (P r = 1) and (Sc = 1). For more
complex cases with varying (P r) and (Sc), the study must rely on numerical methods,
highlighting a significant limitation in the analytical approach. Future work in this
research should focus on developing more generalized analytical methods that can
handle arbitrary Prandtl and Schmidt numbers without necessitating simplification.
This enhancement will allow the model to be applied to a broader range of conditions,
improving its utility in practical scenarios. Additionally, exploring different analytical
approaches could provide a more complete understanding of the complex interactions
in MHD parabolic flows, further refining the predictive capabilities of the model. The
stability criteria of a model’s solutions involve analyzing the behavior of the system
3.1 Introduction 49
under small perturbations. Stability is achieved when small disturbances diminish over
time, ensuring the solution returns to equilibrium.
The sensitivity analysis of key parameters in the study provides insights into how
thermal stratification (γ), mass stratification (ξ), Grashof numbers (Gr and Gc), mag-
netic parameter (M ), and Darcy number (Da) influence fluid dynamics, temperature,
and concentration profiles in a porous medium under magnetohydrodynamic (MHD)
flow. Thermal and mass stratification significantly impact velocity and stability. As γ
and ξ increase, both temperature and concentration gradients steepen, which dampens
fluid motion and stabilizes the flow. Higher stratification suppresses buoyancy forces,
reducing the overall velocity and temperature profiles. Grashof numbers (Gr and Gc)
represent thermal and solutal buoyancy forces. Higher values of Gr lead to increased
velocity due to stronger thermal buoyancy, while higher Gc tends to hinder the flow
by opposing the thermal buoyancy forces. The magnetic parameter (M ) introduces
Lorentz forces that oppose fluid motion, thus reducing velocity and temperature profiles.
This electromagnetic drag transforms kinetic energy into heat, which stabilizes the
system and limits convective heat transfer. On the other hand, the Darcy number
(Da), representing the permeability of the porous medium, has a direct effect on
fluid movement. Higher Da values allow the fluid to pass more freely through the
porous medium, enhancing both velocity and heat transfer by reducing flow resistance.
Thus, the combined effects of these parameters shape the behavior of the fluid system,
optimizing heat and mass transfer efficiency in various engineering applications.
This study introduces several innovative aspects: it employs the Laplace transform
technique to derive an exact solution without any approximation, proving to be a
highly effective method. Previously, there has been no research exploring the combined
effects of thermal and mass stratification on MHD unsteady parabolic flow along an
infinite vertical plate within a porous medium, especially under conditions of Periodic
Temperature Variation and Variable Mass Diffusion. The aim of this research was
to analytically explore the MHD unsteady parabolic flow in the presence of both
thermal and mass stratification, comparing the results of fluid dynamics in this scenario
against the classical case without stratification. This study thoroughly examines
and presents the effects of key physical parameters such as γ, ξ, Gr, Gc, the Magnetic
parameter M , and Darcy number (Da) on the flow profiles through detailed graphs and
tables. The insights gained have wide-ranging applications across various industries,
notably in enhancing the performance of thermal systems by minimizing temperature
discrepancies between the walls and the fluid and accelerating mass diffusion processes.
The real-world applicability of the findings from this study lies in optimizing fluid
behavior in engineering and environmental processes. The effects of thermal and
mass stratification on MHD flow in porous media can enhance the efficiency of heat
and mass transfer systems, such as in industrial cooling, pollutant management, and
energy system design. These findings also provide valuable insights into process control,
improving the prediction and stabilization of fluid motion in varied magnetic and
porous environments.
50 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
• Infinite Vertical Plate: The flow is considered along an infinite vertical plate,
meaning the flow variables depend only on the coordinate perpendicular to the
plate and time, not along the plate length.
• Variable Mass Diffusion: Concentration at the plate varies with time, signifi-
cantly affecting mass transfer characteristics.
• No Side Boundaries: The lack of side boundaries implies that the flow extends
indefinitely in the horizontal directions parallel to the plate.
• Porous Medium: The presence of a porous medium affects fluid flow and
diffusion characteristics, modeled through parameters like the Darcy number.
• MHD Effects: Magnetic fields are considered, affecting the fluid dynamics
through electromagnetic forces, represented by the magnetic parameter.
These assumptions form the basis for the mathematical modeling and analysis
conducted in the study, focusing on how these specific conditions impact flow dynamics,
heat, and mass transfer in the context of engineering and environmental applications.
Energy Equation
∂T ′
" #
ρCp + (⃗q.∇)T
⃗ ′ = α∇2 T ′ (3.2)
∂t′
Concentration Equation
∂C ′
′
+ (⃗q.∇)C
⃗ ′ = D∇2 C ′ (3.3)
∂t
52 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
∂u′ ∂ 2 u′ σB02 u′ ν ′
= gβ(T ′
− T∞
′
) + gβ ∗
(C ′
− C ′
∞ ) + ν − − u (3.4)
∂t′ ∂y ′2 ρ k
′ 2 ′
∂T ∂ T
′
= α ′2 − γ ′ u′ (3.5)
∂t ∂y
′
∂C ∂ 2C ′
= D − ξ ′ u′ (3.6)
∂t′ ∂y ′2
u′ = 0 T ′ = T∞′
C ′ = C∞
′
∀y ′ , t′ ≤ 0
u′ = u0 t′2 T ′ = Tω′ + (Tw′ − T∞
′
)cosω ′ t′ C ′ = C∞
′
+ (Cw′ − C∞
′
)At′ at y ′ = 0, t′ > 0
u′ = 0 T ′ → T∞ ′
C ′ → C∞ ′
as y ′ → ∞, t′ > 0
!1/3
gβ(Tw′ − T∞ ′
) gβ ∗ (Cw′ − C∞ ′
) ν ν ν
Gr = , Gc = , Pr = , Sc = , ω = ω′ ,
(νu0 )1/3 (νu0 ) 1/3 α D u20
!1/3 !1/3
σB02 ν u20 γ ′ν ξ′ν
M= , Da = k , γ= , ξ=
ρ u20 ν4 u0 (Tw′ − T∞
′ ) u0 (Cw′ − C∞
′ )
u2 1/3
A = ν0 is the constant. Then non-dimensional forms of the equations (3.4)-(3.6)
are given by
∂U ∂ 2U 1
= Grθ + GcC + 2
− M+ U (3.7)
∂t ∂y Da
∂θ 1 ∂ 2θ
= − γU (3.8)
∂t P r ∂y 2
∂C 1 ∂ 2C
= − ξU (3.9)
∂t Sc ∂y 2
d2 U 1
2
− sU + Grθ + GcC − M + U =0 (3.11)
dy Da
d2 θ
− sθ = γU (3.12)
dy 2
and
d2 C
− sC = ξU (3.13)
dy 2
Where the Laplace transform parameter is denoted by s. For example, the Laplace
transform of U (y, t) is given by
ˆ ∞
L{U (y, t)} = e−st U (y, t) dt = U (y, s),
0
such that its inverse Laplace transform exists and satisfies
1 Gr
U = {F h3 (F ) − Qh3 (Q)} + {h4 (Q, iω) + h4 (Q, −iω)
F −Q 2(F − Q)
Gr + Gc
−h4 (F, iω) − h4 (F, −iω)} − {h1 (F ) − h1 (Q)} (3.14)
F −Q
1
θ = (1 + E1 − E2 ) erf c(η) + ( + E3 − E4 ) [h5 (iω) + h5 (−iω)]
2
1 1
" (
Gr
+γ {h3 (F ) − h3 (Q)} + h4 (Q, iω)
F −Q 2(F − Q) Q
1 1 1 Gr + Gc
)
+ h4 (Q, −iω) − h4 (F, iω) − h4 (F, −iω) +
Q F F F −Q
1 1
( )#
h1 (Q) − h1 (F ) (3.15)
Q F
2η
" #
2
C = t (1 + 2η )erf c(η) − √ e−η − (D1 − D2 ) [h5 (iω) + h5 (−iω)]
2
π
ξ
+(D3 − D4 )erf c(η) + {h3 (F ) − h3 (Q)} − D1 {h4 (F, iω) + h4 (F, −iω)}
F −Q
+D2 {h4 (Q, iω) + h4 (Q, −iω)} − D3 h1 (F ) + D4 h1 (Q) (3.16)
3.4 Method of Solution 55
Where,
s
1 1
2
y
η= √ , F Q = γGr+ξGc, F +Q = M + , F −Q = M+ − 4(γGr + ξGc)
2 t Da Da
2η
( )
2
Cc = t 1 + 2η 2 erf c(η) − √ e−η (3.17)
π
1
θc = erf c(η) + [h5 (iω) + h5 (−iω)] (3.18)
2
Gr Gr
Uc = h3 (F + Q) − h1 (f + Q) − {h4 (F + Q, iω) + h4
F +Q 2(F + Q)
Gr Gr
(F + Q, −iω)} + erf c(η) + [h5 (iω) + h5 (−iω)] (3.19)
F +Q 2(F + Q)
56 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
3.4.2 Skin-Friction
The non-dimensional Skin-Friction, which is determined as shear stress on the surface,
is obtained by
dU
τ=−
dy y=0
The solution for the Skin-Friction is calculated from the solution of Velocity profile U ,
represented by (3.14), as follows:
√ √ √ −F t
√ √
s
F t erf ( F t) erf ( F t) te
τ= t2 F erf ( F t) + t e−F t − +t √ + √
2(F − Q) π 4F 3/2
F 2F π
Gr
√ q √ q
− e−iωt F − iω erf ( (F − iω)t) + eiωt F + iω erf ( (F + iω)t)
2(F − Q)
2e−F t 2e−Qt
q q #
q q
−e−iωt
Q − iωerf ( (Q − iω)t) − e iωt
Q + iω erf ( (Q + iω)t) + √ − √
πt πt
Gr + Gc √ √ −F t −Qt
" #
e q q e Q
− F erf ( F t) + √ − Q erf ( Qt) − √ −
F −Q πt πt 2(F − Q)
s √ √ √ −Qt
q q
t2 Q erf ( Qt) +
t erf ( Qt) erf ( Qt) te
t e−Qt − +t √ + √
π 4Q 3/2 Q 2Q π
The solution for the Skin-Friction for the classical case is given from the expression
(3.19), which is represented by
q
erf ( (F + Q)t)
s
1 2q q t
τc = t F + Q erf ( (F + Q))t) + t e−(F +Q)t −
2 π 4(F + Q)3/2
q √ −(F +Q)t
erf ( (F + Q)t) te Gr q
q
+t √ + √ − F + Q erf ( (F + Q)t)
F +Q 2(F + Q) π F +Q
e−(F +Q)t 1 √ √ √
#
Gr Gr h
+ √ + √ + e−iωt −iω erf ( −iωt) + eiωt iω
πt F + Q πt 2(F + Q)
√ 2
#
Gr
q
erf ( iωt) + √ − e−iωt F + Q − iω
πt 2(F + Q)
2e−(F +Q)t
q q #
q
erf ( (F + Q − iω)t) + e iωt
F + Q + iω erf ( (F + Q + iω)t) + √
πt
3.4 Method of Solution 57
dθ
Nu = −
dy y=0
The solution for the Nusselt number is calculated from the solution of Temperature
profile θ, represented by (3.15), as follows:
1 1 h √ √
N u = (1 + E1 − E2 ) √ + ( + E3 − E4 ) e−iωt −iω erf ( −iωt)
πt 2
√ √
s
√ √ 2
#
γ t
+eiωt iω erf ( iωt) + √ + t2 F erf ( F t) + t e−F t
πt 2(F − Q) π
√ √ √ −F t #
erf ( F t) erf ( F t) te √
q
− + t √ + √ − E 3 e −iωt
F − iω erf ( (F − iω)t)
4F 3/2 F 2F π
√ 2e −F t √ √ −F t
# " #
q e
+eiωt F + iω erf ( (F + iω)t) + √ − E1 F erf ( F t) + √
πt πt
−Qt
" #
e
q q q q
+E2 Q erf ( Qt) + √ + E4 e−iωt Q − iω erf ( (Q − iω)t)
πt
2e −Qt
#
γ
q q q q
+e iωt
Q + iωerf ( (Q + iω)t) + √ − t2 Q erf ( Qt)
πt 2(F − Q)
s √ √ √ −Qt
t erf ( Qt) erf ( Qt) te
+ t e−Qt − +t √ + √
π 4Q 3/2 Q 2Q π
The solution for the Nusselt number for the classical case is given from the expression
(3.18), which is represented by
2 1h √ √ √ √ i
N uc = √ + e−iωt −iω erf ( −iωt) + eiωt iω erf ( iωt)
πt 2
dC
Sh = −
dy y=0
The solution for the Sherwood number is calculated from the solution of Concentration
profile C, represented by (3.16), as follows:
58 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
s
√ √ √ √ 2
" #
t
Sh = 2 + (D1 − D2 ) e−iωt −iω erf ( −iωt) + eiωt iω erf ( iωt) + √
π πt
√
√ √
s
1 ξ t erf ( F t)
+(D3 − D4 ) √ + t2 F erf ( F t) + t e−F t −
πt 2(F − Q) π 4F 3/2
√ √ −F t #
erf ( F t) te √
q
+t √ + √ − D1 e−iωt
F − iω erf ( (F − iω)t) + eiωt
F 2F π
√ 2e−F t √ √ e−F t
q # " #
F + iω erf ( (F + iω)t) + √ − D3 F erf ( F t) + √ + D4
πt πt
e−Qt
" # q
q q q
Q erf ( Qt) + √ + D2 e −iωt
Q − iω erf ( (Q − iω)t) + eiωt
πt
2e−Qt
#
ξ
q q q q
Q + iω erf ( (Q + iω)t) + √ − t2 Q erf ( Qt)
πt 2(F − Q)
s √ √ √ −Qt
t erf ( Qt) erf ( Qt) te
+ t e−Qt − +t √ + √
π 4Q 3/2 Q 2Q π
The solution for the Sherwood number for the classical case is given from the
expression (3.17), which is represented by
s
t
Shc = 2
π
the fluid flow velocity. Therefore, the presence of both thermal and mass stratification
reduces fluid velocity compared to conditions without stratification, promoting a more
stable stratified flow. As seen in Fig. 3.4, increasing the thermal Grashof number (Gr)
enhances buoyant forces driven by temperature differences, accelerating fluid motion.
Conversely, increasing the mass Grashof number (Gc) generates stronger concentration
gradients, creating opposing buoyancy forces that hinder flow, leading to reduced
3.5 Results and Discussion 61
velocity. In Figure 3.5, increasing M further amplifies Lorentz forces, which oppose the
fluid’s natural motion. These electromagnetic forces enhance resistive drag, reducing
velocity by converting kinetic energy into heat, thus dampening convective flow and
suppressing fluid movement. We can observe from Fig. 3.6 that higher Da values
correspond to increased permeability in porous media, reducing flow resistance and
allowing fluid to pass more freely. This enhances velocity, as the porous medium offers
less impedance, facilitating smoother fluid movement and stronger convective currents.
62 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
In Fig. 3.7, increasing the thermal stratification parameter (γ) or mass stratification
parameter (ξ) enhances the vertical temperature and concentration gradients, leading
to stability in the fluid layers. These gradients reduce buoyancy-driven convection,
causing a suppression in velocity as the stratified fluid resists vertical mixing, thereby
dampening overall fluid motion. This stabilizing effect hinders convective flow, leading
to lower velocities over time.
The effect of thermal and mass stratification on the temperature profile is seen in
Fig. 3.8. As γ increases, it tends to lower the temperature gradient, suggesting that
thermal stratification reduces the rate at which temperature decreases with η. Similarly,
an increase in ξ also lowers the temperature gradient, implying that mass stratification
affects the thermal field by moderating temperature variations. Combinations of high
values of both γ and ξ show the most significant reduction in temperature changes
across η, indicating the compounded effect of both types of stratification. Figure 3.9
64 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
demonstrates that increasing the thermal Grashof number (Gr) reduces the temperature,
reflecting enhanced thermal buoyancy’s cooling effect. Increasing the mass Grashof
number (Gc) similarly lowers the temperature. Maximizing both parameters leads to
the greatest temperature decrease, due to the combined cooling effects of increased
thermal and mass buoyancy. Figure 3.10 illustrates that elevated M values increase
magnetic damping through stronger Lorentz forces, reducing kinetic-to-thermal energy
conversion, thus lowering temperature profiles by limiting convective heat transfer
within the fluid. Figure 3.11 illustrates that higher Darcy number (Da) values enhance
permeability in porous media, reducing flow resistance and improving heat transfer by
facilitating fluid movement. This leads to greater thermal diffusion, resulting in a flatter
temperature profile. Additionally, Figure 3.12 depicts the influence of thermal (γ) and
mass (ξ) stratification parameters on temperature over time. As γ, ξ, or both increase,
the temperature decreases more significantly throughout the timeline, highlighting
the role of stratification in thermal regulation. For example, lines representing higher
values of both γ and ξ show a greater decline in temperature compared to lines with
lower or zero values. This suggests that increased stratification intensifies the cooling
effect over time, leading to lower temperatures as the interaction between thermal and
mass influences becomes more pronounced.
As demonstrated in Figures 3.13, 3.14, 3.15 and 3.16, the impacts of γ and ξ, Gr
and Gc, M , and Da on concentration profiles mirror those observed in temperature
profiles. Fig. 3.17 shows that increasing both the thermal stratification parameter (γ)
and the mass stratification parameter (ξ) influences the concentration profiles over
time. Higher values of γ and ξ, individually or combined, lead to a slower increase
3.5 Results and Discussion 65
time, indicating improved heat and mass transfer, respectively. Higher γ and ξ values
further enhance these transfers, as seen in the steepness of the curves. This suggests
that increased stratification intensifies the heat and mass transfer efficiencies in the
fluid.
Table 3.1 presents the impact of thermal Grashof number (Gr), mass Grashof
number (Gc), magnetic parameter (M ), and Darcy number (Da) on skin friction,
Nusselt number (N u), and Sherwood number (Sh). As Gr and Gc increase, there is a
3.5 Results and Discussion 67
significant rise in skin friction, indicating that enhanced buoyancy forces lead to greater
convective flow and higher surface resistance. This also results in improved heat and
mass transfer, as reflected by the higher N u and Sh values. Conversely, increasing the
magnetic parameter (M ) reduces skin friction, due to the magnetic field’s damping
effect, which slows fluid motion by introducing Lorentz forces. This also slightly
decreases the Nusselt and Sherwood numbers, as the stabilizing effect of the magnetic
68 Dynamic Interaction of Stratifications in a MHD Parabolic Flow...
field reduces convective heat and mass transfer. The Darcy number (Da) shows a
more complex relationship: as permeability increases, skin friction initially decreases
due to restricted flow but rises again at higher Da values as fluid motion becomes less
restricted. Both Nu and Sh increase with Da, reflecting enhanced heat and mass transfer
as permeability increases. Overall, thermal and mass stratification, magnetic fields,
and porous medium characteristics significantly influence fluid dynamics, affecting both
heat and mass transfer efficiency.
3.6 Conclusion 69
Table 3.1 The impact of different parameters on Skin Friction, Nusselt number and
Sherwood number coefficients with
Gr Gc M Da τ Nu Sh
5 5 0.3 0.5 -3.2582 1.2030 1.3793
10 5 -6.1775 1.2240 1.3960
5 10 -5.0863 1.3208 1.4735
10 10 -8.0389 1.3319 1.4824
0.1 -3.3470 1.2080 1.3833
0.5 -3.1721 1.1983 1.3755
1 -2.9683 1.1873 1.3667
0.3 -2.7331 1.1746 1.3566
0.7 -3.5186 1.2180 1.3912
1 -3.7265 1.2310 1.4017
3.6 Conclusion
The study reveals that increasing thermal and mass stratification parameters consis-
tently reduces velocity and temperature, stabilizing fluid flow and thermal environments.
Enhancements in thermal and mass Grashof numbers improve heat and mass trans-
fer efficiency, as indicated by higher Nusselt and Sherwood numbers. The magnetic
parameter and Darcy number significantly impact heat transfer and flow resistance,
respectively. The novel findings demonstrate that combined thermal and mass strati-
fication dampens fluid velocity, enhances heat and mass transfer, and increases skin
friction due to stratification-induced stability. These insights provide valuable guidance
for optimizing fluid dynamics in engineering applications. Future studies could inves-
tigate non-linear temperature and concentration gradients, variable magnetic fields,
and different porous media to further improve process control. Practical applications
include optimizing industrial cooling systems, controlling environmental pollutants,
and designing efficient heat exchangers, highlighting the relevance of these findings in
engineering and energy systems.
The Comprehensive Analysis of
Stratification and Chemical
4
Reaction on MHD Flow along an
Impulsively Started Infinite Vertical
Plate in Presence of Radiation and
Heat Source Through Porous
Medium
Published in Heat Transfer, 2025, ISSN: 2312-4334, E-ISSN: 2688-4542,
(SCOPUS, ESCI)
4.1 Introduction
In fluid dynamics and heat transfer studies, the interaction of chemical reaction, heat
sources, porous media, thermal radiation, and thermal stratification creates a complex
environment where each factor influences the others. Chemical reaction within a flow
generate or consume heat, altering temperature profiles and often interacting with
heat sources that inject additional energy, intensifying thermal gradients and driving
convection processes. This intricate temperature distribution becomes more complex
in a porous medium, which restricts fluid movement while enhancing thermal and mass
exchange through increased surface area. The laplace transform technique is frequently
employed to simplify these interactions mathematically, converting differential equations
72 The Comprehensive Analysis of Stratification and Chemical Reaction...
the influence of chemical reaction, heat generation, and various stratification effects.
Harikrushna et al. [136] analyzed chemical reaction effects on MHD flow past a vertical
absorbent plate with ramped wall temperature and concentration. Sudarmozhi et
al. [137] investigated Maxwell fluid dynamics under chemical reaction and viscous
dissipation. Appidi et al. [138] examined the influence of thermal radiation and
viscous dissipation on MHD convective flows with micro-polar fluids. Rai and Mishra
[139] explored MHD nanofluid flows with slip conditions and radiation effects on a
moving plate. Ramanjaneyulu and Hari Krishna [140] studied heat absorption in fluids
around inclined porous plates. Gopal et al. [141] assessed second-grade fluid flow in
unsteady MHD convection with chemical reaction. Chandra Sekhar and Vijaya Kumar
[142] investigated chemical impacts on unsteady MHD flows past infinite vertical
plates, while Reddy et al. [143, 144] provided computational analyses on chemically
reactive magneto-nanofluids in oscillating and permeable plates. Sahu and Deka [145–
148, 123, 149, 150] analyzed the dynamic effects of thermal and mass stratification on
MHD parabolic flow with periodic temperature variations, adding valuable insights
into stratification influences in porous media. Khan [151] explored the significance of
ramped temperature in the unsteady dynamics of viscoelastic fluids under Lorentz
force, demonstrating how thermal ramping can impact flow behavior and heat transfer
characteristics. This study is crucial for understanding fluid behavior in magnetic
fields, with applications in industrial processes and biomedical engineering. In a com-
plementary study, Khan [152] investigated the Prabhakar fractional derivative model
for sodium alginate (C6 H9 N aO7 ), highlighting its role in modeling accelerated plate
motions. This work emphasizes the relevance of fractional derivatives in describing the
memory effects in viscoelastic materials. The hemodynamics of blood flow, a critical
aspect of biomedical engineering, has also been addressed. Khan [153] examined blood
flow over an inclined cylinder, contributing to the understanding of blood dynamics
under complex geometrical and external conditions. This study offers insights into
physiological flows, such as those in arterial structures. Khan’s [154] earlier work
introduced the Atangana-Baleanu fractional derivatives to model human blood flow
in nanofluids, shedding light on how fractional calculus can refine the understanding
of rheological properties in biological systems. This model has implications for drug
delivery systems and other biomedical applications. Building on these findings, Khan
[155] analyzed ramped heating in CNTs fractional nanofluids, revealing the impact of
ramped temperature on heat transfer efficiency. This work holds promise for enhancing
thermal management in nanotechnology applications. The entropy generation and
thermal analysis in viscoelastic fluids are vital for optimizing energy systems. Khan
[156] explored the entropy generation in sodium alginate fluids subjected to MHD and
porosity effects, providing insights into energy efficiency and thermal performance in
engineering systems. Furthermore, Khan et al. [157], presented a natural convection
simulation for Prabhakar-like fractional Maxwell fluid on inclined planes, incorporating
generalized thermal flux. Their findings highlight the interplay of fractional derivatives
and natural convection, with implications for materials processing and energy systems.
Khan’s extensive contributions, ranging from fractional derivatives to MHD and ther-
74 The Comprehensive Analysis of Stratification and Chemical Reaction...
mal effects in various fluid systems, underscore the growing importance of advanced
mathematical models in addressing complex fluid behaviors across engineering and
biomedical domains. Kataria and Mittal [158, 159] investigated nanofluid flow over an
oscillating vertical plate. Sheikholeslami et al. [160] examined the effects of radiation on
three-dimensional nanofluid heat transfer, incorporating thermal interfacial resistance
and suspension micromixing. Using a mesoscopic approach, Li et al. [161] analyzed the
MHD flow of a nanofluid through a permeable duct. Additionally, Patel et al. [162]
studied the effect of magnetic fields on the flow of micropolar ferrofluids. Mittal and
Patel [163] investigated the flow of Casson fluid at the stagnation point over an infinite
plate within a media that is permeable in two dimensions, using a mixed convective
MHD framework. Kataria and Mittal [164] analyzed the flow of Casson nanofluid in
the existence of radiation and a magnetic field. Mittal and Kataria [165] studied the
three-dimensional flow of CuO–Water nanofluid through a porous medium within a
rotating system confined between two horizontal parallel plates. Furthermore, Kataria
and Mittal [166] provided a mathematical analysis of three-dimensional nanofluid flow
in a rotating system, accounting for thermal interfacial resistance in suspensions within
porous media. Mittal et al. [167] investigated the influence of phenomena such as
viscous dissipation, mixed convection, and Joule heating on micropolar ferrofluid flow.
Mittal [168] analyzed the mass and heat transfer in a three-dimensional composite MHD
nanofluid flow between two parallel horizontal plates within a rotating system. Mittal
[169] examined the effect of a magnetic field on the flow of a two-dimensional squeezing
nanofluid confined between two equally spaced plates. Kataria et al. [170, 171] studied
the effects of nonlinear radiation on MHD micropolar fluid flow and entropy-optimized
MHD fluid flow, respectively. Patel et al. [172] explored nanofluid flow in unsteady
magnetohydrodynamics with thermo-diffusion and heat-generation effects over an
exponentially accelerated porous plate. The time-dependent flow across a vertical plate
with a stepped surface concentration in a porous medium was investigated by Mittal et
al. [173]. Additionally, Patel et al. [174] analyzed the effects of radiation and magnetic
fields on unsteady mixed convection micropolar nanofluid flow.
This study presents several novel contributions. The Laplace transform technique is
utilized to obtain an exact solution with no approximations, showcasing its effectiveness
as a reliable analytical approach. The main objective of this research was to conduct an
analytical investigation to examine MHD unsteady flow with thermal stratification and
compare the resulting fluid dynamics with the classical case without stratification. Key
physical parameters, including γ, K, Gr, Gc, R, ϕ, Sc, M , and Da, were meticulously
analyzed, with comprehensive graphs and tables illustrating their effects on flow profiles.
The findings offer significant applications in optimizing thermal systems, reducing
temperature disparities between walls and fluids, and improving mass diffusion processes.
These results are invaluable for engineering and environmental applications, improving
heat and mass transfer efficiency in systems like industrial cooling, pollution control,
and energy optimization. Additionally, the study provides critical insights for the
4.2 Mathematical Formulation 75
control of processes, facilitating fluid motion prediction and stabilization in porous and
magnetic mediums.
∂u′ ∂ 2 u′ σB02 u′ ν ′
= gβ(T ′
− T∞
′
) + gβ ∗
(C ′
− C ′
∞ ) + ν − − u (4.1)
∂t′ ∂y ′2 ρ k
∂T ′ k ∂ 2T ′ 1 ∂qr q0
′
= ′2
− ′
− (T ′ − T∞′
) − γ ′ u′ (4.2)
∂t ρCp ∂y ρCp ∂y ρCp
′ 2 ′
∂C ∂ C
′
= D ′2 − K ′ (C ′ − C∞ ′
) (4.3)
∂t ∂y
76 The Comprehensive Analysis of Stratification and Chemical Reaction...
u′ = 0 T ′ = T∞′
C ′ = C∞′
∀y ′ , t′ ≤ 0
u′ = u0 T ′ = Tw′ C ′ = Cw′ at y ′ = 0, t′ > 0
u′ = 0 T ′ → T∞ ′
C ′ → C∞ ′
as y ′ → ∞, t′ > 0
′ ′
Here, the term γ ′ ≡ dT ∞
dx′
+ Cgp , where dT ∞
dx′
represents thermal stratification or
vertical temperature advection. Additionally, Cp denotes the rate at which compression,
g
sometimes referred to as the work of compression, performs reversible work on the fluid
particles. For the purpose of evaluating numerical models, the compression work is
combined with thermal stratification as an additive component.
T′ ∼
4
= 4T∞
′3
T ′ − 3T∞
′4
Consequently, we get
∂T ′ k ∂ 2T ′ 1 h ∗ ∗ ′3 ′ i q0
= 2 − 16a σ T∞ (T − T∞
′
) − (T ′ − T∞
′
) − γ ′ u′
∂t′ ρCp ∂y ′ ρCp ρCp
non-dimensional quantities are presented in the following section:
!1/3
u20 T ′ − T∞
′
C ′ − C∞
′
1/3 1/3
u0 u0
U =u ′
, t=t ′
, y = y′ , θ= , C= ,
ν2 ν ν2 Tw′ − T∞′ Cw′ − C∞′
!1/3
gβ(Tw′ − T∞ ′
) gβ ∗ (Cw′ − C∞′
) µCp ν q0 ν
Gr = , Gc = , Pr = , Sc = , ϕ= ,
(νu0 )1/3 (νu0 )1/3 k D ρCp u20
ν
K = K′
u20
4.3 Solution Method 77
u2 1/3
A = ν0 is the constant. The dimensionless forms of equations (4.1)-(4.3) are
expressed as follows:
∂U ∂ 2U 1
= Grθ + GcC + − M + U (4.4)
∂t ∂y 2 Da
∂θ 1 ∂ 2T 1
= − Rθ − ϕT − γu (4.5)
∂t P r ∂y 2 Pr
∂C 1 ∂ 2C
= − KC (4.6)
∂t Sc ∂y 2
The non-dimensional forms of the initial and boundary conditions are outlined as
follows:
Gc
U = [E4 ξ(A) + E5 ξ(A, Z1 , Z2 ) + E6 ξ(A, Z1 , −Z2 )]
A−B
A − N + Gr + Gc Gc
+ ξ(A) +
A−B Sc − 1
[E1 ξ(B) + E2 ξ(B, Z1 , Z2 ) + E3 ξ(B, Z1 , −Z2 )]
A − N + Gr + Gc Gc
+ ξ(B) −
A−B A−B
Gc
[E4 ξ(B) + E5 ξ(B, Z1 , Z2 ) + E6 ξ(B, Z1 , −Z2 )] −
Sc − 1
[E1 ψ(K) + E2 ψ(K, Z1 , Z2 ) + E3 ψ(K, Z1 , −Z2 )] (4.8)
78 The Comprehensive Analysis of Stratification and Chemical Reaction...
C = ψ(K) (4.10)
Where,
y 1 1
η= √ , AB = M + (R + ϕ) − γGr, A+B =M + + R + ϕ,
2 t Da Da
s
1 1
2
A−B = M+ −R−ϕ − 4γGr, N =M+ , G = R + ϕ,
Da Da
v
2 Sc K − N − G u (N − G)2 − 4γGr Sc K − G
u
Z1 = , Z2 = t , Z3 = ,
2(Sc − 1) 4(Sc − 1)2 Sc − 1
Sc K − B Sc K − N Z3 Z3 − Z1 − Z2
Z4 = , Z5 = , E1 = , E2 = ,
Sc − 1 Sc − 1 Z12 − Z22 2Z2 (Z1 + Z2 )
Z2 − Z1 + Z3 Z3 Z4 (Z3 − Z1 − Z2 )(Z4 − Z1 − Z2 )
E3 = , E4 = , E5 = ,
2Z2 (Z2 − Z1 ) Z12 − Z22 2Z2 (Z1 + Z2 )
(Z2 − Z1 + Z3 )(Z2 − Z1 + Z4 ) γ Gc E4 γ Gc E5
E6 = , λ1 = , λ2 = ,
2Z2 (Z2 − Z1 ) (A − G)(B − G) (A − G)(B − G)
γ Gc E6 γ (N − B − Gr − Gc) γ Gc E1
λ3 = , λ4 = , λ5 = ,
(A − G)(B − G) (A − G)(B − G) (Sc − 1)(B − G)
γ Gc E2 γ Gc E3 γ Gc E1
λ6 = , λ7 = , λ8 = ,
(Sc − 1)(B − G) (Sc − 1)(B − G) Z3 (Sc − 1)2
γ Gc E2 γ Gc E3 γ Gc E4
λ9 = , λ10 = , λ11 = ,
(Sc − 1)2 (Z1 + Z2 − Z3 ) (Sc − 1)2 (Z1 − Z2 − Z3 ) (A − B)(A − G)
γ Gc E5 γ Gc E6 γ Gc E4
λ12 = , λ13 = , λ14 = ,
(A − B)(A − G) (A − B)(A − G) (A − B)(B − G)
γ Gc E5 γ Gc E6 γ (A − N + Gr + Gc)
λ15 = , λ16 = , λ17 = ,
(A − B)(B − G) (A − B)(B − G) (A − B)(B − G)
4.3 Solution Method 79
γ Gc E1 γ Gc E2
λ18 = , λ19 = ,
Z3 (Sc − 1)2 (Sc − 1)2 (Z1 + Z2 − Z3 )
γ Gc E3 Gr Gc
λ20 = , λ21 = , λ22 =
(Sc − 1)2 (Z1 − Z2 − Z3 ) N −G Z5 (Sc − 1)
Additionally, the inverse Laplace transforms, are defined as follows:
√ √
e−y s+A e−y s+A
( ) ( )
ξ(A, Z1 , Z2 ) = L −1
, ξ(A, Z3 ) = L−1 ,
s − (−Z1 − Z2 ) s − (−Z3 )
√ √
e−y Sc(s+K) e−y Sc(s+K)
ψ(K, Z1 , Z2 ) = L−1 , ψ(K, Z3 ) = L−1
s − (−Z1 − Z2 ) s − (−Z3 )
θc = ξ(G) (4.11)
Uc = (1 − λ21 ) ξ(N ) + λ22 [ξ(N ) − ξ(N, Z5 )]
+λ21 ξ(G) − λ22 [ξ(K) − ξ(K, Z5 )] (4.12)
4.3.2 Skin-Friction
The non-dimensional Skin-Friction, representing the shear stress at the surface, is
determined as follows:
dU
τ=−
dy y=0
The Skin-Friction is determined from the Velocity profile U , given by (4.8), and is
expressed as follows:
80 The Comprehensive Analysis of Stratification and Chemical Reaction...
E4 Gc + N − B − Gr − Gc √ √ e−At
" # !
τ = A erf ( At) + √
A−B πt
−At
" !
Gc q q e
+ E5 e(−Z1 −Z2 )t A − Z1 − Z2 erf ( (A − Z1 − Z2 )t) + √
A−B πt
−At
!#
q q e
+E6 e(−Z1 +Z2 )t A − Z1 + Z2 erf ( (A − Z1 + Z2 )t) + √
πt
−E4 Gc + A − N + Gr + Gc √ √ e−Bt
" # !
E1 Gc
+ + B erf ( Bt) + √
Sc − 1 A−B πt
e−Bt
" # !
E2 Gc E5 Gc q q
+ − e (−Z1 −Z2 )t
B − Z1 − Z2 erf ( (B − Z1 − Z2 )t) + √
Sc − 1 A − B πt
−Bt
" # !
E3 Gc E6 Gc q q e
+ − e (−Z1 +Z2 )t
B − Z1 + Z2 erf ( (B − Z1 + Z2 )t) + √
Sc − 1 A − B πt
s
√ √
" !
Gc Sc −Kt
− E1 K Sc erf ( Kt) + e
Sc − 1 πt
s !
q q Sc −Kt
+E2 e (−Z1 −Z2 )t
Sc (K − Z1 − Z2 ) erf ( (K − Z1 − Z2 )t) + e
πt
s !#
q q Sc −Kt
+E3 e (−Z1 +Z2 )t
Sc (K − Z1 + Z2 ) erf ( (K − Z1 + Z2 )t) + e
πt
The Skin-Friction solution for the classical case is obtained from the expression
(4.12) and is expressed as follows:
√ √ e−N t
!
τc = (1 − λ21 + λ22 ) N erf ( N t) + √
πt
e−N t
q q !
−λ22 e N − Z5 erf ( (N − Z5 )t) + √
−Z5 t
πt
√ √ e−Gt
!
+λ21 G erf ( Gt) + √
πt
s
√ √
" !
Sc −Kt
−λ22 K Sc erf ( Kt) + e
πt
s !#
q q Sc −Kt
−e−Z5 t Sc (K − Z5 ) erf ( (K − Z5 )t) + e
πt
4.3 Solution Method 81
dθ
Nu = −
dy y=0
The Nusselt number is derived from the Temperature profile θ, given by (4.9), and is
expressed as follows:
√ √ e−Gt
!
N u = (1 + λ1 − λ4 − λ5 + λ8 ) G erf ( Gt) + √
πt
e−Gt
q q !
+(λ2 − λ6 − λ9 ) e G − Z1 − Z2 erf ( (G − Z1 − Z2 )t) + √
(−Z1 −Z2 )t
πt
e−Gt
q q !
+(λ3 − λ7 − λ10 ) e (−Z1 +Z2 )t
G − Z1 + Z2 erf ( (G − Z1 + Z2 )t) + √
πt
−Gt
!
q q e
+(−λ8 + λ9 + λ10 ) e−Z3 t G − Z3 erf ( (G − Z3 )t) + √
πt
√ √ −At
!
e
+(λ11 + λ4 ) A erf ( At) + √
πt
e−At
q q !
+λ12 e (−Z1 −Z2 )t
A − Z1 − Z2 erf ( (A − Z1 − Z2 )t) + √
πt
−At
!
q q e
+λ13 e (−Z1 +Z2 )t
A − Z1 + Z2 erf ( (A − Z1 + Z2 )t) + √
πt
√ √ −Bt
!
e
+(λ5 − λ14 + λ17 ) B erf ( Bt) + √
πt
e−Bt
q q !
+(λ6 − λ15 ) e (−Z1 −Z2 )t
B − Z1 − Z2 erf ( (B − Z1 − Z2 )t) + √
πt
−Bt
!
q q e
+(λ7 − λ16 ) e (−Z1 +Z2 )t
B − Z1 + Z2 erf ( (B − Z1 + Z2 )t) + √
πt
s !
q q Sc −Kt
+(−λ18 − λ19 − λ20 )e (−Z3 )t
Sc (K − Z3 ) erf ( (K − Z3 )t) + e
πt
s
√ √
!
Sc −Kt
−λ18 E1 K Sc erf ( Kt) + e
πt
s !
q q Sc −Kt
+λ19 e (−Z1 −Z2 )t
Sc (K − Z1 − Z2 ) erf ( (K − Z1 − Z2 )t) + e
πt
s !
q q Sc −Kt
+λ20 e (−Z1 +Z2 )t
Sc (K − Z1 + Z2 ) erf ( (K − Z1 + Z2 )t) + e
πt
82 The Comprehensive Analysis of Stratification and Chemical Reaction...
The Nusselt number solution for the classical case is obtained from the expression
(4.11) and is presented as follows:
√ √ e−Gt
N uc = G erf ( Gt) + √
πt
dC
Sh = −
dy y=0
The Sherwood number is derived from the Concentration profile C, given by (4.10),
and is expressed as follows:
√
√ √ Sc e−Kt
Sh = K Sc erf Kt + √
πt
Figure 4.2 illustrates that the results of the present study closely align with those
of Muthucumaraswamy [175], demonstrating strong agreement; as both studies exhibit
a similar trend. The strong alignment of the curves validates the methodology and
assumptions used in the present study against an established study, strengthening its
reliability.
Table 4.1 shows a remarkable agreement between the Sherwood number (Sh)
derived in this study and those reported by Muthucumaraswamy [175] for K = 0.6 and
Sc = 0.78, confirming the reliability of the present results.
The effect of thermal stratification is depicted in Figure 4.3, which significantly
influences fluid velocity, as shown in the graph. Higher thermal stratification (γ)
reduces velocity due to increased thermal resistance, limiting energy transfer across
layers. As γ increases, the peak velocity decreases, and the velocity profile becomes
more uniform. When γ = 0, the absence of stratification allows maximum velocity,
emphasizing the role of thermal gradients in controlling flow behavior. These findings
align with those of [5]. This indicates a reduction in fluid velocity under thermal
stratification for γ > 0, as higher thermal stratification (γ) diminishes the convective
interaction between the heated plate and the fluid, leading to a reduction in buoyancy
force and, subsequently, flow velocity. As seen in Figure 4.4, chemical reaction signifi-
cantly impact fluid velocity. An increase in (K) leads to a decrease in fluid velocity
because reaction rates suppress flow momentum. When K is low, velocity is higher,
indicating minimal reaction effects. Conversely, higher K leads to reduced velocity,
showcasing the dominant role of chemical reaction in resisting flow and altering the
velocity profile.
Figure 4.5 illustrates the impact of the Grashof number (Gr) and modified Grashof
number (Gc) on the velocity profile (U ). As Gr and Gc increase, the velocity profile
rises, with the peak velocity becoming higher and occurring closer to the boundary
layer. This behavior occurs because larger Gr and Gc enhance the buoyancy effects,
resulting in a stronger upward force and increasing the flow velocity. The velocity
gradually decreases beyond the peak due to viscous effects, eventually stabilizing at
84 The Comprehensive Analysis of Stratification and Chemical Reaction...
the free-stream region. Figure 4.6 depicts the effect of the magnetic parameter (M )
on the velocity profile (U ). As M increases, the velocity profile decreases throughout
the boundary layer, with the peak velocity becoming lower. This reduction occurs
because the presence of a magnetic field induces a Lorentz force, which opposes the
fluid motion and increases the resistive drag in the flow. Consequently, the velocity
diminishes with higher magnetic field strength, demonstrating the suppressive effect of
M on fluid flow.
Figure 4.7 demonstrates the influence of the Darcy number (Da) on the velocity
profile (U ). As Da increases, the velocity profile becomes more pronounced, with higher
4.4 Results and Discussion 85
Figure 4.9 illustrates the effect of the volume fraction (ϕ) on the velocity profile (U ).
As ϕ increases, the velocity decreases due to the higher solid particle concentration in
the fluid, which enhances resistance and dampens the fluid motion. This reduction in
velocity is most prominent near the boundary layer, where the influence of particle
concentration is strongest. Figure 4.10 shows the influence of the Schmidt number (Sc)
on the velocity profile. An increase in Sc leads to a significant reduction in U , as higher
Sc corresponds to a lower molecular diffusivity of the solute. This suppresses the mass
transfer within the fluid, reducing the buoyancy-driven flow velocity. Consequently,
lower Sc values exhibit higher velocities due to enhanced mass diffusion.
4.4 Results and Discussion 87
leading to a thicker thermal boundary layer and higher overall fluid temperatures,
reflecting the impact of intensified chemical activity.
Table 4.2 presents the impact of various parameters on the temperature profile, θ.
An increase in the Grashof number (Gr) decreases the temperature, while a rise in
the modified Grashof number (Gc) also lowers the temperature profile. The magnetic
parameter (M ) exhibits a direct relationship, with higher values resulting in an increase
in θ. The Darcy number (Da) shows an inverse effect, where increasing Da reduces the
temperature. Similarly, the radiation parameter (R) significantly decreases θ, indicating
4.4 Results and Discussion 89
Gr Gc M Da R ϕ Sc θ
5 10 1.2 0.25 1 0.5 0.6 0.3456
15 0.3377
5 0.3862
15 0.2976
0 0.3315
1 0.3551
3 0.3861
5 0.4059
0.2 0.3907
0.4 0.3530
0.7 0.3266
1 0.3127
0 0.4008
3 0.3191
5 0.2815
10 0.2146
0 0.3836
5 0.2733
10 0.2094
15 0.1664
0.22 0.2824
0.6 0.3269
0.78 0.3419
1 0.3575
that higher radiation reduces the thermal effect. The volume fraction (ϕ) has a notable
negative influence on the temperature profile, with higher volume fractions leading to
substantial reductions in θ. In contrast, the Schmidt number (Sc) positively affects
the temperature profile, as larger values correspond to higher θ. These trends reveal
the complex dependencies of θ on the governing parameters.
In Figure 4.13 the influence of the stratification parameter (γ) on skin friction (τ )
over time (t). As γ increases, τ becomes less negative initially and stabilizes more
rapidly, indicating that higher stratification reduces the magnitude of shear stresses
caused by thermal and mass effects. This behavior is attributed to the reduced gradients
in temperature and concentration due to increased stratification. In Figure 4.14, the
effect of γ on the Nusselt number (N u) over time is shown. Higher γ values enhance N u,
especially at early times, due to increased heat transfer caused by stratification. Over
time, N u stabilizes at higher values for greater γ, reflecting the persistent enhancement
90 The Comprehensive Analysis of Stratification and Chemical Reaction...
of thermal energy transfer in stratified conditions. These results suggest that strat-
ification mitigates shear stresses while improving heat transfer efficiency within the flow.
Table 4.3 highlights the effects of various parameters on skin friction (τ ) and the
Nusselt number (N u). An increase in the permeability parameter (K) generally reduces
τ , transitioning from negative to positive values, while N u increases slightly. Variations
in (Gr) and (Gc) significantly influence both τ and N u. Higher Gr and Gc values
tend to decrease τ and slightly increase N u. An increase in the magnetic parameter
(M ) positively impacts τ , moving it from negative to positive, and also enhances N u.
4.4 Results and Discussion 91
Table 4.3 The influence of various parameters on Skin Friction and the Nusselt number
K Gr Gc M Da R ϕ Sc τ Nu
0.2 10 10 0.4 0.5 2 1.5 0.78 -0.4408 2.5466
0.8 -0.2907 2.5825
1.5 -0.0113 2.6151
3 2.6766 2.5475
0.6 5 5 0.4 0.5 2 1.5 0.78 0.6125 2.4069
10 5 0.6428 2.3904
5 10 -0.6125 2.6436
15 15 -1.1417 2.6821
0.6 10 10 0 0.5 2 1.5 0.78 -0.4477 2.5723
1 -0.2249 2.5727
3 0.0757 2.5957
5 0.3136 2.6261
0.6 10 10 0.4 0.4 2 1.5 0.78 -0.2439 2.5721
0.7 -0.4958 2.5739
1 -0.6351 2.5812
0.6 10 10 0.4 0.5 0 1.5 0.78 0.7642 2.4809
1 0.1512 2.5156
5 -1.4337 2.8154
10 -2.4238 3.6441
0.6 10 10 0.4 0.5 2 0 0.78 0.4408 2.4953
5 -1.5734 2.8640
10 -2.4196 3.7093
0.6 10 10 0.4 0.5 2 1.5 0.22 -1.9116 2.9481
0.6 -0.7584 2.6697
0.78 -0.3707 2.5708
1 0.0924 2.4644
The Darcy number (Da) shows a decreasing influence on τ , while its effect on N u is
marginal. Radiation parameter (R) significantly impacts both τ and N u, with higher
R values reducing τ and increasing N u. Finally, increasing the volume fraction (ϕ)
and Schmidt number (Sc) leads to a decrease in τ and an enhancement in N u.
Table 4.4 shows that the Sherwood number (Sh) increases with higher Schmidt
number (Sc) and permeability parameter (K), indicating enhanced mass transfer rates.
A noticeable rise in Sherwood number is observed for both parameters at their respective
higher values.
92 The Comprehensive Analysis of Stratification and Chemical Reaction...
K Sc Sh
0.2 0.78 1.3291
0.8 1.4417
1.5 1.5691
3 1.8284
0.6 0.22 0.7459
0.6 1.2319
0.78 1.4046
1 1.5904
4.5 Conclusion
In conclusion, the study provides a comprehensive analysis of the effects of various
parameters on velocity, temperature, concentration, skin friction, Nusselt number,
and Sherwood number profiles in a fluid flow system. Thermal stratification (γ) and
chemical reaction (K) significantly influence velocity, temperature, and mass transfer,
with γ reducing velocity and temperature due to increased thermal resistance, while
K enhances temperature and mass transfer through exothermic reaction. Buoyancy
forces, represented by thermal and mass Grashof numbers (Gr and Gc), drive the flow
and increase heat and mass transfer, but their effects diminish at higher values due to
reduced gradients. The magnetic parameter (M ) enhances skin friction and Nusselt
number by restricting fluid motion, while the Darcy number (Da) reduces these profiles
by facilitating smoother flow in a porous medium. Radiation (R) and heat source
(ϕ) reduce heat transfer efficiency by dispersing thermal energy, while the Schmidt
number (Sc) enhances Sherwood number and slightly impacts other profiles by limiting
mass diffusivity. These findings emphasize the interplay of thermal, mechanical, and
chemical forces in shaping fluid flow dynamics and provide insights into optimizing
heat and mass transfer in engineering and industrial applications.
Influences of Thermal Stratification
and Chemical Reaction on MHD
5
Free Convective Flow Along an
Accelerated Vertical Plate with
Variable Temperature and
Exponential Mass Diffusion in a
Porous Medium
Published in Heat Transfer, 53(7) (2024), 3643-3666, E-ISSN: 2688-4542,
ISSN: 2688-4534. (SCOPUS, ESCI)
5.1 Introduction
The widespread occurrence of thermal stratification in fluid systems arises from the
inherent tendency of fluids to organize themselves into distinct temperature layers based
on their thermal properties. This temperature stratification within a fluid medium
plays a pivotal role in defining the thermal characteristics of various natural and
man-made environments. Thermal stratification exerts a substantial influence on heat
transfer, energy efficiency, and environmental dynamics, with particular prevalence in
bodies of water, atmospheric conditions, and industrial processes. In the realm of fluid
dynamics, magnetohydrodynamics (MHD) serves as a pivotal field of study, focusing
on how magnetic fields interact with electrically conducting fluids such as plasmas
94 Influences of Thermal Stratification and Chemical Reaction on MHD...
and liquid metals. This interaction is crucial for manipulating fluid behaviors in a
variety of engineering and scientific processes. Concurrently, unsteady flow dynamics,
characterized by temporal variations in flow velocity and pressure, are essential for
understanding fluid behaviors in conditions marked by fluctuating dynamics, like those
observed in pulsatile flows. Moreover, the influence of chemical reaction extends to
modifying key fluid properties—including viscosity and density—which subsequently
shape the flow dynamics critical to the design of chemical reactors and the management
of environmental pollutants. The study of fluid flow through porous medium further
enriches our understanding, as the porous structure inherently modifies flow resistance,
making it indispensable across diverse applications from groundwater remediation to
industrial filtration. Additionally, examining the fluid dynamics near an accelerated
vertical plate sheds light on fluid responses to abrupt boundary movements or velocity
shifts, a phenomenon of significant importance in the engineering of rotating machinery
such as turbines and pumps.
Early work by Siegel [7] and Menold and Yang [8] laid the groundwork for un-
derstanding transient free convection from vertical flat plates. Schetz and Eichhorn
[9] extended this by investigating unsteady natural convection near doubly infinite
vertical plates. Later studies explored diverse aspects of this phenomenon. Park and
Hyun [13] and Park [14] investigated transient buoyant flows in stratified fluids, while
Shapiro and Fedorovich [5] explored convectively driven flows in stably stratified fluids.
Porous media were explored in greater detail. Magyari et al.[16] investigated unsteady
free convection in a fluid-saturated porous medium, and Chaudhary and Jain [72]
considered heat and mass diffusion flow by natural convection in a porous medium.
The influence of radiation and magnetic fields was also examined. Reddy et al. [73]
analyzed MHD boundary layer flow along vertical plates with ramped temperature,
while Reddy et al. [88] considered MHD free convective mass transfer with variable
suction and Soret effect.
Krishna et al. [98] examined the Hall effects on MHD peristaltic flow of Jeffrey
fluid through a porous medium in a vertical stratum, shedding light on the interactions
within complex fluid systems under the influence of a magnetic field and porous
constraints. Krishna and Chamkha [99] studied the Hall effects on MHD squeezing flow
of a water-based nanofluid between two parallel disks. This research focuses on how
magnetic fields and Hall effects can manipulate nanofluid flow dynamics, which is crucial
for applications in cooling systems and devices. Also Krishna et al. [100] explored heat
and mass transfer in the free convective flow of a micropolar fluid through a porous
surface with an inclined magnetic field and Hall effects. Their findings highlight the
complex interactions between thermal and velocity fields in micropolar fluids. At the
same year Krishna et al. [101] investigated the influence of Soret, Joule, and Hall effects
on MHD rotating mixed convective flow past an infinite vertical porous plate. This
study is pivotal for understanding how multiple forces interact in an MHD environment.
Krishna and Chamkha [102] expanded their research to explore the Hall and ion slip
effects on unsteady MHD convective rotating flow of nanofluids, with an emphasis on
biomedical engineering applications. This highlights the potential for using MHD flows
5.1 Introduction 95
in medical devices and treatments. Finally, in the same year, Krishna et al. [77] delved
into the heat and mass transfer on MHD flow of a second-grade fluid through a porous
medium over a semi-infinite vertical stretching sheet. This work provides insights into
the behavior of more complex fluid models under the influence of MHD and stretching
effects.
Chemical reaction and thermal radiation were combined in studies by Muthu-
cumaraswamy and Sivakumar [21] and El-Fayez [177], considering chemical reaction
effects on vertical plates with radiation and thermal stratification. Kandasamy et al.
[37] analyzed the effects of chemical reaction, heat, and mass transfer on MHD flow over
a vertical surface, considering heat source and thermal stratification. They focus on
the impacts of chemical reaction on fluid motion near oscillating plates and stretching
surfaces, under various conditions like heat flux, mass transfer, and magnetic fields,
providing comprehensive insights into the complex interplay of thermal stratification,
heat, and mass transfer in fluid dynamics.
Recent research has delved into the effects of chemical reaction in these flows. Nath
and Deka [82] and Nath et al. [83] studied unsteady flow past vertical plates with
chemical reaction in the presence of thermal stratification. Nath and Deka [178–180]
collectively explored the complex dynamics of magnetohydrodynamic (MHD) nanofluids
in different porous media environments, with an emphasis on the combined effects
of thermal and mass stratification. Lastly, non-Newtonian fluids were considered in
porous media with chemical reaction and thermal stratification by Megahed and Abbas
[181], showcasing the breadth and depth of research in this field.
More recent studies have extensively analyzed MHD flows, focusing on the combined
effects of magnetic fields, chemical reaction, and thermal radiation across various
conditions. Notably, research by Yanala et al. [182] investigate the 3-D radiative
MHD flow over a stretching sheet, highlighting the effects of chemical reaction and
the stretching geometry on Casson fluids. Reddy et al. [183] and Shankar Goud et al.
[184] explored similar themes but focused on the impact of radiation and porosity on
two-dimensional MHD flows at stagnation points with viscous dissipation, and also
considerd the heat transfer characteristics of nanofluids on an inclined spinning disk and
over nonlinearly porous sheets. Kumar and Reddy [84] studied the influence of thermal
radiation and nanoparticles on boundary layer flows, revealing the nuanced interactions
within Maxwell fluids. Additionally, studies by Reddy and Mangamma [185] Yanala
et al. [186] Reddy and Mangamma, [187] and Goud and Reddy [188] delved into
specific scenarios such as nanofluid flow over a nonlinearly porous stretching sheet
with nonuniform heat source, flow past a non-linear stretching Riga plate, flow past
a porous shrinking sheet, and flow over a permeable wedge, all considering elements
like heat generation/absorption and radiation. The research is further enriched by
computational analyses, such as those by Reddy et al. [189] who examined the
independence of computational meshes in their models. The study by Suhasini et
al. [190] delves into the dynamics of magneto-micropolar fluid flow over a porous
stretching surface, addressing the complex interactions influenced by a non-uniform
heat source/sink and the presence of Soret-Dufour effects. Their research focuses on
96 Influences of Thermal Stratification and Chemical Reaction on MHD...
the nuanced behavior of micropolar fluids, which are fluids with microstructure and
are capable of exhibiting microrotational effects, under the influence of magnetic fields.
The following points highlight the novelties of the current study: we achieve an
exact solution through the Laplace transform approach without any error, proving
it to be a highly efficient method for obtaining precise solutions. The concurrent
influence of thermal stratification and chemical reaction on the flow past an accelerated
vertical plate has not been previously investigated, making this study the pioneering
endeavor in this field. The goal of this research was to analytically explore MHD
unsteady flow across an accelerated vertical plate, considering the effects of thermal
stratification and chemical reaction. Subsequently, results for the stratified fluid are
contrasted with those from the classical case without stratification. The effects of
physical parameters including γ, Gr, Gc, Sc, Kc , Magnetic parameter M , Chemical
reaction Kc , and Darcy number Da on the profiles, are examined and presented in
graphical form. These insights deepen our understanding of natural convection under
diverse fluid and environmental conditions, offering profound implications for industries
like chemical manufacturing, polymer production, and crystal separation.
• It is assumed that both the diffusing species and the fluid engage in a first-order
chemical reaction, influencing the flow dynamics.
• The physical setup of the problem is depicted in figure 5.1, where the x′ axis is
aligned with the flow along the plate and the y ′ axis is perpendicular to it.
• At a time t′ > 0, the temperature of the plate is increased following the formula
2
′ u0 ′
′
T∞ + (Tw′ − T∞ ) ν t , and similarly, the concentration is adjusted to C∞
′
+ (Cw′ −
u2
0 ′
′
C∞ )e ν t .
referred to as the work of compression, causes reversible work to be done on the fluid
particles. For the purpose of evaluating numerical models, compression work is added
to thermal stratification as an additive component.
u′ = 0, T ′ = T∞
′
, C ′ = C∞
′
,∀y ′ , t′ ≤ 0
2 u2
′ u0 ′ 0 ′
t′ > 0 : u′ = u0 t′ a′ , T ′ = T∞
′
+ (Tw′ − T∞ ) t , C ′ = C∞
′
+ (Cw′ − C∞
′
)e ν t aty ′ = 0
ν
u = 0, T → T∞ , C → C∞ , as y ′ → ∞
′ ′ ′ ′ ′
Together, these boundary and initial conditions are crucial for modeling the prob-
lem accurately. They provide necessary constraints that help define the fluid flow,
temperature, and concentration fields under the influence of plate movement, thermal
98 Influences of Thermal Stratification and Chemical Reaction on MHD...
stratification, and chemical reaction, ensuring that the physical behavior of the system
is realistically and comprehensively captured.
The subsequent dimensionless variables are introduced:
u0 y ′ u2 t′ T ′ − T∞
′
u′ a′ ν ν
y= , t= 0 , T = ′ , u = , a = 2 , Pr =
ν ν Tw − T∞ ′ u0 u0 α
′ ′ ′ 2 ′
Kν C − C∞ νq0 uK ν
K= 2 , C= ′ , ϕ= 2
, Da = 0 2 p , Sc =
u0 Cw − C∞ ′ ρCp u0 ν D
2
σB0 ν γν′
gβν(Tw − T∞ )
′ ′ ′
gβ ′ ν(Cw′ − C∞ )
M= , γ= , Gr = , Gc =
ρu0 2
u0 (Tw − T∞ )
′ ′ u03
u03
∂u ∂ 2u
= 2 + GrT + GcC − N u (5.4)
∂t ∂y
∂T 1 ∂ 2T
= − ϕT − γu (5.5)
∂t P r ∂y 2
∂C 1 ∂ 2C
= − KC (5.6)
∂t Sc ∂y 2
Where,
1
N =M+
Da
The boundary and initial conditions are expressed in dimensionless terms as follows:
u = 0, T = 0, C=0 ∀y and t ≤ 0
u = at, T = t, C = et at y = 0, t > 0 (5.7)
u → 0, T → 0, C→0 as y → ∞, t > 0
C = η1 (K) (5.10)
Where,
s
1 1
2
y
η= √ , A−B = M+ −ϕ − 4γGr, A+B =M + + ϕ,
2 t Da Da
v
1 2 Sc K − N − ϕ u (N − ϕ)2 − 4γGr
u
N =M+ , Z1 = , Z2 = t ,
Da 2(Sc − 1) 4(Sc − 1)2
Sc K − ϕ Sc K − B Sc K − N
Z3 = , Z4 = , Z5 = ,
Sc − 1 Sc − 1 Sc − 1
1 + Z3 Z3 − Z1 − Z2 Z2 − Z1 + Z3
E1 = , E2 = , E3 = ,
(1 + Z1 )2 − Z22 2Z2 (1 + Z1 + Z2 ) 2Z2 (Z2 − Z1 − 1)
(1 + Z3 ) (1 + Z4 ) (Z3 − Z1 − Z2 )(Z4 − Z1 − Z2 )
E4 = , E5 = ,
(1 + Z1 )2 − Z22 2Z2 (1 + Z1 + Z2 )
(Z2 − Z1 + Z3 )(Z2 − Z1 + Z4 )
E6 = ,
2Z2 (Z2 − Z1 − 1)
Gc(E4 − 1) Gc E5 Gc E6 a(N − B) − Gr
ϵ1 = , ϵ2 = , ϵ3 = , ϵ4 = ,
A−B A−B A−B A−B
100 Influences of Thermal Stratification and Chemical Reaction on MHD...
Gc E1 Gc E2 Gc E3 a(A − N ) + Gr
ϵ5 = , ϵ6 = , ϵ7 = , ϵ8 = ,
Sc − 1 Sc − 1 Sc − 1 A−B
γ Gc (E4 − 1) γ Gc E5 γ Gc E6
λ1 = , λ2 = , λ3 = ,
(A − B)(A − ϕ) (A − B)(A − ϕ) (A − B)(A − ϕ)
γ (a(N − B) − Gr)γ γ Gc E1 γ Gc E2
λ4 = , λ5 = , λ6 = ,
(A − B)(A − ϕ) (Sc − 1)(B − ϕ) (Sc − 1)(B − ϕ)
γ Gc E3 γ Gc (E4 − 1) γ Gc E5
λ7 = , λ8 = , λ9 = ,
(Sc − 1)(B − ϕ) (A − B)(A − ϕ) (A − B)(A − ϕ)
γ Gc E6 γ (a(A − N ) + Gr)γ γ Gc E1
λ10 = , λ11 = , λ12 = ,
(A − B)(A − ϕ) (A − B)(B − ϕ) (1 + Z3 ) (Sc − 1)2
γ Gc E2 γ Gc E3
λ13 = , λ14 = ,
(Sc − 1)2 (Z1 − Z2 − Z3 ) (Sc − 1)2 (Z1 − Z2 − Z3 )
γ Gc E1 γ Gc E2
λ15 = , λ16 = ,
(Sc − 1)2 (1 + Z3 ) (Sc − 1)2 (Z3 − Z2 − Z1 )
γ Gc E3
λ17 =
(Sc − 1)2 (Z3 − Z2 − Z1 )
Additionally, the inverse Laplace transforms, defined as follows:
√ √
e−y s+A e−y s+A
( ) ( )
η1 (A) = L −1
, η2 (A, Z1 , Z2 ) = L −1
,
s−1 s + Z1 + Z2
√ √
e−y s+A e−y s+A
( ) ( )
η3 (A) = L −1
, η4 (A, Z5 ) = L −1
s2 s + Z5
!
Gr Gc
uc = a− η3 (N ) + [η1 (N ) − η4 (N, Z5 )]
N −ϕ (Sc − 1)(Z5 + 1)
Gr Gc
+ η3 (ϕ) − [η1 (K) − η4 (K, Z5 )] (5.11)
N −ϕ (Sc − 1)(Z5 + 1)
Tc = η3 (ϕ) (5.12)
5.3 Method of Solution 101
5.3.2 Skin-Friction
The dimensionless Skin-Friction, defined as the surface shear stress, is calculated by:
du
τ =−
dy
y=0
Now from equation (5.8), we get the expression for skin friction as:
√ e−At
q !
τ = ϵ1 e A + 1 erf ( (A + 1)t) + √
t
πt
e−At
q q !
+ϵ2 e(−Z1 −Z2 )t A − Z1 − Z2 erf ( (A − Z1 − Z2 )t) + √
πt
−At
!#
q q e
+ϵ3 e (−Z1 +Z2 )t
A − Z1 + Z2 erf ( (A − Z1 + Z2 )t) + √
πt
√
√ √
s
t erf ( At)
+ϵ4 t A erf ( At) + e−At + √
π 2 A
√ e−Bt
q !
+(ϵ5 − ϵ1 ) e B + 1 erf ( (B + 1)t) + √
t
πt
e−Bt
q q !
+(ϵ6 − ϵ2 ) e (−Z1 −Z2 )t
B − Z1 − Z2 erf ( (B − Z1 − Z2 )t) + √
πt
e−Bt
q q !
+(ϵ7 − ϵ3 ) e (−Z1 +Z2 )t
B − Z1 + Z2 erf ( (B − Z1 + Z2 )t) + √
πt
√
√ √
s
t erf ( Bt)
+ϵ8 t B erf ( Bt) + e−Bt + √
π 2 B
s
q q Sc −Kt
−ϵ5 et Sc(K + 1) erf ( (K + 1)t) + e
πt
s !
q q Sc −Kt
−ϵ6 (−Z1 −Z2 )t
e Sc (K − Z1 − Z2 ) erf ( (K − Z1 − Z2 )t) + e
πt
s !#
q q Sc −Kt
−ϵ7 (−Z1 +Z2 )t
e Sc (K − Z1 + Z2 ) erf ( (K − Z1 + Z2 )t) + e
πt
For the classical case, the Skin-Friction solution is derived from the expression
(5.11), which is represented by
102 Influences of Thermal Stratification and Chemical Reaction on MHD...
√
Gr √ √
s
t −N t erf N t
τc = a+ t N erf N t + e + √
ϕ−N π 2 N
Gc √ q e−N t q
+ et N + 1 erf ( (N + 1)t) + √ − e−Z5 t N − Z5
(Sc − 1)(Z5 + 1) πt
s √
e−N t (
!
q Gr q q t erf ϕt)
erf ( (N − Z5 )t) − √ + t ϕ erf ( ϕt) + e−ϕt + √
πt N −ϕ π 2 ϕ
√
Gc q q e−Kt Sc
− et (K + 1)Sc erf ( (K + 1)t) + √
(Sc − 1)(Z5 + 1) πt
!
q q Sc
−e−Z5 t Sc(K − Z5 ) erf ( (K − Z5 )t) − √ e−Kt
πt
dT
Nu = −
dy
y=0
Now from equation (5.9), we get the expression for the Nusselt number as:
5.3 Method of Solution 103
s √
q q t −ϕt erf ϕt
N u = (1 − λ4 − λ11 )t ϕerf ϕt + e + √
π 2 ϕ
e−ϕt
q q !
+(−λ1 − λ5 + λ8 + λ15 ) e ϕ + 1 erf ( (ϕ + 1)t) + √
t
+ (−λ2 − λ6
πt
e−ϕt
q q !
+λ9 + λ16 ) e (−Z1 −Z2 )t
ϕ − Z1 − Z2 erf ( (ϕ − Z1 − Z2 )t) + √ + (−λ3
πt
e−ϕt
q q !
−λ7 + λ10 − λ17 ) e (−Z1 +Z2 )t
ϕ − Z1 + Z2 erf ( (ϕ − Z1 + Z2 )t) + √
πt
−ϕt
!
q q e
+(−λ15 − λ16 + λ17 ) e−Z3 t ϕ − Z3 erf ( (ϕ − Z3 )t) + √
πt
√ −At
!
q e
+λ1 et A + 1 erf ( (A + 1)t) + √
πt
e−At
q q !
+λ2 e (−Z1 −Z2 )t
A − Z1 − Z2 erf ( (A − Z1 − Z2 )t) + √
πt
−At
!#
q q e
+λ3 e (−Z1 +Z2 )t
A − Z1 + Z2 erf ( (A − Z1 + Z2 )t) + √
πt
√
√ √
s
t erf ( At)
+λ4 t A erf ( At) + e−At + √
π 2 A
√ e−Bt
q !
+(λ5 − λ8 ) e B + 1 erf ( (B + 1)t) + √
t
πt
e−Bt
q q !
+(λ6 − λ9 ) e (−Z1 −Z2 )t
B − Z1 − Z2 erf ( (B − Z1 − Z2 )t) + √
πt
−Bt
!
q q e
+(λ7 − λ10 ) e (−Z1 +Z2 )t
B − Z1 + Z2 erf ( (B − Z1 + Z2 )t) + √
πt
√
√ √
s
t erf ( Bt)
+λ11 t B erf ( Bt) + e−Bt + √
π 2 B
s
q q Sc −Kt
+λ12 et Sc(K + 1) erf ( (K + 1)t) + e
πt
s !
q q Sc −Kt
+λ13 (−Z1 −Z2 )t
e Sc (K − Z1 − Z2 ) erf ( (K − Z1 − Z2 )t) + e
πt
s !#
q q Sc −Kt
+λ14 (−Z1 +Z2 )t
e Sc (K − Z1 + Z2 ) erf ( (K − Z1 + Z2 )t) + e
πt
104 Influences of Thermal Stratification and Chemical Reaction on MHD...
For the classical case, the Nusselt number solution is derived from the expression
(5.12), which is represented by
s √
q q t −ϕt erf ϕt
N uc = t ϕ erf ϕt + e + √ (5.13)
π 2 ϕ
dC
Sh = −
dy
y=0
The calculation of the Sherwood number is derived from the Concentration profile
C, denoted by (5.10), in the following manner:
√
q q Sc e−Kt
Sh = e t
Sc(1 + K) erf (1 + K)t + √ (5.14)
πt
The velocity profiles for various values of γ, Gr, Gc, M, Da, Kc and (Sc) are illus-
trated in Figures 5.4 through 5.9, both with and without stratification. In Figure
5.4, observations indicate that thermal stratification reduces the flow velocity. As the
thermal stratification parameter (γ) increases, there is a decline in the efficiency of
convective potential between the hot plate and the adjacent fluid. This decline in the
buoyancy force, in turn, leads to a reduction in the velocity of the flow. Figures 5.5 to
5.8 illustrate that higher values of Gr and Gc correlate with an increase in fluid velocity,
indicating enhanced buoyancy forces that drive the fluid motion. Similarly, increased
values of Sc and Kc lead to a velocity decrease, suggesting more effective momentum
diffusion. Conversely, higher M values reduce the velocity due to the damping effect
of the magnetic field, which opposes the fluid’s motion by inducing a Lorentz force
that acts as a brake. In Figure 5.9, the velocity profile increases along with Da as the
porosity of the fluid increases, allowing the fluid particles to flow more freely.
Figure 5.10 illustrates the influence of thermal stratification (γ) on temperature
distribution. A decrease in the temperature differential between the plate’s surface
and the fluid around it is the result of thermal stratification (γ). The result is a
decrease in temperature due to the thickening of the thermal boundary layer. It has
been shown that temperature falls as there is an increase in the thermal stratification
parameter (γ). Figures 5.11 to 5.13 show that higher Kc and Sc values correlate
with increased temperature, possibly due to enhanced reaction rates and reduced
diffusive transport of heat. Conversely, a rise in ϕ increases temperature, likely due to
endothermic reactions absorbing heat. As Gr and Gc increasees, temperature decreases,
suggesting that enhanced buoyancy-driven fluid movement enhances heat dissipation
110 Influences of Thermal Stratification and Chemical Reaction on MHD...
away from the heated surface. As observed in Figures 5.14 and 5.15, an increase in the
Magnetic parameter (M) corresponds to an elevation in temperature, while an upsurge
in the Darcy number (Da) results in a decrease in temperature, showcasing the inverse
relationship between Da and thermal behavior.
Figures 5.16 and 5.18 provide a detailed depiction of how thermal stratification
influences the velocity and temperature of a fluid as time progresses. In scenarios
where stratification is not present, there is a noticeable and steep increase in both the
fluid’s temperature and velocity over time. However, this trend undergoes a significant
transformation with the introduction of thermal stratification. Once stratification is
applied, both temperature and velocity begin to exhibit a pattern of gradual stabiliza-
tion, diverging from the rapid increases observed in the absence of stratification. This
process of stabilization results in lower temperatures and velocities when compared
to the standard scenario where γ = 0. This nuanced understanding underscores the
critical role of thermal stratification in accurately modeling and interpreting fluid
dynamics over time. Furthermore, as depicted in Figures 5.17 and 5.19, an increase
in the thermal Grashof number (Gr) and the mass Grashof number (Gc) correlates
with an increased fluid velocity due to stronger buoyancy forces. However, as these
numbers rise, temperature decreases, indicating that the enhanced fluid motion effec-
tively disperses heat, leading to more efficient cooling and a lower temperature near
the heated surface.
Figure 5.20 shows how the fluid’s concentration is affected by the parameters Sc
and Kc . An observation is made that a rise in the values of Sc and Kc corresponds
to a decrease in fluid concentration. This occurs because higher Sc value slows mass
5.4 Results and Discussion 113
diffusion, and elevated Kc enhances chemical reaction that consume the diffusing
species, thereby reducing the concentration. Figures 5.21 and 5.22 depict how thermal
stratification affects variations in skin friction and nusselt number. Over time, a trend is
observed where skin friction diminishes while the Nusselt number experiences an upward
trajectory. This pattern starkly contrasts with the distinct scenario observed in the
absence of stratification. Remarkably, the introduction of thermal stratification, which
leads to an enhancement of both skin friction and the Nusselt number, indicating a
significant shift in thermal and frictional properties under the influence of stratification.
In Figure 5.23 an increase in the Schmidt number Sc and the chemical reaction rate
constant Kc corresponds to a higher Sherwood number, signifying an improvement in
the mass transfer rate during the chemical process.
5.5 Conclusion
The prominent outcomes of the present work are as follows:
• The movement of the fluid is hindered by the Lorentz force, which originates
electromagnetically, while it is aided by the permeability of the porous medium.
• The buoyancy force acts to accelerate the fluid motion.
• In this study, results are juxtaposed with those from the classical scenario, where
there is no stratification.
• As γ, M Kc and Sc values increase, fluid velocity decreases, whereas it rises with
an increase in Gr, Gc and Da.
• This study’s incorporation of thermal stratification, leading to lower velocity and
temperature compared to the classical case (γ = 0), enhances its relevance over
previous studies.
• Temperature decreases with an increase in γ, Gr, Gc, and Da, whereas it rises
with an increase in ϕ, M, Kc and Sc.
• The presence of stratification elevates both skin friction and the Nusselt number
compared to scenarios lacking stratification.
• It’s noteworthy that the destructive reaction does not contribute to an increase
in the distribution of concentration.
• The changes in rates at the boundary surface are crucial in governing the overall
distribution of flow, heat, and mass.
Numerical study of MHD Ternary
Nanofluid
6
(Ag − CoF e2O4 − ZnO/C2H6O2 + H2O):
Effects of Thermal Stratification on
a Vertically Stretching Cylinder in a
Porous Environment
Published in Biointerface Research in Applied Chemistry, 14, 6 (2024),
ISSN: 2069-5837. (SCOPUS)
6.1 Introduction
This study explores the numerical investigation of an MHD ternary nanofluid con-
sisting of silver (Ag), cobalt ferrite (CoF e2 O4 ), and zinc oxide (ZnO) nanoparticles,
uniformly dispersed in a mixture of ethylene glycol (C2 H6 O2 ) and water (H2 O). This
distinctive ternary nanofluid offers remarkable properties, making it suitable for diverse
applications. Incorporating silver (Ag) nanoparticles enhances thermal conductivity,
while cobalt ferrite (CoF e2 O4 ) and zinc oxide (ZnO) nanoparticles improve heat
transfer performance and stability. Such a nanofluid can effectively serve in heat
exchangers, cooling systems, and electronic devices, optimizing thermal management
and enhancing heat dissipation. Furthermore, silver nanoparticles exhibit antibacterial
traits, while ZnO nanoparticles provide photocatalytic action against bacteria. This
120 Numerical study of MHD Ternary Nanofluid ...
ternary nanofluid formulation (Ag − CoF e2 O4 − ZnO) is also valuable for creating
antimicrobial coatings on various surfaces, including medical devices, textiles, and
food storage containers. These coatings effectively restrict bacterial growth, ensur-
ing enhanced hygiene. Zinc oxide (ZnO) nanoparticles further contribute to organic
pollutant degradation and water purification through photocatalytic properties. The
integration of Ag − CoF e2 O4 − ZnO nanoparticles demonstrates potential in improving
thermal conductivity, heat transfer efficiency, antibacterial properties, and photocat-
alytic activity, making it ideal for cooling systems, antimicrobial coatings, and water
treatment solutions to improve water quality and pollutant removal.
The concept of nanofluid, initially introduced by scientists, transformed fluid
dynamics through the incorporation of nanoparticles. The concept of nanofluid, first
introduced by Choi and Eastman,[6], suggests that heat transfer fluids containing
metallic nanoparticles could revolutionize heat transfer fluids. In groundbreaking
studies, Das and Jana [60] and Das et al. [61] explored the natural convective flow of
nanofluids with radiation in a moving vertical plate and a vertical channel, respectively.
Both investigations examined water-based nanofluids composed of titanium dioxide,
aluminum oxide, and copper. Addressing heat production or absorption scenarios,
Rashidi et al. [192] sought a Lie group solution for the flow of nanofluid over a
chemically reacting horizontal plate. The entropy analysis by Abolbashari et al. [62]
employed the Homotopy Analysis Method (HAM) to examine a nanofluid composed
of water and one of four distinct types of nanoparticles: T iO2 , Al2 O3 , Cu, and CuO,
passing across a stretchable permeable surface. In a numerical study, Motsumi and
Makinde [193] focused on the boundary layer movement of nanofluids over a moving
flat plate, specifically the impacts of thermal radiation, viscous dissipation, and thermal
diffusion. Considering velocity slip and temperature jump conditions, Turkyilmazoglu
[194] conducted an analytical study on MHD nanofluid flow using various water-based
nanoparticles across a continuously stretching/shrinking permeable sheet. Further,
Jeevandhar et al. [24] explored the MHD boundary layer nanofluid flow under thermal
radiation over an exponentially accelerating vertical plate. In addition, Cheng [39, 40]
investigated double-diffusive natural convection and combined heat and mass transfer
in non-Newtonian fluid-saturated porous media with thermal and mass stratification.
Recent studies by Sahu and Deka [123, 149, 146], Das and Deka [195], Nath et al. [180],
Nath and Deka [196, 178] focused on thermal and mass stratification effects on MHD
flow with variable temperature and mass diffusion also including chemical reaction, heat
sources, and thermal radiation in their theoretical analyses. The impact of nanoparticle
shape on heat transfer enhancement in magnetized ternary nanofluid was explored by
Shanmugapriya et al. [197]. Arif et al. [198] investigated heat transfer in a radiator
with differently shaped nanoparticles in a water-based ternary nanofluid, utilizing a
fractional model. Finally, Patil et al. [199] explored the flow of tangent hyperbolic
ternary nanofluid over a rough-yawed cylinder due to impulsive motion. Mahmood et al.
[200] examined the effects of suction and a heat source on MHD stagnation point flow
of ternary nanofluid over a convectively heated cylinder. Additionally, Cao et al. [201]
simulated the dynamics of water-based colloidal mixtures with various nanoparticles,
6.1 Introduction 121
∂(ru) ∂(rv)
+ = 0 (6.1)
∂x ∂r
µmnf 1 ∂ (ρβT )mnf
!
∂u ∂u ∂u σmnf 2
u +v = r + g(T − T∞ ) − B u
∂x ∂r ρmnf r ∂r ∂r ρmnf ρmnf 0
µmnf u
− (6.2)
ρmnf k ′
kmnf 1 ∂
!
∂T ∂T ∂T Q0
u +v = r + (T − T∞ ) (6.3)
∂x ∂r (ρcp )mnf r ∂r ∂r (ρcp )mnf
r2 − r02 T − T∞ (x)
s
a aνf
r
η= , ψ= xr0 f (η), θ=
2r0 νf l l Tw (x) − T0
and the non-dimensional quantities are provided as follows:
νf (µcp )f
s
lσf B02 lνf lνf Gr B Q0 l
M= , K= , γ= , λ= , δ= , Q= , Pr =
aρf ak ar02 Re2x A a(ρcp )f kf
6.2 Mathematical Analysis 123
where,
x4 x6
a1 = x 1 x 2 , a2 = x2 x3 , a3 = x2 x5 , a4 =
, a5 = Qx4
Pr
Here, the variables µmnf , (βT )mnf , ρmnf , (ρCp )mnf , kmnf , σmnf refer to the viscosity,
thermal expansion coefficient, electrical conductivity, heat capacity, thermal con-
ductivity, and density of the ternary nanofluid, respectively. In a similar manner,
µf , (βT )f , ρf , (ρCp )f , kf , σf correspond to the viscosity, thermal expansion coefficient,
electrical conductivity, heat capacity, thermal conductivity, and density of the base
fluid, respectively.
Here is how the transformed boundary conditions are defined:
(σs2 + 2σnf ) − 2ϕ2 (σnf − σs2 ) (σs3 + 2σhnf ) − 2ϕ3 (σhnf − σs3 )
" # " #
σhnf = σnf , σmnf = σhnf
(σs2 + 2σnf ) + ϕ2 (σnf − σs2 ) (σs3 + 2σhnf ) + ϕ3 (σhnf − σs3 )
Here, ϕ1 , ϕ2 , and ϕ3 signify the volume fractions of Ag (Silver), CoF e2 O4 (Cobalt
ferrite), and ZnO (Zinc Oxide) nanoparticles, respectively. The designations mnf, hnf,
nf, f, s1, s2, and s3 refer to ternary nanofluid, hybrid nanofluid, nanofluid, base fluid,
and the dense nanoparticles of Silver (Ag), Cobalt ferrite (CoF e2 O4 ), and Zinc Oxide
(ZnO), respectively.
The definitions of the skin friction coefficient and local Nusselt number are as
follows:
1 kmnf ′
x =
Cf Re1/2 f ′′ (0) and N ux Re−1/2 =− θ (0)
(1 − ϕ1 )2.5 (1 − ϕ2 ) (1 − ϕ3 )
2.5 2.5 x
kf
Table 6.1 Thermal and Physical Properties of Ethylene glycol+ water and nanoparticles
(Paul et al. [1, 2])
y(1)
y(2)
y(2)
y(3)
d
y(3) =
y(2)2 −y(1)y(3)−a2 λy(4)+(a3 M +a1 K)y(2)−2a1 γy(3)
(1+2γη)a1
dη
y(4) y(5)
y(5) y(2)(y(4)+δ)−y(1)y(5)−2a4 γy(5)−a5 y(4)
(1+2γη)a 4
transfer through the fluid. Consequently, the rate of heat transfer diminishes, the
thermal boundary layer narrows, and the fluid temperature decreases.
Figure 6.11 shows that a rise in the percentage of solid volume of Ag nanoparticles,
with the volume fractions of CoF e2 O4 and ZnO held constant, causes an increase in
the temperature profile. Furthermore, Figure 6.12 demonstrates that the temperature
also increases as the percentage of volume that CoF e2 O4 rises, while keeping the
128 Numerical study of MHD Ternary Nanofluid ...
Ag and ZnO fractions constant. Similarly, an elevated volume fraction of ZnO (ϕ3 ),
as depicted in Figure 6.13, results in a higher temperature profile, indicating that
increased concentrations of these nanoparticles improve the fluid’s thermal properties.
Tables 6.3 detail the effects of various flow parameters on skin friction and the
local Nusselt number. With a rise in the thermal stratification parameter (δ), both
skin friction and the local Nusselt number decrease. Higher values of the M and the
6.4 Results and Discussion 129
Table 6.3 Skin friction Coefficient and Local Nusselt number for different values of
δ, M, K, λ, P r, Q, γ
K also lead to a decrease in both skin friction and the Nusselt number. In contrast,
an increase in the value of λ results in higher skin friction and Nusselt number. For
the P r, while skin friction decreases, the Nusselt number increases. Similarly, a rise
in the Q decreases the heat transfer rate, but skin friction increases. Lastly, as the γ
increases, skin friction decreases, but the Nusselt number rises.
6.4 Results and Discussion 131
Tables 6.4 and 6.5 present comparisons of the absolute skin friction (−Cf Re1/2
x ) and
the local Nusselt number (N ux Rex ) across nanofluid, hybrid nanofluid, and ternary
−1/2
for −Cf Re1/2x . The Nusselt number for the ternary nanofluid shows a significant
improvement over that of the conventional nanofluid. Moreover, an increase in thermal
stratification (δ) augments the difference in the Nusselt number between the ternary
nanofluid and the regular nanofluid for N ux Re−1/2
x .
6.5 Conclusion 133
6.5 Conclusion
An extensive study was conducted on the consequences of the thermal stratification
effect regarding the movement of a ternary nanofluid that demonstrates magnetohydro-
dynamics characteristics through a vertically stretchable cylinder. This study took into
account the existence of a thermal buoyancy effect and a heat source within a porous
medium. The study also looks at the flow features and how they affect temperature
and velocity profiles, as well as the Nusselt number and friction factor. Notable findings
from this inquiry include the following: As the parameters δ, K, and M grow, the
velocity profile shows a decreasing trend, while it shows a rising pattern with higher
values of λ, and γ. With rising values of δ and P r, the temperature declines, whereas
with rising values of ϕ1 , ϕ2 , ϕ3 , P r, γ and Q, it rises. Parameters δ, M, K, Q, γ, and P r
cause a decrease in absolute skin friction, whereas λ cause it to increase. An upward
trend is noted in the local Nusselt number as the values of P r and γ increase, whereas
a downward trend is observed as δ, M, K, λ, and Q increase. In comparison to the
hybrid nanofluid, the ternary nanofluid exhibits higher absolute skin friction. Similarly,
when compared to ordinary nanofluids, hybrid nanofluids exhibit higher skin friction.
A ternary nanofluid has a greater heat transfer rate than a hybrid nanofluid, which is
also higher than that of standard nanofluids.
MHD Natural Convection of
Ternary Nanofluid in a Porous
7
Medium Incorporating
Stratification and Radiation
Published in ZAMM - Journal of Applied Mathematics and Mechanics,
ISSN: 0044-2267, E-ISSN: 1521-4001, (SCOPUS, SCIE)
7.1 Introduction
Ternary nanofluids are specialized fluids comprising three distinct types of nanoparticles
dispersed within a base fluid. This paper presents a numerical investigation into an
MHD ternary nanofluid comprised of silver (Ag), cobalt ferrite (CoF e2 O4 ), and zinc
oxide (ZnO) nanoparticles, which are evenly distributed in a blend of ethylene glycol
(C2 H6 O2 ) and water (H2 O). This ternary nanofluid possesses unique properties that
render it appropriate for various uses. The inclusion of silver (Ag) nanoparticles
in the nanofluid has been demonstrated to improve thermal conductivity, while the
inclusion of cobalt ferrite (CoF e2 O4 ) and zinc oxide (ZnO) nanoparticles enhances
heat transfer effectiveness and stability. This nanofluid can be used in multiple areas
including heat exchangers, cooling systems, and electronic devices, aiming to enhance
heat dissipation and improve thermal management. Additionally, silver nanoparticles
exhibit antibacterial properties, while ZnO nanoparticles possess photocatalytic action
against bacteria. The use of the ternary nanofluid, comprising Ag − CoF e2 O4 − ZnO,
136 MHD Natural Convection of Ternary Nanofluid in a Porous Medium...
∂(ru) ∂(rv)
+ = 0 (7.1)
∂x ∂r
µmnf 1 ∂ (ρβT )mnf
!
∂u ∂u ∂u σmnf 2
u +v = r + g(T − T∞ ) − B u
∂x ∂r ρmnf r ∂r ∂r ρmnf ρmnf 0
µmnf u
− (7.2)
ρmnf k ′
kmnf 1 ∂ 1 1 ∂ ∂qr′
! !
∂T ∂T ∂T
u +v = r − r
∂x ∂r (ρcp )mnf r ∂r ∂r (ρcp )mnf r ∂r ∂r
Q0
+ (T − T∞ ) (7.3)
(ρcp )mnf
7.2 Mathematical Analysis 139
The radiative heat flux qr′ in Eqn. (7.3) is estimated using the Rosseland approx-
imation Rosseland [219], under the assumption that the nanofluid is optically thick.
According to this approximation, the radiative heat flux is expressed as:
4σ ∗ ∂T 4
qr′ = − ,
3k ∗ ∂r
where σ ∗ is the Stefan-Boltzmann constant and k ∗ is the mean absorption coefficient.
Under the assumption that temperature differences within the fluid are relatively
small, T 4 is linearized using a Taylor series expansion around T∞ , yielding T 4 ≈
4T∞3
T − 3T∞4
. Substituting this into the radiative heat flux term and differentiating,
the energy equation is modified to account for radiative contributions. This method
ensures the inclusion of radiative heat transfer effects in the thermodynamic model
without introducing excessive computational complexity. Following the studies by
Turkyilmazoglu [194] and Ghadikolaei, Gholinia [220], Eqn. (7.3) is reformulated as
follows:
16σ ∗ T∞
′3
1 ∂
! !
∂T ∂T kf kmnf ∂T Q0
u +v = + r + (T − T∞ ) (7.5)
∂x ∂r (ρcp )mnf kf 3kf k r ∂r ∂r (ρcp )mnf
r2 − r02 T − T∞ (x)
s
a aνf
r
η= , ψ= xr0 f (η), θ= (7.6)
2r0 νf l l Tw (x) − T0
s
lσf B02 lνf lνf Gr B Q0 l
M= , K= , γ= , λ= , δ= , Q=
aρf ak ar02 Re2x A a(ρcp )f
16σ ∗ T∞
′3
νf (µcp )f
R= , Pr =
3kf k kf
where δ represents thermal stratification, γ curvature, M magnetism, λ thermal
buoyancy, K porosity, Q heat source/sink, R radiation and P r the Prandtl number.
140 MHD Natural Convection of Ternary Nanofluid in a Porous Medium...
To meet the requirements of the continuity equation (7.1), the stream function ψ is
introduced, with u = 1r ∂ψ
∂r
and v = − 1r ∂ψ
∂x
. Consequently, the dimensionless versions of
the transformed equations are displayed as follows:
where,
x4 x6 + R
a1 = x 1 x 2 , a2 = x2 x3 , a3 = x2 x5 , a4 =
, a5 = Qx4
Pr
Here, the variables µmnf , (βT )mnf , ρmnf , (ρCp )mnf , kmnf , σmnf refer to the viscosity,
thermal expansion coefficient, electrical conductivity, heat capacity, thermal con-
ductivity, and density of the ternary nanofluid, respectively. In a similar manner,
µf , (βT )f , ρf , (ρCp )f , kf , σf correspond to the viscosity, thermal expansion coefficient,
electrical conductivity, heat capacity, thermal conductivity, and density of the base
fluid, respectively.
Here is how the transformed boundary conditions are defined:
Here, ϕ1 , ϕ2 , and ϕ3 represent the volume fractions of Silver (Ag), Cobalt Ferrite
(CoF e2 O4 ), and Zinc Oxide (ZnO) nanoparticles, respectively. The terms mnf, hnf,
nf, f, s1, s2, and s3 denote ternary nanofluid, hybrid nanofluid, nanofluid, base fluid,
and the dense nanoparticles of Silver (Ag), Cobalt Ferrite (CoF e2 O4 ), and Zinc Oxide
(ZnO), respectively, as presented in Table 7.1.
The definitions of the skin friction coefficient and local Nusselt number are as
follows:
1
x =
Cf Re1/2 f ′′ (0) (7.10)
(1 − ϕ1 )2.5 (1 − ϕ2 )2.5 (1 − ϕ3 )2.5
kmnf ′
N ux Re−1/2
x =− θ (0) (7.11)
kf
Table 7.1 Thermal and Physical Properties of Ethylene glycol+ water and nanoparticles
(Paul et al. [1, 2], and Ghadikolaei et al. [3, 4])
y(1)
y(2)
y(2)
y(3)
d
y(3) =
y(2)2 −y(1)y(3)−a2 λy(4)+(a3 M +a1 K)y(2)−2a1 γy(3)
(1+2γη)a1
dη
y(4) y(5)
y(5) y(2)(y(4)+δ)−y(1)y(5)−2a4 γy(5)−a5 y(4)
(1+2γη)a 4
thermal boundary layer and a corresponding decrease in fluid velocity. The interplay
between these factors highlights the critical role of thermal stratification in modulating
the flow dynamics, emphasizing its impact on the overall efficiency of heat transfer and
fluid motion in the system. Illustrated in Fig. 7.3, the fluid velocity f ′ (η) falls as the
magnetic parameter (M ) increases due to the Lorentz force, which acts as a resistive
force opposing the fluid motion. This force, generated by the interaction between the
magnetic field and the moving charged particles in the nanofluid, enhances the viscous
drag, reducing momentum transfer and thinning the boundary layer. Fig. 7.4 shows
how the porosity parameter (K) affects the nanofluid velocity f ′ (η). With an increase
in (K), the velocity decreases. There exists an inverse correlation between K and the
porous medium’s diameter, where a higher K corresponds to a decreasing diameter.
This smaller diameter restricts fluid move via the porous medium, reducing velocity
due to the resistance posed by the porosity parameter (K).
As shown in Fig. 7.5, a rise in the value of R results in increased fluid velocity f ′ (η),
as the elevated thermal radiation parameter enhances energy transfer through radiative
heat flux. This intensification reduces the temperature gradient between the fluid
and the surrounding environment, promoting thermal buoyancy effects. Physically,
the added energy accelerates the fluid particles near the boundary, decreasing viscous
resistance and expanding the thermal boundary layer, thereby increasing velocity.
This phenomenon demonstrates how radiation augments convective flow by coupling
thermal and momentum fields, significantly influencing the fluid’s dynamic behavior.
Fig. 7.6 illustrates the impact of λ on fluid velocity, where it is noted that the
velocity increases in tandem with increasing λ. This rise in λ, representing the thermal
buoyancy parameter, enhances the temperature-induced density differences within the
fluid. Physically, these density variations create upward buoyancy forces that overcome
viscous resistance, accelerating the fluid flow. The increase in velocity reflects the
144 MHD Natural Convection of Ternary Nanofluid in a Porous Medium...
dominance of thermal buoyancy over inertial forces, effectively driving the convective
motion. This interaction underscores the critical role of λ in intensifying the coupling
between thermal and momentum transport, resulting in a more vigorous fluid motion
and improved heat transfer efficiency.
Figure 7.7 demonstrates that as thermal stratification (δ) rises, the fluid temperature
declines due to the diminishing temperature gradient between the heated wall and
the adjacent fluid. Physically, thermal stratification creates layered temperature
distributions, restricting heat transfer between these layers. This reduced gradient
7.4 Results and Discussion 145
limits the convective energy transfer from the wall to the fluid, causing a decrease in
thermal energy within the fluid. Additionally, the thickening of the thermal boundary
layer acts as a thermal insulator, further impeding heat diffusion into the bulk fluid.
This interplay of reduced convection and increased insulation leads to a lower overall
temperature in the fluid, highlighting the stabilizing effect of thermal stratification on
thermal transport dynamics. Conversely, Fig. 7.8 shows that the fluid temperature
grows with rising values of R, attributed to the enhancement of thermal radiative flux.
Physically, higher R represents increased radiative heat transfer, which supplements
146 MHD Natural Convection of Ternary Nanofluid in a Porous Medium...
conductive and convective heat transfer within the fluid. The radiation acts as an
additional energy source, intensifying the overall thermal energy available in the system.
This results in elevated fluid temperatures, particularly near the heated surfaces, as
the energy transfer mechanisms synergistically amplify. The direct correlation between
R and temperature underscores the critical role of radiation in augmenting thermal
transport, expanding the thermal boundary layer, and enhancing energy distribution
within the fluid.
7.4 Results and Discussion 147
Figure 7.9 illustrates the effect of Q on the temperature profile, showing that as
Q increases, the fluid temperature rises due to the heat source term represented by
Q. Physically, a higher Q indicates an intensified internal heat generation within
the fluid, which directly adds thermal energy to the system. This increased energy
elevates the temperature of the fluid, particularly in regions closer to the source. The
additional heat enhances the convective transport of energy, thickening the thermal
boundary layer and elevating the overall temperature distribution. This relationship
highlights how internal heat generation drives temperature augmentation, reinforcing
148 MHD Natural Convection of Ternary Nanofluid in a Porous Medium...
the thermal energy balance in the fluid. Figure 7.10 presents the effect of the P r on
the fluid’s temperature profile, revealing that an increase in P r leads to a decrease in
fluid temperature. A higher P r suggests lower thermal conductivity, which reduces
the efficiency of heat transfer through the fluid. Consequently, the rate of heat transfer
diminishes, the thermal boundary layer narrows, and the fluid temperature decreases.
Figures 7.11, 7.12, and 7.13 demonstrate that increasing the volume fractions of
Ag, CoF e2 O4 , and ZnO nanoparticles individually enhances the temperature profile
of the fluid. Physically, each nanoparticle type contributes to improving the thermal
7.4 Results and Discussion 149
conductivity of the nanofluid due to their superior intrinsic heat transfer properties.
Silver (Ag) nanoparticles, known for their exceptional thermal conductivity, enable
rapid heat conduction, while CoF e2 O4 enhances the fluid’s heat retention due to
its magnetic properties, aiding in the uniform distribution of thermal energy. Zinc
oxide (ZnO) adds photocatalytic capabilities that can absorb and redistribute energy
efficiently. The cumulative effect of increased nanoparticle volume fractions reduces
thermal resistance and augments heat transfer within the fluid, raising the overall
temperature profile and improving thermal management capabilities. This highlights
the synergistic role of these nanoparticles in optimizing the thermal behavior of the
ternary nanofluid.
Table 7.3 detail the effects of various flow parameters on skin friction and the
local Nusselt number, offering insights into their physical implications. A rise in the
thermal stratification parameter (δ) reduces the buoyancy forces and temperature
gradient, thereby decreasing skin friction and limiting heat transfer, which lowers the
Nusselt number. Higher values of M (magnetic parameter) introduce Lorentz forces
that oppose fluid motion, reducing momentum transfer and consequently lowering
both skin friction and the Nusselt number. Similarly, higher K (porosity parameter)
increases resistance within the porous medium, dampening flow and thermal transfer.
In contrast, increases in λ (thermal buoyancy) enhance temperature-induced density
differences, driving stronger flow and improving heat transfer efficiency, which raises
both skin friction and the Nusselt number. For P r (Prandtl number), a higher value
indicates lower thermal diffusivity, leading to a stronger thermal boundary layer and
improved heat transfer (higher Nusselt number) despite reduced flow momentum (lower
skin friction). An increase in R (radiation parameter) adds thermal energy, raising
fluid velocity and skin friction but reducing the heat transfer rate due to a diminished
temperature gradient. Similarly, higher Q (heat source) increases internal heating,
150 MHD Natural Convection of Ternary Nanofluid in a Porous Medium...
thickening the boundary layer and raising skin friction while decreasing heat transfer
efficiency. Finally, an increase in γ (curvature parameter) weakens flow resistance,
reducing skin friction, while simultaneously enhancing convective heat transfer, which
raises the Nusselt number. These relationships highlight the intricate balance between
flow dynamics and heat transfer mechanisms in the system.
δ M K λ Pr R Q γ Cf Rex1/2 N ux Rex−1/2
0 -3.0085 2.4733
0.15 0.5 1 1.5 6.2 2 0.06 0.3 -3.1192 2.3167
0.25 -3.1919 2.2079
0 -2.7456 2.4952
0.05 0.6 1 1.5 6.2 2 0.06 0.3 -3.1026 2.4084
1.2 -3.4268 2.3306
0 -2.2325 2.6210
0.05 0.5 0.5 1.5 6.2 2 0.06 0.3 -2.6659 2.5147
1.5 -3.3879 2.3398
0 -3.5708 2.2467
0.05 0.5 1 1 6.2 2 0.06 0.3 -3.2131 2.3751
2 -2.8834 2.4630
0.6 -2.7941 0.8328
0.05 0.5 1 1.5 6.2 2 0.06 0.3 -3.0457 2.4221
21 -3.2519 4.8521
0 -3.1632 3.5724
0.05 0.5 1 1.5 6.2 1.5 0.06 0.3 -3.0681 2.6106
3 -3.0094 2.1393
-0.1 -3.0792 2.6734
0.05 0.5 1 1.5 6.2 2 0 0.3 -3.0598 2.5225
0.05 -3.0482 2.4395
0 -2.7907 2.2332
0.05 0.5 1 1.5 6.2 2 0.06 0.25 -3.0037 2.3898
0.50 -3.2111 2.5519
Tables 7.4 and 7.5 highlight the performance of ternary nanofluid formulations
in terms of skin friction (−Cf Re1/2 x ) and the local Nusselt number (N ux Rex
−1/2
),
demonstrating their superior thermal and flow characteristics compared to conventional
and hybrid nanofluids. The enhanced skin friction of the ternary nanofluid can be
attributed to its higher viscosity and density, resulting from the synergistic properties of
the nanoparticles (Ag, CoF e2 O4 , ZnO), which increase resistance to flow but improve
heat retention. As thermal stratification (δ) intensifies, the temperature gradient
diminishes, narrowing the disparity in skin friction between the ternary nanofluid
and the standard nanofluid by reducing buoyancy-driven momentum transfer. The
7.5 Conclusion 151
significant improvement in the Nusselt number for the ternary nanofluid is due to its
superior thermal conductivity, enabled by the optimal dispersion of nanoparticles, which
facilitates efficient heat transfer. Additionally, as δ increases, the thermal conductivity
advantage of the ternary nanofluid becomes more pronounced, further amplifying the
heat transfer rate compared to conventional nanofluids. This indicates the ability
of ternary nanofluids to maintain high thermal performance even under stratified
conditions, underscoring their potential for advanced heat transfer applications.
7.5 Conclusion
A comprehensive investigation was conducted to analyze the effects of thermal stratifica-
tion on the flow of a ternary nanofluid with magnetohydrodynamics (MHD) properties
along a vertically stretchable cylinder. The study incorporates the influence of thermal
buoyancy, heat source/sink, radiation, and the presence of a porous medium. The
analysis examines the flow characteristics, focusing on their impact on velocity and
temperature profiles, skin friction, and the local Nusselt number. The key findings of
this study are summarized as follows:
• The absolute skin friction is increased for parameters λ and R while it decreases
for parameters δ, M, K, Q, P r and γ.
152 MHD Natural Convection of Ternary Nanofluid in a Porous Medium...
• The ternary nanofluid shows greater absolute skin friction when compared to
the hybrid nanofluid. Likewise, the skin friction of hybrid nanofluids shows an
increase in comparison to that of regular nanofluids.
• The ternary nanofluid shows a higher heat transfer rate in comparison to the
hybrid nanofluid, and the hybrid nanofluid exhibits a higher heat transfer rate
when compared to traditional nanofluids.
Future Work
8
The research presented in this thesis has explored various aspects of heat and mass
transfer in stratified fluids under the influence of magnetohydrodynamics (MHD) and
porous media. While significant insights have been gained, there are several avenues
for future research to further enhance the understanding and application of these
phenomena. Below are some potential directions for future work:
• Advanced Nanofluids
Ternary and Hybrid Nanofluids: While the current study explores the be-
havior of ternary nanofluids, future work could investigate the effects of different
nanoparticle combinations, shapes, and concentrations on heat and mass transfer.
Additionally, the impact of nanoparticle aggregation, stability, and sedimenta-
tion in nanofluids could be studied to optimize their performance in various
applications.
• Experimental Validation
Laboratory and Field Experiments: While the current study provides theo-
retical and numerical insights, future work could include experimental validation
of the results. Laboratory experiments and field studies could be conducted to
validate the theoretical models and provide empirical data for further refinement
of the models.
• Multi-Physics Coupling
Coupled Thermal, Mechanical, and Electromagnetic Effects: Future
research could explore the coupling of thermal, mechanical, and electromagnetic
effects in stratified fluids. This would involve studying the interactions between
thermal gradients, fluid flow, and electromagnetic fields in multi-physics systems,
which are relevant in applications such as energy storage, geothermal systems,
and advanced manufacturing processes.
The future work outlined above aims to build on the foundational research presented
in this thesis, addressing some of the limitations and expanding the scope of the study.
By exploring these directions, future research can further enhance our understanding of
heat and mass transfer in stratified fluids, leading to improved designs and optimizations
in various engineering and environmental applications.
80_Recommendation
“I would rather have questions that can’t be answered than answers that
can’t be questioned.” -Richard P. Feynman
This thesis offers a detailed and methodical study of heat and mass transfer in
MHD flows with stratification and chemical reactions through porous media. It
effectively integrates analytical methods (Laplace transform) and numerical approaches
(MATLAB’s bvp4c solver) to solve complex physical models. Each chapter builds
coherently on the last, deepening the understanding of thermal and mass stratification
effects alongside radiation, chemical reactions, and nanofluid dynamics. The innovative
modeling of MHD ternary nanofluid flow adds significant novelty, making this work
relevant for modern thermal engineering applications.
While the thesis demonstrates depth and breadth, there is room for enhance-
ment. Complex nonlinear scenarios could benefit from experimental or high-fidelity
numerical validation to improve real-world relevance. The presentation, though clear,
may be refined by reducing repetition and better integrating insights across chapters.
Additionally, the future directions in Chapter 8—such as data-driven modeling and
CFD simulation—highlight the potential for this research to expand into more ap-
plied, interdisciplinary fields, strengthening its impact in both academic and practical
contexts.
[1] Ashish Paul, Bhagyashri Patgiri, and Neelav Sarma. Darcy-forchheimer flow of
ag–zno–cofe2o4/h2o casson ternary hybrid nanofluid induced by a rotatory disk
with emhd. International Journal of Ambient Energy, 45(1):2313697, 2024.
[2] Ashish Paul, Neelav Sarma, and Bhagyashri Patgiri. Numerical assessment of
mhd thermo-mass flow of casson ternary hybrid nanofluid around an exponentially
stretching cylinder. BioNanoScience, pages 1–16, 2024.
[3] SS Ghadikolaei and M Gholinia. 3d mixed convection mhd flow of go-mos2
hybrid nanoparticles in h2o–(ch2oh) 2 hybrid base fluid under the effect of h2
bond. International Communications in Heat and Mass Transfer, 110:104371,
2020.
[4] SS Ghadikolaei, M Gholinia, ME Hoseini, and DD Ganji. Natural convection
mhd flow due to mos2–ag nanoparticles suspended in c2h6o2h2o hybrid base fluid
with thermal radiation. Journal of the Taiwan Institute of Chemical Engineers,
97:12–23, 2019.
[5] Alan Shapiro and Evgeni Fedorovich. Unsteady convectively driven flow along a
vertical plate immersed in a stably stratified fluid. Journal of Fluid Mechanics,
498:333–352, 2004.
[6] S US Choi and Jeffrey A Eastman. Enhancing thermal conductivity of fluids
with nanoparticles. Technical report, Argonne National Lab.(ANL), Argonne, IL
(United States), 1995.
[7] Robert Siegel. Transient free convection from a vertical flat plate. Transactions
of the American Society of Mechanical Engineers, 80(2):347–357, 1958.
[8] ER Menold and Kwang-Tzu Yang. Asymptotic solutions for unsteady laminar
free convection on a vertical plate. J. Appl. Mech., 1962.
[9] JA Schetz and R Eichhorn. Unsteady natural convection in the vicinity of a
doubly infinite vertical plate. J. Heat Transfer, 1962.
[10] MA Hossain and HS Takhar. Radiation effect on mixed convection along a vertical
plate with uniform surface temperature. Heat and Mass transfer, 31(4):243–248,
1996.
[11] UN Das, RK Deka, and VM Soundalgekar. Transient free convection flow past
an infinite vertical plate with periodic temperature variation. J. Heat Transfer.,
1999.
References References
[12] MA Hossain, MA Alim, and D Andrew S Rees. The effect of radiation on free
convection from a porous vertical plate. International Journal of Heat and Mass
Transfer, 42(1):181–191, 1999.
[13] Jun Sang Park and Jae Min Hyun. Technical note transient behavior of vertical
buoyancy layer in a stratified fluid. International journal of heat and mass
transfer, 41(24):4393–4397, 1998.
[14] Jun Sang Park. Transient buoyant flows of a stratified fluid in a vertical channel.
KSME international journal, 15:656–664, 2001.
[15] UN Das, R Deka, and VM Soundalgekar. Effects of mass transfer on flow past an
impulsively started infinite vertical plate with constant heat flux and chemical
reaction. Forschung im Ingenieurwesen, 60(10):284–287, 1994.
[16] E Magyari, I Pop, and B Keller. Unsteady free convection along an infinite
vertical flat plate embedded in a stably stratified fluid-saturated porous medium.
Transport in porous media, 62:233–249, 2006.
[17] Anjana Bhattacharya and RK Deka. Theoretical study of chemical reaction
effects on vertical oscillating plate immersed in a stably stratified fluid. Research
Journal of Applied Sciences, Engineering and Technology, 3(9):887–898, 2011.
[18] V Rajesh and S Vijaya Kumar Varma. Radiation and mass transfer effects
on mhd free convection flow past an exponentially accelerated vertical plate
with variable temperature. ARPN Journal of Engineering and Applied Sciences,
4(6):20–26, 2009.
[19] AJ Chamkha, MC Raju, T Sudhakar Reddy, and SVK Varma. Unsteady mhd
free convection flow past an exponentially accelerated vertical plate with mass
transfer, chemical reaction and thermal radiation. International Journal of
Microscale and Nanoscale Thermal and Fluid Transport Phenomena, 5(1):57,
2014.
[20] Suresh Babu Ramakrishna, Ramesh NL, and Sravan Kumar Thavada. Effects of
chemical reaction, soret and lorentz force on casson fluid flow past an exponentially
accelerated vertical plate: A comprehensive analysis. Heat Transfer, 51(2):2237–
2257, 2022.
[21] R Muthucumaraswamy and P Sivakumar. Hydro magnetic effects on parabolic
flow past an infinite isothermal vertical plate with variable mass diffusion in the
presence of thermal radiation and chemical reaction. International Journal of
Recent Technology and Engineering, 4(2):5–10, 2015.
[22] A Selvaraj, S Dilip Jose, R Muthucumaraswamy, and S Karthikeyan. Mhd-
parabolic flow past an accelerated isothermal vertical plate with heat and mass
diffusion in the presence of rotation. Materials Today: Proceedings, 46:3546–3549,
2021.
References 161
[23] Rajesh Vemula, Lokenath Debnath, and Sridevi Chakrala. Unsteady mhd free
convection flow of nanofluid past an accelerated vertical plate with variable
temperature and thermal radiation. International Journal of Applied and Com-
putational Mathematics, 3(2):1271–1287, 2017.
[24] Shobha Poovani Jeevandhar, Vidya Kedla, Neeraja Gullapalli, and Sravan Kumar
Thavada. Natural convective effects on mhd boundary layer nanofluid flow over
an exponentially accelerating vertical plate. Biointerface Research in Applied
Chemistry, 11(6):13790–13805, 2021.
[25] EM Sparrow and JL Gregg. Laminar-free-convection heat transfer from the outer
surface of a vertical circular cylinder. Transactions of the American Society of
Mechanical Engineers, 78(8):1823–1828, 1956.
[26] Kwang-Tzu Yang. Possible similarity solutions for laminar free convection on
vertical plates and cylinders. J. Applied Mechanics, 1960.
[27] Richard J Goldstein and DG Briggs. Transient free convection about vertical
plates and circular cylinders. J. Heat Transfer-T ASME, 1964.
[28] P Ganesan and P Loganathan. Effects of mass transfer and flow past a moving
vertical cylinder with constant heat flux. Acta mechanica, 150(3):179–190, 2001.
[29] P Ganesan and P Loganathan. Radiation and mass transfer effects on flow of
an incompressible viscous fluid past a moving vertical cylinder. International
Journal of Heat and Mass Transfer, 45(21):4281–4288, 2002.
[30] P Ganesan and P Loganathan. Unsteady natural convective flow past a moving
vertical cylinder with heat and mass transfer. Heat and Mass Transfer, 37(1):59–
65, 2001.
[31] Harmindar S Takhar, Ali J Chamkha, and Girishwar Nath. Combined heat and
mass transfer along a vertical moving cylinder with a free stream. Heat and Mass
transfer, 36:237–246, 2000.
[32] HP Rani. Transient natural convection along a vertical clyinder with variable
surface temperature and mass diffusion. Heat and mass transfer, 40(1):67–73,
2003.
[33] Rudra Kanta Deka and Ashish Paul. Convectively driven flow past an infinite
moving vertical cylinder with thermal and mass stratification. Pramana, 81:641–
665, 2013.
[34] Ashish Paul and Rudra Kanta Deka. Unsteady natural convection flow past
an infinite cylinder with thermal and mass stratification. Int. J. Eng. Math,
2017:8410691, 2017.
[35] Jae Min Hyun. Transient process of thermally stratifying an initially homoge-
neous fluid in an enclosure. International Journal of Heat and Mass Transfer,
27(10):1936–1938, 1984.
References References
[36] Wenxian Lin and SW Armfield. Natural convection cooling of rectangular and
cylindrical containers. International Journal of Heat and Fluid Flow, 22(1):72–81,
2001.
[37] R Kandasamy, K Periasamy, and KK Sivagnana Prabhu. Chemical reaction,
heat and mass transfer on mhd flow over a vertical stretching surface with heat
source and thermal stratification effects. International Journal of Heat and Mass
Transfer, 48(21-22):4557–4561, 2005.
[38] MA Mansour, NF El-Anssary, and AM Aly. Effects of chemical reaction and
thermal stratification on mhd free convective heat and mass transfer over a
vertical stretching surface embedded in a porous media considering soret and
dufour numbers. Chemical engineering journal, 145(2):340–345, 2008.
[39] Ching-Yang Cheng. Double-diffusive natural convection along a vertical wavy
truncated cone in non-newtonian fluid saturated porous media with thermal and
mass stratification. International Communications in Heat and Mass Transfer,
35(8):985–990, 2008.
[40] Ching-Yang Cheng. Combined heat and mass transfer in natural convection
flow from a vertical wavy surface in a power-law fluid saturated porous medium
with thermal and mass stratification. International Communications in Heat and
Mass Transfer, 36(4):351–356, 2009.
[41] D Srinivasacharya and O Surender. Double stratification effects on mixed
convection along a vertical plate in a non-darcy porous medium. Procedia
Engineering, 127:986–993, 2015.
[42] Darbhasayanam Srinivasacharya and O Surender. Non-darcy mixed convection
induced by a vertical plate in a doubly stratified porous medium. Journal of
Porous Media, 17(3), 2014.
[43] D Srinivasacharya and O Surender. Mixed convection over a vertical plate in a
doubly stratified fluid-saturated non-darcy porous medium with cross-diffusion
effects. Heat Transfer—Asian Research, 43(8):677–690, 2014.
[44] D Srinivasacharya, Sandile S Motsa, and O Surender. Numerical study of
free convection in a doubly stratified non-darcy porous medium using spectral
quasilinearization method. International Journal of Nonlinear Sciences and
Numerical Simulation, 16(3-4):173–183, 2015.
[45] D Srinivasacharya, Ali Chamkha, O Surender, and AM Rashad. Natural convec-
tion on a porous vertical plate in a doubly stratified non-darcy porous medium.
Frontiers in Heat and Mass Transfer (FHMT), 6(1), 2015.
[46] D Srinivasacharya and O Surender. Mixed convection boundary layer flow of a
nanofluid past a vertical plate in a doubly stratified porous medium. Journal of
Computational and Theoretical Nanoscience, 11(8):1853–1862, 2014.
[47] D Srinivasacharya and O Surender. Effect of double stratification on natural
convective boundary layer flow over a vertical plate in a porous medium saturated
with nanofluid. Journal of Nanofluids, 4(3):318–327, 2015.
References 163
[59] Hakan F Oztop and Eiyad Abu-Nada. Numerical study of natural convection
in partially heated rectangular enclosures filled with nanofluids. International
journal of heat and fluid flow, 29(5):1326–1336, 2008.
[60] S Das and RN Jana. Natural convective magneto-nanofluid flow and radiative
heat transfer past a moving vertical plate. Alexandria Engineering Journal,
54(1):55–64, 2015.
[61] S Das, RN Jana, and OD Makinde. Transient natural convection in a vertical
channel filled with nanofluids in the presence of thermal radiation. Alexandria
Engineering Journal, 55(1):253–262, 2016.
[62] Mohammad Hossein Abolbashari, Navid Freidoonimehr, Foad Nazari, and Mo-
hammad Mehdi Rashidi. Entropy analysis for an unsteady mhd flow past a
stretching permeable surface in nano-fluid. Powder Technology, 267:256–267,
2014.
[63] P Chandra Reddy, MC Raju, and GSS Raju. Free convective heat and mass
transfer flow of heat-generating nanofluid past a vertical moving porous plate in
a conducting field. Special Topics & Reviews in Porous Media: An International
Journal, 7(2), 2016.
[64] B Mahanthesh, BJ Gireesha, and Rama Subba Reddy Gorla. Heat and mass
transfer effects on the mixed convective flow of chemically reacting nanofluid past
a moving/stationary vertical plate. Alexandria engineering journal, 55(1):569–581,
2016.
[65] Umair Rashid, Ali Akgül, and Dianchen Lu. Impact of nanosized particles on
hybrid nanofluid flow in porous medium with thermal slip condition. Numerical
Heat Transfer, Part B: Fundamentals, pages 1–14, 2023.
[66] Jianfeng Wang, Shan Ali Khan, Sumeira Yasmin, Mohammad Mahtab Alam,
Haihu Liu, Umar Farooq, Ali Akgül, and Ahmed M Hassan. Central composite
design (ccd)-response surface methodology (rsm) for modeling and simulation of
mwcnt-water nanofluid inside hexagonal cavity: Application to electronic cooling.
Case Studies in Thermal Engineering, 50:103488, 2023.
[67] Sawan Kumar Rawat, Moh Yaseen, Umair Khan, Manoj Kumar, Amal Abdul-
rahman, Sayed M Eldin, Samia Elattar, Ahmed M Abed, and Ahmed M Galal.
Insight into the significance of nanoparticle aggregation and non-uniform heat
source/sink on titania–ethylene glycol nanofluid flow over a wedge. Arabian
Journal of Chemistry, 16(7):104809, 2023.
[68] Moh Yaseen, Sawan Kumar Rawat, Umair Khan, Anup Singh Negi, Manoj
Kumar, El-Sayed M Sherif, Ahmed M Hassan, and Ioan Pop. Inspection of
unsteady buoyancy and stagnation point flow incorporated by ag-tio2 hybrid
nanoparticles towards a spinning disk with hall effects. Case Studies in Thermal
Engineering, 44:102889, 2023.
References 165
[69] Moh Yaseen, Sawan Kumar Rawat, Umair Khan, Ioannis E Sarris, Humera
Khan, Anup Singh Negi, Arshad Khan, El-Sayed M Sherif, Ahmed M Hassan,
and Aurang Zaib. Numerical analysis of magnetohydrodynamics in an eyring–
powell hybrid nanofluid flow on wall jet heat and mass transfer. Nanotechnology,
34(48):485405, 2023.
[70] Ayesha Riaz, Muhammad Imran Asjad, Musawa Yahya Almusawa, Sayed M
Eldin, Nejla Mahjoub Said, and Muhammad Ijaz Khan. New solutions of
fractional maxwell fluid with ternary-hybrid nanoparticles. Case Studies in
Thermal Engineering, 47:103091, 2023.
[71] M Muralidharn and R Muthucumaraswamy. Parabolic started flow past an
infinite vertical plate with uniform heat flux and variable mass diffusion. Int.
Journal of Math. Analysis, 8(26):1265–1274, 2014.
[72] RC Chaudhary and Arpita Jain. Mhd heat and mass diffusion flow by natural
convection past a surface embedded in a porous medium. Theoretical and Applied
Mechanics, 36(1):1–27, 2009.
[73] Yanala Dharmendar Reddy, B Shankar Goud, and M Anil Kumar. Radiation
and heat absorption effects on an unsteady mhd boundary layer flow along an
accelerated infinite vertical plate with ramped plate temperature in the existence
of slip condition. Partial Differential Equations in Applied Mathematics, 4:100166,
2021.
[74] A Selvaraj and E Jothi. Heat source impact on mhd and radiation absorption
fluid flow past an exponentially accelerated vertical plate with exponentially
variable temperature and mass diffusion through a porous medium. Materials
Today: Proceedings, 46:3490–3494, 2021.
[75] Subhrajit Sarma and Nazibuddin Ahmed. Thermal diffusion effect on unsteady
mhd free convective flow past a semi-infinite exponentially accelerated vertical
plate in a porous medium. Canadian Journal of Physics, 100(10):437–451, 2022.
[76] M Veera Krishna and Ali J Chamkha. Hall and ion slip effects on mhd rotating flow
of elastico-viscous fluid through porous medium. International Communications
in Heat and Mass Transfer, 113:104494, 2020.
[77] M Veera Krishna, Kamboji Jyothi, and Ali J Chamkha. Heat and mass transfer
on mhd flow of second-grade fluid through porous medium over a semi-infinite
vertical stretching sheet. Journal of Porous media, 23(8), 2020.
[78] M Veera Krishna, N Ameer Ahamad, and Ali J Chamkha. Hall and ion slip
effects on unsteady mhd free convective rotating flow through a saturated porous
medium over an exponential accelerated plate. Alexandria Engineering Journal,
59(2):565–577, 2020.
[79] M Veera Krishna and Ali J Chamkha. Hall and ion slip effects on magnetohy-
drodynamic convective rotating flow of jeffreys fluid over an impulsively moving
vertical plate embedded in a saturated porous medium with ramped wall temper-
ature. Numerical Methods for Partial Differential Equations, 37(3):2150–2177,
2021.
References References
[80] M Veera Krishna and Ali J Chamkha. Hall and ion slip effects on mhd rotating
boundary layer flow of nanofluid past an infinite vertical plate embedded in a
porous medium. Results in Physics, 15:102652, 2019.
[81] Nitul Kalita, Rudra Kanta Deka, and Rupam Shankar Nath. Unsteady flow past
an accelerated vertical plate with variable temperature in presence of thermal
stratification and chemical reaction. East European Journal of Physics, 3:441–450,
2023.
[82] Rupam Shankar Nath and Rudra Kanta Deka. The effects of thermal stratification
on flow past an infinite vertical plate in presence of chemical reaction. East
European Journal of Physics, 3:223–232, 2023.
[83] Rupam Shankar Nath, Rudra Kanta Deka, and Himangshu Kumar. The effect
of thermal stratification on unsteady parabolic flow past an infinite vertical plate
with chemical reaction. East European Journal of Physics, 4:77–86, 2023.
[84] Himangshu Kumar and Rudra Kanta Deka. Thermal and mass stratification
effects on unsteady flow past an accelerated infinite vertical plate with variable
temperature and exponential mass diffusion in porous medium. East European
Journal of Physics, 4:87–97, 2023.
[85] Subhrajit Sarma and Nazibuddin Ahmed. Dufour effect on unsteady mhd flow
past a vertical plate embedded in porous medium with ramped temperature.
Scientific Reports, 12(1):13343, 2022.
[86] M. Abramowitz and IA Stegun. Handbook of mathematical functions with For-
mulas, Graphs, and Mathematical Tables” edited by Dover Publications. Inc.,
New York, Ninth Printing, 1970.
[87] Richard B Hetnarski. An algorithm for generating some inverse laplace transforms
of exponential form. Zeitschrift für angewandte Mathematik und Physik ZAMP,
26:249–253, 1975.
[88] GV Ramana Reddy, Ch V Ramana Murthy, and N Bhaskar Reddy. Unsteady
mhd free convective mass transfer flow past an infinite vertical porous plate with
variable suction and soret effect. Int. J. of Appl. Math and Mech, 7(21):70–84,
2011.
[89] G Kumaran and N Sandeep. Thermophoresis and brownian moment effects on
parabolic flow of mhd casson and williamson fluids with cross diffusion. Journal
of Molecular Liquids, 233:262–269, 2017.
[90] P Mohan Krishna, N Sandeep, and Ram Prakash Sharma. Computational
analysis of plane and parabolic flow of mhd carreau fluid with buoyancy and
exponential heat source effects. The European Physical Journal Plus, 132:1–15,
2017.
[91] T Gowri and A Selvaraj. Rotational effect of unsteady mhd-parabolic flow past a
vertical plate through porous medium with uniform temperature mass diffusion.
International Journal of Mechanical Engineering, 6:2468–2473, 2021.
References 167
[92] Mehari Fentahun Endalew and Subharthi Sarkar. Incidences of aligned magnetic
field on unsteady mhd flow past a parabolic accelerated inclined plate in a porous
medium. Heat Transfer, 50(6):5865–5884, 2021.
[93] Poli Chandra Reddy, Mummadisetty Umamheswar, S Harinath Reddy, AB Mad-
humohana Raju, and Malaraju Chengal Raju. Numerical study on the parabolic
flow of mhd fluid past a vertical plate in a porous medium. Heat Transfer,
51(4):3418–3430, 2022.
[94] PK Sharma and RC Chaudhary. Effect of variable suction on transient free
convection viscous incompressible flow past a vertical plate with periodic temper-
ature variations in slip-flow regime. Emirates Journal of Engineering Research,
8(2):33–38, 2003.
[95] Shusaku Goto, Makoto Yamano, and Masataka Kinoshita. Thermal response of
sediment with vertical fluid flow to periodic temperature variation at the surface.
Journal of Geophysical Research: Solid Earth, 110(B1), 2005.
[96] Ming-Cong Li, Ying Song, Wei Shen, Cong Wang, Wei-Kang Qi, Yongzhen Peng,
and Yu-You Li. The performance of an anaerobic ammonium oxidation upflow
anaerobic sludge blanket reactor during natural periodic temperature variations.
Bioresource Technology, 293:122039, 2019.
[97] Yikun Jia, Ningzhe Xi, Rongrong Li, Jiazhen Fan, Wenyan Zhao, Tao Zhang,
Jianhong Li, and Junbing Pu. Aquatic photosynthetic carbon pump driven by
periodic temperature variations affects dissolved inorganic carbon and hydro-
chemical characteristics in a groundwater-fed reservoir. Environmental Research,
252:119041, 2024.
[98] M Veera Krishna, K Bharathi, and Ali J Chamkha. Hall effects on mhd peristaltic
flow of jeffrey fluid through porous medium in a vertical stratum. Interfacial
Phenomena and Heat Transfer, 6(3), 2018.
[99] M Veera Krishna and Ali J Chamkha. Hall effects on mhd squeezing flow of
a water-based nanofluid between two parallel disks. Journal of Porous media,
22(2), 2019.
[100] M Veera Krishna, PVS Anand, and Ali J Chamkha. Heat and mass transfer on
free convective flow of amicropolar fluid through a porous surface with inclined
magnetic field and hall effects. Special Topics & Reviews in Porous Media: An
International Journal, 10(3), 2019.
[101] M Veera Krishna, BV Swarnalathamma, and Ali J Chamkha. Investigations of
soret, joule and hall effects on mhd rotating mixed convective flow past an infinite
vertical porous plate. Journal of Ocean Engineering and Science, 4(3):263–275,
2019.
[102] M Veera Krishna and Ali J Chamkha. Hall and ion slip effects on unsteady mhd
convective rotating flow of nanofluids—application in biomedical engineering.
Journal of the Egyptian Mathematical Society, 28(1):1, 2020.
References References
[103] M Veera Krishna. Hall and ion slip impacts on unsteady mhd free convective
rotating flow of jeffreys fluid with ramped wall temperature. International
Communications in Heat and Mass Transfer, 119:104927, 2020.
[104] M Veera Krishna. Hall and ion slip effects on radiative mhd rotating flow of
jeffreys fluid past an infinite vertical flat porous surface with ramped wall velocity
and temperature. International communications in Heat and Mass transfer,
126:105399, 2021.
[105] Masood Khan, Mahnoor Sarfraz, Jawad Ahmed, Latif Ahmad, and Awais Ahmed.
Viscoelastic nanofluid motion for homann stagnation-region with thermal radia-
tion characteristics. Proceedings of the Institution of Mechanical Engineers, Part
C: Journal of Mechanical Engineering Science, 235(21):5324–5336, 2021.
[106] Mahnoor Sarfraz and Masood Khan. Entropy generation analysis of cnt-based
nanofluid flows induced by a moving plate. ZAMM-Journal of Applied Mathe-
matics and Mechanics/Zeitschrift für Angewandte Mathematik und Mechanik,
104(1):e202200555, 2024.
[107] Mahnoor Sarfraz and Masood Khan. Significance of ethylene glycol-based cnt
homann nanofluid flow over a biaxially stretching surface. Waves in Random
and Complex Media, pages 1–15, 2022.
[108] Mahnoor Sarfraz, Masood Khan, A Al-Zubaidi, and S Saleem. Insights into the
thermodynamic efficiency of homann-agrawal hybrid nanofluid flow. Alexandria
Engineering Journal, 82:178–185, 2023.
[109] Jawad Ahmed, S Bourazza, Mahnoor Sarfraz, MA Orsud, Sayed M Eldin, Nadia A
Askar, and Mohamed Abdelghany Elkotb. Heat transfer in jeffrey fluid flow over
a power law lubricated surface inspired by solar radiations and magnetic flux.
Case Studies in Thermal Engineering, 49:103220, 2023.
[110] Masood Khan, Mahnoor Sarfraz, and Rida Zehra. Energy transport near homann
stagnation point flow over a spiraling disk with cattaneo–christov theory. Inter-
national Journal of Modern Physics B, 36(25):2250171, 2022.
[111] Khursheed Muhammad, Zahid Nisar, Abdullah Alhuthali, and Mohamed Hussien.
Statistical and numerical analysis of electrically conducting hybrid nanomaterial
near the stagnation region. Journal of Thermal Analysis and Calorimetry, pages
1–10, 2024.
[112] Mahnoor Sarfraz and Masood Khan. Rheology of gyrotactic microorganisms in
jeffrey fluid flow: a stability analysis. Modern Physics Letters B, 38(05):2450003,
2024.
[113] V Puneeth, Farhan Ali, M Riaz Khan, M Shoaib Anwar, and N Ameer Ahammad.
Theoretical analysis of the thermal characteristics of ree–eyring nanofluid flowing
past a stretching sheet due to bioconvection. Biomass Conversion and Biorefinery,
14(7):8649–8660, 2024.
References 169
[114] Muhammad Irfan, Kiran Rafiq, Muhammad Shoaib Anwar, Masood Khan,
Waqar Azeem Khan, and Kaleem Iqbal. Evaluating the performance of new mass
flux theory on carreau nanofluid using the thermal aspects of convective heat
transport. Pramana, 95:1–9, 2021.
[115] Zakir Hussain, Saima Muhammad, and Muhammad Shoaib Anwar. Effects of
first-order chemical reaction and melting heat on hybrid nanoliquid flow over a
nonlinear stretched curved surface with shape factors. Advances in Mechanical
Engineering, 13(4):1687814021999526, 2021.
[116] Muhammad Shoaib Anwar, Taseer Muhammad, Mumtaz Khan, and V Puneeth.
Mhd nanofluid flow through darcy medium with thermal radiation and heat
source. International Journal of Modern Physics B, 38(28):2450386, 2024.
[117] MS Anwar. Modeling and numerical simulations of some fractional nonlinear
viscoelastic flow problems. PhD diss. Lahore University of Management of
Science, 2019.
[118] Muhammad Shoaib Anwar, Metib Alghamdi, Taseer Muhammad, M Hussain,
and V Puneeth. Analysis of nonlinear convection and diffusion in viscoelastic
fluid flow with variable thermal conductivity and thermal radiations. Modern
Physics Letters B, 38(22):2450146, 2024.
[119] Majid Hussain, Aqsa Lubna, Mehwish Ashraf, MS Anwar, Qasim A Ranjha, and
Akhtar Ali. Ohmically dissipated mhd mixed convective flow of williamson fluid
over a penetrable stretching convective wedge with thermal radiations. Numerical
Heat Transfer, Part B: Fundamentals, 85(8):1026–1040, 2024.
[120] Majid Hussain, Muhammad Shoaib, Qasim A. Ranjha, Muhammad Shoaib Anwar,
Zahoor Ahmad, and Mustafa Inc. Numerical solution to flow of casson fluid
via stretched permeable wedge with chemical reaction and mass transfer effects.
Modern Physics Letters B, 38(16):2341008, 2024.
[121] Ali Haider, MS Anwar, Yufeng Nie, and MS Alqarni. Exploring the role of
nanoparticles in enhancing thermal performance of fractional carreau fluid flow
under magnetic field and hall current. Case Studies in Thermal Engineering,
59:104537, 2024.
[122] Muhammad Irfan, Muhammad Shoaib Anwar, Metib Alghamdi, Masood Khan,
and Taseer Muhammad. Modeling heat-mass transport for mhd bio-convection
carreau nanofluid with joule heating containing both gyrotactic microbes and
nanoparticles diffusion. ZAMM-Journal of Applied Mathematics and Mechan-
ics/Zeitschrift für Angewandte Mathematik und Mechanik, 104(10):e202400234,
2024.
[123] Digbash Sahu and Rudra Kanta Deka. Thermal and mass stratification effects
on mhd flow past an accelerated vertical plate with variable temperature and
exponential mass diffusion embedded in a porous medium. East European Journal
of Physics, 2:161–171, 2024.
References References
[124] GS Seth, SM Hussain, and S Sarkar. Hydromagnetic natural convection flow with
heat and mass transfer of a chemically reacting and heat absorbing fluid past an
accelerated moving vertical plate with ramped temperature and ramped surface
concentration through a porous medium. Journal of the Egyptian Mathematical
Society, 23(1):197–207, 2015.
[125] Kanayo K Asogwa, Sardar M Bilal, Isaac L Animasaun, and Fateh M Mebarek-
Oudina. Insight into the significance of ramped wall temperature and ramped
surface concentration: The case of casson fluid flow on an inclined riga plate with
heat absorption and chemical reaction. Nonlinear Engineering, 10(1):213–230,
2021.
[126] Dharmendar Reddy Yanala, Anil Kumar Mella, Srinivasa Rao Vempati, and
B Shankar Goud. Influence of slip condition on transient laminar flow over
an infinite vertical plate with ramped temperature in the presence of chemical
reaction and thermal radiation. Heat Transfer, 50(8):7654–7671, 2021.
[127] HK Mandal, DK Maiti, and RN Jana. Unsteady mhd free convection in a
radiating fluid flow past a vertically time-dependent moving plate with ramped
double-diffusive condition. Journal of Engineering Thermophysics, 33(3):598–621,
2024.
[128] Rengasamy Rajaraman and Rajamanickam Muthucumaraswamy. Impact of
chemical reaction, viscous dissipation, and thermal radiation on mhd flow past
an oscillating plate: Original scientific paper. Chemical Industry & Chemical
Engineering Quarterly, 30(3):223–230, 2024.
[129] A Aravind and J Ravikumar. Effects of thermal radiation in a rotating fluid
on an oscillating vertical plate with variable temperature and mass diffusion.
IAENG International Journal of Applied Mathematics, 54(1), 2024.
[130] Zia Ullah and Mohammed Alkinidri. Effect of variable viscosity on oscillatory
heat and mass transfer in mixed convective flow with chemical reaction along
inclined heated plate under reduced gravity. Alexandria Engineering Journal,
77:539–552, 2023.
[131] Rajdeep Bordoloi, Dipunja Gohain, Nazibuddin Ahmed, and Ali J Chamkha. An
exact analysis of radiation absorption and dufour effect on mhd convective flow
of cu-water nanofluid with heat generation and chemical reaction. International
Journal of Modern Physics B, 38(27):2450371, 2024.
[132] Rajdeep Bordoloi, Kalyan Chamuah, and Nazibuddin Ahmed. Free convective
mhd radiative flow past a porous vertical plate in a porous medium with chemical
reaction. Biointerface Research in Applied Chemistry, 13(3):1–15, 2023.
[133] Vasa Vijaya Kumar, Mamidi Narsimha Raja Shekar, and Shankar Goud Bejawada.
Heat and mass transfer significance on mhd flow over a vertical porous plate in
the presence of chemical reaction and heat generation. CFD Letters, 16(5):9–20,
2024.
References 171
[134] B Shankar Goud, Y Dharmendar Reddy, and Kanayo Kenneth Asogwa. Inspection
of chemical reaction and viscous dissipation on mhd convection flow over an
infinite vertical plate entrenched in porous medium with soret effect. Biomass
Conversion and Biorefinery, 14(6):7683–7694, 2024.
[135] S Abhilash Anand Kumar, S Sreedhar, and M Veera Krishna. Heat and mass
transfer on unsteady mhd convective flow through porous medium between two
vertical plates with chemical reaction. Proceedings of the Institution of Mechanical
Engineers, Part E: Journal of Process Mechanical Engineering, 238(4):1665–1675,
2024.
[136] M Harikrushna, G Prabhakara Rao, and A Mallikarjuna Reddy. Chemical reaction
effects on magnetohydrodynamic convective flow past a vertical absorbent plate
with ramped wall temperature and concentration. Heat Transfer, 52(6):3914–3935,
2023.
[137] K Sudarmozhi, D Iranian, and Ilyas Khan. Chemical reaction effect across the
moving flat plate with heat generation and mhd flow of maxwell fluid with viscous
dissipation. International Journal of Modern Physics B, page 2550027, 2024.
[138] Lakshmi Appidi, P Pramod Kumar, Sweta Matta, and Bala Siddulu Malga.
Effects of chemical reaction and thermal radiation on mhd free convective flow
of micro polar fluid past an infinite moving vertical porous plate with viscous
dissipation. Heat Transfer, 52(3):2922–2939, 2023.
[139] Purnima Rai and Upendra Mishra. Unveiling the behavior of mhd mixed con-
vective nanofluid slip flow over a moving vertical plate with radiation, chemical
reaction, and viscous dissipation. CFD Letters, 15(10):93–109, 2023.
[140] Thammisetty Ramanjaneyulu and Lingampalli Hari Krishna. Chemical reaction
and magnetohydrodynamics effects on heat absorbing fluid past an inclined
porous plate. Heat Transfer, 52(5):3595–3614, 2023.
[141] Ch Hk Gopal, M Sudhakar, S Masthan Rao, and M Veera Krishna. Chemical
reaction on unsteady mhd convection flow of second grade fluid over an un-
bounded perpendicular absorbent plate. Advances in Materials and Processing
Technologies, 10(3):2263–2285, 2024.
[142] B Chandra Sekhar and P Vijaya Kumar. Chemical reaction impacts the unsteady
mhd convective flow of an incompressible viscous fluid past an infinite vertical
porous plate. Heat Transfer, 52(2):1474–1496, 2023.
[143] B Prabhakar Reddy, MD Shamshuddin, SO Salawu, and L Joseph Sademaki. Com-
putational analysis of transient thermal diffusion and propagation of chemically
reactive magneto-nanofluid, brinkman-type flow past an oscillating absorbent
plate. Partial Differential Equations in Applied Mathematics, page 100761, 2024.
[144] B Prabhakar Reddy, PM Matao, and JM Sunzu. Hydromagnetic mixed convection
flow of heat consuming dissipative casson reacting and radiating fluid with
suction/blowing past a permeable oscillating plate in a porous medium: Fem
study. Special Topics & Reviews in Porous Media: An International Journal,
16:39–61, 2025.
References References
[145] Digbash Sahu and Rudra Kanta Deka. Dynamic interaction of stratifications
in a mhd parabolic flow with periodic temperature variation and variable mass
diffusion in a porous medium. International Communications in Heat and Mass
Transfer, 159:108203, 2024.
[146] Digbash Sahu and Rudra Kanta Deka. Influences of thermal and mass strat-
ification on unsteady magnetohydrodynamics parabolic flow along an infinite
vertical plate with periodic temperature variation and exponential mass diffusion
in porous medium. Heat Transfer, 53(8):4822–4850, 2024.
[147] Digbash Sahu and Rudra Kanta Deka. Influences of thermal stratification and
chemical reaction on mhd free convective flow along an accelerated vertical plate
with variable temperature and exponential mass diffusion in a porous medium.
Heat Transfer, 53(7):3643–3666, 2024.
[148] Digbash Sahu and Rudra Kanta Deka. Combined impacts of thermal and mass
stratification on unsteady mhd parabolic flow along an infinite vertical plate with
periodic temperature variation and variable mass diffusion. Heat Transfer, 2024.
[149] Digbash Sahu and Rudra Kanta Deka. Thermal stratification and chemical
reaction effects on mhd flow through oscillatory vertical plate in a porous medium
with temperature variation and exponential mass diffusion. East European Journal
of Physics, 2:209–218, 2024.
[150] Digbash Sahu and Rudra Kanta Deka. Numerical study of mhd ternary hybrid
nanofluid (ag cofe2o4-zno/c2h6o2+h2o): Effects of thermal stratification on a
vertical stretching cylinder in a porous environment. biointerface research in
applied chemistry, 14, 2024.
[151] Ilyas Khan. Significance of ramped temperature in the dynamics of unsteady
viscoelastic fluid subjected to lorentz force. Frontiers in Physics, 10:924910, 2022.
[152] Ilyas Khan. Prabhakar fractional derivative model of sodium alginate (c6h9nao7)
for accelerated plate motions. Frontiers in Energy Research, 10:1013829, 2022.
[153] Ilyas Khan. Hemodynamics of blood flow over an inclined cylinder. Waves in
Random and Complex Media, pages 1–12, 2022.
[154] Ilyas Khan. New idea of atangana and baleanu fractional derivatives to human
blood flow in nanofluids. Chaos: An Interdisciplinary Journal of Nonlinear
Science, 29(1), 2019.
[155] Ilyas Khan. Ramped heating in cnts fractional nanofluids. Case Studies in
Thermal Engineering, 45:102836, 2023.
[156] Ilyas Khan. Entropy generation and thermal analysis in quadratic translation
of sodium alginate with mhd and porosity effects. Journal of Thermal Analysis
and Calorimetry, pages 1–13, 2024.
References 173
[157] Ilyas Khan, Ali Raza, M Aslam Shakir, Amnah S Al-Johani, Amjad Ali Pasha,
and Kashif Irshad. Natural convection simulation of prabhakar-like fractional
maxwell fluid flowing on inclined plane with generalized thermal flux. Case
Studies in Thermal Engineering, 35:102042, 2022.
[158] Hari R Kataria and Akhil S Mittal. Velocity, mass and temperature analysis of
gravity-driven convection nanofluid flow past an oscillating vertical plate in the
presence of magnetic field in a porous medium. Applied Thermal Engineering,
110:864–874, 2017.
[159] Hari R Kataria and Akhil S Mittal. Mathematical model for velocity and
temperature of gravity-driven convective optically thick nanofluid flow past an
oscillating vertical plate in presence of magnetic field and radiation. Journal of
the Nigerian Mathematical Society, 34(3):303–317, 2015.
[160] M Sheikholeslami, Hari R Kataria, and Akhil S Mittal. Radiation effects on
heat transfer of three dimensional nanofluid flow considering thermal interfacial
resistance and micro mixing in suspensions. Chinese journal of physics, 55(6):2254–
2272, 2017.
[161] Zhixiong Li, M Sheikholeslami, Akhil S Mittal, Ahmad Shafee, and Rizwan-
ul Haq. Nanofluid heat transfer in a porous duct in the presence of lorentz
forces using the lattice boltzmann method. The European Physical Journal Plus,
134:1–10, 2019.
[162] Harshad R Patel, Akhil S Mittal, and Rakesh R Darji. Mhd flow of micropolar
nanofluid over a stretching/shrinking sheet considering radiation. International
Communications in Heat and Mass Transfer, 108:104322, 2019.
[163] Akhil S Mittal and Harshad R Patel. Influence of thermophoresis and brownian
motion on mixed convection two dimensional mhd casson fluid flow with non-
linear radiation and heat generation. Physica A: Statistical Mechanics and its
Applications, 537:122710, 2020.
[164] Hari R Kataria and Akhil S Mittal. Analysis of casson nanofluid flow in presence
of magnetic field and radiation. Mathematics Today, 33(1):99–120, 2017.
[165] Akhil S Mittal and Hari R Kataria. Three dimensional cuo–water nanofluid flow
considering brownian motion in presence of radiation. Karbala International
Journal of Modern Science, 4(3):275–286, 2018.
[166] HR Kataria and AS Mittal. Mathematical analysis of three dimensional nanofluid
flow in a rotating system considering thermal interfacial resistance and brownian
motion in suspensions through porous medium. Mathematics Today, 34(A):7–24,
2018.
[167] Akhil S Mittal, Harshad R Patel, and Rakesh R Darji. Mixed convection
micropolar ferrofluid flow with viscous dissipation, joule heating and convective
boundary conditions. International Communications in Heat and Mass Transfer,
108:104320, 2019.
References References
[168] Akhil S Mittal. Analysis of water-based composite mhd fluid flow using ham.
International Journal of Ambient Energy, 42(13):1538–1550, 2021.
[169] Akhil S Mittal. Study of radiation effects on unsteady 2d mhd al2o3-water flow
through parallel squeezing plates. International Journal of Ambient Energy,
43(1):653–660, 2022.
[170] Hari R Kataria, Mital Mistry, and Akhil Mittal. Influence of nonlinear radiation
on mhd micropolar fluid flow with viscous dissipation. Heat Transfer, 51(2):1449–
1467, 2022.
[171] Hari Kataria, Akhil S Mittal, and Mital Mistry. Effect of nonlinear radiation
on entropy optimised mhd fluid flow. International Journal of Ambient Energy,
43(1):6909–6918, 2022.
[172] Harshad Patel, Akhil Mittal, and Tejal Nagar. Fractional order simulation for
unsteady mhd nanofluid flow in porous medium with soret and heat generation
effects. Heat Transfer, 52(1):563–584, 2023.
[173] Akhil Mittal, Harshad Patel, Ramesh Patoliya, and Vimalkumar Gohil. Effects
of magnetic field and chemical reaction on a time dependent casson fluid flow.
Applications and Applied Mathematics: An International Journal (AAM), 19(3):8,
2024.
[174] Harshad Patel, Akhil Mittal, and Tejal Nagar. Effect of magnetic field on unsteady
mixed convection micropolar nanofluid flow in the presence of non-uniform heat
source/sink. International Journal of Ambient Energy, 45(1):2266748, 2024.
[175] R Muthucumaraswamy. Chemical reaction effects on vertical oscillating plate with
variable temperature. Chemical Industry and Chemical Engineering Quarterly,
16(2):167–173, 2010.
[176] Milton Abramowitz and Irene A Stegun. Handbook of mathematical functions
with formulas, graphs, and mathematical tables, volume 55. US Government
printing office, 1968.
[177] Fayza Mohammed Nasser El-Fayez. Effects of chemical reaction on the unsteady
free convection flow past an infinite vertical permeable moving plate with variable
temperature. Journal of surface engineered materials and advanced technology,
2(2):100–109, 2012.
[178] Rupam Shankar Nath and Rudra Kanta Deka. Thermal and mass stratification
effects on unsteady parabolic flow past an infinite vertical plate with expo-
nential decaying temperature and variable mass diffusion in porous medium.
ZAMM-Journal of Applied Mathematics and Mechanics/Zeitschrift für Ange-
wandte Mathematik und Mechanik, page e202300475, 2024.
[179] Rupam Shankar Nath and Rudra Kanta Deka. Thermal and mass stratification
effects on mhd nanofluid past an exponentially accelerated vertical plate through
a porous medium with thermal radiation and heat source. International Journal
of Modern Physics B, page 2550045, 2024.
References 175
[180] Rupam Shankar Nath and Rudra Kanta Deka. Theoretical study of thermal and
mass stratification effects on mhd nanofluid past an exponentially accelerated
vertical plate in a porous medium in presence of heat source, thermal radiation
and chemical reaction. International Journal of Applied and Computational
Mathematics, 10(2):92, 2024.
[181] Ahmed M Megahed and W Abbas. Non-newtonian cross fluid flow through
a porous medium with regard to the effect of chemical reaction and thermal
stratification phenomenon. Case Studies in Thermal Engineering, 29:101715,
2022.
[182] Dharmendar Reddy Yanala, M Anil Kumar, Shankar Goud Bejawada, Kot-
takkaran Sooppy Nisar, R Srinivasa Raju, and V Srinivasa Rao. Exploration of
heat and mass transfer on 3-d radiative mhd casson fluid flow over a stretching
permeable sheet with chemical reaction. Case Studies in Thermal Engineering,
51:103527, 2023.
[183] Nalivela Nagi Reddy, Dharmendar Reddy Yanala, B Shankar Goud, and Srini-
vasa Rao Vempati. Impact of porosity and radiation on two-dimensional unsteady
magnetohydrodynamics heat transfer stagnation point flow with viscous dissipa-
tion. Heat Transfer, 52(5):3538–3556, 2023.
[184] B Shankar Goud, Dharmendar Reddy Yanala, and Abderrahim Wakif. Numerical
analysis on the heat and mass transfer mhd flow characteristics of nanofluid
on an inclined spinning disk with heat absorption and chemical reaction. Heat
Transfer, 52(5):3615–3639, 2023.
[185] Yanala Dharmendar Reddy and Ippa Mangamma. Significance of radiation and
chemical reaction on mhd heat transfer nanofluid flow over a nonlinearly porous
stretching sheet with nonuniform heat source. Numerical Heat Transfer, Part A:
Applications, 85(18):2940–2966, 2024.
[186] Dharmendar Reddy Yanala, Shankar Goud Bejawada, and Kottakkaran Sooppy
Nisar. Influence of chemical reaction and heat generation/absorption on unsteady
magneto casson nanofluid flow past a non-linear stretching riga plate with
radiation. Case Studies in Thermal Engineering, 50:103494, 2023.
[187] Y Dharmendar Reddy and Ippa Mangamma. Numerical approach of fe3o4-
ethylene glycol heat and mass transfer magneto nanofluid flow past a porous
shrinking sheet with chemical reaction and thermal radiation. Journal of Thermal
Analysis and Calorimetry, 148(22):12639–12668, 2023.
[188] B Shankar Goud, Y Dharmendar Reddy, and ADNAN. Numerical investigation
of the dynamics of magnetized casson fluid flow over a permeable wedge subject to
dissipation and thermal radiations. Surface Review and Letters, 31(07):2450054,
2024.
[189] Yanala Dharmendar Reddy, Shankar Goud Bejawada, Nek Muhammad Katbar,
Ahmed M Zidan, Wasim Jamshed, Mohamed R Eid, and Fayza Abdel Aziz
ElSeabee. Heat generating impact on radiative nanofluid flow via exponential
References References
[200] Zafar Mahmood, Zahoor Iqbal, Maryam Ahmed Alyami, Bader Alqahtani,
Mansour F Yassen, and Umar Khan. Influence of suction and heat source
on mhd stagnation point flow of ternary hybrid nanofluid over convectively
heated stretching/shrinking cylinder. Advances in Mechanical Engineering,
14(9):16878132221126278, 2022.
[201] Wenhao Cao, IL Animasaun, Se-Jin Yook, VA Oladipupo, and Xianjun Ji.
Simulation of the dynamics of colloidal mixture of water with various nanoparti-
cles at different levels of partial slip: Ternary-hybrid nanofluid. International
Communications in Heat and Mass Transfer, 135:106069, 2022.
[202] Rupam Shankar Nath and Rudra Kanta Deka. A numerical investigation of the
mhd ternary hybrid nanofluid (cu-al2o3-tio2/h2o) past a vertically stretching
cylinder in a porous medium with thermal stratification. Journal of Advanced
Research in Fluid Mechanics and Thermal Sciences, 116(1):78–96, 2024.
[203] Rupam Shankar Nath and Rudra Kanta Deka. A numerical study on the mhd
ternary hybrid nanofluid (cu-al2o3-tio2/h2o) in presence of thermal stratification
and radiation across a vertically stretching cylinder in a porous medium. East
European Journal of Physics, 1:232–242, 2024.
[204] Ashish Paul, Neelav Sarma, and Bhagyashri Patgiri. Mixed convection of shear-
thinning hybrid nanofluid flow across a radiative unsteady cone with suction and
slip effect. Materials Today Communications, 37:107522, 2023.
[205] Ashish Paul, Neelav Sarma, and Bhagyashri Patgiri. Thermal and mass transfer
analysis of casson-maxwell hybrid nanofluids through an unsteady horizontal
cylinder with variable thermal conductivity and arrhenius activation energy.
Numerical Heat Transfer, Part A: Applications, pages 1–26, 2023.
[206] Lawrence F Shampine, Jacek Kierzenka, Mark W Reichelt, et al. Solving
boundary value problems for ordinary differential equations in matlab with bvp4c.
Tutorial notes, 2000:1–27, 2000.
[207] Jacek Kierzenka and Lawrence F Shampine. A bvp solver based on residual
control and the maltab pse. ACM Transactions on Mathematical Software
(TOMS), 27(3):299–316, 2001.
[208] NP Hale. A sixth-order extension to the matlab bvp4c software of j. kierzenka
and l. shampine. Department of Mathematics, Imperial College London, 2006.
[209] Iskandar Waini, Anuar Ishak, and Ioan Pop. Mixed convection flow over an expo-
nentially stretching/shrinking vertical surface in a hybrid nanofluid. Alexandria
Engineering Journal, 59(3):1881–1891, 2020.
[210] A Ishak and R Nazar. Laminar boundary layer flow along a stretching cylinder.
European Journal of Scientific Research, 36(1):22–29, 2009.
[211] EMA Elbashbeshy, TG Emam, MS El-Azab, and KM Abdelgaber. Laminar
boundary layer flow along a stretching cylinder embedded in a porous medium.
International Journal of Physical Sciences, 2012.
References References