Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameter 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 217 records · Page 12

Evaluation of the Tau-Omega Model Over Bare and Wheat-Covered Flat and Periodic Soil Surfaces at P- and L-band

It has been over ten years since the successful launch of the first-ever dedicated satellite for global soil moisture monitoring; Soil Moisture and Ocean Salinity (SMOS). Looking towards the future, P-band (0.3-1 GHz) is a promising technique to replace or enhance the L-band (1.4 GHz) SMOS and SMAP (Soil Moisture Active Passive) missions because of an expected reduction in roughness and vegetation impact, leading to an improved soil moisture accuracy over rougher soil surfaces and more densely vegetated areas. Accordingly, this investigation evaluated the tau-omega model at P-band (0.75 GHz) using a tower-based experiment in Victoria, Australia, where brightness temperature observations were collected concurrently at P- and L-band over bare and wheat-covered flat and periodic soil surfaces. The potential to retrieve soil moisture without discriminating periodic and flat surfaces was investigated by applying the roughness and vegetation parameters calibrated for flat soil to retrieve the moisture of periodic soil. Results showed that P-band had a comparable RMSE across different roughness configurations (variations less than 0.016 m3/m3) for both bare and wheat-covered soil, while the L-band RMSE was only comparable for wheat-covered soil, indicating that periodic surfaces did not need to be discriminated in such scenarios. Conversely, a difference of 0.022 m3/m3 was observed for L-band with bare soil. A reduced vegetation impact was also demonstrated at P-band, with an RMSE of 0.029 m3/m3 achieved when completely ignoring the wheat existence with under 4-kg/m2 vegetation water content, whereas at L-band the RMSE increased to 0.063 m3/m3. This study therefore paves the way for a successful P-band radiometer mission for obtaining more accurate global soil moisture information.

P-band↗

Non-contact velocity measurement instruments and systems, and related methods

Methods for calibration of non-contact velocity measurements and systems for implementing the same are described. Generally, the method comprises inducing a shock wave into a sample at a stress intensity that varies across the sample's elastic limit, which corresponds to the elastic-plastic state transition of the sample. That transition state may be at the sample's Hugoniot elastic limit. The velocity of the sample is measured using a non-contact velocity measurement instrument such as a velocimeter. The measurement may be compared to a predicted velocity or a velocity measurement made by another system to determine calibration parameters.

47 OTHER INSTRUMENTATION↗

Gaussian process for calibration and control of GlueX Central Drift Chamber

We have developed and implemented a machine learning based system to calibrate and control the GlueX Central Drift Chamber at Jefferson Lab, VA, in near real-time. The system monitors environmental and experimental conditions during data taking and uses those as inputs to a Gaussian process (GP) with learned prior. The GP predicts calibration constants in order to recommend a high voltage (HV) setting for the detector that maintains consistent detector performance (gain and resolution) throughout data taking. This approach is in stark contrast to traditional detector operations in which the detector operates at fixed HV and its calibration parameters vary quite considerably with time. Additionally, the ML based system utilizes uncertainty quantification to correct the recommended control parameters when appropriate. We will present results from the ML system autonomously during the Charged Pion Polarizability (CPP) experiment conducted in Hall D at Jefferson Lab.

McSpadden, Helen↗

Calibration of parallel bond parameters in bonded particle models via physics-informed adaptive moment optimisation

This study proposes an automated calibration procedure for bond parameters in bonded discrete element modelling. By exploring the underlying physical correlations between microscopic parameters of bonds and macroscopic strength parameters of the continuum to be modelled, the microscopic shear strength and tensile strength are identified as independent variables for calibration purpose. Then a physics-informed iterative scheme is proposed to automatically approximate the bond parameters by viewing the micro-macro relation as an implicitly defined mathematical mapping function. As a result of highly non-convex features of this implicit mapping, the adaptive moment estimation (Adam), which is especially suitable for problems with noisy gradients, is adopted as the basic iterative scheme, in conjunction with other numerical techniques to approximately evaluate the partial derivatives involved. The whole procedure offers a simple and effective framework for bond parameter calibration. Finally, a numerical example of SiC ceramic is provided for validation. By compared with some existing calibration methods, the proposed method shows significant advantages in terms of calibration efficiency and accuracy.

42 ENGINEERING↗

Improving 3D reconstruction quality for root phenotyping: assessing the impact of camera calibration and imaging parameters

Arate 3D reconstruction is essential for high-throughput plant phenotyping, particularly for studying complex structures such as root systems. While photogrammetry and Structure from Motion (SfM) techniques have become widely used for 3D root imaging, the camera settings used are often underreported in studies, and the impact of camera calibration on model accuracyccu remains largely underexplored in plant science. In this study, we systematically evaluate the effects of focus, aperture, exposure time, and gain settings on the quality of 3D root models made with a multi-camera scanning system. We show through a series of experiments that calibration significantly improves model quality, with focus misalignment and shallow depth of field (DoF) being the most important factors affecting reconstruction accuracy. Our results further show that proper calibration has a greater effect on reducing noise than filtering it during post-processing, emphasizing the importance of optimizing image acquisition rather than relying solely on computational corrections. This work improves the repeatability and accuracy of 3D root imaging for phenotyping pipelines by giving useful calibration guidelines. This leads to better trait quantification for use in crop research and plant breeding in downstream analysis.

3D reconstruction↗

Optimizing process-based models to predict current and future soil organic carbon stocks at high-resolution

From hillslope to small catchment scales (< 50 km 2 ), soil carbon management and mitigation policies rely on estimates and projections of soil organic carbon (SOC) stocks. Here we apply a process-based modeling approach that parameterizes the MIcrobial-MIneral Carbon Stabilization (MIMICS) model with SOC measurements and remotely sensed environmental data from the Reynolds Creek Experimental Watershed in SW Idaho, USA. Calibrating model parameters reduced error between simulated and observed SOC stocks by 25%, relative to the initial parameter estimates and better captured local gradients in climate and productivity. The calibrated parameter ensemble was used to produce spatially continuous, high-resolution (10 m 2 ) estimates of stocks and associated uncertainties of litter, microbial biomass, particulate, and protected SOC pools across the complex landscape. Here, subsequent projections of SOC response to idealized environmental disturbances illustrate the spatial complexity of potential SOC vulnerabilities across the watershed. Parametric uncertainty generated physicochemically protected soil C stocks that varied by a mean factor of 4.4 × across individual locations in the watershed and a – 14.9 to + 20.4% range in potential SOC stock response to idealized disturbances, illustrating the need for additional measurements of soil carbon fractions and their turnover time to improve confidence in the MIMICS simulations of SOC dynamics.

54 ENVIRONMENTAL SCIENCES↗

The final COS-B database: In-flight calibration of instrumental parameters

A method for the determination of temporal variation of sensitivity is designed to find a set of parameters which lead to maximum consistency between the intensities derived from different observation periods. This method is briefly described and the resulting sensitivity and background variations presented.

Strong, A. W.↗

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗

A rigorous physics-based enhanced parameter estimation (EPE) methodology for calibration of building energy simulations

Buildings rarely perform as designed/simulated and there are numerous tangible benefits if this gap is reconciled. A new scientifically rigorous yet pragmatic methodology for calibrating building energy simulations - called Enhanced Parameter Estimation (EPE) - is proposed that allows physically relevant parameter estimation rather than a blind force-fit to energy use data. Starting with a rapidly created simulation model, calibration is performed in two stages: (a) building shell calibration with the HVAC system replaced by an ideal system that meets the loads (b) HVAC system calibration with the building shell and all internal loads replaced by a box with only process loads. In the first stage, EPE identifies a small number of high-level heat flows in the energy balance, calculates them with specifically tailored individual driving functions, introduces physically significant parameters to best accomplish energy balance, and, estimates the parameters and their uncertainty bounds. Calibration is thus done with corrective heat flows without any arbitrary tuning of input parameters. Calibration accuracy is enhanced by machine learning of the residual errors. The EPE methodology is demonstrated by means of: a synthetic building and an actual 75,000 sq. ft. building in Pennsylvania. Future work needed for widespread application is discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Determining Levels of Detail for Simulators of Parallel and Distributed Computing Systems via Automated Calibration

There are two sources of inaccuracy when simulating parallel and distributed computing systems: (i) a simulator implemented at an insufficient level of detail; and (ii) incorrectly calibrated simulation parameter values. Increasing the simulator’s level of detail can improve accuracy, but at the cost of higher space, time, and/or software complexity. Furthermore, evaluating the intrinsic accuracy of a simulator requires that its parameters be well-calibrated. Making decisions regarding the level of detail is thus challenging. We propose a methodology for instantiating the simulation calibration process and a framework for automating this process, which makes it possible to pick appropriate levels of detail for any simulator. We demonstrate the usefulness of our approach via two case studies for two different domains.

McDonald, Jessie [University of Hawaii at Manoa, H↗

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗

Long-term measurements in a mixed-grass prairie reveal a change in soil organic carbon recalcitrance and its environmental sensitivity under warming

Soil respiration, the major pathway for ecosystem carbon (C) loss, has the potential to enter a positive feedback loop with the atmospheric CO2 due to climate warming. For reliable projections of climate-carbon feedbacks, accurate quantification of soil respiration and identification of mechanisms that control its variability are essential. Process-based models simulate soil respiration as functions of belowground C input, organic matter quality, and sensitivity to environmental conditions. However, evaluation and calibration of process-based models against the long-term in situ measurements are rare. Here, we evaluate the performance of the Terrestrial ECOsystem (TECO) model in simulating total and heterotrophic soil respiration measured during a 16-year warming experiment in a mixed-grass prairie; calibrate model parameters against these and other measurements collected during the experiment; and explore whether the mechanisms of C dynamics have changed over the years. Calibrating model parameters against observations of individual years substantially improved model performance in comparison to pre-calibration simulations, explaining 79–86% of variability in observed soil respiration. Interannual variation of the calibrated model parameters indicated increasing recalcitrance of soil C and changing environmental sensitivity of microbes. Overall, we found that (1) soil organic C became more recalcitrant in intact soil compared to root-free soil; (2) warming offset the effects of increasing C recalcitrance in intact soil and changed microbial sensitivity to moisture conditions. Furthermore, these findings indicate that soil respiration may decrease in the future due to C quality, but this decrease may be offset by warming-induced changes in C cycling mechanisms and their responses to moisture conditions.

59 BASIC BIOLOGICAL SCIENCES↗