Next Article in Journal
Fully Attentional Network for Low-Resource Academic Machine Translation and Post Editing
Previous Article in Journal
Deep Learning for Intelligent Human–Computer Interaction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temperature-Dependent Density and Magnetohydrodynamic Effects on Mixed Convective Heat Transfer along Magnetized Heated Plate in Thermally Stratified Medium Using Keller Box Simulation

1
Department of Mathematics and Statistics, Sargodha-Campus, The University of Lahore, Sargodha 40100, Pakistan
2
Rare Earth Elements Application and Research Center, Munzur University, 62000 Tunceli, Turkey
3
Department of Mathematics, Faculty of Science, University of Tabuk, P.O. Box 741, Tabuk 71491, Saudi Arabia
4
Center of Research, Faculty of Engineering, Future University in Egypt, New Cairo 11835, Egypt
5
Department of Mathematical Sciences, College of Science, Princess Nourahbint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
6
Department of Mathematics, Division of Science and Technology, University of Education, Lahore 54000, Pakistan
7
Department of Mathematics, Faculty of Science, University of Sargodha, Sargodha 40100, Pakistan
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2022, 12(22), 11461; https://doi.org/10.3390/app122211461
Submission received: 10 October 2022 / Revised: 9 November 2022 / Accepted: 10 November 2022 / Published: 11 November 2022
(This article belongs to the Topic Computational Fluid Dynamics (CFD) and Its Applications)

Abstract

:
The heat transmission properties along the non-magnetized geometries have been numerically obtainedby various researchers. These mechanisms are less interesting in engineering and industrial processes because of excessive heating. According to current studies, the surface is magnetized and the fluid is electrically conductive, which helps to lessen excessive surface heating. The main objective of the current analysis is to numerically compute the temperature-dependent density effect on magnetohydrodynamic convective heat-transfer phenomena of electrical-conductive fluid flow along the vertical magnetized and heated plate placed in a thermally stratified medium. For the purpose of numerical analysis, the theoretical process governing heat and magnetic intensity along a vertical magnetic plate is examined. By using suitable and well-known similarity transformations for integration, the non-linear coupled PDEs for the aforementioned electrical-conductive fluid flow mechanism are changed and subsequently converted into non-similar formulation. The Keller Box method is used to numerically integrate the final non-similar equations. The MATLAB software program plots the transformed algebraic equations graphically and quantitatively. The behavior of the physical quantities such asvelocity graph, magnetic field graph, and temperature plot along with their slopes that arerate of skin friction, the rate of heat transfer, and the rate of magnetic intensity for different parameters included in the flow model. The novelty of the current work is to compute the magneto-thermo analysis of electrically conducting flow along the vertical symmetric heated plate. First, we secure the numerical solution for steady part and then these results are used to find skin friction, heat transfer, and magnetic intensity. In the current work, the fluid becomes electrically conducing due to a magnetized surface which insulates heat during the mechanism and reduces the excessive heating. The results are excellent and accurate because they are satisfied by its given boundary conditions. Additionally, the current problems have a big impact on the production of polymer materials, glass fiber, petroleum, plastic films, polymer sheets, heat exchangers, catalytic reactors, and electronic devices.

1. Introduction and Literature Review

In metallurgy and polymer technology, the heat, magnetic intensity, and momentum transfer in the laminar boundary layer flow on a magnetized heated plate immersed in thermally stratified medium are significant from both a theoretical and practical standpoint. Engineers can benefit greatly from the laws governing heat transfer and magnetohydrodynamics flows in a variety of fields, including heat exchangers, sunspot theory, intercontinental ballistic missiles, interstellar gas motion, liquid metals, cooling of nuclear reactors, plasma confinement, and geophysics and astrophysics. The dynamic study of electrically conducting fluid is known as magnetohydrodynamics, established by Maxwell for electromagnetism. The fundamental idea behind an electrically conducting fluid is that a magnetic field has the potential to cause currents to flow through conductive fluid that is flowing, which in turn exerts forces on the fluid and affects the magnetic field. In literature, most researchers have analyzed the heat transfer characteristics of mixed convection flow along different geometries. However, in current analysis, the aligned magnetic field is used to make a surface magnetized and to reduce the excessive heating along the proposed surface. Some related literature of the current problem is mentioned here to understand the theoretical and numerical points. Furthermore, the current issues have significant implications for the polymer industries, including paper production, glass-fiber production, liquid crystal solidification, petroleum production, production of unusual lubricants, suspension solutions, wire drawing, continuous cooling, fiber spinning, plastic film production, polymer sheet extraction, heat exchangers, petroleum resource recovery, fault zones, catalytic reactors, and the production of electronic devices.
Nima and Ferdows [1] performed multiple solutions of microorganisms across a non-isothermal inclined disk with gyro-tactic particles. Mahanthesh et al. [2] developed an analysis on squeezing flow along a magnetized channel with hall current, heat source, and thermal radiation effects numerically. They found the increasing behavior in concentration and temperature by reducing value of squeezing parameter. Bansood and Ambedkar [3] studied free forced convective boundary layer flow over a semi-infinite plate in a porous medium. They analyzed heat and mass transfer in boundary layer flow region for opposing flow. Hazarika and Ahmed [4] investigated the numerical development of the MHD natural convective nanofluid flow through a vertical cone with heat source effects. They discovered adecreasing behavior of concentration profile when the Brownian motion parameter is raised. The laminar mixed convection flow from a vertical surface was studied by Chaudhry et al. [5]. The melting effects on mixed convection heat and mass transfer in a non-Newtonian nanofluid submerged in a porous medium were carried out by Ram Reddy and Kairi in [6]. They discovered that when the buoyancy ratio grows, the heat transfer coefficient decreases in the entire range of porous geometry. Hussain et al. [7] developed the unsteady convection flow of two phase nanofluid around a cone embedded in porous material numerically. They found that under viscous dissipation impact, the heat transfer rate decreased with maximum viscous heat generation.
The consequences of temperature-dependent density on a steady mixed Newtonian fluid convective flow near to a vertically permeable heated plate in a stratified material has great importance in industries and engineering. In the manufacturing sector, double diffuse free forced convection over vertical bodies is crucial for the development of dependable machinery, nuclear reactors, satellites, and spacecraft. Using convective boundary conditions and induced magnetic field, Mahdy [8] studied the nanofluid including gyro-tactic microorganisms with specified density motile microorganisms for natural convection flow across a vertical flat sheet. They found that as time dilation, heat conduction, and mixed convection parameter increases, the Nusselt number decreases. In order to simulate mixed convection flow of viscous fluid from a vertical cone with the effect of chemical reaction, Pullepu et al. [9] used numerical methods. They discovered that as skin friction increases, the velocity of steady flow is raised. Srinivasacharya and Surender [10] investigated impact of dual stratification on free forced convective boundary layer flow of nanofluid over a vertical shape embedded in a porous media. They obtained a considerable decrease in heat transmission and an increase in mass transfer rates with thermophoresis effect. Ullah et al. [11,12,13,14] elaborated impacts of MHD, velocity slip, reduced gravity, and stratification on oscillatory fluid along a non-conducting circular heated shape numerically. The baffle effects on enclosure filled with nanofluid for convective heat mechanism were analyzed in [15]. The entropy analysis on pseudo plastic flow along circular shape with MHD by using Keller Box scheme was investigated numerically in [16]. Sofiadis [17] numerically analyzed the Reynolds effects on micro polar flow in a channel.
Mehmood et al. [18] performed free forced convection flow in thermally and solutal stratified magneto nanofluid with non uniform heat source. They examined the volume fraction profile of lowered Brownian motion parameter and solute parameter increases. In the presence of internal heat generation, Elbashbeshy et al. [19] investigated the effects of thermal radiation and magnetic field on unsteady mixed convection flow and heat transfer across an exponentially extending surface with suction. They discovered that the radiative heat transfer increased as the thermal radiation parameter rose. Daniel et al. [20] developed a heat energy problem from a nanofluid flow across a stretched sheet with MHD and a chemical reaction. They observed that velocity and temperature profile rise byincrements in the electrical parameter.
The impact of a chemical reaction on the flow of MHD viscoelastic fluid over vertical stretching sheet with heat source effects was described by Jena et al. in [21]. They found that the chemical reaction and Prandtl number increases the rate of heat transfer along the plate. Yasin et al. [22] performed boundary layer flow with free forced convection embedded in thermally stratified porous medium statured by nanofluid numerically. Rawi et al. [23] introduced the influence of copper nanoparticles on Jeffrey fluid of free forced convection flow generated by g-Jitter. Alzahrani et al. [24] developed a third grade hydro magnetic free/forced convection problem for nonmaterial including gyro tactic microorganisms across a stretched surface. The effects of magnetohydrodynamics and velocity-slip on convective flow along the non-conducting geometry were studied by Ullah et al. in [25]. The opposing flow through a vertical cylinder filled with a Darcy porous medium containing gyro tactic microorganisms with double diffusive mixed convection was examined by Alsenafi and Ferdows in [26]. Babu et al. [27] examined the impact of thermal conductivity and variable viscosity on heat and mass transfer of nanofluid along a vertical cone in the presence of a chemical reaction. Variable density effects on hydro magnetic mixed convective non-Newtonian fluid flow over a moving vertical plate were presented by Salem et al. in [28]. Iqbal et al. [29] investigated the forced convective flow along a cone with heat and mass flux model by using implicit finite difference method graphically. Abulnaja et al. [30] performed unsteady free/forced convection flow filled with nanofluid across the straight disc placed in a porous material in the presence of variable surface temperature. Maranna et al. [31] analyzed the mass transpiration phenomena with MHD and radiation effects on viscous fluid past a stretched shape. Mabood et al. [32] performed Jeffrey nanofluid phenomena with entropy and MHD effects along the stretched geometry for various engineering application.
It can be seen that all of the aforementioned investigations deal with a minimum temperature difference between the fluid at the surface and the ambient fluid. However, there are instances where there is a significant temperature difference. The Boussinesq approximation is no longer appropriate in this case (see Kakac and Yenar [33]). To notice these circumstances, the temperature-dependent density function is suggested in [33] for the specific cases where alarge temperature difference occurs. The theoretical mechanism on heat and magnetic intensity along a vertical magnetic plate has been investigated for numerical analysis. Using the idea from previous research [33], the temperature-dependent density effect on magnetohydrodynamic convective heat-transfer phenomena of electrical-conductive fluid flow along the vertical magnetized and heated plate placed in thermally stratified medium was explored numerically. The non-similar equations werenumerically integrated by employing the Keller Box method. The transformed algebraic equations were plotted graphically and numerically on the MATLAB software package. The heat transfer characteristics along the non-magnetized shapes were performed in various previous studies numerically. Due to excessive heating, these mechanisms are less interesting in engineering and industrial processes. In the current analysis, the surface wasmagnetized and the fluid waselectrically conducting which is responsible forreducing excessive heating along the surface. The novelty of the present work is in the analysis of convective heat transfer analysis of viscous fluid flow with temperature dependent density and thermal stratification effects along the vertical symmetric heated plate numerically. First, the velocity, magnetic field, and temperature profiles weresecured and then those profiles used to compute skin friction, heat transfer, and magnetic intensity along the heated plate. The obtained results are in excellent agreement by satisfying the given boundary conditions.

2. Flow Problem and Mathematical Description

The problem is devoted tonumerical solutions of magneto-thermo analysis of convective heat transfer along thermally and magnetized geometry immersed in thermally stratified medium. The extended issue will change into arrangement of partial differential conditions and then transform into ordinary ones using the stream-function forms. The molded problem will be resolved by using FDM method conjunct with the Keller Box method. The numerical result for considered material properties will be depicted in graphical and tabular form.
Consider two dimensional steady and incompressible electrical conducting fluid. In Figure 1, where x and y are axes along and perpendicular to the vertical magnetized plate immersed in thermally stratified medium, respectively, u ,   H x and v ,   H y are the velocity-fields in x and y orders, respectively. The H 0 is the magnetic field strength, H x = H 0   indicates that the magnetic field is exact at the surface, T is the temperature, T is the free stream temperature, κ is the thermal conductivity of the fluid, and C p is the specific heat. The ρ is the fluid density, μ o is the megnatic permeability, ν = μ / ρ is the kinematic fluid viscosity, β * is the volumetric coefficient of thermal expansion, and gravity acceleration is g . The following is how the equations for continuity, momentum, magnetic force, and energy are expressed where high temperature difference arises with temperature dependent density function:
ρ = ρ exp [ β T ( T T ) ] , β T = 1 ρ ( ρ T ) p
It should be emphasized that the Boussinesq approximation holds true for fluids with β T of the magnitude of 10 3 when the temperature differential between the fluid and the surface is of the order of O ( 10 2 ) . Consequently, leaving out higher-order terms and expanding Equation (1a) into a Taylor series, the above equation becomes ρ = ρ [ 1 β T ( T w T ) ] .
( ρ u ) x + ( ρ v ) y = 0
ρ H 1 x + ρ H 2 y = 0
u u x + v u y = μ y ( 1 ρ u y ) + g ρ ( ρ ρ ) + μ 0 4 π ρ ( H 1 H 1 x + H 2 H 1 y )
u H 1 x + v H 1 y = ν m 2 H 1 y 2 + ( H 1 u x + H 2 u y )
u T x + v T y = κ C p y ( 1 ρ T y )
The boundary conditions of the present model are,
u = 0 ;   v = 0 ;   H 1 = H 0 ,   H 2 = 0 ;   T = T w a t y = 0   at   y = 0
U U ,   T T   ,   H 1 0   as   y .
The proper boundary conditions for the temperature with temperature-dependent density and the thermal stratification components are described above. The Reynolds number is R e x = U x ν , U being the free steam velocity, T w = T o + b x L is the variable temperature, and T = T o + c x L is ambient variable temperature of the thermally and magnetized heated plate, b ,   c are positive constants and T o is the temperature constant along the thermally and magnetizedheated plate.

3. The Similarity Transformation and Stream Function Formulation

The proper dimensionless similarity and stream function to convert PDEs into ODEs are given as:
u = 1 ρ ¯ ψ y , v = 1 ρ ¯ ψ x , H 1 = 1 ρ ¯ δ φ δ y , H 2 = 1 ρ ¯ δ φ δ x , θ = T T T w T .
The stream function ϕ and ψ are for magnetic and velocity fields, θ is dimensilized temperature and similarity parameter is η .
u = u f ,                         v = 1 2 ρ ¯ v u x ,
H 1 = g H 0 , H 2 = 1 2 ρ ¯ H 0 U ν U x ( g η g ) ,
η = 0 y ρ ¯ d y u ν x , ψ = ν x u f ( η ) , φ = H 0 U ν x u g
The converted form of the ODEs in Equations (1)–(6) by applying the Equations (7) and (8), we have non linear ODEs given below;
1 2 f f + e n θ f n θ e n θ f + λ ( 1 e n θ 1 e n ) + e n θ M g g = 0
P m ( e n θ g n θ e n θ g ) + f g 2 g f 2 = 0
e n θ θ n e n θ θ 2 + P r ( 1 2 f θ f θ S t f ) = 0
where P r = ν α is Prandtl parameter, P m = ν m ν is magnetic Prandtl parameter, M = H 0 2 L 2 μ 0 ν 2 4 π ρ is a magnetic parameter, λ = g β * T 0 U 2 is the mixed convection parameter, η is similarity variable, θ is dimensionless temperature, ν = μ ρ is kinematic fluid viscosity, and S t = c b is stratification number. The transformed boundary conditions in (6) then becomes,
f = 0 , f = 0 , g = 0 , θ = 1 S t , g = 1     at   η = 0
f 1 , θ 0 , g 0 ,   at   η .
The formula for skin friction is C f x = τ w ρ U 2 , the formula for the Nusselt number is N u x = x q w κ ( T w T ) , and the formula for the magnetic intensity coefficient is M g x = j w U 2 . The values of τ w , q w , and j w are listed below:
τ w = μ ( u y ) y = 0 , q w = κ ( T y ) y = 0 , j w = ν m ( H y ) y = 0
Consequently, the values of magnetic intensity, Nusselt number, and skin friction are provided by,
R e x 1 / 2 C f x = f ( 0 ) , R e x 1 / 2 N u x = θ ( 0 ) , R e x 1 / 2 M g x = g ( 0 ) .

4. The Solution and Computing Techniques

A suitable stream function formulation is used to convert the linked partial differential equations for the aforementioned model into a set of ordinary differential equations. The iterative Keller Box scheme is used to solve the simplified nonlinear ordinary differential Equations (9)–(11) and boundary conditions in (12). We begin by introducing additional independent variables p ( η ) ;   q ( η )   ,   u ( η ) , and v ( η ) by using Equation (13),
f = p ,   f = p = q ,   f = q , g = u ,
g = u = v   , g = v   , θ = l   , θ = l = m .
Equations (9)–(12) get simpler to solve in order to get around the problems that they present. Thus, the more simple forms are,
f = p f p = 0
p = q p q = 0
g = u g u = 0
u = v u v = 0
θ = l θ l = 0
f q 2 + ( 1 n θ ) q n l ( 1 n θ ) q + λ θ + ( 1 + n θ ) M g v = 0
( 1 n θ ) [ P m v n l v ] + f v 2 g q 2 = 0
( 1 n θ ) m n ( 1 n θ ) l 2 + P r ( f l 2 p θ S t p ) = 0
The reduced boundary conditions are
f = 0 ,   g = 0 ,   p = 0 ,   u = 1 , θ = 1 S t           a t           η = 0 ,
p 1 , u 0 , θ 0 , a s             η .
We now consider the segment η n 1 , η n with η n 1 2 as midpoint given in Equation (23) below;
η 0 = 0 ,           η n = η n 1 + h n ,                 η n = η .
This is how the central difference form and average form are presented:
f = f n f n 1   h n ,                       f = f n + f n 1   2 = f n 1 2 ,
and
f n f n 1   1 2 h n ( p n + p n 1   ) = 0 ,
p n p n 1   1 2 h n ( q n + q n 1   ) = 0 ,
g n g n 1 1 2 h n ( u n + u n 1 ) = 0 ,
u n u n 1 1 2 h n ( v n + v n 1   ) = 0 ,
θ n θ n 1   1 2 h n ( l n + l n 1 ) = 0 ,
By applying the above Equations (24)–(29), the governing Equations (14)–(21) become
( 1 8 ( f n + f n 1 ) ) ( q n + q n 1   ) ) + ( q n q n 1 h n ) + n [ 3 4 ( q n + q n 1 ) ( l n + l n 1 ) + 1 2 ( p n + p n 1 ) + 1 16 ( f n + f n 1 ) ( p n + p n 1 ) ( l n + l n 1 ) ] + n 2 [ 1 4 ( p n + p n 1 ) ( l n + l n 1 ) + 1 8 ( f n + f n 1 ) ( l n + l n 1 ) 2 ] + ( 1 n ) ( θ n + θ n 1   ) [ λ ( θ n + θ n 1 ) M 4 ( g n + g n 1 ) ( v n + v n 1 ) ] = 0
1 n 2 ( θ n + θ n 1   ) [ P m h n ( v n v n 1 ) n 4 ( l n + l n 1 ) ( v n + v n 1 ) ] + 1 8 ( f n + f n 1 ) ( v n + v n 1 ) 1 8 ( g n + g n 1 ) ( q n + q n 1 ) = 0
1 n 2 ( θ n + θ n 1   ) [ m n + m n 1 2 n ( l n + l n 1 2 ) 2 ] + P r [ 1 8 ( f n + f n 1 ) ( l n + l n 1 ) 1 4 ( p n + p n 1 ) ( θ n + θ n 1 ) S t 2 ( p n + p n 1 ) ] = 0
along with boundary conditions
f 0 = 0 ,   g 0 = 0 , P 0 = 0 , u 0 = 1 , θ 0 = 1 S t , a t           η = 0
p n 1 , θ n 0 , u n 0 , a s           η
Applying the below-mentioned iterative Newton–Raphson approach for smooth algorithm;
f n k + 1   = f n k + δ f n k , p n k + 1   = p n k + δ p n k
q n k + 1   = q n k + δ q n k , θ n k + 1   = θ n k + δ θ n k
u n k + 1 = u n k + δ u n k , g n k + 1 = g n k + δ g n k
v n k + 1   = v n k + δ v n k , l n k + 1 = l n k + δ l n k
By ignoring all appearances of powers of greater than the first power, just like in the standard Newton–Raphson method, equations become;
δ f n δ f n 1 1 2 h n ( δ p n + δ p n 1 ) = ( r 1 ) n
δ p n δ p n 1 1 2 h n ( δ q n + δ q n 1 ) = ( r 2 ) n
δ g n δ g n 1 1 2 h n ( δ v n + δ v n 1 ) = ( r 3 ) n
δ u n δ u n 1 1 2 h n ( δ v n + δ v n 1 ) = ( r 4 ) n
δ θ n δ θ n 1 1 2 h n ( δ l n + δ l n 1 ) = ( r 5 ) n
The simplified form of the equations is provided below, again utilizing Equations (35)–(39) in Equations (30)–(33);
( a 1 ) n δ f n + ( a 2 ) n δ f n 1 + ( a 3 ) n δ q n + ( a 4 ) n δ q n 1 + ( a 5 ) n δ l n + ( a 6 ) n δ l n 1 + ( a 7 ) n δ θ n + ( a 8 ) n δ θ n 1 + ( a 9 ) n δ g n + ( a 10 ) n δ g n 1 = ( r 6 ) n
( b 1 ) n δ v n + ( b 2 ) n δ v n 1 + ( b 3 ) n δ θ n + ( b 4 ) n δ θ n 1 + ( b 5 ) n δ l n + ( b 6 ) n δ l n 1 + ( b 7 ) n δ g n + ( b 8 ) n δ g n 1 + ( b 9 ) n δ v n + ( b 10 ) n δ v n 1 = ( r 7 ) n
( c 1 ) n δ m n + ( c 2 ) n δ m n 1 + ( c 3 ) n δ θ n + ( c 4 ) n δ θ n 1 + ( c 5 ) n δ l n + ( c 6 ) n δ l n 1 + ( c 7 ) n δ p n + ( c 8 ) n δ p n 1 + ( c 3 ) n δ f n + ( c 4 ) n δ f n 1 = ( r 8 ) n
The specific boundary conditions that can be met without iteration are recalled. As a result, to ensure that these correct values are maintained across all iterations, we take
δ f     0 = 0 , δ g 0 = 0 , δ p     0 = 0 , δ u 0 = 1 , δ θ     0 = 1 S t
δ p     n = 1 , δ θ     n = 0 , δ u n = 0

5. The Arrangement of Difference Equations in Vector Notation/Matrix

The next essential step is to organize the aforementioned difference equations in matrix form. If it is done wrong, either the matrix solution method fails because there is a singular matrix (determinant = 0) or sub-matrix, or the approach becomes extremely inefficient since there is no discernible structure in the matrix. We have,
A δ = r ,
[ A ] = [ [ A 1 ] [ C 1 ] [ B 2 ] [ A 2 ] [ C 2 ] [ B n 1 ] [ A n 1 ] [ C n 1 ] [ B n ] [ A n ] ] , [ δ ] = [ [ δ 1 ] [ δ 2 ] [ δ n 1 ] [ δ n ] ] , [ r ] = [ [ r 1 ] [ r 2 ] [ r n 1 ] [ r n ] ] ,

6. Analysis and Discussion of Results

The present analysis addresses the temperature-dependent density effect on magnetohydrodynamic convective heat transfer phenomena of electrical-conductive fluid flow along the vertical magnetized and heated plate placed in thermally stratified medium. The non-linear coupled PDE equations for the above fluid flow mechanism are formulated and then converted into non-similar formulation by applying an appropriate and well-known similarity transformation for integration. The final non-similar forms are integrated numerically by employing Keller Box scheme method. The transformed algebraic equations are plotted graphically and numerically on the MATLAB software package. The behavior of physical quantities such as velocity graph, magnetic field graph, and temperature graph along with their slopes, that is skin friction, magnetic intensity, and heat transfer, under the effect of different parameters included in the flow model is discussed. The impact of physical parameters such as Prandtl number P r , mixed convective parameter λ , stratification parameter S t , density number m , magnetic-force number M , and magnetic Prandtl P m is interpreted graphically and numerically.
Figure 2a–c are presented for the behavior of velocity graph, temperature plot, and magnetic field for different values of temperature dependent density m = 0.0 ,   1.0 ,   2.0 ,   and   3.0 with some fixed parameters. It was deduced that the velocity profile was enhanced at maximum value of m = 3.0 but similar behavior wasnoted for other values of m in Figure 2a. This happens because an increasing value of density parameter tends to increase the density variation with temperature, which ultimately increases the buoyancy force. The magnetic field wasincreased at the lower value of m = 0.0 but decreased at the maximum value of m = 3.0 in Figure 3b. In Figure 3c, the temperature distribution wasincreased at lower value of m = 0.0 but decreased at maximum value of m = 3.0 . The good variations in each profile werenoted for each value of temperature-density parameter m . Hence, as m increased, the fluid was under two forces: the first force increased the velocity of the fluid due to increase in the buoyancy forces and the second force decreased the velocity of the fluid due to the decrease in temperature. Figure 3a–c are plotted against the various values of stratification parameter S t = 0.0 ,   0.1 ,   0.2 ,   0.3 for velocity graph, temperature graph, and magnetic field. It was examined that the velocity profile increased at lower value of S t = 0.0 but the minimum value wasnoted at the higher value of S t = 0.3 in Figure 3a. It wasdeduced that the magnetic field wasenhanced at the large value of S t = 0.3 but decreased at the small value of S t = 0.0 in Figure 3b. This happens because magnetic diffusion is reduced by increasing the magnetic-Prandtl number which is responsible for above-said phenomena. This behavior of the aforementioned variables supports the physical hypothesis that as the magnetic field grows stronger, more resistance is formed inside the fluid flow domain, which reduces the velocity distribution and improves the temperature profile. The temperature distribution increased at the lower value of S t = 0.0 but decreased at the maximum value of S t = 3.0 in Figure 3c. Figure 4a–c depict the behavior of velocity f ( η ) , temperature θ ( η ) , and magnetic field g ( η ) for different quantities of λ = 0.1 ,   0.5 ,   1.5 ,   and   2.0 with some fixed parameters. It was deduced that the velocity profile wasenhanced at larger λ = 2.0 but decreased at small value of λ = 0.1 in Figure 4a. It was examined that the magnetic field increased at lower value of λ = 0.1 while minimum value wasnoted at higher value of λ = 2.0 in Figure 4b. It wasalso noted that the temperature distribution increased at lower value of λ = 0.1 but decreased at maximum value of λ = 2.0 in Figure 4c. Physically, it was expected because larger values of λ correspond to stronger buoyancy forces which lead to an increase in the acceleration of fluid flow. Due to conducting phenomena, magnetic effects are strongly observed exactly at the surface, but far from the surface, each value is zero. The good variations in each profile werenoted for each value of mix convection parameter λ . Figure 5a–c examined the behavior of velocity f ( η ) , temperature θ ( η ) , and magnetic field g ( η ) for different quantities of Prandtl number P r = 1.0 ,   4.0 ,   7.0 ,   10.0 with some fixed parameters. It was shown that the velocity field increased at the highest value of P r = 10.0 and decreased at P r = 1.0 in Figure 5a. The magnetic field distribution increased at the lowest value of P r = 1.0 but minimum value wasnoted at the highest value of P r = 10.0 in Figure 5b. In Figure 5c, the temperature distribution wasincreased by decreasing the value of P r = 1.0 but the temperature distribution wasdecreased by increasing the value of P r = 10.0 . This phenomenon is expected because an increase in Pr tends to increase density variation with temperature which enhances the buoyancy force.The thermal boundary layer thickness wasreduced due to increase in Pr. The prominent thermal slip with good variations are obtained in temperature graph at each value of P r . The good variations in each profile wereexamined for various values of Prandtl number P r . Figure 6a–c are demonstrated for velocity f ( η ) , temperature θ ( η ) , and magnetic field g ( η ) for different quantities of magnetic Prandtl number P m = 0.1 ,   0.7 ,   1.5 ,   3.0 with some fixed parameters along the surface of a vertical geometry. In Figure 6a, it wasdeduced that the fluid velocity along the surface wasincreased for the large value of P m = 3.0 but the minimum value wasobtained at a lower value of P m = 1.0 with good response. It was deduced that the magnetic graph increased at the large value of P m = 3.0 but decreased at the small value of P m = 1.0 in Figure 6b. In Figure 6c, it is depicted that the temperature distribution increased at the small value of P m = 1.0 but the minimum value wasnoted at the large value of P m = 3.0 . It wasalso noted that temperature plots are in good form and excellent agreement for the given phenomenon. As the magnetic Prandtl number increased, viscosity of fluid increased and fluid became thicker more and consequently the boundary layer thickness decreased. Figure 7a–c are plotted for velocity, magnetic field, and fluid temperature against the various values of magnetic force parameter M = 0.1 ,   5.0 ,   10.0 ,   and   20.0 with some constant parameters long the surface of the vertical geometry. In Figure 7a, it is deduced that the mixed convection effect in the velocity graph increased at the large value of M = 20.0 but minimum value wasobtained at the lower value of M = 0.1 by satisfying the given boundary conditions asymptotically. This result was expected because an increase in M means enhancement in the Lorentz forces which opposes the flow and velocity of the fluid decreases. In Figure 7b, it wasexamined that the magnetic field profile along the surface wasenhanced at the large value of M = 20.0 but the minimum value was obtained at small value of M = 0.1 . Thishappens because the magnetic force parameter is the ratio of magnetic energy to kinetic energy, so with the increase in the magnetic force parameter P m , the magnetic energy is increased, while the kinetic energy is reduced. An excellent and favorable phenomenon is observed in magnetic field graph. In Figure 7c, it is deduced that the temperature profile decreased at the large value of M = 20.0 but increased at the small value of M = 0.1 . This given mechanism is performed with good agreement.
From Table 1, it is concluded that the skin friction f ( 0 ) increased at m = 0.5 but the minimum skin friction wasobtained at small m = 0.0 in the presence of buoyancy forces. It should be noted that m = 0 corresponds to the case of a small temperature difference between the surface and fluid, i.e., the buoyancy approximation becomes valid for this particular case, whereas for a high temperature difference, a considerable enhancement in the numerical values of the skin friction and heat transfer rate wasinspected and can be clearly seen in Table 1. It wasdeduced that the magnetic intensity g ( 0 ) wasmaximum at larger m = 0.5 but the smaller quantity of magnetic intensity wasexamined at smaller m = 0.0 in the presence of strong magnetic field. This result was expected because an increase in M means enhancement in the Lorentz forces which opposes the flow and velocity of the fluid decreases. It was also expected because as magnetic Prandtl number increases, the viscosity of fluid increases and fluid becomes thicker and consequently the boundary layer thickness decreases. It is also depicted that the heat transfer wasincreased at maximum m = 0.5 but the lower quantity of heat transfer is depicted at minimum m = 0.0 in the presence of magnetic force M . With the increase in value of M , the magnetic field became stronger along the surface which is a clear indication that the Lorentz force is more effective in this case. In Table 2, the influence of stratification S t number with various values along a vertical plate is shown to check the behavior of physical properties of the rate of skin friction, the rate of magnetic intensity, and the rate of heat transfer. From Table 2, it is computed that the skin friction is maximum at smaller value of S t = 0.1 but the minimum quantity of skin friction is depicted at a large value of S t = 0.4 in the presence of Prandtl number P r = 7.0 . Physically, it is possible because thermal conductivity decreases as Pr increases, which has a lower magnitude of frictional forces between the viscous layers. It is depicted that the magnetic intensity is decreased at larger S t = 0.4 but the maximum value of magnetic intensity is explored at small S t = 0.1 in the presence of λ = 6.5 . However, due to the strong buoyancy number λ which acts like a pressure gradient and dominated over the resistance, skin friction is increased and slight changes in heatand magnetic intensityare noted. It is also mentioned that the heat transfer θ ( 0 ) wasmaximum at the lower value of S t = 0.1 but the minimum quantity of heat transfer wasnoticed at larger S t = 0.4 in the presence of magnetic Prandtl number P m = 8.3 . The reason behind this is that increasing the value of P m is equivalent to decreasing magnetic diffusivity and consequently the strength of magnetic field becomes loose.

7. Conclusions

The current physical problem addressed the temperature-dependent density effect on magnetohydrodynamic convective heat-transfer phenomena of electrical-conductive fluid flow along the vertical magnetized and heated plate placed in thermally stratified medium. In previous studies, the heat transfer characteristics along the non-magnetized shapes were performed numerically. Due excessive heating, these mechanisms are less interesting in engineering and industrial processes. To overcome this issue in current analysis, the surface wasmagnetized and fluid waselectrically conducting which is responsible forreducing excessive heating along the surface. The non-linear coupled PDE equations for the above fluid flow mechanism wereformulated and then converted into non-similar formulation by applying an appropriate and well-known similarity transformation for integration. The final non-similar forms werenumerically integrated by employing the Keller Box scheme. The transformed algebraic equations wereplotted graphically and numerically on the MATLAB software package. The main finding of the current analysis wasto compute physical quantities such asvelocity graph, magnetic field graph, and temperature plot along with their slopes that areskin friction, magnetic intensity and heat transfer for different parameters included in the flow model. The comparison of Keller Box analysis wascomputed with previous existing numerical results. Furthermore, the current issues have significant implications for the polymer industries, glass-fiber production, petroleum production, fiber spinning, plastic film production, polymer sheet extraction, heat exchangers, catalytic reactors, and the production of electronic devices. The main findings are given below:
It wasdeduced that the velocity profile increased for maximum temperature-dependent density but similar behavior wasnoted for other values of variable density. This happens because an increasing value of density parameter tends to increase the density variation with temperature, which ultimately increases the buoyancy force.
It wasalso noted that the temperature distribution increased for strong buoyancy forces but decreased for lower buoyancy forces.Physically, it was expected because larger values of λ correspond to stronger buoyancy forces which leads to increase the acceleration of fluid flow.
It wasexamined that the magnetic field profile along the surface wasenhanced for strong magnetic forces but minimum value wasobtained at lower magnetic forces. Thishappens because the magnetic force parameter is the ratio of magnetic energy to kinetic energy, so with the increase in the magnetic force parameter P m , the magnetic energy is increased, while the kinetic energy is reduced.
It wascomputed that the skin friction wasmaximum for small thermal stratification but the minimum quantity of skin friction wasdepicted for large stratification in the presence of Prandtl number Pr = 7.0 . Physically, it is possible because thermal conductivity decreases as Pr increases, which has a lower magnitude of frictional forces between the viscous layers.
It wasalso mentioned that the heat transfer is increased for lower stratification, but the minimum quantity of heat transfer is noticed in the presence of magnetic Prandtl number. The reason behind this is that increasing the value of Pm is equivalent to decreasing magnetic diffusivity and consequently the strength of magnetic field becomes loose.

Author Contributions

Conceptualization, Z.U.; Data curation, N.A.; Formal analysis, M.A.; Funding acquisition, N.A.; Methodology, M.J.; Project administration, H.F.A.; Funding acquisition, A.M.A.; Visualization, A.H.; Writing—original draft, Z.U.; Funding acquisition, A.H.; Software, S.M.E.; Writing—review and editing, Z.U.; Funding acquisition, H.F.A.; Validation, A.M.A.; Funding acquisition, S.M.E.; Investigation, M.J.; Methodology, M.A. All authors have contributed in the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Princess Nourahbint Abdulrahman University researchers supporting project number (PNURSP2022R52) by Aisha M. Alqahtani, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Princess Nourahbint Abdulrahman University researchers supporting project number (PNURSP2022R52), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nima, N.I.; Ferdows, M. Dual Solutions in Mixed Convection Flow Along Non-Isothermal Inclined Cylinder Containing Gyrotactic Microorganism. J. Adv. Res. Fluid Mech. Therm. Sci. 2021, 87, 51–63. [Google Scholar] [CrossRef]
  2. Mahanthesh, B.; Gireesha, B.J.; Gorla, R.S.R. Mixed convection squeezing three-dimensional flow in a rotating channel filled with nanofluid. Int. J. Numer. Methods Heat Fluid Flow 2016, 26, 1460–1485. [Google Scholar] [CrossRef]
  3. Bansod, V.J.; Ambedkar, B. An integral treatment for combined heat and mass transfer by mixed convection along vertical surface in a saturated porous medium. J. Mech. Eng. Res. 2011, 3, 138–151. [Google Scholar]
  4. Hazarika, S.; Ahmed, S. Numerical Investigation of Endothermic/Exothermic Reaction in MHD Natural Convective Nano-Fluid Flow over a Vertical Cone with Heat Source/Sink. Walailak J. Sci. Technol. 2021, 18, 22834. [Google Scholar] [CrossRef]
  5. Chaudhary, K.; Sharma, A.; Jha, A. Laminar mixed convection flow from a vertical surface with induced magnetic field and convective boundary. Int. J. Appl. Mech. Eng. 2018, 23, 307–326. [Google Scholar] [CrossRef] [Green Version]
  6. Ramreddy, C.; Kairi, R. The effect of melting on mixed convection heat and mass transfer in non-Newtonian nanofluid SATURATEED in porous medium. Front. Heat Mass Transf. 2015, 6, 6. [Google Scholar] [CrossRef] [Green Version]
  7. Hussain, A.; Haider, Q.; Rehman, A.; Abdussattar, A.; Malik, M.Y. A New Heat Dissipation Model and Convective Two-Phase Nanofluid in Brittle Medium Flow over a Cone. Math. Probl. Eng. 2021, 2021, 6688747. [Google Scholar] [CrossRef]
  8. Mahdy, A. Free convection flow over a vertical flat plate in nanofluid porous media containing gyrotactic microorganisms with prescribed density motile microorganisms flux. Sci. Technol. Asia 2017, 22, 71–85. [Google Scholar]
  9. Pullepu, B.; Sambath, P.; Rani, M.S.; Chamkha, A.J.; Viswanathan, K.K. Numerical Solutions of Free Convective Flow from a Vertical Cone with Mass Transfer under the Influence of Chemical Reaction and Heat Generation/Absorption in the Presence of UWT/UWC. J. Appl. Fluid Mech. 2016, 9, 343–356. [Google Scholar] [CrossRef]
  10. Srinivasacharya, D.; Surender, O. Effect of double stratification on mixed convection boundary layer flow of a nanofluid past a vertical plate in a porous medium. Appl. Nanosci. 2014, 5, 29–38. [Google Scholar] [CrossRef] [Green Version]
  11. Ullah, Z.; Ashraf, M.; Rashad, A.M. Magneto-thermo analysis of oscillatory flow around a non-conducting horizontal circular cylinder. J. Therm. Anal. 2020, 142, 1567–1578. [Google Scholar] [CrossRef]
  12. Ullah, Z.; Ashraf, M.; Zia, S.; Ali, I. Surface temperature and freestream velocity oscillation effects on mixed convention slip flow from surface of a horizontal circular cylinder. Therm. Sci. 2020, 24 (Suppl. 1), 13–23. [Google Scholar] [CrossRef]
  13. Ullah, Z.; Ashraf, M.; Zia, S.; Chu, Y.-M.; Khan, I.; Nisar, K.S. Computational analysis of the oscillatory mixed convection flow along a horizontal circular cylinder in thermally stratified medium. Comput. Mater. Contin. 2020, 65, 109–123. [Google Scholar] [CrossRef]
  14. Ullah, Z.; Ashraf, M.; Sarris, I.E.; Karakasidis, T.E. The Impact of Reduced Gravity on Oscillatory Mixed Convective Heat Transfer around a Non-Conducting Heated Circular Cylinder. Appl. Sci. 2022, 12, 5081. [Google Scholar] [CrossRef]
  15. Al-Farhany, K.; Al-Muhja, B.; Loganathan, K.; Periyasamy, U.; Ali, F.; Sarris, I.E. Analysis of Convection Phenomenon in Enclosure Utilizing Nanofluids with Baffle Effects. Energies 2022, 15, 6615. [Google Scholar] [CrossRef]
  16. Al-Mdallal, Q.; Prasad, V.R.; Basha, H.T.; Sarris, I.; Akkurt, N. Keller box simulation of magnetic pseudoplastic nano-polymer coating flow over a circular cylinder with entropy optimisation. Comput. Math. Appl. 2022, 118, 132–158. [Google Scholar] [CrossRef]
  17. Sofiadis, G.; Sarris, I. Reynolds number effect of the turbulent micropolar channel flow. Phys. Fluids 2022, 34, 075126. [Google Scholar] [CrossRef]
  18. Mehmood, K.; Hussain, S.; Sagheer, M. Mixed convection flow with non-uniform heat source/sink in a doubly stratified magnetonanofluid. AIP Adv. 2016, 6, 065126. [Google Scholar] [CrossRef] [Green Version]
  19. Elbashbeshy, E.M.; Emam, T.; Abdelgaber, K. Effects of thermal radiation and magnetic field on unsteady mixed convection flow and heat transfer over an exponentially stretching surface with suction in the presence of internal heat generation/absorption. J. Egypt. Math. Soc. 2012, 20, 215–222. [Google Scholar] [CrossRef] [Green Version]
  20. Daniel, Y.S.; Aziz, Z.A.; Ismail, Z.; Salah, F. Thermal radiation on unsteady electrical MHD flow of nanofluid over stretching sheet with chemical reaction. J. King Saud Univ. -Sci. 2019, 31, 804–812. [Google Scholar] [CrossRef]
  21. Jena, S.; Dash, G.; Mishra, S. Chemical reaction effect on MHD viscoelastic fluid flow over a vertical stretching sheet with heat source/sink. Ain Shams Eng. J. 2016, 9, 1205–1213. [Google Scholar] [CrossRef] [Green Version]
  22. Yasin, M.H.M.; Arifin, N.M.; Nazar, R.; Ismail, F.; Pop, I. Mixed convection boundary layer flow embedded in a thermally stratified porous medium saturated by a nanofluid. Adv. Mech. Eng. 2013, 121943. [Google Scholar] [CrossRef]
  23. Rawi, N.A.; Ilias, M.R.; Jiann, L.Y.; Isa, Z.M.; Shafie, S. The effect of copper nanoparticles on mixed convection flow of jeffrey fluid induced by g-jitter. J. Nanofluids 2018, 7, 156–162. [Google Scholar] [CrossRef]
  24. Alzahrani, E.O.; Shah, Z.; Dawar, A.; Malebary, S.J. Hydromagnetic mixed convective third grade nanomaterial containing gyrotactic microorganisms toward a horizontal stretched surface. Alex. Eng. J. 2019, 58, 1421–1429. [Google Scholar] [CrossRef]
  25. Ullah, Z.; Ehsan, M.; Ahmad, H.; Ilyas, A. Combined effects of MHD and slip velocity on oscillatory mixed convective flow around a non-conducting circular cylinder embedded in a porous medium. Case Stud. Therm. Eng. 2022, 38, 102341. [Google Scholar] [CrossRef]
  26. Alsenafi, A.; Ferdows, M. Dual solution for double-diffusive mixed convection opposing flow through a vertical cylinder saturated in a Darcy porous media containing gyrotactic microorganisms. Sci. Rep. 2021, 11, 19918. [Google Scholar] [CrossRef] [PubMed]
  27. Babu, S.R.C.; Venkateswarlu, S.; Lakshmi, K.J. Effect of variable viscosity and thermal conductivity on heat and mass transfer flow of nanofluid over a vertical cone with chemical reaction. Int. J. Appl. Eng. Res. 2018, 13, 13989–14002. [Google Scholar]
  28. Salem, A.M.; El-Aziz, M.A.; Abo-Eldahab, E.M.; Abd-Elfatah, I. Effect of variable density on hydromagnetic mixed convection flow of a non-Newtonian fluid past a moving vertical plate. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 1485–1493. [Google Scholar] [CrossRef]
  29. Iqbal, M.S.; Khan, W.; Mustafa, I.; Ghaffari, A. Numerical study of natural convection flow of nanofluid past a circular cone with Cattaneo–Christov heat and mass flux models. Symmetry 2019, 11, 1363. [Google Scholar] [CrossRef] [Green Version]
  30. Abualnaja, K.M.; Elgendy, M.S.; Ibrahim, F.S. Unsteady Mixed Convection Flow along Symmetric Wedge with Variable Surface Temperature Embedded in a Porous Medium Saturated with a Nanofluid. J. Appl. Math. Phys. 2021, 09, 101–126. [Google Scholar] [CrossRef]
  31. Maranna, T.; Sneha, K.N.; Mahabaleshwar, U.S.; Sarris, I.E.; Karakasidis, T.E. An Effect of Radiation and MHD Newtonian Fluid over a Stretching/Shrinking Sheet with CNTs and Mass Transpiration. Appl. Sci. 2022, 12, 5466. [Google Scholar] [CrossRef]
  32. Mabood, F.; Fatunmbi, E.O.; Benos, L.; Sarris, I.E. Entropy Generation in the Magnetohydrodynamic Jeffrey Nanofluid Flow Over a Stretching Sheet with Wide Range of Engineering Application Parameters. Int. J. Appl. Comput. Math. 2022, 8, 98. [Google Scholar] [CrossRef]
  33. Kakac, S.; Yener, Y.; Pramuanjaroenkij, A. Convective Heat Transfer; CRC Press: Boca Raton, FL, USA, 1995; pp. 312–314, 399. [Google Scholar]
Figure 1. Coordinate system and flow geometry.
Figure 1. Coordinate system and flow geometry.
Applsci 12 11461 g001
Figure 2. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against m .
Figure 2. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against m .
Applsci 12 11461 g002
Figure 3. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against S t .
Figure 3. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against S t .
Applsci 12 11461 g003
Figure 4. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against λ .
Figure 4. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against λ .
Applsci 12 11461 g004
Figure 5. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against P r .
Figure 5. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against P r .
Applsci 12 11461 g005
Figure 6. (ac) The geometric profiles of f ( η ) ,     g ( η ) and θ ( η ) against P m .
Figure 6. (ac) The geometric profiles of f ( η ) ,     g ( η ) and θ ( η ) against P m .
Applsci 12 11461 g006
Figure 7. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against M .
Figure 7. (ac) The geometric profiles of f ( η ) ,     g ( η ) , and θ ( η ) against M .
Applsci 12 11461 g007
Table 1. Numerical results for f ( 0 ) ,   g ( 0 ) and for θ ( 0 ) for variousvalues of m = 0.0 ,   0.1 ,   0.3 ,   0.5 , while other parameters are fixed.
Table 1. Numerical results for f ( 0 ) ,   g ( 0 ) and for θ ( 0 ) for variousvalues of m = 0.0 ,   0.1 ,   0.3 ,   0.5 , while other parameters are fixed.
m f ( 0 ) g ( 0 ) θ ( 0 )
0.0 0.6956310646362440.8003916595769380.330556513906720
0.1 0.7284428705608600.8432079660300590.352119633819730
0.3 0.8091060156522430.9394448973139970.400936354502018
0.5 0.9134579742535701.0518444908519580.458533654900388
Table 2. Numerical results for f ( 0 ) ,   g ( 0 ) and for θ ( 0 ) for various values of S t = 0.1 ,   0.2 ,   0.3 ,   0.4 , while other parameters are fixed.
Table 2. Numerical results for f ( 0 ) ,   g ( 0 ) and for θ ( 0 ) for various values of S t = 0.1 ,   0.2 ,   0.3 ,   0.4 , while other parameters are fixed.
S t f ( 0 ) g ( 0 ) θ ( 0 )
0.1 0.9840428504402670.4705787799404170.386755927577512
0.2 0.9738346477763950.4674628686214910.379378708437730
0.3 0.9637459023627560.4643690509223720.372022244909067
0.4 0.9537753196620170.4612971604711970.364686553985701
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ullah, Z.; Akkurt, N.; Alrihieli, H.F.; Eldin, S.M.; Alqahtani, A.M.; Hussanan, A.; Ashraf, M.; Jabeen, M. Temperature-Dependent Density and Magnetohydrodynamic Effects on Mixed Convective Heat Transfer along Magnetized Heated Plate in Thermally Stratified Medium Using Keller Box Simulation. Appl. Sci. 2022, 12, 11461. https://doi.org/10.3390/app122211461

AMA Style

Ullah Z, Akkurt N, Alrihieli HF, Eldin SM, Alqahtani AM, Hussanan A, Ashraf M, Jabeen M. Temperature-Dependent Density and Magnetohydrodynamic Effects on Mixed Convective Heat Transfer along Magnetized Heated Plate in Thermally Stratified Medium Using Keller Box Simulation. Applied Sciences. 2022; 12(22):11461. https://doi.org/10.3390/app122211461

Chicago/Turabian Style

Ullah, Zia, Nevzat Akkurt, Haifaa F. Alrihieli, Sayed M. Eldin, Aisha M. Alqahtani, Abid Hussanan, Muhammad Ashraf, and Mah Jabeen. 2022. "Temperature-Dependent Density and Magnetohydrodynamic Effects on Mixed Convective Heat Transfer along Magnetized Heated Plate in Thermally Stratified Medium Using Keller Box Simulation" Applied Sciences 12, no. 22: 11461. https://doi.org/10.3390/app122211461

APA Style

Ullah, Z., Akkurt, N., Alrihieli, H. F., Eldin, S. M., Alqahtani, A. M., Hussanan, A., Ashraf, M., & Jabeen, M. (2022). Temperature-Dependent Density and Magnetohydrodynamic Effects on Mixed Convective Heat Transfer along Magnetized Heated Plate in Thermally Stratified Medium Using Keller Box Simulation. Applied Sciences, 12(22), 11461. https://doi.org/10.3390/app122211461

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop