Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inverse 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 37 records · Page 2

Quantitative analysis of surface characteristics and morphology in Death Valley, California using AIRSAR data

The Jet Propulsion Laboratory Airborne Synthetic Aperture Radar (JPL-AIRSAR) is used to collect full polarimetric measurements at P-, L-, and C-bands. These data are analyzed using the radar analysis and visualization environment (RAVEN). The AIRSAR data are calibrated using in-scene corner reflectors to allow for quantitative analysis of the radar backscatter. RAVEN is used to extract surface characteristics. Inversion models are used to calculate quantitative surface roughness values and fractal dimensions. These values are used to generate synthetic surface plots that represent the small-scale surface structure of areas in Death Valley. These procedures are applied to a playa, smooth salt-pan, and alluvial fan surfaces in Death Valley. Field measurements of surface roughness are used to verify the accuracy.

Kierein-Young, K. S.↗

Analytical modeling of structure-soil systems for lunar bases

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Macari-Pasqualino, Jose Emir↗

Photosynthetic fractionation of 13C and concentrations of dissolved CO2 in the central equatorial Pacific during the last 255,000 years

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C [is equivalent to epsilon(p)] in a central equatorial Pacific sediment core that spans the last approximately 255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic composition of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon(p), derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 [is equivalent to c(e)] based on a new, site-specific calibration of the relationship between epsilon(p) and c(e). The calibration was based on reassessment of existing epsilon(p) versus c(e) data, which support a physiologically based model in which epsilon(p) is inversely related to c(e). Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index U(K/37). Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon(p) and 1/c(e). These are discussed in detail and it is concluded that the observed record of epsilon(p) most probably reflects significant variations in delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from approximately 110 microatmospheres during glacial intervals (ocean > atmosphere) to approximately 60 microatmospheres during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p<0.01) inverse correlations of delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

NASA Discipline Exobiology↗

Orion EFT-1 Cavity Heating Tile Experiments and Environment Reconstruction

Developing aerothermodynamic environments for deep cavities, such as those produced by micrometeoroids and orbital debris impacts, poses a great challenge for engineers. In order to assess existing cavity heating models, two one-inch diameter cavities were flown on the Orion Multi-Purpose Crew Vehicle during Exploration Flight Test 1 (EFT1). These cavities were manufactured with depths of 1.0 in and 1.4 in, and they were both instrumented. Instrumentation included surface thermocouples upstream, downstream and within the cavities, and additional thermocouples at the TPS-structure interface. This paper will present the data obtained, and comparisons with computational predictions will be shown. Additionally, the development of a 3D material thermal model will be described, which will be used to account for the three-dimensionality of the problem when interpreting the data. Furthermore, using a multi-dimensional inverse heat conduction approach, a reconstruction of a time- and space-dependent flight heating distribution during EFT1 will be presented. Additional discussions will focus on instrumentation challenges and calibration techniques specific to these experiments. The analysis shown will highlight the accuracies and/or deficiencies of current computational techniques to model cavity flows during hypersonic re-entry.

Salazar, Giovanni↗

Stochastic Reconstruction of Thermal Protection Material Properties from Arc-Jet Experiments

Material response models are used to assess reliability using variances in the bond-line temperature predictions based on uncertainties in trajectory, aerothermal environment, and material properties. A key deficiency in the current approach is that input uncertainties are too often subjective, empirical, or ad-hoc, and are not rigorously linked to the arc-jet test data used to develop the TPS material model. While materials such as PICA are well understood, future missions may require more novel materials such as HEEET where unknown uncertainties have real consequences on the ability to assess reliability. A quantifiable estimate of reliability requires an iterative methodology where the parameters driving the variance in bond-line temperature (for example) are systematically identified. A test campaign to collect data or develop new models can then be identified to reduce those input uncertainties. A Bayesian inference loop defines these connections mathematically, i.e., prior knowledge about uncertainty is updated based on observation. While these concepts are well known (and often applied intuitively in a non-rigorous approach), only recent advances in reduced-order modelling have made them computationally viable methods for engineering. By replacing deterministic inverse methods with stochastic approaches, the hope is new materials proposed for future missions can more rapidly be developed with a greater understanding of the TPS material reliability. Two additional steps for the analysis of arc jet test data are discussed. The first is ability to construct a reduced-order model using material response simulations (Icarus/US3D) of the arc-jet test articles, and the second is the inclusion of this surrogate model in the Bayesian inversion process. Both capabilities will be demonstrated using prior PICA arc-jet test data. The quality of a surrogate model will be investigated and the variances on the calibrated material properties will be compared to our current understanding of the PICA material model.

Material response↗

MODIS Solar Diffuser Degradation Determination and Its Spectral Dependency

This study presents a modeling approach to improve solar diffuser (SD) degradation determination from SD stability monitor (SDSM) measurements. The MODIS instrument uses a SDto calibrate its reflective solar bands (RSBs) on-orbit. Due to the imperfectly designed SDSM sun view screen, the SD reflectance tracked by SDSM has large noise. The SDSM measurements noise is spectrally coherent and can be minimized by normalizing measurements to the least degradeddetector 9 (936 nm). In this study, a SD degradation model is used to determine the SDdegradation's wavelength dependency and the detector 9 degradation is estimated by the model solution.The results show the SD degradations measured at 6 SDSM detectors (554 _ 936 nm) have stable relationships, where the degradation is inversely proportion to 1/wavelength^4. The model estimated SD degradation at SDSM detector 9 wavelength (936 nm) is ~0.9% from 2002 to 2018.Based on the SD degradation model solution, the SD degradation at short/mid wave bands are estimated to improve short/mid wave bands calibration. The model can also be used to improve interpolating SD degradation at SDSM detectors to RSB wavelengths. Compared to linear interpolation, bands 9 and 10 show the largest differences of up to 0.3 and 0.4% respectively. These differences directly impact the calibration coefficients of these bands.

radiometric calibration↗

Predicting cell death and mutation frequency for a wide spectrum of LET and across cell lines from DNA break clustering inside repair domains

Cosmic radiation, which is composed of high charged and energy (HZE) particles, is responsible for cell death and mutation, which may be involved in cancer induction. Mutations are consequences of mis-repaired DNA breaks – especially double-strand breaks (DSBs) – that induce inter- and intra-chromosomal rearrangements (translocations, deletions, inversion). In this study, a computer simulation model is used to investigate the clustering of DSBs in repair domains, previously evidenced by our group in human breast cells. This model is calibrated with experimental data measuring persistent 53BP1 radiation-induced foci (RIF) and is used to explain the high relative biological effectiveness (RBE) of HZE for both cell death and DNA mutation frequencies. We first validate our DSB cluster model using a new track structure model deployed on a simple geometrical configuration for repair domains in the nucleus; then we extend the scope from cell death to mutation induction. This work suggests that mechanism based on DSB repair process can explain several biological effects induced by HZE particles on different type of living cells.

Ianik Plante↗

Bayesian Inference and the Effects of Varying Uncertainty Models in Charring Ablator Calibration and Uncertainty Quantification Problems

The Mars Science Laboratory (MSL) vehicle utilized a heat shield constructed from NASA’s Phenolic-Impregnated Carbon Ablator (PICA) material to protect the main structure from the high enthalpy environment encountered during hypersonic atmospheric entry. During the vehicle’s descent through Martian atmosphere, multiple thermocouples embedded within the heat shield captured in-depth material temperature data that allow for studies to be conducted on current material response reconstruction tools. In the present work, material temperature data obtained from thermocouples within the MISP-4 plug (MEDLI (Mars Science Laboratory Entry, Descent, and Landing Instrument) Integrated Sensor Plug) are utilized in the calibration of Theoretical Ablative Composite for Open Testing (TACOT) model parameters in conjunction with NASA’s Porous material Analysis Toolbox (PATO) through Bayesian inference where uncertainty due to parametric, modeling, and experimental sources is simultaneously quantified. Prior to the study, a sensitivity analysis is performed through computation of the robust Sobol indices in an effort to study the relationship between input space and model response and to reduce the dimensionality of the statistical inverse problem. The Bayesian inference methodology necessitates an a-priori choice to be made for the uncertainty model for which numerous possibilities are available. Across most works, however, only basic additive or multiplicative models are utilized with pre-defined magnitudes of uncertainty based on a-priori knowledge or to-be-calibrated multipliers of static covariance matrix structures. The present effort explores the effects of informed uncertainty models, ones with temporal dependence that are simultaneously calibrated through Bayesian inference, on calibrated results for parameters that make up the uncertain input space.

Sensitivity Analysis↗

Predicting Cell Death and Mutation Frequency for a Wide Spectrum of LET by Assuming DNA Break Clustering Inside Repair Domains

Cosmic radiation, which is composed of high charged and energy (HZE) particles, is responsible for cell death and mutation, which may be involved in cancer induction. Mutations are consequences of mis-repaired DNA breaks – especially double-strand breaks (DSBs) – that induce inter- and intra-chromosomal rearrangements (translocations, deletions, inversion). In this study, a computer simulation model is used to investigate the clustering of DSBs in repair domains, previously evidenced by our group in human breast cells [1]. This model is calibrated with experimental data measuring persistent 53BP1 radiation-induced foci (RIF) and is used to explain the high relative biological effectiveness (RBE) of HZE for both cell death and DNA mutation frequencies. We first validate our DSB cluster model using a new track structure model deployed on a simple geometrical configuration for repair domains in the nucleus; then we extend the scope from cell death to mutation induction. This work suggests that mechanism based on DSB repair process can explain several biological effects induced by HZE particles on different type of living cells

Plante, Ianik↗

Validity of the Temperature Reconstruction from Water Isotopes in Ice Cores

Well-documented present-day distributions of stable water isotopes (HDO and others) show the existence, in middle and high latitudes, of a linear relationship between the mean annual isotope content of precipitation (SD and 51"0) and the mean annual temperature at the precipitation site. Paleoclimatologists have used this relationship, which is particularly well obeyed over Greenland and Antarctica, to infer paleotemperatures from ice core data. There is, however, growing evidence that spatial and temporal isotope/ surface temperature slopes differ, thus complicating the use of stable water isotopes as paleothermometers. In this paper we review empirical estimates of temporal slopes in polar regions and relevant information that can be inferred from isotope models: simple, Rayleigh-type distillation models and (particularly over Greenland) general circulation models (GCMS) fitted with isotope tracer diagnostics. Empirical estimates of temporal slopes appear consistently lower than present-day spatial slopes and are dependent on the timescale considered. This difference is most probably due to changes in the evaporative origins of moisture, changes in the seasonality of the precipitation, changes in the strength of the inversion layer, or some combination of these changes. Isotope models have not yet been used to evaluate the relative influences of these different factors. The apparent disagreement in the temporal and spatial slopes clearly makes calibrating the isotope paleothermometer difficult. Nevertheless, the use of a (calibrated) isotope paleothermometer appears justified; empirical estimates and most (though not all) GCM results support the practice of interpreting ice core isotope records in terms of local temperature changes.

Jouzel, J.↗

Fine Particulate Matter Predictions Using High Resolution Aerosol Optical Depth (AOD) Retrievals

To date, spatial-temporal patterns of particulate matter (PM) within urban areas have primarily been examined using models. On the other hand, satellites extend spatial coverage but their spatial resolution is too coarse. In order to address this issue, here we report on spatial variability in PM levels derived from high 1 km resolution AOD product of Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm developed for MODIS satellite. We apply day-specific calibrations of AOD data to predict PM(sub 2.5) concentrations within the New England area of the United States. To improve the accuracy of our model, land use and meteorological variables were incorporated. We used inverse probability weighting (IPW) to account for nonrandom missingness of AOD and nested regions within days to capture spatial variation. With this approach we can control for the inherent day-to-day variability in the AOD-PM(sub 2.5) relationship, which depends on time-varying parameters such as particle optical properties, vertical and diurnal concentration profiles and ground surface reflectance among others. Out-of-sample "ten-fold" cross-validation was used to quantify the accuracy of model predictions. Our results show that the model-predicted PM(sub 2.5) mass concentrations are highly correlated with the actual observations, with out-of- sample R(sub 2) of 0.89. Furthermore, our study shows that the model captures the pollution levels along highways and many urban locations thereby extending our ability to investigate the spatial patterns of urban air quality, such as examining exposures in areas with high traffic. Our results also show high accuracy within the cities of Boston and New Haven thereby indicating that MAIAC data can be used to examine intra-urban exposure contrasts in PM(sub 2.5) levels.

aerosol optical depth↗

Helioseismic Measurements of Convective Power in Solar Cycle 24

Constraining the parameters under which convection in the solar interior operates has important implications for describing how energy is transported by plasma motions, and various physical models have been employed to provide some theoretical estimates on the expected power spectrum. Past attempts to measure the convective power distributed among large spatial scales have, however, found differing and incompatible values. Here, we present measurements of the convective power spectrum in the upper convection zone for Carrington rotations in Solar Cycle 24 obtained from the helioseismic signal corresponding to East-West flows. We also perform calibration on synthetic data using the global acoustic GALE code to make assessments of the flow velocities at various length scales without the need for performing inversions. This allows us to derive the convective power spectrum from the flow maps and to compare with the results from the raw travel times. These results are compared against predictions of convective power in global models of convection produced by the EULAG code. Finally, we show how the steps in our analysis procedure (for example data segmentation, filtering, etc.) affect our estimates of the convective power by comparing with the synthetic data from the GALE code.

Heliophysics↗

Middle infrared multispectral aircraft scanner data - Analysis for geological applications

Multispectral middle IR (8-13 microns) data were acquired with an aircraft scanner over Utah. Because these digital image data were dominated by temperature, all six channels were highly correlated. Extensive processing was required to allow geologic photointerpretation based on subtle variations in spectral emittance between rock types. After preliminary processing, ratio images were produced and color ratio composites created from these. Sensor calibration and an atmospheric model allowed determination of surface brightness, temperature, emittance, and color composite emittance images. The best separation of major rock types was achieved with a principal component transformation, followed by a Gaussian stretch, followed by an inverse transformation to the original axes.

Kahle, A. B.↗

The spectral energy distribution of Zeta Puppis and HD 50896

The ultraviolet spectral energy distribution of the O5f star Zeta Pup and the WN5 star HD 50896 are derived from OAO-2 observations with the calibration of Bless, Code, and Fairchild (1976). An estimate of the interstellar reddening (0.12 magnitude) of the Wolf-Rayet star is determined from the size of the characteristic interstellar extinction bump at 4.6 inverse microns. After correction for extinction, both stars show a flat energy distribution in the ultraviolet. The distribution of HD 50896 from 1100 A to 2 microns is in good agreement with results of extended model atmospheres, but some uncertainty remains because of the interstellar-extinction correction. The absolute energy distribution of Zeta Pup is fitted by a 42,000-K plane-parallel model if the model's flux is adjusted for the effects of electron scattering in the stellar wind and for UV line blanketing that was determined empirically from high-resolution Copernicus satellite observations. To achieve this fit, it is necessary to push both the spectroscopically determined temperature and the ultraviolet calibration to the limits of their probable errors.

Holm, A. V.↗

Analysis and calibration of a three-invariant plasticity model for granular materials

This paper briefly reviews issues related to the characterization of properties of granular materials subjected to microgravity and one-gravity conditions at very low effective stress levels. The development of a three-invariant plasticity model that resembles the model devised by Lade (1979) is described. An inverse-identification scheme where the analysis tools are used to extract constitutive model parameters from experiments is also discussed.

Sture, S.↗

Thermal and Pressure Characterization of a Wind Tunnel Force Balance Using the Single Vector System

Wind tunnel research at NASA Langley Research Center s 31-inch Mach 10 hypersonic facility utilized a 5-component force balance, which provided a pressurized flow-thru capability to the test article. The goal of the research was to determine the interaction effects between the free-stream flow and the exit flow from the reaction control system on the Mars Science Laboratory aeroshell during planetary entry. In the wind tunnel, the balance was exposed to aerodynamic forces and moments, steady-state and transient thermal gradients, and various internal balance cavity pressures. Historically, these effects on force measurement accuracy have not been fully characterized due to limitations in the calibration apparatus. A statistically designed experiment was developed to adequately characterize the behavior of the balance over the expected wind tunnel operating ranges (forces/moments, temperatures, and pressures). The experimental design was based on a Taylor-series expansion in the seven factors for the mathematical models. Model inversion was required to calculate the aerodynamic forces and moments as a function of the strain-gage readings. Details regarding transducer on-board compensation techniques, experimental design development, mathematical modeling, and wind tunnel data reduction are included in this paper.

Lynn, Keith C.↗

Mapping From an Instrumented Glove to a Robot Hand

An algorithm has been developed to solve the problem of mapping from (1) a glove instrumented with joint-angle sensors to (2) an anthropomorphic robot hand. Such a mapping is needed to generate control signals to make the robot hand mimic the configuration of the hand of a human attempting to control the robot. The mapping problem is complicated by uncertainties in sensor locations caused by variations in sizes and shapes of hands and variations in the fit of the glove. The present mapping algorithm is robust in the face of these uncertainties, largely because it includes a calibration sub-algorithm that inherently adapts the mapping to the specific hand and glove, without need for measuring the hand and without regard for goodness of fit. The algorithm utilizes a forward-kinematics model of the glove derived from documentation provided by the manufacturer of the glove. In this case, forward-kinematics model signifies a mathematical model of the glove fingertip positions as functions of the sensor readings. More specifically, given the sensor readings, the forward-kinematics model calculates the glove fingertip positions in a Cartesian reference frame nominally attached to the palm. The algorithm also utilizes an inverse-kinematics model of the robot hand. In this case, inverse-kinematics model signifies a mathematical model of the robot finger-joint angles as functions of the robot fingertip positions. Again, more specifically, the inverse-kinematics model calculates the finger-joint commands needed to place the fingertips at specified positions in a Cartesian reference frame that is attached to the palm of the robot hand and that nominally corresponds to the Cartesian reference frame attached to the palm of the glove. Initially, because of the aforementioned uncertainties, the glove fingertip positions calculated by the forwardkinematics model in the glove Cartesian reference frame cannot be expected to match the robot fingertip positions in the robot-hand Cartesian reference frame. A calibration must be performed to make the glove and robot-hand fingertip positions correspond more precisely. The calibration procedure involves a few simple hand poses designed to provide well-defined fingertip positions. One of the poses is a fist. In each of the other poses, a finger touches the thumb. The calibration subalgorithm uses the sensor readings from these poses to modify the kinematical models to make the two sets of fingertip positions agree more closely.

Goza, Michael↗

Rover mast calibration, exact camera pointing, and camara handoff for visual target tracking

This paper presents three technical elements that we have developed to improve the accuracy of the visual target tracking for single-sol approach-and-instrument placement in future Mars rover missions. An accurate, straightforward method of rover mast calibration is achieved by using a total station, a camera calibration target, and four prism targets mounted on the rover. The method was applied to Rocky8 rover mast calibration and yielded a 1.1-pixel rms residual error. Camera pointing requires inverse kinematic solutions for mast pan and tilt angles such that the target image appears right at the center of the camera image. Two issues were raised. Mast camera frames are in general not parallel to the masthead base frame. Further, the optical axis of the camera model in general does not pass through the center of the image. Despite these issues, we managed to derive non-iterative closed-form exact solutions, which were verified with Matlab routines. Actual camera pointing experiments aver 50 random target image paints yielded less than 1.3-pixel rms pointing error. Finally, a purely geometric method for camera handoff using stereo views of the target has been developed. Experimental test runs show less than 2.5 pixels error on high-resolution Navcam for Pancam-to-Navcam handoff, and less than 4 pixels error on lower-resolution Hazcam for Navcam-to-Hazcam handoff.

camera handoff↗