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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2022

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

DOE OSTI2024

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Novel strategies for modal-based structural material identification

Here, we present modal-based methods for model calibration in structural dynamics, and address several key challenges in the solution of gradient-based optimization problems with eigenvalues and eigenvectors, including the solution of singular Helmholtz problems encountered in sensitivity calculations, non-differentiable objective functions caused by mode swapping during optimization, and cases with repeated eigenvalues. Unlike previous literature that relied on direct solution of the eigenvector adjoint equations, we present a parallel iterative domain decomposition strategy (Adjoint Computation via Modal Superposition with Truncation Augmentation) for the solution of the singular Helmholtz problems. For problems with repeated eigenvalues we present a novel Mode Separation via Projection algorithm, and in order to address mode swapping between inverse iterations we present a novel Injective mode ordering metric. We present the implementation of these methods in a massively parallel finite element framework with the ability to use measured modal data to extract unknown structural model parameters from large complex problems. A series of increasingly complex numerical examples are presented that demonstrate the implementation and performance of the methods in a massively parallel finite element framework [7], [5], using gradient-based optimization techniques in the Rapid Optimization Library (ROL) [21].

36 MATERIALS SCIENCE↗

DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗

DOE OSTI2020

(U) PRAD0697 & PRAD0698: Complex Loading of CeO₂ Powder

Cerium(IV) oxide (CeO₂) powder is shock compressed using the Precision High Energy-density Liner Implosion eXperiment (PHELIX) platform. Experimental results are compared against several modeling approaches. Compaction behavior is best captured with a P-∝ model, which calculates CeO₂ powder bulk densities within 80-99% of experimental values but overpredicts densi cation at the cylindrical target's outer radius and center by up to 20%. Preliminary calculations suggest that accuracy could be increased with the inclusion of a coupled strength model. Several common computational modeling approaches for the shock compression response of granular materials and the magnetohydrodynamic (MHD) force upon the impactor/liner in pulsed power compression experiments are investigated and analyzed for their validity. The Bi-linear Ramp, P-∝ PACXP, and P-∝ Menikoff-Kober continuum compaction models are calibrated to planar impact Hugoniot data for CeO₂ powder and used to predict the powder's shock compaction response under non-planar shock wave compression. MHD calculations of the PHELIX pulsed power driver are performed using an idealized resistor-inductor-capacitor (RLC) circuit calibrated to previous experiments. All simulations are performed using the LANL code FLAG. Two validation experiments are computationally designed using the calibrated compaction and circuit models, executed using the PHELIX platform on CeO₂ targets with initial porous densities of 3.95 and 4.03 g/cm³, measured with proton radiography, and analyzed against the model predictions. The two P-∝ models more accurately describe CeO₂ powder densi cation than the Bi-linear Ramp model. However, the two P-∝ models overpredict bulk density of the shock compressed CeO₂ powder by up to 20% when the appropriate impact velocities are applied. MHD calculations for both validation experiments underpredict liner impact velocities by 4-11% when using the idealized RLC circuit model calibrated to previous experiments. Compensating underpredictions of impact velocity and overpredictions of powder densication lead to a false accuracy in pre-shot calculations compared to experimental data. To improve correlation between simulations and experiments, the following improvements are suggested: 1. A coupled strength model for CeO₂ powder that updates strength as a function of porosity and applied stress. 2. An improved MHD circuit model that more accurately captures the PHELIX machine.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Can Restoring Tidal Wetlands Reduce Estuarine Nuisance Flooding of Coasts Under Future Sea‐Level Rise?

Wetland restoration is an increasingly popular nature‐based method for flood risk mitigation in coastal communities. In this study, we present a novel method using hydrodynamic modeling and harmonic analysis to quantify wetlands' ability to reduce future nuisance flooding. The method leverages a hydrodynamic model calibrated to present day data and was run for a range of future sea‐level rise (SLR) and wetland restoration scenarios to quantify changes to tidal harmonic amplitudes and phases. The harmonic constituents are used to generate water surface elevations over a time period of interest (e.g., one year) and compared to critical exceedance thresholds such as levee elevations. Then, changes to nuisance flooding are calculated by counting the number of hours critical thresholds are exceeded under different SLR and wetland restoration scenarios. We applied the method to Coos Bay, Oregon, USA as a test case. We found restoration reduces the number of hours nuisance flooding occurs in downtown Coos Bay from 15 hr (present day conditions) to 0 hr (fully restored condition) under median SLR (82 cm by 2100). Restoration had spatially variable impacts on reducing peak flood elevations with minimal impacts near the estuary mouth and greatest impact 32 km inland. The effectiveness of restoration was heavily dependent on future SLR. Restoration was maximally effective in 2050 under all SLR scenarios, less effective in 2100 under median SLR, and not effective under high SLR. Modeling results suggest increased tidal prism and accommodation space are driving restoration‐associated reductions in tidal amplitudes.

Brand, Matthew W. [Louisiana State Univ., Baton Ro↗

DOE OSTI2022

Calibration of a mesoscale tritium transport model for ceramic breeder materials in TMAP8 using experimental data

Due to the scarcity of the long-term external tritium supply, lithium-containing breeder materials are used in the blanket of fusion reactors to produce tritium faster than the deuterium-tritium fusion reaction consumes it. A cellular breeder is a promising material with higher lithium density and thermal conductivity than conventional ceramic pebbles and enhanced tritium extraction with interconnected pores. Although several experimental and modeling efforts have improved our understanding of tritium transport, the exact mechanisms governing tritium release from breeder materials are still unknown. As a result, current models cannot reliably assess the tritium breeding capabilities of cellular breeder materials. This work presents a tritium transport multiphysics model and calibrates it using experimental data of deuterium absorption and desorption from cellular breeder materials at different temperatures and pressures. This model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. After calibrating the model and performing sensitivity analysis, we discuss pre-dominant mechanisms governing tritium release from cellular breeder materials. The model is part of the Tritium Migration Analysis Program [TMAP8], a multiscale, multiphysics framework for tritium transport based on the finite element multiphysics framework MOOSE. This study demonstrates some of TMAP8’s capabilities and provides insight into the mechanisms governing tritium transport in cellular ceramic breeder materials.

08 HYDROGEN↗

DOE OSTI2024

Inverter Model Validation and Calibration Using Phasor Measurement Unit Data

As the penetration of inverter-based renewable energy resources increases in the power grid, especially at the distribution and microgrid levels, the need to accurately represent them in planning studies increases as well. However, due to the lack of well-established standard procedures, and vendor reluctance towards the detailed sharing of proprietary models, automated dynamic model validation and parameter calibration tools for inverter based resources (IBRs) remain scarce. This work presents a model validation and parameter calibration platform for representing IBRs with generic phasor-domain models. Phasor measurements of power system events are used for continuous validation using the data playback method, and model parameters are re-calibrated if a significant mismatch between measurements and model response is observed. Unique features of the proposed platform include- (a) an iterative Bayesian optimization approach towards parameter calibration to address a possible mismatch between the structures of generic models implemented in simulation softwares and actual commercial inverters, (b) error metrics designed to account for a possible mismatch between the time resolution of simulation and measurements, and (c) analysis of the measurement-simulation mismatch to provide guidance to engineering personnel regarding model shortcomings. The performance of the platform has been illustrated using both simulated data and field measurements to validate/calibrate inverter models in GridLAB-D.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Mapping use cases and dataset needs for benchmarking buildings data

A perennial challenge in buildings research is the lack of high-quality datasets that can be relied upon for a wide array of tasks, including model calibration and improving energy efficiency and load flexibility. Instrumenting a building for data collection is resource intensive, so it is important to be methodical in the approach and ensure that resulting data are flexible and useful for a broad range of analyses. This study aims to fill the gaps in characterizing potential use cases for buildings datasets and mapping them to dataset needs using a well-defined data infrastructure. Here, we have developed a systematic mapping strategy between buildings dataset needs and use cases to help streamline the processes of efficiently targeting datasets, designing building sensing systems, and determining buildings research use cases. We selected 14 prospective use cases and 11 refined buildings data categories for developing the preliminary dataset-needs-to-use-cases mapping matrix (‘DN-UC mapping matrix’) with generic ‘Tags’—a detailed sub-level of data categories extracted by justifying the needs of an aspect of the datasets to use cases. We present two example applications of the developed mapping matrix to demonstrate use of the mapping matrix and its effectiveness.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Rapid and high-throughput determination of sorghum ( Sorghum bicolor ) biomass composition using near infrared spectroscopy and chemometrics

Compositional characterization of biomass is vital for the biofuel industry. Traditional wet chemistry-based methods for analyzing biomass composition are laborious, time-consuming, and require extensive use of chemical reagents as well as highly skilled personnel. In this study, near-infrared (NIR) spectroscopy was used to quickly assess the composition of above-ground vegetative biomass from 113 diverse, photoperiod-sensitive, biomass-type sorghum (Sorghum bicolor) accessions cultivated under field conditions in Central Illinois. Biomass samples were analyzed using NIR spectra collected in the spectral range of 867–2536 nm, with their chemical compositions determined following the National Renewable Energy Laboratory (NREL) protocol. Advanced spectral pre-treatment and band selection techniques were utilized to develop calibration models using partial least squares regression (PLSR). The models’ effectiveness was assessed through cross-validation and independent data tests. The predictions for moisture, ash, extractives, glucan, xylan, acid-soluble lignin (ASL), acid-insoluble lignin (AIL), and total lignin were accurate and reliable, demonstrating the capability of NIR spectroscopy to provide rapid and precise characterization of sorghum biomass. The results demonstrated that NIR spectroscopy is an efficient tool for rapidly characterizing sorghum biomass, making it a sustainable option for screening desirable feedstock for biofuel or bioproduct production.

09 BIOMASS FUELS↗

DOE OSTI2022

Operability thresholds for thermally damaged EBW detonators

Operability thresholds that differentiate between functional RP-87 exploding bridge wire (EBW) detonators and nonfunctional RP-87 EBW detonators (duds) were determined by measuring the time delay between initiation and early wall movement (function time). The detonators were inserted into an externally heated hollow cylinder of aluminum and fired with current flow from a charged capacitor using an exploding bridge wire (EBW initiated). Functioning detonators responded like unheated pristine detonators when the function time was 4 μs or less. The operability thresholds of the detonators were characterized with a simple decomposition cookoff model calibrated using a modified version of the Sandia Instrumented Thermal Ignition (SITI) experiment. These thresholds are based on the calculated state of the PETN when the detonators fire. The operability threshold is proportional to the positive temperature difference (ΔT) between the maximum temperature within the PETN and the onset of decomposition (~406 K). Furthermore, the temperature difference alone was not sufficient to define the operability threshold. The operability threshold was also proportional to the time that the PETN had been at elevated temperatures. That is, failure was proportional to both temperature and reaction rate. The reacted gas fraction is used in the current work for the reaction correlation. Melting of PETN also had a significant effect on the operability threshold. Detonator failure occurred when the maximum temperature exceeded the nominal melting point of PETN (414 K) for 45±5 s or more.

42 ENGINEERING↗

DOE OSTI2022

Effect of lot microstructure variations on detonation performance of the triaminotrinitrobenzene (TATB)-Based insensitive high explosive PBX 9502

PBX 9502 is an important insensitive high explosive due to its combination of safety properties and detonation performance. It is a polymer-bonded formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) as the high explosive crystal, bound with Kel F-800 (FK-800), a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Two different types are used, one known as virgin that uses only pristine manufactured TATB, and a second known as recycled that has 50 wt.% of its TATB reclaimed from machining scraps of previously pressed virgin PBX 9502. Recycled lots have a higher percentage of fine particles compared to virgin lots, due to the fracturing and damage sustained by TATB crystals during pressing. We examine the influence that this TATB microstructure difference between virgin and recycled PBX 9502 has on detonation performance properties. New rate-stick geometry diameter effect, detonation front shape and cylinder expansion test data are obtained for two previously uncharacterized virgin lots and one recycled lot of PBX 9502. This is combined with previously published data for one virgin and one recycled lot for evaluation purposes. Detonation shock dynamics model calibrations are conducted on each of the five lots to provide an assessment of the detonation timing characteristics of virgin versus recycled PBX 9502 lots. For two of the virgin lots, detonations propagate slower in the rate-stick geometry than those in the two recycled lots. However, the other virgin lot tested has propagation rates comparable to that of the recycled lots for larger diameter rate-sticks. The latter result, though, is shown to be geometry dependent and depends on the range of detonation curvatures accessed in different geometries. New copper-confined, cylinder expansion tests are conducted on each of the three virgin and two recycled lots to obtain detonation product Jones-Wilkins-Lee equations of state, enabling an assessment of the metal push capabilities for each of the five lots. We find that the metal push capabilities, characterized by the evolution of the heat of detonation with volume, are similar between the virgin and recycled lots. Thus, changes in microstructure between different PBX 9502 lots seemingly affect the rate of reaction, but not the overall energy content.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Developing and tuning a community scale energy model for a disadvantaged community

This work describes the development of a community-scale energy model for a mixed-use low-income community located in Huntington Beach, CA. An accurate community-scale energy model is useful for evaluating the use of limited capital resources used to invest in clean energy technologies. This work lays out the process of developing such a model while relying primarily on publicly available data and highlighting critical partnerships necessary for model development success. The primary contribution of this work is the demonstration of the process used to develop an accurate energy model for a disadvantaged community when minimal building and energy use data is available. The heart of the model is the physics-based community scale energy modeling platform URBANopt. Using a bottom-up load modeling approach, energy simulated energy use falls within 3% or less of aggregate annual utility data, and within 10% or less aggregate monthly utility data. The demonstrated model development and tuning process can be used by others to characterize other atypical communities, which may differ significantly from prototypical models.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2023

The physics of desiccation cracks 1: Ductile fracturing and dependence on relative humidity

Ductile deformation is ubiquitously found in the shrinkage of geomaterials. The existence of ductility requires elastoplastic mechanics when analyzing the structure deformation under external stress. Here, in this work, we explore the physics of ductile fracturing based on the results from desiccation experiments and triaxial tests. By using the digital image correlation (DIC) method to generate the strain maps of samples undergoing desiccation cracking under different relative humidities, we obtain results showing the previously postulated Cnoidal Wave ductile failure patterns propagating under atmospheric condition-controlled crack velocities. We then correlate these observations with rate-dependent plasticity models calibrated through triaxial tests undergoing several unloading–reloading cycles and velocity stepping. This work demonstrates the necessity to consider ductility in soil cracking, indicating that the formation of crack patterns in soil desiccation is a slow and predictable process.

58 GEOSCIENCES↗

DOE OSTI2021

Low global warming potential (GWP) refrigerant supermarket refrigeration system modeling and its application

As an environmentally friendly low global warming potential (GWP) refrigerant, Carbon dioxide (CO 2 ) has continuously gained popularity and research attention as alternative refrigerant for supermarket refrigeration system. In this paper, to fulfill the increasing need of accurate Low-GWP supermarket refrigeration models for development of supervisory level control and optimization strategies, a high fidelity model is developed for CO 2 transcritical supermarket refrigeration system which includes compressor rack of low temperature (LT) compressors and medium temperature (MT) compressors, air-cooled gas cooler, evaporator, expansion valves and other auxiliary equipment. A resistance-capacity model structure is proposed to simulate the display cases. Semi-thermodynamic models are proposed to estimate reciprocating compressors volumetric efficiency and power consumption. Furthermore the zone modeling approach is used for evaporator simulation, and air-cooled gas cooler is modeled with distributed modeling method. The expansion valve simulation is based on orifice flow model. To calibrate these models, both manufacture performance data and experimental data are used. The experiments are conducted with a full instrumental CO 2 supermarket refrigeration system installed in Oak Ridge National Lab Flexible Research Platform (FRP). The simulation model can predict the system performance, including power consumption, cooling capacity, mass flow rate, temperature, and pressure, with high accuracy (within ±4%) compared to experimental data. In addition, this developed model has been used to create the system optimum high side pressure for high side expansion valve control, and to generate wide operating range simulation data for developing the batter-equivalent commercial refrigeration system model which can be used in grid interactive control development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Non-destructive measurement and real-time monitoring of apple hardness during ultrasonic contact drying via portable NIR spectroscopy and machine learning

Portable near-infrared spectrometer in the spectral range of 900–1700 nm was evaluated for the first time to assess and monitor apple hardness in real-time during ultrasonic drying. Calibration models were developed using PLS and ANN, and their performances were evaluated by internal leave-one-out cross-validation and an external dataset. Several pre-treatments including standard normal variate (SNV), multiplicative scatter correction (MSC), Savitzky–Golay first and second derivatives were employed to examine the effects of spectral variations in hardness prediction. Seven important wavelengths were selected using weighted regression coefficients to develop a simple MLR model to facilitate the model interpretation and circumvent noise. The models using PLS, MLR, and ANN with selected wavelengths predicted the apple hardness with R 2 p of 0.91, 0.91, 0.95, and RMSEP of 14.78, 14.85, and 12.46 N, respectively. Finally, the results indicate that portable NIR spectrometers are quite promising for real-time monitoring of apple hardness during ultrasonic drying.

47 OTHER INSTRUMENTATION↗

DOE OSTI2021

Geometry-based, Gaussian profile model for optical knife-edge displacement sensor

This article describes an analytical model for the well-known optical knife-edge displacement sensor that predicts the relationship between the knife-edge displacement and photointerrupter output. The model includes: 1) the geometry of a rectangular knife-edge profile penetrating a circular laser beam perpendicular to the beam axis; and 2) the Gaussian intensity profile for the beam cross-section. A comparison between the model predictions and measured photointerrupter for a linear range of knife-edge displacements is provided. It is shown that good agreement is obtained without model calibration.

42 ENGINEERING↗