Next Article in Journal
Application of the k-Prototype Clustering Approach for the Definition of Geostatistical Estimation Domains
Next Article in Special Issue
Hydrothermal Mixed Convection in a Split-Lid-Driven Triangular Cavity Suspended by NEPCM
Previous Article in Journal
A Combinatorial Approach to Study the Nordhaus–Guddum-Type Results for Steiner Degree Distance
Previous Article in Special Issue
Implicit Finite-Difference Scheme for a Duffing Oscillator with a Derivative of Variable Fractional Order of the Riemann-Liouville Type
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Buoyancy Force on an Unsteady Thin Film Flow of Al2O3/Water Nanofluid over an Inclined Stretching Sheet

1
Department of Mathematical Sciences, Faculty of Science and Technology, Universiti Kebangsaan Malaysia (UKM), Bangi 43600, Selangor, Malaysia
2
Department of Mathematics, Faculty of Science, Qassim University, Buraydah 52571, Saudi Arabia
3
Fakulti Teknologi Kejuruteraan Mekanikal dan Pembuatan, Universiti Teknikal Malaysia Melaka, Hang Tuah Jaya, Durian Tunggal 76100, Melaka, Malaysia
4
Department of Mathematics, Faculty of Science, King Khalid University, Abha 62529, Saudi Arabia
5
Department of Mathematics, Faculty of Science, South Valley University, Qena 83523, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(3), 739; https://doi.org/10.3390/math11030739
Submission received: 24 December 2022 / Revised: 17 January 2023 / Accepted: 28 January 2023 / Published: 1 February 2023
(This article belongs to the Special Issue Computational Mathematics and Mathematical Modelling)

Abstract

:
The present study looks at the heat transfer and the unsteady thin film flow of Al2O3 water nanofluid past an inclined stretching sheet having a buoyancy force effect. The boundary value problem solver (bvp4c) package in Matlab is utilized in solving the converted set of ordinary differential equations (ODEs). The multi-shape Al2O3 nanoparticles’ impact with respect to the flow as well as heat transfer characteristics are studied and visually displayed for certain governing parameter values, which include the mixed convection, inclination angle, magnetic, slip, and Biot number. Thus, the skin friction coefficient and the local Nusselt number are also determined. Here, the platelet shape of Al2O3 nanoparticles possesses a high heat transfer and flow rate based on the outcomes. In addition, increasing the slip and magnetic parameters improves the temperature, whereas increasing the buoyancy and inclination angle parameters has reverse effects. The results also show that increasing the unsteadiness parameter and the magnetic parameter reduces the film thickness.

1. Introduction

Over the last few decades, the topic of improving heat transfer in industrial applications has been of great relevance. The majority of applications use pure fluids, such as water, oil, and ethylene, as cooling liquids. However, these fluids have low heat transfer rates. To improve heat transfer rates, nanometer-sized particles are scattered into the base fluid to improve its thermal performance and rate of heat transfer [1]. Choi and Eastman were the first to refer to this combination as nanofluid [2]. The addition of nanoparticles into the base fluid significantly improves the heat transfer and thermal conductivity of these fluids. Since this research, nanofluids have been extensively employed in heating and cooling systems, nuclear applications, power generation, medicine, electronics, automobiles, and other applications [3,4,5,6]. Consequently, numerous numerical and experimental studies have been performed to investigate the nanofluids’ properties and behaviors in various aspects [7,8,9,10,11,12,13]. Sajjadi et al. [14] scrutinized the MHD natural convection flow of a nanofluid in a 3-D cavity where the sinusoidal temperature is distributed on a side wall. They concluded that the augmentation of the volumetric fraction of nanoparticles accelerates the rate of heat transfer. Afterwards, Atashafrooz et al. [15] studied the entropy generation due to an inclined step using two types of nanofluids. According to the predictions, they found that the Al2O3/water nanofluid is better suited for the hydrothermal system’s ideal design.
The boundary layer flow originated from a stretching sheet is of great significance as a result of its critical role in industry and engineering. Crane [16] was the pioneer in investigating the 2D boundary layer flow brought on by a stretching sheet. Furthermore, Gupta and Gupta [17] extended earlier research by considering the case of stretching surfaces under the effect of suction or blowing. Since then, several authors have been working to solve the boundary layer flow issue with respect to a stretching surface by considering various nanofluids under varied physical conditions, such as Sharma et al. [18], Zokri et al. [19], Kho et al. [20], Daniel et al. [21], Ahmed et al. [22], Kho et al. [23], and Ibrahim and Negera [24]. Numerical solutions of MHD Maxwell nanofluid flow induced by stretching surface was scrutinized by Prasannakumara [25]. He reported that the surface drag force is enhanced by increasing the volume fraction. Hazarika et al. [26] used a numerical method to resolve the stretching sheet problem for three different types of nanoparticles, including Ferrous Ferric Oxide ( F e 3 O 4 ) , Silver ( A g ) , and Copper ( C u ) . They claimed that F e 3 O 4 -water nanofluid always overshoots the velocity, followed by C u and A g -water nanofluid, although the temperature profiles exhibit the opposite behaviours.
Furthermore, the horizontal and vertical positions of the flat plate have been utilized in numerous research studies. The inclined flat plate, on the other hand, has received less attention. Thus, it is beneficial to examine the fluid flow induced by an inclined flat plate at diverse angles as are used in many engineering processes. In addition, Gupta et al. [27] examined the radiation and chemical reactions on the laminar magnetohydrodynamics (MHD) stagnation point flow, including heat transfer in a non-Newtonian nanofluid caused by an inclined stretching surface. Usman et al. [28] investigated the combination effects of thermal and velocity slips on Casson nanofluid flow past a permeable inclined stretching cylinder. They discovered that the thermophoresis and Brownian motion impacts play a prominent role in the concentration and thermal boundary layers. Moreover, the unsteady MHD flow using Eyring-Powell nanofluid due to an inclined permeable stretching surface was numerically examined by Kumar and Srinivas [29]. They considered the influences of Joule heating and thermal radiation. Afterwards, a numerical computation of MHD hybrid nanofluid flow driven by an inclined sheet was scrutinized by Soomro et al. [30] who found that Al2O3- C u /water hybrid nanofluid increased the heat transfer rate in contrast to the basic C u /water nanofluid.
Many scholars have recently become interested in problems with unsteady boundary layers because of their significant role in a variety of engineering issues. Start-up processes, for instance, denote the motions at rest in transit between one steady flow and another, as well as periodic working fluid motion [31]. Because of the extra time-based factors that impeded the separation of boundary layer and fluid motion pattern, unsteady boundary layer flow activity seems to behave differently than a steady flow [32,33]. The unsteady mixed convection flow, as well as heat transfer caused by a stretching vertical sheet, were examined by Ishak et al. [34]. Moreover, Ishak et al. [35] then looked at how MHD affected the unsteady flow, as well as heat transfer past a stretching surface, and deduced that the rate of heat transfer rises with a rise of the unsteadiness parameter. In contrast, an increment in the magnetic parameter results in a decrease in the heat transfer rate. Daniel et al. [36] performed an investigation on the unsteady MHD natural convection flow problem with respect to a viscous nanofluid due to a permeable shrinking surface which took into account the influences of chemical reaction, electric field, thermal radiation, and viscous dissipation. Dzulkifli et al. [37] also looked at the effect of slip on the unsteady stagnation-point flow and heat transfer caused by a porous exponentially shrinking/stretching surface. They claimed that there are dual solutions and noted that discovering the dual solutions in the case of suction is much simpler than in the case of injection. In addition, Waini et al. [1] examined the 2D unsteady heat transfer and hybrid nanofluid flow numerically concerning a shrinking/stretching surface. Dual solutions were found for a given range of the unsteadiness parameter. Thereafter, Waini et al. [38] employed a hybrid A l 2 O 3 - C u /water nanofluid produced by a permeable rigid sheet to conduct a temporal stability analysis of an unsteady stagnation point flow. A numerical analysis of the unsteady hybrid nanofluid flow brought on by a permeable stretching/shrinking cylinder was also carried out by Zainal et al. [39]. They stated that including the unsteadiness parameter significantly aided in the heat transfer deterioration.
Given the vital uses in industry and engineering--for instance, paper manufacture, polymer and foodstuff processing, wire thinning, crystal growth, and so on--the thin film flow has recently been a demanding research issue. A thin layer of liquid flows past a solid surface because of an external force acting on it, such as gravity, shear stress, or other factors known as film flow. For example, Wang [40] was the first to investigate the thin film flow problem due to an unsteady stretching surface. Later, Andersson et al. [41] improved the previous work by considering the analysis of heat transfer. In addition, Khan et al. [42] provided an analytical solution to the heat transfer problem and the second-grade fluid’s thin film flow across a stretched surface in a porous medium. The effect of film thickness was explored, and it was discovered that, when film thickness decreases, the skin friction coefficient increases. In addition, Naganthran [43] reported the impacts of injection and thermocapillarity on the Carreau thin film flow, including heat transfer in relation to an unsteady stretching surface. It has been demonstrated that increasing the injection intensity and unsteadiness parameter reduces the film thickness, resulting in a rise in the convective heat transfer rate and skin friction coefficient. Afterwards, Naganthran [44] performed an investigation on the flow of the Carreau thin hybrid nanofluid film caused by an accelerating surface together with the melting heat transfer influence. They concluded that melting heat transfer slows the pace of heat transfer rate while having no effect on the liquid film thickness. Moreover, Ali et al. [45] utilized Al2O3 nanofluid past a stretching sheet with convective boundary conditions in examining the unsteady thin film flow dynamic and heat transfer.
A fair amount of substantial analysis on the thin film flow due to a stretching sheet has been carried out by researchers. However, the thin film flow on an inclined stretching plate is rarely examined, motivating the current research to investigate the thin film flow on an inclined stretching surface. Under the multi-shape nanoparticles’ effect, the primary goal of this research is to investigate the unsteady thin film flow as well as heat transfer with respect to the inclined stretching sheet in Al2O3/water nanofluid. Moreover, the convective boundary conditions and buoyancy force effects have also been considered. By employing the similarity transformation, the provided system is transformed to a system of ordinary differential equations (ODEs) which are later numerically solved using the boundary value problem solver (bvp4c) available in Matlab software. This research has investigated the impacts of the distinct factors with regard to the temperature fields as well as the velocity fields. In addition, the acquired numerical findings and the published results available in the literature have been compared for the purpose of validation. The above literature review reveals that no studies of this kind have been performed by any other researchers.

2. Mathematical Formulation

Consider a thin film flow and heat transfer in a nanofluid across an unsteady stretching sheet. Here, the matching physical model is portrayed in Figure 1, with the x -axis parallel to the plate and the y -axis perpendicular to it. Also, the current study examines the effects of four different nanoparticle shapes, including the platelet, cylinder, brick, and sphere-shaped nanoparticles. It is assumed that the base fluid and the multi-shape nanoparticles are in thermal equilibrium. Further, U w = b x 1 a t represents the stretching surface velocity in the x direction, where t is time and a , b are dimensional constants. The uniform magnetic field is assumed as B ( t ) = B 0 1 a t , and the temperature distribution at the wall is represented by T w = T + T 0 ( x ( 1 a t ) 2 ) , where T and T 0 are ambient temperature and constant reference temperature, see Ref. [45].
From all of the considerations above and using the nanofluid model introduced by Tiwari and Das [46], the boundary layer governing equations for the present problem are given by [45,47].
u x + v y = 0 ,
u t + u u x + v u y = μ n f ρ n f 2 u y 2 σ n f ρ n f B ( t ) 2 u + ( ρ β * ) n f ρ n f g ( T T ) cos α ,
T t + u T x + v T y = k n f ( ρ c p ) n f 2 T y 2 ,
subject to the boundary conditions [45].
y = 0 : u = U w + A ν f u y ,   v = 0 ,   k T y = h f [ T 0 T ] , y = h ( t ) : u y = 0 , T y = 0 , v = d h d t .
Here, the nanofluid velocity components along x and y directions are represented by u and v , respectively. Moreover, T is the fluid temperature, g is the gravitational acceleration, α is the inclination angle, h ( t ) is the thickness of the film, A = A 0 1 a t represent the velocity slip factor with constant A 0 , and h f = h 0 1 a t denotes the convective heat transfer with h 0 being constant. Furthermore, ρ n f , μ n f , σ n f , ( ρ β * ) n f , k n f , and ( ρ c p ) n f denote the density, dynamic viscosity, electrical conductivity, diffusivity, thermal expansion coefficient, thermal conductivity, and the heat capacity of the nanofluid, respectively, which are defined as the following [45,48]:
ρ n f = ( 1 φ ) ρ f + φ ρ s , ( ρ c p ) n f = ( 1 φ ) ( ρ c p ) f + φ ( ρ c p ) s , μ n f = μ f ( 1 + A 1 φ + A 2 φ 2 ) , σ n f = σ f ( 1 φ ) + φ σ s , ( ρ β * ) n f = ( 1 φ ) ( ρ β * ) f + φ ( ρ β * ) s , k n f k f = k s + ( m 1 ) k f + ( m 1 ) ( k s k f ) φ k s ( m 1 ) k f ( k s k f ) φ ,
where φ is the volume fraction of the nanoparticle, and the viscosity enhancement heat capacitance coefficients are symbolized by A 1 and A 2 . Meanwhile, the base fluid, nanofluid and nanoparticles are signified as the subscripts f , n f and s , respectively, while m implies the shape factor. Furthermore, it is important to note that the convective heat transfer property of nanofluids depends on the thermophysical characteristics of the suspended particles and the base fluid; the flow structure and pattern; the dimensions and shape of nanoparticles; and the volumetric fraction of nanoparticles. The thermophysical properties, shape factor of various geometry in this study, and coefficients of viscosity for Al2O3 nanoparticles along with water are presented in Table 1 and Table 2 [49,50,51].
Next, the following similarity transformations are employed [45]:
θ ( η ) = T T T w T , η = b ν f ( 1 a t ) y , ψ = b ν f 1 a t x f ( η ) ,
where ψ indicates the stream function which is defined as u = ψ y , and v = ψ x , ν f denotes the kinematic viscosity and β is the dimensionless film thickness, expressed as
β = b ν f ( 1 a t ) h ( t )
which will eventually result in
d h d t = a β 2 ν f b ( 1 a t ) ,
implying that in the similarity variables (6) in Equations (1)–(4), Equation (1) is identically fulfilled while the remaining equations become
ε 1 f ( η ) M ε 3 f ( η ) S [ f ( η ) + η 2 f ( η ) ] f 2 ( η ) + f ( η ) f ( η ) + ε 4 λ θ cos α = 0 ,
ε 2 P r θ ( η ) S 2 ( 4 θ ( η ) + η θ ( η ) ) θ ( η ) f ( η ) + θ ( η ) f ( η ) = 0 ,
with the transformed boundary conditions:
η = 0 : f ( 0 ) = 0 ,   f ( 0 ) = 1 + K f ( 0 ) ,   θ ( 0 ) = k f k n f B i ( 1 θ ( 0 ) ) , η = β : f ( β ) = S β 2 ,   f ( β ) = 0 ,   θ ( β ) = 0 .
Here, M is the magnetic parameter, S is the unsteadiness parameter, P r is the Prandtl number, K is the slip parameter, B i is the Biot number, and λ is the buoyancy or mixed convection parameter, which are defined as
M = σ f B 0 2 ρ n f b ,   S = a b ,   P r = ( ρ c p ) f ν f k f ,   K = A ν f U w x = A 0 ν f b , B i = h f k f x ν f U w = h 0 k f ν f b ,   λ = G r R e x 2 = T 0 g β f * b 2 ,
where G r = g β f * ( T w T ) x 3 ν f 2 and R e x = U w x ν f represent, respectively, the Grashof number and the local Reynold number [52]. Notice that λ > 0 and λ < 0 respectively refer to the assisting and opposing flows. Further, ε 1 , ε 2 , ε 3 , and ε 4 are constants and are given by
ε 1 = 1 + A 1 φ + A 2 φ 2 1 φ + φ ( ρ s ρ f ) ,   ε 2 = k n f k f 1 φ + φ ( ρ c p ) s ( ρ c p ) f ,   ε 3 = 1 φ + φ ( σ s σ f ) 1 φ + φ ( ρ s ρ f ) ,   ε 4 = ( 1 φ ) ( ρ β * ) f + φ ( ρ β * ) s ( ( 1 φ ) ρ f + φ ρ s ) β f * .
The skin friction coefficient and the Nusselt number are the physical quantities of interest which can be respectively expressed as
C f = τ w ρ f U w 2 ,   N u = q w x k f ( T w T ) ,
where
τ w = μ n f ( u y ) y = 0 ,   q w = k n f ( T y ) y = 0 .
By using (5) and inserting (15) into (14), the following relations are obtained:
C f R e 1 / 2 = ( 1 + A 1 φ + A 2 φ 2 ) f ( 0 ) ,    N u R e 1 / 2 = k n f k f θ ( 0 ) .

3. Numerical Method

The package bvp4c in Matlab is utilized for solving the Equations (9) and (10) with the boundary conditions (11), numerically. The bvp4c solver, discussed by Shampine et al. [53], is based on the finite-difference method which uses the 3-stage Lobatto IIIa formula. Lobatto IIIa is a Runge-Kutta method that uses collocation and implicit trapezoidal rules. As a result, the collocation method is indicated by the bvp4c solver, which produces a C1-continuous solution with fourth-order accuracy uniformly in the interval in which the function is integrated. To begin using the bvp4c method, Equations (9) and (10) must be reduced to a set of first-order ordinary differential equations (ODEs), along with their boundary conditions (11), as follows:
f = y ( 1 ) ,   f = y ( 2 ) ,   f = y ( 3 ) ,   θ = y ( 4 ) ,   θ = y ( 5 ) , f = 1 ε 1 ( μ ε 3 y ( 2 ) + S [ y ( 2 ) + η 2 y ( 3 ) ] + y ( 2 ) y ( 2 ) y ( 1 ) y ( 3 ) ε 4 λ y ( 4 ) c o s α ) ,
θ = P r ε 2 ( S 2 ( 3 y ( 4 ) + η y ( 5 ) ) + 2 y ( 2 ) y ( 4 ) y ( 1 ) y ( 5 ) ) .
The boundary conditions (11) become
y a ( 1 ) ,   y a ( 2 ) 1 K y a ( 3 ) ,   y a ( 5 ) + k k n f B i ( 1 y a ( 4 ) ) ,   y b ( 1 ) s β 2 ,   y b ( 3 ) ,   y b ( 5 ) .
Equations (17) and (18) form a set of five first-order ODEs together with six boundary conditions (19). Here, the crucial parameters are highlighted: S the unsteadiness parameter and β the film thickness parameter. There is a relationship between β and S that can lead to the reduction of the six boundary conditions to five. To find the relationship of β and S , Equations (17) and (19) are solved first, and the initial guess value for β is provided by the bvp4c function in Matlab. The value of β is adjusted so that the condition f ( β ) = S β 2 holds. It is achieved using trial-and-error processes. Equations (17)–(19) are solved utilizing the bvp4c solver for known values of S and β .

4. Findings and Discussion

The boundary value problem solver (bvp4c) in Matlab software was applied in solving the mathematical model given in Equations (9)–(11) numerically for the chosen values of the inclination angle parameter α , magnetic parameter M , slip parameter K , mixed convection parameter λ , and Biot number B i . It is essential to note that the range with respect to multi-shape nanoparticles in a thin film past a stretching surface of A l 2 O 3 /water nanofluid affects both temperature profiles and velocity. The parameters’ impact is depicted in tables and graphs. A comparative study of the results obtained by Abel et al. [54] and Ali et al. [45] are given in Table 3. The outcome shows a favorable agreement and ensures that the bvp4c scheme is highly accurate and effective. Moreover, the skin friction coefficient R e 1 / 2 C f for multi-shape nanoparticles is depicted in Table 4. The table demonstrates that, as the mixed convection parameter λ and slip parameter K increase, the skin friction R e 1 / 2 C f decreases. However, when the volume fraction parameter φ and magnetic parameter M increase, the skin friction R e 1 / 2 C f increases. In addition, the heat transfer rate R e 1 / 2 N u is determined and presented in Table 5. Furthermore, it may also be noted that R e 1 / 2 N u rises with the rising values of the inclination angle, Biot number, and mixed convection parameters. Meanwhile, the rising values of the slip and magnetic parameters tend to reduce the values of Nusselt number R e 1 / 2 N u .
Figure 2 illustrates the effect of multi-shape nanoparticles with respect to the film thickness β of A l 2 O 3 nanofluids. In this case, the remaining physical parameters remain constant. Provided that the nanoparticles’ shape changes, the film thickness also drastically changes. Furthermore, the velocity profile f ( η ) and film thickness β accelerate for the platelet-shaped nanoparticles and decelerate for the brick, cylinder, and sphere-shaped nanoparticles. Figure 3 shows the impact of the slip parameter K on the temperature profile θ ( η ) for multi-shape nanoparticles. The temperature profile θ ( η ) for each multi-shape nanoparticle increases as the slip parameter K rises. The enhancement of the slip parameter has resulted in acceleration of the fluid within the boundary layer. Additionally, the current study’s reflection of the convective heating effect, which controls the temperature of the stretched surface, may have contributed to this phenomenon. In reality, increased convection causes higher surface temperatures, allowing the thermal impact to go deeper into the nanofluid. Next, Figure 4 illustrates the magnetic parameter M effect on the temperature profile θ ( η ) . Having higher levels of the magnetic parameter M , the temperature θ ( η ) becomes greater for all types of nanoparticle shapes. The Lorentz force, which acts as a reverse force when magnetic effect is present, causes this. This force increases the fluid’s friction by opposing the fluid’s motion. It is a resistive force; hence, fluid friction causes the temperature to rise. Figure 5 illustrates the effect with respect to the mixed convection parameter λ towards the temperature profile θ ( η ) . Furthermore, note that the temperature θ ( η ) reduces provided the buoyancy parameter boosts (from opposing to assisting). However, it is important to point out that, for platelet and cylinder shapes, the temperature receives less impact. Near the wall, the assisting flow has a lower temperature than the opposing flow. The assisting flow, in theory, transfers heat from the hot plate to the cool fluid, eventually bringing the wall temperature down.
The temperature profile θ ( η ) behavior for different inclination angle parameters α values are shown in Figure 6. It is discovered that, as the angle of inclination parameter α increases, the temperature θ ( η ) for all multi-shape nanosize particles decreases. The buoyancy force increases when the inclination to the horizontal increases. Physically, when the sheet is oriented vertically, the strongest gravitational forces have an impact on the flow. However, if the plate is turned from vertical to horizontal (and as the value of α increases), the buoyancy force increases. Thus, the thermal boundary layer reduces. It is noticed that the problem is reduced to the vertical plate if α = 0 ° , whereas when α = 90 ° the problem is lowered to the horizontal plate. It is noted that for α = 30 ° , 45 ° , 60 ° , the problem reduces to the inclined plate. Figure 7 shows the Biot number B i effect on the temperature profile θ ( η ) . Moreover, it is discovered that raising B i increases the fluid temperature. Additionally, B i is calculated by dividing convection at the surface by conduction inside the body surface. Convection will become more potent near the surface as a consequence of the increase in B i , leading to an increase in the thickness of the thermal boundary layer. Additionally, as illustrated in Figure 8, the temperature rises as the film’s thickness reduces as the shape of the nanoparticles shifts from platelet to spherical. It should be noted that the shape factor has a greater impact on platelet shape than other factors. In Figure 9, the magnetic parameter M and unsteadiness parameter S impacts on the film thickness β are depicted. The decreasing trends in film thickness β are shown in this figure for every multi-shape nanoparticle when the unsteadiness parameter S increases, coupled with a rise in the magnetic parameter M . Also, when S 0 , as shown in Figure 9, it represents the liquid film thickness of an infinite thickness ( β ) . For the particular values S 0 of S , no solutions can be obtained, while S S 0 represents the case of an infinitesimal thick fluid layer ( β 0 ) .

5. Conclusions

A numerical study has been done on the unsteady thin film flow and the heat transfer of A l 2 O 3 /water nanofluid steered by an inclined stretching sheet. Applying similarity transformation, the governing partial differential equations (PDEs) are reduced to a set of ordinary differential equations (ODEs), which are then numerically solved via the boundary value problem solver (bvp4c) in the Matlab software. The findings are listed below:
  • Provided that the magnetic parameter M and the slip parameter K rise, the film thickness β reduces.
  • The temperature θ(η) enhances as the Biot number Bi boosts. For platelet-shaped particles, the film thickness is high, whereas, for cylinder, brick, and sphere-shaped nanoparticles, it ultimately diminishes.
  • On the platelet and sphere shape nanoparticles, the velocity of Al2O3 nanofluid achieves its maximum and minimum, whereas a different pattern can be seen in the temperature.
  • The skin friction coefficient lessens when the slip and mixed convection parameters rise. In contrast, incrementally the volume fraction and magnetic parameters increase the skin friction coefficient.
  • The local Nusselt number increases as the Biot number Bi, inclination angle, and mixed convection parameters increase, while it declines with an augmentation in the slip and magnetic parameters.
  • The assisting buoyancy flow (λ > 0) appears to have a greater local Nusselt number than the opposing buoyancy flow (λ < 0), despite the fact that the skin friction coefficient produces the opposite outcome.
  • A decrease in film thickness is produced by a rise in the unsteadiness parameter, as well as the magnetic parameter.

Author Contributions

Conceptualization, A.I., S.E.A.; funding acquisition, A.I.; methodology, S.A., I.W.; Project administration, A.I.; supervision, A.I., S.A.B., S.E.A.; validation, I.W., S.A.B., A.I., S.E.A.; writing—original draft, S.A.; writing—review and editing, A.I., S.A.B., I.W., S.E.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universiti Kebangsaan Malaysia (Project Code: DIP-2020-001).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a : b dimensional constants ( s 1 )
B i Biot number
B 0 uniform magnetic field
c p specific heat of the fluid ( J   k g 1 K 1 )
g acceleration due to gravity ( m   s 2 )
G r Grashof number
h ( t ) film thickness ( m )
k thermal conductivity of the fluid ( W   m 1 K 1 )
M magnetic parameter
N u Nusselt number
P r Prandtl number
R e Reynolds number
S unsteadiness parameter
t time ( s )
T fluid temperature ( K )
T w surface temperature at the wall ( K )
T 0 reference temperature ( K )
T ambient temperature ( K )
U w surface velocity ( m   s 1 )
u , v velocity components in the x and y directions ( m   s 1 )
x , y Cartesian coordinates ( m )
Greek symbols
α inclination angle
β dimensionless film thickness
β * thermal expansion coefficient of the fluid ( K 1 )
η similarity variable
λ mixed convection parameter
μ dynamic viscosity of the fluid ( k g   m 1 s 1 )
ν kinematic viscosity of the fluid ( m 2   s 1 )
φ dimensionless nanoparticle volume fraction
ψ stream function ( m 2   s 1 )
ρ density of the fluid ( k g   m 3 )
( ρ c p ) heat capacity of the fluid ( J   K 1 m 3 )
σ electrical conductivity of the fluid ( s   m 1 )
θ dimensionless temperature
Subscripts
f base fluid
n f Nanofluid
s solid nanoparticle
w condition at the wall
Superscript
differentiation with respect to η

References

  1. Waini, I.; Ishak, A.; Pop, I. Unsteady flow and heat transfer past a stretching/shrinking sheet in a hybrid nanofluid. Int. J. Heat Mass Transf. 2019, 136, 288–297. [Google Scholar] [CrossRef]
  2. Choi, S.U.S.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Argonne National Lab.: Lemont, IL, USA, 1995. [Google Scholar]
  3. Choi, S.U.S. Nanofluids: From vision to reality through research. J. Heat Transfer. 2009, 131, 033106. [Google Scholar] [CrossRef]
  4. Senthilraja, S.; Karthikeyan, M.; Gangadevi, R. Nanofluid applications in future automobiles: Comprehensive review of existing data. Nano-Micro Lett. 2010, 2, 306–310. [Google Scholar] [CrossRef]
  5. Saidur, R.; Leong, K.Y.; Mohammed, H.A. A review on applications and challenges of nanofluids. Renew. Sustain. Energy Rev. 2011, 15, 1646–1668. [Google Scholar] [CrossRef]
  6. Sheremet, M.A. Applications of nanofluids. Nanomaterials 2021, 11, 1716. [Google Scholar] [CrossRef]
  7. Waini, I.; Ishak, A.; Pop, I. Nanofluid flow on a shrinking cylinder with Al2O3 nanoparticles. Mathematics 2021, 9, 1612. [Google Scholar] [CrossRef]
  8. Hamid, K.A.; Azmi, W.H.; Mamat, R.; Sharma, K. V Heat transfer performance of TiO2–SiO2 nanofluids in a tube with wire coil inserts. Appl. Therm. Eng. 2019, 152, 275–286. [Google Scholar] [CrossRef]
  9. Safwa Khashi’ie, N.; Md Arifin, N.; Hafidzuddin, E.H.; Wahi, N. Dual stratified nanofluid flow past a permeable shrinking/stretching sheet using a non-Fourier energy model. Appl. Sci. 2019, 9, 2124. [Google Scholar] [CrossRef]
  10. Anwar, M.I.; Rafique, K.; Misiran, M.; Khan, I. Numerical solution of casson nanofluid flow over a non-linear inclined surface with soret and dufour effects by keller-box method. Front. Phys. 2019, 7, 139. [Google Scholar]
  11. Bachok, N.; Ishak, A.; Pop, I. Boundary-layer flow of nanofluids over a moving surface in a flowing fluid. Int. J. Therm. Sci. 2010, 49, 1663–1668. [Google Scholar] [CrossRef]
  12. Zaimi, K.; Ishak, A.; Pop, I. Flow past a permeable stretching/shrinking sheet in a nanofluid using two-phase model. PLoS ONE 2014, 9, e111743. [Google Scholar] [CrossRef]
  13. Salleh, S.N.A.; Bachok, N.; Arifin, N.M.; Ali, F.M. Numerical analysis of boundary layer flow adjacent to a thin needle in nanofluid with the presence of heat source and chemical reaction. Symmetry 2019, 11, 543. [Google Scholar] [CrossRef]
  14. Sajjadi, H.; Delouei, A.A.; Atashafrooz, M.; Sheikholeslami, M. Double MRT Lattice Boltzmann simulation of 3-D MHD natural convection in a cubic cavity with sinusoidal temperature distribution utilizing nanofluid. Int. J. Heat Mass Transf. 2018, 126, 489–503. [Google Scholar] [CrossRef]
  15. Atashafrooz, M.; Sajjadi, H.; Amiri Delouei, A.; Yang, T.-F.; Yan, W.-M. Three-dimensional analysis of entropy generation for forced convection over an inclined step with presence of solid nanoparticles and magnetic force. Numer. Heat Transf. Part A Appl. 2021, 80, 318–335. [Google Scholar] [CrossRef]
  16. Crane, L.J. Flow past a stretching plate. Zeitschrift für Angew. Math. und Phys. 1970, 21, 645–647. [Google Scholar] [CrossRef]
  17. Gupta, P.S.; Gupta, A.S. Heat and mass transfer on a stretching sheet with suction or blowing. Can. J. Chem. Eng. 1977, 55, 744–746. [Google Scholar] [CrossRef]
  18. Sharma, R.; Ishak, A.; Pop, I. Partial slip flow and heat transfer over a stretching sheet in a nanofluid. Math. Probl. Eng. 2013, 2013, 724547. [Google Scholar] [CrossRef]
  19. Zokri, S.M.; Arifin, N.S.; Salleh, M.Z.; Kasim, A.R.M.; Mohammad, N.F.; Yusoff, W. MHD Jeffrey nanofluid past a stretching sheet with viscous dissipation effect. J. Phys. Conf. Ser. 2017, 890, 12002. [Google Scholar] [CrossRef]
  20. Kho, Y.B.; Hussanan, A.; Mohamed, M.K.A.; Sarif, N.M.; Ismail, Z.; Salleh, M.Z. Thermal radiation effect on MHD Flow and heat transfer analysis of Williamson nanofluid past over a stretching sheet with constant wall temperature. J. Phys. Conf. Ser. 2017, 890, 12034. [Google Scholar] [CrossRef]
  21. Daniel, Y.S.; Aziz, Z.A.; Ismail, Z.; Salah, F. Thermal stratification effects on MHD radiative flow of nanofluid over nonlinear stretching sheet with variable thickness. J. Comput. Des. Eng. 2018, 5, 232–242. [Google Scholar] [CrossRef]
  22. Ahmed, S.E.; Mohamed, R.A.; Abd Elraheem, M.A.; Soliman, M.S. Magnetohydrodynamic Maxwell nanofluids flow over a stretching surface through a porous medium: Effects of non-linear thermal radiation, convective boundary conditions and heat generation/absorption. World Acad. Sci. Eng. 2019, 13, 436–443. [Google Scholar]
  23. Kho, Y.B.; Hussanan, A.; Sarif, N.M.; Ismail, Z.; Salleh, M.Z. Thermal radiation effects on MHD with flow heat and mass transfer in Casson nanofluid over a stretching sheet. In Proceedings of the MATEC Web of Conferences, Penang, Malaysia, 6–7 December 2017; EDP Sciences: Ulis, France, 2018; Volume 150, p. 6036. [Google Scholar]
  24. Ibrahim, W.; Negera, M. MHD slip flow of upper-convected Maxwell nanofluid over a stretching sheet with chemical reaction. J. Egypt. Math. Soc. 2020, 28, 7. [Google Scholar] [CrossRef] [Green Version]
  25. Prasannakumara, B.C. Numerical simulation of heat transport in Maxwell nanofluid flow over a stretching sheet considering magnetic dipole effect. Partial Differ. Equ. Appl. Math. 2021, 4, 100064. [Google Scholar] [CrossRef]
  26. Hazarika, S.; Ahmed, S.; Chamkha, A.J. Investigation of nanoparticles Cu, Ag and Fe3O4 on thermophoresis and viscous dissipation of MHD nanofluid over a stretching sheet in a porous regime: A numerical modeling. Math. Comput. Simul. 2021, 182, 819–837. [Google Scholar] [CrossRef]
  27. Gupta, S.; Kumar, D.; Singh, J. MHD mixed convective stagnation point flow and heat transfer of an incompressible nanofluid over an inclined stretching sheet with chemical reaction and radiation. Int. J. Heat Mass Transf. 2018, 118, 378–387. [Google Scholar] [CrossRef]
  28. Usman, M.; Soomro, F.A.; Haq, R.U.; Wang, W.; Defterli, O. Thermal and velocity slip effects on Casson nanofluid flow over an inclined permeable stretching cylinder via collocation method. Int. J. Heat Mass Transf. 2018, 122, 1255–1263. [Google Scholar] [CrossRef]
  29. Kumar, B.; Srinivas, S. Unsteady hydromagnetic flow of Eyring-Powell nanofluid over an inclined permeable stretching sheet with joule heating and thermal radiation. J. Appl. Comput. Mech. 2020, 6, 259–270. [Google Scholar]
  30. Soomro, F.A.; Usman, M.; El-Sapa, S.; Hamid, M.; Haq, R.U. Numerical study of heat transfer performance of MHD Al2O3-Cu/water hybrid nanofluid flow over inclined surface. Arch. Appl. Mech. 2022, 92, 2757–2765. [Google Scholar] [CrossRef]
  31. Schlichting, H.; Gersten, K. Boundary-Layer Theory; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2003; ISBN 3540662707. [Google Scholar]
  32. Smith, F.T. Steady and unsteady boundary-layer separation. Annu. Rev. Fluid Mech. 1986, 18, 197–220. [Google Scholar] [CrossRef]
  33. White, F.M.; Majdalani, J. Viscous Fluid Flow; McGraw-Hill: New York, NY, USA, 2006. [Google Scholar]
  34. Ishak, A.; Nazar, R.; Pop, I. Boundary layer flow and heat transfer over an unsteady stretching vertical surface. Meccanica 2009, 44, 369–375. [Google Scholar] [CrossRef]
  35. Ishak, A. Unsteady MHD flow and heat transfer over a stretching plate. J. Appl. Sci. 2010, 10, 2127–2131. [Google Scholar] [CrossRef]
  36. Daniel, Y.S.; Aziz, Z.A.; Ismail, Z.; Salah, F. Slip effects on electrical unsteady MHD natural convection flow of nanofluid over a permeable shrinking sheet with thermal radiation. Eng. Lett. 2018, 26, 1–13. [Google Scholar]
  37. Dzulkifli, N.F.; Bachok, N.; Yacob, N.A.; Md Arifin, N.; Rosali, H. Unsteady stagnation-point flow and heat transfer over a permeable exponential stretching/shrinking sheet in nanofluid with slip velocity effect: A stability analysis. Appl. Sci. 2018, 8, 2172. [Google Scholar] [CrossRef]
  38. Waini, I.; Ishak, A.; Pop, I. Unsteady hybrid nanofluid flow on a stagnation point of a permeable rigid surface. ZAMM-J. Appl. Math. Mech. für Angew. Math. und Mech. 2021, 101, e202000193. [Google Scholar] [CrossRef]
  39. Zainal, N.A.; Nazar, R.; Naganthran, K.; Pop, I. Stability Analysis of Unsteady Hybrid Nanofluid Flow Past a Permeable Stretching/Shrinking Cylinder. J. Adv. Res. Fluid Mech. Therm. Sci. 2021, 86, 64–75. [Google Scholar] [CrossRef]
  40. Wang, C.Y. Liquid film on an unsteady stretching surface. Q. Appl. Math. 1990, 48, 601–610. [Google Scholar] [CrossRef]
  41. Andersson, H.I.; Aarseth, J.B.; Dandapat, B.S. Heat transfer in a liquid film on an unsteady stretching surface. Int. J. Heat Mass Transf. 2000, 43, 69–74. [Google Scholar] [CrossRef]
  42. Khan, N.S.; Islam, S.; Gul, T.; Khan, I.; Khan, W.; Ali, L. Thin film flow of a second grade fluid in a porous medium past a stretching sheet with heat transfer. Alexandria Eng. J. 2018, 57, 1019–1031. [Google Scholar] [CrossRef]
  43. Naganthran, K.; Hashim, I.; Nazar, R. Triple Solutions of Carreau Thin Film Flow with Thermocapillarity and Injection on an Unsteady Stretching Sheet. Energies 2020, 13, 3177. [Google Scholar] [CrossRef]
  44. Naganthran, K.; Nazar, R.; Siri, Z.; Hashim, I. Entropy Analysis and Melting Heat Transfer in the Carreau Thin Hybrid Nanofluid Film Flow. Mathematics 2021, 9, 3092. [Google Scholar] [CrossRef]
  45. Ali, R.; Shahzad, A.; us Saher, K.; Elahi, Z.; Abbas, T. The thin film flow of Al2O3 nanofluid particle over an unsteady stretching surface. Case Stud. Therm. Eng. 2022, 29, 101695. [Google Scholar] [CrossRef]
  46. Tiwari, R.K.; Das, M.K. Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf. 2007, 50, 2002–2018. [Google Scholar] [CrossRef]
  47. Alabdulhadi, S.; Waini, I.; Ahmed, S.E.; Ishak, A. Hybrid Nanofluid Flow and Heat Transfer Past an Inclined Surface. Mathematics 2021, 9, 3176. [Google Scholar] [CrossRef]
  48. Nasir, S.; Shah, Z.; Islam, S.; Bonyah, E.; Gul, T. Darcy Forchheimer nanofluid thin film flow of SWCNTs and heat transfer analysis over an unsteady stretching sheet. AIP Adv. 2019, 9, 15223. [Google Scholar] [CrossRef] [Green Version]
  49. Vajjha, R.S.; Das, D.K. Experimental determination of thermal conductivity of three nanofluids and development of new correlations. Int. J. Heat Mass Transf. 2009, 52, 4675–4682. [Google Scholar] [CrossRef]
  50. Oztop, H.F.; Abu-Nada, E. Numerical study of natural convection in partially heated rectangular enclosures filled with nanofluids. Int. J. Heat Fluid Flow 2008, 29, 1326–1336. [Google Scholar] [CrossRef]
  51. Timofeeva, E.V.; Routbort, J.L.; Singh, D. Particle shape effects on thermophysical properties of alumina nanofluids. J. Appl. Phys. 2009, 106, 14304. [Google Scholar] [CrossRef]
  52. Schlichting, H.; Gersten, K. Boundary Layer Theory; Springer: Berlin/Heidelberg, Germany, 2017. [Google Scholar]
  53. Shampine, L.F.; Gladwell, I.; Thompson, S. Solving ODEs with MATLAB; Cambridge University Press: Cambridge, UK, 2003; ISBN 9780521824040. [Google Scholar]
  54. Abel, M.S.; Mahesha, N.; Tawade, J. Heat transfer in a liquid film over an unsteady stretching surface with viscous dissipation in presence of external magnetic field. Appl. Math. Model. 2009, 33, 3430–3441. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic diagram of the present study.
Figure 1. Schematic diagram of the present study.
Mathematics 11 00739 g001
Figure 2. Impact of shape factor on velocity profile for S = 0.4, M = 1.0, K = 0.5 and φ = 0.02.
Figure 2. Impact of shape factor on velocity profile for S = 0.4, M = 1.0, K = 0.5 and φ = 0.02.
Mathematics 11 00739 g002
Figure 3. Impact of the slip parameter K on the temperature profile for S = 0.4, φ = 0.02, M = 1.0, Bi = 3.0, λ = −2.0, α = π/4 and Pr = 8.0.
Figure 3. Impact of the slip parameter K on the temperature profile for S = 0.4, φ = 0.02, M = 1.0, Bi = 3.0, λ = −2.0, α = π/4 and Pr = 8.0.
Mathematics 11 00739 g003
Figure 4. Impact of the magnetic parameter M on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, Bi = 3.0, λ = −2.0, α = π/4 and Pr = 8.0.
Figure 4. Impact of the magnetic parameter M on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, Bi = 3.0, λ = −2.0, α = π/4 and Pr = 8.0.
Mathematics 11 00739 g004aMathematics 11 00739 g004b
Figure 5. Impact of the mixed convection parameter λ on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, M = 1.0, Bi = 3.0, α = π/4 and Pr = 8.0.
Figure 5. Impact of the mixed convection parameter λ on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, M = 1.0, Bi = 3.0, α = π/4 and Pr = 8.0.
Mathematics 11 00739 g005
Figure 6. Impact of the inclination angle parameters α on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, M = 1.0, Bi = 3.0, λ = −2.0 and Pr = 8.0.
Figure 6. Impact of the inclination angle parameters α on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, M = 1.0, Bi = 3.0, λ = −2.0 and Pr = 8.0.
Mathematics 11 00739 g006
Figure 7. Impact of the Biot number Bi on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, M = 1.0, λ = −2.0, α = π/4 and Pr = 8.0.
Figure 7. Impact of the Biot number Bi on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, M = 1.0, λ = −2.0, α = π/4 and Pr = 8.0.
Mathematics 11 00739 g007
Figure 8. Impact of shape factor on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, M = 1.0, Bi = 3.0, λ = −2.0, α = π/4 and Pr = 8.0.
Figure 8. Impact of shape factor on the temperature profile for S = 0.4, K = 0.5, φ = 0.02, M = 1.0, Bi = 3.0, λ = −2.0, α = π/4 and Pr = 8.0.
Mathematics 11 00739 g008
Figure 9. The film thickness β against the unsteadiness parameter S for diverse values of M when K = 0.5 and φ = 0.02.
Figure 9. The film thickness β against the unsteadiness parameter S for diverse values of M when K = 0.5 and φ = 0.02.
Mathematics 11 00739 g009
Table 1. Thermophysical properties of Al2O3 and water.
Table 1. Thermophysical properties of Al2O3 and water.
Physical PropertiesAl2O3H2O
ρ ( k g m 3 ) 3970997.1
k ( W   m 1 K 1 ) 400.613
c p ( J   k g 1 K 1 ) 7654179
σ ( s   m 1 ) 16.55.50
β * × 10 5 ( K 1 ) 0.8521
Table 2. Viscosity and shape factor values of nanoparticles.
Table 2. Viscosity and shape factor values of nanoparticles.
Nanoparticle ShapesPlatelets
Mathematics 11 00739 i001
Cylinder
Mathematics 11 00739 i002
Brick
Mathematics 11 00739 i003
Sphere
Mathematics 11 00739 i004
Parameters
A137.113.51.92.5
A2612.6904.4471.46.5
Shape factors ( m )5.724.823.723.0
Table 3. Comparison of numerical results when M = K = φ = 0.
Table 3. Comparison of numerical results when M = K = φ = 0.
Abel et al. [54]Ali et al. [45]Present Result
Sβf″(0)βf″(0)βf″(0)
0.44.981455−1.1340984.981468−1.13409574.981466−1.1340953
0.63.131710−1.1951283.131711−1.19512523.131662−1.1951202
0.82.151990−1.2458052.152021−1.24580642.152011−1.2458089
1.01.543617−1.2777691.543615−1.27776931.543456−1.2777079
1.21.1227780−1.2791711.127779−1.27917181.127490−1.2789750
1.40.8221033−1.2335450.8210317−1.23354530.821016−1.2335280
1.60.576176−1.1149410.5761743−1.11493680.576061−1.1147574
1.80.356390−0.8674160.3563871−0.867410490.356236−0.8670686
Table 4. Numerical values of skin-friction coefficient of multi-shape nanoparticles.
Table 4. Numerical values of skin-friction coefficient of multi-shape nanoparticles.
Physical ParametersPlateletsCylinderBrickSphere
K φ M λ R e 1 / 2 C f
0.30.021.0−2.01.6812682661.4895362981.2478317611.133526940
0.51.4333466761.2586675261.0396798010.937343360
0.71.2536315221.0938017400.8938778000.802628774
0.50.021.0−2.01.4333466761.2586675261.0396798010.937343360
0.042.0890891071.8984502861.3783065370.986750837
0.062.8180981302.6941129951.8512882721.039162880
0.50.020.0−2.01.1638710571.0288904320.8596964350.778968598
0.51.3185075531.1611618070.9613120980.869720683
1.01.4333466761.2586675261.0396798010.937343360
0.50.021.0−2.01.4333466761.2586675261.0396798010.937343360
0.01.3360314661.1654732570.9520853410.853015117
2.01.2438477951.0780978240.8714607770.776276122
Table 5. Numerical values of Nusselt number for multi-shape nanoparticles.
Table 5. Numerical values of Nusselt number for multi-shape nanoparticles.
Physical ParametersPlateletsCylinderBrickSphere
B i α K M λ R e 1 / 2 N u
2.0 π / 4 0.51.0−2.01.4711844041.4334862921.3842835011.353043680
2.51.6706693801.6281712111.5721617751.536832098
3.01.8364417221.7899560211.7281639601.689397686
3.0 π / 6 0.51.0−2.01.8333108381.7862653591.7234156331.683998156
π / 4 1.8364417221.7899560211.7281639601.689397686
π / 3 1.8404193121.7946186921.7340992221.696099546
3.0 π / 4 0.31.0−2.01.8783817031.8333163581.7735331271.735593784
0.51.8364417221.7899560211.7281639601.689397686
0.71.8020961371.7548179891.6922742491.652735779
3.0 π / 4 0.50.0−2.01.8752384101.8301205811.7701306601.732319575
0.51.8533059271.8073473471.7468796531.708144155
1.01.8364417221.7899560211.7281639601.689397686
3.0 π / 4 0.51.0−2.01.8364417221.7899560211.7281639601.689397686
0.01.8495845491.8052601391.7474069421.710956457
2.01.8616246341.8190469551.7642273141.729451205
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alabdulhadi, S.; Abu Bakar, S.; Ishak, A.; Waini, I.; Ahmed, S.E. Effect of Buoyancy Force on an Unsteady Thin Film Flow of Al2O3/Water Nanofluid over an Inclined Stretching Sheet. Mathematics 2023, 11, 739. https://doi.org/10.3390/math11030739

AMA Style

Alabdulhadi S, Abu Bakar S, Ishak A, Waini I, Ahmed SE. Effect of Buoyancy Force on an Unsteady Thin Film Flow of Al2O3/Water Nanofluid over an Inclined Stretching Sheet. Mathematics. 2023; 11(3):739. https://doi.org/10.3390/math11030739

Chicago/Turabian Style

Alabdulhadi, Sumayyah, Sakhinah Abu Bakar, Anuar Ishak, Iskandar Waini, and Sameh E. Ahmed. 2023. "Effect of Buoyancy Force on an Unsteady Thin Film Flow of Al2O3/Water Nanofluid over an Inclined Stretching Sheet" Mathematics 11, no. 3: 739. https://doi.org/10.3390/math11030739

APA Style

Alabdulhadi, S., Abu Bakar, S., Ishak, A., Waini, I., & Ahmed, S. E. (2023). Effect of Buoyancy Force on an Unsteady Thin Film Flow of Al2O3/Water Nanofluid over an Inclined Stretching Sheet. Mathematics, 11(3), 739. https://doi.org/10.3390/math11030739

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