Next Article in Journal
Photosynthetic Efficiency is Higher in Asymmetric Leaves than in Symmetric Leaves of the Same Plant
Next Article in Special Issue
Peristaltic Pumping of Nanofluids through a Tapered Channel in a Porous Environment: Applications in Blood Flow
Previous Article in Journal
Online Visual Tracking of Weighted Multiple Instance Learning via Neutrosophic Similarity-Based Objectness Estimation
Previous Article in Special Issue
Heat Transfer of Oil/MWCNT Nanofluid Jet Injection Inside a Rectangular Microchannel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modified MHD Radiative Mixed Convective Nanofluid Flow Model with Consideration of the Impact of Freezing Temperature and Molecular Diameter

1
Department of Mathematics and Statistics, Hazara University, Mansehra 21120, Pakistan
2
Department of Mathematics, Mohi-ud-Din Islamic University Nerian Sharif, Azad Jammu & Kashmir 12080, Pakistan
3
Department of Mathematics, Faculty of Sciences, HITEC University Taxila Cantt, Punjab 47080, Pakistan
4
Department of Mathematics, College of Science Al-Zulfi, Majmaah University, Al-Majmaah 11952, Saudi Arabia
5
Department of Mathematics, COMSATS University Islamabad, Abbottabad 22010, Pakistan
6
Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City 72915, Vietnam
7
Department of Mechanical Engineering, College of Engineering, Prince Mohammad Bin Fahd University, Al Khobar 31952, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(6), 833; https://doi.org/10.3390/sym11060833
Submission received: 13 March 2019 / Revised: 9 April 2019 / Accepted: 16 April 2019 / Published: 25 June 2019
(This article belongs to the Special Issue Symmetry and Fluid Mechanics)

Abstract

:
Magnetohydrodynamics (MHD) deals with the analysis of electrically conducting fluids. The study of nanofluids by considering the influence of MHD phenomena is a topic of great interest from an industrial and technological point of view. Thus, the modified MHD mixed convective, nonlinear, radiative and dissipative problem was modelled over an arc-shaped geometry for Al2O3 + H2O nanofluid at 310 K and the freezing temperature of 273.15 K. Firstly, the model was reduced into a coupled set of ordinary differential equations using similarity transformations. The impact of the freezing temperature and the molecular diameter were incorporated in the energy equation. Then, the Runge–Kutta scheme, along with the shooting technique, was adopted for the mathematical computations and code was written in Mathematica 10.0. Further, a comprehensive discussion of the flow characteristics is provided. The results for the dynamic viscosity, heat capacity and effective density of the nanoparticles were examined for various nanoparticle diameters and volume fractions.

1. Introduction

The liquids regularly used in heat transfer applications, such as water, propylene glycol, ethylene glycol, kerosene oil, engine oil and transformer oil, are extensively used in industry and in thermal power plants. Due to their reduced thermal conductivity, these liquids do not have effective heat transfer characteristics. However, for a great deal of industrial production, remarkable amounts of heat are required. The thermal conductivity of the solid materials, such as different metals and oxides, is very high in comparison with regular liquids. Thus, scientists and engineers hypothesized that the heat transfer in working fluids could be enhanced by mixing in the nanoparticles in above regular liquids. Finally, Choi [1] developed a colloidal composition that has effective heat transfer characteristics, as compared to regular fluids. Choi [1] unlocked a new innovative research area and researchers, scientists and engineers focused on the analysis of nanofluids. Before the development of nanofluids, heat transfer was a major problem from an industrial point of view, since considerable amounts of heat transfer were required for a great deal of technological and industrial production, and regular liquids failed to provide the desired amount of heat. Thermal conductivity plays a major role in the heat transfer rate of nanofluids. Thus, several theoretical models were proposed for thermal conductivity. The thermal conductivity model was developed by Maxwell in 1873 [2], who considered nanosized particles, and this can be considered the origin of the concept of the nanofluids.
Several theoretical models based on nanoparticle characteristics that take into account the effects of temperature and the shape and diameter of nanoparticles, as well as Brownian motion, have been presented. In 1935, Bruggemann [3] constructed a thermal conductivity correlation for spherical nanoparticles that was limited to high concentration patterns. The behavior of thermal conductivity was developed by Hamilton [4] in 1962, who explored the effects of nanoparticle shape. In 1996, Lu and Lin [5] proposed a model incorporating the effects of Brownian dynamics. The thermal conductivity model for the interaction between the nanoparticles and their surrounding liquid was developed by Koo and Kleinstreuer [6,7]. Xue [8] developed a thermal conductivity model for carbon nanotubes. In 2005, Prasher et al. [9] found a correlation by considering the influence of convection on the thermal conductivity of nanoparticles. In 2006, Li [10] made apparent the influence of temperature on thermal conductivity and outlined a correlation for Al2O3/H2O and CuO/H2O nanofluids. In 2011, Corcione [11] developed a model for Al2O3 + H2O nanofluids by incorporating the effects of freezing temperatures. By incorporating thermal conductivity models, the above researchers presented various models and described the heat transfer enhancement due to thermal conductivity. Some useful studies for nanofluids are described in [12,13,14,15,16].
Similar to nanoparticles, carbon nanotubes also have high thermal conductivity and unique mechanical and chemical properties. Carbon nanotubes are subcategorized as either single or multiple walled carbon nanotubes. The concept of colloidal suspension in relation to carbon nanotubes was presented by Iijima [17]. After the development of the thermal conductivity correlation for carbon nanotubes, many studies were presented outlining thermal enhancement due to suspended carbon nanotubes. Recently, Ahmed et al. [18] explored the influence of thermal radiation and viscous dissipation on the flow of water suspended by carbon nanotubes. The effect of thermophysical characteristics of the nanotubes on the heat transfer enhancement water over a curved surface and non-parallel walls was described in [19,20], respectively.
Recently, the flow over an arc-shaped geometry has become a point of interest. Reddy et al. [21] recently modified the curve-shaped flow model for nonlinear radiative heat flux. They also examined the impact of the cross-diffusion phenomenon on heat and mass transfer. Another useful mechanism to enhance the fluid temperature is ohmic heating, which produces extra heat in the conductor, with electrons supplying energy to the atoms of the conductor through collisions. In 2017, Hayat et al. [22] examined the effects of resistive heating on the curve surface flow.
A literature review revealed that there have not yet been any studies of the impact of freezing temperatures and the diameter of nanoparticles on the flow of incompressible fluids due to the effects of nonlinear radiative heat flux, viscous dissipation, mixed convection and Lorentz forces. This study is presented to cover this significant gap. The nanofluids Al2O3 + H2O were used to study the characteristics of the flow and other effective thermophysical properties, such as effective density, heat capacity and thermal conductivity. The results for shear stress and local heat transfer are also described and discussed comprehensively. Finally, major findings of the study is presented.

2. Model Formulation

We considered the laminar time independent and the incompressible flow of the Al2O3 + H2O nanofluid by taking into account the influence of a nonlinear, radiative heat flux and the imposed variable magnetic field over an arc geometry situated in the curvilinear frame r and s. Further, the r-axis was perpendicular, and the arc was placed in the direction of s. The velocity and magnetic field were functions of s and mathematically described as below:
u w ^ ( s ) = b s m , B ^ ( s ) = B 0 ( s 0.5 ( m 1 ) ) 1 ,
It was assumed that the induced magnetic field was inconsequential, and therefore was not taken into consideration. The temperatures at the arched and the free surface were T ^ w and T ^ , respectively. The value of m > 1 represented the flow over an arc shaped, which was nonlinearly stretched. m = 1 was for the flow of a linearly stretching geometry. Figure 1 presents the flow description in a curvilinear frame. The flow chart of the study presented in Figure 2.
The Al2O3 + H2O nanofluid flow, incorporating the phenomena of Lorentz forces, viscous dissipation and nonlinear radiative heat flux, is described by the following system [21]:
r ( v   ^ ( r + a ) ) = 0 ,
u   ^ 2 = p r ( r + a ) ( ρ n f ) 1 ,
ρ n f ( v   ^ u   ^ r + ( a r + a ) u   ^ u   ^ s + u   ^ v ^ ( r + a ) ) = a ( r + a ) p s + μ n f ( 2 u   ^ r 2 + ( 1 r + a ) u   ^ r ( 1 r + a ) 2 u   ^   ) σ n f B 0 2 u   ^ + g ( ρ β ) n f ( T T ) ,
a ( r + a ) u   ^ T   ^ s + v   ^ T   ^ r = k n f ( ρ C p ) n f ( 1 ( r + a ) T ^ r + 2 T   ^ r 2 ) + μ n f ( ρ C p ) n f ( u   ^ r u   ^ ( r + a ) ) 2 1 ( ρ C p ) n f ( r + a ) r ( q R ( r + a ) ) ,
with boundary conditions at the curve and far from the curve of:
At r = 0 :
u   ^ = b s m v   ^ = 0 T ^ = T w } .
At r :
u   ^ 0 u   ^ r 0 T   ^ T } .
The expression for nonlinear radiative heat flux is described as:
q R = 4 σ 3 k ^ ^   r ( T ^ 4 ) = 16 3 σ k ^ ^   T   ^ 3 T ^   r ,
Equations (1)–(4) are the conservation of mass, momentum and energy, respectively. Further, the Stefan–Boltzmann law and adsorption coefficient are σ ^   and k   ^ , respectively. To enhance the thermal and physical characteristics of the model, the following effective models for thermophysical characteristics were used [23]:
ρ n f = [ ( 1 ϕ ) + ϕ ρ p ρ f ] ρ f ,
( ρ C p ) n f = [ ( 1 ϕ ) + ϕ ( ρ C p ) p ( ρ C p ) f ] ( ρ C p ) f ,
μ n f = μ f ( 1 34.87 ( d p a r t i c l e d f l u i d ) 0.3 ϕ 1.03 ) 1 .
k n f = k f ( 1 + 4.4 R e b 0.4 Pr 0.66 ( T T f r e e z i n g ) 10 ( k p k f ) 0.03 ϕ 0.66 ) ,
( ρ β ) n f = [ ( 1 ϕ ) + ϕ ( ρ β ) p ( ρ β ) f ] ( ρ β ) f ,
σ n f = σ f [ 1 + 3 ( Θ 1 ) ϕ ( Θ + 2 ) ( Θ 1 ) ϕ ] ,   where ,   Θ = σ p σ f .
In Equation (11), R e b shows the Reynolds number due to Brownian motion, and is described in the following pattern:
R e b ( μ f ) = d p ρ f u b ,
In Equation (14), the velocity of Brownian motion is calculated by the formula:
u b = 2 T k b ( π d p 2 μ f ) ,
Here, k b is the Stefan–Boltzmann coefficient and its value is 1.380648 × 10 23 (JK−1). d p represents the molecular diameter which is calculated by the expression [24]:
d f = 6 M * ( N * ρ f π ) 1 ,
The molecular weight of the regular liquid and Avogadro number are denoted by M * and N * , respectively. Further, the value of d f is calculated as:
d f = ( 6 × 0.01801528 998.62 × ( 6.022 × 10 23 ) × π ) 1 3 = 3.85 × 10 10   nm ,
Table 1 shows the data for thermal conductivity, effective density, thermal expansion coefficient, and effective dynamic viscosity [23]:
The similarity transformations are described in the following set of equations [21]:
η = b ν b f s 0.5 ( m 1 ) r u   ^ = b s m L p = ρ b f b 2 s 2 m P v   ^ = a ( r + a )   b ν b f s 0.5 ( m 1 ) { 0.5 ( m + 1 ) L + 0.5 ( m 1 ) η L } N = T T T w T }
To analyze the phenomena of nonlinear radiative heat flux, the following expression was used:
T   ^ = T ( 1 + ( β w 1 ) N )
The ratio of wall and free surface temperature was denoted by β w .
The following model was attained after incorporating the similarity transformations and partial derivatives in Equations (1)–(4):
1 ( 1 34.87 ( d p d f ) 0.3 ϕ 1.03 ) [ ( η + K ) 3 L + 2 ( η + K ) 2   L ( η + K ) L + L ] [ 1 + 3 ( Θ 1 ) ϕ ( Θ + 2 ) ( Θ 1 ) ϕ ] M [ ( η + K ) 3 L + ( η + K ) 2 L ] + ( ( 1 ϕ ) + ϕ ( ρ β ) p ( ρ β ) f ) α ( ( η + K ) 3 N + ( η + K ) 2 L ) + { ( 1 ϕ ) + ϕ ρ p ρ f } K [ 0.5 ( 1 3 m ) ( η + K ) ( L ) 2 + 0.5 ( m + 1 ) ( η + K ) L L 0.5 ( m + 1 ) L L + 0.5 ( m + 1 ) ( η + K ) 2 L L + ( η + K ) 2 0.5 ( 1 3 m ) L L ] = 0   ,
R d [ 3 ( N ) 3 ( 1 + ( β w 1 ) N ) 2 ( β w 1 ) ( η + K ) + ( 1 + ( β w 1 ) N ) 3 ( ( η + K ) N + N ) ] + [ 1 + 4.4 R e b 0.4 Pr 0.66 ( T T f r e e z i n g ) 10 ( k p k f ) 0.03 ϕ 0.66 ] ( ( η + K ) N + N ) + P r E c ( 1 34.87 ( d p d f ) 0.3 ϕ 1.03 ) ( η + K ) ( ( η + K ) L L ) 2 = 0 .
In Equations (20) and (21), L and N are the functions of η . After solving, L and N provide the velocity and temperature distributions, respectively.
L ( η ) = 0 ,     L ( η ) = 1 ,       N ( η ) = 1 ,
At η :
L ( η ) 0 , L ( η ) 0 , N ( η ) 0 ,
Dimensionless quantities were described by the following formulas:
K = a   b ν b f , E c = b 2 s 2 m ( C p ) f ( T w T ) , R d = 16 σ ^ T 3 3 k ^ k f , M = σ f B 0 2 ρ f b , P r = ( c p ) f μ f k f , α = G r s ( R e s 2 ) 1 , R e s = b s 2 ( ν f ) 1 , G r s = s 3 ( T w T ) g β ( ν f ) 2 .
Moreover, the dimensional formula for skin friction and local Nusselt number were defined as:
C f = τ r s ( ρ f   u   ^ w 2   ) 1 , N u s = s q w ( k f ( T w T ) ) 1 ,
where,
τ r s = ( u   ^ r u   ^ ( a + r ) ) r = 0 , q w = k n f T   ^ r r = 0 ,
Drawing the values from Equation (26) into Equation (25), the following dimensionless formulas were obtained:
C f ( R e s ) 1 2 = 1 ( 1 34.87 ( d p d f ) 0.3 ϕ 1.03 ) ( L ( 0 ) L ( 0 ) K 1 ) ,
N u s ( R e s ) 1 2 = ( ( 1 + 4.4 R e b 0.4 Pr 0.66 ( T T f r e e z i n g ) 10 ( k p k f ) 0.03 ϕ 0.66 ) + R d β w 3 ) N ( 0 ) .

3. Mathematical Analysis

Highly nonlinear and coupled systems of differential equations usually possess no closed-form solution. Our flow model for Al2O3 + H2O is a highly nonlinear fourth-order model defined at a semi-infinite domain. Therefore, the model was tackled numerically using the Runge-Kutta scheme [13,16,25], as the RK scheme is used for the first order initial value problem (IVP). First the following substitutions were made and the model was reduced into first order IVP.
h 1 = L ,   h 2 = L ,   h 3 = L ,   h 4 = L ,   h 5 = N ,   h 6 = N .
After the successful transformation into the first order IVP, Mathematica 10.0 was used and the system was solved successfully.

4. Graphical Results and Discussion

This section emphasizes the flow and thermophysical characteristics of the fluid phase and the nanoparticles of Al2O3. The values for the thermophysical characteristics were calculated at 310 K [23]. The results for the shear stress and heat transfer rate are elaborated using bar charts and are discussed comprehensively.

4.1. Velocity and Temperature Distribution

Magnetic field phenomena are of a great significance from an industrial point of view. Many industrial productions contain impurities that need to be removed. However, the magnetic parameter opposes the fluid motion, and the impurities remain at the bottom and the nanofluid velocity drops. The influence of Lorentz forces on the velocity distribution of Al2O3 + H2O nanofluids is elaborated in Figure 3a. It was shown that the applied magnetic field opposed the nanofluid motion, and the velocity of the Al2O3 + H2O nanofluid dropped. The velocity declined more slowly for a weaker magnetic field, and a rapid decrement in the nanofluid velocity was observed for a stronger magnetic field. Near the arched surface, variations in the velocity ( L ( η ) ) were almost negligible. This behavior of the velocity distribution was due to the friction between the surface and the nanolayer of Al2O3 + H2O. In the successive nanolayers, the velocity field was altered significantly. These influences became negligible far from the curve and showed an asymptotic pattern of velocity distribution at the free surface. Figure 3b shows the velocity distribution of the parameter m. The velocity of the Al2O3 + H2O nanofluid dropped rapidly for m in comparison with M. As the values for parameter m became larger, the velocity decreased promptly.
The effects of surface curvature on the velocity ( L ( η ) ) are elucidated in Figure 4a. Altering the surface curvature caused the velocity to increase. For a smaller curvature, the velocity increased slowly and then vanished asymptotically far away from the surface. Similar behavior of the velocity of the Al2O3 + H2O nanofluid is depicted in Figure 4b. For α , a prominent behavior of the velocity was noticed in 0.5 η 1.5 . Besides this, the velocity ( L ( η ) ) was almost inconsequential.
The temperature distribution ( N ( η ) ) for the radiation parameter ( R d ) and Eckert number is highlighted in Figure 5a,b, respectively. In Figure 5a, we can see rapid drops in the temperature that were investigated by altering the radiation parameter ( R d ). The temperature ( N ( η ) ) dropped rapidly for a stronger radiation parameter, however at the free surface this was almost negligible and vanished asymptotically. The Eckert number, which appeared due to viscous dissipation, played a vibrant role in the heat transfer enhancement. These effects are elucidated in Figure 5b. It was obvious that the temperature of Al2O3 + H2O nanofluid grew for the more dissipative nanofluid. For the larger Ec, the temperature distribution rose rapidly.

4.2. Streamlines and Isotherms

This subsection is devoted to analyzing the behavior of the streamlines and isotherm patterns by altering different pertinent flow parameters. Figure 6 presents the streamline pattern by a varying magnetic parameter (M). For a smaller magnetic parameter, the streamlines were more curved near the surface and for a stronger M, the streamlines assumed a less curved shape. At the free stream, these streamlines became straight, since, as already mentioned for the parameter m > 1 , it showed that it was a nonlinear stretching curved surface. The streamline pattern versus m is elaborated in Figure 7. It was noticed that by decreasing the parameter m, the streamlines assumed a more curved pattern. The curvature parameter showed a fascinating pattern for the streamlines, in comparison with M and m. These alterations are illustrated in Figure 8. Figure 9 depicts the flow pattern for varying α . For a higher α , the streamlines shrank and became almost straight at the top. Figure 10 elaborates the isotherm pattern for the radiation parameter (Rd). When there was more radiative nanofluid, the isotherms increased, and vice versa. Further, a 3D scenario of the isotherms is depicted in Figure 11.

4.3. Thermophysical Characteristics

This subsection describes the impacts of the volume fraction factor ( ϕ ) on the effective characteristics of the nanofluids, and the behavior of the shear stress and local heat transfer rate by varying embedded flow parameters.
Figure 12 describes the effects of the volume fraction ( ϕ ,) and diameter of the nanoparticles on the effective dynamic viscosity of Al2O3 + H2O. The volume fraction of Al2O3 showed a vibrant role in enhancing the dynamic viscosity of the nanofluid. The observed high dynamic viscosity corresponded to a greater volume fraction. On the other hand, the nanoparticle diameter ( d p ) induced inverse variations in the dynamic viscosity. Increasing the diameter of the nanoparticles caused the dynamic viscosity to drop. This means that nanoparticles with a smaller diameter are important to enhance the dynamic viscosity of nanofluids.
Figure 13 highlights the effective density ( ρ n f ) and heat capacity ( ( ρ c p ) n f ) of the nanofluid versus ϕ . ϕ and ρ n f were in direct proportion to each other, and the effective heat capacity dropped when ϕ increased. Therefore, smaller values of ϕ enhanced the effective heat capacity. Due to the high volume fraction, the colloidal suspension Al2O3 + H2O became denser, which enhanced the effective density ( ρ n f ). Figure 14 highlights that the volume fraction ( ϕ ) and the effective electrical conductivity were in inverse proportion to each other.

4.4. Skin Fraction and Heat Transfer Rate

The shear stress and the local rate of heat transfer are very interesting and have attained great significance from an industrial point of view. The radiation parameter, Eckert number, and the curvature parameter all play a significant role in the shear stress and local Nusselt number. Figure 15 describes the heat transfer behavior for Rd, Ec, K and the volume fraction of the nanoparticles. It was noted that the more radiative nanofluids favored the heat transfer. On the other hand, a smaller amount of the heat transfer was noticed for a higher Eckert number. Therefore, the less dissipative fluids similarly favored the heat transfer. For a more curved surface, a large curvature worked against the heat transfer. At greater volume concentrations the heat transfer rate grew slowly.
Figure 16 elucidates the shear stress behavior for the mixed convection parameter ( α ), magnetic parameter (M), curvature of the surface, and m. For a more convective fluid, the heat transfer increased rapidly at the surface. The magnetic parameter (M) highlighted the reverse behavior of the shear stress. With an increase of the parameter M, the shear stress dropped quickly, and an almost negligible influence of m on the shear stresses was also observed. Significant alterations were pointed out for surface curvature, as the curve was along the circle loop of radius a. Therefore, for a smaller radius, shear stress declined promptly. Increasing the radius of the loop caused the surface curvature to become larger, which favored the shear stress.

5. Conclusions

The mixed convective laminar flow of water, composed by Al2O3 nanofluids in the presence of Lorentz forces and nonlinear radiative heat flux, was examined over an arc-shaped geometry. To enhance the heat transfer rate, a thermal conductivity model that considered the impact of freezing temperature and molecular diameter was used. It was found that the nanofluid velocity ( L ( η ) ) dropped for a stronger magnetic parameter (M), which is very significant from an industrial point of view. Further, the velocity of the nanofluid increased due to the mixed convection and larger curvature of the surface. The temperature N ( η ) intensified for more dissipative fluid, and an inverse relationship between the temperature and the thermal radiation parameter was found. The dynamic viscosity of the nanofluid increased with the volume fraction, and the diameter of the nanoparticles showed reverse alterations for dynamic viscosity. The nanofluid became denser for a high volume fraction, and the electrical conductivity dropped. It was found that the heat transfer reduced the surface of a smaller curvature and intensified for larger curvatures. A more convective fluid and a larger curvature better opposed the shear stress. Finally, when considering the influence of the freezing temperature and the molecular diameter, the nanofluid flow model is very useful for heat transfer in comparison with existing studies.

Author Contributions

Conceptualization, U.K. and A.A.; methodology, Adnan, N.A.; software, S.N.; validation, writing—review and editing, I.K. and W.A.K.; supervision, S.T.M.-D; Funding acquisition, S.O.A. All the authors contributed to the manuscript equally.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

u Component of the velocity
a Radius
T Temperature
T Temperature far from the surface
k f Thermal conductivity of the host fluid
k s Thermal conductivity of the nanoparticles
k n f Effective thermal conductivity of the nanofluid
( C p ) f Heat capacity of the host fluid
σ s Electrical conductivity of the nanoparticles
σ n f Electrical conductivity of the nanofluid
μ ) _ f   Dynamic viscosity of the fluid
M * Molecular weight
d f Molecular diameter
β w Temperature ratio parameter
P r Prandtl number
M Hartmann number
L ( η ) Dimensionless velocity
N u Nusselt number
v Component of the velocity
p Pressure
T w Temperature at the surface
ρ f Density of the host fluid
ρ s Density of the nanoparticles
ρ n f Effective density of the nanofluid
( C p ) _ s Heat capacity of the nanoparticles
( C p ) n f Heat capacity of the nanofluid
σ f Electrical conductivity of the host fluid
μ n f Effective dynamic viscosity
ϕ Volume fraction of the nanoparticles
N * Avogadro number
k b Stefan Boltzmann constant
R d Radiation parameter
E c Eckert number
K Curvature parameter
N ( η ) Dimensionless temperature
C f Skin fraction coefficient

References

  1. Choi, S. Enhancing thermal conductivity of fluids with nanoparticles in developments and applications of non-newtonians flows. ASME J. Heat Transf. 1995, 66, 99–105. [Google Scholar]
  2. Clerk, M.J. Treatise on Electricity and Magnetism; Oxford University Press: Oxford, UK, 1873. [Google Scholar]
  3. Bruggeman, D.A.G. Berechnung verschiedener physikalischer konstanten von heterogenen substanzen, I—dielektrizitatskonstanten und leitfahigkeiten der mischkorper aus isotropen substanzen. Ann. Phys. Leipz. 1935, 24, 636–679. [Google Scholar] [CrossRef]
  4. Hamilton, H.L.; Crosser, O.K. Thermal conductivity of heterogeneous two-component systems. Ind. Eng. Chem. Fundam. 1962, 1, 187–191. [Google Scholar] [CrossRef]
  5. Lu, S.; Lin, H. Effective conductivity of composites containing aligned spherical inclusions of finite conductivity. J. Appl. Phys. 1996, 79, 6761–6769. [Google Scholar] [CrossRef]
  6. Koo, J.; Kleinstreuer, C. A new thermal conductivity model for nanofluids. J. Nanopart. Res. 2004, 6, 577–588. [Google Scholar] [CrossRef]
  7. Koo, J.; Kleinstreuer, C. Laminar nanofluid flow in micro-heat sinks. Int. J. Heat Mass Transf. 2005, 48, 2652–2661. [Google Scholar] [CrossRef]
  8. Xue, Q.Z. Model for thermal conductivity of carbon nanotube-based composites. Phys. B Phys. Condens. Matter. 2005, 368, 302–307. [Google Scholar] [CrossRef]
  9. Prasher, R.; Bhattacharya, P.; Phelan, P.E. Thermal conductivity of nanoscale colloidal solutions (nanofluids). Phys. Rev. Lett. 2005, 92, 25901. [Google Scholar] [CrossRef]
  10. Li, C.H.; Peterson, G.P. Experimental investigation of temperature and volume fraction variations on the effective thermal conductivity of nanoparticle suspensions (nanofluids). J. Appl. Phys. 2006, 99. [Google Scholar] [CrossRef]
  11. Corcione, M. Rayleigh–Be´nard convection heat transfer in nanoparticle suspensions. Int. J. Heat Fluid Flow 2011, 32, 65–77. [Google Scholar] [CrossRef]
  12. Sheikholeslami, M.; Li, Z.; Shamlooei, M. Nanofluid MHD natural convection through a porous complex shaped cavity considering thermal radiation. Phys. Lett. A 2018, 382, 1615–1632. [Google Scholar] [CrossRef]
  13. Ahmed, N.; Khan, A.U.; Mohyud-Din, S.T. Influence of an effective prandtl number model on squeezed flow of γAl_2 O_3-H_2 O and γAl_2 O_3-C_2 H_6 O_2 nanofluids. J. Mol. Liq. 2017, 238, 447–454. [Google Scholar] [CrossRef]
  14. Sheikholeslami, M.; Zia, Q.M.Z.; Ellahi, R. Influence of induced magnetic field on free convection of nanofluid considering Koo-Kleinstreuer-Li (KKL) correlation. Appl. Sci. 2016, 6, 324. [Google Scholar] [CrossRef]
  15. Asadullah, A.M.; Khan, U.; Naveed, A.; Mohyud-Din, S.T. Analytical and numerical investigation of thermal radiation effects on flow of viscous incompressible fluid with stretchable convergent/divergent channels. J. Mol. Liq. 2016, 224, 768–775. [Google Scholar]
  16. Khan, U.; Naveed, A.A.; Mohyud-Din, S.T. 3D squeezed flow of γAl_2 O_3-H_2 O and γAl_2 O_3-C_2 H_6 O_2 nanofluids: A numerical study. Int. J. Hydrog. Energy 2017, 42, 24620–24633. [Google Scholar] [CrossRef]
  17. Iijima, S. Helical microtubules of graphitic carbon. Nature 1991, 354, 56–58. [Google Scholar] [CrossRef]
  18. Naveed, A.A.; Khan, U.; Mohyud-Din, S.T. Influence of thermal radiation and viscous dissipation on squeezed flow of water between two riga plates saturated with carbon nanotubes. Colloids Surf. A Physciochem. Eng. Asp. 2017, 522, 389–398. [Google Scholar]
  19. Saba, F.; Naveed, A.; Hussain, S.; Khan, U.; Mohyud-Din, S.T.; Darus, M. Thermal analysis of nanofluid flow over a curved stretching surface suspended by carbon nanotubes with internal heat generation. Appl. Sci. 2018, 8, 395. [Google Scholar] [CrossRef]
  20. Khan, U.; Naveed, A.; Mohyud-Din, S.T. Heat transfer effects on carbon nanotubes suspended nanofluid flow in a channel with non-parallel walls under the effect of velocity slip boundary condition: A numerical study. Neural Comput. Appl. 2017, 28, 37–46. [Google Scholar] [CrossRef]
  21. Reddy, J.V.R.; Sugunamma, V.; Sandeep, N. Dual solutions for nanofluid flow past a curved surface with nonlinear radiation, soret and dufour effects. J. Phys. Conf. Ser. 2018, 1000, 12152. [Google Scholar] [CrossRef]
  22. Hayat, T.; Qayyum, S.; Imtiaz, M.; Alsaedi, A. Double stratification in flow by curved stretching sheet with thermal radiation and joule heating. J. Therm. Sci. Eng. Appl. 2017, 10. [Google Scholar] [CrossRef]
  23. Alsabery, A.I.; Sheremet, M.A.; Chamkha, A.J.; Hashim, I. MHD convective heat transfer in a discretely heated square cavity with conductive inner block using two-phase nanofuid model. Sci. Rep. 2018, 8, 1–23. [Google Scholar] [CrossRef]
  24. Corcione, M. Empirical correlating equations for predicting the efective thermal conductivity and dynamic viscosity of nanofuids. Energy Convers. Manag. 2011, 52, 789–793. [Google Scholar] [CrossRef]
  25. Naveed, A.A.; Khan, U.; Mohyud-Din, S.T. Unsteady radiative flow of chemically reacting fluid over a convectively heated stretchable surface with cross-diffusion gradients. Int. J. Therm. Sci. 2017, 121, 182–191. [Google Scholar]
Figure 1. Flow description in a curvilinear frame.
Figure 1. Flow description in a curvilinear frame.
Symmetry 11 00833 g001
Figure 2. Flowchart of the study.
Figure 2. Flowchart of the study.
Symmetry 11 00833 g002
Figure 3. Impacts of (a) M and (b) m on the velocity distribution ( L ( η ) ).
Figure 3. Impacts of (a) M and (b) m on the velocity distribution ( L ( η ) ).
Symmetry 11 00833 g003
Figure 4. Impacts of (a) K and (b) α on the velocity distribution ( L ( η ) ).
Figure 4. Impacts of (a) K and (b) α on the velocity distribution ( L ( η ) ).
Symmetry 11 00833 g004
Figure 5. Impacts of (a) Rd and (b) Ec on the temperature distribution ( N ( η ) ).
Figure 5. Impacts of (a) Rd and (b) Ec on the temperature distribution ( N ( η ) ).
Symmetry 11 00833 g005
Figure 6. Streamlines for (a) M = 6 (b) M = 2.
Figure 6. Streamlines for (a) M = 6 (b) M = 2.
Symmetry 11 00833 g006
Figure 7. Streamlines for (a) m = 4 and (b) m = 2.
Figure 7. Streamlines for (a) m = 4 and (b) m = 2.
Symmetry 11 00833 g007
Figure 8. Streamlines for (a) k = 8 and (b) k = 5.
Figure 8. Streamlines for (a) k = 8 and (b) k = 5.
Symmetry 11 00833 g008
Figure 9. Streamlines for (a) α = 6 and (b) α = 2 .
Figure 9. Streamlines for (a) α = 6 and (b) α = 2 .
Symmetry 11 00833 g009
Figure 10. Isotherms for (a) Rd = 6 and (b) Rd = 2.
Figure 10. Isotherms for (a) Rd = 6 and (b) Rd = 2.
Symmetry 11 00833 g010
Figure 11. 3D scenario of Figure 10; (a) Rd = 6 and (b) Rd = 2.
Figure 11. 3D scenario of Figure 10; (a) Rd = 6 and (b) Rd = 2.
Symmetry 11 00833 g011
Figure 12. The impact of the volume fraction ( ϕ ) and the nanoparticle diameter ( d p ) on the dynamic viscosity.
Figure 12. The impact of the volume fraction ( ϕ ) and the nanoparticle diameter ( d p ) on the dynamic viscosity.
Symmetry 11 00833 g012
Figure 13. The impact of ϕ on the effective density and heat capacity.
Figure 13. The impact of ϕ on the effective density and heat capacity.
Symmetry 11 00833 g013
Figure 14. The impact of ϕ on the effective electrical conductivity.
Figure 14. The impact of ϕ on the effective electrical conductivity.
Symmetry 11 00833 g014
Figure 15. The impact of various flow parameters on the local heat transfer.
Figure 15. The impact of various flow parameters on the local heat transfer.
Symmetry 11 00833 g015
Figure 16. The impact of various parameters on the shear stress.
Figure 16. The impact of various parameters on the shear stress.
Symmetry 11 00833 g016
Table 1. Thermal and physical characteristics of the fluid phase and nanoparticles at T = 310 K [23].
Table 1. Thermal and physical characteristics of the fluid phase and nanoparticles at T = 310 K [23].
Properties d p   ( nm ) ρ (kg/m3) β (1/k) c p (J/Kg K) μ f (kg/ms) k (W/mk) σ (S/m)
H2O 0.385 993 36.2 × 10 5 4178 695 × 10 6 0.628 0.005
Al2O333 3970 0.85 × 10 5 765 - 40 0.05 × 10 6

Share and Cite

MDPI and ACS Style

Khan, U.; Abbasi, A.; Ahmed, N.; Alharbi, S.O.; Noor, S.; Khan, I.; Mohyud-Din, S.T.; Khan, W.A. Modified MHD Radiative Mixed Convective Nanofluid Flow Model with Consideration of the Impact of Freezing Temperature and Molecular Diameter. Symmetry 2019, 11, 833. https://doi.org/10.3390/sym11060833

AMA Style

Khan U, Abbasi A, Ahmed N, Alharbi SO, Noor S, Khan I, Mohyud-Din ST, Khan WA. Modified MHD Radiative Mixed Convective Nanofluid Flow Model with Consideration of the Impact of Freezing Temperature and Molecular Diameter. Symmetry. 2019; 11(6):833. https://doi.org/10.3390/sym11060833

Chicago/Turabian Style

Khan, Umar, Adnan Abbasi, Naveed Ahmed, Sayer Obaid Alharbi, Saima Noor, Ilyas Khan, Syed Tauseef Mohyud-Din, and Waqar A. Khan. 2019. "Modified MHD Radiative Mixed Convective Nanofluid Flow Model with Consideration of the Impact of Freezing Temperature and Molecular Diameter" Symmetry 11, no. 6: 833. https://doi.org/10.3390/sym11060833

APA Style

Khan, U., Abbasi, A., Ahmed, N., Alharbi, S. O., Noor, S., Khan, I., Mohyud-Din, S. T., & Khan, W. A. (2019). Modified MHD Radiative Mixed Convective Nanofluid Flow Model with Consideration of the Impact of Freezing Temperature and Molecular Diameter. Symmetry, 11(6), 833. https://doi.org/10.3390/sym11060833

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