Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHOCK TESTING MACHINE”

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.

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

At least 37 records · Page 2

arcjetCV: automating recession extraction from video

Arc jet Computer Vision (arcjetCV)[1][2] is a software application built to automate time-resolved recession tracking of edges in test videos, specifically for tracking material recession and the shock-material standoff. This provides a new capability to resolve and validate new physics associated with non-linear processes and an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcticCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.

Recession tracking↗

ArcjetCV: Automating Recession Tracking

Arc jet Computer Vision (arcjetCV) is a software application built to automate time-resolved recession tracking of edges in test videos, specifically for tracking material recession and the shock-material standoff. This provides a new capability to resolve and validate new physics associated with non-linear processes and an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcticCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.

Recession tracking↗

ArcjetCV: Automating Arc Jet Analysis

Arc jet Computer Vision (arcjetCV) is a software application built to automate time-resolved recession tracking of edges in test videos, specifically for tracking material recession and the shock-material standoff. This provides a new capability to resolve and validate new physics associated with non-linear processes and an essential step to reduce testing, modeling, and validation uncertainties for heatshield material performance. ArcticCV uses several types of machine learning (convolutional neural net: CNN, decision tree: DT, k-means unsupervised clustering: KM) to automate the video processing pipeline. These include inferring the start/stop of time segments of interest (1D CNN), measuring the time-dependent 2D recession of the material samples (2D CNN, DT), measuring the time-dependent shock standoff distance (2D CNN, DT), and post-processing cleaning of the recession data (KM). The software also provides a graphical user interface for ease of use. The results of using this tool on arc jet videos show non-linear time-dependent effects can be important for certain materials and characterizing certain failure modes.

Recession tracking↗

Durability of YSZ Coated Ti2AlC in 1300 °C Mach 0.3 Burner Rig Tests

A thermal barrier coating system survived burner rig testing at 1300 °C for 500 h. A 160 µm thick yttria stabilized zirconia (YSZ) coating was applied to a Ti2AlC MAX phase bar sample by plasma spray physical vapor deposition (PS-PVD) and tested face-on in an atmospheric Mach 0.3 jet fuel burner, using 5-h thermal cycles. No thermal barrier coating (TBC) spallation or recession was observed, only a 2.4 mg/cm2 mass gain. The modest weight gain precluded severe volatility losses under high velocity burner conditions. The coating surface exhibited colonies of (111)flourite fiber-textured columns separated by craze patterns, with no visible moisture attack. The metastable tetragonal t' YSZ phase was obtained initially, transitioning to equilibrium teq and cubic YSZ, but with little detrimental monoclinic. The thickness of the alumina TGO was ~21 to 23 μm under the heated YSZ face and ~13 to 15 μm on the uncoated, cooler backside. The backside exhibited removal of initial transient TiO2 nodules and partial etching of the underlying Al2O3 scale by volatile hydroxides formed in high temperature, high velocity water vapor. Aerodynamic forces produced some bending of the cantilevered sample via creep. The test indicated exceptional stability of YSZ coatings on Ti2AlC under turbine conditions, with thermal expansion matching playing a key role. The purpose of this study was to demonstrate long term durability of YSZ/MAX phase system in aggressive high temperature burner rig testing. MAX phases have been keenly studied because of their unique crystal structure and intriguing properties (Refs. 1 and 2). Having Mn+1(Al,Si)(C,N)n general composition, they are defined as ceramics, but possess unusual desirable attributes such as high conductivity, thermal shock resistance, easy machinability, and deformation tolerance. The mechanical properties derive from weak M-(Al,Si) bonding in the basal plane that leads to sliding and kinking in preference to catastrophic crack growth. Like most ceramics they are phase stable at high temperatures, generally up to 1500 °C. High temperature oxidation resistance is excellent for alumina-forming Ti3AlC2, Ti2AlC, and Cr2AlC, as reviewed by Tallman, et al. (Ref. 3). Compatibility with α-Al2O3 scales is further enhanced in cyclic exposures by a close matching of thermal expansion coefficients, (Ref. 4) i.e., (~9.3, 10.2, 11.3×10–6/K for Al2O3, Ti2AlC, and YSZ, to be discussed). Turbine environments generally contain 10 percent water vapor in the combustion gases, therefore moisture effects can be a concern for some materials (Ref. 5). Furnace tests of MAX phases in high temperature steam generally showed little effect on Al2O3 scale growth (Ref. 6). However, high velocity and high pressure gas can influence scale losses by the formation of volatile reaction products, such as TiO(OH)2 and Al(OH)3 (Refs. 7 to 10). This phenomenon had been discussed for 1100 to 1300 °C high pressure burner rig tests of Ti2AlC (Ref. 11). A single cubic growth rate parameter kcubic was measurably lower than comparable furnace TGA data, but it could be matched reasonably well if corrected for a slight volatility term. In general, a two-parameter cubic-linear growth-volatility law was believed to apply. Corresponding scale volatility loss rates, directly measured at 1300 °C on a pre-oxidized sample, were moderate (0.012 mg/cm2/h) and largely attributed to removal of the initial TiO2 transient scale. A related CH4 burner study of high purity Cr2AlC MAX phase demonstrated 1200 °C durability after 500 rapid (5 min. heat and 2 min. cool) thermal shock cycling (29 h hot time) (Ref. 12). Heating and cooling rates were ~1000 and 500 °C per minute, with a gas velocity of 5 m/s, producing a 75 °C/mm gradient. A 7 μm Al2O3 surface scale and a 13 μm Cr7C3 depletion zone formed with no signs of failure. No evidence of scale volatility was evident, although weight change was not provided, the velocity was moderate, and the total hot time was not extensive. The same high gradient BRT was used to produce 1400 °C surface temperatures for a YSZ/Cr2AlC/IN738 system in the first study of MAX phases used as bond coats for thermal barrier coatings (TBC) (Ref. 25). Here TBC failure was reported after 745 cycles, with only a 1.5 μm Al2O3 scale entrained within a porous, Cr7C3 bondcoat depletion phase. YSZ thermal barrier coatings have been considered to be a compatible complement to Al-MAX phases because of thermal expansion matching and extremely low volatility in water vapor. Initial studies showed superior oxidative stability up to 1300 °C, for long times (at least 500 h) for Ti2AlC substrates and less (268 h) for Cr2AlC, while withstanding large alumina TGO scale thickness (~35 to 40 μm) (Refs. 13 and 14). By comparison, typical superalloy systems can only survive 1150 °C maximum interface temperatures for extended periods, with a maximum sustained TGO below 10 μm (Ref. 15). High temperature SiC based systems are known to form slow-growing SiO2 scales. But these are subject to rate enhancement and volatile Si(OH)4 products in the presence of water vapor, as described comprehensively by Opila, et al. (Refs. 5, 16 to 19). Net weight losses are generally observed in high velocity, high pressure burner rig studies (e.g., 0.084 mg/cm2/h at 1300 °C) (Ref. 20). Furthermore, the loss rates have been shown from chemical physics to scale with v1/2 and pH2O 2 (Ref. 16). Low activity, moisture-resistant environmental barrier coatings (EBC), such as rare earth silicates, are needed to prevent substrate recession under turbine conditions (Refs. 21 to 23).

Smialek, James L.↗

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗

Experimental studies on the tripping behavior of narrow T-stiffened flat plates subjected to hydrostatic pressure and underwater shock

An experimental investigation was conducted to determine the static and dynamic responses of a specific stiffened flat plate design. The air-backed rectangular flat plates of 6061-T6 aluminum with an externally machined longitudinal narrow-flanged T-stiffener and clamped boundary conditions were subjected to static loading by water hydropump pressure and shock loading from an eight pound TNT charge detonated underwater. The dynamic test plate was instrumented to measure transient strains and free field pressure. The static test plate was instrumented to measure transient strains, plate deflection, and pressure. Emphasis was placed upon forcing static and dynamic stiffener tripping, obtaining relevant strain and pressure data, and studying the associated plate-stiffener behavior.

Budweg, H. L.↗

An Explainable Machine-Learning Model for Compensatory Reserve Measurement: Methods for Feature Selection and the Effects of Subject Variability

Tracking vital signs accurately is critical for triaging a patient and ensuring timely therapeutic intervention. The patient’s status is often clouded by compensatory mechanisms that can mask injury severity. The compensatory reserve measurement (CRM) is a triaging tool derived from an arterial waveform that has been shown to allow for earlier detection of hemorrhagic shock. However, the deep-learning artificial neural networks developed for its estimation do not explain how specific arterial waveform elements lead to predicting CRM due to the large number of parameters needed to tune these models. Alternatively, we investigate how classical machine-learning models driven by specific features extracted from the arterial waveform can be used to estimate CRM. More than 50 features were extracted from human arterial blood pressure data sets collected during simulated hypovolemic shock resulting from exposure to progressive levels of lower body negative pressure. A bagged decision tree design using the ten most significant features was selected as optimal for CRM estimation. This resulted in an average root mean squared error in all test data of 0.171, similar to the error for a deep-learning CRM algorithm at 0.159. By separating the dataset into sub-groups based on the severity of simulated hypovolemic shock withstood, large subject variability was observed, and the key features identified for these sub-groups differed. This methodology could allow for the identification of unique features and machine-learning models to differentiate individuals with good compensatory mechanisms against hypovolemia from those that might be poor compensators, leading to improved triage of trauma patients and ultimately enhancing military and emergency medicine.

60 APPLIED LIFE SCIENCES↗

Experimental and modeling study of the autoignition behavior of a saturated heterocycle: Pyrrolidine

Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.

Autoignition↗

Oxidation and pyrolysis of methyl propyl ether

The ignition, oxidation, and pyrolysis chemistry of methyl propyl ether (MPE) was probed experimentally at several different conditions, and a comprehensive chemical kinetic model was constructed to help understand the observations, with many of the key parameters computed using quantum chemistry and transition state theory. Experiments were carried out in a shock tube measuring time variation of CO concentrations, in a flow tube measuring product concentrations, and in a rapid compression machine (RCM) measuring ignition delay times. The detailed reaction mechanism was constructed using the Reaction Mechanism Generator software. Sensitivity and flux analyses were used to identify key rate and thermochemical parameters, which were then computed using quantum chemistry to improve the mechanism. Validation of the final model against the 1-20 bar 600-1500 K experimental data is presented with a discussion of the kinetics. The model is in excellent agreement with most of the shock tube and RCM data. Strong non-monotonic variation in conversion and product distribution is observed in the flow-tube experiments as the temperature is increased, and unusually strong pressure dependence and significant heat release during the compression stroke is observed in the RCM experiments. These observations are largely explained by a close competition between radical decomposition and addition to O2 at different sites in MPE; this causes small shifts in conditions to lead to big shifts in the dominant reaction pathways. The validated mechanism was used to study the chemistry occurring during ignition in a diesel engine, simulated using Ignition Quality Test (IQT) conditions. At the IQT conditions, where the MPE concentration is higher, bimolecular reactions of peroxy radicals are much more important than in the RCM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Data-scarce surrogate modeling of shock-induced pore collapse process

Understanding the mechanisms of shock-induced pore collapse is of great interest in various disciplines in sciences and engineering, including materials science, biological sciences, and geophysics. However, numerical modeling of the complex pore collapse processes can be costly. To this end, a strong need exists to develop surrogate models for generating economic predictions of pore collapse processes. Here, in this work, we study the use of a data-driven reduced-order model, namely dynamic mode decomposition, and a deep generative model, namely conditional generative adversarial networks, to resemble the numerical simulations of the pore collapse process at representative training shock pressures. Since the simulations are expensive, the training data are scarce, which makes training an accurate surrogate model challenging. To overcome the difficulties posed by the complex physics phenomena, we make several crucial treatments to the plain original form of the methods to increase the capability of approximating and predicting the dynamics. In particular, physics information is used as indicators or conditional inputs to guide the prediction. In realizing these methods, the training of each dynamic mode composition model takes only around 30 s on CPU. In contrast, training a generative adversarial network model takes 8 h on GPU. Moreover, using dynamic mode decomposition, the final-time relative error is around 0.3% in the reproductive cases. We also demonstrate the predictive power of the methods at unseen testing shock pressures, where the error ranges from 1.3 to 5% in the interpolatory cases and 8 to 9% in extrapolatory cases.

97 MATHEMATICS AND COMPUTING↗

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dynamic loads on twin jet exhaust nozzles due to shock noise

Acoustic near field data were collected with model single and twin jet nozzles to determine if closely spaced nozzles produce higher acoustic loading than do single nozzles. The tests were spurred by structural failure of the B-1 exhaust nozzle external flaps and similar damage on the F-15. The test was performed using two 5/8 in. ID pipes machined and placed side-by-side to mimic B-1 nozzles. A microphone mounted on the internozzle fairing measured acoustic levels near the nozzle exit plane. The nozzles oscillated significantly more than did a single nozzle over a wide range of nozzle pressure ratios. Acoustic levels in the dual jets exceeded single jet noise by as much as 20 dB, making acoustic resonance a definite candidate for structural damage in the twin jet configuration.

Norum, T. D.↗

Multi-dimensional high order essentially non-oscillatory finite difference methods in generalized coordinates

The nonlinear stability of compact schemes for shock calculations is investigated. In recent years compact schemes were used in various numerical simulations including direct numerical simulation of turbulence. However to apply them to problems containing shocks, one has to resolve the problem of spurious numerical oscillation and nonlinear instability. A framework to apply nonlinear limiting to a local mean is introduced. The resulting scheme can be proven total variation (1D) or maximum norm (multi D) stable and produces nice numerical results in the test cases. The result is summarized in the preprint entitled 'Nonlinearly Stable Compact Schemes for Shock Calculations', which was submitted to SIAM Journal on Numerical Analysis. Research was continued on issues related to two and three dimensional essentially non-oscillatory (ENO) schemes. The main research topics include: parallel implementation of ENO schemes on Connection Machines; boundary conditions; shock interaction with hydrogen bubbles, a preparation for the full combustion simulation; and direct numerical simulation of compressible sheared turbulence.

Shu, Chi-Wang↗

Measurement of material recession and shock standoff in plasma windtunnel using neural nets

Arcjets are plasma wind tunnels used to test the performance of heatshield materials for spacecraft atmospheric entry. These facilities present an extremely harsh flow environment with heat fluxes up to 10^9 W/m^2 for up to 30 minutes. The plasma is low-temperature (~1 eV) but high pressure (> 10 kPa) creating high-enthalpy supersonic flows similar to atmospheric entry conditions. Typically, material samples are measured before and after a test to characterize the total recession. However, this does not capture time-dependent effects such as material expansion and non-linear recession. This work will present new analysis of arcjet test videos which measure both the time-dependent 2D recession of the material samples and the shock standoff distance. New results showing non-linear material erosion rates will be highlighted. The material and shock edges are extracted from the videos by training and applying a convolutional neural network. Due to the consistent camera settings, the machine learning model achieves high accuracy (~99%) on new data with only a small number of training frames (~80). The new results will be discussed in the context of temperature dependent plasma-surface interaction.

Magnus A Haw↗

Measurement of Material Recession and Shock Standoff in Plasma Windtunnel using Neural Nets

Arcjets are plasma wind tunnels used to test the performance of heatshield materials for spacecraft atmospheric entry. These facilities present an extremely harsh flow environment with heat fluxes up to 109 W/m2 for up to 30 minutes. The plasma is low-temperature (∼1 eV) but high pressure (> 10 kPa) creating high-enthalpy supersonic flows similar to atmospheric entry conditions. Typically, material samples are measured before and after a test to characterize the total recession. However, this does not capture time-dependent effects such as material expansion and non-linear recession. This work will present new analysis of arcjet test videos which measure both the time-dependent 2D recession of the material samples and the shock standoff distance. The results show non-linear time-dependent effects are present for some conditions. The material and shock edges are extracted from the videos by training and applying a convolutional neural network. Due to the consistent camera settings, the machine learning model achieves high accuracy (± 2 px) relative to manually segmented images with only a small number of training frames (80).

Neural network↗

Toward a machine-guided approach to energetic material discovery

In this article, we trained a machine learning (ML) model to connect microstructural details of an energetic material formulation to its performance for the purpose of guiding the discovery of new explosive formulations. Our hypothesis was that the algorithm would robustly learn the training data and produce an accurate surrogate model. Specifically, the algorithm learned the relationship between details of the void size distribution (VSD), initiating shock pressure, and the energetic material performance. We used realistic constraints on the VSD and a range of cases were ingested by a physically informed reactive flow model working within a hydrodynamic solver running on high-performance computing resources. The ML algorithm produced a surrogate model that accurately predicted known test points around the parameter space. In addition to the utility of the model and the process used for its development, we noted interesting comparisons between what we, the authors—subject matter experts, would heuristically conclude from the training data and the surrogate model predictions. We detected nuanced details that were missed by the surrogate model; however, these details are not important to an energetic material formulator. We concluded that the algorithm did indeed robustly learn the training data and produce an accurate surrogate model. We further concluded that the surrogate model is a powerful tool to guide the formulator in the absence of subject matter experts and limited-access computing resources.

42 ENGINEERING↗

Machine learning changes the rules for flux limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations has seen recent interest for reducing the computational load for fluid simulations. Here, we focus on determining a flux-limiter for shock capturing methods. Focusing on flux limiters provides a specific plug-and-play component for existing numerical methods. Since their introduction, an array of flux limiters has been designed. Using the coarse-grained Burgers' equation, we show that flux-limiters may be rank-ordered in terms of their log-error relative to high-resolution data. We then develop a theory to find an optimal flux-limiter and present flux-limiters that outperform others tested for integrating Burgers' equation on lattices with [Formula: see text], and 2x, 3x, 4x, and 8x coarse-grainings. We train a continuous piecewise linear limiter by minimizing the mean-squared misfit to six-grid point segments of high-resolution data, averaged over all segments. While flux limiters are generally designed to have an output of φ(r) = 1 at a flux ratio of r = 1, our limiters are not bound by this rule and yet produce a smaller error than standard limiters. Here we find that our machine learned limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we use our theory to learn flux-limiters that outperform standard limiters across a range of values (as opposed to at a specific fixed value) of coarse-graining, number of discretized bins, and diffusion parameter. This demonstrates the ability to produce flux limiters that should be more broadly useful than standard limiters for general applications.

97 MATHEMATICS AND COMPUTING↗