No Arabic abstract
Ramp compression experiment are used to deduce the relation between compression and normal stress in a material, by measuring how a compression wave evolves as it propagates through different thicknesses of the sample material. The compression wave is generally measured by Doppler velocimetry from a surface that can be observed with optical or near-optical photons. For high-pressure ramp loading, the reflectivity of a free surface often decreases as it is accelerated by the ramp wave, and window materials transparent to the probing photons are used to keep the surface flatter and preserve its reflectivity. We previously described a method of analyzing ramp-wave data measured at the free surface which did not require numerical iteration. However, this method breaks down when the pressure at the surface changes and hence cannot be used for data taken with a finite-impedance window. We have now generalized this non-iterative analysis method to apply to measurements taken through a window. Free surfaces can be treated seamlessly, and the need for sampling at uniform intervals of velocity has been removed. These calculations require interpolation of partially-released states using the partially-constructed stress-compression relation, making them slower than the previous free-surface scheme, but they are still much more robust and fast than iterative analysis.
In the canonical ramp compression experiment, a smoothly-increasing load is applied to the surface of the sample, and the particle velocity history is measured at two or more different distances into the sample, at interfaces where the surface of the sample can be probed. The velocity histories are used to deduce a stress-density relation, usually using iterative Lagrangian analysis to account for the perturbing effect of the impedance mismatch at the interface. In that technique, a stress- density relation is assumed in order to correct for the perturbation, and is adjusted until it becomes consistent with the deduced stress-density relation. This process is subject to the usual difficulties of nonlinear optimization, such as the existence of local minima (sensitivity to the initial guess), possible failure to converge, and relatively large computational effort. We show that, by considering the interaction of successive characteristics reaching the interfaces, the stress-density relation can be deduced directly by recursion rather than iteration. This calculation is orders of magnitude faster than iterative analysis, and does not require an initial guess. Direct recursion may be less suitable for very noisy data, but it was robust when applied to trial data. The stress-density relation deduced was identical to the result from iterative Lagrangian analysis.
Diamond is used extensively as a component in high energy density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first principles theoretical EOS models, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Greuneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. The resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.
In this paper, a fast synthetic iterative scheme is developed to accelerate convergence for the implicit DOM based on the stationary phonon BTE. The key innovative point of the present scheme is the introduction of the macroscopic synthetic diffusion equation for the temperature, which is obtained from the zero- and first-order moment equations of the phonon BTE. The synthetic diffusion equation, which is asymptomatically preserving to the Fouriers heat conduction equation in the diffusive regime, contains a term related to the Fouriers law and a term determined by the second-order moment of the distribution function that reflects the non-Fourier heat transfer. The mesoscopic kinetic equation and macroscopic diffusion equations are tightly coupled together, because the diffusion equation provides the temperature for the BTE, while the BTE provides the high-order moment to the diffusion equation to describe the non-Fourier heat transfer. This synthetic iterative scheme strengthens the coupling of all phonons in the phase space to facilitate the fast convergence from the diffusive to ballistic regimes. Typical numerical tests in one-, two-, and three-dimensional problems demonstrate that our scheme can describe the multiscale heat transfer problems accurately and efficiently. For all test cases convergence is reached within one hundred iteration steps, which is one to three orders of magnitude faster than the traditional implicit DOM in the near-diffusive regime.
Tomography has made a radical impact on diverse fields ranging from the study of 3D atomic arrangements in matter to the study of human health in medicine. Despite its very diverse applications, the core of tomography remains the same, that is, a mathematical method must be implemented to reconstruct the 3D structure of an object from a number of 2D projections. In many scientific applications, however, the number of projections that can be measured is limited due to geometric constraints, tolerable radiation dose and/or acquisition speed. Thus it becomes an important problem to obtain the best-possible reconstruction from a limited number of projections. Here, we present the mathematical implementation of a tomographic algorithm, termed GENeralized Fourier Iterative REconstruction (GENFIRE). By iterating between real and reciprocal space, GENFIRE searches for a global solution that is concurrently consistent with the measured data and general physical constraints. The algorithm requires minimal human intervention and also incorporates angular refinement to reduce the tilt angle error. We demonstrate that GENFIRE can produce superior results relative to several other popular tomographic reconstruction techniques by numerical simulations, and by experimentally by reconstructing the 3D structure of a porous material and a frozen-hydrated marine cyanobacterium. Equipped with a graphical user interface, GENFIRE is freely available from our website and is expected to find broad applications across different disciplines.
We study the solution of the system of equations describing the dynamical evolution of spontaneous ruptures generated in a prestressed elastic-gravitational deforming body and governed by rate and state friction laws. We propose an iterative coupling scheme based on a weak formulation with nonlinear interior boundary conditions, both for continuous time and with implicit discretization (backward Euler) in time. We regularize the problem by introducing viscosity. This guarantees the convergence of the scheme for solutions of the regularized problems in both cases. We also make precise the conditions on the relevant coefficients for convergence to hold.