Nanofluids modeling and energy transformation under the influence of mixed convection: Role of aluminum and γ-aluminum nanoparticles

Introduction

A lot of researchers studied temperature behavior in the base fluids such as ethylene glycol, kerosene oil, and H 2 O . These fluids have limited applications due to their low thermal performance, as for many industrial uses high thermal transport is required to finish the processes of different products. Therefore, researchers focused on the development of a new class of fluids with high thermal performance than the regular liquids. These fluids gained much intentions of the researches and industrialists and strengthen their roots in the modern world.

Choi ¹ potentially contributed in the development of nanofluids. Finally, his untiring efforts opened a new lock in the modern world. They explored that thermal performance in the regular liquids could be intensify by adding nano additives of various metals and oxides. These particles are saturated in the liquid uniformly and thermally in equilibrium. The addition of nanoparticles in the regular liquid boosts thermal performance of the fluid. These became very popular among the researcher’s community and resolved many problems of the industrialist which they faced from many decades in the context of heat transport. Therefore, the investigation of thermal transport in nanofluids prepared by various metals oxides by dispersing in various regular liquids are reported by Ashraf et al.^2,3and Chu et al. ⁴

The study related to nanofluids containing Al 2 O 3 and γ − Al 2 O 3 nanomaterials is described by Maciver et al. ⁵ The effective thermal conductivity of nanofluids is imperative and playing significant role to boosts the internal energy of the nanofluids which ultimately results in rise in the fluid temperature. Therefore, numerous thermal conductance correlations were introduced and to analyze their thermal capability, researchers and engineers utilized those correlations in the modeling of nanofluids and examined their influences in the fluid thermal transport ability. The fluid dynamists made their efforts day by day and reported the nanofluid temperature behavior prepared by various metals and regular liquids. The significant studies in this regard under the impacts of various physical quantities like viscus dissipation, thermal radiations, combined convection effects, nonlinear thermal radiations, chemical reaction, and internal heat generation and absorption are reported in Hamad, ⁶ Vajravelu etal., ⁷ and Ashraf et al. ⁸ Synthetization process of nanofluids given in Figure 1(a).

OPEN IN VIEWER

Figure 1.

Synthetization of nanofluid in (a) and physical description of nanofluids flow in (b).

Thermal conductivity and thickness reproductions are expenditure for Al 2 O 3 . Individually investigative in addition mathematical investigation are supported available. The flow above a stretching surface is an important problem in lots of industrial procedures with uses in manufacturing such as melt-spinning, extrusion, the boiling developing, wire illustration, glass–fiber manufacture, rubber sheets and refrigeration of a big metal plate in immersion, which may be an electrolyte, etc. The production, polymer pieces, production of plastic, and fibers are industrial in unbroken extrusion of polymer commencing die to a mechanical wave, which is placed at a limited coldness. Considered the fixed 2D incompressible flow of a Newtonian fluid instigated using the widening of an elastic smooth sheet which transfers in its specific horizontal velocity changing linearly as the distance from a steady idea due to the use of a constant pressure. Problematic is mostly stimulating since a particular result of the 2D Navier–Stokes equations has been developed by Tiwari and Das. ⁹

Later on, Crane ¹⁰ considered an effective substance to study the heat transfer mechanism and unsteady flow of nanofluid past a stretching sheet. They modeled the problem for nanofluid and examined the behavior of thermal transport under the impacts of various flow governing quantities. The important investigations of nanofluids upshot temperature under physical phenomena like thermal radiations and convective conditions etc. reports in Ashraf et al.^11,12 The increase of the unusual kind of flow due to reduction was first explored by Wang ¹³ and considered the performance of fluid over an unsteady stretching sheet. They studied the mixed convection heat transfer on a non-linear stretching geometry with inconstant fluid done by Prasad et al. ¹⁴ Yazdi et al. ¹⁵ conversed the slip flow and heat transfer over a non-linear porous stretching surface. The investigation of the nanofluid dynamics saturated by aluminum alloys past a stretching surface was described in Zaidi et al. ¹⁶ Recently, Rana and Bhargava ¹⁷ studied the boundary layer flow of nanofluid flow on a non-linearly stretching sheet. Just, Das ¹⁸ considered heat transfer convective of nanofluids over a stretching sheet in the existence of partial slip-up and thermal radioactivity. First idea on show in their investigated that the stretching of the sheet may not automatically linear by Gupta and Gupta. ¹⁹

Thermal enhancement and entropy^{20 –23} study in nanofluids convinced by thermal radiations behavior made place in the scientist’s community. The forced scientists from various applied areas of research to investigate and explore fascinating insights in thermal and enhancement studies. Carbon nanotubes which are found to be very useful for enhanced heat transport due to their unique structure compelled the researchers to examine the dynamics of nanofluids. The Fluid Dynamists made several attempts to study the role of CNTs in the heat as well as entropy of the fluid. The studies conducted by adding variety of body forces and under various conditions specified at the boundary domain, such studies available in the literature.^{24 –27} Among the distinct heat transport mode, convection with addition of CCHFM (Cattaneo Christov Heat Flux Model) is one of the prominent sources to upshot the heat transfer inside the working nanofluid. Recently, many studies^{28 –31} reported in this regard by taking various sorts of working nanofluids synthesized by different metals and host solvents. The literature revealed that inside temperature of the nanofluids could be enriched by inducing different physical phenomena in the modeling of the heat transfer problem. Further, many fascinating and useful results presented in the studies^{32 –35} and found excellent heat transport which would be helpful to the industrialist and engineers to cope their issues in the current time.

From the above literature, it is found that the novel comparative analysis of heat transfer in Al₂O₃-H₂O, γAl₂O₃-H₂O, Al₂O₃-C₂H₆O₂, and γAl₂O₃- C₂H₆O₂ nanofluids is absent which is significant research gap. The γAl₂O₃ nanoparticles are very famous due to their progressive applications in the world specifically in biomedical engineering. Therefore, in this study the behavior of Al₂O₃-H₂O, γAl₂O₃-H₂O, Al₂O₃-C₂H₆O₂, and γAl₂O₃-C₂H₆O₂ heat transfer with outstanding thermophysical attributes of two base solvents and updated nanofluid models for γAl₂O₃ is presented. Further, the model will be upgraded with new insights of combined convection which is an important factor in the analysis of nanofluids. This study will fulfill the above gap and help the researchers in the selection of parameters for better heat transfer.

The study will subject to the following limitations:

Formulation of the problem

Considered an unsteady incompressible flow of nanofluids composed by Al 2 O 3 , γ − Al 2 O 3 nanomaterials and the regular liquids namely C 2 H 6 O 2 and H 2 O above a stretching sheet. The flow is configurated due to a stretching surface in two-dimensions. The velocity of the sheet is u w = a x along the x-axis and temperature at the stretchable surface is T w = T ∞ + b x 2 . Where a, b are constants. The configuration of nanofluid flow over a stretchable surface is depicted in Figure 1.

The basic physical laws that govern the flow of Al₂O₃-H₂O, γAl₂O₃-H₂O, Al₂O₃-C₂H₆O₂ and γAl₂O₃-C₂H₆O₂ over unsteady stretching surface are described as below in vector form ³⁶ :

∇ . V = 0 ρ n f [ ∂ V ∂ t + ( V . ∇ ) V ] = μ n f [ ∇ 2 V ] + S 1 ( ρ c p ) n f [ ∂ T ∂ t + ( V . ∇ ) T ] = k n f ∇ 2 T }

Where; S 1 incorporates the effects of combined convection.

∂ u ∂ x + ∂ v ∂ y = 0

(1)

ρ n f ( ∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y ) = μ n f ∂ 2 u ∂ y 2 + g ( β ρ ) n f ( T − T ∞ ) .

(2)

( ρ C p ) n f ( ∂ T ∂ t + u ∂ T ∂ x + v ∂ T ∂ y ) = K n f ∂ 2 T ∂ y 2

(3)

The flow is subject to the following flow conditions:

at y = 0 :

U = U w , v = 0 , T w = T ∞ + b x 2 ,

(4)

as y → ∞ :

U → 0 , T = T ∞ .

(5)

Thermal physical values

The effective correlations which work particularly for γ A l 2 O 3 nanoparticles and by considering the base solvents are listed in equations (7)–(12). These properties of the nanofluids will be utilized in the development of Al₂O₃-H₂O, γAl₂O₃-H₂O, Al₂O₃-C₂H₆O₂, and γAl₂O₃-C₂H₆O₂ thermal transport models after incorporating in the basic governing laws ³⁷ :

ρ n f = ( 1 − ϕ ) ρ f + ϕ ρ s ( ρ C p ) n f = ( 1 − ϕ ) ( ρ C p ) f + ϕ ( ρ C p ) s ( ρ β ) n f = ( 1 − ϕ ) ( ρ β p ) f + ϕ ( ρ β p ) s }

(6)

Wherever ϕ represent the nanoparticles volumetric fraction.

μ n f μ f = 1 ( 1 − ϕ ) 2 . 5 . ( f o r A l 2 O 3 − H 2 O , and A l 2 O 3 − C 2 H 6 O 2 ) .

(7)

μ n f μ f = 306 ϕ 2 − 0 . 19 ϕ + 1 . ( f o r γ − A l 2 O 3 − C 2 H 6 O 3

(8)

μ n f μ f = 123 ϕ 2 − 7 . 3 ϕ + 1 , ( f o r γ − A l 2 O 3 − H 2 O )

(9)

Thermal conductivity of the nanofluid is given by the relations:

K n f = K f [ K s + 2 K f − 2 ϕ ( K f − K s ) ] [ K s + 2 K f + ϕ ( K f − K s ) ] , ( f o r A l 2 O 3 − C 2 H 6 O 2 , a n d A l 2 O 3 − C 2 H 6 O 2 ) .

(10)

K n f = K f [ 28 . 905 ϕ 2 + 2 . 8273 ϕ + 1 ] , ( f o r γ − A l 2 O 3 − C 2 H 6 O 2 ) .

(11)

K n f = K f [ 4 . 97 ϕ 2 + 2 . 72 ϕ + 1 ] , ( f o r γ − A l 2 O 3 − H 2 O ) .

(12)

Now, we introduce the following dimensionless functions F and θ , and the similarity variable η as:

η = a υ f ( 1 − c t ) y ψ = a υ f ( 1 − c t ) x F ( η ) θ ( η ) = T − T ∞ T w − T ∞ }

(13)

Here, ψ is the stream function, u and v are the velocity components and defined as:

u = ∂ ψ ∂ y , v = − ∂ ψ ∂ x .

(14)

The transformed momentum and energy equations combine with the boundary conditions are given by equations (2)–(4), can be written as:

F ‴ ( η ) + ( 1 − ϕ + ϕ ρ s ρ f ) ( 1 − ϕ ) 2 . 5 { F ( η ) F ″ ( η ) − ( F ′ ( η ) ) 2 − A ( F ′ ( η ) + 1 2 η F ″ ( η ) ) } + ( 1 − ϕ ) 2 . 5 ( 1 − ϕ + ϕ ( ρ β ) s ( ρ β ) f ) ( γ θ ( η ) ) = 0 , ( for A l 2 O 3 − H 2 O a n d A l 2 O 3 − C 2 H 6 O 2 ) .

(15)

F ‴ ( η ) + ( 1 − ϕ + ϕ ρ s ρ f ) ( 306 ϕ 2 − 0 . 19 ϕ + 1 ) { F ( η ) F ″ ( η ) − ( F ′ ( η ) ) 2 − A ( F ′ ( η ) + 1 2 η F ″ ( η ) ) } + ( 1 − ϕ + ϕ ( ρ β ) s ( ρ β ) f ) ( 306 ϕ 2 − 0 . 19 ϕ + 1 ) ( γ θ ( η ) ) = 0 , ( f o r γ − A l 2 O 3 − C 2 H 6 O 2 ) .

(16)

F ‴ ( η ) + ( 1 − ϕ + ϕ ρ s ρ f ) ( 123 ϕ 2 − 7 . 3 ϕ + 1 ) { F ( η ) F ″ ( η ) − ( F ′ ( η ) ) 2 − A ( F ′ ( η ) + 1 2 η F ″ ( η ) ) } + ( 1 − ϕ + ϕ ( ρ β ) s ( ρ β ) f ) ( 123 ϕ 2 − 7 . 3 ϕ + 1 ) ( γ θ ( η ) ) = 0 . ( f o r γ − A l 2 O 3 − H 2 O ) .

(17)

The boundary conditions are:

F ( 0 ) = 0 , F ′ ( 0 ) = 1 , a n d F ′ ( ∞ ) = 0 , θ ( 0 ) = 1 , θ ( ∞ ) = 0 .

(18)

So, the energy equation becomes:

K s − 2 K f + 2 ϕ ( K f − K s ) K s − 2 K f + ϕ ( K f − K s ) P r ( 1 − ϕ + ϕ ( ρ C p ) s ( ρ C p ) f ) θ ″ ( η ) − A 1 2 η θ ′ ( η ) + θ ′ ( η ) F ( η ) = 0 , ( f o r A l 2 O 3 − H 2 O a n d A l 2 O 3 − C 2 H 6 O 2 ) .

(19)

1 P r ( ( 1 − ϕ ) + ϕ ( ρ C P ) s ( ρ C p ) f ) ( 4 . 97 ϕ 2 + 2 . 72 ϕ + 1 ) θ ″ ( η ) − 1 2 A η θ ′ ( η ) + F ( η ) θ ′ ( η ) = 0 , ( f o r γ − A l 2 O 3 − H 2 O ) 1 P r ( ( 1 − ϕ ) ( ρ C p ) f + ϕ ( ρ C P ) s ) ( 28 . 905 ϕ 2 + 2 . 8273 ϕ + 1 ) θ ″ ( η ) − 1 2 A η θ ′ ( η ) + F ( η ) θ ′ ( η ) = 0 ( f o r γ − A l 2 O 3 − C 2 H 6 O 2 )

(20)

The dimensionless parameters appeared in the final model due to the induction of combined convection and unsteady effects are A = c a (unsteady parameter), P r = ( ν f ) ( ρ c p ) f k f (Prandtl number) and γ = g ( ρ β ) f ρ f b a 2 (buoyancy parameter).

Procedure for solution

The numerical technique based on shooting algorithm is adopted to tackle the model. This technique usually based on a set of first order ordinary differential equations which are obtained from the nonlinear flow model by using some sort of substitutions. For under consideration model, these substitutions for the velocity and temperature equations are made in the following way and the method implementation is also given in the chart:

Further,

( y 1 , y 2 , y 3 , y 3 ′ , y 5 , y 6 , y 6 ′ ) = ( F , F ′ , F ″ , F ‴ , θ , θ ′ , θ ″ ) y ′ 3 = [ ( 1 − ϕ ) 2 . 5 { ( 1 − ϕ + ϕ ρ s ρ f ) ( A ( y 1 + 1 2 η y 3 ) + y 2 2 − y 1 y 3 ) − ( 1 − ϕ + ϕ ( ρ β ) s ( ρ β ) f ) γ y 5 } ] , ( f o r A l 2 O 3 − H 2 O a n d A l 2 O 3 − C 2 H 6 O 2 )

(21)

y ′ 3 = [ 1 ( 306 ϕ 2 − 0 . 19 ϕ + 1 ) { ( 1 − ϕ + ϕ ρ s ρ f ) ( A ( y 1 + 1 2 η y 3 ) + y 2 2 − y 1 y 3 ) − ( 1 − ϕ + ϕ ( ρ β ) s ( ρ β ) f ) γ y 5 } ] , ( f o r γ A l 2 O 3 − C 2 H 6 O 2 )

(22)

y ′ 3 = [ 1 ( 123 ϕ 2 − 7 . 3 ϕ + 1 ) { ( 1 − ϕ + ϕ ρ s ρ f ) ( A ( y 1 + 1 2 η y 3 ) + y 2 2 − y 1 y 3 ) − ( 1 − ϕ + ϕ ( ρ β ) s ( ρ β ) f ) γ y 5 } ] , ( f o r γ A l 2 O 3 − H 2 O )

(23)

y 6 ′ = − ( 1 − ϕ + ϕ ( ρ C p ) s ( ρ C p ) f ) K s − 2 K f + 2 ϕ ( K f − K s ) K s − 2 K f + ϕ ( K f − K s ) P r [ y 1 y 6 − 1 2 A η y 6 ] , ( f o r A l 2 O 3 − H 2 O a n d A l 2 O 3 − C 2 H 6 O 2 ) .

(24)

y 6 ′ = − ( 1 − ϕ + ϕ ( ρ C p ) s ( ρ C p ) f ) ( 4 . 97 ϕ 2 + 2 . 72 ϕ + 1 ) P r [ y 1 y 6 − 1 2 A η y 6 ] . ( f o r γ A l 2 O 3 − H 2 O )

(25)

y 6 ′ = − ( 1 − ϕ + ϕ ( ρ C p ) s ( ρ C p ) f ) ( 4 . 97 ϕ 2 + 2 . 72 ϕ + 1 ) P r [ y 1 y 6 − 1 2 A η y 6 ] . ( f o r γ A l 2 O 3 − H 2 O )

(26)

y 6 ′ = − ( 1 − ϕ + ϕ ( ρ C p ) s ( ρ C p ) f ) ( 28 . 905 ϕ 2 + 2 . 8273 ϕ + 1 ) P r [ y 1 y 6 − 1 2 A η y 6 ] . ( f o r γ A l 2 O 3 − C 2 H 6 O 2 )

(27)

Quantities of engineering concern

The analysis of skin friction coefficient and local heat transport mechanism is significant from industrial and technological point of view. Further, thermophysical values described in Table 1. For under consideration model, the dimensional mathematical relations define as:

C f = − 2 μ n f ( ∂ u ∂ y ) y = 0 ρ f u w 2 a n d N u x = x q w K f ( T w − T ∞ )

(28)

Where q w is the heat transfer from sheet and can be written as:

q w = − K n f ( ∂ T ∂ y ) y = 0

(29)

OPEN IN VIEWER

Table 1.

Thermo-physical properties.^38,39.

Properties	Pure water H 2 O	Al 2 O 3	Cu	C 2 H 6 O 2
ρ ( k g / m 3 )	997 . 1	3970	8933	1116 . 6
C p ( J / k g k )	4179	765	385	2382
k ( W m − 1 k − 1 )	0 . 6130	40	401	0 . 249
β × 10 5 k − 1	21	0.850	1 . 670	6.5
P r	6.96	-----	-----	204

Applying the dimensional transformation (12) we obtain:

1 2 R e x 1 2 C f = − ( 1 − ϕ ) − 2 . 5 F ″ ( 0 ) . ( f o r A l 2 O 3 − H 2 O a n d A l 2 O 3 − C 2 H 6 O 2 )

(30)

1 2 R e x 1 2 C f = − ( 123 ϕ 2 + 7 . 3 ϕ + 1 ) F ″ ( 0 ) . ( f o r γ A l 2 O 3 − H 2 O )

(31)

1 2 R e x 1 2 C f = − ( 306 ϕ 2 − 0 . 19 ϕ + 1 ) F ‴ ( 0 ) . ( f o r γ A l 2 O 3 − C 2 H 6 O 2 )

(32)

Where 1 2 R e x − 1 2 N u x is the reduced Nusselt number:

1 2 R e x − 1 2 N u x = − ( K s + 2 K f − 2 ϕ ( K f − K s ) K s + 2 K f + ϕ ( K f − K s ) ) θ ′ ( 0 ) . ( f o r A l 2 O 3 − C 2 H 6 O 2 )

(33)

1 2 R e x − 1 2 N u x = − ( K s + 2 K f − 2 ϕ ( K f − K s ) K s + 2 K f + ϕ ( K f − K s ) ) θ ′ ( 0 ) ( f o r A l 2 O 3 − H 2 O a n d A l 2 O 3 − C 2 H 6 O 2 )

(34)

1 2 R e x − 1 2 N u x = − ( 123 ϕ 2 + 7 . 3 ϕ + 1 ) θ ′ ( 0 ) ( f o r γ A l 2 O 3 − H 2 O )

(35)

1 2 R e x − 1 2 N u x = − ( 4 . 9 ϕ 2 + 2 . 72 ϕ + 1 ) θ ′ ( 0 ) , ( f o r γ A l 2 O 3 − C 2 H 6 O 2 )

(36)

Results and discussion

In the study of nanofluids, the behavior of working fluid regarding the momentum and energy transport is of huge significance which make them prominent from the common fluids. In this research, variety of nanoliquids (Al₂O₃-H₂O, γAl₂O₃-H₂O, Al₂O₃-C₂H₆O₂, and γAl₂O₃-C₂H₆O₂) which prepared by two distinct base solvents (H₂O, C₂H₆O₂) and important nanomaterial (Al₂O₃, γAl₂O₃) are taken in the existence of combined convection and unsteady effects. These nanomaterials became much popular and can be utilized significantly in industrial applications due to their superior chemical properties. The comparative alterations in the velocity ( F ′ ( η ) ) and temperature ( θ ( η ) ) profiles under varying flow quantities are demonstrated in this section.

The velocity F ′ ( η ) variations for Al₂O₃-H₂O, γAl₂O₃-H₂O, due to unsteady effects A are captured in Figure 2(a) and (b). Keen observations of the plotted results reveal that the motion of the both nanoliquids decays under stronger unsteady effects. The movement of nanoliquids based on simple Al₂O₃ promptly decays than γAl₂O_3. Physically, the nanoliquids Al₂O₃-H₂O and Al₂O₃-C₂H₆O₂ are more-denser and in the view of their thermophysical characteristics which opposes the motion. The viscosity inside the nanoliquids intensifies which resists the particles motion. Figure 2(c) to (f) demonstrating thermal distribution (θ (η)) for the nanofluids against increasing unsteady effects A . These results reveal that nanoliquids are better for thermal transport in which thermal conductivity of Al₂O₃ and γAl₂O₃ are have important role.

OPEN IN VIEWER

Figure 2.

Variations in the velocity field against A (a and b) and in temperature via A in (c–f).

The velocity trends against increasing ϕ are depicted in Figure 3(a) and (b) for Al₂O₃-H₂O, γAl₂O₃-H₂O, Al₂O₃-C₂H₆O₂, and γAl₂O₃-C₂H₆O₂ nanoliquids, respectively. Physically, increment in the velocity occurs because for higher fraction ϕ , viscosity and density are decreases and the fluid particles move freely within the domain. Figure 3(c) to (f) furnished to investigate the temperature behavior by enhancing volume concentration. Thermal conductivity of nanoparticles playing significant role in the temperature enhancing of the nanoliquids. Physically, additional thermal conductivity of Al₂O₃ γAl₂O₃ is a key tool for thermal enhancement because these are better conductor and intensify the heat transport inside the fluid.

OPEN IN VIEWER

Figure 3.

The velocity via ϕ (a and b) and temperature via ϕ (c–f).

Figure 4(a) and (b) furnished for the analysis of F ′ ( η ) due to increasing buoyancy effects. It is noticed that the velocity grows almost inconsequentially by strengthen γ . Further, momentum boundary layer portion rises for γAl₂O₃-H₂O and γAl₂O₃-C₂H₆O₂ and it decays for Al₂O₃-H₂O and Al₂O₃-C₂H₆O₂. Figure 4(c) to (f) demonstrating impacts of imposed mixed convection effects ( γ ) on the temperature distribution θ ( η ) . These results showing that the temperature drops due to upsurges in γ . Physically, the motion declines due to stronger γ as a result the particles collide slowly and friction between the successive layer of the fluid reduces. This reduction in collision affects the temperature and hence it increases very slowly.

OPEN IN VIEWER

Figure 4.

The velocity via γ (a and b) and temperature via γ (c–f).

Conclusions

The flow of nanofluids comprising γ Al 2 O 3 , Al 2 O 3 nanomaterials is conducted under unsteady effects. The model is handled numerically and demonstrated the results against the significant physical quantities. From the graphical representation of the results, it is examined that the nanofluid motion declines abruptly for the nanofluid composed by Al 2 O 3 than γ Al 2 O 3 . The volumetric fraction of the metallic particles is a key ingredient for thermal enhancement in the nanofluid and showing dominant behavior for γ -nanomaterials based nanofluid. Dual behavior of the temperature is investigated against combined convection effects. Finally, it is concluded that the nanofluids are reliable for thermal enhancement and beneficial for many industrial and technological purposes as compared to that of base liquids which have low thermal performance characteristics.