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 145 records · Page 8

A study of shock initiation experiments for the explosive PBX 9502 using three reactive burn models

Shock to detonation transition (SDT) experiments are essential in calibrating and validating reactive burn models for explosives. This work investigates the large collection of SDT test data for the explosive PBX 9502 at ambient temperature that was presented by Gustavsen, Sheffield, and Alcon [Journal of Applied Physics, 99, 114907 (2006)]. We first analyze the experimental data and compare two different methods of determining the shock transition time/distance (namely, the bilinear method and the single-curve method). This reveals some of the uncertainty in estimating shock transition points, which contributes to scatter in Pop-plot data. Next, we compare the WSD, AWSD, and SURFplus reactive burn models for a collection of approximately 20 experimental shots using the FLAG hydrocode. Error estimates are used to quantify how well each reactive burn model (and their respective parameter calibrations) performs at predicting the SDT process for a range of loading conditions. Additionally, the importance of mesh resolution and numerical dissipation in SDT simulations will be assessed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Observational magnetometer calibration with the Hubble Space Telescope's new magnetometers

The two magnetometers recently replaced on the Hubble Space Telescope during the STS-61 Servicing Mission are now being used successfully for Coarse Attitude Determination during spacecraft vehicle safemode recovery operation. The magnetometer alignments relative to the spacecraft's vehicle's reference frame and the magnetic coupling of the sensors to the four magnetic torquer bars were determined. Coarse Attitude determination errors are now reduced to an average of 0.6 deg. Magnetometer Sensing System calibration and Coarse Attitude determination testing with the new calibration parameters is a geometrical problem. Telemetered earth magnetic field data was collected at twenty-six different vehicle attitudes. The spacecraft attitudes selected were distributed as widely apart as possible throughout the Geocentric Inertial Coordinate reference frame. It is also desirable to sample the Earth's magnetic field over as many different locations of the spacecraft's passage over the Earth as possible, within the limitation of the +/- 28.5 deg orbital inclination. A full range of magnetic moment outputs from the torquer bars needs to be sampled, +/- 3600 ampmeters squared, as well as data when the torquer bars have zero current. Graphic utilities were also developed to visually aid in optimizing the data collection process. Finally, a brief discussion of a method for collecting data for future calibrations is suggested.

Broude, Sidney M.↗

Reliable vision-guided grasping

Automated assembly of truss structures in space requires vision-guided servoing for grasping a strut when its position and orientation are uncertain. This paper presents a methodology for efficient and robust vision-guided robot grasping alignment. The vision-guided grasping problem is related to vision-guided 'docking' problems. It differs from other hand-in-eye visual servoing problems, such as tracking, in that the distance from the target is a relevant servo parameter. The methodology described in this paper is hierarchy of levels in which the vision/robot interface is decreasingly 'intelligent,' and increasingly fast. Speed is achieved primarily by information reduction. This reduction exploits the use of region-of-interest windows in the image plane and feature motion prediction. These reductions invariably require stringent assumptions about the image. Therefore, at a higher level, these assumptions are verified using slower, more reliable methods. This hierarchy provides for robust error recovery in that when a lower-level routine fails, the next-higher routine will be called and so on. A working system is described which visually aligns a robot to grasp a cylindrical strut. The system uses a single camera mounted on the end effector of a robot and requires only crude calibration parameters. The grasping procedure is fast and reliable, with a multi-level error recovery system.

Nicewarner, Keith E.↗

Reliable vision-guided grasping

Automated assembly of truss structures in space requires vision-guided servoing for grasping a strut when its position and orientation are uncertain. This paper presents a methodology for efficient and robust vision-guided robot grasping alignment. The vision-guided grasping problem is related to vision-guided 'docking' problems. It differs from other hand-in-eye visual servoing problems such as tracking in that the distance from the target is a relevant servo parameter. The methodology described in this paper is a hierarchy of levels in which the vision/robot interface is decreasingly 'intelligent', and increasingly fast. Speed is achieved primarily by information reduction. This reduction exploits the use of region-of-interest windows in the image plane and feature motion prediction. These reductions invariably require stringent assumptions about the image. Therefore, at a higher level, these assumptions are verified using slower, more reliable methods. This hierarchy provides for robust error recovery in that when a lower-level routine fails, the next-higher routine will be called and so on. A working system is described which visually aligns a robot to grasp a cylindrical strut. The system uses a single camera mounted on the end effector of a robot and requires only crude calibration parameters. The grasping procedure is fast and reliable, with a multi-level error recovery system.

Nicewarner, Keith E.↗

PACE OCI Polarization Sensitivity Based on Pre-launch Testing

The Ocean Color Instrument (OCI) is a sensor on the upcoming Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission, scheduled for launch in early 2024. OCI is a grating spectrometer with hyperspectral coverage from the ultraviolet (about 310 nm) to near-infrared (about 900 nm), with additional filtered channels in the short-wave infrared (940 nm – 2260 nm). This instrument will provide ocean color science data to continue the data sets collected by heritage sensors MODIS, SeaWiFs, and VIIRS, but with increased spectral coverage and improved accuracy. In order to achieve the high levels of accuracy demanded by the science community, a rigorous ground test program was conducted to calibrate the instrument and ensure that the calibration can be transferred to on-orbit operations. Some calibration parameters can only be measured during pre-launch testing; one such parameter is the polarization sensitivity. Polarization testing measured the Mueller matrix components needed to determine the polarization sensitivity for all spectral bands for a series of telescope scan angles covering the expected on-orbit scan range. Results indicate that the sensitivity is below 0.6 % except at the shortest wavelengths (less than 340 nm) and was characterized to better than 0.1 % above 340 nm. This indicates that any polarized scenes measured on orbit can be corrected for with a high degree of confidence.

PACE↗

PACE OCI Polarization Sensitivity Based on Pre-launch Testing

The Ocean Color Instrument (OCI) is a sensor on the upcoming Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission, scheduled for launch in early 2024. OCI is a grating spectrometer with hyperspectral coverage from the ultraviolet (about 310 nm) to near-infrared (about 900 nm), with additional filtered channels in the short-wave infrared (940 nm – 2260 nm). This instrument will provide ocean color science data to continue the data sets collected by heritage sensors MODIS, SeaWiFs, and VIIRS, but with increased spectral coverage and improved accuracy. In order to achieve the high levels of accuracy demanded by the science community, a rigorous ground test program was conducted to calibrate the instrument and ensure that the calibration can be transferred to on-orbit operations. Some calibration parameters can only be measured during pre-launch testing; one such parameter is the polarization sensitivity. Polarization testing measured the Mueller matrix components needed to determine the polarization sensitivity for all spectral bands for a series of telescope scan angles covering the expected on-orbit scan range. Results indicate that the sensitivity is below 0.6 % except at the shortest wavelengths (less than 340 nm) and was characterized to better than 0.1 % above 340 nm. This indicates that any polarized scenes measured on orbit can be corrected for with a high degree of confidence.

PACE↗

Long-Term Evaluation of Landsat 9 Oli Data Product Uniformity at Focal Plane Module-Boundaries Using Earth Scenes

The uniformity performance characteristic in remote sensing data products is a metric that helps evaluate how an image accurately represents radiometric responses to geophysical variability at the Earth’s surface or within the atmosphere. These uniformity attributes are used in Earth remote sensing science applications such as investigating water quality, land-cover/land-use classifications, and crop health assessments. To ensure these applications perform effectively, it is crucial to maintain and regularly assess the uniformity of calibrated remote sensing data products. In the Landsat missions, data from the onboard radiometric reference Solar diffuser device are used as the standard approach to derive the calibration parameters, ensuring the maintenance of both absolute radiance and radiometric uniformity through quarterly updates. This publication presents an alternative approach that leverages Earth scene data to assess the on-orbit uniformity performance of the Landsat 9 (L9) mission. The analysis characterizes uniformity at a specific set of spatial positions (the focal-plane boundary zones) that are intrinsic to the Operational Land Imager (OLI) focal plane design. By employing these Earth data statistics, the article further explores the temporal trends in the radiometric uniformity results. The results shown in this article demonstrate that for pushbroom imaging design systems with staggered sensor chip focal plane architectures, the co-registered pixels, i.e., overlapping measurements of the same ground target viewed by different detectors, offer a self-contained mechanism for assessing product uniformity. Remarkably, this alternative approach achieves a level of radiometric quality comparable to that delivered by a stable, well-characterized onboard Solar diffuser calibration device. The findings demonstrate the potential and capability of such Earth scene statistics in maintaining radiometric uniformity in calibrated data products to better than 1% (1-sigma) throughout the mission’s operational lifetime.

trend↗

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 ↗

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗

Bayesian calibration and uncertainty quantification of a rate-dependent cohesive zone model for polymer interfaces

In this work we present a rate-dependent cohesive zone model for the fracture of polymeric interfaces and performs a Bayesian calibration, an uncertainty quantification, and a sensitivity analysis for the model. The proposed cohesive zone model accounts for both reversible elastic and irreversible rate-dependent separation sliding deformation at the interface. The viscous dissipation due to the irreversible opening at the interface is modeled using elastic-viscoplastic kinematics that incorporates the effects of strain rate. Inverse calibration of parameters for such complex models through trial and error is challenging due to the large number of parameters of the model. Moreover, the calibrated parameter values are often non-unique and uncertain when the available experimental data is limited. To tackle this challenge, we employ a Bayesian calibration approach to identify parameters from experimental data, the resulting parameters significantly enhance the accuracy of the model. To quantify the uncertainty associated with the inverse parameter estimation, a modular Bayesian approach is employed to calibrate the unknown model parameters, accounting for the parameter uncertainty of the cohesive zone model. The advantages of the Bayesian calibration over a deterministic parameter fit are demonstrated. Further, to quantify the model uncertainties, such as incorrect assumptions or missing physics, a discrepancy function is introduced, which significantly improves the model’s prediction. Finally, the total uncertainty of the model is quantified in a predictive setting. A sensitivity analysis is performed to assess how changes in the input variables of the model affect the peak load, facilitating the identification of a concise set of highly influential parameters. The present approach can be used for calibration and uncertainty quantification for other complex computational mechanics models. It should also facilitate the designing of interface materials under uncertainty.

42 ENGINEERING↗

Finding Liebig’s law of the minimum

Abstract Liebig’s law of the minimum (LLM) is often used to interpret empirical biological growth data and model multiple substrates co‐limited growth. However, its mechanistic foundation is rarely discussed, even though its validity has been questioned since its introduction in the 1820s. Here we first show that LLM is a crude approximation of the law of mass action, the state of art theory of biochemical reactions, and the LLM model is less accurate than two other approximations of the law of mass action: the synthesizing unit model and the additive model. We corroborate this conclusion using empirical data sets of algae and plants grown under two co‐limiting substrates. Based on our analysis, we show that when growth is modeled directly as a function of substrate uptake, the LLM model improperly restricts the organism to be of fixed elemental stoichiometry, making it incapable of consistently resolving biological adaptation, ecological evolution, and community assembly. When growth is modeled as a function of the cellular nutrient quota, the LLM model may obtain good results at the risk of incorrect model parameters as compared to those inferred from the more accurate synthesizing unit model. However, biogeochemical models that implement these three formulations are needed to evaluate which formulation is acceptably accurate and their impacts on predicted long‐term ecosystem dynamics. In particular, studies are needed that explore the extent to which parameter calibration can rescue model performance when the mechanistic representation of a biogeochemical process is known to be deficient.

54 ENVIRONMENTAL SCIENCES↗

An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic—specifically, transversely isotropic—and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semianalytical homogenization of a bilayer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.

58 GEOSCIENCES↗

Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon

Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.

Wang, Chen↗

Hierarchical Bayesian modeling for Inverse Uncertainty Quantification of system thermal-hydraulics code using critical flow experimental data

The best estimate plus uncertainty methodology in nuclear system thermal-hydraulic studies necessitates a comprehensive understanding of uncertainties in system code predictions. The forward uncertainty quantification (UQ) process involves the propagation of input uncertainties through the computational models to obtain uncertainties in the outputs. To this end, achieving an accurate estimation of input uncertainties is important, which is the focus of inverse UQ (IUQ). Traditionally, research in Bayesian IUQ within the nuclear engineering domain has largely relied on single-level Bayesian inference. While being effective for relatively small datasets, this approach encounters limitations for cases with large datasets. The use of a single-level model may prove inefficient, as the resultant posterior distributions can significantly differ when distinct subsets of data are employed. To address this issue, we employ an hierarchical Bayesian model for IUQ. Furthermore, this approach involves organizing observations into different groups based on the test conditions, thereby accommodating varying calibration parameters across these distinct groups. In this study, we developed and implemented a hierarchical Bayesian IUQ method to consider the grouping effect of critical flow measurement data from various geometries. Comparing the outcomes of IUQ under different selections of test data using hierarchical Bayesian IUQ against those obtained from single-level Bayesian IUQ, the forward propagation of hierarchical Bayesian IUQ results demonstrates a notably improved agreement with the experimental data.

42 ENGINEERING↗

Fragmentation analysis of a bar with the Lip-field approach

The Lip-field approach was introduced in Moës and Chevaugeon (2021) as a new way to regularize softening material models. It was tested in 1D quasistatic in Moës and Chevaugeon (2021) and 2D quasistatic in Chevaugeon and Moës (2021): this paper extends it to 1D dynamics, on the challenging problem of dynamic fragmentation. The Lip-field approach formulates the mechanical problem to be solved as an optimization problem, where the incremental potential to be minimized is the non-regularized one. Spurious localization is prevented by imposing a Lipschitz constraint on the damage field. Here, the displacement and damage field at each time step are obtained by a staggered algorithm, that is the displacement field is computed for a fixed damage field, then the damage field is computed for a fixed displacement field. Indeed, these two problems are convex, which is not the case of the global problem where the displacement and damage fields are sought at the same time. The incremental potential is obtained by equivalence with a cohesive zone model, which makes material parameters calibration simple. A non-regularized local damage equivalent to a cohesive zone model is also proposed. It is used as a reference for the Lip-field approach, without the need to implement displacement jumps. These approaches are applied to the brittle fragmentation of a 1D bar with randomly perturbed material properties to accelerate spatial convergence. Both explicit and implicit dynamic implementations are compared. Favorable comparison to several analytical, numerical and experimental references serves to validate the modeling approach.

36 MATERIALS SCIENCE↗

Developing an experimental-computational framework to investigate the deformation mechanisms and mechanical properties of Al-8Ce-10Mg alloys at micro and macroscales

There is a promising future for the use of aluminum-cerium-magnesium alloys in a broad range of applications, including devices that operate at high temperatures. With cerium currently considered a waste product of rare earth mining validation studies of possible applications are essential to reduce the environmental waste. Here, a computational-experimental framework is developed to investigate the role both the intermetallic and matrix have on the mechanical properties of these alloys. A set of experiments, including SEM/EBSD imaging, nanoindentation, in-situ SEM tensile testing, and in-situ SEM-DIC tests are performed to characterize the microstructure and mechanical properties of these alloys. Furthermore, the elastic, plastic, and failure deformation mechanisms of the microstructure, and their correlation with the bulk scale mechanical properties are investigated. Experimental results are also used to calibrate parameters for a crystal plasticity finite element model, by performing a computational framework that minimizes the error between the computational and experimental results. This model is then utilized to investigate how the area percentage of intermetallics and the crystallographic texture control the mechanical properties of the alloy. Simulation results show that an increase in the percentage of intermetallics increase the strength but decrease the ductility of the alloy. Also, a change in material texture improves strength and reduces damage that leads to material failure. The development of the crystal plasticity model, as discussed in this work, opens opportunities for future investigations of similar aluminum-cerium-magnesium alloys.

36 MATERIALS SCIENCE↗

A potential model study of the nucleon's charge and mass radius

We study the charge and mass distributions within a nucleon and compute the associated squared radii based on a potential model approach. Different constituent quark configurations such as Δ, Y, and quark-diquark are considered and compared, with model parameters calibrated by experimental measurements of the proton and neutron charge radius. The results suggest that while the charge radius is dictated by quark dynamics, the mass radius is strongly influenced by nonperturbative QCD contributions to a nucleon's mass that are not sensitive to the constituent quarks. As a result, the mass radius could become substantially different from the charge radius. The obtained nucleon mass distributions of different configurations are further used for simulations of the initial conditions in heavy ion collisions. Furthermore, the computed eccentricities ε 2 and ε 3 are found to demonstrate a considerable sensitivity to the input nucleon profiles, especially to the mass radius in the peripheral region as well as for systems with fewer participants.

Gallimore, Daniel [Indiana University, Bloomington↗

Insights into Waterflooding in Hydrocarbon-Bearing Nanochannels of Varying Cross Sections from Mesoscopic Multiphase Flow Simulations

Waterflooding is one of the geotechnique used to recover fuel sources from nanoporous geological formations. The scientific understanding of the process that involves the multiphase flow of nanoconfined fluids, however, has lagged, mainly due to the complex nanopore geometries and chemical compositions. To enable benchmarked flow of nanoconfined fluids, the architected geomaterials, such as synthesized mesoporous silica with tunable pore shapes and surface chemical properties, are used for designing and conducting experiments and simulations. This work uses a modified many-body dissipative particle dynamics (mDPD) model with accurately calibrated parameters to perform parametric flow simulations for studying the influences of waterflooding driven power, pore shape, surface roughness, and surface wettability on the multiphase flow in heptane-saturated silica nanochannels. Remarkably, up to 80\% reduction in the effective permeability is found for water-driven heptane flow in a baseline 4.5 nm-wide slit channel, when compared with the Hagen–Poiseuille equation. In the 4.5 nm-wide channels with architected surface roughness, the flow rate is found either higher or lower than the baseline case, depending on the shape and size of cross-sections. High wettability of the solid surface to water is essential for achieving high recovery of heptane, regardless of surface roughness. When the solid surface is less wetting or non-wetting to water, the existence of an optimal waterflooding driven power is found to allow for the highest possible recovery. A detailed analysis on the evolution of the transient water-heptane interface in those nanochannels is presented to elucidate the underlying mechanisms that impact or dictate the multiphase flow behaviors.

58 GEOSCIENCES↗