Search NASASearch

SEARCH · Search NASA

Results for “4 point bend”

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.

612 records · Page 19

Assembly bias and local Primordial non-Gaussianity from DESI DR1 quasars

The analysis of the large-scale clustering of quasars (QSO) observed by the Dark Energy Spectroscopic Instrument (DESI) represents a promising avenue for constraining local Primordial non-Gaussianity (PNG), parameterized by f NL . The signal to be constrained is the scale-dependent bias induced in the 2-point clustering of the considered tracer sample. The resulting constraints on f NL , however, are fully degenerate with the local PNG bias parameter b ϕ , dependent on the assembly bias parameter p. Using IllustrisTNG hydrodynamical simulations, we select a QSO sample reflecting the selection criteria and properties of DESI QSOs, and provide a robust prior for p, and thus for b ϕ , building on the findings of Fondi et al. 2024. We find a distribution with mean p̅ ≃ 1.4 with weak redshift dependence, stable to selection noise and consistent with the expected recent merger history typical of quasar-hosting halos. By comparing with the CAMELS simulations we demonstrate that this prior is robust to astrophysical assumptions and cosmic variance. Finally, applying this prior to the DESI DR1 dataset, we derive updated constraints on local PNG, obtaining f NL = -3.3±9.2.

cosmological parameters from LSS

Understanding and Controlling Hybrid Electric Gas Turbine Engine Transient Dynamics

The electrification of the gas turbine engine is known to increase the flexibility of aircraft architectures by enabling the generation of electrical power to distribute to other electrically based, thrust producing subsystems. It also has potential for direct performance benefits in the gas turbine engine itself, both at steady state and dynamically. Although the design focus of the gas turbine engine performance is primarily at steady state, it is often the instabilities occurring during transients that cause disequilibrium and constrain performance improvements. Instabilities arise due to the disequilibrium of the energy storage mechanisms within the traditional engine system and likewise for the electrified engine system. The primary energy storage mechanisms in the traditional system are the rotational inertia, gas path volumes, the thermal masses that make up the mechanical structure, and now, the electrical power system will provide additional contributions. Understanding the effect each of these energy storage mechanisms has on the others and controlling them appropriately allows for the suppression of state changes within the turbomachinery components to the degree that the components remain near steady state, thus reducing the disequilibrium within the system. The ability to tightly regulate the state changes of turbomachinery components, such as the compressor, minimizes the excursion of the compressor operating point from the operating line (operability), allowing for higher performing, more efficient compressor designs by decreasing the amount of stall margin needed for safe engine power level changes. Preliminary studies with the electrification of the turbine engine have shown that this is possible, and this concept can lead to design trades benefitting engine performance, weight, and volume. This paper explores in detail the energy storage mechanisms and control approaches for coordinating their state changes, and ultimately proposes that a higher performing, more efficient compressor design can result.

transient dynamics

Low thermal budget dopant activation in shallow junctions of implanted Si via Tunable plasma enhanced annealing

Low thermal budget (LTB) annealing has become increasingly critical for shallow junctions as semiconductor devices are scaled down. Plasma Enhanced Annealing (PEA) has emerged as a promising LTB annealing process, however its mechanisms remain unclear. In this study, arsenic and boron implanted silicon wafers were subjected to helium or argon ion bombardment generated by annealing plasmas under controlled ion energies and doses in a home-built annealing reactor. The structural and electrical changes were characterized by four-point probe measurements, secondary ion mass spectrometry, Raman spectroscopy and spectroscopic ellipsometry. A pronounced reduction in sheet resistance, accompanied by a decrease in the damage layer thickness demonstrates the surface selective annealing capability of PEA process. Comparative studies using He and Ar plasmas and related ion bombardment reveal a dependence of activation efficacy on both ion species and the dopant implantation profile. Furthermore, a multi-step annealing mechanism is proposed, consisting of an initial (low dose) structural recovery stage followed by the generation of “extended” defects and concluded at large doses by an enhanced dopant activation. The proposed annealing kinetics are thought to be controlled by ion flux, dose and energy. Finally, these results highlight PEA as an effective, surface-selective, LTB approach for shallow-dopant activation and near surface annealing, and provide mechanistic insights into PEA processes.

dopant activation

Automation of Solid Oxyde Electrolyzer Cell (SOEC) & Stack Assembly

The work performed under this agreement before the No/Go decision amounted to the following: 1. Preliminary manufacturing requirements defined by business sensitivity, baseline processes & risk assessment; 2. Process & materials development, in-line gauge exploration, and preliminary automation work used to reduce risk & further refine equipment specifications; 3. “Request for Quote” issued to multiple suppliers for major equipment; 4. Application specific equipment, hardware, and fixturing will be more beneficial to fabricate inhouse; 5.Current proposals & estimates meet cycle time, capital spend, & direct labor; space is on target but requires awareness; 6. Team will continue to explore opportunities to reduce risk (dry time, traceability); 7. Go / No-Go, Purchase Orders, Equipment build & commissioning next All other future task were not completed because during the Go /No-Go decision it was confirmed this project was a "No-Go" and did not proceed past BP1; these included the following: 1. Preliminary manufacturing requirements defined by business sensitivity, baseline processes & risk assessment; 2. Process & materials development, in-line gauge exploration, and preliminary automation work used to reduce risk & further refine equipment specifications; 3. “Request for Quote” issued to multiple suppliers for major equipment; 4. Application specific equipment, hardware, and fixturing will be more beneficial to fabricate inhouse; 5. Current proposals & estimates meet cycle time, capital spend, & direct labor; space is on target but requires awareness; 6. Team will continue to explore opportunities to reduce risk (dry time, traceability) In final preparation of the termination of the Automation of Solid Oxyde Electrolyzer Cell (SOEC) & Stack Assembly, Cummins has purchased no equipment with government funds, and there is no government-owned property in Cummins possession related to this project. Minimal labor and travel were completed and paid. This acts as the Final technical report and concludes our business with DOE on this grant.

36 MATERIALS SCIENCE

Resistance of Silicon Nitride Turbine Components to Erosion and Hot Corrosion/Oxidation Attack

Silicon nitride turbine components are under intensive development by AlliedSignal to enable a new generation of higher power density auxiliary power systems. In order to be viable in the intended applications, silicon nitride turbine airfoils must be designed for survival in aggressive oxidizing combustion gas environments. Erosive and corrosive damage to ceramic airfoils from ingested sand and sea salt must be avoided. Recent engine test experience demonstrated that NT154 silicon nitride turbine vanes have exceptional resistance to sand erosion, relative to superalloys used in production engines. Similarly, NT154 silicon nitride has excellent resistance to oxidation in the temperature range of interest - up to 1400 C. Hot corrosion attack of superalloy gas turbine components is well documented. While hot corrosion from ingested sea salt will attack silicon nitride substantially less than the superalloys being replaced in initial engine applications, this degradation has the potential to limit component lives in advanced engine applications. Hot corrosion adversely affects the strength of silicon nitride in the 850 to 1300 C range. Since unacceptable reductions in strength must be rapidly identified and avoided, AlliedSignal and the NASA Lewis Research Center have pioneered the development of an environmental life prediction model for silicon nitride turbine components. Strength retention in flexure specimens following 1 to 3300 hour exposures to high temperature oxidation and hot corrosion has been measured and used to calibrate the life prediction model. Predicted component life is dependent upon engine design (stress, temperature, pressure, fuel/air ratio, gas velocity, and inlet air filtration), mission usage (fuel sulfur content, location (salt in air), and times at duty cycle power points), and material parameters. Preliminary analyses indicate that the hot corrosion resistance of NT154 silicon nitride is adequate for AlliedSignal's initial engine applications. Protective coatings and/or inlet air filtration may be required to achieve required ceramic component lives in more aggressive environments.

Nitride

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]

Development, Characterization, and Validation of Elevated-Temperature Constitutive Models: Deformation and Damage

This report provides a brief review of experimentally observed hereditary and nonhereditary material behavior along with background information on standard as well as advanced internal state variable constitutive modeling at elevated temperature. A description of exploratory, characterization, and validation testing is presented along with a detailed outline of what constitutes “sufficient” data content (i.e., quality and quantity) for developing or enhancing, characterizing, and validating a particular sophisticated nonlinear time- and history-dependent (hereditary) class of constitutive models known as GVIPS (generalized viscoplasticity with potential structure). The tests described are necessary to reveal a material’s behavior in the reversible (or viscoelastic) and irreversible (or viscoplastic) regimes, for the identification of both deformation and damage model parameters. Results presented are primarily for metallic materials. In all cases, both uniaxial and multiaxial tests are described, and the linkage between the specific tests and the parameters within the model that can be characterized from the results of these tests are also discussed. Discussion is also provided relative to the role information management must play relative to material data collection, analysis, maintenance, and dissemination. The need for such an information system is particularly important as both analyst and designer move toward utilizations of sophisticated, nonlinear time- and history dependent (hereditary) constitutive models. Lastly, the concept of state space and its utility in understanding and establishing constitutive models is addressed throughout. The intent behind this document is to help both the modeler and experimentalist understand each other’s specific points of view and provide guidance for both model development and characterization. Emphasis has been placed on providing the mechanician with information regarding how tests are performed and what issues to be aware of when interpreting results. It is hoped that experimentalists will take away a new perspective on the types of information that modelers and analysts are looking for from them.

Constitutive Modeling

Enhanced Decoding for the Galileo S-Band Mission

A coding system under consideration for the Galileo S-band low-gain antenna mission is a concatenated system using a variable redundancy Reed-Solomon outer code and a (14,1/4) convolutional inner code. The 8-bit Reed-Solomon symbols are interleaved to depth 8, and the eight 255-symbol codewords in each interleaved block have redundancies 64, 20, 20, 20, 64, 20, 20, and 20, respectively (or equivalently, the codewords have 191, 235, 235, 235, 191, 235, 235, and 235 8-bit information symbols, respectively). This concatenated code is to be decoded by an enhanced decoder that utilizes a maximum likelihood (Viterbi) convolutional decoder; a Reed Solomon decoder capable of processing erasures; an algorithm for declaring erasures in undecoded codewords based on known erroneous symbols in neighboring decodable words; a second Viterbi decoding operation (redecoding) constrained to follow only paths consistent with the known symbols from previously decodable Reed-Solomon codewords; and a second Reed-Solomon decoding operation using the output from the Viterbi redecoder and additional erasure declarations to the extent possible. It is estimated that this code and decoder can achieve a decoded bit error rate of 1 x 10 -7 at a concatenated code signal-to-noise ratio of 0.76 dB. By comparison, a threshold of 1.17 dB is required for a baseline coding system consisting of the same (14,1/4) convolutional code, a (255,223) Reed-Solomon code with constant redundancy 32 also interleaved to depth 8, a one-pass Viterbi decoder, and a Reed Solomon decoder incapable of declaring or utilizing erasures. The relative gain of the enhanced system is thus 0.41 dB. It is predicted from analysis based on an assumption of infinite interleaving that the coding gain could be further improved by approximately 0.2 dB if four stages of Viterbi decoding and four levels of Reed-Solomon redundancy are permitted. Confirmation of this effect and specification of the optimum four-level redundancy profile for depth-8 interleaving is currently being done.

S Dolinar

Speciation of Transition Metal Dissolution in Electrolyte from Common Cathode Materials

Significant capacity loss has been observed across extended cycling of lithium-ion batteries cycled to high potential. One of the sources of capacity fade is transition metal dissolution from the cathode active material, ion migration through the electrolyte, and deposition on the solid-electrolyte interphase on the anode. While much research has been conducted on the oxidation state of the transition metal in the cathode active material or deposited on the anode, there have been limited investigations of the oxidation state of the transition metal ions dissolved in the electrolyte. Here, in this work, X-ray absorption spectroscopy has been performed on electrolytes extracted from cells built with four different cathode active materials (LiMn 2 O 4 (LMO), LiNi 0.5 Mn 1.5 O 4 (LNMO), LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), and (x Li 2 MnO 3 *(1-x) LiNi a Mn b Co c O 2 , with a+b+c=1) (LMRNMC)) that were cycled at either high or standard potentials to determine the oxidation state of Mn and Ni in solution. Inductively coupled plasma-mass spectrometry has been performed on the anodes from these cells to determine the concentration of deposited transition metal ions. While transition metal ions were found dissolved in all electrolytes, the oxidation state(s) of Mn and Ni were determined to be dependent on the cathode material and independent of cycling potential.

25 ENERGY STORAGE

Localizing the non-thermal X-ray emission of PSR J2229+6114 from its multi-wavelength pulse profiles

Context.Pulsars are detected over the whole electromagnetic spectrum, from radio wavelengths up to very high energies, in the GeV-TeV range. While the radio emission site for young pulsars is well constrained to occur at altitudes about several percent of the light-cylinder radius andγ-ray emission is believed to be produced in the striped wind, outside the light cylinder, the non-thermal X-ray production site remains unknown. Aims.The aim of this Letter is to localize the non-thermal X-ray emission region based on a multi-wavelength pulse profile fitting for PSR J2229+6114, which stands as a particularly good candidate given its high X-ray brightness. Methods.Based on the geometry deduced from the joint radio andγ-ray pulse profiles, we fixed the magnetic axis inclination angle and the line-of-sight (LoS) inclination angle. However, we left the region of X-ray emission unlocalized, setting it somewhere between the surface and the light cylinder. We localized this region and its extension by fitting the X-ray pulse profile as observed by the NICER,NuSTAR, and RXTE telescopes in the ranges of 2–7 keV, 3–10 keV, and 9.4–22.4 keV, respectively. Results.We constrained the non-thermal X-ray emission to arise from altitudes between 0.2 r L and 0.55 r L wherer L is the light-cylinder radius. The magnetic obliquity is approximatelyα ≈ 45° −50° and the LoS inclination angle isζ ≈ 32° −48°. Conclusions.This Letter is among the first works to offer a tight constraint on the location of non-thermal X-ray emission from pulsars. We plan to apply this procedure to several other promising candidates to confirm this new result.

Astronomy & Astrophysics

In-suit CO2 Washout Test System (ICWTS) for CO2 Washout Verification in Spacesuits

Since the time of Mercury, Gemini, and Apollo, NASA has been performing CO 2 washout verification of spacesuits with Human in the Loop (HITL) test methods. Classically, the only instrumentation both with the accuracy and response time necessary was that of external mass spectrometers or medical gas analyzers fed by a long capillary line from the oro-nasal region within the suit. The helmet CO 2 washout configuration including the challenges posed by variation of a suited human subject frequently created difficulty in capturing adequate breath by breath data. For example, a nasal cannula could be worn with reasonable comfort however it posed the following challenges: (1) impeded the test subject’s ability to do a proper Valsalva maneuver; (2) frequently provided a corrupted waveform with relative movement; (3) precluded the ability to evaluate the Flow Weighted Average (FWA) at subject interface. A mouthpiece could be used in lieu of a cannula but posed challenges as well: (1) created potential disturbances of the flow field; (2) precluded the test subject from speaking during usage; and (3) presented comfort issues limiting the duration of test points. These challenges made it impractical to use a mouthpiece for anything other than short lab ambient environment tests. Technology development and NASA investment have afforded new capabilities to address some of these measurements providing improved data collection and are implemented in the In-suit CO 2 Washout Test System (ICWTS) which provides suit inlet/outlet pressure, temperature, relative humidity, partial pressure of oxygen, partial pressure of CO 2 , and flow measurement coupled with the acquisition of the inspired PPCO 2 including the ability for FWA measurement. This test system, coupled with use of the In-suit Respiration Mannequin Assembly (IRMA) and Computational Fluid Dynamics (CFD) combine to provide a full picture of the in-suit oro-nasal flow field and effectiveness of CO 2 washout in a spacesuit design.

Colin Campbell

In-suit CO 2 Washout Test System (ICWTS) for CO 2 Washout Verification in Spacesuits

Since the time of Mercury, Gemini, and Apollo, NASA has been performing CO 2 washout verification of spacesuits with Human in the Loop (HITL) test methods. Classically, the only instrumentation both with the accuracy and response time necessary was that of external mass spectrometers or medical gas analyzers fed by a long capillary line from the oro-nasal region within the suit. The helmet CO 2 washout configuration including the challenges posed by variation of a suited human subject frequently created difficulty in capturing adequate breath by breath data. For example, a nasal cannula could be worn with reasonable comfort however it posed the following challenges: (1) impeded the test subject’s ability to do a proper Valsalva maneuver; (2) frequently provided a corrupted waveform with relative movement; (3) precluded the ability to evaluate the Flow Weighted Average (FWA) at subject interface. A mouthpiece could be used in lieu of a cannula but posed challenges as well: (1) created potential disturbances of the flow field; (2) precluded the test subject from speaking during usage; and (3) presented comfort issues limiting the duration of test points. These challenges made it impractical to use a mouthpiece for anything other than short lab ambient environment tests. Technology development and NASA investment have afforded new capabilities to address some of these measurements providing improved data collection and are implemented in the In-suit CO 2 Washout Test System (ICWTS) which provides suit inlet/outlet pressure, temperature, relative humidity, partial pressure of oxygen, partial pressure of CO 2 , and flow measurement coupled with the acquisition of the inspired PPCO 2 including the ability for FWA measurement. This test system, coupled with use of the In-suit Respiration Mannequin Assembly (IRMA) and Computational Fluid Dynamics (CFD) combine to provide a full picture of the in-suit oro-nasal flow field and effectiveness of CO 2 washout in a spacesuit design.

Colin Campbell

Mars Sample Retrieval Lander Thermal Protection System Design

The Mars Sample Return (MSR) Sample Retrieval Lander (SRL) was slated to launch in 2031 and enter Mars’ atmosphere after a 15-month cruise. After several iterations of vehicle design and trajectories, the latest mission architecture included a 4.72-meter diameter entry capsule entering the Martian atmosphere at 8 km/s, marking a 3.5x increase in kinetic energy when entering the atmosphere in comparison to Mars 2020. Designed to be the heaviest and fastest vehicle to enter Mars’ atmosphere, SRL required slim mass margins, and as such, the thermal protection system (TPS) design methodology needed to be tailored appropriately. The entry capsule aeroshell is composed of a heatshield, backshell, parachute support cone, and parachute lid, as shown in Fig. 1.

Hannah Alpert

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

stagnation heat flux

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory

CatTestHub: A benchmarking database of experimental heterogeneous catalysis for evaluating advanced materials

The ability to quantitatively compare newly evolving catalytic materials and technologies is hindered by the widespread availability of catalytic data collected in a consistent manner. While certain catalytic chemistries have been widely studied across decades of scientific research, quantitative comparisons based on literature information is hindered by variability in reaction conditions, types of reported data, and reporting procedures. Here, we present CatTestHub, an open-access database dedicated to benchmarking experimental heterogeneous catalysis data. Combining systematically reported catalytic activity data for selected probe chemistries, with relevant material characterization and reactor configuration information, the database provides a collection of catalytic benchmarks for distinct classes of active site functionality. Through key choices in data access, availability, and traceability, CatTestHub seeks to balance the fundamental information needs of chemical catalysis and the FAIR data design principles. Details of the database architecture and the means through which to navigate it are presented, highlighting examples of catalytic insights readily drawn from the available benchmarking data. In its current iteration, CatTestHub spans over 250 unique experimental data points, collected over 24 solid catalysts, that facilitated the turnover of 3 distinct catalytic chemistries. Here, a roadmap is presented through which to expand the open-access platform that serves as a community wide benchmark, primarily through continuous addition of kinetic information on select catalytic systems by members of the heterogeneous catalysis community at large.

Benchmark

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC