Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “input estimation”

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 1,693 records · Page 94

Southern Idaho Health & Air Quality II: Evaluating Atmospheric Mixing Height Estimations in the Western United States

Wildfires in the western United States have caused immense infrastructure damage and loss of human life in recent years. Wildfire smoke, which travels far from its original source, is also harmful to human health. Mixing height, which acts as a lid and prevents smoke from rising above a certain altitude in the lower troposphere, is a critical input in smoke dispersion and air quality models used by agencies that monitor wildfires. These models, coupled with forecaster expertise, are also used to decide when it is safe to execute a prescribed burn. The DEVELOPID team partnered with the National Weather Service (NWS) Fire Weather Program, Bureau of Land Management (BLM), National Interagency Fire Center (NIFC), and National Park Service (NPS) Fire Management Program Center (FMPC) to help improve reliability and confidence in mixing height estimations, and therefore the burn prescription decision-making process. To that end, the team developed a toolbox for measuring smoke-related aerosol mixing heights using Cloud-Aerosol LiDAR and Infrared Pathfinder Satellite Observations (CALIPSO) Vertical Feature Mask granules. CALIPSO mixing heights and NWS estimations covaried significantly and positively. However, substantial disagreement between the methods stymied the team’s attempts to quantify systematic bias in a meaningful way. The relative error between the methods was especially large at low mixing heights, which suggests that this method of validation may only be suitable at higher altitudes.

Dean Berkowitz↗

Polaris-PARCS Sensitivity Study on LWR Fuel Cycles: Polaris Input Options

This study is the first of a multi-phase effort to assess the sensitivity of light-water reactor (LWR) core-level prediction biases to changes in lattice-level calculation parameters. Prediction bias is the measured-to-predicted difference in a core-level quantity of interest (QOI) which can be estimated by comparing the simulation results with the plant-measured data for key nuclear parameters. The LWR two-step neutronics codes employed herein are the SCALE–Polaris lattice physics code (v6.3.1) and the Purdue Advanced Reactor Core Simulator (PARCS) nodal diffusion simulator (v3.4.2), both funded and used for confirmatory analysis to support licensing by the US Nuclear Regulatory Commission (NRC). Polaris–PARCS is used to model Watts Bar Unit 1 cycles 1–3 and Peach Bottom Unit 2 cycles 1–3. This study focuses on the impact of changes to Polaris input options such as scattering treatment or quadrature settings and how these input options induce changes in core-level quantities of interest (QOIs)bias. The report documents multiple bias assessments for different modeling choices and compares the bias magnitude to the QOI measurement uncertainties. Future companion reports will investigate the sensitivity of core-level LWR prediction bias to Polaris input options and Polaris-computed QOIs such as few-group assembly-homogenized cross sections to gain an understanding of the key drivers of prediction bias at lattice and core levels for application of a two-step LWR neutronics procedure in a licensing scenario.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Automated reaction kinetics and network exploration (Arkane): A statistical mechanics, thermodynamics, transition state theory, and master equation software

The open-source statistical mechanics software described here, Arkane–Automated Reaction Kinetics and Network Exploration–facilitates computations of thermodynamic properties of chemical species, high-pressure limit reaction rate coefficients, and pressure-dependent rate coefficient over multi-well molecular potential energy surfaces (PES) including the effects of collisional energy transfer on phenomenological kinetics. Arkane can use estimates to fill in information for molecules or reactions where quantum chemistry information is missing. The software solves the internal energy master equation for complex unimolecular reaction systems. Inputs to the software include converged electronic structure computations performed by the user using a variety of supported software packages (Gaussian, Molpro, Orca, TeraChem, Q-Chem, Psi4). The software outputs high-pressure limit rate coefficients and pressure-dependent phenomenological rate coefficients, as well as computed thermodynamic properties (enthalpy, entropy, and constant pressure heat capacity) with added energy corrections. Some of the key features of Arkane include treatment of 1D, 2D or ND hindered internal rotation modes, treatment of free internal rotation modes, quantum tunneling effect consideration, transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficient computations, master equation solution with four implemented methods, inverse-Laplace transform of high-pressure limit rate coefficients into the energy domain, energy corrections based on bond-additivity or isodesmic reactions, automated and efficient PES exploration, and PES sensitivity analysis. The present work describes the design of Arkane, how it should be used, and refers to the theory that it employs. Arkane is distributed via the RMG-Py software suite (https://github.com/ReactionMechanismGenerator/RMG-Py).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Techno-economic prospects for producing Fischer-Tropsch jet fuel and electricity from lignite and woody biomass with CO 2 capture for EOR

This work explores the prospective techno-economic performance of facilities that produce low- and net-negative-carbon liquid transportation fuels and electricity with CO 2 capture for enhanced oil recovery. The lignite and biomass-to-jet fuel process is based on KBR’s TRIG gasifier, Rectisol (for sulfur removal and CO 2 capture), fixed-bed low temperature Fischer-Tropsch synthesis of liquid fuels, and Brayton/Rankine combined cycles to convert synthesis/refining off-gases and waste heat to electricity. This work leverages a recent, highly-detailed assessment of a prospective first-of-a-kind (FOAK) demonstration facility to develop highly detailed Aspen Plus process simulations for nine prospective N th -of-a-kind (NOAK) plant equipment configurations. Component-level capital costs from the FOAK study are scaled and adjusted to reflective prospective learning-by-doing to estimate capital costs for the NOAK designs. NOAK plant economic performance is found to be largely insensitive to variations in plant configurations and electricity output fraction, but biomass input fraction significantly affects profitability. Facilities that consume only carbon–neutral biomass, with no lignite co-feed, have significantly net-negative carbon emissions and the most favorable prospective economics when carbon emissions are priced. For these facilities, the crude oil price required for plant economic viability falls rapidly from $\$$100/bbl as carbon emission prices increase above $\$$120/tonne CO 2eq . In general, plants that co-fire lignite with biomass are less profitable (than 100% biomass plants) due to their higher net greenhouse gas emissions.

09 BIOMASS FUELS↗

ROADRUNNER uranium nitride MiniFuel: Experimental design, fabrication and pre-irradiation baseline characterization for accelerated burnup testing

Uranium nitride (UN) is a promising fuel candidate for advanced reactor systems owing to its high uranium density and thermal conductivity; however, its qualification remains constrained by the scarcity of well-controlled irradiation performance data. Here, to address this limitation, the ROADRUNNER (Research On ADvancing the peRformance of UraNium Nitrides in Extreme enviRonments) campaign employs the MiniFuel platform in the High Flux Isotope Reactor (HFIR) to enable accelerated burnup irradiation testing under tightly controlled and largely isothermal conditions. This paper presents the experimental design, fuel fabrication, and pre-irradiation baseline characterization of the ROADRUNNER UN MiniFuel campaign. Thirty-six UN minidisc specimens were fabricated with systematically varied as-fabricated density (86–96% of theoretical density), carbon impurity content (961–5240 ppm), oxygen content (≤ ∼2000 ppm), and grain size (2.5–24 μm). The irradiation matrix spans nominal fuel temperatures of 873 K, 1173 K, and 1473 K and target burnups of 3.75%, 6.0%, and 7.5% fissions per initial metal atom (FIMA). Neutronic and thermal analyses were performed to define specimen-specific burnup accumulation and temperature histories, establishing the boundary conditions for subsequent in-pile behavior. Comprehensive pre-irradiation characterization—including dimensional metrology, density verification, impurity analysis, X-ray diffraction, Raman spectroscopy, scanning electron microscopy, X-ray computed tomography, and confocal profilometry—provides a detailed baseline for post-irradiation examination. Pre-irradiation data were further used to generate predictive estimates of fission gas release and swelling using existing empirical correlations. This quantitative comparison reveals substantial inter-model divergence at intermediate and elevated temperatures that exceeds propagated input uncertainties, highlighting structural gaps in the historical irradiation database. The ROADRUNNER irradiation campaign is currently underway in HFIR, with initial firs cycle completed in late 2025 and remaining targets scheduled through 2027. The experimental design and baseline dataset presented here establish the framework needed to interpret forthcoming post-irradiation measurements and to provide discriminating data for the validation and refinement of physics-based UN fuel performance models.

Lopes, Denise Adorno [Oak Ridge National Laborator↗

Transverse single-spin asymmetries in $γ$ SIDIS as a direct probe of quark-gluon-quark longitudinal momentum structure

Transverse single-spin asymmetries in the semi-inclusive deep-inelastic production of isolated photons ($\gamma$SIDIS), $A_{UT}^{\gamma {\rm SIDIS}}$, provide an unprecedented opportunity to extract the quark-gluon-quark correlators $F_{FT}(x,x')$ and $G_{FT}(x,x')$ point-by-point in their full support $x,x'$. We utilize realistic models for these functions, based on input from the Sivers transverse momentum dependent parton distribution function and imposing constraints from the $d_2$ matrix element calculated in lattice QCD, in order to provide numerical estimates for $A_{UT}^{\gamma {\rm SIDIS}}$ at the Electron-Ion Collider (EIC). We thoroughly explore the EIC phase space in order to isolate in which regions the asymmetry can be sizable, finding it can be as much as $10\%$ or larger for certain kinematics. Given that $F_{FT}(x,x')$ and $G_{FT}(x,x')$ are basically unknown, $A_{UT}^{\gamma {\rm SIDIS}}$ will be an important future measurement to learn about multi-parton correlations in the nucleon.

Electron-Ion Collider↗

Conditioned Simulation of Ground-Motion Time Series at Uninstrumented Sites Using Gaussian Process Regression

Ground-motion time series are essential input data in seismic analysis and performance assessment of the built environment. Because instruments to record free-field ground motions are generally sparse, methods are needed to estimate motions at locations with no available ground-motion recording instrumentation. In this study, given a set of observed motions, ground-motion time series at target sites are constructed using a Gaussian process regression (GPR) approach, which treats the real and imaginary parts of the Fourier spectrum as random Gaussian variables. Model training, verification, and applicability studies are carried out using the physics-based simulated ground motions of the 1906 Mw 7.9 San Francisco earthquake and Mw 7.0 Hayward fault scenario earthquake in northern California. Additionally, the method’s performance is further evaluated using the 2019 Mw 7.1 Ridgecrest earthquake ground motions recorded by the Community Seismic Network stations located in southern California. These evaluations indicate that the trained GPR model is able to adequately estimate the ground-motion time series for frequency ranges that are pertinent for most earthquake engineering applications. The trained GPR model exhibits proper performance in predicting the long-period content of the ground motions as well as directivity pulses.

58 GEOSCIENCES↗

Quality Assurance Plan for Federal Guidance Report 16

Federal Guidance Report No. 16 will provide risk coefficients for use in projecting the probability of cancer morbidity or mortality from exposure to environmental radionuclides. The derivation of nuclide-specific risk coefficients involves multiple steps, with each step relying on the accuracy of the previous step as well as the numerical data utilized within each step. Thus, to ensure the quality of the derived risk coefficients, it is necessary to evaluate and confirm the quality and proper implementation of input data at each step of the derivation. The Oak Ridge National Laboratory (ORNL) dosimetry research team currently has two largely independent software codes for deriving radiation dose estimates, dose coefficients, and cancer risk coefficients for internally deposited radionuclides: DCAL and QCAL. The agreement of derived values based on alternate codes is an important step in quality assurance (QA) of the risk coefficients produced in this project. DCAL is an integration of the dosimetric software and numerical databases developed at ORNL over the past 25 years for the EPA’s Office of Radiation and Indoor Air (ORIA). QCAL was developed in 2015 to serve as a substitute for DCAL while DCAL was undergoing revision and to provide QA for the updated version of DCAL. The complex methodology employed in the derivation of cancer risk coefficients presents a challenge for QA, and as such, this report describes the QA steps used by the ORNL dosimetry team in this process as implemented in DCAL and QCAL for Federal Guidance Report No. 16.

61 RADIATION PROTECTION AND DOSIMETRY↗

Evaluated 238 U(n,f) Average Prompt Fission Neutron Multiplicities Including the CGMF Model

This report documents an evaluation of the average prompt fission neutron multiplicity, $\overline{v}_p$, of 238 U from 800 keV to MeV. This evaluation had to be re-done from “scratch” as the input to previous $\overline{v}_p$ evaluations, specifically ENDF/B-VIII.0, was not found. That means that all available experimental data were re-analyzed and uncertainties were re-estimated. The new evaluated 238 U $\overline{v}_p$ based on only experimental data differs distinctly from ENDF/B-VIII.0 $\overline{v}_p$ from 2 to 4.5 MeV, and from 6 to 7 MeV, and is otherwise similar. The difference from 2 to 4.5 MeV stems from the fact that ENDF/B-VIII.0 was tweaked in this energy range to data of Frehaut, while two other, equally trustworthy, data sets would indicate an evaluated 238 U $\overline{v}_p$ that is up to 2% higher. Also, second chance fission in ENDF/B-VIII.0 was smoothed over from 6–7 MeV. Another major difference to ENDF/B-VIII.0 is that one of the evaluations presented here includes model information from the Hauser-Feshbach fission fragment decay code CGMF, while ENDF/B-VIII.0 is based purely on experimental data. CGMF links several fission quantities with each other; $\overline{v}_p$ is predicted by assumptions made on, e.g., pre-neutron emission yields as a function of mass, the total kinetic energy, or spin and parity of fission fragments. This allows to validate the new 238 U $\overline{v}_p$ by using CGMF parameters obtained from fitting to experimental 238 U $\overline{v}_p$ to predict yields as a function of mass, the average total kinetic energy, or the mean energy of the prompt fission neutron spectrum. These model-predicted values can then be compared to experimental and evaluated data. The model-predicted fission-observable values using evaluated parameters obtained here are reasonably close to experimental data for some observables, but are farther away from experimental data related to TKE observables. In addition to that, the evaluated 238 U(n,f) $\overline{v}_p$ shows similar deviations from ENDF/B-VIII.0 as for the evaluation with only experimental data. This difference is expected to lead to changes in simulated effective neutron multiplication factor, $k_{eff}$ of ICSBEP critical assemblies that are sensitive to 238 U in the fast range (BigTen, Flattop, Flattop-Pu). These changes in $k_{eff}$ need to be counter-balanced. Chi-Nu PFNS experimental data are expected to be released in the next few months that might lead to the needed changes in the PFNS. Until then, we hold off in benchmarking the new 238 U(n,f) $\overline{v}_p$ as well as submitting it to ENDF/B-VIII.1. Also, new high-precision 238 U $\overline{v}_p$ are expected to be measured by the CEA in the next two years that will shed further light on question on 238 U $\overline{v}_p$ from 2–4.5 and 6–7 MeV.

238U↗

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗

Space shuttle orbital maneuvering system failure detection and identification software requirements (uncontrolled)

Candidate designs and their software implementation are presented for the Orbital Maneuvering System (OMS) Failure Detection and Identification (FDI) algorithms in the Redundance Management (RM) module of the Space Shuttle Guidance, Navigation, and Control (GN&C) software. The OMS engine FDI algorithm monitors OMS engine thrust performance, and the OMS actuator FDI algorithm monitors OMS gimbal actuator performance. The software functional requirements of the algorithms are described along with the objective of each algorithm. A list of the assumptions which have governed its design, input/output requirements, a functional description of the algorithm (including a functional block diagram), and input interface requirements are given. The HAL (the language of the space shuttle flight computer) software formulation of the algorithms is considered including structured flowcharts of the procedures, estimates of flight computer core storage and CPU time, and processing requirements. A glossary of the symbols used to define the software requirements and formulations is included.

Damario, L. A.↗

Estimation of characteristics and stochastic control of an aircraft flying in atmospheric turbulence

An adaptive control technique to improve the flying qualities of an aircraft in turbulence was investigated. The approach taken was to obtain maximum likelihood estimates of the unknown coefficients of the aircraft system and then, using these estimates along with the separation principle, to define the stochastic optimal control. The maximum likelihood estimation technique that accounted for the effects of turbulence provided good estimates of the unknown coefficients and of the turbulence. The assessment of the stochastic optimal control based on the maximum likelihood estimates showed that the desired effects were attained for the regulator problem of minimizing pitch angle and the tracking problem of requiring normal acceleration to follow the pilot input.

Iliff, K. W.↗

A parameter estimation subroutine package

Linear least squares estimation and regression analyses continue to play a major role in orbit determination and related areas. A library of FORTRAN subroutines were developed to facilitate analyses of a variety of estimation problems. An easy to use, multi-purpose set of algorithms that are reasonably efficient and which use a minimal amount of computer storage are presented. Subroutine inputs, outputs, usage and listings are given, along with examples of how these routines can be used. The routines are compact and efficient and are far superior to the normal equation and Kalman filter data processing algorithms that are often used for least squares analyses.

Bierman, G. J.↗

A comprehensive method for preliminary design optimization of axial gas turbine stages. II - Code verification

The present effort represents an extension of previous work wherein a calculation model for performing rapid pitchline optimization of axial gas turbine geometry, including blade profiles, is developed. The model requires no specification of geometric constraints. Output includes aerodynamic performance (adiabatic efficiency), hub-tip flow-path geometry, blade chords, and estimates of blade shape. Presented herein is a verification of the aerodynamic performance portion of the model, whereby detailed turbine test-rig data, including rig geometry, is input to the model to determine whether tested performance can be predicted. An array of seven (7) NASA single-stage axial gas turbine configurations is investigated, ranging in size from 0.6 kg/s to 63.8 kg/s mass flow and in specific work output from 153 J/g to 558 J/g at design (hot) conditions; stage loading factor ranges from 1.15 to 4.66.

Jenkins, R. M.↗

Derated ion thruster design issues

Preliminary activities to develop and refine a lightweight 30 cm engineering model ion thruster are discussed. The approach is to develop a 'derated' ion thruster capable of performing both auxiliary and primary propulsion roles over an input power range of at least 0.5 to 5.0 kilo-W. Design modifications to a baseline thruster to reduce mass and volume are discussed. Performance data over an order of magnitude input power range are presented, with emphasis on the performance impact of engine throttling. Thruster design modifications to optimize performance over specific power envelopes are discussed. Additionally, lifetime estimates based on wear test measurements are made for the operation envelope of the engine.

Patterson, Michael J.↗

Aberrations of a horizontal-vertical depolarizer

Ray-trace equations for uniaxial birefringent materials are used here to derive third-order estimates for aberrations that are produced in imaging through uniaxial plates and horizontal-vertical (HV) depolarizers. An HV depolarizer is a spatial pseudodepolarizer; it converts a uniform input polarization state into a continuum of spatially varying polarization states in an output beam. An HV depolarizer consists of two birefringent wedges whose crystal axes are crossed at 90 deg. The interface between the wedges is included, which leads to a spatially varying retardance that provides the spatial pseudodepolarization. In HV depolarizers, spherical aberration, astigmatism, and image doubling are the principal aberrations for on-axis objects. Only spherical aberration occurs in isotropic plates, while the presence of birefringent wedges introduces astigmatism and image doubling. It is shown that image separation is proportional to the magnitude of the retardance variation.

Mcclain, Stephen C.↗

Software Fault Tolerance: A Tutorial

Because of our present inability to produce error-free software, software fault tolerance is and will continue to be an important consideration in software systems. The root cause of software design errors is the complexity of the systems. Compounding the problems in building correct software is the difficulty in assessing the correctness of software for highly complex systems. After a brief overview of the software development processes, we note how hard-to-detect design faults are likely to be introduced during development and how software faults tend to be state-dependent and activated by particular input sequences. Although component reliability is an important quality measure for system level analysis, software reliability is hard to characterize and the use of post-verification reliability estimates remains a controversial issue. For some applications software safety is more important than reliability, and fault tolerance techniques used in those applications are aimed at preventing catastrophes. Single version software fault tolerance techniques discussed include system structuring and closure, atomic actions, inline fault detection, exception handling, and others. Multiversion techniques are based on the assumption that software built differently should fail differently and thus, if one of the redundant versions fails, it is expected that at least one of the other versions will provide an acceptable output. Recovery blocks, N-version programming, and other multiversion techniques are reviewed.

Torres-Pomales, Wilfredo↗

Self-Tuning of Design Variables for Generalized Predictive Control

Three techniques are introduced to determine the order and control weighting for the design of a generalized predictive controller. These techniques are based on the application of fuzzy logic, genetic algorithms, and simulated annealing to conduct an optimal search on specific performance indexes or objective functions. Fuzzy logic is found to be feasible for real-time and on-line implementation due to its smooth and quick convergence. On the other hand, genetic algorithms and simulated annealing are applicable for initial estimation of the model order and control weighting, and final fine-tuning within a small region of the solution space, Several numerical simulations for a multiple-input and multiple-output system are given to illustrate the techniques developed in this paper.

Lin, Chaung↗