Strain localization in planar shear of granular media: the role of porosity and boundary conditions.

Shear strain localization into shear bands is associated with velocity weakening instabilities and earthquakes. Here, we simulate steady-state plane-shear flow of numerical granular material (gouge), confined between parallel surfaces. Both constant shear stress and constant strain-rate boundary conditions are tested, and the two types of boundary conditions are found to yield distinct velocity profiles and friction laws. The inertial number, I, exerts the largest control on the layers' behavior, but additional dependencies of friction on normal stress and thickness of the layer are observed under constant stress boundary condition. We find that shear-band localization, which is present in the quasistatic regime ([Formula: see text]) in rate-controlled shear, is absent under stress-controlled loading. In the latter case, flow ceases when macroscopic friction coefficient approaches the quasistatic friction value. The inertial regime that occurs at higher inertial numbers ([Formula: see text]) is associated with distributed shear, and friction and porosity that increase with shear rate (rate-strengthening regime). The finding that shear under constant stress boundary condition produces the inertial, distributed shear but never quasistatic, localized deformation is rationalized based on low fluctuations of shear forces in granular contacts for stress-controlled loading. By examining porosity within and outside a shear band, we also provide a mechanical reason why the transition between quasistatic and inertial shear coincides with the transition between localized and distributed strain.

Revealing isomerism in sodium-water clusters: Photoionization spectra of Na(H2O)n (n = 2-90).

Soft ionization of sodium tagged polar clusters is increasingly used as a powerful technique for sizing and characterization of small aerosols with possible application, e.g., in atmospheric chemistry or combustion science. Understanding the structure and photoionization of the sodium doped clusters is critical for such applications. In this work, we report on measurements of photoionization spectra for sodium doped water clusters containing 2-90 water molecules. While most of the previous studies focused on the ionization threshold of the Na(H2O)n clusters, we provide for the first time full photoionization spectra, including the high-energy region, which are used as reference for a comparison with theory. As reported in previous work, we have seen an initial drop of the appearance ionization energy with cluster size to values of about 3.2 eV for n<5. In the size range from n = 5 to n = 15, broad ion yield curves emerge; for larger clusters, a constant range between signal appearance (∼2.8 eV) and signal saturation (∼4.1 eV) has been observed. The measurements are interpreted with ab initio calculations and ab initio molecular dynamics simulations for selected cluster sizes (n≤ 15). The simulations revealed theory shortfalls when aiming at quantitative agreement but allowed us identifying structural motifs consistent with the observed ionization energy distributions. We found a decrease in the ionization energy with increasing coordination of the Na atom and increasing delocalization of the Na 3s electron cloud. The appearance ionization energy is determined by isomers with fully solvated sodium and a highly delocalized electron cloud, while both fully and incompletely solvated isomers with localized electron clouds can contribute to the high energy part of the photoionization spectrum. Simulations at elevated temperatures show an increased abundance of isomers with low ionization energies, an entropic effect enabling size selective infrared action spectroscopy, based on near threshold photoionization of Na(H2O)n clusters. In addition, simulations of the sodium pick-up process were carried out to study the gradual formation of the hydrated electron which is the basis of the sodium-tagging sizing.

Unsteady granular flows down an inclined plane.

The continuum description of granular flows is still a challenge despite their importance in many geophysical and industrial applications. We extend previous works, which have explored steady flow properties, by focusing on unsteady flows accelerating or decelerating down an inclined plane in the simple shear configuration. We solve the flow kinematics analytically, including predictions of evolving velocity and stress profiles and the duration of the transient stage. The solution shows why and how granular materials reach steady flow on slopes steeper than the angle of repose and how they decelerate on shallower slopes. The model might facilitate development of natural hazard assessment and may be modified in the future to explore unsteady granular flows in different configurations.

Mutual diffusion in the ternary mixture of water + methanol + ethanol and its binary subsystems.

Mutual diffusion is investigated by means of experiment and molecular simulation for liquid mixtures containing water + methanol + ethanol. The Fick diffusion coefficient is measured by Taylor dispersion as a function of composition for all three binary subsystems under ambient conditions. For the aqueous systems, these data compare well with literature values. In the case of methanol + ethanol, experimental measurements of the Fick diffusion coefficient are presented for the first time. The Maxwell-Stefan diffusion coefficient and the thermodynamic factor are predicted for the ternary mixture as well as its binary subsystems by molecular simulation in a consistent manner. The resulting Fick diffusion coefficient is compared to present measurements and that obtained from the classical simulation approach, which requires experimental vapor-liquid equilibrium or excess enthalpy data. Moreover, the self-diffusion coefficients and the shear viscosity are predicted by molecular dynamics and are favorably compared to experimental literature values. The presented ternary diffusion data should facilitate the development of aggregated predictive models for diffusion coefficients of polar and hydrogen-bonding systems.

Determination of the distance-dependent viscosity of mixtures in parallel slabs using non-equilibrium molecular dynamics.

We generalize a technique for determination of the shear viscosity of mixtures in planar slabs using non-equilibrium computer simulations by applying an external force parallel to the surface generating Poiseuille flow. The distance-dependent viscosity of the mixture, given as a function of the distance from the surface, is determined by analysis of the resulting velocity profiles of all species. We present results for a highly non-ideal water + methanol mixture in the whole concentration range between rutile (TiO(2)) walls. The bulk results are compared to the existing equilibrium molecular dynamics and experimental data while the inhomogeneous viscosity profiles at the interface are interpreted using the structural data and information on hydrogen bonding.