Predicting Physical Properties of Nanofluids by Computational

Dec 27, 2016 - Under this investigation, a new database of thermal conductivity and ... This contribution aims at the creation of a unified expert sys...
0 downloads 0 Views 850KB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Article

Predicting Physical Properties of Nanofluids by Computational Modeling Natalia Sizochenko, Michael Syzochenko, Agnieszka Gajewicz, Jerzy Leszczynski, and Tomasz Puzyn J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.6b08850 • Publication Date (Web): 27 Dec 2016 Downloaded from http://pubs.acs.org on December 28, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Predicting Physical Properties of Nanofluids by Computational Modeling Natalia Sizochenko1,2, Michael Syzochenko2,3, Agnieszka Gajewicz1, Jerzy Leszczynski2, Tomasz Puzyn1* 1

Laboratory of Environmental Chemometrics, Faculty of Chemistry, University of Gdansk, Wita

Stwosza 63, 80-308, Gdansk, Poland 2

Interdisciplinary Center for Nanotoxicity, Department of Chemistry, Jackson State University,

1400 J. R. Lynch Street, P. O. Box 17910, 39217, Jackson, MS, USA 3

Department of Computer Science, Jackson State University, 1400 J. R. Lynch Street, P. O. Box

17910, 39217, Jackson, MS, USA

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 23

ABSTRACT: The focal point of the current contribution was to develop global quantitative structure-property relationship (QSPR) models for nanofluids. Two target properties: thermal conductivity and viscosity of nanofluids were thoroughly investigated. Under this investigation a new database of thermal conductivity and viscosity of nanofluids (more than 150 data points) was created. A hierarchical system of molecular representation reflecting features of nanoparticle’s structure at the different levels of organization was introduced. Also sizedependent, volume-dependent and intensive parameters were calculated. The model for thermal conductivity is characterized with determination coefficient R2 = 0.81 and root mean squared error RMSE = 0.055; model for viscosity: R2 = 0.79 and RMSE = 0.234. Developed models are in agreement with modern theories of nanofluids behavior. Size- and concentration-related behavior of target properties were discussed. Findings suggest that the increase of surface area ratio, interfacial layer thickness and decrease of nanoparticles size lead to thermal conductivity and viscosity increase. Thermal conductivity and viscosity increase with an increase of weighted fraction-dependent parameters. Up to date, reliable theoretical models were created only for a single type of nanoparticles. In the current article, developed models can simultaneously predict the thermal conductivity and viscosity in an effective way using both size and volume concentration of nanofluid.

2 ACS Paragon Plus Environment

Page 3 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction Nanofluids are engineered colloids of nanoparticles (NPs) dispersed homogenously in a base fluid.1 Transport, thermophysical and other related properties of nanofluids differ from properties of base fluids.2-4 For instance, addition of small fractions of NPs can cause significant enhancement of thermal conductivity. Many experimental and theoretical works reported equations for thermal conductivity and viscosity changes.2-3, 5-15,

16-18

Existing equations explain

the influence of nanofluid’s viscosity and thermal conductivity using different parameters, e.g. type of base fluids, NP volume fraction, particle size, particle shape, temperature, surface charge, pH value, particle nature, Brownian motion of NPs, effect of clustering, monolayer and dispersion techniques.1,

4, 19-20

Almost all available studies report size-dependent and/or

concentration-dependent thermal conductivity and viscosity focusing on one type of nanofluids. However, there is a lot of controversy over experimental findings and theoretical models. So far, no satisfactory mathematical models have been proposed to describe the thermal response of a nanofluid. One of the most promising and popular techniques to develop mathematical dependencies between target properties and features of the chemical structure is the Quantitative StructureProperty/Activity Relationships (QSPR or QSAR) approach.21 In the case of nano-sized compounds, this technique is often referred to as “nano-QSPR” or ”nano-QSAR”.22-23 A large number of nano-QSAR/QSPR models has been developed for inorganic nanoparticles since 2009.23 Despite the urgency of this task and promising results of the existing nano-QSAR/QSPR models, there is no universal methodology that provides a unified description for a variety of nanoscale structures. Moreover, direct description of dispersed-phase particles is possible only

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 23

using molecular dynamics models,8 but these methods are not suitable for fast predictions for screening purposes. There are several studies on the influence of the chemical structure on thermal conductivity and viscosity of nanofluids by using the QSPR approach.24-25 All of them employed simplified representation of the nanofluid structure. For example, in our previous contribution, nanofluids structures were represented by so-called “liquid drop” model.25 However, applicability of the model lacked (1) diversity (only 5 types of NPs were considered) and (2) volume fraction (NPs concentration lays within 0.01 - 0.55 %). Also only thermal conductivity was considered as a target property.25 The current study is aimed at a detailed analysis of the thermal conductivity and viscosity using hierarchical combination of newly developed descriptors, which reflect nanofluid structure at the different levels of organization. This contribution aims at the creation of a unified expert system, which could predict thermal conductivity and viscosity values for different types of nanofluids at different concentrations and sizes. In this paper, the values of thermal conductivities and viscosities for nanofluids, predicted with the proposed computational model are compared with the experimental data reported in the literature.

Experimental Data As it was discussed in previous contributions,25-26 physical properties of nanofluids may depend on the type of the base fluid as well as on the size of NPs. The original experimental data on thermal conductivity and viscosity were gathered from literature (see Supporting Information, Table S1) in the framework of NanoBRIDGES project (FP7-PEOPLE-2011-IRSES, grant agreement #295128). The gathered dataset was presented in current form as a part of

4 ACS Paragon Plus Environment

Page 5 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

NanoBRIDGES report. In the current investigation applicability domain of the model is restricted to spherical and near-spherical NPs at temperature range of 21-25 0C (the most popular temperature range among collected research articles). We have collected the largest available dataset, consisting of 100 nanofluids for which the experimentally measured thermal conductivity data were available (k) and 69 nanofluids for which the experimentally measured viscosity data were available (η). The new datasets are larger than the dataset, collected in our previous contribution (thermal conductivity, 23 data points).25 In our previous contribution the database consisted of nanofluids which contained metals (Cu, Fe) and metal oxides (Al2O3, TiO2, ZrO2) NPs. In current work, we report a highly diverse database which consists of experimental data for different types of nanofluids: metals (Au, Ag, Cu, Fe), metal oxides (Al2O3, CuO, Fe2O3, TiO2, ZrO2) and silica-derived (SiO2, SiC) nanofluids dispersed in water. The thermal conductivity ratio (k) was defined as the ratio of nanofluid’s thermal conductivity to thermal conductivity of pure fluid (water). The viscosity ratio (η) was defined as the ratio of nanofluid viscosity to viscosity of pure fluid (water). Additionally, the experimental data on the average NP size (SA) and the volume fraction of NPs ( ) for each nanofluid was extracted. Structures of NPs, target properties (the endpoints) and related experimental values are presented in SI (Tables S2 and S3).

Theoretical Methods Descriptors Several levels of the structure representation were applied. First, chemical diversity of the collected database was reflected by simple descriptors based on the periodic table (PTD): heat capacity (heat), mass density (ρ), molecular mass (M), indicator (1 for yes/0 for no) for metal

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 23

oxides of 4th or 5th period metal oxides (PTD-1), number of cations in metal oxides composition (PTD-2 for 4th period, PTD-3 for 5th period, PTD-5 for 3th period).31-32 Next, size-dependent descriptors derived from the “liquid drop” model (LDM) were utilized to describe the features of pristine nanoparticles.26 The scheme of LDM calculations was discussed in our previous contributions.25-26 It was demonstrated, that LDM encodes possible interactions between NPs using simplified geometrical representations of the particles at the nano-level. In LDM, the basic elements are densely packed and the density of the “liquid drop” is equal to the mass density (so-called Wigner-Seitz radius, eq. 1): 

 = 

 

(1)

where: M - molecular mass, ρ - mass density. LDM approximation assumes, that probable NP’s shape is a sphere. Based on this assumption, the number of elementary particles in a nanoaggregate is defined as follows (eq. 2): 

n =  





(2)

where: r0 – the mean radius of a NP. Using equation 2, nanoparticle’s surface area (F, eq. 3) and the surface-area-to-volume ratio (SV, eq. 4) could be calculated: F = 4 SV =





  

(3)

 ! #$%! %!

= #$%!

%! & '$%#!



(4)

Unique properties of nanofluids are affected by the layering of liquid molecules at the solid-particle surface (so-called, interfacial layer or nanolayer).27-28 Interfacial layer plays the role of a bridge between a solid NP and the base liquid. NPs have highly ordered structures, whereas base liquids are less organized. The interfacial layer is more ordered than the bulk 6 ACS Paragon Plus Environment

Page 7 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

liquid, but less ordered than the solid structure of NP. Thus, to reflect, specific adsorbed layer between surface of NPs and fluid, the interfacial nanolayer model (INM) was used (Figure 1).2729

Several models were developed to describe features of interfacial nanolayer, for example

thickness. Thickness of this layer (h, nm) is highly dependent on the temperature and can be described by the following equation (eq. 5):29 ( = 0.01 ∙ - . 273 ∙ 23.4

(5)

where 2 – size of single NP [nm].

Figure 1. Structure of nanoparticle in the base fluid. Molecules of the liquid (2) can form specific, highly ordered layer (3) near the nanoparticle surface (1).

Let us refer to the LDM- and INM-based descriptors as unweighted descriptors. Unfortunately, both types of descriptors are not able to express the effects related to changes of the volume fraction of NPs.3 When the value of the target property changes with the concentration of its components, a weighting procedure could be applied to describe concentration-dependent changes.30 Weighted descriptors are calculated as concentrationweighted sums using the descriptor value and volume fraction of NP (eq. 6): A6 =  ∙ 7,

(6)

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 23

where Aw - is the weighted descriptor,  – volume fraction of NPs, A – unweighted descriptor. It should be noted that in this case, each NP is defined by several points (several property values) corresponding to the different concentration values. Experimentally measured structure characteristics, such as size and concentration of NPs were also considered as descriptors. All calculated descriptors are provided in the electronic SI (Tables S4, S5, S6).

Model development and validation procedures For the purpose of validation, the initial dataset was split between training and test sets.21 Each training set covers ~ 80% of the initial dataset; related test set covers the remaining ~ 20% (see SI, Tables S7, S8). Nanofluids from the training set were structurally diverse enough to cover the whole descriptors space of the overall data set. Values of thermal conductivity were ranged in order to split datasets evenly. Each type of nanofluid (metal, metal oxide, etc.) was included in both training and test sets.21 Chemical diversities of nanofluids in the training and test sets were similar. After calculating an initial pool of descriptors, the descriptors with zero-variance and highly internally correlated (|r| > 0.80) ones were eliminated and the procedure of z-scaling (standardization) was performed. The models were developed using the Weka software package33 and its integration workflow plan for KNIME 2.11. For the purposes of QSPR modeling the M5P classifier was applied.33-34 M5P methodology combines a conventional decision tree technique with the possibility of including linear regression functions at the nodes. In the training set, each NP is specified by the fixed set of descriptors and has an associated target value (here: thermal

8 ACS Paragon Plus Environment

Page 9 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

conductivity or viscosity). Tree-based M5P models are constructed by the divide and conquer method. At the first step tree model is generated and specific tree nodes are specified. At the next step, multivariate linear model is constructed for the training samples at each node obtained at first step. The statistical fit of the QSPR models and their predictive ability were assessed using different internal and external validation strategies. The goodness-of-fit of the models was estimated by the quality metric R2, as well as external validation by metric R2test.35 Additional rm2 2

metrics for the training and test set data ( rm and ∆rm2) were estimated to determine the closeness between the values of predicted and corresponding observed.35 Moreover, for each model we calculated the root-mean-square error (RMSE) both for test and training sets. One of the OECD principles for model validation requires defining the Applicability Domain (AD) for the QSPR models.21 Reliable predictions are impossible only for chemicals structurally similar to the training compounds used to build the QSPR model. As in present work we used several types of chemically equal compounds, we suggest, that current models are applicable only for similar compounds in the same experimental conditions.

Results and Discussion Thermal conductivity model 2

The model is characterized by R2 = 0.81, RMSE = 0.055, rm = 0.76, ∆rm2 = 0.15 for the 2

training set; R2 = 0.70, rm = 0.67, ∆rm2 = 0.20 and RMSE = 0.077 for the test set. Dataset splitting, experimental and predicted values for the thermal conductivity model are summarized in SI (Table S7). Satisfactory values of the correlation coefficients and low values of the prediction errors confirm both accuracy of fitting, and predictive ability of the developed model. 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 23

Graphical representation of the correlation plot is presented in Figure 2.

Figure 2. Observed (experimental) vs. predicted values of thermal conductivity (k). Dots represent training set, squares – test set. It can be seen that some values are predicted with lower accuracy than others. Thermal conductivity ratios for nanofluids which contained larger NPs than usual, or high concentrations (higher than 1 %) were predicted with > 10 % error. For instance, ZrO2 (20 nm, 1 - 2 %, Table S2, nanofluids № 13 and № 14), TiO2 (50 nm, 1 - 2 %, Table S2, № 13 and № 14), Al2O3 (33 nm, 4 – 5 %, Table S2, № 49 and № 50), Cu (75 nm, Table S2, № 31) CuO (5 nm, 5 %, Table S2, № 55; and 23 nm, 10 %, Table S2, № 79). Nanofluid which contained CuO of 23.9 nm, 2.5 % (Table S2, № 63) was highly overestimated and thus excluded from final model. Mathematically, M5P model for thermal conductivity of nanofluids could be presented as 10 ACS Paragon Plus Environment

Page 11 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

a series of if-then-else rules (Figure 3).

Figure 3. M5P tree for thermal conductivity (k in %). As one can see (Figure 3), the developed model contains nine descriptors, from which the collection of rules included four pure descriptors: density (ρ), molecular weight (M), number of 5th period cations in chemical formula of metal oxide NP (PTD-3), presence/absence of metal

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 23

oxides of periods 4 and 5 (PTD-1); two weighted descriptors: number of NPs per volume fraction (nw) and heat capacity per volume fraction (heatw). Node equations included two unweighted descriptors – volume fraction of NPs ( ) and number of 4th period cations in chemical formula of metal oxide NP (PTD-2) and two weighted descriptors number of NPs per volume fraction (nw) and INM per volume fraction (hw). Relative importance values (RI), namely – the percentage of involvement of each descriptor – are presented in Table 1. Table 1. List of descriptors and their relative importance (RI, in %) in the thermal conductivity model Descriptor

Symbol

Descriptor type RI, %

volume fraction of NPs



experimental

3.5

molecular weight

M

PTD

14.0

presence/absence of metal oxides of periods 4 and 5

PTD-1

PTD

15.0

number of 4th period cations in chemical formula of metal oxide NP

PTD-2

PTD

4.0

number of 5th period cations in chemical formula of metal oxide NP

PTD-3

PTD

4.0

density

ρ

PTD

30.0

number of NPs per volume fraction

nw

weighted LDM

8.5

heatw

weighted PTD

5.0

hw

weighted INM

15

heat capacity per volume fraction interfacial layer thickness per volume fraction

Technical model node for KNIME is presented in SI. Installation process is presented on official website (https://www.knime.org/). Potential users can install this node and use it to make own prediction. Alternatively, if-then-else rules could be applied directly. For instance, for NPs, which mass density (ρ) is less than 4580, the first equation from if-then-else model should be applied (eq. 7): 12 ACS Paragon Plus Environment

Page 13 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

k = 7.6103 · PTD-1 + 2.8122 · hw + 0.9347 (7) In case if mass density (ρ) is higher than 4580, and molecular mass is lower than 59.75, the second equation from if-then-else model should be applied (eq. 8): k = -1.4912 ·  + 15.0476

(8)

This procedure should be repeated until suitable rules are found.

Viscosity model

The developed model for viscosity is characterized with quite good statistical parameters: 2

2

R2 = 0.79, rm = 0.73, ∆rm2 = 0.088 and RMSE = 0.234 for the training set; R2 = 0.78, rm = 0.72, ∆rm2 = 0.14 and RMSE = 0.244 for the test set. Dataset splitting, experimental and predicted values for the viscosity model are summarized in SI (Table S8). The plot of experimental versus predicted values is shown on Figure 4.

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 23

Figure 4. Observed (experimental) vs. predicted values of viscosity (η). Dots represent test set, squares – training set. Viscosity of nanofluids which contained Al2O3 at high concentrations (higher than 5 %) were usually overestimated or underestimated with > 23 % error. For instance, nanofluids (Table S3) № 5 (5 %), № 6 (6 %), № 12 (9 %), № 15 (8.5 %) and № 16 (12.2 %). Next, in the case of SiC (Table S3, № 66-69) predictions were similar. This issue could be related to non-metallic behavior of selected nanofluids (All other NPs were metals or metal oxides). Similar situation was observed for nanofluids with CuO at 33 nm (Table S3, № 53-59). At the same time, predictions for other nanofluids with CuO (Table S3, № 60-65) were characterized with high variation (Pearson correlation coefficient 0.93).

14 ACS Paragon Plus Environment

Page 15 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Mathematically, M5P model for viscosity of nanofluids could be presented as a series of if-then-else rules (Figure 5).

Figure 5. M5P tree for viscosity. The developed model for viscosity contains seven parameters. Collection of if-then-else rules included two pure descriptors: number of 3rd period cations in chemical formula of NP (PTD-4) and interfacial nanolayer thickness (h); three weighted descriptors (number of NPs per volume fraction, nmix; heat capacity per volume fraction, heatw; surface ratio, Fw). Node equations included pure descriptors of presence/absence of TiO2, (PTD-5) and three weighted descriptors: density per volume fraction (ρw), number of NPs per volume fraction (nw) and heat 15 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 23

capacity per volume fraction (heatw). Relative variable importance of descriptors is presented in Table 2. Technical information about KNIME node are presented in SI. Table 2. List of descriptors and their relative importance (RI in %) in the viscosity model Descriptor

Symbol Descriptor type RI, %

number of 3rd period cations in chemical formula of NP

PTD-4

PTD

29.0

presence/absence of TiO2

PTD-5

PTD

8.0

heat capacity per volume fraction

heatw

weighted PTD

8.0

density per volume fraction

ρ

weighted PTD

8.0

surface area ratio per volume fraction

Fw

weighted LDM

19.0

number of NPs per volume fraction

nw

weighted LDM

29.0

interfacial layer thickness

hw

INM

15.0

Interpretation of developed models

The models developed within this study represent physical properties of nanofluids as functions of experimentally measured and calculated descriptors. As it is shown in Tables 1-2, there are several ways of combining important descriptors into groups. Let us refer to these groups with fundamental physical concepts. One important physical parameter, which plays a major role in nanofluid thermal properties and stability, is the size of NP.2-5 In most cases, thermal conductivity and viscosity could linearly or non-linearly increase as particle size decreases.4 In the developed models this outcome is supported by the size-dependent LDM- and INM-based descriptors, namely surface area ratio, interfacial layer thickness and number of NPs. As is was previously mentioned “liquid drop” model-derived and interfacial nanolayer model derived descriptors are size-dependent. As developed M5P model is non-linear, both effects of potential linear and non-linear decrease and

16 ACS Paragon Plus Environment

Page 17 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

increase are reflected. Cumulative impact of these descriptors is 23.5% for the thermal conductivity model and 63% for the viscosity model. The other important empirical finding is NPs fraction-property response of nanofluids.1-10 Thermal conductivity and viscosity increase with volume concentration. As weighted descriptors are fraction-dependent, these fundamental relationships are clearly reflected in the developed models. Cumulative relative influence for the fraction-related descriptors is 32% for the thermal conductivity model and 64% for the viscosity model. On the other hand, with the decrease in volume fraction, effective volume of the cluster increases, hence the enhancement of thermal conductivity and viscosity. Clustering of NPs can create paths of lower thermal resistance and enhance thermal conductivity significantly. If one has an aggregated cluster of closely packed spherical particles, ~ 25% of the volume is filled with liquid in the space between particles, which increases the effective volume of the highly conductive region by ~ 30%. The last group of properties – intensive properties of NPs (here: PTD) contribute 72% for the thermal conductivity model and 53% for the viscosity model. An intensive property is a bulk property, meaning that it is a physical property of a system that does not depend on the system size or the amount of material in the system (here: periodic-table based). Comparing our results with existing models (here: models for nanofluids with water as a base fluid) we have found that every research group compared their data with existing theoretical models only for single type of NPs (using different volume concentration).1,3,4 For instance, Duangthongsuk et al.,6 Lee et al.,13 Li et al.,15 and Paul et al.19 measured the target properties and

built relationships of them with volume fraction or compared their own results with results of other authors. Kang et al.,11, Lee et al.,14 and Xie et al.,17 measured target properties and compared their results with the Hamilton-Crosser model.9 Kang et al.,11 found that the Hamilton-

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 23

Crosser model predicts properties of Ag nanofluids with lower values in comparison to the experimental data; this model works well for SiO2 nanofluids. However, Kang et al.,11 does not compare Ag and SiO2 nanofluids simultaneously. Lee et al.,14 demonstrated that his experimental data is in agreement with the Hamilton-Crosser model for Al2O3 nanofluids, but gave inadequate results for CuO nanofluids. Beck et al. 5 used a simple equation (k = 4.4134  ) to measure thermal conductivity for

Al2O3 nanofluids. R2 for this equation was between 0.90 and 0.93. Unfortunately, this equation was not tested for other types of nanofluids and is not applicable to NPs of different sizes. Therefore, our models are the only models, which can simultaneously predict the thermal conductivity and viscosity in an effective way using both size and volume concentration of nanofluid.

Conclusions In the present study we have successfully applied the computational modeling to build theoretical models of thermal conductivity and viscosity of nanofluids. We introduced a new hierarchical system of structure representation which allows reflecting features of nanofluids’ at the different levels of organization: from intensive properties (as period of metal, density, etc.) to nanoparticles size, concentration and thickness of interfacial layer. Existing models of thermal conductivity and viscosity of nanofluids can predict only sizeor concentration-dependent behavior. On the contrary, models developed in the framework of hierarchical approach, make an allowance for predictions of thermal conductivity and viscosity for nanofluids that contain NPs of different sizes, concentrations and types (metal, silica- or metal oxide). KNIME nodes for developed models are available for public use.

18 ACS Paragon Plus Environment

Page 19 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

We have developed statistically significant and externally predictive quantitative structure-property relationships models. Models were obtained using the largest available dataset of structurally diverse nanofluids, including a variety of concentration and size classes. From this study it is clear that the best models included a combination of different types of descriptors. Developed models are in agreement with the fact, that viscosity and thermal conductivity depend on many parameters such as particle volume fraction, particle size, particle size distribution and thickness of interfacial layer. From the performed calculations we found the descriptors that are highly correlated to the observed viscosity and thermal conductivity. These descriptors are mainly related to the volume fraction and size of the nanomaterials. It was found that size-dependent “liquid drop” model and the interfacial nanolayer model derived descriptors (surface area ratio, interfacial layer thickness, number of NPs) have positive influence on thermal conductivity. Weighted fraction-dependent descriptors reflected an increase in thermal conductivity and viscosity with increase of volume concentration. Applicability of the models is limited to spherical nanoparticles. We believe, that the obtained models reflect substantial factors concerning the physical properties of nanofluids. The developed models could help in the future with the development of suitable candidate materials.

AUTHOR INFORMATION Corresponding Author *Phone: + 48 (058) 523-5248 Dr. Tomasz Puzyn, e-mail: [email protected] Author Contributions N.S. and M.S. contributed equally. N.S. and M.S. collected data, carried out calculations, processed the data, developed code and analyzed results. N.S., M.S. and A.G. drafted and revised 19 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 23

the manuscript. J.L. and T.P. supervised the research, coordinated the study and revised the manuscript. All authors have read and have given approval to the final version of the manuscript. Notes The authors declare no competing financial interests. ACKNOWLEDGMENTS The authors thank the support from the NSF CREST Interdisciplinary Nanotoxicity Center: grants EPSCoR 362492-190200-01\NSFEPS-0903787 and N00014-13-1-0501. This research has received funding from the European Union Seventh Framework Program, NanoPUZZLES project (FP7/2007-2013, grant agreement #309837) and the European Commission through the Marie Curie IRSES program, NanoBRIDGES project (FP7-PEOPLE-2011-IRSES, grant agreement #295128). SUPPORTING INFORMATION Collected experimental data, values of calculated descriptors and KNIME code are provided in the electronic Supporting Information.

ABBREVIATIONS NP nanoparticle; LDM “liquid drop” model; QSAR Quantitative Structure-Activity Relationships; QSPR Quantitative Structure-Property Relationships; RMSE the root-meansquare error.

REFERENCES 1.

Saidur, R.; Leong, K. Y.; Mohammad, H. A. A Review on Applications and Challenges 20 ACS Paragon Plus Environment

Page 21 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

of Nanofluids. Renew. Sust. Energ. Rev. 2011, 15, 1646-1668. 2. Chandrasekar, M.; Suresh, S.; Senthilkumar, T. Mechanisms Proposed through Experimental Investigations on Thermophysical Properties and Forced Convective Heat Transfer Characteristics of Various Nanofluids – a Review. Renew. Sust. Energ. Rev. 2012, 16, 39173938. 3. Kleinstreuer, C.; Feng, Y. Experimental and Theoretical Studies of Nanofluid Thermal Conductivity Enhancement: A Review. Nanoscale Res. Lett. 2011, 6, 1-13. 4. Özerinç, S.; Kakaç, S.; Yazıcıoğlu, A. Enhanced Thermal Conductivity of Nanofluids: A State-of-the-Art Review. Microfluid Nanofluid 2010, 8, 145-170. 5. Beck, M.; Yuan, Y.; Warrier, P.; Teja, A. The Effect of Particle Size on the Thermal Conductivity of Alumina Nanofluids. J. Nanopart. Res. 2009, 11, 1129-1136. 6. Duangthongsuk, W.; Wongwises, S. Measurement of Temperature-Dependent Thermal Conductivity and Viscosity of TiO2-Water Nanofluids. Exp. Therm. Fluid Sci. 2009, 33, 706714. 7. Eastman, J. A.; Choi, U. S.; Li, S.; Thompson, L. J.; Lee, S. Enhanced Thermal Conductivity through the Development of Nanofluids. Mater. Res. Soc. Symp. Proc. 1996, 457, null-null. 8. Ghosh, M. M.; Roy, S.; Pabi, S. K.; Ghosh, S. A Molecular Dynamics-Stochastic Model for Thermal Conductivity of Nanofluids and Its Experimental Validation. J. Nanosci. Nanotechno. 2011, 11, 2196-2207. 9. Hamilton, R. L.; Crosser, O. K. Thermal Conductivity of Heterogeneous TwoComponent Systems. Ind. Eng. Chem. Fund. 1962, 1, 187-191. 10. Hwang, Y.; Park, H. S.; Lee, J. K.; Jung, W. H. Thermal Conductivity and Lubrication Characteristics of Nanofluids. Curr. Appl. Phys. 2006, 6, Supplement 1, e67-e71. 11. Kang, H. U.; Kim, S. H.; Oh, J. M. Estimation of Thermal Conductivity of Nanofluid Using Experimental Effective Particle Volume. Exp. Heat Transfer 2006, 19, 181-191. 12. Koo, J.; Kleinstreuer, C. A New Thermal Conductivity Model for Nanofluids. J. Nanopart. Res. 2004, 6, 577-588. 13. Lee, J.-H.; Hwang, K. S.; Jang, S. P.; Lee, B. H.; Kim, J. H.; Choi, S. U. S.; Choi, C. J. Effective Viscosities and Thermal Conductivities of Aqueous Nanofluids Containing Low Volume Concentrations of Al2O3 Nanoparticles. Int. J. Heat Mass Tran. 2008, 51, 2651-2656. 14. Lee, S.; Choi, S. U. S.; Li, S.; Eastman, J. A. Measuring Thermal Conductivity of Fluids Containing Oxide Nanoparticles. J Heat Tran. 1999, 121, 280-289. 15. Li, Q.; Xuan, Y.; Wang, J. Experimental Investigations on Transport Properties of Magnetic Fluids. Exp. Therm. Fluid Sci. 2005, 30, 109-116. 16. Wei, W. A Comprehensive Model for the Enhanced Thermal Conductivity of Nanofluids. J. Adv. Res. Phys. 2013, 3. 17. Xie, H.; Wang, J.; Xi, T.; Liu, Y. Thermal Conductivity of Suspensions Containing Nanosized Sic Particles. Int. J. Thermophys. 2002, 23, 571-580. 18. Zhu, H.; Zhang, C.; Liu, S.; Tang, Y.; Yin, Y. Effects of Nanoparticle Clustering and Alignment on Thermal Conductivities of Fe3O4 Aqueous Nanofluids. Appl. Phys. Lett. 2006, 89, 023123. 19. Paul, G.; Chopkar, M.; Manna, I.; Das, P. K. Techniques for Measuring the Thermal Conductivity of Nanofluids: A Review. Renew. Sust. Energ. Rev. 2010, 14, 1913-1924. 20. Zhong, S. J. n. A. A. S.-K. n. A. W.-H. Enhancement of Fluid Thermal Conductivity by the Addition of Single and Hybrid Nano-Additives. Thermochim. Acta 2007. 21 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 23

21. OECD, Guidance Document on the Validation of (Quantitative) Structure-Activity Relationship [(Q)Sar] Models; OECD Publishing. 22. Puzyn, T.; Leszczynska, D.; Leszczynski, J. Toward the Development of "Nano-Qsars": Advances and Challenges. Small 2009, 5, 2494-2509. 23. Sizochenko, N.; Leszczynski, J. Review of Current and Emerging Approaches for Quantitative Nanostructure-Activity Relationship Modeling - the Case of Inorganic Nanoparticles. Journal of Nanotoxicology and Nanomedicine 2016, 1, 1-16. 24. Toropov, A. A.; Leszczynska, D.; Leszczynski, J. Predicting Thermal Conductivity of Nanomaterials by Correlation Weighting Technological Attributes Codes. Mater. Lett. 2007, 61, 4777-4780. 25. Sizochenko, N., Jagiello, K.; Leszczynski, J.; Puzyn, T. How “Liquid Drop” Approach Could Be Applied for Quantitative Structure-Property Modeling of Nanofluids. J. Phys. Chem. C. 2015, 119, 25542–25547. 26. Sizochenko, N.; Rasulev, B.; Gajewicz, A.; Kuz'min, V.; Puzyn, T.; Leszczynski, J. From Basic Physics to Mechanisms of Toxicity: The "Liquid Drop" Approach Applied to Develop Predictive Classification Models for Toxicity of Metal Oxide Nanoparticles. Nanoscale 2014, 6, 13986-13993. 27. Yu, W.; Choi, S. U. S. The Role of Interfacial Layers in the Enhanced Thermal Conductivity of Nanofluids: A Renovated Maxwell Model. J. Nanopart. Res. 2003, 5, 167-171. 28. Leong, K. C.; Yang, C.; Murshed, S. M. S. A model for the thermal conductivity of nanofluids – the effect of interfacial layer. J. Nanopart. Res. 2006, 8, 245-254. 29. Sitprasert, C.; Dechaumphai, P.; Juntasaro, V. A thermal conductivity model for nanofluids including effect of the temperature-dependent interfacial layer. J. Nanopart. Res. 2009, 11, 1465-1476. 30. Oprisiu, I.; Novoratskyi, S.; Tetko, I. V. Modeling of non-additive mixture properties using the Online CHEmical database and Modeling environment (OCHEM). J. Chemoinform. 2013, 5, 4-10. 31. Kar, S.; Gajewicz, A.; Puzyn, T.; Roy, K.; Leszczynski, J. Periodic Table-Based Descriptors to Encode Cytotoxicity Profile of Metal Oxide Nanoparticles: A Mechanistic QSTR Approach. Ecotoxicol. Environ. Saf. 2014, 107c, 162-169. 32. Sizochenko, N.; Gajewicz, A.; Leszczynski, J.; Puzyn, T. Causation or only correlation? Application of causal inference graphs for evaluating causality in nano-QSAR models. Nanoscale 2016, 8, 7203-7208. 33. Witten, I. H.; Frank, E.; Trigg, L.; Hall, M.; Holmes, G.; Sally J. C. Weka: Practical Machine Learning Tools and Techniques with Java Implementations; Hamilton, New Zealand: University of Waikato, Department of Computer Science., 1999. 34. Quinlan, J. R. Learning with continuous classes. In Processings AI’92, 343-359, Singapore: World Scientific, 1992. 35. Roy, K.; Chakarborty, P.; Mirta, M.; Ojha, P. K.; Kar, S.; Das, R. N. Some Case Studies on Application of ‘‘rm2’’ Metrics for Judging Quality of Quantitative Structure–Activity Relationship Predictions: Emphasis on Scaling of Response Data. J. Comp Chem. 2012, 34, 1071-1082.

22 ACS Paragon Plus Environment

Page 23 of 23

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table of Contents

23 ACS Paragon Plus Environment