Numerical Investigation of the Microscopic Heat Current inside a

3 hours ago - The curves fit to data of the probability distribution of the microscopic heat current are established, and the parameters including the...
0 downloads 7 Views 425KB Size
Subscriber access provided by Grand Valley State | University

Numerical Investigation of the Microscopic Heat Current inside a Nanofluid System based on Molecular Dynamics Simulation and Wavelet Analysis Tao Jia, and Di Gao Anal. Chem., Just Accepted Manuscript • DOI: 10.1021/acs.analchem.7b05350 • Publication Date (Web): 19 Mar 2018 Downloaded from http://pubs.acs.org on March 19, 2018

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 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 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.

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

Analytical Chemistry

Numerical Investigation of the Microscopic Heat Current inside a Nanofluid System based on Molecular Dynamics Simulation and Wavelet Analysis Tao Jia1

Di Gao2

1. Department of Thermal Engineering, Taiyuan University of Technology, Taiyuan, 030024, China 2. Department of Statistics, North Dakota State University, Fargo, 58102, United States

Abstract Molecular dynamics simulation is employed to investigate the microscopic heat current inside an argon-copper nanofluid. Wavelet analysis of the microscopic heat current inside the nanofluid system is conducted. The signal of the microscopic heat current is decomposed into two parts: one is the approximation part; the other is the detail part. The approximation part is associated with the low frequency part of the signal, and the detail part is associated with the high frequency part of the signal. Both the probability distributions of the high frequency and the low frequency parts of the signals demonstrate Gaussian-like characteristics. The curves fit to data of the probability distribution of the microscopic heat current are established, and the parameters including the mean value and the standard deviation in the mathematical formulas of the curves show dramatic changes for the cases before and after adding copper nanoparticles into the argon base fluid.

1

Corresponding author: [email protected]

ACS Paragon Plus Environment

Analytical 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

Keywords: Nanofluid; Molecular Dynamics Simulation; Wavelet analysis; microscopic heat current

Introduction Nanofluid is a homogenous mixture of nanoparticles and conventional fluid which is often termed base fluid, the nanoparticles are uniformly dispersed with low concentration in the base fluid, and the diameter of the nanoparticle is usually less than 100nm. Heat transfer is a crucial issue in the thermal management in industry. A traditional method to augment heat transfer was to enlarge the surface of heat transfer. However, once the equipment is built and installed on site, the surface is almost impossible to be changed. Nanofluid provides an effective and easy-toimplement solution to enhance the heat transfer without increasing the heat transfer surface. The thermal performance of nanofluid composed of carbon nanotubes and water was experimentally investigated, and it was found that the thermal conductivity of the nanofluid increased as the length of the nanotube decreased [1]. ZnO nanoparticle and propylene glycol were mixed to construct a nanofluid, and the technologies of ultraviolet–visible spectroscopy and light scattering were employed to study the nanoparticle size distribution and the stability of the suspension of the nanoparticle in the base fluid [2]. The effects of temperature and nanoparticle concentration on the viscosity of naofluids were investigated, and the result showed that the viscosity decreased with temperature and increased with the nanoparticle concentration [3]. In the base fluid composed of 60 mas% propylene glycol and 40 mas% water, different nanoparticles made of Al2O3, ZnO, TiO2, and SiO2 were injected, the nanofluids exhibited lower surface tension, compared to the surface tension of the base fluid [4]. The thermal performance of water-silver nanofluids were experimentally studied. The nanofluid was formed by the twostep method, and the experiment results demonstrated that the thermal performance increased

ACS Paragon Plus Environment

Page 2 of 23

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

Analytical Chemistry

about 28% by using the nanofluid [5]. Ultrasound-based technology was employed to detect the concentration characterization of the water-copper nanofluid. The ultrasonic velocity in the nanofluid was greater than that in the base fluid [6]. Experimental investigation of the thermal performance of the nanofluid flowing in an equilateral triangular duct was carried out. Nanoparticles made of Al2O3 and SiO2 were added to distilled water, and then were dispersed with ultrasonic vibration method. Finally a significant increase of heat transfer was found under the condition of the employment of the nanofluids [7]. The effect of using water-silver nanofluid on the thermal performance of an inclined heat pipe was experimentally investigated. Four different inclination angles including 0°, 30°, 60° and 90° were considered. As the concentration of the silver nanoparticle increased the thermal resistance of the heat pipe decreased. The thermal conductivity of the heat pipe reached maximum at the angle of 60° [8]. Al2O3 nanoparticles were distributed inside ethylene glyol to form a nanofluid, and an improved transient-hot-wire method was used to measure the thermal conductivity of the nanofluid. The influences of the nanoparticle volume fraction and size distribution on the thermal conductivity were investigated. The experiment showed that the influence of particle size was over that of particle volume fraction on the thermal conductivity under the condition of volume fraction less than 0.25% [9]. The character of TiO2-water nanofluid was experimentally investigated under the conditions of laminar and mixed flows. Brownian motion and thermophoresis were assumed to be responsible for the free movement of the nanoparticle in the base fluid [10]. Contrast-enhanced video microscopy technology was used to find the growth of fingering patterns in dewetting nanofluids [11]. Both temperature-sensitive and non-temperature-sensitive magnetic nanofluids were discussed, and the particle rotation and the influence of the torque on the particle by an externally magnetic field were thought to be the explanation of the effective viscosity of the

ACS Paragon Plus Environment

Analytical 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

nanofluid [12]. Experiments showed that the employment of highly thermophilic particles resulted in an active control of the heat transferred by smart nanofluids [13]. It was found that there was a dramatic increase in heat transfer by using nanfluid under the condition of turbulent flow in a numerical investigation based on Euler and Lagrangian methods. Two different kinds of nanofluids including Cu and Al2O3 nanoparticles were investigated [14]. The shadowgraph images were used to discover the pattern of power law scaling in the early stage of solutal convection in a nanofluid under conditions far from equilibrium [15]. Nanofluid has been widely used in boiling heat transfer which is a high efficient heat transfer technique [16, 17]. The investigation of the thermal performance of nanofluids in convective boiling flows is experimentally conducted. It was found that adding alumina and copper oxide nanoparticles and multi-walled carbon nanotubes to water did enhance the whole thermal performance of the boiling flow system [16]. Experimental investigation of the heat transfer coefficient of Titana nano-fluids was carried out in the boiling heat transfer process. It was shown that the heat transfer coefficient of the nanofluid system is higher than that of the base fluid, and the nucleation site density was increased after adding the nanoparticles [18]. It was found that the heat transfer coefficient of a copper-made heat sink with rectangular microchannel was increased in the regime of laminar flow after using the nanofluid composed of silver nanoparticles and deionized water as the coolant [19]. Green synthesis method is employed to produce silver nanoparticles. The produced nanoparticles were shown to have good stability when they are dispersed in deionized water [20]. The influences of stirring, sonication, and surface active agents on the stability of cupric oxide nanoparticles dipersed in water, ethylene glycol, and mixture of water and ethylene glycol were experimentally investigated. The experimental results showed that the ethylene glycol worked best to disperse the cupric oxide nanoparticles [21]. The

ACS Paragon Plus Environment

Page 4 of 23

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

Analytical Chemistry

thermal performances of three kinds of coolants: gallium, nanofluid composed of cupric oxide nanoparticles and water, and water in experiments to cool CPU were investigated. It was found that gallium had the best performance [22]. Carbon nano tube aqueous nanofluid was used as a coolant to flow inside a rectangular microchannel of a copper-made sink, and it was found that as the mass concentration of the nanoparticle increased, the heat transfer coefficient increased significantly [23]. An experimental study of the thermal performance of the nanofluid composed of alumina nanoparticles and water under the process of flow boiling heat transfer was conducted. The experimental results were that the heat transfer coefficient of the thermal system increased during short time period between 0 and 60 minutes, and it decreased in the extended time period between 60 and 1000 minutes [24]. Molecular dynamics simulation serves as an effective method to explore the microscopic behavior. The information of each individual microscopic particle including its velocity and momentum can be calculated based on Verlet algorithm [25-27]. The microscopic heat current is defined based on the information of all particles, and it is directly linked to the thermal conductivity of the nanofluid. So the investigation of the stochastic characteristics of the microscopic heat current can help us understand the microscopic thermal fluctuation inside the nanofluid system [25-26]. Wavelets-based technology [29-32] has found its application in a wide range of engineering fields, and it provides us a powerful tool to qualitatively represent a variety of signals. The basic idea of wavelet is to decompose the original signal into approximate and detail parts. The approximation part reflects the general trend of the original signal, and it is associated with high scale and low frequency. The detail part reflects local difference in the original signal, and it is associated with low scale and high frequency. The characteristics of the signal of the

ACS Paragon Plus Environment

Analytical 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

microscopic heat current inside a nanofluid system is a direct representation of the microscopic thermal transport inside the nanofluid system [25-27]. In this paper, the wavelet analysis of the microscopic heat current signal is conducted with the aim to find some certain patterns as the nanoparticle volume fraction changes. SIMULATION The numbers of the argon atoms are 4000, 3859, 3799, 3775, 3751, 3748, and the numbers of copper atoms are 137, 381, 456, 528, 601 for the cases of copper nanoparticle volume fraction of 0%, 2%, 2.5%, 3%, 3.5%, 4% respectively. On each side of the simulation box, periodical boundary condition is applied to ensure the configuration that the nanoparticles are uniformly distributed in the base fluid. At the beginning of the simulation, the copper nanoparticle is arranged at the center of the simulation box. Figure 1 illustrates that copper nanoparticles are distributed in the argon base fluid. The large particles with black color represent the copper nanoparticle and the smaller particles with blue color represent the argon atoms. Different cases of nanoparticle volume fractions are investigated. As to the cases of nanoparticle volume fraction of 0%, 2%, 2.5%, 3%, 3.5%, 4%, the numbers of argon atoms are 4000, 3859, 3799, 3775, 3751, 3748, and the numbers of copper atoms are 137, 381, 456, 528, 601 respectively. To describe the interaction between the atoms including the interaction between argon atoms, the interaction between copper atoms, and the interaction between argon and copper atoms, Lennard-Jones potential function [27, 28] is employed. The potential well depth is 1.67 × 10−21 J and the characteristic length is 0.3405 nm for the interactions between argon atoms. The potential well depth is 65.63 × 10−21 J and the characteristic length is 0.233 nm for the interaction between copper atoms. The potential well depth is 10.415 × 10−21 J and the characteristic length is 0.287 nm for the interaction between argon atom and copper atom based on Berthlot mixing rule [27, 28]. To get

ACS Paragon Plus Environment

Page 6 of 23

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

Analytical Chemistry

the thermal conductivity of the nanofluid, Green-Kubo relation [27, 28] is employed. GreenKubo relation is widely used to calculate the macroscopic thermal transport coefficient from the microscopic quantities. The underlying mechanism is the fluctuation-dissipation theorem [27, 28] which is first discovered by Einstein to mathematically relate the macroscopic transport coefficient of a material to the microscopic properties inside the material. To describe the transport inside a system, the phenomenological relationship is employed [28] : Flux = −TransportCoefficient × Gradient

(1)

Based on the above, we have the following

Ex& = − D

∂E ∂x

(2)

dx is the derivative of the dt position with respect to time, that is the velocity at the position x , D is the diffusion coefficient, Ex& is the flux, which is a measure of the magnitude and direction of the flow of the energy . where E = E ( x, t ) is the energy at the position x and at time t . x& =

The conservation of energy leads to the following:

∂E ∂ ( Ex& ) + =0 ∂t ∂x

(3)

Combining the above two equations, we get the following diffusion equation:

∂E ∂2 E =D 2 ∂t ∂x

(4)

If originally, the energy E0 is located at the position x = 0 at the time t = 0 , then the solution of Eq (4) is the following:

E ( x, t ) =

E0 2 π Dt

e

− x2 4 Dt

(5)

To obtain the link between the microscopic property and the macroscopic property, the mean square displacement of the energy is used as the following:

ACS Paragon Plus Environment

Analytical 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

〈 [ E ( x, t ) − E ( x, 0)]2 〉 =

1 L [ E ( x, t ) −E ( x, 0)]2 dx ∫ 0 L

Page 8 of 23

(6)

where L is the length of the one-dimensional domain. The combination of Eq (5) and Eq (6) results in the following: 〈 [ E ( x, t ) − E ( x, 0)]2 〉 = 2 Dt

(7)

As t → ∞ , we have: 〈[ E ( x, t ) − E ( x, 0)]2 〉 =D lim t →∞ 2t The above is Einstein relation [27, 28] which calculates the macroscopic transport coefficient based the time average of the microscopic properties. Green-Kubo relation can be obtained from Einstein relation; the two are proved to be equivalent [27, 28]. The calculation of the thermal conductivity through Green-Kubo relation is the following [27, 28]: k=

1 3VkBT 2





J (0) J (t ) dt

(8)

0

where k is the thermal conductivity, kB is Boltzmann constant ( 1.38 × 10−23 J/K ), T and V are the temperature and the volume of the system respectively, and J is the vector the microscopic heat current. To make the argon be in a liquid state, the nanofluid system temperature is kept as 86 K during the simulation. The microscopic heat current is determined as the following:

J=

d N ∑ ri Ei dt i =1

(9)

where ri is atom position vector, Ei is the site energy of atom i , and N is the total atom number: 1 1 2 Ei = mi vi + ∑ uij 2 2 j

ACS Paragon Plus Environment

(10)

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

Analytical Chemistry

where mi is the atom mass, vi is the atom velocity, and uij is the potential energy between two atoms. The microscopic heat current in the nanofluid system, is determined as the following [26]: kl Nk β Nk β  1 1 β β N k Nl  ∂φ J = ∑∑ mik (vik ) 2 v ik − ∑∑∑∑ rijkl ijkl − φijkl I  ⋅ v ik − ∑ hk ∑ v ik 2 k =α l =α i =1 j =1  ∂rij k =α i =1 2 k =α i =1 

(11)

there are two components including argon and copper in the nanofluid system, α and β denote copper and argon atoms respectively, v ki is the velocity of atom i of species k ( k = α or β ), Nα is the number of atoms of specie α , N β is the number of atoms of specie β , hk is the mean partial enthalpy, and I is a unit tensor. In order to simplify the analysis of the microscopic heat current signal, we have conducted the nondimensionalization of the signal. It is rooted on four fundamental quantities including argon atom mass, the characteristic length of argon atom, and the potential well depth of argon atom.

J * = J ma / ε aa / (ε aaσ aa2 )

(12)

where ma is the argon atom mass, σ aa is the characteristic length of the argon atom, and ε aa is the potential well depth of argon atom. The time-evolution of the signal of the nondimensionalized microscopic heat current is shown in Figure 2. The approximation and detail parts of the signals of the microscopic heat currents inside the nanofluids with different nanoparticle volume fractions are plotted in Figure 3 and Figure 4 respectively. Verlet algorithm [27] is chosen to conduct the simulation. It is rooted on Taylor series which decomposes a function into an infinite sum of the function's derivatives at a reference

ACS Paragon Plus Environment

Analytical 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

point. As to the position of one atom, r , the Taylor series r (t + ∆t ) and r (t − ∆t ) are expressed as follows: dr (t ) 1 d 2 r (t ) 2 1 d 3 r (t ) 3 r (t + ∆t ) = r (t ) + ∆t + ∆t + ∆t + ο (∆t 4 ) 2 3 dt 2 dt 3! dt

(13)

dr (t ) 1 d 2 r (t ) 2 1 d 3 r (t ) 3 ∆t + ∆t − ∆t + ο (∆t 4 ) dt 2 dt 2 3! dt 3

(14)

r (t − ∆t ) = r (t ) −

where t is time and ∆t is the time-step in the simulation. The sum of Eq (10) and Eq (11) results in the following:

r (t + ∆t ) = 2r (t ) − r (t − ∆t ) +

d 2 r (t ) 2 ∆t + ο (∆t 4 ) dt 2

(15)

Based on the first-order central difference estimator, we have the velocity of the atom at time t as the following:

v(t ) =

r (t + ∆t ) − r (t − ∆t ) 2∆t

(16)

and at the next time-step the velocity of the atom is determined as the following : d 2r v(t + ∆t ) = v(t ) + 2 ∆t dt

(17)

Haar wavelet is employed here to analyze the microscopic heat current. The mother wavelet function of Haar wavelet is mathematically described as the following:

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

Analytical Chemistry

 1 0 ≤ x < 0.5  ϕ ( x) = −1 0.5 ≤ x < 1  0 otherwise 

(18)

and the scaling function of Haar wavelet is mathematically described as the following:

1 0 ≤ x < 1 0 otherwise

φ ( x) = 

(19)

Haar wavelet has been extensively used in the approximation of signal. The basic idea of the Haar-wavelet-based analysis is to decompose the original signal into two parts: one is to represent the trend of the signal which manifestates the global character of the original signal; the other is to represent the detail in local the domain which serves as the local signature of the original signal. We give a simple example here. As to a signal composed of two numbers, {10, 15}, the number 12.5=(10+15)/2 represents the trend of the signal, and the number -2.5 = (1015)/2 represents local detail of the signal..

RESULTS AND DISCUSSION The curve fitting of the microscopic heat current signal is conducted based on least-square method. It is found that both the probability distribution of the approximation and the probability distribution of detail parts of the microscopic heat current signal can be described as the following:

y = I × exp(( x − c) 2 /(2 × σ 2 ))

ACS Paragon Plus Environment

(20)

Analytical 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

where the parameter I represents the amplitude of the signal, the parameter c represents the position of the line x = c with respect to which the curve is symmetric, the parameter σ is the standard deviation, the larger is σ the wider range does the probability distribution covers. Table 1 and Table 2 give the data of the above parameters under the different conditions of nanoparticle volume fractions. The R squared analysis for the data in Table 1 and Table 2 has been conducted. One main source of the uncertainty of the data obtained from the molecular dynamics simulation is that the Verlet integration method used to solve the mathematical equations describing the motion of the atoms is precise to third-order; that means the local trunction error varies as (∆t ) 4 [27,28] when conducting the numerical integration. For the approximation part of the microscopic heat current, the value of the parameter I is the minimum under the condition that the nanoparticle volume fraction is zero, which means there is no nanoparticle in the base fluid, and it is the maximum in the case that the nanoparticle volume fraction is 4%. The same trend appears for the value of the parameter σ . The parameter

σ increases with the increase of the nanoparticle volume fraction. As to the detail part of the microscopic heat current, the value of the parameter I has it minimum value in the case that there is no nanoparticle in the base fluid, and reaches its maximum under the condition that the nanoparticle volume fraction is 4%. As the nanoparticle volume fraction increases, the parameter σ changes accordingly.

CONCLUSION The microscopic thermal transport inside an argon-copper nanofluid is studied based on equilibrium molecular dynamics simulation. The signal of the microscopic heat current is analyzed based on Haar wavelet, and the signal is decomposed into an approximation part and

ACS Paragon Plus Environment

Page 12 of 23

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

Analytical Chemistry

detail part. The probability distributions of the approximation and the detail parts of microscopic heat current signal are investigated and it is found that the probability distributions show Gaussian-like character. The parameter of the standard deviation of the Gaussian-like curve of the probability distribution increases with the increase of the nanoparticle volume fraction. The signal of the microscopic heat current is a microscopic parameter that reflects the microscopic thermal fluctuation inside the nanofluid system, and the nanoparticle volume fraction is a macroscopic parameter of the nanofluid system. The relationship between the microscopic heat current and the nanoparticle volume fraction reflects the relationship between the microscopic and macroscopic characters of the nanofluid system.

Notes The authors declare no competing financial interest.

References (1) Farbod, M.; Ahangarpour, A.; Etemadb, S. G. Particuology. 2015, 22 ,59-65. (2) Taghi, M.; Z.Moattar, Shekaaria, H.; Rast,R.M.; Cegincara, R. M. Fluid Phase Equilibria. 2015,403 136–144. (3) Singh, R.; Sanchez, O.; Ghosh, S.; Kadimcherla,N.; Sen,S.; Balasubramanian, G. Physics Letters A. 2015, 379, 2641-2644. (4) Chinnam, J.; Das, D.K.; Vajjha, R.S.; Jagannadha, R. S. International Journal of Thermal Sciences. 2015, 98, 68-80. (5) Parametthanuwat,T.; Bhuwakietkumjohn, N.; Rittidech, S.; Ding, Y. International Journal of Heat and Fluid Flow. 2015, 56, 80-90. (6) Sahu, R. K.; Hiremath,S. S.; Manivannan, P.V. Powder Technology. 2015, 284, 429-436. (7) Ahmed, H. E.; Ahmed, M.I.; Yusoff, M.Z.; Hawlader, M.N.A.; Ani, H. A. International Journal of Heat and Mass Transfer. 2015, 90, 1197-1206. (8) Ghanbarpour, M.; Nikkam, N.; Khodabandeh, R.; Toprak, M. S. International Communications in Heat and Mass Transfer. 2015, 67, 14-20.

ACS Paragon Plus Environment

Analytical 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

(9) Lee, J.; Lee, H.; Baik,Y. J.; Koo, J. International Journal of Heat and Mass Transfer. 2015, 89,116-12. (10) Colla, L.; Fedele, L.; Buschmann, M.H. International Journal of Thermal Sciences. 2015, 97, 26-40. (11) Vaujour, E. P.; Stannard, A.; Martin, C. P.; Blunt, M. O.; Notingher, I.; Moriarty, P. J. Physical Review Letters. 2008, 100, 176102. (12) Patel, R. Physical Review E. 2012, 85, 026316. (13) Bernardin, M.; Comitani, F.; Vailati, A. Physical Review E. 2012, 85, 066321. (14) Kondaraju, S.; Jin, E. K.; Lee, J. S. Physical Review E. 2010, 81, 016304. (15) Giavazzi , F.; Vailati, A. Physical Review E. 2009, 80, 015303R. (16) Sarafraz, M. M.; Hormozi, F. Powder Technology. 2016, 287, 412-430. (17) Sarafraz, M. M. Heat Mass Transfer. 2012, 48, 611–619.

(18) Salari, E.; Peyghambarzadeh, S. M.; Sarafraz, M. M.; Hormozi, F. Periodica Polytechnica Chemical Engineering. 2016, 60, 106-122. (19) Sarafraz, M. M.; Nikkhah, V.; Nakhjavani, M . Experimental Thermal & Fluid Science. 2018, 59, 509-519. (20) Nakhjavani , M.; Nikkhah, V.; Sarafraz, M. M.; Shoja, S.; Sarafraz, M. Heat & Mass Transfer. 2017, 53, 1-9. (21) Kamalgharibi, M.; Hormozi , F.; Zamzamian, S.; Sarafraz, M. M. Heat & Mass Transfer. 2016, 52, 55-62. (22) Sarafraz, M. M.; Arya, A.; Hormozi ,F.; Nikkhah, V. Applied Thermal Engineering. 2017, 112, 1373-1381. (23) Sarafraz, M. M.; Nikkhah, V.; Nakhjavani , M.; Arya, A. Applied Thermal Engineering. 2017, 123, 29-39. (24) Salari, E.; Peyghambarzadeh, M.; Sarafraz, M. M.; Hormozi, F. Periodica Polytechnica Chemical Engineering, 2016, 60, 252-258. (25) Jia, T. ; Gao, D. Physical Chemistry Chemical Physics. 2016,18, 21930.

ACS Paragon Plus Environment

Page 14 of 23

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

Analytical Chemistry

(26) Jia, T.; Gao, D.; Ying, Y. Y. Journal of Nanofluids. 2014, 3, 1–9. (27) Rapaport, D. C. The Art of Molecular Dynamics Simulation, Cambridge University Press

1995. (28) Haile, J. M. Molecular Dynamics Simulation: Elementary Methods, John Wiley &Sons, Inc

1992. (29) Patel, V. K.; Singh, S.; Singh, V. K. Journal of Computational and Applied Mathematics.

2017, 317, 307–330. (30) Sahu, P. K.; Ray, S. S. Fuzzy Sets and Systems. 2017, 309,131–144. (31) Akbari, Y.; Nouri, K.;

Sadri, J.; Djeddi, C.; Siddiqi, I. Image and Vision Computing.

2017, 59, 17–30. (32) Postnikov, E. B.; Stiukhina, E. S.; Postnov, D. E. Applied Mathematics and Computation.

2017, 305,251–26.

ACS Paragon Plus Environment

Analytical 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

Table 1 Parameters for Curve Fitting to Approximate Part of Microscopic Heat Current

NVF

I (Approximation part)

Sig(Approximation part)

C(Approximation part)

R2

0%

0.055

0.12

0

0.9927

2%

0.07

0.28

0

0.9916

2.5%

0.065

0.29

0

0.9980

3%

0.063

0.3

0

0.9969

3.5%

0.063

0.31

0

0.9971

4%

0.075

0.32

0

0.9973

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

Analytical Chemistry

Table 2 Parameters for Curve Fitting to Detail Part of Microscopic Heat Current C( Detail part )

R2

NVF

I (Detail part)

Sig( Detail part )

0%

0.06

0.002

0

0.9947

2%

0.07

0.015

0

0.9952

2.5%

0.062

0.017

0

0.9974

3%

0.062

0.019

0

0.9985

3.5%

0.07

0.02

0

0.9958

4%

0.073

0.024

0

0.9962

ACS Paragon Plus Environment

Analytical 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

Figure caption Figure 1 Copper nanoparticles distributed in argon fluid

Figure 2 Microscopic heat currents inside nanofluids

Figure 3 Approximation parts of the signals of the microscopic heat currents inside the nanofluids

Figure 4 Detail parts of the signals of the microscopic heat currents inside the nanofluids

ACS Paragon Plus Environment

Page 18 of 23

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

Analytical Chemistry

Figure 1

ACS Paragon Plus Environment

0

1

2 3 t ( second )

4 x 10

1 0 -1

0

1

2 3 t ( second )

4

0

1

2 3 t ( second )

-1

0

1

4

2 3 t ( second )

4 -10

x 10

1 0 -1

0

1

x 10

0

-1

0

-10

1

Page 20 of 23

1

-10

J* ( NVF = 3% )

0

J* ( NVF = 4% )

J* ( NVF = 2.5% )

-0.5

J* ( NVF = 2% )

0.5

J* ( NVF = 3.5% )

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

J* ( NVF = 0% )

Analytical Chemistry

2 3 t ( second )

4 -10

x 10

2

0

-2

0

-10

x 10

1

2 3 t ( second )

Figure 2

ACS Paragon Plus Environment

4 -10

x 10

Page 21 of 23

NVF = 0 %

NVF = 2 % Probability

0.05 0 -0.5

Probability

0.1

0 J* approximation part NVF = 2.5 %

0.1 0.05 0 -2

-1 0 1 J* approximation part NVF = 3.5 %

2

-1 0 1 J* approximation part NVF = 3 %

2

-1 0 1 J* approximation part NVF = 4 %

2

-1 0 1 J* approximation part

2

0.1 0.05 0 -2

0.1 Probability

0.1 0.05 0 -2

0.05 0 -2

0.5

Probability

Probability

0.1

Probability

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

Analytical Chemistry

-1 0 1 J* approximation part

2

0.05 0 -2

Figure 3

ACS Paragon Plus Environment

Analytical Chemistry

NVF = 0 %

NVF = 2 % Probability

0.1

0.05 0 -0.01 -0.005 0 0.005 J* detail part NVF = 2.5 % 0.1 0.05 0 -0.1

-0.05 0 0.05 J* detail part NVF = 3.5 %

0.1

-0.05 0 0.05 J* detail part NVF = 3 %

0.1

-0.05 0 0.05 J* detail part NVF = 4 %

0.1

-0.05 0 0.05 J* detail part

0.1

0.1 0.05 0 -0.1

0.1 Probability

0.1 0.05 0 -0.1

0.05 0 -0.1

0.01

Probability

Probability

Probability

0.1

Probability

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

-0.05 0 0.05 J* detail part

0.1

0.05 0 -0.1

Figure 4

ACS Paragon Plus Environment

Page 23 of 23 1 2 3 4 5 6 7 8 9 10 11

Analytical Chemistry

ACS Paragon Plus Environment