Heat Transfer Capability of (Ethylene Glycol + Water)-Based Nanofluids Containing Graphene Nanoplatelets: Design and Thermophysical Profile

This research aims at studying the stability and thermophysical properties of nanofluids designed as dispersions of sulfonic acid-functionalized graphene nanoplatelets in an (ethylene glycol + water) mixture at (10:90)% mass ratio. Nanofluid preparation conditions were defined through a stability analysis based on zeta potential and dynamic light scattering (DLS) measurements. Thermal conductivity, dynamic viscosity, and density were experimentally measured in the temperature range from 283.15 to 343.15 K and nanoparticle mass concentrations of up to 0.50% by using a transient plate source, a rotational rheometer, and a vibrating-tube technique, respectively. Thermal conductivity enhancements reach up to 5% without a clear effect of temperature while rheological tests evidence a Newtonian behavior of the studied nanofluids. Different equations such as the Nan, Vogel-Fulcher-Tamman (VFT), or Maron-Pierce (MP) models were utilized to describe the temperature or nanoparticle concentration dependences of thermal conductivity and viscosity. Finally, different figures of merit based on the experimental values of thermophysical properties were also used to compare the heat transfer capability and pumping power between nanofluids and base fluid.


Background
As energy transport processes are integrated in most industrial areas, the development of more efficient and compact heat transfer equipment has focused a lot of attention during the last decades. Once different approaches such as material modifications, the use of extended surfaces, or the optimization of process parameters have been extensively exploited, many research activities are focused now on improving the weak heat-transfer abilities of typical conventional fluids such as water (W), ethylene glycol (EG), or engine oils [1,2]. In this sense, the enhancement of thermal transport properties of heat transfer fluids by dispersing high-thermal conductive nanomaterials, also known as nanofluids, has become a prominent research avenue [3,4].
A wide variety of nanoadditives, including metallic, organic, and non-organic materials, has been utilized to engineer nanofluids. Among them, carbon nanostructures are those which seem to exhibit the highest potential [5]. Within the family of graphite carbon allotropes, the remarkable mechanical, structural, thermal, and electrical properties of graphene [6][7][8] have raised great interest since this material was discovered by Novoselov et al. [9] in 2004. Graphene structure is ideally envisaged as a single-atom-thick sheet of hexagonally arranged, sp 2bonded carbon atoms tightly packed into a honeycomb lattice [10,11]. This two-dimensional material is available commercially in the form of several-layer stacks (usually more than 10 layers) which are known as graphene nanoplatelets or nanosheets (GnPs).
GnPs combine the advantages of both single-layer graphene and highly ordered graphitic carbon [12,13]. Thus, the thermal conductivity of graphene nanoplatelets has been reported to be much higher than the value presented by other carbon allotropes like multi-wall nanotubes (MWNT), single-wall nanotubes (SWNT), or diamond [8,14]. In addition, graphene nanosheets possess larger surfaces areas compared to nanotubes or other nanoparticles, allowing bigger contact areas/interfaces with the base fluid. This superior contact could lead to a reduction of Kapitza resistance at the graphene-fluid interface and, consequently, help to improve the effective thermal conductivity of the nanofluid [15]. Another upside of graphene nanomaterials is the relatively easy and cost-effective production at a large scale [13,16,17]. Hence, graphene appears to meet all the conditions to develop new nanofluids with improved thermal properties. Nevertheless, unfortunately, graphene is hydrophobic and consequently it cannot be dispersed in water or some polar organic solvents for a long time without agglomerating [18,19]. In order to increase the dispersibility, the interactions between graphene nanostructures and aqueous/organic solvents can be improved by means of non-covalent (using surfactants) or covalent (adding hydrophobic or hydrophilic groups on the high energy features such as edges of graphene nanosheets) functionalizations [16].
In water-and ethylene glycol-based nanofluids, the surfactants predominating in the literature for non-covalent options are sodium dodecyl benzene sulphonate (SDBS) [16,[20][21][22], polyvinyl alcohol (PVA) [23], and gum Arabic [14,22], while the oxidation following the original version [24,25] or modifications [26] of Hummers method are the most common with covalent functionalizations. The main downside of using non-covalent methods is that the addition of a surfactant considerably increases the viscosity of the nanofluid, leading to higher pressure drops in thermal equipment. As an example, Amiri et al. [16] found that the viscosity enhancement of non-covalent graphene/ water nanofluids synthesized using SDBS doubles the values of covalent functionalized nanofluids at the same nanoparticle concentration. Regarding covalent functionalizations, although Paredes et al. [27] reported that good dispersions of graphene oxide (GO) can be obtained not only in water but also in other organic solvents such as ethylene glycol, acid treatment can lead to the formation of defect sites within the systematically arranged conjugated graphene structure of the nanosheets [28][29][30]. As these defects could reduce the thermal conductivity of graphene by an order of magnitude or more below its intrinsic value, some reaction conditions such as acid content, temperature, or time of exposure must be optimized in order to minimize overoxidation and subsequent defect formation [11]. Another peculiar property of graphene oxide is the high natural acidity of its aqueous solutions due to the generation of hydrogen cations resulting from GO-water interactions [31,32]. Thus, a special attention must be paid to pH control in order to avoid possible damages to metallic elements of the installation [33].
During the last decades, thermal conductivity (k) of nanofluids has received a lot of attention due to the high influence of this property on the heat transfer performance of thermal facilities [34,35]. However, other thermophysical properties, such as dynamic viscosity (η), density (ρ), or even specific heat capacity (c p ), should also be taken into account in order to assess whether the replacement of the conventional fluid with the new nanofluid would be really beneficial, as well as to make technical calculations of thermal facilities. Dynamic viscosity or the rheological characteristics in general have a critical effect on the type of flow and consequently on the heat transfer and the necessary pumping power, for example [36]. Although the complex behavior of nanofluids precludes a generalization of their thermophysical properties, the addition of nanoparticles usually leads to higher thermal conductivities, viscosities, and densities as well as to lower specific heat capacities. The improvement of some properties and the worsening of others hinders the choice of the fluid with the best features, and so it is necessary to use figures of merit (FoMs) such as the Mouromtseff number (Mo) [37] in order to select the fluid with better heat transfer capabilities.
Over the last years, an important number of studies about graphene nanofluids have been performed using water [1,3,13,[38][39][40][41] or ethylene glycol [12,20,31,39,40] as base fluids. However, many industrial facilities do not utilize these pure compounds as heat transfer media but their mixtures in order to combine the advantages of both compounds, such as the better thermal conductivity of water or the lower freezing points of glycols. To our knowledge, only Kole and Dey [15], Amiri et al. [42], and Ijam et al. [43] have analyzed the thermophysical properties of graphene nanofluids based on (ethylene glycol + water) mixtures. Kole and Dey [15] studied the thermal conductivity and rheological behavior of functionalized hydrogen exfoliated graphene in EG + W (70:30 vol.%) at nanoparticle concentrations between 0.041 and 0.395 vol.%. They found that both properties increase with the addition of graphene, maximum increases reaching 15% for thermal conductivity and 100% for dynamic viscosity. Amiri et al. [42] used a mixture of EG + W (60:40 vol.%) as base fluid to prepare nanofluids at graphene mass concentrations from 0.01 to 0.2%. Thermal conductivity, density, and viscosity were reported to rise up to (65, 0.6, and 4.9)%, respectively, for the highest concentration. Contrarily, specific heat capacity decreases up to 5%. Ijam et al. [43] studied dispersions of graphene oxide nanosheets in EG + W (40:60 vol.%) finding thermal conductivity and viscosity increases at the highest concentration (0.10 wt.%) of up to 10.5 and 35%, respectively.
In this work, sulfonic acid-functionalized graphene nanoplatelets were dispersed in a (ethylene glycol + water) mixture at (10:90)% mass ratio to prepare nanofluids at (0.10, 0.25, and 0.50)% nanoparticle mass fractions, which correspond to 0.00045, 0.00112, and 0.00225 nanoparticle volume fractions, respectively. Volume fractions were calculated by using a density value of 2.25 g cm −3 [28] for graphene oxide and the experimental densities obtained in this work for base fluid. Preparation conditions such as sonication time or pH value were optimized by studying the zeta potential and variation of the average nanoparticle size with time. Thermal conductivity, dynamic viscosity, and density of nanofluids and base fluid were obtained experimentally, and the influences of temperature and nanoparticle concentration on these three properties were analyzed. Finally, the modifications of heat transfer performance and pumping power were assessed from the thermophysical properties here obtained by means of different figures of merit.

Materials
Sulfonic acid-functionalized graphene oxide nanoplatelets (GOnPs) were provided by NanoInnova Technologies S.L. (Madrid, Spain). The base fluid consists of a (ethylene glycol + water) mixture at (90:10)% mass ratio. Ethylene glycol was purchased from Sigma-Aldrich with a mass purity of 99.5% while Milli-Q Grade water was produced with a resistivity of 18.2 MΩ cm at 298 K by means of a Millipore system (Billerica, MA, USA). An ammonium hydroxide solution (30-33% NH 3 in water) from Sigma-Aldrich was used for the modification of the pH value.
A Sartorius analytical balance with an uncertainty of 0.0001 g was utilized to weigh the reagents.

Powder Characterization, Nanofluid Preparation, and Stability Analysis
The morphology of the dry powder was studied through scanning electron microscopy analyses conducted on a JEOL JSM-6700F field emission gun-SEM (JEOL, Tokyo, Japan) working at an accelerator voltage of 20 kV in backscattering electron image (Yttrium Aluminum Garnet type detector). SEM samples were prepared by depositing a drop of a GOnP dispersion in analytical grade methanol on the top of a silica support and drying it under atmospheric conditions. A typical SEM image of GOnP additive is presented in Fig. 1. GO particles exhibit a plate-like shape of up to some micrometers, rough surfaced with tiny wrinkles over the whole surface and mountainous peaks similar to those presented by Geng et al. [44]. The edges are not flat and smooth but rounded, which may be due to the synthesis route as also pointed out in [44]. A further characterization of the studied GOnP nanoplatelets can be found in Agromayor et al. [45].
Nanofluids were prepared through a two-step process. Thus, a predetermined amount of graphene nanopowder was first added to the mass of liquid necessary to obtain the desired nanoparticle concentration and then stirred for 120 min. Afterwards, samples were sonicated by means of a CP104 Ultrasonic Bath (CEIA, Italy) working with a maximum sonication power of 200 W and a sonication frequency of 40 kHz. As it was pointed out, one of the disadvantages of graphene oxide nanoplatelets is the high acidity of its aqueous solutions due to the presence of surface acidic groups [31,32]. Thus, in our study, pH values between 2.2 and 2.8 were measured for the nanofluids once the nanomaterial was added. In order to modify the pH conditions, the pH value of dispersions was raised by adding ammonium hydroxide. pH measurements were performed with a PHM 210 standard pH meter (Radiometer Analytical S.A., France) with a pH electrode code 5208 (Crison Instruments SA, Spain). With the aim of optimizing nanofluid design, the influences of the pH value and sonication procedure on nanofluid stability were evaluated through zeta potential and size measurements by using a Zetasizer Nano (Malvern, UK) [46].
Firstly, zeta potential was studied for 0.10 and 0.25% nanoparticle mass concentrations in the pH range from 2.2 to 10. No zeta potential measurements were performed for the highest mass concentration since the intensity detected by the device was outside the optimum working range. Figure 2 shows the pH influence on the zeta potential of graphene nanofluids. As can be observed, zeta potential remains constant around 40-43 mV for pH values up to 5 and then decreases as pH value rises. In order to reach a compromise between high zeta potentials to ensure strong electrical repulsion charges around the particles and appropriate pH values to avoid corrosion issues, a pH = 5 was selected to prepare the nanofluids.
The nanoparticle size distribution of the nanofluids was determined by measuring the random changes in the intensity of light scattered from the dispersions at 298.15 K and with a scattering angle of 173°. This technique was used to select the optimum sonication time determining "apparent" size measurements for nanofluids prepared with 0.10 wt.% GOnP concentration and sonication times from 0 to 300 min. It should be pointed out that dynamic light scattering (DLS) measurements are based on the assumption that particles are spherical while the studied nanoadditives are sheet-like shaped. The DLS results indicate that the graphene nanoplatelets of the fresh nanofluids exhibit a trimodal distribution with apparent sizes from some nanometers to 4 μm.
This polydispersity of graphene nanofluids was previously reported in literature [24,47]. Width and medium value of the different peaks of the size distribution reduce as sonication time increases. However, this reduction mainly occurs at times between 0 and 240 min. while only slight variations were found between samples prepared at 240 and 300 min. Thus, a sonication time of 300 min was selected to prepare the nanofluids at the three studied nanoadditive concentrations. Subsequently, the evolution of nanofluid stability with the time elapsed after preparation was studied following the procedure defined by Fedele et al. [46]. Thus, two samples of each nanoparticle concentration were put in two different cuvettes and their apparent sizes were studied by DLS for a month. The sample of the first cuvette was measured almost every day without shaking the fluid (in static conditions) in order to evaluate the changes in the size distribution due to natural sedimentation, while the second sample was measured after being manually shaken for some seconds to measure apparent size distribution after mechanically recovering settled particles. The apparent size distributions of the 0.10 wt.% graphene concentration measured just after preparation and the 28th day are plotted in Fig. 3. As it can be observed, in the case of the static samples, the main peak slightly moves to the left and the 4000-nm peak disappears while the 50-nm and 400-nm peaks move to the right for the shaken sample. This indicates a partial agglomeration of nanoadditives and sedimentation (especially of the largest nanoplatelets) under static conditions.

Experimental Methods
Thermal conductivities were determined at temperatures ranging from 293 to 343 K by means of a TPS 2500 S Hot Disk Thermal Analyzer® (Hot Disk AB, Sweden), which is appropriate to perform measurements of different materials in the range from 0.005 to 500 W m −1 K. This instrument was used together with a 7877 sensor based on the transient plate source (TPS) and consisting of a double spiral of 2-mm radius made of thin nickel wire. The probe, which works as both heat source and temperature sensor, was vertically immersed in the sample so that heat could freely diffuse in all directions. In order to ensure a uniform initial temperature and remove thermal gradients, the box which contains both sample and probe was placed in a thermostatic bath and measurements of thermal conductivity were performed at least with a delay of 15 min between them. Tests were carried out using low thermal powers, 40-55 mW, and a short power input time, 4 s, in order to minimize convection effects. With the aim of checking the measuring procedure in the studied temperature range, the thermal conductivity of bidistilled water was measured and the results were compared with the NIST database [38]. Maximum deviations between the here reported thermal conductivities and literature data [38] are lower than 1.5%, well within the 5% accuracy reported by the manufacturer. Additional information about the experimental device and measuring procedure can be found in Fedele et al. [48].
Rheological behavior was studied by using an AR-G2 magnetic bearing rotational rheometer (TA Instruments, New Castle, USA) with a cone-plate geometry of 1°cone angle and 40-mm diameter. This device allows controlling torques from 0.1 μN m to 200 mN m and normal forces between 0.005 and 50 N. The sample temperature was regulated in the range between 293 and 343 K by using a Peltier plate and an upper heated plate (UHP). Before tests, the rheometer was carefully calibrated at each temperature as further described by Bobbo et al. [49]. The geometry was imposed to a gap of 30 μm and an amount of 0.34 cm 3 was considered optimal for the experiments. Flow curve tests were carried out at constant temperature and shear rates ranging from 80 to 1200 s −1 . In order to assess measurement accuracy, water (a fluid of well-known viscosity) was studied in the entire temperature range and the values were compared with those of the NIST database [38]. Average Absolute Deviations, AADs%, between both data sets are less than 2% which is well within the uncertainty reported by the manufacturer for viscosity measurements, i.e., 5%.
Densities were experimentally obtained with an Anton Paar DMA 4500 vibrating U-tube densimeter, which relates sample density and oscillation period. Temperature was controlled with a resolution of 0.01 K in the range from 293 to 343 K by means of a solid-state thermostat. Calibration was performed using water and air, and uncertainty was estimated to be 5 × 10 −4 g cm −3 [50].

Thermal Conductivity
The thermal conductivity of the base fluid and the three designed nanofluids was measured at ambient pressure in the temperature range from 293 to 344 K with steps of 10 K. The temperature dependence of experimental values is depicted in Fig. 4. To our knowledge, only Melinder [33] reported the thermal conductivities of (ethylene glycol + water) mixture at (10:90)% mass concentration. In addition, Sun et al. [51] and Bohne et al. [52] studied different concentrations of (ethylene glycol + water) system and proposed thermal conductivity correlations as a function of composition and temperature. The obtained values for the base fluid were compared with these existing literature data. AADs% of (2.6, 1.9, and 4.0)% between our experimental values and those reported or calculated by using the equations proposed by Melinder [33], Sun et al. [51], and Bohne et al. [52] were obtained, respectively. For each fluid, the temperature dependence was described by using a second-order polynomial: where the a i adjustable coefficients were obtained with AADs% lower than 0.3% by using a least-square fitting and are gathered in Table 1. Thermal conductivity increases with temperature up to 14.7% throughout the entire studied range. It also increases up to 5% with the addition of nanoparticles, and the enhancements have no temperature dependence as showed in the inset of Fig. 4. The reported temperature-independent thermal conductivity enhancements are in agreement with those found by Yu et al. [53] and Hadadian et al. [31] for GOnP/EG nanofluids or by Kole and Dey [15] for hydrogen exfoliated graphene-(EG + W) nanofluids. During the last decades, several theoretical or semiempirical equations, based on Maxwell model [54] for spherical and well-dispersed particles, have been developed in order to predict or describe the thermal conductivity of nanofluids. A good agreement was found in literature between experimental data and the values obtained using Nan model [55] for different nanofluids in general [56] and for graphene nanoplatelet dispersions in particular [15,20,31,57,58]. Nan et al. [55] generalized Maxwell equation including the effects of particle geometry and finite interfacial resistance by the following expression: where L ii are the geometrical factors which take a value of L 11 = 0 and L 33 = 1 in the case of nanoplatelets [15,20], ϕ is the volumetric fraction of particles, and β ii coefficients are defined as: Based on that model, in-plane thermal conductivity of the nanomaterial, k np , can be obtained by using a least squares fitting of experimental data of the nanofluids, k nf , and base fluid, k bf . Thus, the thermal conductivity obtained by using Eqs. (2)(3) for the studied graphene nanoplatelets is 17 W m −1 K, which is slightly higher than the values reported by Kole and Dey [15], Yu et al. [20], or Hadadian et al. [31] when studied graphene oxide dispersions in (ethylene glycol + water), ethylene glycol and water, respectively. These differences can be  due to the density of defects, size, and/or roughness of graphene oxide, as it was theoretically demonstrated by Nika et al. [59]. The values provided by Nan model are plotted together with the experimental values in Fig. 4, exhibiting a deviation of 0.6%, which is well within the experimental uncertainty.

Dynamic Viscosity
The influence of GOnP concentration on the rheological behavior was analyzed at temperatures ranging from 293.15 to 343.15 K with steps of 10 K. Figure 5a shows the flow curves of the base fluid and different studied nanofluids at 293.15 K. Dynamic viscosities obtained for the base fluid were compared with previous literature data [33,[60][61][62] for the (ethylene glycol + water) mixture at (10:90)% mass ratio. The AADs% between these values range from 1.5 to 3.6%, showing a good agreement.
On the other hand, the linear relationship between shear stress and shear rate indicates that both base fluid and nanofluids exhibit a Newtonian behavior over the studied conditions. This Newtonian behavior coincides with the results found by Kamatchi et al. [57] and Mehrali et al. [58] for graphene oxide-water nanofluids at shear rates higher than 80 s −1 or by Ma et al. [63] for nanofluids based on dimethyl silicone oil, for example. The temperature dependence of dynamic viscosity measurements is plotted in Fig. 5b. As can be observed, dynamic viscosity decreases considerably with temperature and a modification of Andrade's equation, which is also known as the three-coefficient Vogel-Fulcher-Tamman (VFT) model, was used to describe this temperature dependence: where η 0 , D, and T 0 are the adjustable coefficients. D and T 0 parameters are also known as the Angell Strength and Vogel temperature, respectively. The values of these coefficients as well as AADs% and standard deviations between experimental values and data fitted by using this equation are gathered in Table 1. The good description of the temperature dependence, with standard deviations lower than 0.02 mPa s, is also shown in Fig. 5b. The base fluid and studied nanofluids exhibit small values of Angell Strength parameter, which indicates that the fluids present a rapid breakdown of their configurational structure with increasing temperature near and above the glass transition [64,65]. Figure 6 shows the influence of GOnP concentration on viscosity ratio, i.e., η nf /η bf , at several temperatures. As expected, the dynamic viscosity of nanofluids increases with the addition of graphene oxide nanoplatelets. These increases rise with increasing temperature, especially for the 343-K isotherm reaching a maximum increment of 12.6% for the 0.50 wt.% GOnP mass concentration.
The relative viscosity of non-interacting and stationary hard particles in infinitely dilute suspensions can be described by using different equations derived from the pioneering Einstein model [66]. In this sense, Maron and Piece (MP) [67] proposed an equation with the same functional form than that of Krieger-Dougherty relationship [68]. Nevertheless, unlike Krieger-Dougherty model, MP equation does not require the knowledge of intrinsic viscosity, which is not always an easy task when particles are not spherical [69,70]. The Maron and Piece model [67] has already been successfully utilized to model viscosity enhancements of particles with different aspect rations such as fibers [71,72] or platelets/flakes [73,74] and can be written as follows: where ϕ a is the effective volume fraction of nanoparticles which can be reduced to the volume fraction of nanoparticles in absence of aggregates, while ϕ m is the maximum volume packing fraction and can be used as a fitting parameter [70]. In this case, an ADD% of 1.3% and a maximum deviation of 2.8% between experimental and correlated values were obtained finding a value of 0.048 for ϕ m . This maximum volume packing fraction is similar to those values obtained by Fisa [74] and Utracki et al. [75] when studied flaky mica particles.

Density
Density was studied for the base fluid and the three designed nanofluids at temperatures between 293.15 and 343.15 K. Experimental results are depicted in Fig. 7a. A comparison between the values here presented for the base fluid and previous density data [33,60] shows AADs% less than 0.05%. As usually happens with nanofluids, density increases with nanoadditive concentration. These increases slightly rise with temperature ranging from 0.28 to 0.30% for the 0.50% mass concentration, as shown in the inset of the Fig. 7a. For each fluid, the temperature dependence of density was described by using a second-order polynomial: where the b i adjustable coefficients were obtained by using a least-square fitting and are gathered in Table 1.
As shown in Fig. 7a, this equation provides a good density data description with standard deviations lower than 2 × 10 −4 g cm −3 . In addition, isobaric thermal expansivities, α p = -(1/ρ) · (δρ/δT) p , were numerically obtained through the derivatives of the polynomial fits of density. Figure 7b shows α p values for the base fluid and the three studied nanofluids at different temperatures. As can be observed, this property slightly decreases with concentration, diminutions reaching up to 2% for 0.50% mass concentration at 333.15 K. This reduction in α p with nanoparticle concentration agrees with previous nanofluid literature [76,77].

Nanofluid Comparison Based on Thermophysical Properties
As pointed out, different figures of merit based on the thermophysical properties can be used to compare the heat transfer performance or pumping power of nanofluids to those of the base fluid. In the case of fully developed internal laminar flow conditions, a merit criterion to assess the thermal capability is the ratio between viscosity and thermal conductivity enhancements [78][79][80]: Thus, the increase in dynamic viscosity must be lower than four times the thermal conductivity improvement so that the replacement is beneficial in terms of thermal energy. Regarding the turbulent flow conditions, the heat transfer performance can be evaluated through the Mouromtseff number (Mo) [37] which is defined as: Mo ¼ ρ 0:8 ⋅k 0:67 ⋅c 0:33 Higher Mo numbers indicate higher heat transfer capabilities and, consequently, a Mo nf /Mo bf ratio higher than one is desirable.
In a circular tube with a uniform flux at the wall and considering that the mass flow rate of both base fluid and nanofluids is the same, the increment in pumping power can be assessed by using the following expressions for laminar and turbulent flow conditions, respectively [81]: Moderate Ẇ nf /Ẇ bf ratios are needed to avoid increasing the energy consumption due to fluid pumping.
In order to perform these analyses, the experimental values of thermal conductivity, dynamic viscosity, and density were utilized, while specific heat capacities were obtained by using the following weighted average equation [82]: In order to carry out these calculations, the heat capacity data of graphene oxide and base fluid were obtained from literature [28,33]. Figure 8 depicts the temperature dependence of these figures of merit for the different studied nanofluids. Under laminar flow rate conditions only the 0.50 wt.% GnOP concentration is within the benefit zone in the entire studied temperature range. However, for turbulent flow, this same concentration would obtain an improvement only at 303 K (Fig. 8a). In addition, pumping power rises with the concentration of graphene platelets up to 2.4 and 11.9% for laminar and turbulent flow conditions, respectively. Although pumping power increases with temperature under both flow regimens, increases are more appreciable under laminar conditions for which variations reach up to 5% throughout the analyzed range, as can be seen in Fig. 8b.

Conclusions
New nanofluids consisting of GOnP dispersions in an (ethylene glycol + water) mixture were designed following a two-step process. According to zeta potential and dynamic light scattering analyses, dispersions exhibit zeta potentials of around 40 mV and trimodal apparent size distributions. A careful study of the thermodynamic profile of this nanofluid set shows that thermal conductivity, dynamic viscosity, and density increase with nanoplatelet mass concentration up to (5, 12.6, and 0.3)%, respectively. Temperature has a clear effect on the dynamic viscosity and density increases but not on the thermal conductivity enhancements. The temperature or nanoparticle behaviors of the two transport properties were also described by using Nan, Vogel-Fulcher-Tamman (VFT), or Maron and Pierce (MP) models with AADs% of (0.6, 1.6, and 1.3)%, respectively. Based on an analysis of different figures of merit and using the reported thermophysical property data of the nanofluids and base fluid, no significant enhancements in the heat transfer performance would be expected for the studied GOnP concentrations under turbulent conditions, while the pumping power increases would reach up to 2.4 and 11.9% for turbulent and laminar flow conditions, respectively.