Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer 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.

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Digital Twins for Materials

Digital twins are emerging as powerful tools for supporting innovation as well as optimizing the in-service performance of a broad range of complex physical machines, devices, and components. A digital twin is generally designed to provide accurate in-silico representation of the form (i.e., appearance) and the functional response of a specified (unique) physical twin. This paper offers a new perspective on how the emerging concept of digital twins could be applied to accelerate materials innovation efforts. Specifically, it is argued that the material itself can be considered as a highly complex multiscale physical system whose form (i.e., details of the material structure over a hierarchy of material length) and function (i.e., response to external stimuli typically characterized through suitably defined material properties) can be captured suitably in a digital twin. Accordingly, the digital twin can represent the evolution of structure, process, and performance of the material over time, with regard to both process history and in-service environment. This paper establishes the foundational concepts and frameworks needed to formulate and continuously update both the form and function of the digital twin of a selected material physical twin. The form of the proposed material digital twin can be captured effectively using the broadly applicable framework of n-point spatial correlations, while its function at the different length scales can be captured using homogenization and localization process-structure-property surrogate models calibrated to collections of available experimental and physics-based simulation data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quench process modeling & simulation in the heat-treatment of critical aerospace components

To manufacture light-weight, advanced metal alloy components for gas turbine engines, quench heat-treatment processes are typically used. By quenching the component from elevated tempera-tures, the alloy sometimes undergoes a solid-state phase transformation which produces special microstructures with the required, enhanced mechanical properties. However, the quenching can also lead to cracks forming in the component. Addressing the quench cracking problems adds a significant burden to the cost, schedule, and energy demand of manufacture. Currently, optimizing the quench process to mitigate or avoid the cracking is performed largely by trial-and-error, relying heavily on costly experimental (thermocouple) trials to understand the local thermal gradients which cause the cracks to form. In this first part (Phase 1) of the work, high-performance computing is employed to establish the ability of modern CFD (computational fluid dynamics) to alleviate or wholly replace the experimental quenching trials by virtual testing. A baseline CFD model is defined and its accuracy established to be comparable to (and which usually exceeds) the accuracy of existing HTC (heat-transfer correlation) based simulation methods of quenching. As a first-principles based approach, “calibration” of the Baseline CFD model is independent of the quench process itself, but instead relies on the accuracy of the underlying (modeled), generic two-phase fluid processes which cannot be currently resolved by CFD for large, industrial-scale cases. A novel, high-fidelity DNS capability has been developed and verified to examine and further improve upon the mean-field closure submodels on which the Baseline CFD approach is based.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Neural Density Estimation and Uncertainty Quantification for ChemCam Spectra [Slides]

The ChemCam instrument of Curiosity uses laser-induced breakdown spectroscopy (LIBS). It fires a laser at target and vaporizes rock surfaces, creating a plasma. Three spectrographs divide the plasma light into wavelengths for chemical analysis: ultraviolet, violet, and visible near-infrared. Regression methods (SVR, PCR, CNN) have been employed for calibration (prediction of the elemental composition of samples); however, labeled ChemCam samples are limited. Here, we focus on unsupervised learning and employ generative models from ChemCam analysis. Further, we use labels (supervised) in combination to the generative model to compute uncertainties related to predictions. We report generative modeling can be successfully applied to model real-world data. Normalizing flow models can be efficiently constructed on latent spaces for fast downstream inference. Unsupervised and supervised learning can be combined to form an uncertainty quantification framework.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A permafrost implementation in the simple carbon–climate model Hector v.2.3pf

Abstract. Permafrost currently stores more than a fourth of global soil carbon. A warming climate makes this carbon increasingly vulnerable to decomposition and release into the atmosphere in the form of greenhouse gases. The resulting climate feedback can be estimated using land surface models, but the high complexity and computational cost of these models make it challenging to use them for estimating uncertainty, exploring novel scenarios, and coupling with other models. We have added a representation of permafrost to the simple, open-source global carbon–climate model Hector, calibrated to be consistent with both historical data and 21st century Earth system model projections of permafrost thaw. We include permafrost as a separate land carbon pool that becomes available for decomposition into both methane (CH4) and carbon dioxide (CO2) once thawed; the thaw rate is controlled by region-specific air temperature increases from a preindustrial baseline. We found that by 2100 thawed permafrost carbon emissions increased Hector’s atmospheric CO2 concentration by 5 %–7 % and the atmospheric CH4 concentration by 7 %–12 %, depending on the future scenario, resulting in 0.2–0.25 ∘C of additional warming over the 21st century. The fraction of thawed permafrost carbon available for decomposition was the most significant parameter controlling the end-of-century temperature change in the model, explaining around 70 % of the temperature variance, and was distantly followed by the initial stock of permafrost carbon, which contributed to about 10 % of the temperature variance. The addition of permafrost in Hector provides a basis for the exploration of a suite of science questions, as Hector can be cheaply run over a wide range of parameter values to explore uncertainty and can be easily coupled with integrated assessment and other human system models to explore the economic consequences of warming from this feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗
DOE OSTI2020

Improvements to CTF Closure Models for Modelling of Two-Phase Flow

This report documents the efforts to improve the CTF prediction of void fraction and two-phase pressure drop by improving two-phase closure models. Previous validation activities have revealed that CTF tends to overpredict void fraction and two-phase pressure drop. In response to this, two approaches were taken to improve CTF’s predictive capabilities. Based on findings that the interfacial drag and subcooled boiling models significantly impact void prediction, alternative closure models for these physical effects were found and implemented into the code. An extensive assessment was performed by using the existing and newly added two-phase experimental data, which showed that void prediction, wall temperature, and two-phase pressure-drop results are improved by using the newly implemented models. A second approach for improving CTF modeling accuracy involved the use of a Bayesian calibration process that uses CTF validation data to optimize selected modeling coefficients and closure model multipliers to achieve a more accurate prediction of experimental results. A similar assessment was performed with calibrated models that improved the void and pressure-drop prediction for test cases that were both included and not included in the calibration dataset. Results of this study will be used to change the models used in CTF to achieve more accurate BWR analyses moving forward. This study also identified opportunities for improving additional closure models and other opportunities to use calibration techniques to improve CTF.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Four-year Cosmology Large Angular Scale Surveyor (CLASS) Observations: On-sky Receiver Performance at 40, 90, 150, and 220 GHz Frequency Bands

The Cosmology Large Angular Scale Surveyor (CLASS) observes the polarized cosmic microwave background (CMB) over the angular scales of 1° ≲ θ ≤ 90° with the aim of characterizing primordial gravitational waves and cosmic reionization. We report on the on-sky performance of the CLASS Q-band (40 GHz), W-band (90 GHz), and dichroic G-band (150/220 GHz) receivers that have been operational at the CLASS site in the Atacama desert since 2016 June, 2018 May, and 2019 September, respectively. We show that the noise-equivalent power measured by the detectors matches the expected noise model based on on-sky optical loading and lab-measured detector parameters. Using Moon, Venus, and Jupiter observations, we obtain power to antenna temperature calibrations and optical efficiencies for the telescopes. From the CMB survey data, we compute instantaneous array noise-equivalent-temperature sensitivities of 22, 19, 23, and 71 $\mu {{\rm{K}}}_{\mathrm{cmb}}\sqrt{{\rm{s}}}$ for the 40, 90, 150, and 220 GHz frequency bands, respectively. These noise temperatures refer to white noise amplitudes, which contribute to sky maps at all angular scales. Future papers will assess additional noise sources impacting larger angular scales.

Observational cosmology↗
DOE OSTI2023

Source Characterization of the Declared North Korean Nuclear Tests From Regional Distance Coda Wave Spectral Ratios

Abstract Seismic observations of underground nuclear explosions provide crucial data on source yield and depth that cannot easily be estimated from other geophysical methods. However, it is difficult to obtain reliable yield estimates for test sites for which we do not have direct seismic calibration experiments. To obtain source information from uncalibrated sites and paths, local and regional seismic records of six, proximal, declared underground nuclear explosions in North Korea are used to compute spectral ratios of narrow‐band waveform envelopes of body‐wave coda that remove path and site effects to reveal precise, relative source moment. The yields of these explosions are obtained from the observed source ratios by simultaneously fitting the classical source model of Mueller and Murphy (1971), https://doi.org/10.1785/bssa0610061675 to all event pairs. The source model provides an impressive fit to the observations considering that the P phase coda derived source spectral ratios did not require a prior knowledge of the source or regionally calibrated corrections to be applied to the data. However, the observed corner frequencies from S wave coda spectral ratios are lower than the source model predictions, but are well fit by models calculated using the corner frequency consistent with the Fisk conjecture. The results presented here provide novel constraints on the spectral distributions of the energy radiated by the sources of the DPRK test series, and allows for an independent evaluation of existing estimated source model parameters.

58 GEOSCIENCES↗
DOE OSTI2022

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.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

ReNew100: Reliable Power System Operation with 100% Renewable Generation

The ReNew100 project has developed and demonstrated an operator support system (OSS) to operate power systems with 100% renewable power generation from inverter-based resources (IBRs) like wind and solar that significantly reduces the risk of widespread power outages in a simulated operational environment at Technology Readiness Level (TRL) 6. The OSS achieves resilient operation for power systems under changing operation modes including varying combinations of conventional and renewable generation, including cases with 100% renewable generation from wind and solar IBRs. The OSS continuously monitors the N-1 security of the power system, which refers to the ability of the power systems to survive single outages of power system elements, like the loss of a generator or a power line, without causing widespread power outages beyond those expected and planned for. If the OSS identifies that the system is not N-1 secure, an automatic Controller Parameter Optimization (CoPO) developed within Renew100 is activated to optimize the N-1 security of the system by tuning controller parameters within different assets like batteries, wind, and solar plants. The core innovation of ReNew100 is the development of this CoPO as well as the demonstration of the OSS for a real-time simulation of the Hawai'i Island power system. Moreover, Renew100 has developed new model calibration techniques for dynamic power system models and demonstrated them for Hawai'i Island.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

ML-based calibration and control of the GlueX Central Drift Chamber

The GlueX Central Drift Chamber (CDC) in Hall D at Jefferson Lab, used for detecting and tracking charged particles, is calibrated and controlled during data taking using a Gaussian process. The system dynamically adjusts the high voltage applied to the anode wires inside the chamber in response to changing environmental and experimental conditions such that the gain is stabilized. Control policies have been established to manage the CDC's behavior. These policies are activated when the model's uncertainty exceeds a configurable threshold or during human-initiated tests during normal production running. Finally, we demonstrate the system reduces the time detector experts dedicate to calibration of the data offline, leading to a marked decrease in computing resource usage without compromising detector performance.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Optimization of thermal barrier coating performance and durability over a drive cycle

A methodology to thermo-mechanically optimize a piston thermal barrier coating over a full drive cycle was established. The optimization objective was to minimize the heat transfer to the engine wall while maintaining structural integrity of the coating. Over 800 candidate materials were investigated and the optimization required more than one million non-road transient drive cycle calculations; real materials were investigated to ensure a realizable result and the existence of thermal and mechanical properties. High computational efficiency was achieved using a recently developed analytical heat transfer technique for multilayer engine walls. An uncoupled approach was utilized for the optimization, wherein the gas temperature and heat transfer coefficient profiles from a fully coupled and calibrated baseline model over the 20-min drive cycle were employed. The coating/piston interface temperature was constrained to be below the maximum piston service temperature limit. The durability was assessed using a recently developed analytical coating delamination framework for engine in-cylinder coatings based on the energy release rate when a crack forms. Results are presented for a mechanically unconstrained optimization and for cases constrained to three fixed levels of drive-cycle maximum energy release rate, and also constrained by the individual material’s toughness. The best-performing coating materials identified were verified using the fully coupled system-level model, which compared well to the uncoupled predictions. A study on the effect of adding a sealing layer to some high-performing, but porous, coatings showed a reduction in fuel consumption benefit and an increased exhaust temperature over the cycle, but the system still outperformed the uncoated case. The results of the study elucidate the importance of including engine performance and mechanical failure considerations in thermal barrier coating design.

Engineering↗
DOE OSTI2022

Calibrating constitutive models with full‐field data via physics informed neural networks

Abstract The calibration of solid constitutive models with full‐field experimental data is a long‐standing challenge, especially in materials that undergo large deformations. In this paper, we propose a physics‐informed deep‐learning framework for the discovery of hyperelastic constitutive model parameterizations given full‐field surface displacement data and global force‐displacement data. Contrary to the majority of recent literature in this field, we work with the weak form of the governing equations rather than the strong form to impose physical constraints upon the neural network predictions. The approach presented in this paper is computationally efficient, suitable for irregular geometric domains, and readily ingests displacement data without the need for interpolation onto a computational grid. A selection of canonical hyperelastic material models suitable for different material classes is considered including the Neo–Hookean, Gent, and Blatz–Ko constitutive models as exemplars for general non‐linear elastic behaviour, elastomer behaviour with finite strain lock‐up, and compressible foam behaviour, respectively. We demonstrate that physics informed machine learning is an enabling technology and may shift the paradigm of how full‐field experimental data are utilized to calibrate constitutive models under finite deformations.

Hamel, Craig M.↗
DOE OSTI2021

How Equivalent Are Equivalent Porous Media?

Abstract Geoenergy and geoengineering applications usually involve fluid injection into and production from fractured media. Accounting for fractures is important because of the strong poromechanical coupling that ties pore pressure changes and deformation. A possible approach to the problem uses equivalent porous media to reduce the computational cost and model complexity instead of explicitly including fractures in the models. We investigate the validity of this simplification by comparing these two approaches. Simulation results show that pore pressure distribution significantly differs between the two approaches even when both are calibrated to predict identical values at the injection and production wells. Additionally, changes in fracture stability are not well captured with the equivalent porous medium. We conclude that explicitly accounting for fractures in numerical models may be necessary under some circumstances to perform reliable coupled thermohydromechanical simulations, which could be used in conjunction with other tools for induced seismicity forecasting.

58 GEOSCIENCES↗
DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2021

Simulation and Modeling of Time-Resolved X-Ray Detector for the Saturn Accelerator

In this work, we present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current–voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia’s Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn’s Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16×10 10 and 1.88×10 10 rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV’s non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗