Engineering Papers⌕ Search

Engineering topics

Ma, Xia

Publications and source records attributed to Ma, Xia.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An investigation of the multi-mode Richtmyer-Meshkov instability at a gas/HE interface using Pagosa

Here, in this work, we present a hydrocode Pagosa and explore the Richtmyer-Meshkov Instability (RMI) at an air/high explosive (HE) interface for the first time that is important but has not received much attention yet in the high explosive safety field. Thus, the presented Pagosa can be expected to predict the whole deflagration-to-detonation transition (DDT) process in future. In Pagosa, spatial discretization is implemented on cubic staggered grids by computing different variables at the vertex and the cell center, respectively, a special operator-splitting technique is employed to reduce the computational cost, and an artificial viscosity is added to handle the discontinuous shock waves in our simulations. To quantitatively evaluate the capability of Pagosa to solve these kinds of instabilities, the single mode Rayleigh-Taylor instability (RTI) and the multimode RMI at an air/SF 6 interface are first performed, respectively. The Pagosa results are compared with the related numerical solutions in the existing references and the experimental result. Moreover, a theoretical derivation of growth of RTI is also provided based on our numerical method. Subsequently, we explore the multi-mode RMI at an air/HE interface as well as the effects of several factors using Pagosa. Numerical results show that Pagosa is a powerful toolset to generate the right structures and the amplitude of RTI and RMI at an air/SF 6 interface. The solid HE can be penetrated by a strong shock wave and forms RMI deformations. The RMI at an air/HE interface behaves very different than at an air/SF 6 interface, periodic, decreased oscillation is observed due to material character, and is very sensitive to the initial simulation settings, that is, a tiny change in physical quantities will lead to a remarkable RMI structure, which is also observed in a shock bubble interaction. The findings in this work are significant, and will present a new insight for the high explosive field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Simulation and analysis of the 3D double shock-induced recompression of spall using a hybrid PAGOSA with FLIP+MPM

In this work, an innovative hybrid algorithm PAGOSA with FLIP + MPM is first presented by coupling PAGOSA with the particle FLIP + MPM, and applied to systematically exploring the 3D recompression of spall induced by double-shock waves in ductile materials as the few existing studies are currently limited to the pure longitudinal recompression of spall, and ignore the 3D effects. This algorithm solves the difficulties the grid-based methods encounter when applied to capturing the fracture in material. We first validate the capabilities of PAGOSA with FLIP + MPM to predict different fracture/fragmentation cases by solving two benchmark problems and comparing with the analytical solutions or the experiment results. The convergences and the infinity-norm errors are also investigated. Subsequently, the 3D effects are demonstrated by simulating the spallation in material using PAGOSA with FLIP + MPM and comparing with those in the longitudinal cases. Additionally, some factors that affect the spallation with 3D effects are also discussed. Finally, the recompressions of spall in material with/without 3D effects are systematically simulated and analyzed. The recompression of spall driven by a high explosive is also simulated to show the ability of PAGOSA with FLIP + MPM to handle this kind of problem. The conditions that can lead to a recompression of spall are explored. Also, numerical results show that PAGOSA with FLIP + MPM can accurately predict different fracture cases, that the 3D effects play an important role in fracture/fragmentation, and the recompression of spall is very sensitive to the occurrence conditions, shows great differences when considering 3D effects in real applications. Moreover, the present hybrid method can be easily extended to other grid-based techniques employed for fracture in materials.

3D effects↗
DOE OSTI2021

The simulation of 3D hypervelocity spallation using a hydrocode PAGOSA with FLIP+MPM

n this work, a hydrocode PAGOSA with FLIP+MPM is presented and exercised to investigate the fracture in ductile material. The merit of PAGOSA with FLIP+MPM to solve the advection problem is first illustrated by a solid piston periodically moving in a sealed tube with air. Furthermore, the ability of PAGOSA with FLIP+MPM to capture the fracture in material is shown by a simple stretching fracture in ductile material. In both of two benchmark problems, the PAGOSA results and analytical solutions are also presented for comparison. Then PAGOSA with FLIP+MPM is used to model complex spall in ductile material, which is a challenging problem in engineering applications. The convergences of PAGOSA with FLIP+MPM—based on both the mesh size and the marker density—are investigated by monitoring free surface velocity. To further show the ability of PAGOSA with FLIP+MPM to predict fracture in ductile material, the numerical results are compared with the experimental results and other published numerical results. Moreover, the effect of the spall parameter in PAGOSA with FLIP+MPM on the numerical simulation is also analyzed by investigating the free surface velocity. Finally, the effects of the peak compressive stress, the tensile strain rate, and the loading rate on the spallation are further investigated. The numerical results show that PAGOSA with FLIP+MPM can improve PAGOSA's performance when applied to predicting fracture in ductile materials, and is robust enough to accurately predict spall fracture in ductile materials.

42 ENGINEERING↗
DOE OSTI2020

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗
DOE OSTI2020

Braiding Majorana fermions and creating quantum logic gates with vortices on a periodic pinning structure

We show how vortices that support Majorana fermions when placed on a periodic pinning array can be used for vortex exchange and independent braiding by performing a series of specific moves with a probe tip. Using these braiding operations, we demonstrate realizations of a Hadamard and a CNOT gate. We specifically consider the first matching field at which there is one vortex per pinning site, and we show that there are two basic dynamic operations, move and exchange, from which basic braiding operations can be constructed in order to create specific logic gates. The periodic pinning array permits both control of the world lines of the vortices and freedom for vortex manipulation using a set of specific moves of the probe during which the probe tip strength and height remain unchanged. We measure the robustness of the different moves against thermal effects and show that the three different operations produce distinct force signatures on the moving tip.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Shock Initiation Response of PBX 9502 Considering Rarefaction Wave Effects

The effect of rarefaction waves on the shock initiation processes of impacted explosives is investigated. The Scaled Uniform Reactive Front (SURF) explosives model is convenient for this investigation because mechanisms responsible for ignition, i.e., the physics of hot spots, is based on the shock pressure, while rarefactions are controlled by local pressure and treated separately in SURF. One -dimensional gas gun experiments are used to calibrate the ignition portion of SURF. Impact experiments conducted using thin disk and small diameter rods, where controlled one- and two-dimensional rarefactions are introduced, allow the local pressure-dependence of SURF to be parameterized.The Theory of Singular Surfaces was invoked to elucidate the thermo-mechanical and thermo-chemical behavior. Room and elevated (250°C) temperature conditions were explored. In all cases SURF could be parameterized to match the experimental results. In doing so, the role of rarefactions on the shock initiation process was gleaned.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2017

Relating microstructure, temperature, and chemistry to explosive ignition and shock sensitivity

An analysis is put forth that relates explosive properties including microstructure, temperature, and chemistry to explosive ignition and sensitivity. Unlike approaches that focus only on elucidating mesoscopic mechanisms important to ignition, the present analysis seeks a methodology directly applicable to continuum explosive burn models. Because the Scaled Uniform Reactive Front (SURF) burn model has a microstructural basis, it is chosen as the starting point of the present analysis. To build upon the SURF framework, a literal translation is performed of the high-level conceptual notions for which SURF is based, to concrete ignition–combustion parameters, a statistical description of the microstructure, and other thermo-chemical data acquired by non-shock experiments. The analysis requires a void volume fraction distribution acquired from ultra-small angle neutron scattering (USANS). The shock response of PBX 9502 is used to illustrate the theory. While the analysis requires calibration to a shock experiment (the pop plot, for example), the results are a self-consistent set of realistic physical quantities that contribute to the shock initiation process including, initial temperature, chemical kinetics, a statistical description of the microstructure, and the hot spot size, spacing, and temperature. The utility of having a mesoscale based theory is that once the critical hot spot temperature, as a function of shock pressure, is known for a given explosive and initial porosity distribution, the entire set of meso- and macro-scale results, including the pop plot can be calculated for any other porosity. Thus, one can understand changes in sensitivity at different densities (void size distributions), relying only on the assumption that the hot spot temperature curve would not differ significantly if the morphology of the hots spots were similar. Another important utility of the present analysis is to address the question of the role of initial temperature on the observed shock sensitivity of PBX 9502 and other explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

A separate phase drag model and a surrogate approximation for simulation of the steam assisted gravity drainage (SAGD) process

General ensemble phase averaged equations for multiphase flows have been specialized for the simulation of the steam assisted gravity drainage (SAGD) process. In the average momentum equation, fluid-solid and fluid-fluid viscous interactions are represented by separate force terms. This equation has a form similar to that of Darcy’s law for multiphase flow but augmented by the fluid-fluid viscous forces. Models for these fluid-fluid interactions are suggested and implemented into the numerical code CartaBlanca. Numerical results indicate that the model captures the main features of the multiphase flow in the SAGD process, but the detailed features, such as plumes are missed. We find that viscous coupling among the fluid phases is important. Advection time scales for the different fluids differ by several orders of magnitude because of vast viscosity differences. Numerically resolving all of these time scales is time consuming. To address this problem, we introduce a steam surrogate approximation to increase the steam advection time scale, while keeping the mass and energy fluxes well approximated. This approximation leads to about a 40-fold speed-up in execution speed of the numerical calculations at the cost of a few percent error in the relevant quantities.

02 PETROLEUM↗
DOE OSTI2015

NSR&D FY15 Final Report. Modeling Mechanical, Thermal, and Chemical Effects of Impact

The main goal of this project is to develop a computer model that explains and predicts coupled mechanical, thermal and chemical responses of HE under impact and friction insults. The modeling effort is based on the LANL-developed CartaBlanca code, which is implemented with the dual domain material point (DDMP) method to calculate complex and coupled thermal, chemical and mechanical effects among fluids, solids and the transitions between the states. In FY 15, we have implemented the TEPLA material model for metal and performed preliminary can penetration simulation and begun to link with experiment. Currently, we are working on implementing a shock to detonation transition (SDT) model (SURF) and JWL equation of state.

36 MATERIALS SCIENCE↗
DOE OSTI2014

NSR&D FY14 Final Report: Friction/Impact Modeling

The LANL-developed CartaBlanca code uses advanced techniques, including the Dual Domain Material Point (DDMP) method, to calculate fluid, solid motions and fluid-structure interactions. In the last year we have implemented the ViscoSCRAM material model in CartaBlanca based on Clements’ Abaqus-Explicit Finite Element implementation. In this fiscal year, we study numerical properties of this implementation. To consider coupled mechanical, thermal and chemical effects, we have implemented Henson’s HE decomposition model to couple with the thermal and mechanical packages in CartaBlanca. We have also developed porous media flow modules to consider heat convection effects of the reaction product gas inside pores of HE. The next immediate step is to perform numerical simulations to study interactions among these different mechanisms. We have started such calculation and obtained interesting initial results at the end of this fiscal year.

97 MATHEMATICS AND COMPUTING↗