No Arabic abstract
Criterion for contacting is critically important for the Generalized Interpolation Material Point(GIMP) method. We present an improved criterion by adding a switching function. With the method dynamical response of high melting explosive(HMX) with cavities under shock is investigated. The physical model used in the present work is an elastic-to-plastic and thermal-dynamical model with Mie-Gruneissen equation of state. We mainly concern the influence of various parameters, including the impacting velocity $v$, cavity size $R$, etc, to the dynamical and thermodynamical behaviors of the material. For the colliding of two bodies with a cavity in each, a secondary impacting is observed. Correspondingly, the separation distance $D$ of the two bodies has a maximum value $D_{max}$ in between the initial and second impacts. When the initial impacting velocity $v$ is not large enough, the cavity collapses in a nearly symmetric fashion, the maximum separation distance $D_{max}$ increases with $v$. When the initial shock wave is strong enough to collapse the cavity asymmetrically along the shock direction, the variation of $D_{max}$ with $v$ does not show monotonic behavior. Our numerical results show clear indication that the existence of cavities in explosive helps the creation of ``hot spots.
The collapse of cavities under shock is a key problem in various fields ranging from erosion of material, ignition of explosive, to sonoluminescence, etc. We study such processes using the material-point-method developed recently in the field of solid physics. The main points of the research include the relations between symmetry of collapsing and the strength of shock, other coexisting interfaces, as well as hydrodynamic and thermal-dynamic behaviors ignored by the pure fluid models. In the case with strong shock, we study the procedure of jet creation in the cavity; in the case with weak shock, we found that the cavity can not be collapsed completely by the shock and the cavity may collapse in a nearly isotropic way. The history of collapsing significantly influences the distribution of hot spots in the shocked material. The change in symmetry of collapsing is investigated. Since we use the Mie-Gr% {u}neisen equation of state and the effects of strain rate are not taken into account, the behavior is the same if one magnifies the spatial and temporal scales in the same way.
Melting of orthorhombic boron silicide B6Si has been studied at pressures up to 8 GPa using in situ electrical resistivity measurements and quenching. It has been found that in the 2.6-7.7 GPa range B6Si melts congruently, and the melting curve exhibits negative slope of -31(2) K/GPa that points to a higher density of the melt as compared to the solid phase. At very high temperatures B6Si melt appears to be unstable and undergoes disproportionation into silicon and boron-rich silicides. The onset temperature of disproportionation strongly depends on pressure, and the corresponding low-temperature boundary exhibits negative slope of -92(3) K/GPa which is indicative of significant volume decrease in the course of B6Si melt decomposition.
Direct modeling of porous materials under shock is a complex issue. We investigate such a system via the newly developed material-point method. The effects of shock strength and porosity size are the main concerns. For the same porosity, the effects of mean-void-size are checked. It is found that, local turbulence mixing and volume dissipation are two important mechanisms for transformation of kinetic energy to heat. When the porosity is very small, the shocked portion may arrive at a dynamical steady state; the voids in the downstream portion reflect back rarefactive waves and result in slight oscillations of mean density and pressure; for the same value of porosity, a larger mean-void-size makes a higher mean temperature. When the porosity becomes large, hydrodynamic quantities vary with time during the whole shock-loading procedure: after the initial stage, the mean density and pressure decrease, but the temperature increases with a higher rate. The distributions of local density, pressure, temperature and particle-velocity are generally non-Gaussian and vary with time. The changing rates depend on the porosity value, mean-void-size and shock strength. The stronger the loaded shock, the stronger the porosity effects. This work provides a supplement to experiments for the very quick procedures and reveals more fundamental mechanisms in energy and momentum transportation.
{em Ab initio} techniques based on density functional theory in the projector-augmented-wave implementation are used to calculate the free energy and a range of other thermodynamic properties of liquid iron at high pressures and temperatures relevant to the Earths core. The {em ab initio} free energy is obtained by using thermodynamic integration to calculate the change of free energy on going from a simple reference system to the {em ab initio} system, with thermal averages computed by {em ab initio} molecular dynamics simulation. The reference system consists of the inverse-power pair-potential model used in previous work. The liquid-state free energy is combined with the free energy of hexagonal close packed Fe calculated earlier using identical {em ab initio} techniques to obtain the melting curve and volume and entropy of melting. Comparisons of the calculated melting properties with experimental measurement and with other recent {em ab initio} predictions are presented. Experiment-theory comparisons are also presented for the pressures at which the solid and liquid Hugoniot curves cross the melting line, and the sound speed and Gr{u}neisen parameter along the Hugoniot. Additional comparisons are made with a commonly used equation of state for high-pressure/high-temperature Fe based on experimental data.
We investigate the approach to catastrophic failure in a model porous granular material undergoing uniaxial compression. A discrete element computational model is used to simulate both the micro-structure of the material and the complex dynamics and feedbacks involved in local fracturing and the production of crackling noise. Under strain-controlled loading micro-cracks initially nucleate in an uncorrelated way all over the sample. As loading proceeds the damage localizes into a narrow damage band inclined at 30-45 degrees to the load direction. Inside the damage band the material is crushed into a poorly-sorted mixture of mainly fine powder hosting some larger fragments. The mass probability density distribution of particles in the damage zone is a power law of exponent 2.1, similar to a value of 1.87 inferred from observations of the length distribution of wear products (gouge) in natural and laboratory faults. Dynamic bursts of radiated energy, analogous to acoustic emissions observed in laboratory experiments on porous sedimentary rocks, are identified as correlated trails or cascades of local ruptures that emerge from the stress redistribution process. As the system approaches macroscopic failure consecutive bursts become progressively more correlated. Their size distribution is also a power law, with an equivalent Gutenberg-Richter b-value of 1.22 averaged over the whole test, ranging from 3 down to 0.5 at the time of failure, all similar to those observed in laboratory tests on granular sandstone samples. The formation of the damage band itself is marked by a decrease in the average distance between consecutive bursts and an emergent power law correlation integral of event locations with a correlation dimension of 2.55, also similar to those observed in the laboratory (between 2.75 and 2.25).