No Arabic abstract
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.
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.
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.
A general formulation was developed to represent material models for applications in dynamic loading. Numerical methods were devised to calculate response to shock and ramp compression, and ramp decompression, generalizing previous solutions for scal
Most electronic properties of metals are determined solely by the low-energy states around the Fermi level, and for topological metals/semimetals, these low-energy states become distinct because of their unusual energy dispersion and emergent pseudospin degree of freedom. Here, we propose a class of materials which are termed as quadratic contact point (QCP) semimetals. In these materials, the conduction and valence bands contact at isolated points in the Brillouin zone, around which the band dispersions are quadratic along all three directions. We show that in the absence/presence of spin-orbit coupling, there may exist triply/quadruply-degenerate QCPs that are protected by the crystalline symmetry. We construct effective models to characterize the low-energy fermions near these QCPs. Under strong magnetic field, unlike the usual 3D electron gas, there appear unconventional features in the Landau spectrum. The QCP semimetal phase is adjacent to a variety of topological phases. For example, by breaking symmetries via Zeeman field or lattice strain, it can be transformed into a Weyl semimetal with Weyl and double Weyl points, a Z2 topological insulator/metal, or a Dirac semimetal. Via first-principles calculations, we identify realistic materials Cu2Se and RhAs3 as candidates for QCP semimetals.
Mixtures of fluids and granular sediments play an important role in many industrial, geotechnical, and aerospace engineering problems, from waste management and transportation (liquid--sediment mixtures) to dust kick-up below helicopter rotors (gas--sediment mixtures). These mixed flows often involve bulk motion of hundreds of billions of individual sediment particles and can contain both highly turbulent regions and static, non-flowing regions. This breadth of phenomena necessitates the use of continuum simulation methods, such as the material point method (MPM), which can accurately capture these large deformations while also tracking the Lagrangian features of the flow (e.g. the granular surface, elastic stress, etc.). Recent works using two-phase MPM frameworks to simulate these mixtures have shown substantial promise; however, these approaches are hindered by the numerical limitations of MPM when simulating pure fluids. In addition to the well-known particle ringing instability and difficulty defining inflow/outflow boundary conditions, MPM has a tendency to accumulate quadrature errors as materials deform, increasing the rate of overall error growth as simulations progress. In this work, we present an improved, two-phase continuum simulation framework that uses the finite volume method (FVM) to solve the fluid phase equations of motion and MPM to solve the solid phase equations of motion, substantially reducing the effect of these errors and providing better accuracy and stability for long-duration simulations of these mixtures.