Atomistic methods, such as molecular dynamics and direct simulation Monte Carlo, constitute a powerful and growing set of techniques for fluid-dynamics simulation. The more fundamental nature of such methods, which exhibit nonlinear transport effects and small-scale fluctuations, extends their modelling accuracy to a significantly wider range of scales and regimes than the more traditional Navier-Stokes-based continuum fluid-simulation techniques. In this paper, we describe the current state of the art in atomistic fluid simulation, from both a theoretical and a computational standpoint, and outline the advantages and limitations of such methods. In addition, we present an overview of some recent atomistic-simulation results on fluid instabilities and on the physical scaling of atomistic techniques. Finally, we suggest possible avenues of future research in the field.

A ubiquitous example of fluid mixing is the Rayleigh-Taylor instability, in which a heavy fluid initially sits atop a light fluid in a gravitational field. The subsequent development of the unstable interface between the two fluids is marked by several stages. At first, each interface mode grows exponentially with time before transitioning to a nonlinear regime characterized by more complex hydrodynamic mixing. Unfortunately, traditional continuum modeling of this process has generally been in poor agreement with experiment. Here, we indicate that the natural, random fluctuations of the flow field present in any fluid, which are neglected in continuum models, can lead to qualitatively and quantitatively better agreement with experiment. We performed billion-particle atomistic simulations and magnetic levitation experiments with unprecedented control of initial interface conditions. A comparison between our simulations and experiments reveals good agreement in terms of the growth rate of the mixing front as well as the new observation of droplet breakup at later times. These results improve our understanding of many fluid processes, including interface phenomena that occur, for example, in supernovae, the detachment of droplets from a faucet, and ink jet printing. Such instabilities are also relevant to the possible energy source of inertial confinement fusion, in which a millimeter-sized capsule is imploded to initiate nuclear fusion reactions between deuterium and tritium. Our results suggest that the applicability of continuum models would be greatly enhanced by explicitly including the effects of random fluctuations.atomistic simulation ͉ magnetic levitation ͉ Rayleigh-Taylor A s early as 400 B.C., the Greek philosopher Democritus of Abdera postulated that matter is comprised of indivisible units, and that it therefore can change its character due to the constant motion and temporary clustering of these ''atoms.'' Today, the existence of atoms can be verified directly via electron microscopy. Therefore, when modeling a sample of material, it is desirable to take into account the dynamics of the atoms, which ultimately are responsible for the macroscopic properties of the gas, liquid, or solid of interest. However, such calculations of the underlying atomistic dynamics have not previously been attempted for a complex fluid flow because it was thought that the number of atoms required was too large for even the most powerful computers. The physical problem has usually been approximated by a continuum description, in which the spatially and temporally averaged properties of a large number of atoms are described. Such descriptions are usually adequate when the scales of interest are large compared with the atomistic scale.A continuum obviously describes many physical situations very well, such as the fluid flow around a car or low-altitude aircraft. However, the continuum hypothesis must fail when the mean distance between atoms is not very small compared with the characteristic length of flow features, as is the...

No abstract

Raman spectra from 50 to 3500 cm(-1) and 4-296 K are analyzed for molecular crystal powders of the explosives pentaerythritol tetranitrate (PETN), beta-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX), and 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) and the inert naphthalene. Temperature-dependent Raman spectroscopy is utilized for its sensitivity to anharmonic couplings between thermally populated phonons and higher frequency vibrations relevant to shock up-pumping. The data are analyzed with anharmonic perturbation theory, which is shown to have significant fundamental limitations in application to real data. Fitting to perturbation theory revealed no significant differences in averaged anharmonicities among the three explosives, all of which exhibited larger averaged anharmonicities than naphthalene by a factor of 3. Calculations estimating the multiphonon densities of states also failed to correlate clearly with shock sensitivity. However, striking differences in temperature-dependent lifetimes were obvious: PETN has long lived phonons and vibrons, HMX has long lived phonons but short lived vibrons, while TATB has short lived phonons and vibrons at low temperature. Naphthalene, widely used as a model system, has significantly different anharmonicities and density of states from any of the explosives. The data presented suggest the further hypothesis that hindered vibrational energy transfer in the molecular crystals is a significant factor in shock sensitivity.

Survey spectra of single-crystal HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine), RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine), and PETN (pentaerythritol tetranitrate) were acquired in the region from 10 to 80 cm(-1) using terahertz time-domain spectroscopy. The spectra were taken at temperatures ranging from 8.4 to 300 K. Generally, the spectra show multiple absorption peaks in the range 50-80 cm(-1), with PETN (110) showing strong absorption features at room temperature. RDX (210) is the most notable in the region 10-40 cm(-1), showing multiple spectral features, while HMX (010) shows a very broad absorption at 47.8 cm(-1) with a fwhm of 37.3 cm(-1). Future plans include polarization-dependent investigations for multiple crystallographic orientations over an increased spectral range and higher-level theoretical calculations.

No abstract

No abstract

No abstract