SEARCH · Engineering Papers

Results for “shock capturing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Novel Mechanisms for the Formation of Dislocation Cell Patterns in BCC Metal

In this study, we present the first simulation results of the formation of dislocation cell wall microstructures in tantalum subjected to shock loading. Dislocation patterns and cell wall formation are important to understanding the mechanical properties of the materials in which they spontaneously arise, and yet the processing and self-assembly mechanisms leading to their formation are poorly understood. By employing transmission electron microscopy and discrete dislocation dynamics, we propose a new mechanism involving coplanar dislocations and pseudo-dipole mixed dislocation arrays that is essential to the pattern formation process. Our large-scale 3D DDD simulations demonstrate the self-organization of dislocation networks into cell walls in deformed BCC metal (tantalum) persisting at strain ε = 20%. The simulation analysis captures several crucial aspects of how the dislocation cell pattern affects metal plasticity, as observed in experiments. Although experimental evidence is inconclusive regarding whether cell wall formation takes place at the shock front, after the shock, during release, or when the sample has had enough time to relax post-recovery, our simulations indicate cell wall formation occurs after the shock and before release. The extended Taylor hardening composite model effectively considers the non-uniform dislocation density when cell walls form and accurately describes the corresponding flow stress.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗

DOE OSTI2023

Numerical simulations of the Los Alamos gapstick experiment

The gapstick is a high explosive (HE) sensitivity test recently developed at Los Alamos. The experiment design was motivated by the traditional gap test and consists of a series of HE and inert pellets in a rate stick configuration. The inert pellets are made increasingly longer so that eventually the attenuated shock is unable to initiate the next HE pellet. While the gapstick is a mechanically simple experiment it poses several challenges for numerical simulation. In particular, accurately modeling the HE initiation and detonation phenomena requires a reactive burn model with sufficient mesh resolution to capture the reaction scales. In this work, the Scaled Uniform Reactive Front (SURF) and Arrhenius Wescott-Stewart-Davis (AWSD) burn models are used for simulations of a PBX 9501 (95 wt% HMX, 5% binder) gapstick. Furthermore, the material model for the inert pellet material, Polyvinylidene Fluoride (PVDF), has a direct influence on the shock propagation and a new equation of state (EOS) calibration for PVDF is developed using available Hugoniot data. Despite some challenges, the simulations are able to reasonably predict shock transit velocities and detonation failure in the gapstick.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Strength of porous α-SiO 2 in a shock loaded environment: Calibration via Richtmyer–Meshkov instability and validation via Mach lens

The strength of brittle porous media is of concern in numerous applications, for example, earth penetration, crater formation, and blast loading. Thus, it is of importance to possess techniques that allow for constitutive model calibration within the laboratory setting. The goal of the current work is to demonstrate an experimental technique allowing for strength assessment of porous media subjected to shock loading, which can be implemented into pressure-dependent yield surfaces within numerical simulation schemes. As a case study, the deviatoric response of distended α-SiO 2 has been captured in a tamped Richtmyer–Meshkov instability (RMI) environment at a pressure regime of 4–10 GPa. Hydrocode simulations were used to interpret RMI experimental data, and a resulting pressure-dependent yield surface akin to the often employed modified Drucker–Prager model was calibrated. Simulations indicate that the resulting jet length generated by the RMI is sensitive to the porous media strength, thereby providing a feasible experimental platform capable of capturing the pressurized granular deviatoric response. Furthermore, in efforts to validate the RMI-calibrated strength model, a set of Mach-lens experiments was performed and simulated with the calibrated pressure-dependent yield surface. Excellent agreement between the resulting Mach-lens length in experiment and simulation provides additional confidence to the RMI yield-surface calibration scheme.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Experimentally validated multiphysics modeling of fracture induced by thermal shocks in sintered UO 2 pellets

Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗

DOE OSTI2026

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗

DOE OSTI2025

Magnetically Driven Neutron-rich Ejecta Unleashed: Global 3D Neutrino–General Relativistic Magnetohydrodynamic Simulations of Collapsars Probe the Conditions for r -process Nucleosynthesis

Abstract Collapsars—rapidly rotating stellar cores that form black holes—can power gamma-ray bursts and are proposed to be key contributors to the production of heavy elements in the Universe via the rapid neutron capture process ( r -process). Previous neutrino-transport collapsar simulations have been unable to unbind neutron-rich material from the disk. However, these simulations have not included sufficiently strong magnetic fields and the black hole (BH), both of which are essential for launching mass outflows. We present ν h-amr , a novel neutrino-transport general relativistic magnetohydrodynamic ( ν GRMHD) code, which we use to perform the first 3D global ν GRMHD collapsar simulations. We find a self-consistent formation of a weakly magnetized dense accretion disk, which has sufficient time to neutronize. Eventually, substantial magnetic flux accumulates near the BH, becomes dynamically important, leads to a magnetically arrested disk (MAD), and unbinds some of the neutron-rich material. However, the strong flux also hinders accretion, lowers density, and increases neutrino-cooling timescale, which prevents further disk neutronization. Typical collapsar progenitors with mass accretion rates, M ̇ ∼ 0.1 − 1 M ⊙ s - 1 , do not produce significant neutron-rich ( Y e < 0.25) ejecta. However, we find that MADs at higher mass accretion rates, M ̇ ≳ few M ⊙ s - 1 (e.g., for more centrally concentrated progenitors), can unbind M ej ≲ M ⊙ of neutron-rich ejecta. The outflows inflate a shocked cocoon that mixes with the infalling neutron-poor stellar gas and raises the final outflow Y e ; however, the final r -process yield may be determined earlier at the point of neutron capture freeze-out. Future work will explore under what conditions more typical collapsar engines become r -process factories.

Issa, Danat (ORCID:0009000524787631)↗

DOE OSTI2026

Improving Chirped Fiber Bragg Grating Resolution for Position-Sensitive Sensors in Shock- and Detonation-Driven Experiments

Chirped fiber Bragg gratings (CFBGs) are robust diagnostic sensors that are widely used to track detonation-driven and shock wave propagation. CFBGs are inscribed with a linearly chirped periodic index of refraction changes that alter the Bragg wavelength along the length of the probe. The light return of each individual Bragg element is captured by a detector at a unique time to map the full reflected spectrum. The CFBG spectrum is measured with a dispersive Fourier transform of the reflected light that temporally stretches the spectrum to increase spatial resolution and make a one-to-one map of the wavelength on a time axis. Here, we propose an improvement of CFBG temporal resolution by incorporating two co-linear laser pulses with orthogonal polarization states and a 5 ns time offset. The two separate signals were split and tracked by two separate detectors. An oscilloscope captured good separation in the signals, and two separate spectrograms were generated and interleaved in the post-processing of the data. This novel technique doubled the CFBG temporal resolution and led to a doubled location resolution. As a proof-of-concept of this technique, the resolution improvement was compared between standard CFBG measurements and the two polarization states method on a position-sensitive CFBG sensor. CFBG resolution doubling will advance sensor capabilities and will have a direct impact on improving capture and analysis in dynamic, high-explosive experiments.

42 ENGINEERING↗

DOE OSTI2023

Two-Dimensional Thomson Scattering in Laser-Produced Plasmas

We present two-dimensional (2D) optical Thomson scattering measurements of electron density and temperature in laser-produced plasmas. The novel instrument directly measures ne(x,y) and Te(x,y) in two dimensions over large spatial regions (cm2) with sub-mm spatial resolution, by automatically translating the scattering volume while the plasma is produced repeatedly by irradiating a solid target with a high-repetition-rate laser beam (10 J, ∼1012 W/cm2, 1 Hz). In this paper, we describe the design and motorized auto-alignment of the instrument and the computerized algorithm that autonomously fits the spectral distribution function to the tens-of-thousands of measured scattering spectra, and captures the transition from the collective to the non-collective regime with distance from the target. As an example, we present the first 2D scattering measurements in laser-driven shock waves in ambient nitrogen gas at a pressure of 0.13 mbar.

Zhang, Haiping (ORCID:0009000228133852)↗

DOE OSTI2022

Finding Fast Transients in Real Time Using a Novel Light-curve Analysis Algorithm

The current data acquisition rate of astronomical transient surveys and the promise for significantly higher rates in the next decade necessitate the development of novel approaches to analyze astronomical data sets and promptly detect objects of interest. The Deeper, Wider, Faster (DWF) program is a survey focused on the identification of fast-evolving transients, such as fast radio bursts, gamma-ray bursts, and supernova shock breakouts. It employs multifrequency simultaneous coverage of the same part of the sky over several orders of magnitude. Using the Dark Energy Camera mounted on the 4 m Blanco telescope, DWF captures a 20 s g -band exposure every minute, at a typical seeing of ~1'' and an air mass of ~1.5. These optical data are collected simultaneously with observations conducted over the entire electromagnetic spectrum—from radio to γ -rays—as well as cosmic-ray observations. In this paper, we present a novel real-time light-curve analysis algorithm, designed to detect transients in the DWF optical data; this algorithm functions independently from, or in conjunction with, image subtraction. We present a sample of fast transients detected by our algorithm, as well as a false-positive analysis. Our algorithm is customizable and can be tuned to be sensitive to transients evolving over different timescales and flux ranges.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Multidimensional Numerical Modeling of Combustion Dynamics in a Non-Premixed Rotating Detonation Engine With Adaptive Mesh Refinement

In the present work, a novel computational fluid dynamics (CFD) methodology was developed to simulate full-scale non-premixed rotating detonation engines (RDEs). A unique feature of the modeling approach was the incorporation of adaptive mesh refinement (AMR) to achieve a good trade-off between model accuracy and computational expense. Here, unsteady Reynolds-averaged Navier–Stokes (RANS) simulations were performed for an Air Force Research Laboratory (AFRL) non-premixed RDE configuration with hydrogen as fuel and air as the oxidizer. The finite-rate chemistry model, along with a ten-species detailed kinetic mechanism, was employed to describe the H 2 -Air combustion chemistry. Three distinct operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel/air mass flowrates. For all conditions, the capability of the model to capture essential detonation wave dynamics was assessed. An exhaustive verification and validation study was performed against experimental data in terms of a number of waves, wave frequency, wave height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model was demonstrated to accurately predict the sensitivity of these wave characteristics to the operating conditions, both qualitatively and quantitatively. A comprehensive heat release analysis was also conducted to quantify detonative versus deflagrative burning for the three simulated cases. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

42 ENGINEERING↗

DOE OSTI2023

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Deep Learning-based Parameterization of Complex 3D CO2 Saturation Data in Large-scale Geological Carbon Storage

In deep learning (DL), dimension reduction plays a pivotal role in improving training efficiency and minimizing overfitting, especially when working with complex datasets like three-dimensional (3D) saturation data. In the context of geological carbon storage (GCS), 3D saturation data introduces unique challenges due to its sparse nature and sharp transitions at plume boundaries, known as shock fronts. To tackle these challenges, we developed a novel DL framework that combines dimension reduction with advanced 3D reconstruction techniques. Our approach utilizes latent variables derived from 2D average saturation fields to efficiently capture the essential features of high-dimensional data while reducing the number of variables. This enhances both the robustness and accuracy of DL models, making the framework more practical for real-world applications. By offering a tailored solution for modeling complex 3D saturation dynamics, this framework holds significant potential for environmental monitoring, energy storage, and other geological applications.

Wang, Hongsheng [University of Texas at Austin]↗

DOE OSTI2021

On the numerical accuracy in finite-volume methods to accurately capture turbulence in compressible flows

The goal of the present article is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over the cell volume. Here, third-, fifth- and seventh-order WENO-Z schemes are investigated. On a problem with a smooth solution, the theoretical order of convergence rate for each method is retrieved, and changing the order of the reconstruction at cell faces does not impact the results, whereas for a shock-driven problem all the methods collapse to first-order. Here, study of the decay of compressible homogeneous isotropic turbulence reveals that using a high-order quadrature rule to compute the average over a finite-volume cell does not improve the spectral accuracy and that all methods present a second-order convergence rate. However the choice of the numerical method to reconstruct data at cell faces is found to be critical to correctly capture turbulent spectra. In the context of simulations with finite-volume methods of practical flows encountered in engineering applications, it becomes apparent that an efficient strategy is to perform the average integration with a low-order quadrature rule on a fine mesh resolution, whereas high-order schemes should be used to reconstruct data at cell faces.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗