A molecular dynamics simulation on the convective heat transfer in nanochannels

The understanding of the flow and heat transfer processes for fluid through micro- and nanochannels becomes imperative due to its wide application in micro- and nano-fluidic devices. In this paper, the method to simulate the convective heat transfer process in molecular dynamics is improved based on a previous study. With this method, we simulate a warm dense fluid flowing through a cold parallel-plate nanochannel with constant wall temperature. The characteristics of the velocity and temperature fields are analysed. The temperature difference between the bulk average temperature of fluid and the wall temperature decreases in an exponential form along the flow direction. The Nusselt number for the laminar flow in parallel-plate nanochannel is smaller than its corresponding value at macroscale. It could be attributed to the temperature jump at the fluid–wall interface, which decreases the temperature gradient near the wall. The results also reveal that the heat transfer coefficient is related to the surface wettabilities, which differs from that in the macroscopic condition.     

Influence of water on HFO-1234yf oxidation pyrolysis via ReaxFF molecular dynamics simulation

ABSTRACT

The effect of water molecules on HFO-1234yf oxidation pyrolysis was investigated by ReaxFF-molecular dynamics simulation from 1900 to 4200?K. The initial pyrolysis of HFO-1234yf starts around 2500?K and the water molecules participate in chemical reactions at 2800?K when the reactants pyrolysis reached the highest reaction rate. The primary products including HF, COF₂ and CO₂ are observed at 2600, 2700 and 2900?K, respectively. The influence of water molecules on products is mainly reflected in the promotion activity on the conversion from COF₂ to CO₂ and the generation of HF molecules. Four formation pathways are observed and calculated to further elucidate the procedure of pyrolysis. The main conversion process from H₂O to HF is the ?F?+?H₂O?=?HF+?OH reaction, and the paths from H₂O to ?OH radical and COF₂ to ?CFO radical which are promoted by ?F and ?H radical, respectively, have relatively low energy barriers of 10.44 and 40.29?kJ/mol, and both reaction processes released HF molecules.     

Relaxation dynamics and power spectra of liquid water: a molecular dynamics simulation study

We present molecular dynamics simulations of liquid water at normal and supercooled conditions. Autocorrelation functions (ACFs) of several structural quantities and their fourier transforms are obtained and analysed. Structural correlations and relaxation times increase linearly with degree of supercooling. Power spectra of ACFs show increase in librational motion of liquid water with cooling. These modes intensify with supercooling because of structuring and ordering of water molecules. Overall, liquid water structure is homogenous over the temperatures and pressures studied and undergoes fluctuation–dissipation in its local-density variations [English and Tse, Phys. Rev. Lett. 106, 037801 (2011)].     

Study of mechanical properties of amorphous copper with molecular dynamics simulation

The formation and mechanical properties of amorphous copper are studied using molecular dynamics simulation. The simulations of tension and shearing show that more pronounced plasticity is found under shearing, compared to tension. Apparent strain hardening and strain rate effect are observed. Interestingly, the variations of number density of atoms during deformation indicate free volume creation, especially under higher strain rate. In particular, it is found that shear induced dilatation does appear in the amorphous metal.     

A molecular dynamics simulation of nanoscale convective heat transfer with the effect of axial heat conduction

Effective heat dissipation from nano-fluidic devices is sometimes necessary to ensure their performance and lifespan. In the molecular dynamics simulation of nanoscale convective heat transfer, thermostats cannot be directly applied to the fluid because of the non-uniform temperature distribution. Periodic boundary is typically utilised, but unrealistic axial heat conduction exists when there is a temperature difference between the outlet and images of inlet atoms. In this paper, the effect of axial conduction caused by periodic boundary is investigated through the Péclet number (Pe). Taking viscous dissipation into consideration, the magnitude of outlet thermal diffusion is observed to decrease with increasing Pe. The local average temperature of fluid changes in an exponential form except in the region close to the outlet. Results show that the contribution of outlet axial conduction to the local average temperature is less than 2.0% when Pe > 10. The main reason is that the magnitude of fluid velocity and viscous heat dissipation in nanochannels is much larger than that in macro-channels at the same Péclet number.     

Study of lattice thermal conductivity of alpha-zirconium by molecular dynamics simulation

The non-equilibrium molecular dynamics method is adapted to calculate the phonon thermal conductivity of alphazirconium. By exchanging velocities of atoms in different regions, the stable heat flux and the temperature gradient are established to calculate the thermal conductivity. The phonon thermal conductivities under different conditions, such as different heat exchange frequencies, different temperatures, different crystallographic orientations, and crossing grain boundary (GB), are studied in detail with considering the finite size effect. It turns out that the phonon thermal conductivity decreases with the increase of temperature, and displays anisotropies along different crystallographic orientations. The phonon thermal conductivity in [0001] direction (close-packed plane) is largest, while the values in other two directions of [2īī0] and [01ī0] are relatively close. In the region near GB, there is a sharp temperature drop, and the phonon thermal conductivity is about one-tenth of that of the single crystal at 550 K, suggesting that the GB may act as a thermal barrier in the crystal.     

Physical insights into the sonochemical degradation of recalcitrant organic pollutants with cavitation bubble dynamics

This paper tries to discern the mechanistic features of sonochemical degradation of recalcitrant organic pollutants using five model compounds, viz. phenol (Ph), chlorobenzene (CB), nitrobenzene (NB), p-nitrophenol (PNP) and 2,4-dichlorophenol (2,4-DCP). The sonochemical degradation of the pollutant can occur in three distinct pathways: hydroxylation by OH radicals produced from cavitation bubbles (either in the bubble–bulk interfacial region or in the bulk liquid medium), thermal decomposition in cavitation bubble and thermal decomposition at the bubble–liquid interfacial region. With the methodology of coupling experiments under different conditions (which alter the nature of the cavitation phenomena in the bulk liquid medium) with the simulations of radial motion of cavitation bubbles, we have tried to discern the relative contribution of each of the above pathway to overall degradation of the pollutant. Moreover, we have also tried to correlate the predominant degradation mechanism to the physico-chemical properties of the pollutant. The contribution of secondary factors such as probability of radical–pollutant interaction and extent of radical scavenging (or conservation) in the medium has also been identified. Simultaneous analysis of the trends in degradation with different experimental techniques and simulation results reveals interesting mechanistic features of sonochemical degradation of the model pollutants. The physical properties that determine the predominant degradation pathway are vapor pressure, solubility and hydrophobicity. Degradation of Ph occurs mainly by hydroxylation in bulk medium; degradation of CB occurs via thermal decomposition inside the bubble, degradation of PNP occurs via pyrolytic decomposition at bubble interface, while hydroxylation at bubble interface contributes to degradation of NB and 2,4-DCP.     

First-principles molecular dynamics simulation of biased electrode/solution interface

First-principles molecular dynamics simulations have been carried out for water in contact with Pt(1 1 1) surface. To apply negative bias potential to the water/Pt interface, excess electrons were added to our slab model using the recently developed computational scheme called “effective screening medium (ESM)”. Water molecules located away from the surface reoriented themselves to screen the electric field, but they responded differently near the surface. Water molecules nearest to the surface, forming a distinct layered structure with the hydrogen atom directed to the surface, increased the density with increasing field. On these bases, we discuss microscopic aspects of the electric double layer.     

A molecular dynamics study on iridium

In this study, molecular dynamics simulations are performed by using a modified form of Morse potential function in the framework of the Embedded Atom Method (EAM). Temperature-and pressure-dependent behaviours of bulk modulus, second-order elastic constants (SOEC), and the linear-thermal expansion coefficient is calculated and compared with the available experimental data. The melting temperature is estimated from 3 different plots. The obtained results are in agreement with the available experimental findings for iridium.      

A molecular dynamics study of the swelling patterns of Na/Cs-montmorillonites and the hydration of interlayer cations

We report on a molecular dynamics study of the swelling patterns of an Na-rich/Cs-poor montomorillonite and a Csmontomorillonite.The recently developed CLAYFF force field is used to predict the basal spacing as a function of the water content in the interlayer.The simulations reproduce the swelling patterns of the Na and Cs-montomorillonite,suggesting a mechanism of its hydration different from that of the montomorillonite.In addition,we find that the differences in size and hydration energy of Na and Cs ions have strong implications for the structure and the internal energy of interlayer water.In particular,our results indicate that the hydrate difference in the presence of coexistent Na and Cs has a larger influence on the behavior of a clay-water system.For Na-rich/Cs-poor montomorillonite,the hydration energy values of Na ions and water molecules each have a dramatic increase compared with those in Na-montomorillonite on the interlayer spacing,and the hydration energy values of Cs ions and water molecules decrease somewhat compared with those in Cs-montomorillonite.     

Effect of dissolved gases in water on acoustic cavitation and bubble growth rate in 0.83 MHz megasonic of interest to wafer cleaning

Changes in the cavitation intensity of gases dissolved in water, including H₂, N₂, and Ar, have been established in studies of acoustic bubble growth rates under ultrasonic fields. Variations in the acoustic properties of dissolved gases in water affect the cavitation intensity at a high frequency (0.83 MHz) due to changes in the rectified diffusion and bubble coalescence rate. It has been proposed that acoustic bubble growth rates rapidly increase when water contains a gas, such as hydrogen faster single bubble growth due to rectified diffusion, and a higher rate of coalescence under Bjerknes forces. The change of acoustic bubble growth rate in rectified diffusion has an effect on the damping constant and diffusivity of gas at the acoustic bubble and liquid interface. It has been suggested that the coalescence reaction of bubbles under Bjerknes forces is a reaction determined by the compressibility and density of dissolved gas in water associated with sound velocity and density in acoustic bubbles. High acoustic bubble growth rates also contribute to enhanced cavitation effects in terms of dissolved gas in water. On the other hand, when Ar gas dissolves into water under ultrasound field, cavitation behavior was reduced remarkably due to its lower acoustic bubble growth rate. It is shown that change of cavitation intensity in various dissolved gases were verified through cleaning experiments in the single type of cleaning tool such as particle removal and pattern damage based on numerically calculated acoustic bubble growth rates.     

Diffusion activation energy versus the favourable energy in two-order-parameter model:A molecular dynamics study of liquid Al

In the present work, we find that both diffusion activation energy E_a(D) and E_a(S^ex) increase linearly with pressure and have the same slope (0.022±0.001 eV/GPa) for liquid Al. The temperature and pressure dependence of excess entropy is well fitted by the expression -S^ex(T,P)/k_B=a(P)+b(P)T+c(P)exp(E_f/k_BT), which together with the small ratio of E_f/k_BT leads to the relationship of excess entropy to temperature and pressure, i.e. S^ex≈-cE_f/T, where c is about 12 and E_f (=Δ E-PΔV) is the favourable energy. Therefore, there exists a simple relation between E_a(S^ex) and E_f, i.e. E_a(S^ex)≈cE_f.     

Thermal conductive performance of deposited amorphous carbon materials by molecular dynamics simulation

A series of diamond-like carbon (DLC) films with different microstructure were prepared by depositing carbon atoms on diamond surface with incident energy ranging from 1 to 100 eV. The thermal conductivity of the deposited films and the Kapitza resistance between the film and the diamond substrate were investigated. Results show that the average density, the average fraction of sp³ bonding and the thermal conductivity of the DLC films increase first, reaching a maximum around 20–40 eV before decreasing, while the Kapitza resistance decreases gradually with increased deposition energy. The analysis suggests that the thermal resistance of the interface layer is in the order of 10^?10 m²K/W, which is not ignorable when measuring the thermal conductivity of the deposited film especially when the thickness of the DLC film is not large enough. The fraction of sp³ bonding in the DLC film decreases gradually normal to the diamond surface. However, the thermal conductivity of the film in normal direction is not affected obviously by this kind of structural variation but depends linearly on the average fraction of sp³ bonding in the entire film. The dependence of the thermal conductivity on the fraction of sp³ bonding was analysed by the phonon theory.     

Experimental research into the cavitation noise spectrum in water solutions of high molecular weight polymers

The paper presents cavitation noise spectrum measurements results for water and a water solution of polyacryloamide with the 1.5 × 10⁷ g mol⁻¹ particle mass. Cavitation noise spectrum characteristics were shown for different ultrasound intensity, with different sonification times (at a constant ultrasound intensity) as well as characteristics which enable comparing particular sonification parameters.     

Free energy and scalings for polymer translocation through a nanopore: A molecular dynamics simulation study combined with milestoning

Coarse-grained molecular dynamics simulations combined with milestoning method are used to study the stochastic process of polymer chain translocation though a nanopore. We find that the scalings for polymer translocation process (the chain is initialized with the first monomer in the nanopore) and for polymer escape process (the chain is initialized with the middle monomer in the nanopore) are different. The translocation process is mainly controlled by the entropic barrier, while the polymer escape process is driven by the effective force due to free energy difference.     

A molecular dynamics study of the structural change differences between Au225 and Au369 clusters on MgO surfaces at low temperature

The differences in structural change between Au 225 and Au 369 clusters with their(111) facets supported on MgO(100) surfaces at 5 K are studied by using molecular-dynamics simulations with the atomic interchange potentials of the Au/MgO interface.The parameters are obtained from the ab initio energies using the Chen-Mo¨bius inversion method.Analyses of the pair distribution functions show that the two Au clusters use different deformation processes to adjust the distances between the interface atoms,owing to the misfit between the atom distances among the clusters and the substrates.The local structural changes are identified by atomic density profiles.     

A study on dissociation of sII krypton hydrate and the effect of hydrocarbon guest molecules as stabilizer by molecular dynamics simulation

In this paper, the characteristics of structure II krypton hydrate are studied by molecular dynamics simulation under isobaric-isothermal (NPT) ensemble condition. The dissociation process of the hydrate is simulated and the effect of krypton (Kr) and various types of hydrocarbon guest molecules (HGMs) on the stability of the hydrate structure is investigated during the simulation time of 1 ns. The studied HGMs are propane, isobutane, neopentane, cyclopropane, cyclobutane, cyclopentane and cyclopentene. The structural change of the Kr-hydrate is analyzed with the radial distribution function, mean square displacement and diffusion coefficient. As temperature increases, the obtained results indicate a gradual increase in the Kr-hydrate cell size, which leads to distortion of the hydrate lattice and escaping of the encapsulated Kr molecules from the hydrate structure to form small bubbles of Kr aggregated in the aqueous solution.     

Role of oxygen functional groups for structure and dynamics of interfacial water on low rank coal surface: a molecular dynamics simulation

Molecular dynamics simulations were employed to study the effects of oxygen functional groups for structure and dynamics properties of interfacial water molecules on the subbituminous coal surface. Because of complex composition and structure, the graphite surface modified by hydroxyl, carboxyl and carbonyl groups was used to represent the surface model of subbituminous coal according to XPS results, and the composing proportion for hydroxyl, carbonyl and carboxyl is 25:3:5. The hydration energy with ?386.28 kJ/mol means that the adsorption process between water and coal surface is spontaneous. Density profiles for oxygen atoms and hydrogen atoms indicate that the coal surface properties affect the structural and dynamic characteristics of the interfacial water molecules. The interfacial water exhibits much more ordering than bulk water. The results of radial distribution functions, mean square displacement and local self-diffusion coefficient for water molecule related to three oxygen moieties confirmed that the water molecules prefer to absorb with carboxylic groups, and adsorption of water molecules at the hydroxyl and carbonyl is similar.     

Calculation of the thermodynamic properties of copper by molecular dynamics simulation

It is shown that the thermodynamic characteristics of a system can be accurately described using many-body potentials of the interatomic interaction in molecular dynamics calculations. Original Russian Text ? O.V. Stepanyuk, D.B. Alekseev, A.M. Saletskii, 2009, published in Vestnik Moskovskogo Universiteta. Fizika, 2009, No. 2, pp. 115–116.