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

<jats:p>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.</jats:p>

section_content: 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 nonisothermal 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 electricalconductive 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.
section_name: Introduction and Literature Review

section_content: 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 velocityfields 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:

Appl. Sci. 2022, 12, 11461 5 of 16

The boundary conditions of the present model are,

The proper boundary conditions for the temperature with temperature-dependent density and the thermal stratification components are described above. The Reynolds number is =
section_name: Flow Problem and Mathematical Description

section_content: The proper dimensionless similarity and stream function to convert PDEs into ODEs are given as: 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 ρ

The boundary conditions of the present model are,

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

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.
section_name: The Similarity Transformation and Stream Function Formulation

section_content: The proper dimensionless similarity and stream function to convert PDEs into ODEs are given as:

The stream function φ and ψ are for magnetic and velocity fields, θ is dimensilized temperature and similarity parameter is η. 

The converted form of the ODEs in Equations ( 1 )-( 6 ) by applying the Equations ( 7 ) and ( 8 ), we have non linear ODEs given below;

1e -n + e nθ Mgg = 0 (9)

where

is the mixed convection parameter, η is similarity variable, θ is dimensionless temperature, ν = µ ρ is kinematic fluid viscosity, and St = c b is stratification number. The transformed boundary conditions in (6) then becomes,

The formula for skin friction is C f x = τ w ρU ∞ 2, the formula for the Nusselt number is

, and the formula for the magnetic intensity coefficient is Mg x = -j w U ∞ 2. The values of τ w, q w, and j w are listed below:

Consequently, the values of magnetic intensity, Nusselt number, and skin friction are provided by,
section_name: The Similarity Transformation and Stream Function Formulation

section_content: 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 ),

Equations ( 9 )-( 12 ) get simpler to solve in order to get around the problems that they present. Thus, the more simple forms are,

The reduced boundary conditions are

p → 1, u → 0, θ → 0, as η → ∞. 

We now consider the segment η n-1, η n with η n-1 2 as midpoint given in Equation ( 23 ) below;

This is how the central difference form and average form are presented:

and

By applying the above Equations ( 24 )-( 29 ), the governing Equations ( 14 )-( 21 ) become

along with boundary conditions

Applying the below-mentioned iterative Newton-Raphson approach for smooth algorithm;

By ignoring all appearances of powers of greater than the first power, just like in the standard Newton-Raphson method, equations become;

The simplified form of the equations is provided below, again utilizing Equations ( 35 )-(39) in Equations ( 30 )-( 33 );

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

δp n = 1, δθ n = 0, δu n = 0
section_name: The Solution and Computing Techniques

section_content: 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,
section_name: The Arrangement of Difference Equations in Vector Notation/Matrix

section_content: 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 St, 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 St = 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 St = 0. 0 but the minimum value wasnoted at the higher value of St = 0. 3 in Figure 3a. It wasdeduced that the magnetic field wasenhanced at the large value of St = 0. 3 but decreased at the small value of St = 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 St = 0. 0 but decreased at the maximum value of St = 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. 

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. The good variations in each profile wereexamined for various values of Prandtl number. Figure 6a-c are demonstrated for velocity ′ (), temperature (), and magnetic field ′() for different quantities of magnetic Prandtl number = 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 = 3. 0 but the minimum value wasobtained at a lower value of = 1. 0 with good response. It was deduced that the magnetic graph increased at the large value of = 3. 0 but decreased at the small value of = 1. 0 in Figure 6b. In Figure 6c, it is depicted that the temperature distribution increased at the small value of = 1. 0 but the minimum value wasnoted at the large value of = 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 = 0. 1, 5. 0, 10. 0, and 20. 0with 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 = 20. 0 but minimum value wasobtained at the lower value of = 0. 1 by satisfying the given boundary conditions asymptotically. This result was expected because an increase in 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 = 20. 0 but the minimum value was obtained at small value of = 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, 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 = 20. 0 but increased at the small value of = 0. 1. This given mechanism is performed with good agreement. From Table 1, it is concluded that the skin friction ′′ (0) increased at = 0. 5 but the minimum skin friction wasobtained at small = 0. 0 in the presence of buoyancy forces. It should be noted that = 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 - ′′ (0)wasmaximum at larger = 0. 5 but the smaller quantity of magnetic intensity wasexamined at smaller = 0. 0 in the presence of strong magnetic field. This result was expected because an increase in 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 = 0. 5 but the lower quantity of heat transfer is depicted at minimum = 0. 0 in the presence of magnetic force. With the increase in value of, 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 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 = 0. 1 but the minimum quantity of skin friction is depicted at a large value of = 0. 4 in the presence of Prandtl number = 7. 0. Physically, it is possible because thermal conductivity decreases as Pr increases, which has a lower magnitude of frictional forces 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 St 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 St = 0. 1 but the minimum quantity of skin friction is depicted at a large value of St = 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 St = 0. 4 but the maximum value of magnetic intensity is explored at small St = 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 St = 0. 1 but the minimum quantity of heat transfer wasnoticed at larger St = 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.
section_name: Analysis and Discussion of Results

section_content: 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 MAT-LAB 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 P r = 7. 0. Physically, it is possible because thermal conductivity decreases as P r 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.
section_name: Conclusions

section_name: Acknowledgments:

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

section_content: 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.

section_content: Data Availability Statement: Not applicable.

section_name: Conflicts of Interest:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: