SEARCH · Engineering Papers
Results for “Model calibration”
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.
Quasi-Newton Variational Bayes v.3.0
SAND2024-09052O Quasi-Newton Variational Bayes (QNVB) is a training algorithm that performs high-dimensional variational inference on machine learning models. It provides mechanisms to calibrate and control model uncertainty during training. The PyTorch implementation of projective integral updates for Gaussian mean-fields supports the manuscript, "Projective Integral Updates for High-Dimensional Variational Inference." Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
2014-2015 Puget Sound Regional Travel Study
The 2014-2015 Puget Sound Regional Travel Study collected information about household and individual travel patterns for residents throughout a four-county region in Washington State. Study results were used to update the region's travel and land-use models and to calibrate local traffic and travel models. The study also helped the Puget Sound Regional Council (PSRC) and its regional partners develop plans that accommodate the diverse travel needs and preferences of residents. The Resource Systems Group administered the study on behalf of PSRC. Global positioning system (GPS)-equipped smartphones were used to provide data pertaining to the daily travel of 547 individual participants. Because the region's university students may have been underrepresented in the initial 2014 household travel study, the PSRC added a college-population travel survey in fall 2014. In spring 2015, a second household data collection effort was conducted to increase the frequency of data collection and to collect GPS data as well as a sample of longitudinal data from households that completed the 2014 survey.
Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling
The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3, documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. The overall objective of the saturated zone modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within and downgradient of the Central Plateau of the Hanford Site. Specifically, the purpose of this environmental calculation is to describe the application of the hydraulic property fields and recharge parameters documented in ECF-HANFORD-20-0027, Null Space Monte Carlo Evaluation of the Plateau to River Model, to the CA flow and fate and transport simulation results to quantify the uncertainty in the simulated results due to input parameter selection. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, models lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. The result of NSCM analysis is a set of F&T simulations that provide an estimate of the range of possible outcomes that are used to quantify the uncertainty in simulated concentrations produced using the base case simulations. The simulated concentrations from all simulations will support calculation of the uncertainty of the total dose calculated in a separate calculation.
Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling
The Plateau-to-River (P2R) Model is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by Central Plateau Cleanup Company (CPCCo) at the Hanford Site in Washington State. Figure 1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the composite analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3 documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application.
Predictive Contaminant Transport Simulation with P2R Model for the Composite Analysis Limited Source Sensitivity Case.
The Plateau to River (P2R) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by the Central Plateau Cleanup Company at the Hanford Site in Washington State. Figure 1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3 documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. Simulations conducted to support the dose calculations required by the CA are documented in ECF-HANFORD-19-0119, Predictive Flow Simulation with the P2R Model for the Composite Analysis Base Case and ECF-HANFORD-19-0120, Contaminant Transport Simulation with the P2R Model for the Composite Analysis Base Case.
CMIP6-based Multi-model Streamflow Projections over the Conterminous US
This dataset presents an ensemble of streamflow projections based on the hydroclimate projections dataset supporting the SECURE Water Act Section 9505 Assessment for the US Department of Energy (DOE) Water Power Technologies Office (WPTO). The six-member General Climate Model (GCM) ensemble from the Coupled Models Intercomparison Project phase 6 (CMIP6) downscaled using statistical (i.e., DBCCA) and dynamical (i.e., RegCM) and bias-corrected using two meteorological reference observations (Daymet & Livneh) are driven through two calibrated hydrologic models (VIC & PRMS) to simulate projected future hydrologic responses. This leads to the production of an ensemble of hydroclimate projections, including total runoff (surface runoff and baseflow) for 1980–2019 baseline and 2020–2059 near-term future periods under multiple emission scenarios at 1/24° (~4 km) spatial resolution across the CONUS. The total runoff projections are routed through the Routing Application for Parallel computatIon of Discharge (RAPID) routing model to produce an ensemble of streamflow projections for both periods across 2.7 million NHDPlusV2 stream reaches in the CONUS.
Utilizing remote sensing of thematic mapper data to improve our understanding of estuarine processes and their influence on the productivity of estuarine-dependent fisheries
The land-water interface of coastal marshes may influence the production of estuarine-dependent fisheries more than the area of these marshes. To test this hypothesis, a spatial model was created to explore the dynamic relationship between marshland-water interface and level of disintegration in the decaying coastal marshes of Louisiana's Barataria, Terrebonne, and Timbalier basins. Calibrating the model with Landsat Thematic Mapper satellite imagery, a parabolic relationship was found between land-water interface and marsh disintegration. Aggregated simulation data suggest that interface in the study area will soon reach its maximum and then decline. A statistically significant positive linear relationship was found between brown shrimp catch and total interface length over the past 28 years. This relationship suggests that shrimp yields will decline when interface declines, possibly beginning about 1995.
Dynamical Model for the Zodiacal Cloud and Sporadic Meteors
The solar system is dusty, and would become dustier over time as asteroids collide and comets disintegrate, except that small debris particles in interplanetary space do not last long. They can be ejected from the solar system by Jupiter, thermally destroyed near the Sun, or physically disrupted by collisions. Also, some are swept by the Earth (and other planets), producing meteors. Here we develop a dynamical model for the solar system meteoroids and use it to explain meteor radar observations. We find that the Jupiter Family Comets (JFCs) are the main source of the prominent concentrations of meteors arriving to the Earth from the helion and antihelion directions. To match the radiant and orbit distributions, as measured by the Canadian Meteor Orbit Radar (CMOR) and Advanced Meteor Orbit Radar (AMOR), our model implies that comets, and JFCs in particular, must frequently disintegrate when reaching orbits with low perihelion distance. Also, the collisional lifetimes of millimeter particles may be longer (approx. > 10(exp 5) yr at 1 AU) than postulated in the standard collisional models (approx 10(exp 4) yr at 1 AU), perhaps because these chondrule-sized meteoroids are stronger than thought before. Using observations of the Infrared Astronomical Satellite (IRAS) to calibrate the model, we find that the total cross section and mass of small meteoroids in the inner solar system are (1.7-3.5) 10(exp 11) sq km and approx. 4 10(exp 19) g, respectively, in a good agreement with previous studies. The mass input required to keep the Zodiacal Cloud (ZC) in a steady state is estimated to be approx. 10(exp 4)-10(exp 5) kg/s. The input is up to approx 10 times larger than found previously, mainly because particles released closer to the Sun have shorter collisional lifetimes, and need to be supplied at a faster rate. The total mass accreted by the Earth in particles between diameters D = 5 micron and 1 cm is found to be approx 15,000 tons/yr (factor of 2 uncertainty), which is a large share of the accretion flux measured by the Long Term Duration Facility (LDEF). Majority of JFC particles plunge into the upper atmosphere at <15 km/s speeds, should survive the atmospheric entry, and can produce micrometeorite falls. This could explain the compositional similarity of samples collected in the Antarctic ice and stratosphere, and those brought from comet Wild 2 by the Stardust spacecraft. Meteor radars such as CMOR and AMOR see only a fraction of the accretion flux (approx 1- 10% and approx 10-50%, respectively), because small particles impacting at low speeds produce ionization levels that are below these radars detection capabilities.
Assessment of the Uniqueness of Wind Tunnel Strain-Gage Balance Load Predictions
A new test was developed to assess the uniqueness of wind tunnel strain-gage balance load predictions that are obtained from regression models of calibration data. The test helps balance users to gain confidence in load predictions of non-traditional balance designs. It also makes it possible to better evaluate load predictions of traditional balances that are not used as originally intended. The test works for both the Iterative and Non-Iterative Methods that are used in the aerospace testing community for the prediction of balance loads. It is based on the hypothesis that the total number of independently applied balance load components must always match the total number of independently measured bridge outputs or bridge output combinations. This hypothesis is supported by a control volume analysis of the inputs and outputs of a strain-gage balance. It is concluded from the control volume analysis that the loads and bridge outputs of a balance calibration data set must separately be tested for linear independence because it cannot always be guaranteed that a linearly independent load component set will result in linearly independent bridge output measurements. Simple linear math models for the loads and bridge outputs in combination with the variance inflation factor are used to test for linear independence. A highly unique and reversible mapping between the applied load component set and the measured bridge output set is guaranteed to exist if the maximum variance inflation factor of both sets is less than the literature recommended threshold of five. Data from the calibration of a six{component force balance is used to illustrate the application of the new test to real-world data.
Development of a thermal creep model for aluminum alloy 6061 cladding in U-10Mo monolithic fuel plates
Plate-type fuel elements consisting of a high-density, low-enriched uranium (LEU) U–10Mo-based fuel foil encapsulated in an aluminum alloy (AA) cladding are fabricated using the hot isostatic pressing (HIP) technique. During the HIP process, the fuel plate system is heated to 560 °C, then cooled to room temperature. This heat cycle significantly affects the mechanical properties of the aluminum cladding, and experimental investigations have shown that, post-HIP bonding, the mechanical properties of the aluminum cladding transition from those of AA 6061-T6 to something closer to the O temper. More specifically, the ultimate strength of the cladding decreases while its ductility increases, making it challenging to capture the changes in mechanical behavior and material properties. Understanding the residual stresses generated during the HIP process is critical for assessing the fuel plate’s integrity under various temperature, pressure, and irradiation. To simulate the HIP bonding process, the elastic, plastic, and thermal properties of the cladding are assumed to be similar to those of AA 6061-O temper. However, the primary challenge lies in the lack of available data for the creep model of the AA 6061 cladding during this transient process of HIP. The present study focuses on developing a computational model that predicts the creep behavior of the aluminum cladding in the fuel plates during the HIP process, as cladding creep significantly influences the residual stresses generated in U-10Mo fuel plates during HIP fabrication. Furthermore, as HIP bonding occurs at high temperatures that are nearing the melting point of aluminum, the present work considered a temperature-dependent Arrhenius-type creep model. In particular, a hyperbolic sine creep model is employed to estimate the creep properties of the as-fabricated aluminum cladding. In conclusion, the residual stresses predicted in the U-10Mo fuel when using the newly calibrated creep model closely align with the experimental measurements, validating the model’s accuracy.
Electrical characterization of a Space Station Freedom alpha utility transfer assembly
Electrical power, command signals and data are transferred across the Space Station Freedom solar alpha rotary joint by roll rings, which are incorporated within the Utility Transfer Assembly (UTA) designed and manufactured by Honeywell Space Systems Operations. A developmental Model of the UTA was tested at the NASA Lewis Research Center using the Power Management and Distribution DC test bed. The objectives of these tests were to obtain data for calibrating system models and to support final design of qualification and flight units. This testing marked the first time the UTA was operated at high power levels and exposed to electrical conditions similar to that which it will encounter on the actual Space Station. Satisfactory UTA system performance was demonstrated within the scope of this testing.
Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches
Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.
Modeling Laminar-to-Turbulent Transition in the Panel Test Facility Arcjet
A computational fluid dynamics (CFD) laminar-to-turbulence transition model was developed for the NASA Ames Research Center’s 20 MW Panel Test Facility (PTF). Surface pressure, heat flux to a water-cooled plate, and surface temperature on a tile plate coated with reaction-cured glass were measured across several conditions in the facility and compared with laminar and fully-turbulent CFD simulations. The potential for bypass transition in the PTF nozzle was assessed via application of the Langtry-Menter four-equation transitional shear-stress transport (SST) model. Results from the transition model were inconsistent with measurements. Flow interaction with a boundary conditioning plate feature inside the nozzle was also investigated as a potential source of laminar-to-turbulent transition using two turbulence models with specified transition locations. The Baldwin-Lomax turbulence model was configured to simulate a transition at the upstream edge of the boundary conditioning plate and produced results consistent with the surface pressure measurements but not the cold-wall heat flux. Finally, the SST turbulence model was calibrated to transition at the upstream edge of the boundary conditioning plate and produced results consistent with both the surface pressure and cold-wall heat flux measurements. The SST-based model demonstrated reasonable agreement with surface temperature measurements on the reaction-cured glass tile, albeit with some discrepancies.
Flux Comparison of Master-8 and Ordem 3.1 Modelled Space Debris Population
With ESA’s Meteoroid And Space debris Terrestrial Environment Reference (MASTER-8) model and NASA’s Orbital Debris Engineering Model(ORDEM) 3.1, the two premier orbital debris engineering models have been officially released. The two models come with significant enhancements and now represent state-of-the-art orbital debris modelling for their respective agencies. Both models provide the community with estimates of the space debris environment from low Earth orbit (LEO) up to at least geostationary altitude. The MASTER population is an event-based simulation of all known events that generate debris and objects that are part of the U.S. Space Surveillance Network (SSN) catalog, which provides coverage of objects with diameters down to approximately 10cm in LEO and 1 m in geosynchronous Earth orbit (GEO). Different models are used to simulate the artificial objects and their orbit evolution over time. These models are called “sources” since they assign an origin to each individual object and consist of fragments, solid rocket motor (SRM) remainders, sodium-potassium(NaK) droplets, paintflakes, ejecta, and multi-layer insulation (MLI) fragments. The objects from each source are characterized by having individual release mechanisms, as well as orbital distributions, material composition, size, and mass distributions. Dedicated radar and telescope observation data is used to calibrate the model for objects larger than 1 cm in LEO and larger than 10 cm in GEO. For calibrating the small-sized objects, below 1 cm, impact data from returned surfaces are analyzed. Because the >1cm object population is dominated by fragments, the fragmentation event database was updated to include new events, as well as re-evaluate past events. Special attention was drawn to re-evaluating theFengyun-1C anti-satellite test from 2007 and Cosmos-Iridium collision event from 2009 since these events shape the fragment population because of their severity. After 2009, the two largest fragmentations in terms of number of tracked debris are the Briz-M explosion in 2012 and the NOAA-16explosion in 2015. In total, there are 261 confirmed fragmentations in the database up to November 2016.The baseline population for ORDEM 3.1 is based on the U.S. SSN catalog, and observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller size regions. These regions are not well-covered by the SSN catalog yet may pose the greatest threat to operational spacecraft. The NASA Standard Satellite Breakup Model is used to generate fragments greater than 1mm from collisions and explosions, and these fragment populations are scaled using ground-based radar data. Specific major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos2251 debris clouds, and unique populations, such as NaK droplets, were re-examined, modelled, and added to the ORDEM environment separately. Optical measurement data is used to model the GEO population down to 10 cm. The debris environment is propagated using NASA’s LEO-to-GEO Environment Debris model, and future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically. The environment from a few millimetres down to 10 𝜇m is modelled using a special degradation model where small particles are generated from intact spacecraft and rocket bodies, then the populations are scaled to fit in situ cratering data from Space Shuttle returned surfaces. Fragments smaller than 10 cm are differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to upper stages and spacecraft. This paper will discuss the MASTER and ORDEM approaches for modelling populations and compare fluxes for specific orbits, including sun-synchronous, ISS-altitude, geosynchronous transfer, and GEO. In the end, a conclusion is drawn towards the importance of having multiple fundamentally different, yet validated, models to estimate the space debris population.
MatCal Users Guide: Release 1.3.0
Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.
How Accurately Do Maize Crop Models Simulate the Interactions of Atmospheric CO2 Concentration Levels With Limited Water Supply on Water Use and Yield?
This study assesses the ability of 21 crop models to capture the impact of elevated CO2 concentration [CO2] on maize yield and water use as measured in a 2-year Free Air Carbon dioxide Enrichment experiment conducted at the Thunen Institute in Braunschweig, Germany (Manderscheid et al. 2014). Data for ambient [CO2] and irrigated treatments were provided to the 21 models for calibrating plant traits, including weather, soil and management data as well as yield, grain number, above ground biomass, leaf area index, nitrogen concentration in biomass and grain, water use and soil water content. Models differed in their representation of carbon assimilation and evapotranspiration processes. The models reproduced the absence of yield response to elevated [CO2] under well-watered conditions, as well as the impact of water deficit at ambient [CO2], with 50 percent of models within a range of plus/minus 1 Mg ha(exp. -1) around the mean. The bias of the median of the 21 models was less than 1 Mg ha(exp. -1). However under water deficit in one of the two years, the models captured only 30 percent of the exceptionally high [CO2] enhancement on yield observed. Furthermore the ensemble of models was unable to simulate the very low soil water content at anthesis and the increase of soil water and grain number brought about by the elevated [CO2] under dry conditions. Overall, we found models with explicit stomatal control on transpiration tended to perform better. Our results highlight the need for model improvement with respect to simulating transpirational water use and its impact on water status during the kernel-set phase.
Dynamic calibration of differential equations using machine learning, with application to turbulence models
We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.