1,332 research outputs found
Molecular-level simulations of turbulence and Its decay
We provide the first demonstration that molecular-level methods based on gas kinetic theory and molecular chaos can simulate turbulence and its decay. The direct simulation Monte Carlo (DSMC) method, a molecular-level technique for simulating gas flows that resolves phenomena from molecular to hydrodynamic (continuum) length scales, is applied to simulate the Taylor-Green vortex flow. The DSMC simulations reproduce the Kolmogorov − 5 / 3 law and agree well with the turbulent kinetic energy and energy dissipation rate obtained from direct numerical simulation of the Navier-Stokes equations using a spectral method. This agreement provides strong evidence that molecular-level methods for gases can be used to investigate turbulent flows quantitatively
Gas-kinetic simulation of sustained turbulence in minimal Couette flow
We provide a demonstration that gas-kinetic methods incorporating molecular chaos can simulate the sustained turbulence that occurs in wall-bounded turbulent shear flows. The direct simulation Monte Carlo method, a gas-kinetic molecular method that enforces molecular chaos for gas-molecule collisions, is used to simulate the minimal Couette flow at Re=500. The resulting law of the wall, the average wall shear stress, the average kinetic energy, and the continually regenerating coherent structures all agree closely with corresponding results from direct numerical simulation of the Navier-Stokes equations. These results indicate that molecular chaos for collisions in gas-kinetic methods does not prevent development of molecular-scale long-range correlations required to form hydrodynamic-scale turbulent coherent structures
DSMC simulations of turbulent flows at moderate Reynolds numbers
The Direct Simulation Monte Carlo (DSMC) method has been used for more than 50 years to simulate rarefied gases. The advent of modern supercomputers has brought higher-density near-continuum flows within range. This in turn has revived the debate as to whether the Boltzmann equation, which assumes molecular chaos, can be used to simulate continuum flows when they become turbulent. In an effort to settle this debate, two canonical turbulent flows are examined, and the results are compared to available continuum theoretical and numerical results for the Navier-Stokes equations
Solute effects on edge dislocation pinning in complex alpha-Fe alloys
Reactor pressure vessel steels are well-known to harden and embrittle under neutron irradiation, mainly because of the formation of obstacles to the motion of dislocations, in particular, precipitates and clusters composed of Cu, Ni, Mn, Si and P. In this paper, we employ two complementary atomistic modelling techniques to study the heterogeneous precipitation and segregation of these elements and their effects on the edge dislocations in BCC iron. We use a special and highly computationally efficient Monte Carlo algorithm in a constrained semi-grand canonical ensemble to compute the equilibrium configurations for solute clusters around the dislocation core. Next, we use standard molecular dynamics to predict and analyze the effect of this segregation on the dislocation mobility. Consistently with expectations our results confirm that the required stress for dislocation unpinning from the precipitates formed on top of it is quite large. The identification of the precipitate resistance allows a quantitative treatment of atomistic results, enabling scale transition towards larger scale simulations, such as dislocation dynamics or phase field.Fil: Pascuet, Maria Ines Magdalena. Consejo Nacional de Investigaciones Científicas y Técnicas; Argentina. Comision Nacional de Energia Atomica. Centro Atomico Constituyentes. Departamento de Materiales; ArgentinaFil: Martínez, E.. Los Alamos National High Magnetic Field Laboratory; Estados UnidosFil: Monnet, G.. EDF–R&D; FranciaFil: Malerba, L.. SCK•CEN. Structural Materials Expert Group. Nuclear Materials Institute; Bélgic
Lowest matric potential in quartz: Metadynamics evidence
The lowest matric potential is an important soil property characterizing the strength of retaining water molecules and a key parameter in defining a complete soil water retention curve. However, the exact value of the lowest matric potential is still unclear and cannot be measured due to the limitation of current experimental technology. In this study, a general theoretical framework based on metadynamics was proposed to determine the lowest matric potential in quartz minerals. The matric potential was derived from partial volume free energy and can be further calculated by the difference between the adsorption free energy and self‐hydration free energy. Metadynamics was employed to enhance molecular dynamics for determination of the adsorption free energy. In addition to the water‐mineral interaction, the adsorptive water layer structure was identified as an important mechanism that may lower the free energy of water molecules. The lowest matric potential for quartz mineral was found as low as −2.00 GPa
Segue Between Favorable and Unfavorable Solvation
Solvation of small and large clusters are studied by simulation, considering
a range of solvent-solute attractive energy strengths. Over a wide range of
conditions, both for solvation in the Lennard-Jones liquid and in the SPC model
of water, it is shown that the mean solvent density varies linearly with
changes in solvent-solute adhesion or attractive energy strength. This behavior
is understood from the perspective of Weeks' theory of solvation [Ann. Rev.
Phys. Chem. 2002, 53, 533] and supports theories based upon that perspective.Comment: 8 pages, 7 figure
Dynamical heterogeneities in a supercooled Lennard-Jones liquid
We present the results of a large scale molecular dynamics computer
simulation study in which we investigate whether a supercooled Lennard-Jones
liquid exhibits dynamical heterogeneities. We evaluate the non-Gaussian
parameter for the self part of the van Hove correlation function and use it to
identify ``mobile'' particles. We find that these particles form clusters whose
size grows with decreasing temperature. We also find that the relaxation time
of the mobile particles is significantly shorter than that of the bulk, and
that this difference increases with decreasing temperature.Comment: 8 pages of RevTex, 4 ps figure
Boron Nitride Monolayer: A Strain-Tunable Nanosensor
The influence of triaxial in-plane strain on the electronic properties of a
hexagonal boron-nitride sheet is investigated using density functional theory.
Different from graphene, the triaxial strain localizes the molecular orbitals
of the boron-nitride flake in its center depending on the direction of the
applied strain. The proposed technique for localizing the molecular orbitals
that are close to the Fermi level in the center of boron nitride flakes can be
used to actualize engineered nanosensors, for instance, to selectively detect
gas molecules. We show that the central part of the strained flake adsorbs
polar molecules more strongly as compared with an unstrained sheet.Comment: 20 pages, 9 figure
Time and length scales in supercooled liquids
We numerically obtain the first quantitative demonstration that development
of spatial correlations of mobility as temperature is lowered is responsible
for the ``decoupling'' of transport properties of supercooled liquids. This
result further demonstrates the necessity of a spatial description of the glass
formation and therefore seriously challenges a number of popular alternative
theoretical descriptions.Comment: 4 pages, 4 figs; improved version: new refs and discussion
Fluctuations of water near extended hydrophobic and hydrophilic surfaces
We use molecular dynamics simulations of the SPC-E model of liquid water to
derive probability distributions for water density fluctuations in probe
volumes of different shapes and sizes, both in the bulk as well as near
hydrophobic and hydrophilic surfaces. To obtain our results, we introduce a
biased sampling of coarse-grained densities, which in turn biases the actual
solvent density. The technique is easily combined with molecular dynamics
integration algorithms. Our principal result is that the probability for
density fluctuations of water near a hydrophobic surface, with or without
surface-water attractions, is akin to density fluctuations at the water-vapor
interface. Specifically, the probability of density depletion near the surface
is significantly larger than that in bulk. In contrast, we find that the
statistics of water density fluctuations near a model hydrophilic surface are
similar to that in the bulk
- …