Viscosity calculations at molecular dynamics simulations
Viscosity and diffusion are chosen as an example to demonstrate the universality of diagnostics methods in the molecular dynamics method. To emphasize the universality, three diverse systems are investigated, which differ from each other drastically: liquids with embedded atom method and pairwise interatomic interaction potentials and dusty plasma with a unique multiparametric interparticle interaction potential. Both the Einstein–Helfand and Green–Kubo relations are used. Such a particular process as glass transition is analysed at the simulation of the aluminium melt. The effect of the dust particle charge fluctuation is considered. The results are compared with the experimental data.
Wave processes occurring under the interaction of the Earth's magnetosphere with dusty plasma near the lunar surface are studied. Ion-acoustic waves are shown to be excited in some regions of the magnetosphere due to the development of a linear hydrodynamic instability. This results in the excitation of ion-acoustic turbulence in these regions. Dust-acoustic waves are demonstrated to be generated due to the development of linear kinetic instability in the entire region of magnetotail interaction with dusty plasma near the Moon. Correspondingly, dust-acoustic turbulence can be excited in the entire region of the interaction of the Earth's magnetosphere with dusty plasma near the lunar surface. We discuss magnetic reconnection processes, which are related to the development of plasma turbulence at the Moon.
The quantum nuclear effects are studied in water using centroid molecular dynamics (CMD) method. The aim is the calibration of CMD implementation in LAMMPS. The calculated intramolecular energy, atoms gyration radii and radial distribution functions are shown in comparison with previous works. The work is assumed to be the step toward to solution of the discrepancy between the simulation results and the experimental data of liquid n-alkane properties in our previous works.
It is well recognized that excessive ice accumulation at low-temperature conditions can cause significant damage to civil infrastructure. The passive anti-icing surfaces provide a promising solution to suppress ice nucleation and enhance ice removal. However, despite extensive efforts, it remains a challenge to design anti-icing surfaces with low ice adhesion. Using all-atom molecular dynamics (MD) simulations, we show that surfaces with single-walled carbon nanotube array (CNTA) significantly reduce ice adhesion due to the extremely low solid areal fraction. It was found that the CNTA surface exhibits up to a 45% decrease in the ice adhesion strength in comparison with the atomically smooth graphene surface. The details of the ice detachment from the CNTA surface were examined for different water-carbon interaction energies and temperatures of the ice cube. Remarkably, the results of MD simulations demonstrate that the ice detaching strength depends linearly on the ratio of the ice-surface interaction energy and the ice temperature. These results open the possibility for designing novel robust surfaces with low ice adhesion for passive anti-icing applications.
Molecular-level knowledge of the thermodynamic, structural, and transport properties of water confined by interfaces and nanopores of various materials is crucial for quantitative understanding and prediction of many natural and technological processes, including carbon sequestration, water desalination, nuclear waste storage, cement chemistry, fuel cell technology, etc. Computational molecular modeling is capable to significantly complement the experimental investigations of such systems by providing invaluable atomic-scale information leading to improved understanding of the specific effects of the substrate structure and composition on the structure, dynamics and reactivity of interfacial and nano-confined aqueous solutions. This paper offers a brief overview of recent efforts to quantify some of these effects for individual H2O molecules and hydrated ions confined at the interfaces and in nanopores of several typical hydrophilic and hydrophobic materials. The first molecular layer of aqueous solution at all substrates is often highly ordered, indicating reduced translational and orientational mobility of the H2O molecules. This ordering cannot be simply described as “ice-like”, but rather resembles the behavior of supercooled water or amorphous ice, although with very significant substrate-specific variations.
In this work, we perform coexistence simulations of methane hydrates for pressures up to 5000 bar for different water models. We calculate the kinetic stability boundary of the superheated metastable sI structure and analyze the effects of the heating rate, system size and cage occupancy. We also report molecular dynamics simulation of several possible structuresfor the new hydrogen hydrate clathrate. We show the strength of molecular simulation as a supplement tool for the analysis of experimental data. © 2015 by Nova Science Publishers, Inc. All rights reserved.
Cooling of tokamak boundary plasma owing to radiation of non-fully stripped lithium ions is considered as a promising way for protection of plasma facing elements (PFE) in tokamak. It may be effectively realized when the main part of lithium ions are involved in the closed circuit of migration between plasma and PFE surface. Such an approach may be implemented with the use of lithium device whose hot (500-600 °C) area to be effected by plasma serves as a Li-emitter and the cold part (∼180 °C) as a Li-collector in the shadow. Capillary-pore system (CPS) provides the returning of collected and condensed lithium to emitting zone by capillary forces. The main goals of the last T-11M lithium experiments were investigating Li ions transport in the tokamak scrape of layer (SOL) and their collecting by different kinds of limiters. The design of devices based on lithium CPS with different ratio of emitting/collecting area is the main subject of this paper. © 2015 The Authors.
The dynamics of a two-component Davydov-Scott (DS) soliton with a small mismatch of the initial location or velocity of the high-frequency (HF) component was investigated within the framework of the Zakharov-type system of two coupled equations for the HF and low-frequency (LF) fields. In this system, the HF field is described by the linear Schrödinger equation with the potential generated by the LF component varying in time and space. The LF component in this system is described by the Korteweg-de Vries equation with a term of quadratic influence of the HF field on the LF field. The frequency of the DS soliton`s component oscillation was found analytically using the balance equation. The perturbed DS soliton was shown to be stable. The analytical results were confirmed by numerical simulations.
Radiation conditions are described for various space regions, radiation-induced effects in spacecraft materials and equipment components are considered and information on theoretical, computational, and experimental methods for studying radiation effects are presented. The peculiarities of radiation effects on nanostructures and some problems related to modeling and radiation testing of such structures are considered.