Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter sampling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Bayesian analysis of (3 +1)⁢D relativistic nuclear dynamics with the RHIC beam energy scan data

This work presents a Bayesian inference study for relativistic heavy-ion collisions in the beam energy scan program at the BNL Relativistic Heavy-Ion Collider. The theoretical model simulates event-by-event (3+1)-dimensional [(3+1)⁢D] collision dynamics using hydrodynamics and hadronic transport theory. We analyze the model's 20-dimensional posterior distributions obtained using three model emulators with different accuracy and demonstrate the essential role of training an accurate model emulator in the Bayesian analysis. Our analysis provides robust constraints on the quark-gluon plasma's transport properties and various aspects of (3+1)⁢D relativistic nuclear dynamics. By running full model simulations with 100 parameter sets sampled from the posterior distribution, we make predictions for p T -differential observables and estimate their systematic theory uncertainty. Here, a sensitivity analysis is performed to elucidate how individual experimental observables respond to different model parameters, providing useful physics insights into the phenomenological model for heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Constraining Cosmology with Simulation-based inference and Optical Galaxy Cluster Abundance

We test the robustness of simulation-based inference (SBI) in the context of cosmological parameter estimation from galaxy cluster counts and masses in simulated optical datasets. We construct ``simulations'' using analytical models for the galaxy cluster halo mass function (HMF) and for the observed richness (number of observed member galaxies) to train and test the SBI method. We compare the SBI parameter posterior samples to those from an MCMC analysis that uses the same analytical models to construct predictions of the observed data vector. The two methods exhibit comparable performance, with reliable constraints derived for the primary cosmological parameters, ($\Omega_m$ and $\sigma_8$), and richness-mass relation parameters. We also perform out-of-domain tests with observables constructed from galaxy cluster-sized halos in the Quijote simulations. Again, the SBI and MCMC results have comparable posteriors, with similar uncertainties and biases. Unsurprisingly, upon evaluating the SBI method on thousands of simulated data vectors that span the parameter space, SBI exhibits worsened posterior calibration metrics in the out-of-domain application. We note that such calibration tests with MCMC is less computationally feasible and highlight the potential use of SBI to stress-test limitations of analytical models, such as in the use for constructing models for inference with MCMC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Scan strategies in EBM-printed IN718 and the physics of bulk 3D microstructure development

Three-dimensional (3D) characterization provides opportunities for understanding processing-structure relationships in additively manufactured (AM) materials. Bulk samples of Inconel 718 were fabricated via electron beam melting (EBM) in order to study microstructural development as a function of energy input and beam scan strategy. TriBeam tomography of bulk Inconel 718 microstructures built under steady-state growth conditions reveals the sensitivity of microstructure formation and evolution to machine process parameters. Here, samples manufactured using a narrow range of energy input per unit build area result in varied grain morphologies and crystallographic textures. Using TRUCHAS, a thermal simulation software, the thermal history of bulk scan strategies was predicted, and combined with a calibrated microstructure-processing map to accurately predict bulk grain morphologies. The solidification parameters and the 3D measured nucleation density are used to predict the transition between columnar and equiaxed grain morphologies, providing a process map to guide AM parameter choices to locally control as-printed microstructure. A two-dimensional metric for characterizing bulk grain morphology was also found to agree well with predictions from the process map calibrated by 3D data. Combined with 3D tomography and thermal modelling, the physics of structure development were understood at a new level of detail with respect to the competing processes of grain nucleation and epitaxial growth.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dark Energy Survey Year 3 Results: Constraints on cosmological parameters and galaxy bias models from galaxy clustering and galaxy-galaxy lensing using the redMaGiC sample

We constrain cosmological parameters and galaxy-bias parameters using the combination of galaxy clustering and galaxy-galaxy lensing measurements from the Dark Energy Survey Year-3 data. We describe our modeling framework and choice of scales analyzed, validating their robustness to theoretical uncertainties in small-scale clustering by analyzing simulated data. Using a linear galaxy bias model and redMaGiC galaxy sample, we obtain constraints on the matter content of the universe to be $\Omega_{\rm m} = 0.325^{+0.033}_{-0.034}$. We also implement a non-linear galaxy bias model to probe smaller scales that includes parameterizations based on hybrid perturbation theory, and find that it leads to a 17% gain in cosmological constraining power. Using the redMaGiC galaxy sample as foreground lens galaxies, we find the galaxy clustering and galaxy-galaxy lensing measurements to exhibit significant signals akin to decorrelation between galaxies and mass on large scales, which is not expected in any current models. This likely systematic measurement error biases our constraints on galaxy bias and the $S_8$ parameter. We find that a scale-, redshift- and sky-area-independent phenomenological decorrelation parameter can effectively capture this inconsistency between the galaxy clustering and galaxy-galaxy lensing. We perform robustness tests of our methodology pipeline and demonstrate stability of the constraints to changes in the theory model. After accounting for this decorrelation, we infer the constraints on the mean host halo mass of the redMaGiC galaxies from the large-scale bias constraints, finding the galaxies occupy halos of mass approximately $1.5 \times 10^{13} M_{\odot}/h$.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

Evaluation of GIS Technology in Assessing and Modeling Land Management Practices

There is an increasing concern of land owners to protect and maintain healthy and sustainable agroecosystems through the implementation of best management practices (BMP). The objectives of this study were: (1) To develop and evaluate the use of a Geographic Information System (GIS) technology for enhancing field-scale management practices; (2) evaluate the use of 2-dimensional displays of the landscape and (3) define spatial classes of variables from interpretation of geostatistical parameters. Soil samples were collected to a depth of 2 m at 15 cm increments. Existing data from topographic, land use, and soil survey maps of the Winfred Thomas Agricultural Research Station were converted to digital format. Additional soils data which included texture, pH, and organic matter were also generated. The digitized parameters were used to create a multilayered field-scale GIS. Two dimensional (2-D) displays of the parameters were generated using the ARC/INFO software. The spatial distribution of the parameters evaluated in both fields were similar which could be attributed to the similarity in vegetation and surface elevation. The ratio of the nugget to total semivariance, expressed as a percentage, was used to assess the degree of spatial variability. The results indicated that most of the parameters were moderate spatially dependent Biophysical constraint maps were generated from the database layers, and used in multiple combination to visualize results of the BMP. Understanding the spatial relationships of physical and chemical parameters that exists within a field should enable land managers to more effectively implement BMP to ensure a safe and sustainable environment.

Archer, F.↗
NASA NTRSDate not supplied

Sampling Functions from Gaussian Processes and Structured Covariance Gaussian Networks

When learning aerodynamic models from data, it is critical to incorporate estimates of model uncertainty. This motivates the design of probabilistic aerodynamic databases which can be sampled to generate physically and statistically plausible aerodynamic models. In this talk we discuss how to sample deterministic functions from two different kinds of probabilistic models and demonstrate their use. First, Gaussian Process Regressors (GPRs) are a widely used probabilistic kernel-based model which can be thought of as Gaussian distributions over functions. GPRs are generally trained by maximizing the marginal likelihood of seeing the training data over the kernel parameter space. Sample functions are easily generated by drawing points from the Gaussian distribution at desired input points. However, when the points are not known ahead of time, the classical sampling approach is not possible since successive function samples will generate different function realizations. We present an approach for sampling consistent function evaluations from a GPR over multiple samples. Second, we describe a neural network architecture which learns a conditional Gaussian distribution by maximizing the marginal likelihood at each point in the input space. We then discuss and compare several options for generating sample functions which match this distribution. Finally, we demonstrate the use of these probabilistic aerodynamic models in an atmospheric reentry simulation.

Gaussian process regression↗
DOE OSTI2020

Radiological Monitoring Plan for the Oak Ridge Y-12 National Security Complex: Surface Water

DOE Order 458.1 requires that dose estimates consider contributions from all facilities. In the Y-12 Radiological Monitoring Plan (RMP), surface water is monitored at points that reflect individual facilities, as well as at points that reflect the combined contributions of all facilities. This monitoring plan does not consider other potential routes (i.e., airborne releases and food chains). Thus, a complete determination of total effective dose (TED) cannot be made based on this plan alone. The other routes from Y-12, and all routes from other DOE facilities on the Oak Ridge Reservation (e.g., Oak Ridge National Laboratory (ORNL) and The Heritage Center), must be considered in order to satisfy DOE Order 458.1 requirements. Determination of TED from all sites and pathways is done through the use of dose-assessment models and is documented in the Annual Site Environmental Report. This monitoring plan provides adequate monitoring goals for Y-12 surface water releases to provide input of sufficient sensitivity and accuracy to reliably determine the Y-12 surface water component of the TED. The routine radiological monitoring program is designed to monitor effluents at four types of locations: (1) treatment facilities, (2) other point and area source discharges, (3) instream locations, and (4) production building roof run-off. With this sampling and analysis program, data will be obtained on primary point sources as well as on locations that represent the composite of other potential sources. This plan will be reviewed periodically to determine necessary modifications to the sampling frequencies, parameters, and locations. Modifications, if any, will be based on the analysis of the previous data and its effectiveness in satisfying the objectives of this plan. Appendix A contains graphs of the sum of the DCS fractions for locations and frequencies contained in a previous version of this plan. The data was collected from January 2009 through December 2019. Each sample was analyzed, and each result was divided by the appropriate DCS to compute a DCS fraction. These fractions were summed for all isotopes. According to DOE –STD-1196-2011, the annual average of these sums should be below 1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Model-independent measurement of 𝐷 0 −$\overline{D}$ 0 mixing parameters in 𝐷 0 → 𝐾$^{0}_{𝑆}$⁢𝜋 + ⁢𝜋 − decays at Belle and Belle II

We perform a model-independent measurement of the 𝐷 0 −$\overline{D}$ 0 mixing parameters using samples of 𝑒 + ⁢𝑒 − -collision data collected by the Belle and Belle II experiments that have integrated luminosities of 951 fb −1 and 408 fb −1 , respectively. Approximately 2.05 × 10 6 neutral 𝐷 mesons are reconstructed in the 𝐷 0 → 𝐾$^{0}_{S}$⁢𝜋 + ⁢𝜋 − channel, with the neutral 𝐷 flavor tagged by the charge of the pion in the 𝐷* + → 𝐷 0⁢ 𝜋 + decay. Assuming charge-parity symmetry, the mixing parameters are measured to be 𝑥 = (4.0 ±1.7 ±0.4) × 10 −3 and 𝑦 = (2.9 ±1.4 ±0.3) × 10 −3 , where the first uncertainties are statistical and the second systematic. The results are consistent with previous determinations.

Charmed mesons↗
NASA NTRS2008

Joint Bayesian Component Separation and CMB Power Spectrum Estimation

We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are (1) conditional sampling of foreground spectral parameters and (2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global template amplitudes) given spectral parameters. Given a parametric model of the foreground signals, we estimate efficiently and accurately the exact joint foreground- CMB posterior distribution and, therefore, all marginal distributions such as the CMB power spectrum or foreground spectral index posteriors. The main limitation of the current implementation is the requirement of identical beam responses at all frequencies, which restricts the analysis to the lowest resolution of a given experiment. We outline a future generalization to multiresolution observations. To verify the method, we analyze simple models and compare the results to analytical predictions. We then analyze a realistic simulation with properties similar to the 3 yr WMAP data, downgraded to a common resolution of 3 deg FWHM. The results from the actual 3 yr WMAP temperature analysis are presented in a companion Letter.

numerical↗
DOE OSTI2025

Decayheatml

This code is designed to predict and analyze the decay heat generated in molten salt reactors (MSRs) using a hybrid approach that combines machine learning and segmented polynomial fitting. The accurate prediction of decay heat is essential for reactor safety and the optimization of spent fuel storage. The code operates through several key components: 1) Data Architecture: It incorporates a modular data architecture that handles various MSR-specific operational parameters such as power density, humidity content, and air ingress. These parameters are sampled using Sobol sequences to ensure comprehensive coverage of operational uncertainties. 2) Machine Learning Framework: The code employs a diverse set of machine learning models, including polynomial regression, decision trees, random forests, gradient boosting, support vector regression, k-nearest neighbors, multi-layer perceptrons, and symbolic regression. These models are trained to predict decay heat over a wide temporal range, from immediate shutdown up to 10,000 years. 3) Region-Optimized Training: The temporal domain is divided into multiple regions, each modeled separately to capture distinct decay heat characteristics across different time scales. This approach significantly improves the accuracy and interpretability of predictions. 4) Segmented Polynomial Interpretation (SPI): The SPI method translates machine learning predictions into piecewise polynomial equations. These equations are physically interpretable and can be directly integrated into existing engineering workflows and safety analyses. 5) Front-End Interfaces: The code includes both a Jupyter notebook interface for research development and a Streamlit web application for operational deployment. These interfaces allow users to interactively explore decay heat predictions, adjust operational parameters, and visualize results in real-time. 6) Applications: The framework supports various applications, including safety system validation and spent fuel container optimization. It enables real-time evaluation of worst-case decay heat scenarios, informing the design of passive safety systems and optimizing container designs for long-term storage. Overall, this code provides a robust, accurate, and user-friendly tool for predicting decay heat in MSRs, enhancing reactor safety, and optimizing spent fuel management.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2022

Improvement of two-phase closure models in CTF using Bayesian inference

Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, extensive capabilities have been developed in CTF to analyze light-water reactors (LWRs) for normal operating conditions, departure from nucleate boiling (DNB), and system transients. However, further improvements are required in the modeling and simulation of boiling water reactors (BWRs), which is a focus of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, CTF validation results were used to optimize selected modeling coefficients by calibrating to experimental data using a Bayesian inference approach. Here, calibration studies were conducted to improve (vapor) void fraction prediction without worsening the two-phase pressure drop prediction, as well as to improve the two-phase pressure drop prediction. Calibration was performed for interfacial drag and wall shear models. Surrogates were developed to alleviate the computational expense required for sampling the parameter space using Markov chain Monte Carlo (MCMC). An assessment performed with calibrated models demonstrated an improvement of CTF in its prediction of key parameters such as void fraction and two-phase pressure drop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Quantifying the propagation of parametric uncertainty on flux balance analysis

Flux balance analysis (FBA) and associated techniques operating on stoichiometric genome-scale metabolic models play a central role in quantifying metabolic flows and constraining feasible phenotypes. At the heart of these methods lie two important assumptions: (i) the biomass precursors and energy requirements neither change in response to growth conditions nor environmental/genetic perturbations, and (ii) metabolite production and consumption rates are equal at all times (i.e., steady-state). Despite the stringency of these two assumptions, FBA has been shown to be surprisingly robust at predicting cellular phenotypes. In this paper, we formally assess the impact of these two assumptions on FBA results by quantifying how uncertainty in biomass reaction coefficients, and departures from steady-state due to temporal fluctuations could propagate to FBA results. In the first case, conditional sampling of parameter space is required to re-weigh the biomass reaction so as the molecular weight remains equal to 1 g mmol –1 , and in the second case, metabolite (and elemental) pool conservation must be imposed under temporally varying conditions. Results confirm the importance of enforcing the aforementioned constraints and explain the robustness of FBA biomass yield predictions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

BeyondPlanck X. Bandpass and beam leakage corrections

We discuss the treatment of bandpass and beam leakage corrections in the Bayesian BeyondPlanck CMB analysis pipeline as applied to the Planck LFI measurements. As a preparatory step, we first apply three corrections to the nominal LFI bandpass profiles including removal of a known systematic effect in the ground measuring equipment at 61 GHz; smoothing of standing wave ripples; and edge regularization. The main net impact of these modifications is an overall shift in the 70 GHz bandpass of +0.6 GHz; we argue that any analysis of LFI data products, either from Planck or BeyondPlanck, should use these new bandpasses. In addition, we fit a single free bandpass parameter for each radiometer of the form Δ i = Δ 0 + $δ$ i , where Δ 0 represents an absolute frequency shift per frequency band and $δ$i is a relative shift per detector. The absolute correction is only fitted at 30 GHz with a full $\chi$ 2 -based likelihood, resulting in a correction of Δ 30 = 0.24 ± 0.03 GHz. The relative corrections are fitted using a spurious map approach, fundamentally similar to the method pioneered by the WMAP team, but without introducing many additional degrees of freedom. All bandpass parameters are sampled using a standard Metropolis sampler within the main BeyondPlanck Gibbs chain, and bandpass uncertainties are thus propagated to all other data products in the analysis. In total, we find that our bandpass model significantly reduces leakage effects. For beam leakage corrections, we adopt the official Planck LFI beam estimates without additional degrees of freedom, and only marginalize over the underlying sky model. We note that this is the first-time leakage from beam mismatch has been included for Planck LFI maps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A comparison of Raman and pyrometry dynamic temperature measurements of shocked cyclohexane

Optical pyrometry is a well-developed, broadly applicable method to determine temperature in dynamic compression experiments. However, it measures radiation from only the skin depth of the radiant material, is highly susceptible to spurious backgrounds, and requires an often-unmeasurable free parameter—the sample dynamic emissivity. Raman spectroscopy offers a method to determine temperature directly from fundamental theory, but it is practically limited in applications to a small subset of shock experiments. In this work, we describe an experiment allowing simultaneous application of Raman and pyrometric temperature techniques to benchmark a specific instance of pyrometry. We attempted both measurements on multiple experiments and compare the successful temperature results between the separate but highly repeatable shots. Here, we find that dynamic shock temperatures inferred from Raman spectroscopy and pyrometry generally agree within the uncertainties of each measurement, although those uncertainties are larger than typical for a specially designed pyrometry experiment on a shocked metal sample.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

A Time-Domain Protection Approach for AC Transmission Systems With Grid-Forming Resources

Ac transmission protection must reliably detect, classify, and locate short-circuit faults from voltage and current measurements. At present, these functionalities, which have been classically engineered using phasors approaches, are being challenged by the dynamic behavior and fault-current limits of converter-based generation. This paper tackles these challenges by engineering a time-domain protection approach that leverages the classical Bergeron model in a new manner. Low- and high-impedance faults are detected and classified by ascertaining how well line voltage and current measurements match the Bergeron equations. Faults are located by posing a novel one-variable optimization problem, whereas voltage and current waveforms at the fault location are estimated by unveiling rigorous relationships. The proposed elements are secure against external faults, measurement errors, and variation of line parameters and sampling time. Furthermore, these advances are tested via electromagnetic transient simulations and are significant to satisfy IEEE and North American Electric Reliability Corporation requirements.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1973

Moessbauer spectroscopic studies of iron-doped rutile.

Moessbauer spectra were obtained of single crystal and powdered samples of rutile (TiO2) doped to about one percent by weight in isotopically enriched iron. It is shown that the oxidation state may be reversibly altered in situ. After reduction, the oxygen neighbors of the dopant ion are apparently shifted to accomodate the larger ferrous ion. The agreement of calculated quadrupole splittings with experimental results suggests that impurities and oxygen vacancies are uniformly distributed in powdered samples, giving the dopant ions an 'ideal lattice' local environment. The differences between the single crystal and powder sample hyperfine parameters are attributed to variations in stoichiometry, charge compensation mechanisms, or other diffusion related parameters.

Stampfl, P. P.↗
NASA NTRS1977

Flexible high-speed instrumentation system

Remote multiplexer/demultiplexer digital instrumentation system, which is suitable to both airborne and ground-based data acquisition/process control applications, can be employed in a variety of research and flight test applications where great flexibility to accommodate changes in number of parameters, data sampling rates, and signal conditioning is required.

Bartoli, F.↗