Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Mesoscale Particle-Based Model of Electrophoresis

Here, we develop and evaluate a semi-empirical particle-based model of electrophoresis using extensive mesoscale simulations. We parameterize the model using only measurable quantities from a broad set of colloidal suspensions with properties that span the experimentally relevant regime. With sufficient sampling, simulated diffusivities and electrophoretic velocities match predictions of the ubiquitous Stokes-Einstein and Henry equations, respectively. This agreement holds for non-polar and aqueous solvents or ionic liquid colloidal suspensions under a wide range of applied electric fields.

36 MATERIALS SCIENCE↗
DOE OSTI2015

“Using Statistical Comparisons between SPartICus Cirrus Microphysical Measurements, Detailed Cloud Models, and GCM Cloud Parameterizations to Understand Physical Processes Controlling Cirrus Properties and to Improve the Cloud Parameterizations”

The dual objectives of this project were improving our basic understanding of processes that control cirrus microphysical properties and improvement of the representation of these processes in the parameterizations. A major effort in the proposed research was to integrate, calibrate, and better understand the uncertainties in all of these measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Parameterization of a Cookoff Model for LX-07

A thermal decomposition model for PBX 9501 (95% HMX, 2.5% Estane® binder, 2.5% BDNPA/F nitro-plasticizer) was implemented by Dickson, et. al. The objective in this study is to estimate parameters associated with this kinetics model so it can be applied to carry out thermal ignition predictions for LX-07 (90% HMX, 10% Viton binder). LX-07 thermal ignition experiments have been carried out using the “Sandia Instrumented Thermal Ignition Apparatus”, SITI. The SITI design consists of solid cylinders (1” diameter × 1” height) of high explosive (HE) confined by a cylindrical aluminum case. An electric heater is wrapped around the outer surface of the case. This heater produces a temperature heating ramp on the outer surface of the case. Internal thermocouples measure the HE temperature rise from the center to locations close to the HE-aluminum interface. The energetic material is heated until thermal ignition occurs. A two–dimensional axisymmetric heat conduction finite element model is used to simulate these experiments. The HE thermal decomposition kinetics is coupled to a heat conduction model trough the definition of an energy source term. The parameters used to define the HE thermal decomposition model are optimized to obtain a good agreement with the experimental time to thermal ignition and temperatures. Also, heat capacity and thermal conductivity of the LX-07 mixture were estimated using temperatures measured at the center of the HE before the solid to solid HMX phase transition occurred.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2018

A parameterized ray-tracing model for evaluation of hemispherical lensing of a uniform, telecentric object aperture -- performance improvement estimates for segmented scintillators with miniature lens arrays

A parametric model describing the collimation of radiation from a uniform circular aperture by a hemispherical lens of the same radius has been developed for rapid estimation of the potential benefit of miniature dielectric lens arrays in lens-coupled gamma ray cameras (GRCs) with segmented scintillators. A simplified ray-tracing approach is adopted that makes use of small angle approximations relevant to GRCs used in nuclear science at Los Alamos National Laboratory and elsewhere. Application of this analysis technique in a first-pass design optimization process results in the optimistic estimate that miniature lenses fabricated from typical optical quality glass might boost optical coupling by as much as a factor of 4 (6 dB gain). However, this is not likely to be economically attractive compared to other techniques, such as reduction of the GRC lens f-number.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Biological Mechanisms May Contribute to Soil Carbon Saturation Patterns: Modeling Archive

This Modeling Archive is in support of a TES-SFA publication “Biological Mechanisms May Contribute to Soil Carbon Saturation Patterns” (Craig et al., 2021). We ran and evaluated a multi-assumption soil organic carbon (SOC) model to investigate whether alternative assumptions regarding constraints on soil microbial biomass could lead to soil carbon saturation patterns. We developed this model in the Multi-Assumption Architecture and Testbed (MAAT, https://github.com/walkeranthonyp/MAAT, tag: v1.2.1_Craig2021; Walker et al. 2018). Using MAAT, we embedded three alternative hypotheses in a microbially explicit three-pool SOC model: 1) the efficiency of mineral-associated SOC formation decreases as mineral-associated SOC approaches a maximum value (“Mineral saturation”), 2) the microbial biomass turnover rate increases with increasing microbial biomass (“Density-dependent turnover”), and 3) community carbon use efficiency decreases as microbial biomass increases toward an upper limit (“Density-dependent growth”). We ran a factorial combination of these hypotheses resulting in eight models for three different classes of model (linear decay, Michaelis-Menten decay, or reverse Michaelis-Menten decay), resulting in 24 models, 12 of which are presented or discussed in the related publication. Models were parameterized using values from previous studies with similar models (Wang et al. 2013, Wieder et al. 2014, Li et al. 2014, Georgiou et al. 2017, Hassink and Whitmore 1997) and ran to an approximate steady state (200 years) at six (6) different C input rates corresponding to 0.5, 1, 2, 4, 7, and 10 times the default input value. Further model details are available in the related publication. This archive contains output from three MAAT simulations, and scripts to run these simulations and process and plot the data. Simulations are labeled “lin”, “MM_highKm”, and “RMM_highKm” reflecting factorial runs for linear, Michealis-Menten, and reverse Michaelis-Menten models, respectively. This archive contains: • 3 R scripts prepended with “init_MAAT_” to initialize runs (1 for each simulation), • 1 csv file containing years over which to run simulations (“met_year.csv”), • 1 bash script (.bs) to run MAAT, • 6 XML files that are output from MAAT describing a run (2 for each simulation), • 3 model output csv files prepended by “out_” (1 for each simulation), and • 1 analysis R script for reproducing figures 3 and 4 in Craig et al. 2021. See included user guide (Craig_2021_modeling_archive_20210315.pdf) for file organization details.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Technical Note: An efficient method for accelerating the spin-up process for process-based biogeochemistry models

To better understand the role of terrestrial ecosystems in the global carbon cycle and their feedbacks to the global climate system,process-based biogeochemistry models need to be improved with respect tomodel parameterization and model structure. To achieve these improvements,the spin-up time for those differential equation-based models needs to be shortened. Here, an algorithm for a fast spin-up was developed by finding the exact solution of a linearized system representing the cyclo-stationary state of a model and implemented in a biogeochemistry model, the Terrestrial Ecosystem Model (TEM). With the new spin-up algorithm, we showed that the model reacheda steady state in less than 10 years of computing time, while the original method requires more than 200 years on average of model run. For the test sites with five different plant functional types, the new method saves over 90 % of the original spin-up time in site-level simulations. In North American simulations, average spin-up time savings for all grid cells is 85 % for either the daily or monthly version of TEM. Thus, the developed spin-up method shall be used for future quantification of carbon dynamics at fine spatial and temporal scales.

58 GEOSCIENCES↗
DOE OSTI2020

Simulation and Emulation of X-Ray Diffraction from Dynamic Compression Experiments

Many important aspects of the dynamic thermo-mechanical response of materials occur at the mesoscale, i.e. a physical scale of interactions smaller than what can be adequately described by homogenous behaviors, yet larger than the scale of the atomic lattice. Concurrent advancements in computational power, continuum theory, and experimental diagnostics are enabling unprecedented understanding of such interactions. However, we cannot develop a sufficient level of confidence in such mesoscale capability until the constitutive description of the underlying constituents is reliably representative of their actual physical behavior. Therefore, there is a strong need to combine experimental, modeling, and data-science techniques to validate models of the thermomechanical response of individual single crystals. One experimental diagnostic with high potential impact to shock physics and materials science is in-situ x-ray diffraction. This paper is primarily focused on simulation of x-ray diffraction in shock physics, but with an aim toward quantifying parametric uncertainty of simulation models. Here, we develop and demonstrate a data-science and model-driven approach to constrain the parameterization of continuum models of crystal lattice deformation associated with the shock response of crystalline materials. The framework is built around the connection between continuum hydrodynamic simulations of lattice deformation and a new Bragg diffraction simulation code, BarberShop. The dynamic deformation of a crystal lattice is modeled using the DiscoFlux model within an arbitrary Lagrangian-Eulerian hydrodynamic code, FLAG. These detailed continuum simulations of lattice deformation can be computationally slow, thus a statistical model is used to emulate the evolution of lattice deformation fields in time and across the considered model parameter space. Emulated lattice deformation fields can then be generated rapidly for any combination of physics model parameters. In turn, these fields can be fed into BarberShop to realize a rapid prediction of Bragg diffraction patterns associated with particular values of physics model parameters. The framework enables parameterization of the single crystal model to obtain Bragg diffraction patterns that most closely resemble a corresponding measurement. Furthermore, the framework naturally provides sensitivities of the lattice deformation to the physics parameters. We highlight the utility of this framework through the application to a synthetic closed-loop inverse problem leading to the parameterization of a single crystal material model. As a model problem, we consider the dynamic response of the energetic molecular crystal, cyclotrimethylenetrinitramine (or RDX), under dynamic compression induced by simulated flyer plate impact experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Carbon Organisms Rhizosphere and Protection in Soil Environment model script and input data for soil moisture-respiration responses in tropical forests

Objectives: Climatic drying is predicted for many tropical forests, yet models remain poorly parameterized for tropical forests, hampering predictions of forest-climate feedbacks. We applied an integrated model–experiment approach, parameterizing an ecosystem model Carbon Organisms Rhizosphere and Protection in the Soil Environment (CORPSE) with tropical forest observational data, and comparing model predictions with a field drying manipulation. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We used the field data to parameterize and run tests in the model.Results: Measured CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. We used this data to parameterize the model, which then predicted increased soil CO2 fluxes in wetter and fertile forests with drying, and decreased fluxes in drier, infertile forests. In contrast to model predictions, a chronic throughfall exclusion experiment in the forests initially suppressed soil CO2 fluxes across forests, with sustained suppression after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season), as predicted by the model. The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Code files:CORPSE_array.py: Defines the equations of the CORPSE modelCORPSE_solvers: Functions for running the CORPSE model using either iterative or ordinary differential equation (ODE) solversrun_Panama_sims.py: Read in datasets and run the model simulations for this studyInput data:PanamaGradientEcosystemChem_BT_CPools_20152016CO2_DC_20190615.xlsx: Plot characteristics used in running model simulationsLiCor compiled surface flux only to 2020_03 DC_20200825.xlsx: Surface gas exchange fluxes used in model-data comparisonsPARCHED litterfall data for Ben Sulman LD 20200902.xlsx: Litterfall data used to drive model simulationsInitialization data:state_500y_20190823.csv: Initial state of model pools based on previous spinup runsOutput data:Outputs/prev_moisture_response.csv: Simulations of multiple sites using original model moisture response function.Outputs/updated_moisture_response.csv: Simulations of multiple sites using updated model moisture response function.Outputs/dry15_prev_moisture_response.csv: Simulations with soil moisture reduced by 15%, using original moisture response function.Outputs/dry15_updated_moisture_response.csv: Simulations with soil moisture reduced by 15%, using updated moisture response function.Outputs/dry30_prev_moisture_response.csv: Simulations with soil moisture reduced by 30%, using original moisture response function.Outputs/dry30_updated_moisture_response.csv: Simulations with soil moisture reduced by 30%, using updated moisture response function.Outputs/latestart_prev_moisture_response.csv: Simulations with extended dry season, using original moisture response function.Outputs/latestart_updated_moisture_response.csv: Simulations with extended dry season, using updated moisture response function.Outputs/[site name]_oneyear.csv: One-year simulation for each site in expanded site list using original moisture response function.Outputs/[site name]_oneyear_dried.csv: One-year simulation for each site in expanded site list using original moisture response function, with soil moisture reduced by 25%.Outputs/[site name]_oneyear_updated_moisture_response.csv: One-year simulation for each site in expanded site list using updated moisture response function.Outputs/[site name]_oneyear_updated_moisture_response_dried.csv: One-year simulation for each site in expanded site list using updated moisture response function, with soil moisture reduced by 25%.Field plot location data:There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗
DOE OSTI2016

Global parameterization and validation of a two-leaf light use efficiency model for predicting gross primary production across FLUXNET sites: TL-LUE Parameterization and Validation

We present the first extended validation of satellitemicrowave (MW) liquidwater path (LWP) for low nonprecipitating clouds, from four operational sensors, against ship-borne observations from a three-channel MW radiometer collected along ship transects over the northeast Pacific during May–August 2013. Satellite MW retrievals have an overall correlation of 0.84 with ship observations and a bias of 9.3 g/m2. The bias for broken cloud scenes increases linearly with water vapor path and remains below 17.7 g/m2. In contrast, satelliteMWLWP is unbiased in overcast scenes with correlations up to 0.91, demonstrating that the retrievals are accurate and reliable under these conditions. Satellite MW retrievals produce a diurnal cycle amplitude consistent with ship-based observations (33 g/m2). Observations taken aboard extended ship cruises to evaluate not only satellite MW LWP but also LWP derived from visible/infrared sensors offer a new way to validate this important property over vast oceanic regions.

Zhou, Yanlian↗
DOE OSTI2022

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗
DOE OSTI2023

Multifidelity uncertainty quantification with models based on dissimilar parameters

Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. Here we present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.

42 ENGINEERING↗
DOE OSTI2021

Evaluating Diurnal and Semi-Diurnal Cycle of Precipitation in CMIP6 Models Using Satellite- and Ground-Based Observations

We report the diurnal and semi-diurnal cycle of precipitation simulated from CMIP6 models during 1996-2005 are evaluated globally between 60°S and 60°N, as well as at ten selected locations representing three categories of diurnal cycle of precipitation: (1) afternoon precipitation over land, (2) early morning precipitation over ocean, and (3) nocturnal precipitation over land. Three satellite-based and two ground-based rainfall products are used to evaluate the climate models. Globally, the ensemble mean of CMIP6 models shows a diurnal phase of 3 to 4 hours earlier over land and 1 to 2 hours earlier over ocean, when compared with the latest satellite products. These biases are in line with what were found in previous versions of climate models but reduced compared to the CMIP5 ensemble mean. Analysis at the selected locations complimented with in-situ measurements further reinforces these results. Several CMIP6 models have shown a significant improvement in the diurnal cycle of precipitation compared to their CMIP5 counterparts, notably on delaying afternoon precipitation over land. This can be attributed to the use of more sophisticated convective parameterizations. Most models are still unable to capture the nocturnal peak associated with elevated convection and propagating mesoscale convective systems, with a few exceptions that allow convection to be initiated above the boundary layer to capture nocturnal elevated convection. We also quantify an encouraging consistency between the satellite- and ground-based precipitation measurements despite differing spatiotemporal resolutions and sampling periods, which provides confidence in using them to evaluate the diurnal and semi-diurnal cycle of precipitation in climate models.

58 GEOSCIENCES↗
DOE OSTI2024

A Bayesian Multi-fidelity Neural Network to Predict Nonlinear Frequency Backbone Curves

The use of structural mechanics models during the design process often leads to the development of models of varying fidelity. Often low-fidelity models are efficient to simulate but lack accuracy, while the high-fidelity counterparts are accurate with less efficiency. Here, this paper presents a multi-fidelity surrogate modeling approach that combines the accuracy of a high-fidelity finite element model with the efficiency of a low-fidelity model to train an even faster surrogate model that parameterizes the design space of interest. The objective of these models is to predict the nonlinear frequency backbone curves of the Tribomechadynamics Research Challenge benchmark structure which exhibits simultaneous nonlinearities from frictional contact and geometric nonlinearity. The surrogate model consists of an ensemble of neural networks that learn the mapping between low and high-fidelity data through nonlinear transformations. Bayesian neural networks are used to assess the surrogate model's uncertainty. Once trained, the multi-fidelity neural network is used to perform sensitivity analysis to assess the influence of the design parameters on the predicted backbone curves. Additionally, Bayesian calibration is performed to update the input parameter distributions to correlate the model parameters to the collection of experimentally measured backbone curves.

42 ENGINEERING↗