Next Article in Journal
Development of CO2 Concentration Prediction Tool for Improving Office Indoor Air Quality Considering Economic Cost
Next Article in Special Issue
Interchangeability of Hydrogen Injection in Zhejiang Natural Gas Pipelines as a Means to Achieve Carbon Neutrality
Previous Article in Journal
High Mining Face Flexible Reinforcement to Prevent Coal Wall Spalling by Cuttable Aluminum–Plastic Pipe Pre-Grouting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Friction Factor and Heat Transfer of Giesekus-Fluid-Based Nanofluids in a Pipe Flow

1
School of Media and Design, Hangzhou Dianzi University, Hangzhou 310018, China
2
State Key Laboratory of Fluid Power Transmission and Control, Zhejiang University, Hangzhou 310027, China
3
Laboratory of Impact and Safety Engineering (Ningbo University), Ministry of Education, Ningbo 315201, China
*
Author to whom correspondence should be addressed.
Energies 2022, 15(9), 3234; https://doi.org/10.3390/en15093234
Submission received: 26 March 2022 / Revised: 20 April 2022 / Accepted: 26 April 2022 / Published: 28 April 2022
(This article belongs to the Special Issue Multiphase Flows Related to Energies)

Abstract

:
The friction factor and heat transfer of Giesekus-fluid-based nanofluids in a pipe flow were studied in the ranges of 0.5 ≤ Reynolds number (Re) ≤ 500, 1 ≤ Weissenberg number (Wi) ≤ 8, 0.5% ≤ particle volume concentration (Φ) ≤ 3.0%, 0 ≤ viscosity ratio (β0) ≤ 1, and 0 ≤ mobility parameter (α) ≤ 0.5. Our numerical method was validated by comparing the results with available ones in the literature. The effects of Wi, Φ, β0, Re, and α on the relative friction factor (Cf/CfNew), Nusselt number (Nu), and ratio (PECnf/PECf) of energy performance evaluation criterion for Giesekus-fluid-based nanofluids to those for Giesekus fluid were discussed. The results showed that the values for the Cf/CfNew and Nu of Giesekus-fluid-based nanofluids were larger than those for Newtonian fluid-based nanofluids and those for pure Giesekus fluid. The values for Cf/CfNew increased with increasing Φ and Re, but they increased with decreasing β0 and α. As Wi increased, the values of Cf/CfNew first increased and then decreased. The values of Nu and PECnf/PECf were enhanced with increasing Wi, Φ, Re, and α, but with decreasing β0. It is more effective to use Giesekus-fluid-based nanofluids to improve heat transfer with the conditions of a larger Wi, Φ, Re, and α and a smaller β0. Finally, the correlation formula for PECnf/PECf as a function of Wi, Φ, β0, Re, and α was derived.

1. Introduction

The friction factor and heat transfer of non-Newtonian fluids have aroused extensive concern due to their wide range of applications, e.g., drilling processes, synthesis of composites, paper making, and plastic molding [1,2]. In order to reduce friction and increase the heat transfer of fluid flowing in a pipe and achieve the goal of saving energy, nanofluids, i.e., fluids containing nanoparticles, are used. So far, nanofluids have been extensively used, and their properties have also been widely investigated. For example, adding graphite and carbon nanofibers increases the efficiency of a gear pump-driven hydraulic circuit using ethanol [3]. The thermal efficiency was increased by 13.8%, 1.5%, and 1.3% for nanofluids of water–CuO compared with pure water, water–TiO2, and water–Al2O3, respectively [4]. Low friction factor and high heat transfer were reached when the nanofluids of water-Al2O3 and water-SiO2 were used in a heat exchanger in the range of 10,551 ≤ Re ≤ 17,220 and 17,220 ≤ Re ≤ 31,910, respectively [5].
Newtonian fluid was employed as the base fluid in most of the previous studies. However, there are many occasions where the base fluid is non-Newtonian fluid in applications; only a few studies have focused on non-Newtonian fluid as the base fluid. For instance, the Deborah number had an important impact on convective heat transfer, and the effect of heat conduction of Jeffrey nanofluids was better than that of Oldroyd-B nanofluids [6]. The impacts of fluid relaxation and retardation time parameters on the velocity of Oldroyd-B nanofluids were converse [7]. The heat transfer performance of Jeffrey nanofluids was better than that of the Oldroyd-B nanofluids when taking the thermal radiation, suction, and nonuniform heat source/sink into account [8]. The heat transfer performance was greatly enhanced in the Oldroyd-B nanofluids at higher Biot and Prandtl numbers [9]. The heat transfer rate of EG-Si3N4 was higher than that of EG for various nanoparticle concentrations and flow rates [10]. The heat transfer rate of Oldroyd-B nanofluids was much higher than those of Jeffery and Maxwell nanofluids, while the friction factor of Maxwell nanofluids was less than those of Oldroyd-B and Jeffrey nanofluids over a permeable stretching sheet [11]. Polyalphaolefins-based nanofluids with carbon tubes promoted heat conductivity, followed by exfoliated graphite and heat-treated nanofibers [12]. The precipitation of particles on a wall enhanced the heat transfer rate in polyalphaolefins-based nanofluids with exfoliated graphite fibers [13]. The thermophysical properties of polyalphaolefins-based nanofluids containing Al2O3 nanoparticles were dependent on particle volume fraction, particle aggregation, and dispersion [14].
The non-Newtonian fluids mentioned above involved Oldroyd-B, Jeffrey, and Maxwell fluids. The constitutive equations of these three non-Newtonian fluids can account for elasticity; hence, they are suited to model viscoelastic fluid but cannot account for the shear-thinning property of fluid. The constitutive equations of Jeffrey and Maxwell fluids are linear, and they can only be used to describe some simple flows. In practical applications, non-Newtonian fluid with both viscoelastic and shear-thinning characteristics is also common. To describe both characteristics, the constitutive equation of Giesekus fluid can be used. Hence, Giesekus fluid has attracted wide attention. However, a literature survey revealed that there is little research available on the friction factor and heat transfer of Giesekus-fluid-based nanofluids. Yang et al. [15] investigated the turbulent flow and heat transfer of Giesekus-fluid-based nanofluids in a channel and drew the conclusion that the viscoelasticity of Giesekus fluid resulted in drag reduction and the enhancement of convective heat transfer. Wang et al. [16] studied drag reduction in the gas–liquid two-phase turbulent flow of Giesekus fluids and indicated that it is better to use polymer drag-reducing agents with high concentrations in high shear flow. With the method of using Giesekus-fluid-based nanofluids, the effects of reducing friction and enhancing heat transfer performance are dependent on the Reynolds number, Weissenberg number, particle volume concentration, viscosity ratio, and mobility parameter, and there is still a lack of research on this issue. In addition, in previous studies, the distribution of particle volume concentration is usually assumed to be uniform without taking particle convection, diffusion, coagulation, or breakage into account. In this work, therefore, a numerical simulation was conducted to investigate the effects of the Reynolds number, the Weissenberg number, rheological properties, and particle concentration on the friction factor and heat transfer performance of Giesekus-fluid-based Al2O3 nanofluids in a pipe flow while considering the effects of particle convection, diffusion, coagulation, and breakage.
The rest of this paper is organized as follows. The flow and equations of nanofluids and particle equations are introduced in Section 2 and Section 3, respectively. The numerical method is given in Section 4. In Section 5, the numerical results are presented and discussed. Finally, conclusions are drawn in Section 6.

2. Flow and Equations of Nanofluids

The flow of Giesekus-fluid-based nanofluids in a pipe is shown in Figure 1. For the incompressible flow, the continuity and momentum equations of Giesekus-fluid-based nanofluids are [17]:
u = 0
ρ nf ( u t + u u ) = σ
σ = p I + 2 μ nf D + τ
where is the differential operator; u, p, and ρnf are the velocity, pressure, and density of nanofluids, respectively; σ is the total stress tensor; I is the unit tensor; μnf is the solvent viscosity; τ is the stress tensor; and D is the rate-of-deformation tensor:
D = [ u + ( u ) T ] 2
The constitutive equation of Giesekus fluid is:
λ τ + α λ μ p τ τ + τ = 2 μ p D
where λ is the fluid relaxation time; μp is the contribution of the polymer to viscosity; α is the mobility parameter; and the symbol ∇ is the upper-convected time derivative:
τ τ t + u τ ( u ) T τ τ u
The density and viscosity of nanofluids in Equations (2) and (3) are [18]:
ρ nf = ( 1 Φ ) ρ f + Φ ρ p ,           μ nf = μ f ( 1 + 2.5 Φ + 6.2 Φ 2 )
where Φ is the particle concentration and subscripts f and p denote fluid and particles, respectively.
The energy equation is:
T t + u T = D nf 2 T
where T is the nanofluid temperature and Dnf is the thermal diffusivity coefficient of nanofluids:
D nf = k nf ( ρ C p ) nf
where heat capacitance (ρCp)nf and knf are [19]:
k nf = k f ( k p + 2 k f ) 2 Φ ( k f k p ) ( k p + 2 k f ) + Φ ( k f k p ) ,       ( ρ C p ) nf = ( 1 Φ ) ( ρ C p ) f + Φ ( ρ C p ) p
where k is thermal conductivity.

3. Equations of Particles

Particle concentration Φ in Equations (7) and (10) should be determined before solving Equations (1)–(6). The diffusion, convection, coagulation, and breakage of nanoparticles leads to a nonuniform distribution of Φ in the flow process; hence, it is necessary to solve the equation of particle concentration.

3.1. Dynamics Equation for Nanoparticles

Considering the convection, diffusion, coagulation, and breakage of nanoparticles, the equation of particle volume can be given as follows under a steady condition:
u j n ( v ) x j x j [ D p n ( v ) x j ] = 1 2 0 v β ( v 1 , v v 1 ) n ( v 1 ) n ( v v 1 ) d v 1 0 β ( v 1 , v ) n ( v ) n ( v 1 ) d v 1 + 1 2 0 v β ( v 1 , v v 1 ) n ( v 1 ) n ( v v 1 ) ¯ d v 1 0 β ( v 1 , v ) n ( v ) n ( v 1 ) ¯ d v 1 + v a ( v 1 ) b ( v | v 1 ) n ( v 1 ) d v 1 a ( v ) n ( v )
where uj is the nanofluids’ velocity; the subscript j indicates summation; n(v) is the distribution function of particle volume; Dp = kBT/3πμnfdp is the diffusion coefficient of particles with kB being the Boltzmann constant, μnf the viscosity, and dp the particle diameter; β(v, v1) is the coagulation kernel for two particles of volume v and v1; α(v) is the breakage kernel expressing the breakage frequency of particle v; and b(v|v1) is the breakage distribution function.

3.2. Moment Equation of Nanoparticles

Multiplying Equation (11) with vk and then integrating over v yields:
u j M k x j x j ( D p M k x j ) = 1 2 0 0 [ ( v + v 1 ) k v k v 1 k ] β ( v , v 1 ) n ( v ) n ( v 1 ) d v d v 1 + 0 v k 0 a ( v 1 ) b ( v | v 1 ) n ( v 1 ) d v 1 d v 0 v k a ( v ) n ( v ) d v
where moment Mk is:
M k = 0 v k n ( v ) d v
Taking k = 0 and 1 yields:
M 0 = 0 n ( v ) d v = N ,           M 1 = 0 v n ( v ) d v = V
where M0 and M1 are the total particle number and particle volume, respectively. Based on Equation (14), Φ can be calculated.

3.3. Particle Coagulation and Breakage

The coagulation kernel in Equation (12) is [20]:
β = B ( 1 v 1 / 3 + 1 v 1 1 / 3 ) ( v 1 / 3 + v 1 1 / 3 )
where B = 2kBT/μnf with kB being the Boltzmann constant, T the fluid temperature, and μnf the nanofluids’ viscosity.
The particle breakage function is composed of breakage kernel α(v) and breakage distribution function b(v|v1). The breakage kernel α(v) is [21]:
a ( v ) = k B ( ν nf G τ ) q v p 1 / 3 ( d c ( v ) d p ) 3 / D f
where νnf is the nanofluids’ viscosity; G is the shear rate of flow; τ* is characteristic shear stress and is a strength measure of the coagulated particles; q, vp, and dp are constants dependent on Φ; dc is particle collision diameter and dependent on the number of primary particles; and Df is the fractal dimension of particles (Df = 3 for spherical particles).
The breakage distribution function, b(v|v1), is [22]:
b ( v v 1 ) = 2         if   v 1 = 2 v 0         otherwise
Replacing the second term on the right-hand side of Equation (12) with Equation (17) yields:
0 v k 0 a ( v 1 ) b ( v | v 1 ) n ( v 1 ) d v 1 d v = 0 [ 0 v k b ( v v 1 ) d v ] a ( v 1 ) n ( v 1 ) d v 1 = 0 2 1 k v i k a ( v 1 ) n ( v 1 ) d v 1

3.4. Taylor Series Expansion Moment Method

The Taylor series expansion technique [23,24] is used to solve moment equations. The equations of zero-, first-, and second-order moment can be given by substituting Equations (15)–(18) into Equation (12):
u j M 0 x j x j ( D p M 0 x j ) = B 2 ( 151 M 1 4 + 2 M 2 2 M 0 2 13 M 2 M 1 2 M 0 ) M 0 2 81 M 1 4
u j M 1 x j x j ( D p M 1 x j ) = 0
u j M 2 x j x j ( D p M 2 x j ) = 2 B 2 ( 2 M 2 2 M 0 2 13 M 2 M 1 2 M 0 151 M 1 4 ) 81 M 1 2

4. Numerical Method

4.1. Steps of Numerical Simulation

The numerical simulation included the following steps: (1) solving Equations (1)–(7) with Φ = 0 to get u, p, and τ; (2) solving Equations (11)–(21) to get n and Φ; (3) substituting Φ into Equations (7), (9), and (10) to get ρnf, μnf, Dnf, kt, and (ρCp)nf; (4) substituting Φ, ρnf, μnf, Dnf, kt, and (ρCp)nf into Equations (1)–(8) to get u, p, τ, and T; (5) turning to Step (2) based on the new values of u, p, and τ if necessary; (6) calculating the Fanning friction factor Cf and Nusselt number Nu:
C f = Δ p 4 ρ nf ( L / D ) ( u z , av 2 / 2 ) ,       Nu = h D k nf
where uz,av and ∆p are the mean velocity and pressure drop of the fluid, respectively, and h is the heat transfer coefficient; L and D are shown in Figure 1.

4.2. Numerical Method and Boundary Condition

The finite-volume method was used to solve Equations (2), (5), (8), and (19)–(21). The SIMPLE scheme was employed to handle the term of velocity–pressure coupling, and the power-law scheme was used to handle the convection term. The discretized equations were solved using a staggered mesh system and an alternating direction implicit method. The mesh size was uniform along z and θ directions, while it was refined in r direction. On the wall, the velocity, zero-, first-, and second-order moment were zero; the heat flux was constant. The initial distributions of zero-, first-, and second-order moment and temperature were uniform.

4.3. Parameters

Giesekus fluid has both viscoelasticity and shear-thinning characteristics, and Al2O3 nanoparticles have the features of high energy density and thermal conductivity, safety, and low cost. Therefore, the Giesekus fluid and Al2O3 nanoparticles were selected as the base fluid and disperse phase, respectively, in the present study. The parameters used in the numerical simulation were: ρf = 990 kg/m3, ρp= 3970 kg/m3, kf = 0.65 W/m·K, kp = 40 W/m·K, Cpf = 4.15 kJ/kg·K, Cpp = 0.765 kJ/kg·K, and kB = 1.38 × 10−23 J/K. The Reynolds number was defined as Re = ρnfUD/μnf.

4.4. Validation

The flow field was divided into 128(r) × 32(θ) × 256(z) = 1,048,576 grids, and a test of independence for the grid size was performed, as shown in Table 1. The convergence criterion was that all the residual errors were less than 10−4.
To validate the numerical model used in the present simulation, the present numerical result of the axial velocity profile along the radial direction in the pipe flow of Giesekus fluid was compared with the result given by Vachagina et al. [25], as shown in Figure 2, where uz is the ratio of local velocity to average velocity and the corresponding parameters are: Wi = 4.29 and λ = 84.25 s−1. The present result of the relationship between the Nusselt number and Weissenberg number in the pipe flow of Giesekus fluid was compared with theoretical results [26], as shown in Figure 3, where Br = ηua2/Dqw with η being zero-shear viscosity, ua the average velocity, D the pipe diameter, and qw the heat flux. As can be seen from Figure 2 and Figure 3, our results were in agreement with previous results.

5. Results and Discussion

5.1. Friction Factor

The Fanning friction factor Cf, as shown in Equation (22), was used to express the shear stress exacting on the flow by the wall.

5.1.1. Effect of Weissenberg Number

The relationship between the relative friction factor, Cf/CfNew, and the modified Weissenberg number, Wi*, for different particle concentrations, Φ, is shown in Figure 4, where the results [27] for pure Giesekus fluid are also given for comparison. The parameters in Figure 4 are: CfNew is the friction factor for Newtonian fluid; Wi* = 2Wiτ0/(Reb0L), where Wiτ0 = λuτ2/ν0 (λ is relaxation time, uτ is friction velocity, and ν0 is zero shear-rate kinematic viscosity) is zero shear-rate friction Wi; Reb0 = Dub/ν0 (D is pipe diameter and ub is bulk average velocity), L is the maximum extensibility parameter of polymer molecules; β0 = μs/μ0 = μs/(μs + μp), where μs is solvent viscosity, μ0 is zero shear-rate viscosity, and μp is viscosity from polymer contributions; and α is a mobility parameter.
From Figure 4, we can see that Cf/CfNew is always larger than 1, i.e., the Cf of Giesekus-fluid-based nanofluids is larger than that of Newtonian fluid. Cf/CfNew = 1 when Wi* = 0 because Wi* = 2Wiτ0/(Reb0L) is the ratio of elastic force to inertial force; Wi* = 0 for Newtonian fluid. As Wi* increases, Cf/CfNew increases rapidly when 0 < Wi* < 0.5 and slowly when 0.5 < Wi* < 1, and then it decreases slowly when 1 < Wi*< 7. Wi is also the ratio of the first normal stress difference (FNSD) to shear stress. When 0 < Wi* < 1, the increase in Wi* means an increase in FNSD, so as to produce additional extension resistance and result in an increase in Cf. When 1 < Wi* < 7, the effect of the increase in FNSD is weakened; meanwhile, the effect of shear-thinning is increased, resulting in the decrease in Cf. Friction is the result of the combined action of various forces exacting on fluid. Therefore, there exists a maximum Cf for certain rheological and physical properties of nanofluids and certain flow conditions.

5.1.2. Effects of Particle Concentration and Viscosity Ratio

In Figure 4, Cf/CfNew increases with increasing Φ. The nanoparticles in Giesekus fluid are responsible for energy absorption, leading to a large Cf. In addition, nanofluids’ viscosity, μnf, is proportional to Φ, as shown in Equation (7). μnf increases with increasing Φ, intensifying the viscous sublayer and resulting in an increase in Cf. As Φ increases, the increase in Cf/CfNew is large when 0 < Φ < 0.5% but small when 0.5% < Φ < 3.0% because the variation in Φ has little impact on the viscous sublayer in this range of Φ, showing that the effect of particle concentration on Cf is nonlinear. In previous studies, there were different conclusions on the impact of nanoparticles in Newtonian-fluid-based nanofluids on Cf; most conclusions are that nanoparticles lead to an increase in Cf (e.g., [28,29,30,31]). The results in Figure 4 show that nanoparticles in Giesekus-fluid-based nanofluids also lead to an increase in Cf.
Figure 5 shows the relationship between Cf/CfNew and Wi* for different viscosity ratios, β0. β0 = 0 and 1 corresponds to pure polymer and Newtonian fluid because β0 = μs/(μs + μp). We can see that Cf/CfNew = 1 for Newtonian fluid (β0 = 1). As Wi* increases, Cf/CfNew increases rapidly when 0 < Wi* < 2 and slowly when 2 < Wi* < 7 for the pure polymer (β0 = 0), showing that Cf is directly proportional to Wi* in the absence of solvent viscosity μs. However, Cf/CfNew first increases and then gradually decreases with the increase in Wi* when 0 < β0 < 1. Therefore, the changing trend of Cf with Wi* depends on β0. In addition, Cf/CfNew increases with decreasing β0, i.e., a large β0 should be selected in order to reduce Cf.

5.1.3. Effects of Reynolds Number and Mobility Parameter

The Reynolds number represents the ratio of fluid inertia to viscosity, and its value represents the relative importance of the two. The range in Reynolds numbers in this paper was 0.5–500 because this can more fully reflect the situations of inertia or viscosity dominance. The relationship between Cf/CfNew and Wi* for different Re is shown in Figure 6, where Cf/CfNew is reduced with increasing Re, but the relationship between the reducing rate of Cf/CfNew and the increasing rate of Re is nonlinear, i.e., the change rate for Cf/CfNew is almost the same as the small change rate for Re when 0.5 < Re < 1, and there is a large change rate when 100 < Re < 500. Cf is proportional to viscosity and shear rate, and Re is the ratio of inertia to viscosity. At a small Re, the function of viscosity is large, and the shear rate is large in most regions of flow, so there is a large Cf. At a small Re, the flow is sensitive to changes in Re, so Cf is also sensitive to changes in Re. At a large Re, the inertial force is dominant against the viscous forces with increasing Re, so Cf is less sensitive to changes in Re.
Figure 7 shows the relationship between Cf/CfNew and Wi* for different mobility parameters, α; the degree of shear-thinning is proportional to α. Increasing α leads to a decrease in Cf/CfNew due to an increase in the shear-thinning effect. Cf/CfNew is smaller for Giesekus-fluid-based nanofluids than for Oldroyd-B-based nanofluids (α = 0) under the same conditions. For a definite α and larger Wi*, the change in Cf/CfNew becomes unobvious and Cf/CfNew approaches asymptotic values. Despite the shear-thinning effect, Cf/CfNew is still larger than 1 because there exists a polymer contribution, μp, to viscosity.

5.2. Heat Transfer

In order to understand the heat transfer performance of Giesekus-fluid-based nanofluids, we explored the impacts of Wi*, Φ, β0, Re, and α on Nu, which is shown in Equation (22).

5.2.1. Effect of Weissenberg Number

The relationship between Nu and Wi for different Φ is shown in Figure 8, where the results [32] for pure Giesekus fluid are also given. Nu is at its minimum when Wi = 0 for Φ = 0, i.e., the heat transfer effect of Giesekus fluid is better than that of Newtonian fluid. As Wi increases, Nu increases rapidly when 0 < Wi < 2 and slowly when 2 < Wi < 8. For a definite Φ and larger Wi, Nu approaches asymptotic values, indicating that the fluid elasticity has little impact on heat transfer when Wi exceeds a certain value.

5.2.2. Effects of Particle Concentration and Viscosity Ratio

In Figure 8, the Nu for Φ≠0 is larger than that for Φ = 0, indicating that adding nanoparticles into fluid can improve heat transfer. An increase in Φ usually leads to an increase in viscosity, as shown in Equation (7), and intensifies the viscous sublayer, which causes a decrease in heat transfer. However, the thermal conductivity of nanofluids is higher than that of pure fluid when adding nanoparticles with high thermal conductivity into fluid, which enhances the heat transfer effect. Therefore, the results in Figure 8 imply that the impact of thermal conductivity enhancement overcomes the impact of viscosity increase, which is consistent with previous results, e.g., [33]. Nu increases with increasing Φ because the interaction between particles is weak at a small Φ, meaning a weak heat transfer. At a large Φ, both frequent particle interaction and movement disturb the flow and enhance heat transfer. The increasing rate of Nu decreases with increasing Φ because the increase in viscosity is far greater than the increase in thermal conductivity at a large Φ. In addition, the nanoparticles with high Φ are easier to coagulate, resulting in a reduction in particle numbers, thereby affecting heat transfer performance.
Figure 9 shows the relationship between Nu and Wi for different viscosity ratios, β0. The Nu for β0 < 1 is larger than that for β0 = 1, i.e., the heat transfer performance of Giesekus-fluid-based nanofluids was better than that of Newtonian-fluid-based nanofluids, which is in agreement with previous conclusions that viscoelastic fluid can enhance heat transfer rates compared to Newtonian fluids, e.g., [34]. Nu decreases with increasing β0, and the reducing rate was almost the same when 0 < β0 < 1, showing the enhancement of heat transfer depends linearly on the viscosity ratio.

5.2.3. Effects of Reynolds Number and Mobility Parameters

The relationship between Nu and Wi for different Re is shown in Figure 10, where Nu increases with increasing Re. The relationship between the growth rate of Nu and increasing rate of Re was nonlinear. As the flow was laminar in the range of Re in this paper, there are two reasons for why Nu was proportional to Re: (1) the laminar sublayer becomes thinner and shear stress subsequently increases with increasing Re; (2) more heat is carried by a faster moving fluid than a slower one.
Figure 11 shows the relationship between Nu and Wi for different α. Nu increases with the increase in α, and the growth rate of Nu is almost the same when 0 < α < 0.5, showing the increase in heat transfer depends linearly on α. The larger α is, the more obvious the shear-thinning effect is and the thinner the viscous sublayer is, leading to increases in heat transfer. In Equation (5), α = 0 corresponds to Oldroyd-B fluid; therefore, the heat transfer performance of Giesekus-fluid-based nanofluids is better than that of Oldroyd-B-based nanofluids.

5.3. Performance Evaluation Criterion

As shown above, the heat transfer performance of Giesekus-fluid-based nanofluids is better than that of Giesekus fluid, while the friction factor in the former is larger than that in the latter. Therefore, we need to balance the increase in heat transfer and enhancement of consumed power. An energy performance evaluation criterion (PEC) was employed to express the ratio of heat flow rate to required consumed power [35]:
PEC = ( T o u t T in ) D / 2 D / 2 2 π r u z ( ρ C p ) nf d r Δ P D / 2 D / 2 2 π r u z d r
where Tin and Tout are the temperatures at the inlet and outlet of the pipe, respectively.

5.3.1. Effects of Various Factors on the PEC

The ratio of PEC for Giesekus-fluid-based nanofluids (PECnf) to that for Giesekus fluid (PECf) is expressed as PECnf/PECf. Figure 12, Figure 13, Figure 14 and Figure 15 show PECnf/PECf as a function of Wi for different Φ, β0, Re, and α, respectively. At a small Wi (Wi < 3~6), the difference in Cf/CfNew is larger than that in Nu between Giesekus-fluid-based nanofluids and Giesekus fluid; hence PECnf/PECf is less than 1. On the contrary, PECnf/PECf is larger than 1 at a large Wi because the difference in Nu is larger than that in Cf/CfNew. PECnf/PECf increases with increasing Wi, Φ, Re, and α, and with decreasing β0. Therefore, it is more effective to use nanofluids at a large Wi, Φ, Re, and α, but small β0.

5.3.2. Correlation Model

In Figure 12, Figure 13, Figure 14 and Figure 15, PECnf/PECf is proportional to Wi, Φ, Re, and α, but inversely proportional to β0. In order to describe the impact of related parameters on PECnf/PECf, we built a correlation model relating PECnf/PECf to related parameters. Firstly, related parameters were combined into a dimensionless parameter:
ξ = Wi Φ Re α β 0
We established a relationship between PECnf/PECf and ξ based on formula (24) and numerical data in Figure 12, Figure 13, Figure 14 and Figure 15 as:
PEC nf / PEC f = 0.17632 e ( ξ / 171.88918 ) + 1.12805
Figure 16 shows the numerical data in Figure 12, Figure 13, Figure 14 and Figure 15 and formula (25) on a fitted curve.

6. Conclusions

The friction factor and heat transfer of Giesekus-fluid-based nanofluids in a pipe flow were investigated in order to understand the function of non-Newtonian-fluid-based nanofluids as a heat transfer medium. The impacts of Wi, Φ, β0, Re, and α on Cf/CfNew, Nu, and PECnf/PECf were discussed. The main conclusions are summarized as follows:
(1) Cf/CfNew is larger for Giesekus-fluid-based nanofluids than for Newtonian fluid, and there exists a maximum Cf/CfNew under certain rheological and physical properties of nanofluids and flow conditions. As Wi increases, Cf/CfNew first increases and then decreases. Nanoparticles in Giesekus-fluid-based nanofluids lead to an increase in Cf/CfNew. An increase in Φ results in an increase in Cf/CfNew, and the effect of Φ on Cf/CfNew is nonlinear. Cf/CfNew is proportional to Wi in the absence of solvent viscosity, and the changing trend of Cf/CfNew with Wi depends on β0. A large β0 should be selected in order to reduce friction factor. Cf/CfNew is reduced with increasing Re, but the relationship between the reducing rate of Cf/CfNew and increasing rate of Re is nonlinear. Increasing α leads to a decrease in Cf/CfNew, and Giesekus-fluid-based nanofluids have a smaller Cf/CfNew than Oldroyd-B-based nanofluids.
(2) As Wi increases, Nu increases rapidly when 0 < Wi < 2 and slowly when 2 < Wi < 8, and fluid elasticity has little effect on Nu when Wi exceeds a certain value. Adding nanoparticles into fluid can enhance heat transfer. Nu increases with increasing Φ, and the increasing rate of Nu decreases with the increase in Φ. Giesekus-fluid-based nanofluids have a better heat transfer effect than Newtonian-fluid-based nanofluids. Nu decreases with increasing β0, and the enhancement of Nu depends linearly on β0. Nu increases with increasing Re, and the relationship between the growth rate of Nu and increase rate of Re is nonlinear. Nu increases with an increase in α, and the enhancement of Nu depends linearly on α.
(3) PECnf/PECf increases with increasing Wi, Φ, Re, and α, and with decreasing β0. It is more effective to use nanofluids at a large Wi, Φ, Re, and α, but small β0. Finally, the formula of PECnf/PECf as a function of Wi, Φ, β0, Re, and α was derived based on numerical data.

Author Contributions

Conceptualization, J.L. and W.L.; methodology, W.L. and H.Y.; software, W.L. and H.Y.; validation, W.L. and H.Y.; writing, W.L. and H.Y.; resources W.L. and J.L.; review, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant no. 12132015).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available.

Conflicts of Interest

The authors declare no conflict of interest regarding the publication of this paper.

References

  1. Cruz, D.A.; Coelho, P.M.; Alves, M.A. A simplified method for calculating heat transfer coefficients and friction factors in laminar pipe flow of non-newtonian fluids. J. Heat Transf. Trans. ASME 2012, 134, 091703. [Google Scholar] [CrossRef]
  2. Chang, X.; Zhou, J.; Guo, Y.T.; He, S.M.; Wang, L.; Chen, Y.L.; Tang, M.; Jian, R. Heat transfer behaviors in horizontal wells considering the effects of drill pipe rotation, and hydraulic and mechanical frictions during drilling procedures. Energies 2018, 11, 2414. [Google Scholar] [CrossRef] [Green Version]
  3. Martorana, P.; Bayer, I.S.; Steele, A.; Loth, E. Effect of graphite and carbon nanofiber additives on the performance efficiency of a gear pump driven hydraulic circuit using ethanol. Ind. Eng. Chem. Res. 2010, 49, 11363–11368. [Google Scholar] [CrossRef]
  4. Peng, Y.P.; Zahedidastjerdi, A.; Abdollahi, A.; Amindoust, A.; Bahrami, M.; Karimipour, A.; Goodarzi, M. Investigation of energy performance in a u-shaped evacuated solar tube collector using oxide added nano-particles through the emitter, absorber and transmittal environments via discrete ordinates radiation method. J. Therm. Anal. Calorim. 2020, 139, 2623–2631. [Google Scholar] [CrossRef]
  5. Tian, Z.; Abdollahi, A.; Shariati, M.; Amindoust, A.; Arasteh, H.; Karimipour, A.; Goodarzi, M.; Bach, Q.V. Turbulent flows in a spiral double-pipe heat exchanger optimal performance conditions using an enhanced genetic algorithm. Int. J. Numer. Methods Heat Fluid Flow 2020, 30, 39–53. [Google Scholar] [CrossRef]
  6. Sandeep, N.; Malvandi, A. Enhanced heat transfer in liquid thin film flow of non-newtonian nanofluids embedded with graphene nanoparticles. Adv. Powder Technol. 2016, 27, 2448–2456. [Google Scholar] [CrossRef]
  7. Sulaiman, M.; Ali, A.; Islam, S. Heat and mass transfer in three-dimensional flow of an Oldroyd-B nanofluid with gyrotactic micro-organisms. Math. Probl. Eng. 2018, 2018, 6790420. [Google Scholar] [CrossRef] [Green Version]
  8. Sandeep, N.; Sulochana, C. Momentum and heat transfer behavior of Jeffrey, Maxwell and Oldroyd-B nanofluids past a stretching surface with non-uniform heat source/sink. Ain Shams Eng. J. 2018, 9, 517–524. [Google Scholar] [CrossRef] [Green Version]
  9. Aziz, A.; Muhammad, T.; Alsaedi, A.; Hayat, T. An optimal study for 3D rotating flow of Oldroyd-B nanofluid with convectively heated surface. J. Braz. Soc. Mech. Sci. Eng. 2019, 41, 236. [Google Scholar] [CrossRef]
  10. Berberovic, E.; Bikic, S. Computational study of flow and heat transfer characteristics of EG-Si3N4 nanofluid in laminar flow in a pipe in forced convection regime. Energies 2020, 13, 74. [Google Scholar] [CrossRef] [Green Version]
  11. Sandeep, N.; Kumar, B.R.; Kumar, M.S.J. A comparative study of convective heat and mass transfer in non-Newtonian nanofluid flow past a permeable stretching sheet. J. Mol. Liq. 2015, 212, 585–591. [Google Scholar] [CrossRef]
  12. Shaikh, S.; Lafdi, K.; Ponnappan, R. Thermal conductivity improvement in Carbon nanoparticle doped PAO oil: An experimental study. J. Appl. Phys. 2007, 101, 064302. [Google Scholar] [CrossRef] [Green Version]
  13. Nelson, I.C.; Banerjee, D.; Ponnappan, R. Flow loop experiments using polyalphaolefin nanofluids. J. Heat Transf. 2009, 23, 752–761. [Google Scholar] [CrossRef]
  14. Yu, L.; Liu, D.; Botz, F. laminar convective heat transfer of alumina-polyalphaolefin nanofluids containing spherical and non-spherical nanoparticles. Exp. Therm. Fluid Sci. 2012, 37, 72–83. [Google Scholar] [CrossRef]
  15. Yang, J.C.; Li, F.C.; Cai, W.H.; Zhang, H.N.; Yu, B. Direct numerical simulation of viscoelastic-fluid-based nanofluid turbulent channel flow with heat transfer. Chin. Phys. B 2015, 24, 084401. [Google Scholar] [CrossRef]
  16. Wang, Y.; Wang, Y.; Cheng, Z. Direct numerical simulation of gas-liquid drag-reducing cavity flow by the voset method. Polymer 2019, 11, 596. [Google Scholar] [CrossRef] [Green Version]
  17. D’Avino, G.; Maffettone, L. Particle dynamics in viscoelastic liquids. J. Non-Newton. Fluid Mech. 2015, 215, 80–104. [Google Scholar] [CrossRef]
  18. Batchelor, G.K. The Effect of brownian motion on the bulk stress in a suspension of spherical particles. J. Fluid Mech. 1977, 83, 97–117. [Google Scholar] [CrossRef]
  19. Zhang, X.; Gu, H.; Fujii, M. Effective thermal conductivity and thermal diffusivity of nanofluids containing spherical and cylindrical nanoparticles. Exp. Therm. Fluid Sci. 2007, 31, 593–599. [Google Scholar] [CrossRef]
  20. Barrett, J.C.; Webb, N.A. A Comparison of some approximate methods for solving the aerosol general dynamic equation. J. Aerosol Sci. 1998, 29, 31–39. [Google Scholar] [CrossRef]
  21. Barthelmes, G.; Pratsinis, S.E.; Buggisch, H. Particle size distributions and viscosity of suspensions undergoing shear-induced coagulation and fragmentation. Chem. Eng. Sci. 2003, 58, 2893–2902. [Google Scholar] [CrossRef]
  22. Marchisio, D.L.; Vigil, R.D.; Fox, R.O. Implementation of the quadrature method of moments in CFD codes for aggregation-breakage problems. Chem. Eng. Sci. 2003, 58, 3337–3351. [Google Scholar] [CrossRef]
  23. Yu, M.Z.; Lin, J.Z.; Chan, T.L. A new moment method for solving the coagulation equation for particles in brownian motion. Aerosol Sci. Technol. 2008, 42, 705–713. [Google Scholar] [CrossRef]
  24. Yu, M.Z.; Lin, J.Z. Binary homogeneous nucleation and growth of water-sulfuric acid nanoparticles using a TEMOM model. Int. J. Heat Mass Transf. 2010, 53, 635–644. [Google Scholar] [CrossRef]
  25. Vachagina, E.; Dushin, N.; Kutuzova, E.; Kadyirov, A. Exact solution for viscoelastic flow in pipe and experimental validation. Polymers 2022, 14, 334. [Google Scholar] [CrossRef]
  26. Khatibi, A.M.; Mirzazadeh, M.; Rashidi, F. Forced convection heat transfer of Giesekus viscoelastic fluid in pipes and channels. Heat Mass Transf. 2010, 46, 405–412. [Google Scholar] [CrossRef]
  27. Housiadas, K.D.; Beris, A.N. On the skin friction coeffificient in viscoelastic wall-bounded flows. Int. J. Heat Fluid Flow 2013, 42, 49–67. [Google Scholar] [CrossRef]
  28. Bao, F.B.; Lin, J.Z. Burnett simulation of gas flow and heat transfer in micro Poiseuille flow. Int. J. Heat Mass Transf. 2008, 51, 4139–4144. [Google Scholar] [CrossRef]
  29. Azmi, W.H.; Sharma, K.V.; Sarma, P.K.; Septiadi, W.N. Experimental determination of turbulent forced convection heat transfer and friction factor with SiO2 nanofluid. Exp. Therm. Fluid Sci. 2013, 51, 103–111. [Google Scholar] [CrossRef] [Green Version]
  30. Lin, J.Z.; Xia, Y.; Ku, X.K. Pressure drop and heat transfer of nanofluid in turbulent pipe flow considering particle coagulation and breakage. J. Heat Transf. 2014, 136, 111701. [Google Scholar] [CrossRef]
  31. Zhang, P.J.; Lin, J.Z.; Ku, X.K. Friction factor and heat transfer of nanofluid in the turbulent flow through a 90° bend. J. Hydrodyn. 2021, 33, 1105–1118. [Google Scholar] [CrossRef]
  32. Melhi, S.A.; Filali, A.; Khezzar, L.; Alshehhi, M. Flow and heat transfer of a Giesekus fluid in plane and 3D ducts. Heat Transf.-Asian Res. 2017, 46, 1380–1398. [Google Scholar] [CrossRef]
  33. Abbasian Arani, A.A.; Amani, J. Experimental investigation of diameter effect on heat transfer performance and pressure drop of TiO2–water nanoflfluid. Exp. Therm. Fluid Sci. 2013, 44, 520–533. [Google Scholar] [CrossRef]
  34. Khezzar, L.; Filali, A.; AlShehhi, A.M. Flow and heat transfer of FENE-P fluids in ducts of various shapes: Effffect of Newtonian solvent contribution. J. Non-Newton Fluid. 2014, 207, 7–20. [Google Scholar] [CrossRef]
  35. Ferrouillat, S.; Bontemps, A.; Poncelet, O.; Soriano, O.; Gruss, J.A. Inflfluence of nanoparticle shape factor on convective heat transfer and energetic performance of water-based SiO2 and ZnO nanoflfluids. Appl. Therm. Eng. 2013, 51, 839–851. [Google Scholar] [CrossRef]
Figure 1. Flow in a pipe and cylindrical coordinate system.
Figure 1. Flow in a pipe and cylindrical coordinate system.
Energies 15 03234 g001
Figure 2. Axial velocity profile along the radial direction. ⸺●⸺: present numerical result;: analytical result [25]; : experimental result [25].
Figure 2. Axial velocity profile along the radial direction. ⸺●⸺: present numerical result;: analytical result [25]; : experimental result [25].
Energies 15 03234 g002
Figure 3. Variation in Nu with Wi (Br = 0.1). ⸺●⸺: present numerical result; : theoretical solution [26].
Figure 3. Variation in Nu with Wi (Br = 0.1). ⸺●⸺: present numerical result; : theoretical solution [26].
Energies 15 03234 g003
Figure 4. Friction factor vs. modified Weissenberg number for different particle concentrations (β0 = 0.2, Re = 2.5, α = 0.2). : Φ = 0% (Giesekus fluid) [27]; : Φ = 0% (Giesekus fluid); ⸺●⸺: Φ = 0.5%; : Φ = 1.0%; : Φ = 2.0%; and :Φ = 3.0%.
Figure 4. Friction factor vs. modified Weissenberg number for different particle concentrations (β0 = 0.2, Re = 2.5, α = 0.2). : Φ = 0% (Giesekus fluid) [27]; : Φ = 0% (Giesekus fluid); ⸺●⸺: Φ = 0.5%; : Φ = 1.0%; : Φ = 2.0%; and :Φ = 3.0%.
Energies 15 03234 g004
Figure 5. Friction factor vs. modified Weissenberg for different viscosity ratios (Φ =1.0%, Re = 2.5, α = 0.2). : β0 =1.0; : β0 = 0.9; : β0 = 0.3; : β0 = 0.2; ⸺: β0 = 0.1; and : β0 = 0.
Figure 5. Friction factor vs. modified Weissenberg for different viscosity ratios (Φ =1.0%, Re = 2.5, α = 0.2). : β0 =1.0; : β0 = 0.9; : β0 = 0.3; : β0 = 0.2; ⸺: β0 = 0.1; and : β0 = 0.
Energies 15 03234 g005
Figure 6. Friction factor vs. modified Weissenberg number for different Reynolds numbers (β0 = 0.2, Φ = 1.0%, α = 0.2). : Re = 500; : Re = 100; ⸺●: Re = 10; : Re = 2.5; : Re = 1; and : Re = 0.5.
Figure 6. Friction factor vs. modified Weissenberg number for different Reynolds numbers (β0 = 0.2, Φ = 1.0%, α = 0.2). : Re = 500; : Re = 100; ⸺●: Re = 10; : Re = 2.5; : Re = 1; and : Re = 0.5.
Energies 15 03234 g006
Figure 7. Friction factor vs. modified Weissenberg number for different mobility parameters (Re = 2.5, Φ = 1.0%, β0 = 0.2). : α = 0.5; : α = 0.4; : α = 0.3; : α = 0.2; : α = 0.1; and : α = 0.
Figure 7. Friction factor vs. modified Weissenberg number for different mobility parameters (Re = 2.5, Φ = 1.0%, β0 = 0.2). : α = 0.5; : α = 0.4; : α = 0.3; : α = 0.2; : α = 0.1; and : α = 0.
Energies 15 03234 g007
Figure 8. Nusselt number vs. Weissenberg number for different particle concentrations (β0 = 0.2, Re = 2.5, α = 0.2). : Φ = 0% (Giesekus fluid) [32]; : Φ = 0% (Giesekus fluid); ⸺●⸺: Φ = 0.5%; : Φ = 1.0%; : Φ = 2.0%; and : Φ = 3.0%.
Figure 8. Nusselt number vs. Weissenberg number for different particle concentrations (β0 = 0.2, Re = 2.5, α = 0.2). : Φ = 0% (Giesekus fluid) [32]; : Φ = 0% (Giesekus fluid); ⸺●⸺: Φ = 0.5%; : Φ = 1.0%; : Φ = 2.0%; and : Φ = 3.0%.
Energies 15 03234 g008
Figure 9. Nusselt number vs. Weissenberg number for different viscosity ratios (Φ = 1.0%, Re = 2.5, α = 0.2). : β0 = 1.0; : β0 = 0.9; : β0 = 0.3; : β0 = 0.2; ⸺●: β0 = 0.1; and : β0 = 0.
Figure 9. Nusselt number vs. Weissenberg number for different viscosity ratios (Φ = 1.0%, Re = 2.5, α = 0.2). : β0 = 1.0; : β0 = 0.9; : β0 = 0.3; : β0 = 0.2; ⸺●: β0 = 0.1; and : β0 = 0.
Energies 15 03234 g009
Figure 10. Nusselt number vs. Weissenberg number for different Reynolds numbers (β0 = 0.2, Φ = 1.0%, α = 0.2). : Re = 0.5; : Re = 1; : Re = 2.5; : Re = 10; : Re = 100; and : Re = 500.
Figure 10. Nusselt number vs. Weissenberg number for different Reynolds numbers (β0 = 0.2, Φ = 1.0%, α = 0.2). : Re = 0.5; : Re = 1; : Re = 2.5; : Re = 10; : Re = 100; and : Re = 500.
Energies 15 03234 g010
Figure 11. Nusselt number vs. Weissenberg number for different mobility parameters (Re = 2.5, Φ = 1.0%, β0 = 0.2). : α = 0; : α = 0.1; : α = 0.2; : α = 0.3; : α = 0.4; and : α = 0.5.
Figure 11. Nusselt number vs. Weissenberg number for different mobility parameters (Re = 2.5, Φ = 1.0%, β0 = 0.2). : α = 0; : α = 0.1; : α = 0.2; : α = 0.3; : α = 0.4; and : α = 0.5.
Energies 15 03234 g011
Figure 12. PECnf/PECf vs. Weissenberg number for different particle concentrations (β0 = 0.2, Re = 2.5, α = 0.2). : Φ = 0% (Giesekus fluid); : Φ = 0.5% (Giesekus fluid); : Φ = 1.0%; : Φ = 2.0%; and : Φ = 3.0%.
Figure 12. PECnf/PECf vs. Weissenberg number for different particle concentrations (β0 = 0.2, Re = 2.5, α = 0.2). : Φ = 0% (Giesekus fluid); : Φ = 0.5% (Giesekus fluid); : Φ = 1.0%; : Φ = 2.0%; and : Φ = 3.0%.
Energies 15 03234 g012
Figure 13. PECnf/PECf vs. Weissenberg number for different viscosity ratios (Φ = 1.0%, Re = 2.5, α = 0.2). : β0 = 1.0; : β0 = 0.9; : β0 = 0.3; : β0 = 0.2; : β0 = 0.1; and : β0 = 0.
Figure 13. PECnf/PECf vs. Weissenberg number for different viscosity ratios (Φ = 1.0%, Re = 2.5, α = 0.2). : β0 = 1.0; : β0 = 0.9; : β0 = 0.3; : β0 = 0.2; : β0 = 0.1; and : β0 = 0.
Energies 15 03234 g013
Figure 14. PECnf/PECf vs. Weissenberg number for different Reynolds numbers (β0 = 0.2, Φ = 1.0%, α = 0.2). : Re = 0.5; : Re = 1; : Re = 2.5; : Re = 10; : Re = 100; and : Re = 500.
Figure 14. PECnf/PECf vs. Weissenberg number for different Reynolds numbers (β0 = 0.2, Φ = 1.0%, α = 0.2). : Re = 0.5; : Re = 1; : Re = 2.5; : Re = 10; : Re = 100; and : Re = 500.
Energies 15 03234 g014
Figure 15. PECnf/PECf vs. Weissenberg number for different mobility parameters (Re = 2.5, Φ =1.0%, β0 = 0.2). : α = 0; : α = 0.1; : α = 0.2; : α = 0.3; : α = 0.4; and : α = 0.5.
Figure 15. PECnf/PECf vs. Weissenberg number for different mobility parameters (Re = 2.5, Φ =1.0%, β0 = 0.2). : α = 0; : α = 0.1; : α = 0.2; : α = 0.3; : α = 0.4; and : α = 0.5.
Energies 15 03234 g015
Figure 16. Relationship between PECnf/PECf and ξ. : numerical data; ⸺⸺: Formula (25).
Figure 16. Relationship between PECnf/PECf and ξ. : numerical data; ⸺⸺: Formula (25).
Energies 15 03234 g016
Table 1. Values of M0 when changing grid points.
Table 1. Values of M0 when changing grid points.
r × θ × SM0r × θ × SM0r × θ × SM0
112 × 32 × 2561.08669128 × 24 × 2561.08664128 × 32 × 2161.08660
120 × 32 × 2561.08651128 × 28 × 2561.08648128 × 32 × 2361.08647
128 × 32 × 2561.08636128 × 32 × 2561.08636128 × 32 × 2561.08636
136 × 32 × 2561.08628128 × 36 × 2561.08630128 × 32 × 2761.08631
144 × 32 × 2561.08622128 × 40 × 2561.08626128 × 32 × 2961.08628
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, W.; Yang, H.; Lin, J. Friction Factor and Heat Transfer of Giesekus-Fluid-Based Nanofluids in a Pipe Flow. Energies 2022, 15, 3234. https://doi.org/10.3390/en15093234

AMA Style

Lin W, Yang H, Lin J. Friction Factor and Heat Transfer of Giesekus-Fluid-Based Nanofluids in a Pipe Flow. Energies. 2022; 15(9):3234. https://doi.org/10.3390/en15093234

Chicago/Turabian Style

Lin, Wenqian, Hailin Yang, and Jianzhong Lin. 2022. "Friction Factor and Heat Transfer of Giesekus-Fluid-Based Nanofluids in a Pipe Flow" Energies 15, no. 9: 3234. https://doi.org/10.3390/en15093234

APA Style

Lin, W., Yang, H., & Lin, J. (2022). Friction Factor and Heat Transfer of Giesekus-Fluid-Based Nanofluids in a Pipe Flow. Energies, 15(9), 3234. https://doi.org/10.3390/en15093234

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