Search NASA⌕ Search

SEARCH · Search NASA

Results for “value of data”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Proxy-based Bayesian inversion of strain tensor data measured during well tests

Recent instrument developments have made it possible to measure the strain tensor caused by injecting or pumping fluid from aquifers or reservoirs, but the full value of these data is limited because the long runtimes of poroelastic forward models makes it impractical to use many inversion schemes. This limits the interpretation of strain data for managing the recovery of resources or storage of wastes in the subsurface. This paper describes a method of inverting deformation data using a poroelastic numerical simulator so the results can be used to manage reservoirs or aquifers. We developed a workflow designed to reduce the number of simulations sufficiently to make it feasible to use DREAMzs, an advanced Bayesian inversion method that translates the uncertainties from different sources into unbiased posterior parameter distributions and uncertainty envelopes around the field data. Using a KNN proxy model for the poroelastic simulator is key to reducing the overall computations, and the workflow includes a strategy for ensuring the proxy model results converge on the results from the simulator. The workflow is tested using an idealized example that verifies the ability to correctly identify parameters and characterize noise used to perturb the data. Field data from an injection test at an oil reservoir near Tulsa, Oklahoma, are also used to evaluate the efficacy of the workflow with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), where the results were used to characterize posterior parameter distribution and evaluate the prediction uncertainty. Furthermore, this workflow is significant because it enables strain tensor, or other geomechanical measurements to be interpreted to guide decision-making during energy and environmental processes in the subsurface.

42 ENGINEERING↗
DOE OSTI2025

WHONDRS laboratory time series moisture manipulative experiment from soil core layers across eastern contiguous US: time series aerobic respiration, geochemistry, and aggregates

This dataset supports a broader study examining the effects of wetting and drying on soil layers across the eastern contiguous United States (CONUS). The dataset provides data generated from a laboratory moisture manipulation experiment. The contents include time series aerobic respiration and moisture; dissolved oxygen; sediment geochemistry data; and field metadata. Samples were collected as part of a collaboration between WHONDRS (Worldwide Hydrobiogeochemistry Observation Network for Dynamic River Systems; https://whondrs.pnnl.gov) and MONet (Molecular Observation Network; https://www.emsl.pnnl.gov/monet). The field samples (soil cores) were labeled as MEL_##_COR and subsequent subsamples begin with MEL_##. Additional subsamples were taken for the laboratory experiment and were labeled as EL_##. The labels from the MEL field samples and the EL subsamples can be mapped directly based on the digits following the prefix and underscore (i.e., EL_01 is a subsample from MEL_01). See the critical details section below for more details on sample naming and experimental design.For details on how to navigate this data package, see this infographic from the River Corridor SFA https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About.In addition to this readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions.This dataset is comprised of (1) a folder containing environmental context photos; (2) file-level metadata; (3) data dictionary; (4) field metadata; (5) readme; (6) international generic sample number (IGSN) mapping file; and (7) a subfolder with soil sample data from field samples and the incubation experiment. The sample data subfolder contains (1) effect size; (2) gravimetric moisture from field samples and incubation experiment; (3) respiration rates, raw dissolved oxygen values, and plots; (4) specific conductance, pH, and temperature from the incubation; (5) soil aggregates; (6) a summary containing median values of each data type for each treatment (wet and dry) in the incubation; (7) a summary containing averages for each data type of each soil layer; and (8) methods codes. All files are .csv, .pdf, .jpeg, or .jpg.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

COMPASS-FME Terrestrial Ecosystem Manipulation to Probe the Effects of Storm Treatments (TEMPEST) Soil Greenhouse Gas Fluxes 04/2020-06/2024

Both raw and processed data measured using a LI-COR 7810 Greenhouse Gas Analyzer in the TEMPEST experiment, part of COMPASS-FME (https://compass.pnnl.gov/FME/COMPASSFME), at Smithsonian Environmental Research Center. This ecosystem-scale experiment probes the effects of saltwater versus freshwater flooding in a coastal deciduous forest.The data consist of both concentrations and fluxes of carbon dioxide and methane measured using static chambers on the soil surface, approximately every two weeks from early 2020 to mid 2024. Some of the measurement points are controls, and some subject to root-exclusion techniques; all measurements are embedded in the TEMPEST control, freshwater, and saltwater plots (see Hopple et al. 2023). These data were generated to understand changing soil greenhouse gas (CO2 and CH4) production and consumption. File types are comma-separated value (.csv) for data, and markdown (.md) for supplementary information.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Elimination LArTPC Simulation Uncertainty

Liquid Argon Time Projection Chambers (LArTPC) are essential for detecting muons and neutrinos by capturing electrons released during particle collisions, which drift toward wire planes under an electric field and induce currents measured to reconstruct particle paths. However, LArTPCs face challenges from effects such as electron-ion recombination, electron diffusion, and electron attenuation, complicating data simulation. The Short Baseline Neutrino (SNB) detector aims to measure neutrinos before oscillation occurs. To bridge the gap between simulation and actual data, we propose modifying the amplitude and width of signals on the TPC wires, addressing uncertainties by adjusting signal characteristics to better match observed data. A Gaussian fit to current waveforms produces hits with associated charge and width, and by comparing data and simulated values, discrepancies highlight areas where the model fails. Initial results indicate the current modification algorithm may increase divergence between simulation and data, necessitating further refinement. A discovered bug in the WireModMakeHist_plug.cpp file, which incorrectly computed simulation and data ratios, underscores the need for precise algorithm adjustments. Future work involves correcting code errors, fine-tuning the model, and conducting multiple simulation runs to enhance statistical confidence and reduce uncertainties, ultimately aiming for accurate LArTPC operation and reliable neutrino detection.

Mkrtchyan, Ka'ren↗
NASA NTRS1995

The Atmospheric Data Acquisition And Interpolation Process For Center-TRACON Automation System

The Center-TRACON Automation System (CTAS), an advanced new air traffic automation program, requires knowledge of spatial and temporal atmospheric conditions such as the wind speed and direction, the temperature and the pressure in order to accurately predict aircraft trajectories. Real-time atmospheric data is available in a grid format so that CTAS must interpolate between the grid points to estimate the atmospheric parameter values. The atmospheric data grid is generally not in the same coordinate system as that used by CTAS so that coordinate conversions are required. Both the interpolation and coordinate conversion processes can introduce errors into the atmospheric data and reduce interpolation accuracy. More accurate algorithms may be computationally expensive or may require a prohibitively large amount of data storage capacity so that trade-offs must be made between accuracy and the available computational and data storage resources. The atmospheric data acquisition and processing employed by CTAS will be outlined in this report. The effects of atmospheric data processing on CTAS trajectory prediction will also be analyzed, and several examples of the trajectory prediction process will be given.

Jardin, M. R.↗
NASA NTRS2022

NASA’s Global Precipitation Measurement Mission: Leveraging Stakeholder Engagement & Applications Activities to Inform Decision-making

The application of satellite precipitation estimates from NASA’s Global Precipitation Measurement (GPM) Mission for decision-making has been a focus for the mission since launch. As a result, GPM data have enabled a range of applications that address societal needs, including water resource management, crop forecasting, ecological monitoring, disaster response, public health, aviation, weather forecasting, and climate modeling, among others. GPM applications activities have continued to focus on user engagement through in person trainings and interviews, workshops, webinars, and educational outreach activities. The goals of these efforts are to synthesize community data needs in order to effectively support and enable decision-making across agencies, academia and the global community. While these efforts have helped the GPM mission establish a large stakeholder community that encompasses federal and state partners, academic institutions, nd private and nonprofit companies, there remains difficulties associated with accessing, processing, and applying the data to support or enable applications. In this article, we present GPM applications strategies and approaches used to enhance the applications value of GPM data, and most importantly, demonstrate how these efforts have and can inform different decision-making contexts. This work also provides a discussion on key lessons learned from the user community and how this information can be utilized to help better support and shape applications approaches for future NASA Earth Science missions.

Satellite precipitation↗
DOE OSTI2015

MO-F-CAMPUS-T-03: Data Driven Approaches for Determination of Treatment Table Tolerance Values for Record and Verification Systems

Purpose: To determine appropriate couch tolerance values for modern radiotherapy linac R&V systems with indexed patient setup. Methods: Treatment table tolerance values have been the most difficult to lower, due to many factors including variations in patient positioning and differences in table tops between machines. We recently installed nine linacs with similar tables and started indexing every patient in our clinic. In this study we queried our R&V database and analyzed the deviation of couch position values from the acquired values at verification simulation for all patients treated with indexed positioning. Mean and standard deviations of daily setup deviations were computed in the longitudinal, lateral and vertical direction for 343 patient plans. The mean, median and standard error of the standard deviations across the whole patient population and for some disease sites were computed to determine tolerance values. Results: The plot of our couch deviation values showed a gaussian distribution, with some small deviations, corresponding to setup uncertainties on non-imaging days, and SRS/SRT/SBRT patients, as well as some large deviations which were spot checked and found to be corresponding to indexing errors that were overriden. Setting our tolerance values based on the median + 1 standard error resulted in tolerance values of 1cm lateral and longitudinal, and 0.5 cm vertical for all non- SRS/SRT/SBRT cases. Re-analizing the data, we found that about 92% of the treated fractions would be within these tolerance values (ignoring the mis-indexed patients). We also analyzed data for disease site based subpopulations and found no difference in the tolerance values that needed to be used. Conclusion: With the use of automation, auto-setup and other workflow efficiency tools being introduced into radiotherapy workflow, it is very essential to set table tolerances that allow safe treatments, but flag setup errors that need to be reassessed before treatments.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2003

Interleaved Observation Execution and Rescheduling on Earth Observing Systems

Observation scheduling for Earth orbiting satellites solves the following problem: given a set of requests for images of the Earth, a set of instruments for acquiring those images distributed on a collecting of orbiting satellites, and a set of temporal and resource constraints, generate a set of assignments of instruments and viewing times to those requests that satisfy those constraints. Observation scheduling is often construed as a constrained optimization problem with the objective of maximizing the overall utility of the science data acquired. The utility of an image is typically based on the intrinsic importance of acquiring it (for example, its importance in meeting a mission or science campaign objective) as well as the expected value of the data given current viewing conditions (for example, if the image is occluded by clouds, its value is usually diminished). Currently, science observation scheduling for Earth Observing Systems is done on the ground, for periods covering a day or more. Schedules are uplinked to the satellites and are executed rigorously. An alternative to this scenario is to do some of the decision-making about what images are to be acquired on-board. The principal argument for this capability is that the desirability of making an observation can change dynamically, because of changes in meteorological conditions (e.g. cloud cover), unforeseen events such as fires, floods, or volcanic eruptions, or un-expected changes in satellite or ground station capability. Furthermore, since satellites can only communicate with the ground between 5% to 10% of the time, it may be infeasible to make the desired changes to the schedule on the ground, and uplink the revisions in time for the on-board system to execute them. Examples of scenarios that motivate an on-board capability for revising schedules include the following. First, if a desired visual scene is completely obscured by clouds, then there is little point in taking it. In this case, satellite resources, such as power and storage space can be better utilized taking another image that is higher quality. Second, if an unexpected but important event occurs (such as a fire, flood, or volcanic eruption), there may be good reason to take images of it, instead of expending satellite resources on some of the lower priority scheduled observations. Finally, if there is unexpected loss of capability, it may be impossible to carry out the schedule of planned observations. For example, if a ground station goes down temporarily, a satellite may not be able to free up enough storage space to continue with the remaining schedule of observations. This paper describes an approach for interleaving execution of observation schedules with dynamic schedule revision based on changes to the expected utility of the acquired images. We describe the problem in detail, formulate an algorithm for interleaving schedule revision and execution, and discuss refinements to the algorithm based on the need for search efficiency. We summarize with a brief discussion of the tests performed on the system.

Khatib, Lina↗
DOE OSTI2021

Developing an Automated Uncertainty Quantification Tool to Improve Watershed-Scale Predictions of Water and Nutrient Cycling

Managing the flow of water, nutrients, and contaminants in watersheds is vital to addressing pressing issues related to water scarcity, access to clean drinking water, energy production, resilience to natural and anthropogenic perturbations, and ecological restoration. Decisions about the management of watersheds critically depend on the accuracy with which the flow of water and chemicals through the watershed can be predicted by computer models. Prediction uncertainty can be reduced by matching the model to data, which are collected in the field at great expense. The contribution of watershed characterization data to reducing uncertainty of relevant model predictions can be evaluated in a so-called data-worth analysis, which provides transparent, quantitative metrics about a data set’s value for the support of relevant watershed management objectives. To achieve this goal, we developed a software package that implements the data-worth analysis approach for use with state-of-the-art watershed models. The purpose of the proposed data-worth analysis is to help decision-makers allocate resources for watershed characterization such that the uncertainty in model predictions can be significantly reduced, which leads to better, more effective management decisions. At the same time, watershed characterization costs can be reduced. The specific technical objectives of this SBIR/STTR Phase II project were to develop a framework and associated software toolsets that implement the uncertainty quantification and data-worth analysis approach for use with state-of-the-art watershed models. This goal was achieved by (A) developing a user-friendly, robust software package that is accessible to a wide audience, including watershed managers, policy-makers, and public stakeholders; (B) by demonstrating application of the prototype on several use cases that are representative of complex watershed management challenges spanning a range of scales and that consider different open-source, DOE-based codes and other modeling platforms; and (C) by gathering information about the needs and requirements from potential users to help guide future developments, ensuring that the final product will be commercially viable. The developed software consists of a graphical user interface that guides the user through a sequence of analysis steps, supported by toolsets that leverage state-of-the-art computational simulation-optimization capabilities. A prototype of the software runs on multiple platforms (PC, Mac, multi-processor Linux environment), is linked to diverse watershed simulators (e.g., ECOSYS, TOUGH2, TOUGHREACT, Amanzi-ATS), performs multiple analysis tasks (predictive simulations, sensitivity analysis, uncertainty analysis, automatic parameter estimation, and data-worth analysis, multicomponent geothermometry), and is readily extensible to include external simulators and analysis tools. The software is being commercialized and will be continually updated to address user needs.

58 GEOSCIENCES↗
NASA NTRS1989

Testing general relativity with Landers on the Martian satellite Phobos

A planned experiment to obtain range and Doppler data with the Phobos 2 Lander on the surface of the Martian satellite Phobos is described. With the successful insertion on January 29, 1989 of Phobos 2 into Mars orbit, it is anticipated that the Lander will be placed on the surface of Phobos in April 1989. Depending on the longevity of the Lander, range and Doppler data for a period of from one to several years are expected. Because these data are of value in performing solar-system tests of general relativity, the current accuracy of the relevant relativity tests using Deep Space Network data from the Mariner-9 orbiter of Mars in 1971 and from the Viking Landers in 1976-1982 is reviewed. The expected improvement from data anticipated during the Phobos 2 Lander Mission is also discussed; most important will be an improved sensitivity to any time variation in the gravitational 'constant' as measured in atomic units.

Anderson, J. D.↗
NASA NTRS1993

Global Reference Atmospheric Model - 1990

GRAM-90, improvement on original model, provides values of atmospheric parameters, either automatically at positions along linear path or along any set of connected positions specified by user. Parameters produced include monthly total and mean values of pressure, density, temperature, and wind, as well as perturbation values and magnitudes. Also provides thermal wind shear for monthly mean winds, percent deviation from standard atmosphere, mean vertical wind, and perturbation data for each position. Perturbation data consist of values of quasi-biennial and random perturbations. Written in FORTRAN 77.

Justus, C. G.↗
NASA NTRS2017

Oculomotor Assessment of Diurnal Arousal Variations

Saccadic and pupillary responses are reliable indices of arousal decrement (e.g. fatigue), that might be exploited to improve work schedule guidelines. In this study, we tested the sensitivity of a short 30-s oculomotor test to detect diurnal arousal variations. Twelve participants (5 females, 7 males, 37.7+-11.9 years) volunteered to be assessed every hour (66+-20 min) for three consecutive working days, during their regular office-hours. We used a fully automated testing system, the FIT 2000 Fitness Impairment Tester (Pulse Medical Instruments Inc., Rockville, MD, USA), to measure and record saccadic peak velocity, pupil diameter, and latency and amplitude of the pupillary light reflex. In addition, we collected subjective levels of arousal using the Stanford Sleepiness Scale, and body core temperature. We analyzed the data using a linear mixed model approach for longitudinal data. Both saccadic velocity and subjective alertness decreased over the course of a day, while body core temperature increased (all p-values.05). The data also weakly suggested an increase of the pupil diameter (p 07). The findings support the use of oculomotor indices in the assessment of arousal and fatigue in applied settings.

saccadic velocity↗
DOE OSTI2025

Elucidation of Marcus Relationships for Hydride Transfer Reactions Involving Transition Metal Hydrides

The rate of hydride transfer from three Ir hydride complexes of the type Cp*Ir( R bpy)H + (Cp* = C 5 Me 5 ; R bpy = 4,4′-R-2,2′-bipyridine, R = OMe, H, CO 2 Me) to six N-methylacridinium ( R Acr + ) acceptors with electronically different substituents in the 2- or 2,7-positions were measured. Using the thermodynamic hydricity of the donors and the hydride affinity of the acceptors the thermodynamic driving forces for hydride transfer were determined. Brønsted plots, which correlate kinetic and thermodynamic hydricity, demonstrate distinct linear free energy relationships for each complex, with different Brønsted α values. Thus, at the same driving force hydride transfer from Cp*Ir( OMe bpy)H + is faster than for Cp*Ir(bpy)H + or Cp*Ir( CO2Me bpy)H + . Experimental and computational analyses are consistent with a concerted hydride transfer mechanism for all Ir complexes. As the thermodynamic driving force increases an earlier transition state is observed and all transition states also include π-stacking interactions between the donor and acceptor, which likely contribute to the different α values. The experimental data fits well to the Marcus model, enabling the determination of reorganization energies (λ) that range from 58 to 69 kcal mol -1 . These are lower than λ values for hydride transfer reactions involving organic donors and acceptors. This work provides a rare example of the correlation of kinetic and thermodynamic hydricity using only experimental data and shows that hydride transfer reactions involving metal hydrides can follow Marcus theory. Furthermore, the findings offer insight into controlling metal-catalyzed hydride transfer reactions, which is valuable for designing improved systems for a range of transformations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

XCOM intrinsic dimensionality for low-Z elements at diagnostic energies

Purpose: To determine the intrinsic dimensionality of linear attenuation coefficients (LACs) from XCOM for elements with low atomic number (Z = 1-20) at diagnostic x-ray energies (25-120 keV). H{sub 0}{sup q}, the hypothesis that the space of LACs is spanned by q bases, is tested for various q-values. Methods: Principal component analysis is first applied and the LACs are projected onto the first q principal component bases. The residuals of the model values vs XCOM data are determined for all energies and atomic numbers. Heteroscedasticity invalidates the prerequisite of i.i.d. errors necessary for bootstrapping residuals. Instead wild bootstrap is applied, which, by not mixing residuals, allows the effect of the non-i.i.d residuals to be reflected in the result. Credible regions for the eigenvalues of the correlation matrix for the bootstrapped LAC data are determined. If subsequent credible regions for the eigenvalues overlap, the corresponding principal component is not considered to represent true data structure but noise. If this happens for eigenvalues l and l + 1, for any l{<=}q, H{sub 0}{sup q} is rejected. Results: The largest value of q for which H{sub 0}{sup q} is nonrejectable at the 5%-level is q = 4. This indicates that the statistically significant intrinsic dimensionality of low-Z XCOM data at diagnostic energies is four. Conclusions: The method presented allows determination of the statistically significant dimensionality of any noisy linear subspace. Knowledge of such significant dimensionality is of interest for any method making assumptions on intrinsic dimensionality and evaluating results on noisy reference data. For LACs, knowledge of the low-Z dimensionality might be relevant when parametrization schemes are tuned to XCOM data. For x-ray imaging techniques based on the basis decomposition method (Alvarez and Macovski, Phys. Med. Biol. 21, 733-744, 1976), an underlying dimensionality of two is commonly assigned to the LAC of human tissue at diagnostic energies. The finding of a higher statistically significant dimensionality thus raises the question whether a higher assumed model dimensionality (now feasible with the advent of multibin x-ray systems) might also be practically relevant, i.e., if better tissue characterization results can be obtained.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Simple analytical model for fitting QENS data from liquids

The advent of new generation of neutron spectrometers for quasielastic scattering measurements has demonstrated the need for multi-component (at least two-component) fitting of data collected from liquids. The narrowest component is associated with long-range translational diffusion, whereas the broader component originates from localized center-of-mass motions in the transient cage of the nearest neighbors. Intra-particle dynamics (e.g., of molecular side groups) may give rise to additional components besides the center-of-mass motions. We present a simple analytical model for two-component fits of quasielastic neutron scattering data collected from liquids. Here, the use of such models is optional (albeit helpful) in the traditional approach, where fits are performed sequentially and independently for the data collected at separate Q values, but imperative in the global fitting, where the data are fitted simultaneously at all Q values, which may be necessary to rationalize the data from increasingly complex systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Does the presence of planets affect the frequency and properties of extrasolar Kuiper belts? Results from the HErschel DEBRIS and DUNES surveys

The study of the planet–debris disk connection can shed light on the formation and evolution of planetary systems and may help “predict” the presence of planets around stars with certain disk characteristics. In preliminary analyses of subsamples of the Herschel DEBRIS and DUNES surveys, Wyatt et al. and Marshall et al. identified a tentative correlation between debris and the presence of low-mass planets. Here we use the cleanest possible sample out of these Herschel surveys to assess the presence of such a correlation, discarding stars without known ages, with ages <1 Gyr, and with binary companions <100 AU to rule out possible correlations due to effects other than planet presence. In our resulting subsample of 204 FGK stars, we do not find evidence that debris disks are more common or more dusty around stars harboring high-mass or low-mass planets compared to a control sample without identified planets. There is no evidence either that the characteristic dust temperature of the debris disks around planet-bearing stars is any different from that in debris disks without identified planets, nor that debris disks are more or less common (or more or less dusty) around stars harboring multiple planets compared to single-planet systems. Diverse dynamical histories may account for the lack of correlations. The data show a correlation between the presence of high-mass planets and stellar metallicity, but no correlation between the presence of low-mass planets or debris and stellar metallicity. Comparing the observed cumulative distribution of fractional luminosity to those expected from a Gaussian distribution in logarithmic scale, we find that a distribution centered on the solar system’s value fits the data well, while one centered at 10 times this value can be rejected. This is of interest in the context of future terrestrial planet detection and characterization because it indicates that there are good prospects for finding a large number of debris disk systems (i.e., with evidence of harboring planetesimals, the building blocks of planets) with exozodiacal emission low enough to be appropriate targets for an ATLAST-type mission to search for biosignatures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

$J/\psi $-p scattering length from the total and differential photoproduction cross sections

The J/ψ-p scattering length, αJ/ψp, can be extracted from the J/ψ photoproduction cross section near threshold using the Vector Meson Dominance (VMD) model to relate the reaction γp→J/ψp to J/ψp→J/ψp. Such estimates based on experimental data result in values for |αJ/ψp|, which are much lower than most of the theoretical predictions. In this work, we study the relations between the different results, depending on the use of the total or the differential cross sections, and the method of extrapolating the data to threshold in the case of a low-statistics data sample, such as the near threshold J/ψ photoproduction dataset. We estimate a range for |αJ/ψp| of 0.003–0.025 fm as extracted from experimental data within the VMD model and discuss possible reasons for such lower values compared to the theoretical results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Landmark-Warped Emulators for Models with Misaligned Functional Response

Many computer models output functional data, and in some cases, these functional data have similar, but misaligned, shape characteristics. In this paper, we introduce a general approach for building emulators for computer models that output misaligned functional data when key values in the functional response (landmarks) can be easily identified. This approach has two main parts: modeling the aligned (using the landmarks) functional data, and modeling the functions that map the misaligned data to the aligned space (warping functions). As the warping functions are required to be monotonic, we give special attention to modeling monotonic functional response data. We discuss how our approach can be easily applied for a variety of typical emulators, such as Gaussian processes, Bayesian multivariate adaptive regression splines, and Bayesian additive regression trees, and how sensitivity analysis can be performed. We demonstrate our approach by building emulators for two applications: (1) a high-energy-density physics computer model used to simulate inertial confinement fusion ignition experiments, where model outputs are highly misaligned, and (2) a multiphysics continuum hydrocode used to simulate high-velocity impact experiments, where model outputs are only slightly misaligned. In case (1) traditional methods cannot be applied, while in (2) they can be applied, but the proposed method performs significantly better.

97 MATHEMATICS AND COMPUTING↗