Search NASASearch

SEARCH · Search NASA

Results for “Discrete Element Method”

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.

417 records · Page 20

ClassNMSW- a real-time classification approach for non-recycled municipal solid waste using hyperspectral imaging

Real-time classification of non-recycled municipal solid waste (NMSW) is essential for efficient valorization. This study introduces ClassNMSW, a comprehensive framework for classifying 22 NMSW subclasses under industrial constraints by using hyperspectral imaging (HSI). A primary innovation of this work is the development of a variance-controlled spectral extraction algorithm. Unlike traditional methods that rely on simple averaging, this approach systematically investigates the extent of pixel extraction to minimize the loss of critical chemical information while maximizing data reduction thus ensuring high spectral fidelity with low computational cost. The approach developed in this work integrates automated, computer-vision-based background removal, eliminating the need for the manual thresholding common in current literature. To resolve ambiguities among chemically similar subclasses, a tiered classification and multi-camera fusion strategy (NIR17 and NIR22) is implemented. Results demonstrate that ClassNMSW achieves an object-wise weighted accuracy of 98.70% for single-sensor configurations and 100% under sensor fusion. A novel rolling-window strategy satisfies desired end-to-end latency of <2 s, satisfying the strict deterministic requirements of high-speed industrial sorting environments. The ClassNMSW framework provides a scalable foundation for advancing circularity and resource recovery in large-scale waste valorization operations.

99 - GENERAL AND MISCELLANEOUS

High-Intensity UV Exposure for the Rapid Screening of Silicon Photovoltaic Architectures

Advanced Si photovoltaic architectures incorporate different materials and processing pathways that influence degradation modes. Ultraviolet-induced degradation (UVID) is an understudied degradation mode for advanced cell architectures and is of increasing concern to industry due to growing adoption of UV-transparent encapsulation and bifacial technologies. In order to adopt new and evolving technologies confidently, novel component materials and processing techniques must be evaluated and designed for long-term stability, in addition to the conventional design focus on efficiency. In this work, a study protocol framework is presented for the rapid screening of unencapsulated devices against UVID. Unencapsulated passivated emitter rear contact (PERC) and tunnel oxide passivated contact (TOPCon) devices were aged under different UV irradiance intensities and measured via conventional nondestructive electrical characterization methods to assess performance degradation. Based on the results, protocol efficacy and recommendations for further study are discussed. As a result, this work is part of a broader effort to develop rapid screening processes that cut across architectures and exposure conditions to aid module manufacturers in vetting new materials choices for long-term stability.

Accelerated exposure

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

corrosion

Coercivity enhancement in hematite/permalloy heterostructures across the Morin transition

Interfacial effects between antiferromagnetic (AFM) and ferromagnetic (FM) materials have long been a center of magnetism studies. Aside from the exchange bias occurring at the AFM/FM interface, controlling the coercivity is another significant topic in magnetic recordings. The coercivity of FM materials is often determined through varying grain size, alloy composition, density of defects, etc., which is set during material growth and offers limited room for modification after growth. Hematite (α - Fe 2 O 3 ) is an AFM material that undergoes a temperature-controlled spin-flip transition, the so-called Morin transition. This transition gives an extra degree of freedom making hematite an intriguing component to study the exchange coupling when interfaced with an FM material. Here, in this work, changes in the magnetic properties of soft magnetic permalloy (Ni 81 Fe 19 , or Py) thin films grown on hematite were studied across the Morin transition. Surprisingly, these samples showed a remarkable change in coercivity during the Morin transition. We attribute this effect to the magnetic domain mixture of hematite during the Morin transition. Our findings present a novel method of controlling the coercivity of plain ferromagnetic thin films.

Antiferromagnet

Rinse-Free, Sodium-Efficient Synthesis of O3-Type Layered Oxide Materials Enabled by Acetate Precursors

Sodium-ion batteries (SIBs) are a sustainable alternative to lithium-ion systems for global electrification, with O3-type layered oxide cathodes offering high specific capacity and feasibility of scalable synthesis. Industrial co-precipitation synthesis of these cathodes typically uses transition metal sulfates, requiring extensive water rinsing to remove Na 2 SO 4 impurities, a process that consumes significant water and risks residual inactive phases if incomplete. Here, this work introduces a rinse-free, resource-efficient approach using metal acetate precursors. Residual sodium acetate in non-rinsed precursors decomposes during sintering to generate Na 2 CO 3 in situ, partially substituting an external sodium resource (e.g., NaOH and Na 2 CO 3 ) and reducing its consumption by ∼18–20%. Phase-pure O3-Na 1.0 Ni 1/3 Fe 1/3 Mn 1/3 O 2 (NFM111) cathodes synthesized via this method exhibit microstructure and electrochemical performance comparable to rinsed sulfate-derived counterparts, with initial capacities of 141 mAh g −1 at C/20 (7.5 mA g −1 ). By eliminating rinsing and minimizing sodium reagent use, this acetate-based route enhances sustainability and scalability of layered oxide production for SIBs.

acetate vs. sulfate

Vanishing dynamic strength measured in a transient high-pressure phase of tin

This work reports two independent experimental estimations of the dynamic strength of the body-centered tetragonal (bct) 𝛾 phase of tin, which occurs at pressures above about 9 GPa and transforms promptly back to the ambient 𝛽 phase upon pressure release. Measuring strength in such a transient high-pressure phase is challenging. One measurement used free-surface Richtmyer-Meshkov instabilities generated with gas gun impact. Another set of measurements used ramp-release loading in a pulsed-power facility. Strength estimations came from comparison to simulations using a comprehensive multiphase modeling framework that includes equation of state, shear moduli, and strength separately for each phase, and treats mixed-phase regions. Both experimental methods found the 𝛾-phase strength to be very low, within experimental uncertainty of zero. The existing literature on other materials, by contrast, almost universally reports higher strength in transient high-pressure phases compared with ambient phases. Recent Molecular-Dynamics simulations on tin in the literature showed almost zero deviatoric strength during a deformation-induced bct → bct transformation in which one of the 𝛾-phase 𝑎 axes “flips” to the 𝑐 axis. This reorientation currently provides the most plausible explanation for the low observed strength in 𝛾-phase tin.

36 MATERIALS SCIENCE

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen

Why Perovskite Thermal Stress is Unaffected by Thin Contact Layers

Metal halide perovskite photovoltaics have emerged as a high efficiency, low-cost alternative that can potentially rival or enhance conventional silicon technology. Despite exceptional initial power conversion efficiencies, achieving compliance with international standards and widespread adoption requires further enhancements to their operational stability. Notably, addressing mechanical strain and stress in brittle perovskites has emerged as a pivotal approach to mitigate chemical degradation and improve reliability during thermal cycling. Here, in this study, a popularized strain engineering strategy is investigated in which a high coefficient of thermal expansion (CTE) hole transport layer (i.e., PDCBT) is cast onto inorganic perovskite (CsPbI 2 Br) at 100 °C. Contrary to previously published results, the X-ray diffraction (XRD):Sin 2 ψ and substrate curvature measurement techniques show that the hole transport layer has no discernible impact on perovskite strain. The accuracy of the XRD:Sin 2 ψ method for measuring strain is highlighted in contrast to an analysis based on shifts of single XRD peaks which can be influenced by multiple artifacts. The findings in this study are in accordance with mechanics theory: thin layers are unable to induce significant strain changes in perovskite thin films as the force they apply is negligible compared to that applied by a thick and stiff substrate.

14 SOLAR ENERGY

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew

Patchy nanoparticles by atomic stencilling

Stencilling, in which patterns are created by painting over masks, has ubiquitous applications in art, architecture and manufacturing. Modern, top-down microfabrication methods have succeeded in reducing mask sizes to under 10 nm, enabling ever smaller microdevices as today’s fastest computer chips. Meanwhile, bottom-up masking using chemical bonds or physical interactions has remained largely unexplored, despite its advantages of low cost, solution-processability, scalability and high compatibility with complex, curved and three-dimensional (3D) surfaces. Here we report atomic stencilling to make patchy nanoparticles (NPs), using surface-adsorbed iodide submonolayers to create the mask and ligand-mediated grafted polymers onto unmasked regions as ‘paint’. We use this approach to synthesize more than 20 different types of NP coated with polymer patches in high yield. Polymer scaling theory and molecular dynamics (MD) simulation show that stencilling, along with the interplay of enthalpic and entropic effects of polymers, generates patchy particle morphologies not reported previously. These polymer-patched NPs self-assemble into extended crystals owing to highly uniform patches, including different non-closely packed superlattices. We propose that atomic stencilling opens new avenues in patterning NPs and other substrates at the nanometre length scale, leading to precise control of their chemistry, reactivity and interactions for a wide range of applications, such as targeted delivery, catalysis, microelectronics, integrated metamaterials and tissue engineering.

36 MATERIALS SCIENCE

The POINTER Imaging baseline cohort: Associations between multimodal neuroimaging biomarkers, cardiovascular health, and cognition

Abstract INTRODUCTION The U.S. Study to Protect Brain Health Through Lifestyle Intervention to Reduce Risk (U.S. POINTER) is evaluating lifestyle interventions in older adults at risk for cognitive decline and dementia. Here we characterize the baseline data set of the POINTER Imaging ancillary study. METHODS Participants underwent health and cognitive assessments and neuroimaging with multimodal positron emission tomography (PET) (beta‐amyloid [Aβ] and tau) and magnetic resonance imaging (MRI). Framingham risk score (FRS) was used to quantify cardiovascular disease (CVD) risk. RESULTS A total of 1052 participants (31% from underrepresented ethnoracial groups) were enrolled. Compared to Aβ−, Aβ+ (29%) participants were older, had higher apolipoprotein E (APOE) ε4 carriage rate and white matter hyperintensity volume, and greater temporal tau. FRS was related to MRI measures, but not AD biomarkers. FRS and tau had independent effects on cognition. DISCUSSION In this heterogenous, at‐risk cohort, CVD risk was related to more abnormal brain structure and poorer cognition, representing a putative non‐AD (Alzheimer's disease) pathway to brain injury and cognitive decline. Highlights The U.S. Study to Protect Brain Health Through Lifestyle Intervention to Reduce Risk (U.S. POINTER) cohort is enriched for cardiovascular disease (CVD) and poor lifestyle POINTER Imaging collected multimodal neuroimaging data in this unique, at‐risk cohort Amyloid burden was related to age, apolipoprotein E (APOE) ε4 carriage, and measures of disease progression Associations between amyloid and tau, and tau and cognition, were relatively weak CVD risk and tau pathology were independently related to memory

Neurosciences & Neurology

AutoSourceID-Classifier: Star-galaxy classification using a convolutional neural network with spatial information

Aims.Traditional star-galaxy classification techniques often rely on feature estimation from catalogs, a process susceptible to introducing inaccuracies, thereby potentially jeopardizing the classification’s reliability. Certain galaxies, especially those not manifesting as extended sources, can be misclassified when their shape parameters and flux solely drive the inference. We aim to create a robust and accurate classification network for identifying stars and galaxies directly from astronomical images. Methods.The AutoSourceID-Classifier (ASID-C) algorithm developed for this work uses 32x32 pixel single filter band source cutouts generated by the previously developed AutoSourceID-Light (ASID-L) code. By leveraging convolutional neural networks (CNN) and additional information about the source position within the full-field image, ASID-C aims to accurately classify all stars and galaxies within a survey. Subsequently, we employed a modified Platt scaling calibration for the output of the CNN, ensuring that the derived probabilities were effectively calibrated, delivering precise and reliable results. Results.We show that ASID-C, trained on MeerLICHT telescope images and using the Dark Energy Camera Legacy Survey (DECaLS) morphological classification, is a robust classifier and outperforms similar codes such as SourceExtractor. To facilitate a rigorous comparison, we also trained an eXtreme Gradient Boosting (XGBoost) model on tabular features extracted by SourceExtractor. While this XGBoost model approaches ASID-C in performance metrics, it does not offer the computational efficiency and reduced error propagation inherent in ASID-C’s direct image-based classification approach. ASID-C excels in low signal-to-noise ratio and crowded scenarios, potentially aiding in transient host identification and advancing deep-sky astronomy.

Astronomy & Astrophysics

Sintering effects on NpO 2 grain size and morphology: The role of precursor Np phase

Neptunium dioxide (NpO 2 ) is a key phase in nuclear material processing as a target material for the production of plutonium-238 ( 238 Pu) and has been historically synthesized via the calcination of a Np oxalate precursor. Alternative synthesis methods for NpO 2 are now more prevalent, necessitating their study and comparison with the more common oxalate route. The purpose of this work was to investigate the microstructural properties of NpO 2 synthesized from a nitrate-based Np precursor phase via the assessment of NpO 2 particle size and morphology as a function of calcination temperature and time. Scanning electron microscopy (SEM) was used to probe the primary grain size and morphology of NpO 2 after calcination at temperatures ranging from 700 to 1100 °C and hold times ranging from 1 to 10 h. Post-image analysis using ImageJ software enabled the quantification of mean particle diameter. This analysis indicated that particle diameter increases with both increasing calcination temperature and hold time. Primary particles were shown to be clumped in irregular patterns into the overall rough, blocky aggregates, but this macroscopic morphology was not affected by calcination time or temperature. Although trends in primary grain size of NpO 2 were consistent with available literature from other Np precursor phases, the macroscopic morphology of the NpO 2 aggregates was quite different than reported for other precursors. Through comparison with historical literature on Np oxalate, this work emphasizes the importance of Np precursor on the physical properties of NpO 2 .

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA

Validating a Dynamic PWR Safety and Security Model?

Nuclear power plants (NPPs) are assessed for safety and security using separate models that cannot capture how an attacker's decisions and a plant's response unfold together in real time, leaving regulators and operators without a complete picture of true plant vulnerability. Traditional probabilistic risk assessment (PRA) methods treat adversarial events as fixed initiators with predetermined outcomes, and are structurally incapable of representing the time-dependent interplay between physical security events, safety system response, and operator mitigative actions. At Idaho National Laboratory (INL), I contributed to the development and validation of Modeling and Analysis for Safety and Security using the Dynamic EMRALD Framework (MASS-DEF). Where static PRA relies on event-tree logic that cannot evolve mid-scenario, MASS-DEF couples a time-dependent dynamic PRA tool EMRALD (Event Modeling Risk Assessment using Linked Diagrams) with attack simulation software, allowing attacker behavior, plant system states, and operator actions to interact across time. My work focused on validating a general Pressurized Water Reactor (PWR) model. I traced model logic against PWR plant to identified errors in logic and confirm accuracy. I then built and tested attack scenarios against a general PWR model to verify that the model produced expected outcomes across all logical pathways. I also contributed a section to a related technical paper applying the same EMRALD platform to radiation dose modeling. Results show that MASS-DEF can quantitatively demonstrate that many plants exceed their regulatory security thresholds. This demonstrated margin provides a technically defensible basis for reducing the number of guards without compromising regulatory compliance. Physical security costs represent roughly 10% of annual operating budgets, making such reductions directly meaningful to INL's mission of sustaining existing commercial NPPs. This internship strengthened my understanding of nuclear systems, probabilistic modeling, and technical writing, and has solidified my pursuit of a career at a national laboratory.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL

Influence of functional additives, fillers, and pigments on thermal and catalytic pyrolysis of polyethylene for waste plastic upcycling

Pyrolysis offers a relatively green and economical method to convert waste plastics into valuable chemicals and fuels without the need for harmful solvents, toxic chemicals, or costly high-pressure reactors. Despite its popularity among chemical upcycling technologies, industrial adoption suffers from feedstock heterogeneity, low-quality products, and catalyst deactivation. Most plastics in our daily lives are formulated with functional additives, fillers, and colorants. These additives remaining in end-of-life waste streams increase feedstock heterogeneity, creating a challenging issue in recycling plastics. Still, the potential impacts of additives on the chemical upcycling of plastics have been poorly understood. In this study, polyethylene compounded with a range of widely used additives (antioxidants, stabilizers, pigments, fillers, slip agents, and flame retardants) was subjected to both thermal pyrolysis and catalytic pyrolysis in different catalyst-to-feedstock contact modes. It showed that many inorganic additives, such as talc, kaolin, CaCO 3 , TiO 2 , carbon black, and zinc stearate, facilitated polymer decomposition during pyrolysis, increasing light hydrocarbons while also promoting aromatic and carbon residue formation. Conversely, antioxidants and stabilizers inhibited depolymerization, favoring heavier hydrocarbons. During catalytic pyrolysis with HZSM-5 zeolite, additives strongly enhanced aromatic and catalytic coke formation, especially when there was direct contact between plastics and catalysts. Although certain additives seem beneficial in the short term by promoting polymer cracking and improving the selectivity of aromatics, the transport of the additives and their degradation products and increased carbon coking can contaminate products, deactivate or modify catalysts, and foul reactors. These findings address a critical knowledge gap in effectively converting waste plastics via a greener route.

42 ENGINEERING

Polycarbonate‐Based Solid‐Polymer Electrolytes for Solid‐State Sodium Batteries

Solid-polymer electrolytes comprised of polypropylene carbonate (PPC) and varied sodium bis(fluorosulfonyl)imide (NaFSI) salt concentrations are investigated for implementation as a conductive solid polymer electrolyte into solid-state cathode composites utilizing a sodium-layered oxide active material. The ionic conductivity generally increases with NaFSI salt content, reaching ≈1 mS cm −1 at 80 °C at the highest salt concentration (PPC:NaFSI = 0.5:1). Through an all-in-one slurry casting method, Na 2/3 Ni 1/3 Mn 2/3 O 2 cathode composites are fabricated in which the dispersed PPC electrolyte acts as the primary binder. Enabled by a bilayer polymer electrolyte system, cycling performance with the PPC cathode electrolyte is optimized with respect to salt concentration and anode material. The best cyclability is achieved with a moderate salt concentration electrolyte (PPC:NaFSI = 5:1), showcasing an initial capacity of 83 mA h g −1 with a remarkable 80% capacity retention after 150 cycles at C/5 rate and 60 °C. The superior performance of the lower salt concentration electrolyte is attributed to better electrochemical stability, as confirmed by linear sweep voltammetry and online electrochemical mass spectrometry measurements. In conclusion, these results underscore the potential of carbonate-based polymer electrolytes and the importance of balancing electrolyte conductivity and stability in cell design.

25 ENERGY STORAGE

3D-Printed mmWave Quasi-Holographic Antenna for 2-D Beamforming

This paper presents the design and manufacturing of a novel 2D-scanning antenna that integrates a 3D-printed Rotman lens with a quasi-holographic leaky-wave antenna (HLWA). The proposed design achieved beam-scanning capabilities by leveraging the beamforming of the Rotman lens and the high-gain directive properties of the quasi-HLWA. The Rotman lens (RL) enables beam steering in the elevation plane by switching between input ports. The quasi-HLWA, designed using holographic principles, achieves frequency-controlled beam scanning in the azimuth plane. The entire antenna structure was fabricated using additive manufacturing with an Ink1092 substrate and silver ink for the conductive traces. This approach provides greater control over material placement and design freedom compared to traditional methods. A 25° transmission linear substrate taper was used to ensure good impedance matching between the Rotman lens and the quasi-HLWA, allowing greater gain while maintaining a good scanning range. The experimental results validate the 2D scanning capability of the proposed antenna. The antenna system provides coverage from −54° to 54° in the elevation θ plane and −28° to 28° in the azimuth plane ϕ , with a maximum measured gain of 21.3 dBi at 28 GHz with an average radiation efficiency η=60 %. The fabricated prototype was tested, and the performance was in good agreement with the simulated performance.

Rotman lens

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