SEARCH · Engineering Papers

Results for “experimental data”

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

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Plant-wide modeling and techno-economic analysis of a direct non-oxidative methane dehydroaromatization process via conventional and microwave-assisted catalysis

Direct non-oxidative methane dehydroaromatization (DHA) process via conventional and microwave (MW)-assisted thermo-catalytic catalysis is studied. Rate models for methane DHA reactions, including the effect of catalyst deactivation, are developed by using the in-house experimental data. Model results for gas concentration profile and catalyst deactivation are in good agreement with the experimental data. This rate model is then used for the development of dynamic multi-scale, multi-physics commercial-scale reactor models. Total number of fixed bed reactors desired for a cyclic steady state process is estimated. Plant-wide models are then developed for conventional and MW-assisted processes for producing products of desired specifications. Techno-economic analysis of the methane DHA process is undertaken. Economics of these methane DHA processes are compared with the typical multi-step natural gas to aromatics production process via methanol synthesis. Sensitivity of internal rate of return (IRR) and net present value (NPV) to various economic and process parameters such as plant scale, desired rate of return, reactor cost, feedstock and utility cost, catalyst variable cost, and MW reactor cost is studied. Here, electric equivalent efficiency of the conventional methane DHA process is found to be 69.2 % and 67.3 % at 750 °C and 800 °C, respectively, while the MW-assisted methane DHA process has the electric equivalent efficiency of 48.9 % at 800 °C. IRRs of the conventional methane DHA process at 750 °C and 800 °C, and MW-assisted process are 15.2 %, 17.5 %, and 18.8 %, respectively for a methane feed flowrate of 19,782 kg/h, while the IRR of the multi-step natural gas to aromatics production process is estimated to be 0 % for the same plant scale. Impact of change in the methane price, electricity price, and catalyst cost is found to be considerable on the process economics, while the cost of the MW reactor is found to have negligible impact.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

A RELAP5-3D Model of the Lobo Lead Loop

The goal of our research is to build upon the capability of RELAP5-3D to model molten lead systems. Molten lead has several potential uses in future advanced reactors like the lead fast reactor or fusion reactors that utilize dual-coolant lead lithium blankets. This potential for use in future generations of reactors highlights the necessity to develop molten lead models to ensure that they can accurately predict the thermohydraulic behavior. We have developed a RELAP5-3D model of the Lobo Lead Loop facility, located at the University of New Mexico, to verify the accuracy of RELAP5-3D via comparison to existing computational fluid dynamics results and analytical calculations. It was found that RELAP5-3D accurately calculated radiative heat transfer (within < 1%) when compared to theoretical calculations. In addition, pressure drop calculations done in RELAP5-3D demonstrated reasonable agreement within 20 kPa, mostly within ~7-15%, when compared to the computational fluid dynamics model of the facility developed by the University of New Mexico, and captured the dependence of pressure drop on flow velocity accurately. Finally, a hypothetical loss of flow transient was imposed on the RELAP5-3D model to determine the feasibility of performing a similar experiment with the Lobo Lead Loop. It was found that such an experiment could be possible as the RELAP5-3D model indicated that the temperatures of the fluid would not exceed limiting temperatures of the structure (1658 K) nor the maximum temperature of the electromagnetic pump inlet (823 K). Although there is not experimental data to begin validation, the model will be readily available for future validation studies when the experimental data is generated, especially as the model continues to evolve over time. Furthermore, the results so far demonstrate a promising first step in the verification/validation of the RELAP5-3D model of the Lobo Lead Loop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

An experimental and kinetic modeling study of the pyrolysis of isoprene, a significant biogenic hydrocarbon in naturally occurring vegetation fires

Isoprene dominates the carbon flux emitted by vegetation and constitutes 40% of non-methane biogenic emissions worldwide. Despite pyrolysis experiments at temperatures above 1000 K showing a link between isoprene combustion and aromatic species formation, comprehensive mechanistic research on isoprene is scarce in the literature. Here, in this work, we carry out an experimental and theoretical study to build, for the first time, a chemical kinetic model describing isoprene pyrolysis. The formation of polycyclic aromatic hydrocarbon (PAH) precursor species, often observed in vegetation fire plumes, is partially explained by isoprene pyrolysis experiments and theoretical modeling. Molecular dynamics (MD) simulations unveil reaction pathways from allylic isoprenyl radicals to allene and cyclopentadiene (CPD) intermediates, two relevant species detected in the experiments. Rate constants for these identified pathways are calculated using variational transition state theory to update the kinetic model, which is validated against single-pulse shock tube (SPST), and jet-stirred reactor (JSR) experimental data in the temperature range of 850–1690 K. The kinetic model presents satisfactory agreement with the SPST experimental data, and a reaction pathway analysis shows that association of propargyl radicals results in benzene formation. The JSR pathway analysis also identifies the prominent reactions for CPD, benzene, styrene, and toluene formation. Our model does not reproduce the CPD experimental profiles, indicating that additional studies are necessary. Overall, our findings advance the understanding of isoprene pyrolysis and its related atmospheric pollutants in naturally occurring vegetation fires where smoldering and oxygen-deficient combustion processes are present.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Microstructure homogenization of concrete used in nuclear power plants

Almost all nuclear power plants in the United States are operating past their intended lifetimes or are requesting lifetime extensions. Therefore, understanding changes to the concrete containment structure over time is crucial to evaluate the structure’s continued viability. Concrete materials are heterogeneous particulate composites that exhibit viscoelastic material properties, which can lead to slow deformation over time, causing stress redistribution and the potential for creep cracking. A code to generate random, three dimensional (3D) concrete microstructures has been developed and paired with finite element analysis to predict the long-term viscoelastic properties of concrete. Data from these simulations are used to develop constitutive equations for the viscoelastic behavior of the homogenized concrete. The codes in this work are used to virtualize laboratory experiments, to obtain long-term creep data in a faster, cheaper manner. To validate this work, the simulated creep behavior of concrete is compared to 800 d of experimental data that has been extended to 27 y of data using the Time-Temperature superposition (TTS) principal. Excellent agreement between the simulation results and experimental data is seen.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Experimental workflow to estimate model parameters for evaluating long term viscoelastic response of CO2 storage caprocks

Understanding the time-dependent behavior of reservoir and sealing formations is critical to assessing risks associated with geological carbon storage since time-dependent deformation strongly influences mechanical responses of some rock types. Many studies have evaluated the risk of CO2 leakage and induced seismicity by assuming poroelastic rheology in sealing formations. Few have considered viscoelastic or other time-dependent responses, where the existing literature adopts 1D models to represent long-term time-dependent responses. This is primarily because to date, the general form of a reasonable 3D time-dependent model for rocks remains unclear. In this paper, we address this unclear issue by proposing a new workflow to select constitutive modeling parameters to evaluate if a 3D viscoelastic model is reasonable using several-hour-long experimental data and a power-law response to extrapolate to the decades-long time frames of interest in geologic carbon storage. To provide experimental data, we conducted multi-level loading/unloading triaxial relaxation tests with four rock types. The experimental results showed that the maximum load relaxation observed is approximately 49%, with some rock types showing as little as 1.4%. Using a simple linear viscoelastic model, parameters were chosen such that a maximum deviation of 1.5 MPa in axial stress and 7 MPa in radial stress was attained with the extrapolated 30-year data. We found that a reasonable parameter range for the normalized elastic modulus is 0.1~2 for rocks with significant time-dependent responses and 0.01~0.06 for those with small time-dependent responses. No matter how significant time-dependent responses are for rocks considered, our results showed that the relaxation time has a general range of 1~10^10 s, whose time scale can be one or two orders higher than a time frame typically envisioned for CO2 injection projects.

Stress relaxation, 3D time-dependent model, viscoe↗

DOE OSTI2025

Utility of a hybrid approach to the hadronic vacuum polarization contribution to the muon anomalous magnetic moment

An accurate calculation of the leading-order hadronic vacuum polarization (LOHVP) contribution to the anomalous magnetic moment of the muon ( a μ ) is key to determining whether a discrepancy, suggesting new physics, exists between the Standard Model and experimental results. This calculation can be expressed as an integral over Euclidean time of a current-current correlator G ( t ) , where G ( t ) can be calculated using lattice QCD or, with dispersion relations, from experimental data for e + e − → hadrons . The BMW/DMZ collaboration recently presented a hybrid approach in which G ( t ) is calculated using lattice QCD for most of the contributing t range, but using experimental data for the largest t (lowest energy) region. Here we study the advantages of varying the position t = t 1 separating lattice QCD from data-driven contributions. The total LOHVP contribution should be independent of t 1 , providing both a test of the experimental input and the robustness of the hybrid approach. We use this criterion and a correlated fit to show that Fermilab/HPQCD/MILC lattice QCD results from 2019 strongly favor the CMD-3 cross section data for e + e − → π + π − over a combination of earlier experimental results for this channel. Further, the resulting total LOHVP contribution obtained is consistent with the result obtained by BMW/DMZ, and supports the scenario in which there is no significant discrepancy between the experimental value for a μ and that expected in the Standard Model. We then discuss how improved lattice results in this hybrid approach could provide a more accurate total LOHVP across a wider range of t 1 values with an uncertainty that is smaller than that from either lattice QCD or data-driven approaches on their own. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Inference of strength and phase transition kinetics in dynamically-compressed tin

Dynamic compression experiments in condensed matter are of interest in part because they provide opportunities to examine material response under extreme conditions; however, the inference of material behavior from dynamic experiments is challenging in the presence of phase transitions exhibiting kinetic processes. Demonstrating an approach to quantitative interpretation of such dynamic experiments, here we present a Bayesian model calibration of strength and phase transformation parameters to data drawn from pulsed power and gas gun shot experiments. The posterior predictions of the Bayesian model capture the experimental measurements and account for the various uncertainties in the experimental configurations. This holistic approach to model calibration utilizing multiple types of experimental data identifies important cross correlations among kinetics, strength, and the phase boundary. Improved insight into potential sources of current model form error is provided by comparing the differences between calibrations against different subsets of the experimental data.

42 ENGINEERING↗

DOE OSTI2024

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗

DOE OSTI2023

Templates of expected measurement uncertainties for average prompt and total fission neutron multiplicities

In this paper, we provide templates of measurement uncertainty sources expected to appear for average prompt- and total-fission neutron multiplicities, $ \overline\nu_p $ and $ \overline\nu_t $, for the following measurement types: absolute manganese-bath experiments for $ \overline\nu_t $, absolute and ratio liquid-scintillator measurements for $ \overline\nu_p $. These templates also suggest a typical range of these uncertainties and their correlations based on a survey of available experimental data, associated literature, and feedback from experimentalists. In addition, the information needed to faithfully include the associated experimental data into the nuclear-data evaluation process is provided.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

239 Pu(n,f) Neutron Multiplicity Evaluation with CGMF-Very First Release Candidate [Slides]

239 Pu(n,f) nu-bar evaluated with CGMF from E inc = 0.1-20 MeV. Changes compared to VIII.0: Prior: CGMF model included via Kalman and sensitivities of CGMF model parameters to nu-bar. Evaluation technique: Kalman including correction for PPP. Experimental data: Nearly all data that Phil took into account (I rejected: Huanqiao, Johnstone, Leroy, Nesterov, Smirenkin); New UQ for all experimental data; Marini; No correlations between unc. of different exp., except for Cf-252(sf) nu-bar uncertainty cross-correlating all uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

239 Pu(n,f) Neutron Multiplicity Evaluation with CGMF-Very First Release Candidate. [Slides]

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Elastic scattering in the 12 N + 197 Au system at 𝐸 lab = 70 MeV

Halo structure is an interesting exotic configuration developed in some light weakly bound nuclei, where a valence particle orbits a nuclear core. Signatures of halo structure can be observed in the angular distributions of the elastic scattering induced by these nuclei at energies around the Coulomb barrier. There are some well-studied reactions with neutron-rich halo nuclei, such as 6 He and 11 Li . However, the information is scarce on the proton-rich side. Recent works confirm the halo structure in the 8B nuclei but still lack more experimental studies for other proton-halo candidates, such as 12 N and 17 F. In this work, we report experimental data for the elastic scattering of 12 N on 197 Au target at 𝐸 lab = 70 MeV. The 12 N is a proton-rich nucleus with proton separation energy 𝑆 𝑝 = 600 keV, which is higher than the 8 B (𝑆 𝑝 = 137 keV) and almost the same as the 17 F (𝑆 𝑝 = 601 keV). Data were obtained at the Cyclotron Institute of Texas A&M University where the 12 N radioactive beam was produced by the momentum achromatic recoil spectrometer. The calculation of the optical model was used to fit the measured angular distribution for the elastic scattering and to obtain the reaction cross section 𝜎 𝑅 . We also performed continuum discretized coupled-channel calculations to compare with the experimental data. The angular distribution of the elastic cross sections exhibits a suppression of the Fresnel peak. From the fitting of the optical model, we obtain the total reaction cross section, 𝜎 𝑅 = 1269 ± 41 mb. The agreement between the CDCC calculation and the experimental elastic scattering data is limited and the breakdown does not exhaust the measured 𝜎 𝑅 . The resulting reduced reaction cross section 𝜎 Red for 12 N + 197 Au is large and comparable to the one obtained for the 8 B + 208 Pb system. This suggests a strong decoupling of the valence proton from the core because of the low binding energy and a dynamic polarization effect. Here, more research is required to estimate the contribution of core excitation.

Elastic scattering reactions↗

DOE OSTI2023

Developing and applying quantifiable metrics for diagnostic and experiment design on Z

This project applies methods in Bayesian inference and modern statistical methods to quantify the value of new experimental data, in the form of new or modified diagnostic configurations and/or experiment designs. We demonstrate experiment design methods that can be used to identify the highest priority diagnostic improvements or experimental data to obtain in order to reduce uncertainties on critical inferred experimental quantities and select the best course of action to distinguish between competing physical models. Bayesian statistics and information theory provide the foundation for developing the necessary metrics, using two high impact experimental platforms on Z as exemplars to develop and illustrate the technique. We emphasize that the general methodology is extensible to new diagnostics (provided synthetic models are available), as well as additional platforms. We also discuss initial scoping of additional applications that began development in the last year of this LDRD.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Discovering equations that govern experimental materials stability under environmental stress using scientific machine learning

Abstract While machine learning (ML) in experimental research has demonstrated impressive predictive capabilities, extracting fungible knowledge representations from experimental data remains an elusive task. In this manuscript, we use ML to infer the underlying differential equation (DE) from experimental data of degrading organic-inorganic methylammonium lead iodide (MAPI) perovskite thin films under environmental stressors (elevated temperature, humidity, and light). Using a sparse regression algorithm, we find that the underlying DE governing MAPI degradation across a broad temperature range of 35 to 85 °C is described minimally by a second-order polynomial. This DE corresponds to the Verhulst logistic function, which describes reaction kinetics analogous to self-propagating reactions. We examine the robustness of our conclusions to experimental variance and Gaussian noise and describe the experimental limits within which this methodology can be applied. Our study highlights the promise and challenges associated with ML-aided scientific discovery by demonstrating its application in experimental chemical and materials systems.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Predicting transport of intra-articularly injected growth factor fusion proteins into human knee joint cartilage

There are no drugs or treatment methods known to prevent the development of post-traumatic osteoarthritis (PTOA), a type of osteoarthritis (OA) that is triggered by traumatic joint injuries and accounts for ~12% of the nearly 600 million OA cases worldwide. Lack of effective drug delivery techniques remains a major challenge in developing clinically effective treatments, but cationic delivery carriers can help overcome this challenge. Scaling up treatments that are effective in in vitro models to achieve success in preclinical in vivo models and clinical trials is also a challenging problem in the field. Here we use a cationic green fluorescent protein (GFP) as a carrier to deliver Insulin-Like Growth Factor 1 (IGF-1), a drug considered as a potential therapeutic for PTOA. GFP-IGF-1 conjugates were first synthesized as fusion proteins with different polypeptide linkers, and their transport properties were characterized in human cartilage explants. In vitro experimental data were used to develop a predictive mathematical transport model that was validated using an independent in vitro experimental data set. In this study, the model was used to predict the transport of these fusion proteins upon intra-articular injection into human knee joints. The predictions included results for the rate and extent of fusion protein penetration into cartilage, and the maximum levels of fusion proteins that would escape into systemic circulation through the joint capsule. Together, our transport measurements and model set the stage for translation of such explant culture studies to in vivo preclinical studies and potentially clinical application.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Elucidating the temperature and density dependence of silver chloride hydration numbers in high-temperature water vapor: A first-principles molecular simulation study

Hydration numbers of metal complexes in low-density aqueous solutions are required for developing geochemical models for ore-forming metals and for designing supercritical desalination processes. In this study, we investigate the temperature and density dependence of the hydration numbers of silver chloride at 623 K, 673 K, and 713 K and densities of 10–100 kg/m 3 . Experimental estimates of the hydration number in the literature for AgCl at these conditions are inconclusive and possibly contradictory as to the temperature and density dependence. First-principles molecular simulation presents an attractive alternative to experimental measurements. Specifically, recent work shows that machine-learning-accelerated nested Monte Carlo simulations provide reliable estimates for the hydration numbers of CuCl at 623 K from 10 to 100 kg/m 3 . Using the same technique, we find a monotonic temperature dependence, with the hydration number decreasing slightly with increasing temperature. In addition, the simulation-predicted hydration numbers steadily increase with increasing density. These temperature and density trends are in agreement with certain experimental data sets. Therefore, this work demonstrates how first-principles Monte Carlo simulations assist in resolving discrepancies between experimental data sets. Our simulation results also correctly predict that the hydration number, and thus also the solubility, of AgCl is lower than that of CuCl under the same conditions. Furthermore, the bond length and angle formed between the water complex and AgCl differ from those for CuCl, consistent with the lower solubility of AgCl.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Impact of Surface Recombination on Single-Event Charge Collection in an SOI Technology

Semiconductor–insulator interfaces play an important role in the reliability of integrated devices; however, the impact of these interfaces on the physical mechanisms related to single-event effects has not been previously reported. We present experimental data that demonstrate that single-event charge collection can be impacted by changes in interface quality. The experimental data, combined with simulations, show that single-event response may depend on surface recombination at interface defects. The effect depends on strike location and increases with increasing linear energy transfer (LET). Surface recombination can affect single-event charge collection for interfaces with a surface recombination velocity (SRV) of 1000 cm/s and is a dominant charge collection mechanism with SRV >10 5 cm/s.

42 ENGINEERING↗

DOE OSTI2023

Large Eddy Simulation of Gasoline Sprays in a Lagrangian–Eulerian Framework Using the High-Order Spectral Element Method

Predicting the spray evolution using simulations requires accurate modeling of the turbulent gas-phase flow field. Here, in this study, the high-order spectral-element method (SEM), implemented in the code Nek5000, was used to provide highly resolved solutions to the turbulent flow field. Spray modeling capabilities were implemented into the Nek5000 code. The spray is modeled in a Lagrangian–Eulerian (LE) framework, where the liquid is represented by discrete parcels of droplets. The method for coupling liquid and gas in the context of SEM is described, which allows for very fine meshes to be used without affecting the stability of the solution. Large-eddy simulations (LES) of the eight-hole ECN Spray G gasoline injector were conducted. Numerical results are compared against experimental data for liquid penetration, droplet size and gas velocity. The morphology of the multiplume spray is compared against experimental data. The effect of different spray injection inputs is analyzed. It was found that using a plume direction of 33 deg and an injection cone angle of 30 deg produced the best results overall. This work shows the applicability of SEM for spray modeling applications, where use of a high-order flow solver can help us understand the multiplume spray aerodynamics and how it leads to plume collapse under certain conditions. Results also highlight the need for tuning spray input parameters in the LE framework, even when high-fidelity gas flow solutions are possible.

33 ADVANCED PROPULSION SYSTEMS↗