On the dependence of the Navier Stokes equations on the distribution of moleular velocities

In this work we introduce a completely general Chapman Enskog procedure in which we divide the local distribution into an isotropic distribution with anisotropic corrections. We obtain a recursion relation on all integrals of the distribution function required in the derivation of the moment equations. We obtain the hydrodynamic equations in terms only of the first few moments of the isotropic part of an arbitrary local distribution function. The incompressible limit of the equations is completely independent of the form of the isotropic part of the distribution, whereas the energy equation in the compressible case contains an additional contribution to the heat flux. This additional term was also found by Boghosian and by Potiguar and Costa in the derivation of the Navier Stokes equations for Tsallis thermostatistics, and is the only additional term allowed by the Curie principle

Degeneracy Algorithm for Random Magnets

It has been known for a long time that the ground state problem of random magnets, e.g. random field Ising model (RFIM), can be mapped onto the max-flow/min-cut problem of transportation networks. I build on this approach, relying on the concept of residual graph, and design an algorithm that I prove to be exact for finding all the minimum cuts, i.e. the ground state degeneracy of these systems. I demonstrate that this algorithm is also relevant for the study of the ground state properties of the dilute Ising antiferromagnet in a constant field (DAFF) and interfaces in random bond magnets.Comment: 17 pages(Revtex), 8 Postscript figures(5color) to appear in Phys. Rev. E 58, December 1st (1998

Spinodal Decomposition in Binary Gases

We carried out three-dimensional simulations, with about 1.4 million particles, of phase segregation in a low density binary fluid mixture, described mesoscopically by energy and momentum conserving Boltzmann-Vlasov equations. Using a combination of Direct Simulation Monte Carlo(DSMC) for the short range collisions and a version of Particle-In-Cell(PIC) evolution for the smooth long range interaction, we found dynamical scaling after the ratio of the interface thickness(whose shape is described approximately by a hyperbolic tangent profile) to the domain size is less than ~0.1. The scaling length R(t) grows at late times like t^alpha, with alpha=1 for critical quenches and alpha=1/3 for off-critical ones. We also measured the variation of temperature, total particle density and hydrodynamic velocity during the segregation process.Comment: 11 pages, Revtex, 4 Postscript figures, submitted to PR

A scaling theory of 3D spinodal turbulence

A new scaling theory for spinodal decomposition in the inertial hydrodynamic regime is presented. The scaling involves three relevant length scales, the domain size, the Taylor microscale and the Kolmogorov dissipation scale. This allows for the presence of an inertial "energy cascade", familiar from theories of turbulence, and improves on earlier scaling treatments based on a single length: these, it is shown, cannot be reconciled with energy conservation. The new theory reconciles the t^{2/3} scaling of the domain size, predicted by simple scaling, with the physical expectation of a saturating Reynolds number at late times.Comment: 5 pages, no figures, revised version submitted to Phys Rev E Rapp Comm. Minor changes and clarification

Observation of off-Hugoniot shocked states with ultrafast time resolution

We apply ultrafast single shot interferometry to determine the pressure and density of argon shocked from up to 7.8 GPa static initial pressure in a diamond anvil cell. This method enables the observation of thermodynamic states distinct from those observed in either single shock or isothermal compression experiments, and the observation of ultrafast dynamics in shocked materials. We also present a straightforward method for interpreting ultrafast shock wave data which determines the index of refraction at the shock front, and the particle and shock velocities for shock waves in transparent materials. Based on these methods, we observe shocked thermodynamic states between the room temperature isotherm of argon and the shock adiabat of cryogenic argon at final shock pressures up to 28 GPa

Percolation in three-dimensional random field Ising magnets

The structure of the three-dimensional random field Ising magnet is studied by ground state calculations. We investigate the percolation of the minority spin orientation in the paramagnetic phase above the bulk phase transition, located at [Delta/J]_c ~= 2.27, where Delta is the standard deviation of the Gaussian random fields (J=1). With an external field H there is a disorder strength dependent critical field +/- H_c(Delta) for the down (or up) spin spanning. The percolation transition is in the standard percolation universality class. H_c ~ (Delta - Delta_p)^{delta}, where Delta_p = 2.43 +/- 0.01 and delta = 1.31 +/- 0.03, implying a critical line for Delta_c < Delta <= Delta_p. When, with zero external field, Delta is decreased from a large value there is a transition from the simultaneous up and down spin spanning, with probability Pi_{uparrow downarrow} = 1.00 to Pi_{uparrow downarrow} = 0. This is located at Delta = 2.32 +/- 0.01, i.e., above Delta_c. The spanning cluster has the fractal dimension of standard percolation D_f = 2.53 at H = H_c(Delta). We provide evidence that this is asymptotically true even at H=0 for Delta_c < Delta <= Delta_p beyond a crossover scale that diverges as Delta_c is approached from above. Percolation implies extra finite size effects in the ground states of the 3D RFIM.Comment: replaced with version to appear in Physical Review

Transport properties of dense fluid argon

We calculate using molecular dynamics simulations the transport properties of realistically modeled fluid argon at pressures up to $\simeq 50GPa$ and temperatures up to $3000K$. In this context we provide a critique of some newer theoretical predictions for the diffusion coefficients of liquids and a discussion of the Enskog theory relevance under two different adaptations: modified Enskog theory (MET) and effective diameter Enskog theory. We also analyze a number of experimental data for the thermal conductivity of monoatomic and small diatomic dense fluids.Comment: 8 pages, 6 figure

A study of tantalum pentoxide Ta2O5 structures up to 28 GPa

Tantalum pentoxide Ta2O5 with the orthorhombic L-Ta2O5 structure has been experimentally studied up to 28.3 GPa (at ambient temperature) using synchrotron angle-dispersive powder X-ray diffraction (XRD). The ambient pressure phase remains stable up to 25 GPa where with increased pressure a crystalline to amorphous phase transition occurs. A detailed equation of state (EOS), including pressure dependent lattice parameters, is reported. The results of this study were compared with a previous high-pressure XRD study by Li et al. A clear discrepancy between the ambient-pressure crystal structures and, consequently, the reported EOSs between the two studies was revealed. The origin of this discrepancy is attributed to two different crystal structures used to index the XRD patterns

Isentropic Compression with a Rectangular Configuration for Tungstene and Tantalum, Computations and Comparison with Experiments

Isentropic compression experiments and numerical simulations on metals are performed at Z accelerator facility from Sandia National Laboratory and at Lawrence Livermore National Laboratory in order to study the isentrope, associated Hugoniot and phase changes of these metals [1]. 3D configurations have been calculated here to benchmark the new beta version of the electromagnetism package coupled with the dynamics in Ls-Dyna and compared with the ICE Z shots 1511 and 1555. The electromagnetism module is being developed in the general-purpose explicit and implicit finite element program LS-DYNA{reg_sign} in order to perform coupled mechanical/thermal/electromagnetism simulations. The Maxwell equations are solved using a Finite Element Method (FEM) for the solid conductors coupled with a Boundary Element Method (BEM) for the surrounding air (or vacuum). More details can be read in the reference [2], [3]