SEARCH · Engineering Papers

Results for “SHOCK”

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 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

NASA NTRS2014

High Order Numerical Methods for LES of Turbulent Flows with Shocks

Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten's subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.

High Order↗

NASA NTRS2016

Relativistic Electrons Produced by Foreshock Disturbances Observed Upstream of Earth's Bow Shock

Charged particles can be reflected and accelerated by strong (i.e., high Mach number) astrophysical collisionless shock waves, streaming away to form a foreshock region in communication with the shock. Foreshocks are primarily populated by suprathermal ions that can generate foreshock disturbances-largescale (i.e., tens to thousands of thermal ion Larmor radii), transient (approximately 5-10 per day) structures. They have recently been found to accelerate ions to energies of several keV. Although electrons in Saturn's high Mach number (M > 40) bow shock can be accelerated to relativistic energies (nearly 1000 keV), it has hitherto been thought impossible to accelerate electrons beyond a few tens of keV at Earth's low Mach number (1 ≤M <20) bow shock. Here we report observations of electrons energized by foreshock disturbances to energies up to at least approximately 300 keV. Although such energetic electrons have been previously observed, their presence has been attributed to escaping magnetospheric particles or solar events. These relativistic electrons are not associated with any solar or magnetospheric activity. Further, due to their relatively small Larmor radii (compared to magnetic gradient scale lengths) and large thermal speeds (compared to shock speeds), no known shock acceleration mechanism can energize thermal electrons up to relativistic energies. The discovery of relativistic electrons associated with foreshock structures commonly generated in astrophysical shocks could provide a new paradigm for electron injections and acceleration in collisionless plasmas.

SLAMS↗

NASA NTRS2015

High Order Numerical Methods for LES of Turbulent Flows with Shocks

Simulation of turbulent flows with shocks employing explicit subgrid-scale (SGS) filtering may encounter a loss of accuracy in the vicinity of a shock. In this work we perform a comparative study of different approaches to reduce this loss of accuracy within the framework of the dynamic Germano SGS model. One of the possible approaches is to apply Harten’s subcell resolution procedure to locate and sharpen the shock, and to use a one-sided test filter at the grid points adjacent to the exact shock location. The other considered approach is local disabling of the SGS terms in the vicinity of the shock location. In this study we use a canonical shock-turbulence interaction problem for comparison of the considered modifications of the SGS filtering procedure. For the considered test case both approaches show a similar improvement in the accuracy near the shock.

LES↗

NASA NTRS2017

Solar Flare Termination Shock and the Synthetic Fe XXI 1354.08 Å Line

Solar flares are one of the most energetic phenomena occurred in the solar system. In the standard solar flare model, a fast mode shock, which is often referred to as the flare termination shock (TS), can exist above the loop-top source of hard X-ray emissions. The existence of the termination shock has been recently related to spectral hardening of flare hard X-ray spectrum at energies > 300 keV. Observations of the Fe XXI 1354.08 Å line during solar flares by the IRIS spacecraft have found significant redshift with >100 km/s, which is consistent with a reconnection downflow. The ability to identify such a redshift by IRIS is made possible by IRIS's high time resolution, high spatial resolution, high sensitivity and cadence spectral observations. The ability to identify such a redshift by IRIS suggests that one may be able to use IRIS observations to identify flare termination shocks. Using a MHD simulation to model magnetic reconnection of a solar flare and assuming the existence of a TS in the downflow of the reconnection plasma, we model the synthetic emission of the Fe XXI 1354.08 Å line in this work. We show that the existence of the TS in the solar flare may manifest itself from the Fe XXI 1354.08 Å line.Plain Language SummaryA solar flare is a brief eruption of intense high-energy radiation from the sun's surface, associated with sunspots and causing electromagnetic disturbances on the earth, as with radio frequency communications and power line transmissions. The standard solar flare model proposes that magnetic reconnection drives solar flares and fast mode shock, which is often referred to as the flare termination shock (TS), can exist at the reconnection outflows. The existence of the termination shock has been recently related to spectral hardening of flare hard X-ray spectrum.

Guo, Lijia↗

NASA NTRS2019

Time Accurate Simulation of Nonequilibrium Flow Inside the NASA Ames Electric Arc Shock Tube

Numerical investigations of the flowfield inside NASA Ames' Electric Arc Shock Tube have been performed. The focus is to simulate the experiments designed to reproduce shock layer radiation layer relevant to Earth re-entry conditions. This paper assess the current computational capability in simulating time-accurate unsteady nonequilibrium flows in the presence of strong shock waves with state-of-the-art physical models. The technical approach is described with preliminary results presented for one specific flow condition. It was found that the axisymmetric source term generates a numerical instability that appears as shock bending. This instability is time dependent which greatly affects the shock speed. Post-shock conditions are discussed and compared to CEA equilibrium prediction and good agreement was obtained close to the test-section and just behind the shock.

Bensassi, Khalil↗

NASA NTRSDate not supplied

Deciphering Extreme Mineral Records; Microstructural Phase Heritage of Shocked Materials

High-pressure minerals such as coesite, stishovite, and reidite are often used as indicators of shock metamorphism. However, high-pressure and high-temperature phases are metastable at ambient conditions and, therefore, often revert to a more energetically favorable phase. During solid-state transformations, the crystallographic orientations of the stable polymorph are controlled by the transformation pathway from the parent. Microstructural orientation analysis using electron backscatter diffraction (EBSD) of the stable phase can reveal systematic intercrystalline orientation relationships (OR) diagnostic of solid-state transformations from these high-pressure or high-temperature polymorph phases (Cayron et al., 2006). Knowledge of polymorph stability fields for a given system can thus be used to infer minimum shock P and/or T conditions that rocks have experienced, which are otherwise unavailable using traditional thermobarometers such as element partition coefficients. Using orientation relationships to infer the former presence of mineral phases has been termed microstructural “phase heritage”, and has successfully been used to probe for evidence of extreme pressures and temperatures resulting from hypervelocity impacts. The reconstruction of polymorphs of ZrSiO4 (zircon, reidite) and ZrO2 (baddeleyite) have further elucidated the behavior of these phases under hypervelocity impact conditions (Timms et al., 2017a; White et al., 2018), and OR analysis of shocked monazite has identified a previously unrecognized tetragonal high-P polymorph (Erickson et al., 2019). Zircon will convert to the low tetragonal polymorph reidite at shock stresses above ~21 GPa. However, upon decompression at temperatures above 1200 °C reidite will revert to zircon (Kusaba et al., 1985). Recrystallized zircon neoblasts encased in impact melt often show systematic misorientation relationships of 90° about <110> indicative of reversion from the reidite (Fig. 1), and preserving cryptic evidence of the high-P history that is often erased by the post shock thermal spike (Cavosie et al., 2016). Microstructural OR analysis of a dissociated zircon corona composed of baddeleyite (monoclinic-ZrO2) has identified the former presence of cubic zirconia within impact melt glass from the Mistastin Lake impact structure, Labrador, CA. Based on stability fields in the ZrO2-SiO2 binary phase diagram, cubic zirconia records thermal conditions in excess of 2370 °C at ambient P, indicating that the Mistastin impact melt is the hottest rock identified on Earth’s surface (Timms et al., 2017b). Shock-deformed monazite (monoclinic La,Ce,ThPO4) containing lamellae comprised of interlocking laths has been identified from the Haughton Dome impact structure, Nunavut, Canada and Nördlinger-Ries Crater, southern Germany. The lath-structured lamellae are composed of four systematic orientation variants, which OR analyses indicate originate from a tetragonal parent whereby the [010]monoclinic aligns with either [100] or [010]tetragonal (Erickson et al., 2019). These results highlight utility of microstructural phase heritage analyses for the identification of unstable high-P or high-T polymorphs that uniquely record the extremely transient shock conditions produced by hypervelocity impacts.

Timmons Erickson↗

NASA NTRS2021

MMS Observations of Energized He+ Pickup Ions at Quasiperpendicular Shocks

Accelerated He+ PUIs, downstream of quasi-perpendicular shocks, are studied as 10a function of the fast mode Mach number (Mf) and shock obliquity (𝜃Bn).We analyze quasi-perpendicular shocks with Mach numbers in the range [1, 7] observed by the Magnetospheric MultiScale (MMS) mission, and compare upstream and downstream He+ velocity distribution functions (VDFs). For each shock event, we characterize the upstream PUI distribution and derive reduced 1D velocity distributions for the selected upstream and downstream intervals. We also compare the upstream to downstream ratio of spectral indices, computed from the He+ perpendicular distributions, to Mf and 𝜃Bn. We find a positive correlation between this spectral index ratio and Mf, which suggests that perpendicular energization of He+ PUIs is enhanced as the shock becomes stronger. These results inform modeling efforts of PUIs and shock acceleration processes, particularly those taking place at the termination shock.

Michael J. Starkey↗

NASA NTRSDate not supplied

MARS 2020 Backshell Radiative Heating Measurement and Shock Tube Verification

In February 2021, the Perseverance rover was brought to the surface of Mars by the Mars 2020 mission. A feature of the Mars 2020 capsule was instrumentation to measure its entry, descent and landing (EDL) with the so-called Mars EDL Instrumentation 2 (MEDLI2) [1]. The MEDLI2 introduced, among other things, backshell instrumentation, including a broadband radiometer. The radiometer was mounted on the leeside of the vehicle next to thermocouple plugs and a heat flux gauge. The data returned by the leeside MEDLI2 heat flux gauge is largely analogous to that measured by the COMARS gauge flown on the ExoMars Schiaparelli [2] entry in that it is measuring heat flux in an area that is entirely dominated by radiative heating. The COMARS measurement provided excellent validation of backshell radiative heating models for Mars entry, albeit at a limited number of points [3]. The MEDLI2 heat flux gauge measurement effectively confirmed the quality of the prediction, extended over the full trajectory, as will be presented in this paper. The MEDLI2 radiometer,however, was blocked by ablation products and suffered a loss of half of its signal. The second backshell heat flux gauge installed on the windside of the vehicle was also well predicted, although the heating had both radiative and convective contributions. It was desired to reproduce the conditions of the Mars 2020 entry via ground testing in the Electric Arc Shock Tube (EAST) at NASA Ames. Tests to verify stagnation line heating were previously reported in EAST, confirming the presence of shock layer radiation as the major discrepancy in heatshield temperature modelling [4]. Therefore, tests for stagnation line heating were not repeated. Instead, the shock tube informed bias method [5, 6] was used to identify test conditions that may produce similarity to streamlines that pass around the backside of the vehicle and are responsible for the radiation observed at the two heat flux gauge locations. This method was used to identify a range of velocities and densities in the shock tube that are relevant for confirming the radiative environment encountered. This paper reports the results obtained in the 10 cm diameter EAST shock tube, corresponding to later trajectory points at ambient pressures of 1.1-2.0 Torr and velocities from 1.2-3.5 km/s. The test series employed two primary diagnostics: emission spectroscopy and tunable diode laser absorption spectroscopy (TDLAS). The emission spectroscopy performed broadband measurements of the radiative emission of the 4.3 and 2.7 m bands of CO2 at flight similar conditions, obtaining both spectral and spatial data corresponding to the relaxation behind the shock front. The TDLAS measured the absorption of several lines of CO and CO2 and obtains species number densities and temperatures as a function of time behind the shock front. This paper will review highlights of this test series and analyses of the emission and absorption data. While the datasets generally show good agreement with predictions, a few discrepancies and items for additional investigation are identified and will be discussed

Brett A Cruden↗

DOE OSTI2019

Ultra-High Fidelity Laser-Induced Air Shock from Energetic Materials

Recent interest in the implementation of the Laser-induced Air Shock from Energetic Materials (LASEM) technique prompted an investigation using ultra-high-speed imaging diagnostics to provide early-time (sub-μs) shock-wave-radius-versus-time data necessary to accurately determine the characteristic laser-induced-shock velocity (y-intercept) for high explosives (HE) of interest. Ultra-high-speed focused-shadowgraphy images were collected from nanosecond-pulsed laser-ablated samples of HE similar to those in the published literature. Shadowgraphy images were collected using interframe times ranging from 50–750 ns, with exposure times of 5 ns. Acquired shock-wave-radius-versus-time data permitted a high-fidelity assessment of the shock wave velocity produced at the characteristic radius of the generated laser-plasma pulse. The resulting data from eight different HE indicated that early-time energy contributions and the resulting laser-induced shock velocity profiles produced from ablation of the HE material were indistinguishable amongst the explosives tested for the majority of the time domain characterized (0–12 μs).

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗

DOE OSTI2020

High-precision shock equation of state measurements for metallic fluid carbon between 15 and 20 Mbar

Diamond is an efficient ablator material to convert the energy of high-power giant lasers into ablation pressure with applications for High-Energy-Density (HED) science, planetary science, and Inertial Confinement Fusion (ICF) research at the National Ignition Facility (NIF). Unfortunately, current theoretical equation of state models cannot reproduce all the observed experimental data in the multi-megabar regime particularly relevant for HED and ICF research. New experimental data on the behavior of carbon at extreme pressures and temperatures are, therefore, essential to improve our predictive capability to design and analyze dynamic compression experiments for HED or ICF research and build improved equation of state models in the future. In this paper, we report high-precision laser-driven shock compression measurements on diamond single crystals at the Omega Laser Facility. Using ultrafast Doppler optical Velocimetry Interferometer System for Any Reflector (VISAR) to track the leading shock front and a quartz plate as an in situ reference, we obtain relative pressure-density shock equation-of-state measurements between 15 and 20 Mbar with an impedance-matching procedure. We also report shock-and-release measurements in a spherical geometry at the NIF. The new data provide tight constraints on the compressibility of warm dense carbon along the Hugoniot of full density diamond, allowing us to discriminate between existing theoretical equation-of-state models. We find that both LLNL LEOS 9061 and LANL Sesame 7835 models capture well the shock compressibility in the explored range. LANL Sesame 7835 also reproduces well the observed shock-and-release behavior of diamond near 10–20 Mbar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Experimental observations of shock-wave-induced bubble collapse and hot-spot formation in nitromethane liquid explosive

We studied the collapse of individual helium gas bubbles in the homogeneous explosive nitromethane (NM) to investigate effects of hot-spot formation on the detonation process. A bubble was injected into a NM sample, and a shock wave from an explosive detonator compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and the pre- and post-bubble collapse processes. An infrared camera image showed the residual radiance temperature distribution after the bubble collapse, and an optical fiber pyrometer measured the time-resolved thermal radiance. We measured the optical spectra of light emitted from detonating NM without a bubble and from a collapsing bubble in shocked, undetonated NM. We estimated temperatures of the detonation fronts and of the hot spots formed by bubble collapse. To study the incipient detonation process, we performed all bubble collapse experiments at pressures below the threshold for creating a sustained detonation. Where the bubble collapsed, we observed an opaque, thermally emissive region believed to be chemical reaction products. Chemical reactions in NM can be produced with lower shock pressures (~1 GPa) when a helium bubble is present than without a bubble (~10 GPa). We used hydrodynamic modeling to predict shock wave propagation, extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. Simulations using a temperature-dependent Arrhenius burn model gave much better results than reactive burn models that depend only on pressure and density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

The role of pre-existing heterogeneities in materials under shock and spall

There has been a challenge for many decades to understand how heterogeneities influence the behavior of materials under shock loading, eventually leading to spall formation and failure. Experimental, analytical, and computational techniques have matured to the point where systematic studies of materials with complex microstructures under shock loading and the associated failure mechanisms are feasible. This is enabled by more accurate diagnostics as well as characterization methods. As interest in complex materials grows, understanding and predicting the role of heterogeneities in determining the dynamic behavior becomes crucial. Early computational studies, hydrocodes, in particular, historically preclude any irregularities in the form of defects and impurities in the material microstructure for the sake of simplification and to retain the hydrodynamic conservation equations. Contemporary computational methods, notably molecular dynamics simulations, can overcome this limitation by incorporating inhomogeneities albeit at a much lower length and time scale. This review discusses literature that has focused on investigating the role of various imperfections in the shock and spall behavior, emphasizing mainly heterogeneities such as second-phase particles, inclusions, and voids under both shock compression and release. Pre-existing defects are found in most engineering materials, ranging from thermodynamically necessary vacancies, to interstitial and dislocation, to microstructural features such as inclusions, second phase particles, voids, grain boundaries, and triple junctions. This literature review explores the interaction of these heterogeneities under shock loading during compression and release. Systematic characterization of material heterogeneities before and after shock loading, along with direct measurements of Hugoniot elastic limit and spall strength, allows for more generalized theories to be formulated. Further, continuous improvement toward time-resolved, in situ experimental data strengthens the ability to elucidate upon results gathered from simulations and analytical models, thus improving the overall ability to understand and predict how materials behave under dynamic loading.

36 MATERIALS SCIENCE↗

DOE OSTI2013

SHOCK HEATING OF THE MERGING GALAXY CLUSTER A521

A521 is an interacting galaxy cluster located at z = 0.247, hosting a low-frequency radio halo connected to an eastern radio relic. Previous Chandra observations hinted at the presence of an X-ray brightness edge at the position of the relic, which may be a shock front. We analyze a deep observation of A521 recently performed with XMM-Newton in order to probe the cluster structure up to the outermost regions covered by the radio emission. The cluster atmosphere exhibits various brightness and temperature anisotropies. In particular, two cluster cores appear to be separated by two cold fronts. We find two shock fronts, one that was suggested by Chandra and that is propagating to the east, and another to the southwestern cluster outskirt. The two main interacting clusters appear to be separated by a shock-heated region, which exhibits a spatial correlation with the radio halo. The outer edge of the radio relic coincides spatially with a shock front, suggesting that this shock is responsible for the generation of cosmic-ray electrons in the relic. The propagation direction and Mach number of the shock front derived from the gas density jump, M = 2.4 {+-} 0.2, are consistent with expectations from the radio spectral index, under the assumption of Fermi I acceleration mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRSDate not supplied

Application of Flux Conserved Modeling to an Unsteady Combustion Driven Pseudo-Shock

Quasi-steady correlation-based and flux conserved models of the pseudo-shock pressure profile are experimentally assessed in the Air Force Research Lab RC-18 Sub-Scale DirectConnect Isolator Combustion Facility. A bifurcated normal pseudo-shock is established in the 2.67 aspect ratio isolator test section by a downstream cavity-stabilized combustor. The Waltrup-Billig correlation, Smart flux conserved model, and Fiévet data driven flux conserved models are compared against measured pressure profiles. The literature models are found to perform poorly in a high-amplitude forced environment of RC-18, where shock speeds upwards of 20% the isolator Mach number occur. Two modifications to the quasi-steady Fiévet data driven flux conserved model are proposed including incorporating an estimate of pseudo-shock static pressure recovery and exercising the model in a shock stationary reference frame. It is shown that these modifications allow the quasi-steady flux conserved modeling approach to accurately predict pseudo-shock pressure profiles over the discrete phases of high-amplitude forced dynamics.

Shock Train↗

DOE OSTI2022

Shock physics in warm dense matter: A quantum hydrodynamics perspective

Abstract Warm dense matter (WDM), an exotic, highly compressed state of matter between solid and plasma phases, is of high current interest, in particular for astrophysics and inertial confinement fusion. For the latter, in particular the propagation of compression shocks is crucial. The main unknown in the shock propagation in WDM is the behaviour of the electrons since they are governed by correlations, quantum and spin effects that need to be accounted for simultaneously. Here we describe the shock dynamics of the warm dense electron gas using a quantum hydrodynamic model. From the numerical hydrodynamic simulations, we observe that the quantum Bohm pressure induces shear force that weakens the formation and strength of the shock. In addition, the Bohm pressure induces an electron density response that takes the form of oscillations. This is confirmed by the theoretical analysis of the early stage of the shock formation. Our theoretical and numerical analyses allow us to identify characteristic dimensionless shock propagation parameters at which the effect of the Bohm force is important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Off-center blast in a shocked medium

When multiple blasts occur at different times, the situation arises in which a blast wave is propagating into a medium that has already been shocked. Determining the evolution in shape of the second shock is not trivial, as it is propagating into air that is not only non-uniform, but also non-stationary. To accomplish this task, we employ the method of Kompaneets to determine the shape of a shock in a non-uniform media. We also draw from the work of Korycansky [1] on an off-center explosion in a medium with radially varying density. Extending this to treat non-stationary flow, and making use of approximations to the Sedov solution for the point blast problem, we are able to determine an analytic expression for the evolving shape of the second shock. Specifically, we consider the case of a shock in air at standard ambient temperature and pressure, with the second shock occurring shortly after the original blast wave reaches it, as in a sympathetic detonation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

A MUSCL-SCNI approach for meshfree modeling of shock waves in fluids

Here, a stable and nodally integrated meshfree formulation for modeling shock waves in fluids is developed. The reproducing kernel approximation is employed to discretize the conservation equations for compressible flow, and a flux vector splitting approach is applied to allow proper numerical treatments for the advection and pressure parts, respectively, based on the characteristics of each flux term. To capture the essential shock physics in fluids, including the Rankine–Hugoniot jump conditions and the entropy condition, local Riemann enrichment is introduced under the stabilized conforming nodal integration (SCNI) framework. Meanwhile, numerical instabilities associated with the advection flux are eliminated by adopting a modified upwind scheme. To further enhance accuracy, a MUSCL-type method is introduced in conjunction with an oscillation limiter to avoid Gibbs phenomenon and ensure monotonic piecewise linear reconstruction in the smooth region. The present meshfree formulation is free from tunable artificial parameters and is capable of capturing shock and rarefaction waves without over/undershoots. Finally, several numerical examples are analyzed to demonstrate the effectiveness of the proposed MUSCL-SCNI approach in meshfree modeling of complex shock phenomena, including shock diffraction, shock–vortex interaction, and high energy explosion processes.

42 ENGINEERING↗