Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model Inference”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research.Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra: Modeling Archive

This dataset is a model archive of the paper Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra (in prep) to support a modeling study investigating how projected increases in Arctic temperature and precipitation will jointly influence hydrologic conditions in ice-rich tundra landscapes. With this dataset, this study is to address the research question: Will Arctic tundra landscapes become wetter or drier with increasing precipitation and temperature in the future when thaw-induced ground subsidence and associated microtopographic evolution are represented? The simulations focus on ice-wedge polygon tundra, a widespread form of ice-rich permafrost terrain that is highly sensitive to thaw-driven landscape change. This dataset contains model input and output data for four study watersheds in Alaska: Anaktuvuk, Utqiagvik (formerly Barrow), Brooks Foothills, and Prudhoe Bay. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.5), a physics-rich integrated surface–subsurface hydrologic model. For each watershed, ten modeling cases were performed representing two landscape evolution conditions (with subsidence and without subsidence) combined with five climate forcing scenarios derived from Shared Socioeconomic Pathways (SSP5, SSP5 with precipitation trend, SSP2, SSP2 with precipitation trend, and SSP2 with double precipitation trend). Particularly, for each watershed under the forcing SSP2 with precipitation trend, there are two additional simulations considering spatially heterogeneous subsidence distributions: one assumes randomly distributed scaling and the other includes elevation dependent distribution scaling. These simulations span 1980–2099 and include spin-up runs (1980–2009) followed by transient projections (2010–2099). To facilitate reproducibility of simulations, all datasets are organized by watershed. For each study watershed, the dataset contains: (1) Pre-partitioned mesh files for 32-core modeling (.par.32.XX), located in EACH_WATERSHED/mesh/basin; and also a non-partitioned mesh file (.exo) located in EACH_WATERSHED/mesh; (2) Climate forcings corresponding to the five SSP scenarios (.h5), located in EACH_WATERSHED/data; (3) Final states (.h5) from column spin-up modeling used to initialize historical watershed-scale spin-up runs from 1980 to 2009, located in EACH_WATERSHED/PreSpinupHistorical; (4) Final states (.h5) of historical watershed-scale spin-up runs from 1980 to 2009 used to initialize projection runs, located in EACH_WATERSHED/Spinup_daymetERA5; (5) ATS modeling input files (.xml), located in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/inputfiles; (6) ATS modeling output files (.dat), located in in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/combined_obs; (7) For the Brooks Foothills watershed, additional spatial model outputs are provided (.h5) for selected years (2033 and 2093) used to generate spatial figures in this study, located in Brooksfoothills/EACH_SIMULATION_SCENARIO/results-WITH/WITHOUT_SUBSIDENCE-year2033/2093. All data files with suffix .h5 can be accessible through Python h5py, and all data files with suffix of .dat can be imported by Python pandas. Mesh file with .exo can be visualized through Paraview or read by Python netCDF. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION↗
DOE OSTI2026

Multi-resolution Arctic Shrub Cover Dataset Derived from UAS and Airborne SfM and LiDAR (2013-2025)

We synthesized 177 unoccupied aerial system flights and 77 airborne flights across the Arctic and created a multi-resolution benchmark data of low-to-tall shrub fractional cover leveraging Structure-from-Motion and Light Detection and Ranging. The resulting dataset covered a total of 1899 km2 across Alaska, Western Canada, Sweden, and Siberian Arctic, including key sites from the Oro Arctic to the High Arctic. The dataset is organized into 6 primary data collection directories (“Abisko,” “AWI,” “ERE,” “Fairbanks,” “NGEE,” “Toolik”), each containing site and flight subdirectories. Flight directories include shrub cover rasters (*.tifs) at 1 m, 5 m, and 30 m resolution, the canopy height model at 1 m resolution (*.tifs), and a bounding box *.kml file. For the AWI, Abisko, NGEE, and Fairbanks collections, we also include the GCC raster at 1 m resolution (*.tif). Files are organized by Collection > Site > Flight Name > Data Files. Flight rasters are in the local UTM zone and the .kml files are in the geographic coordinate system EPSG 4326. We also include a .csv file that details the source datasets for every flight. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

canopy height model↗
DOE OSTI2026

BatFIT (Battery Feature Inference Toolbox) [SWR-26-034]

This package implements several data-based techniques for parameter fitting in Li-ion battery models. It uses BATMODS-lite to generate the data . The repository contains the code that is used for the paper "Neural posterior estimation is accurate, tractable and scalable for inverse parameter inference in Li-ion batteries", M. Hassanaly, C. R. Randall, P. J. Weddle, P. J. Gasper, C. Kelly, T. R. Tanim, K. Smith.

Hassanaly, Malik [National Laboratory of the Rocki↗
DOE OSTI2015

A variance-based decomposition and global sensitivity index method for uncertainty quantification: Application to retrieved ice cloud properties

This study develops a novel uncertainty quantification (UQ) method for cloud microphysical property retrievals using variance-based decomposition and global sensitivity index. In this UQ framework, empirical orthogonal function (EOF) analysis is applied to the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) ground-based observations, which are the inputs for the cloud retrieval studied here. The principal components (PCs) in the EOF expansion are parameterized as random input variables, and hence, the input dimension is greatly reduced (up to a factor of 50), allowing large ensemble of random samplings. The EOF expansion improves the accuracy of the uncertainty estimation by taking into account the cross correlations in the input data profiles. Here, this method enables a probabilistic representation of a retrieval process by adding normally distributed perturbations into PCs of sample means of input data profiles within a time window. Therefore, it effectively facilitates objective validation of climate models against cloud retrievals under a probabilistic framework for rigorous statistical inferences. Moreover, the variance-based global sensitivity index analysis, part of this method, attributes the output uncertainties to each individual source, thus providing directions for improving retrieval algorithms and observation strategies. For demonstration, we apply this method to quantify the uncertainties of the ARM program's baseline cloud retrieval algorithm for an ice cloud case observed at the Southern Great Plains site on 9 March 2000.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Probabilistic Mixture Model-Based Spectral Unmixing

Spectral unmixing attempts to decompose a spectral ensemble into the constituent pure spectral signatures (called endmembers) along with the proportion of each endmember. This is essential for techniques like hyperspectral imaging (HSI) used in environment monitoring, geological exploration, etc. Several spectral unmixing approaches have been proposed, many of which are connected to hyperspectral imaging. However, most extant approaches assume highly diverse collections of mixtures and extremely low-loss spectroscopic measurements. Additionally, current non-Bayesian frameworks do not incorporate the uncertainty inherent in unmixing. We propose a probabilistic inference algorithm that explicitly incorporates noise and uncertainty, enabling us to unmix endmembers in collections of mixtures with limited diversity. We use a Bayesian mixture model to jointly extract endmember spectra and mixing parameters while explicitly modeling observation noise and the resulting inference uncertainties. We obtain approximate distributions over endmember coordinates for each set of observed spectra while remaining robust to inference biases from the lack of pure observations and the presence of non-isotropic Gaussian noise. As a direct impact of our methodology, access to reliable uncertainties on the unmixing solutions would enable robust solutions to noise, as well as informed decision-making for HSI applications and other unmixing problems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

AI and extreme scale computing to learn and infer the physics of higher order gravitational wave modes of quasi-circular, spinning, non-precessing black hole mergers

We use artificial intelligence (AI) to learn and infer the physics of higher order gravitational wave modes of quasi-circular, spinning, non precessing binary black hole mergers. We trained AI models using 14 million waveforms, produced with the surrogate model NRHybSur3dq8, that include modes up to $\ell$ ≤ 4 and (5,5), except for (4,0) and (4,1), that describe binaries with mass-ratios $\textit{q}$ ≤ 8, individual spins $s^z_{\{1,2\}} \in$[–0.8,0.8], and inclination angle $θ \in$ [0,π]. Our probabilistic AI surrogates can accurately constrain the mass-ratio, individual spins, effective spin, and inclination angle of numerical relativity waveforms that describe such signal manifold. We compared the predictions of our AI models with Gaussian process regression, random forest, k-nearest neighbors, and linear regression, and with traditional Bayesian inference methods through the PyCBC Inference toolkit, finding that AI outperforms all these approaches in terms of accuracy, and are between three to four orders of magnitude faster than traditional Bayesian inference methods. Our AI surrogates were trained within 3.4 hours using distributed training on 1,536 NVIDIA V100 GPUs in the Summit supercomputer.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Infall-driven protostellar accretion and the solution to the luminosity problem

We investigate the role of mass infall in the formation and evolution of protostars. To avoid ad hoc initial and boundary conditions, we consider the infall resulting self-consistently from modeling the formation of stellar clusters in turbulent molecular clouds. We show that infall rates in turbulent clouds are comparable to accretion rates inferred from protostellar luminosities or measured in pre-main-sequence stars. They should not be neglected in modeling the luminosity of protostars and the evolution of disks, even after the embedded protostellar phase. We find large variations of infall rates from protostar to protostar, and large fluctuations during the evolution of individual protostars. In most cases, the infall rate is initially of order 10{sup –5} M {sub ☉} yr{sup –1}, and may either decay rapidly in the formation of low-mass stars, or remain relatively large when more massive stars are formed. The simulation reproduces well the observed characteristic values and scatter of protostellar luminosities and matches the observed protostellar luminosity function. The luminosity problem is therefore solved once realistic protostellar infall histories are accounted for, with no need for extreme accretion episodes. These results are based on a simulation of randomly driven magnetohydrodynamic turbulence on a scale of 4 pc, including self-gravity, adaptive-mesh refinement to a resolution of 50 AU, and accreting sink particles. The simulation yields a low star formation rate, consistent with the observations, and a mass distribution of sink particles consistent with the observed stellar initial mass function during the whole duration of the simulation, forming nearly 1300 sink particles over 3.2 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Conductivity and dissociation in liquid metallic hydrogen and implications for planetary interiors

Liquid metallic hydrogen (LMH) is the most abundant form of condensed matter in our solar planetary structure. The electronic and thermal transport properties of this metallic fluid are of fundamental interest to understanding hydrogen’s mechanism of conduction, atomic or pairing structure, as well as the key input for the magnetic dynamo action and thermal models of gas giants. Here, we report spectrally resolved measurements of the optical reflectance of LMH in the pressure region of 1.4–1.7 Mbar. We analyze the data, as well as previously reported measurements, using the free-electron model. Fitting the energy dependence of the reflectance data yields a dissociation fraction of 65 ± 15%, supporting theoretical models that LMH is an atomic metallic liquid. We determine the optical conductivity of LMH and find metallic hydrogen’s static electrical conductivity to be 11,000–15,000 S/cm, substantially higher than the only earlier reported experimental values. Here, the higher electrical conductivity implies that the Jovian and Saturnian dynamo are likely to operate out to shallower depths than previously assumed, while the inferred thermal conductivity should provide a crucial experimental constraint to heat transport models.

08 HYDROGEN↗
DOE OSTI2015

Investigating Hα, UV, and IR star-formation rate diagnostics for a large sample of z ∼ 2 galaxies

We use a sample of 262 spectroscopically confirmed star-forming galaxies at redshifts 2.08⩽z⩽2.51 to compare Hα, ultraviolet (UV), and IR star formation rate (SFR) diagnostics and to investigate the dust properties of the galaxies. At these redshifts, the Hα line shifts to the K{sub s} band. By comparing K{sub s}-band photometry to underlying stellar population model fits to other UV, optical, and near-infrared data, we infer the Hα flux for each galaxy. We obtain the best agreement between Hα- and UV-based SFRs if we assume that the ionized gas and stellar continuum are reddened by the same value and that the Calzetti attenuation curve is applied to both. Aided with MIPS 24 μm data, we find that an attenuation curve steeper than the Calzetti curve is needed to reproduce the observed IR/UV ratios of galaxies younger than 100 Myr. Furthermore, using the bolometric SFR inferred from the UV and mid-IR data (SFR{sub IR}+SFR{sub UV}), we calculated the conversion between the Hα luminosity and SFR to be (7.5±1.3)×10{sup −42} for a Salpeter initial mass function, which is consistent with the Kennicutt conversion. The derived conversion factor is independent of any assumption of the dust correction and is robust to stellar population model uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

AN EXPLANATION FOR THE DIFFERENT X-RAY TO OPTICAL COLUMN DENSITIES IN THE ENVIRONMENTS OF GAMMA RAY BURSTS: A PROGENITOR EMBEDDED IN A DENSE MEDIUM

We study the {approx}> 10 ratios in the X-ray to optical column densities inferred from afterglow spectra of gamma ray bursts (GRBs) due to gas surrounding their progenitors. We present time-evolving photoionization calculations for these afterglows and explore different conditions of their environment. We find that homogenous models of the environment (constant density) predict X-ray columns similar to those found in the optical spectra, with the bulk of the opacity being produced by neutral material at large distances from the burst. This result is independent of gas density or metallicity. Only models assuming a progenitor immersed in a dense ({approx}10{sup 2-4} cm{sup -3}) cloud of gas (with radius {approx}10 pc), with a strong, declining gradient of density for the surrounding interstellar medium (ISM) are able to account for the large X-ray to optical column density ratios. However, to avoid an unphysical correlation between the size of this cloud and the size of the ionization front produced by the GRB, the models also require that the circumburst medium is already ionized prior to the burst. The inferred cloud masses are {approx}< 10{sup 6} M{sub Sun }, even if low metallicities in the medium are assumed (Z {approx} 0.1 Z{sub Sun }). These cloud properties are consistent with those found in giant molecular clouds and our results support a scenario in which the progenitors reside within intense star formation regions of galaxies. Finally, we show that modeling over large samples of GRB afterglows may offer strong constraints on the range of properties in these clouds, and the host galaxy ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Bayesian merging of numerical modeling and remote sensing for saltwater intrusion quantification in the Vietnamese Mekong Delta

Saltwater intrusion has become one of the most concerning issues in the Vietnamese Mekong Delta (VMD) due to its increasing impacts on agriculture and food security of Vietnam. Reliable estimation of salinity plays a crucial role to mitigate the impacts of saltwater intrusion. Here, this study developed a hybrid technique that merges satellite imagery with numerical simulations to improve the estimation of salinity in the VMD. The salinity derived from Landsat images and by numerical simulations was fused using the Bayesian inference technique. The results indicate that our technique significantly reduces the uncertainties and improves the accuracy of salinity estimates. The Nash–Sutcliffe coefficient is 0.74, which is much higher than that of numerical simulation (0.63) and Landsat estimation (0.6). The correlation coefficient between the ensemble and measured salinity is relatively high (0.88). The variance of the ensemble salinity errors (5.0 ppt 2 ) is lower than that of Landsat estimation (10.4 ppt 2 ) and numerical simulations (9.6 ppt 2 ). The proposed approach shows a great potential to combine multiple data sources of a variable of interest to improve its accuracy and reliability wherever these data are available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008

The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedure of time-delay cosmography by comparing independent and blind time-delay predictions of the system WGD 2038-4008 from two teams using two different software programs: GLEE and LENSTRONOMY. The predicted time delays from the two teams incorporate the stellar kinematics of the deflector and the external convergence from line-of-sight structures. The un-blinded time-delay predictions from the two teams agree within 1.2σ, implying that once the time delay is measured the inferred Hubble constant will also be mutually consistent. However, there is a ~4σ discrepancy between the power-law model slope and external shear, which is a significant discrepancy at the level of lens models before the stellar kinematics and the external convergence are incorporated. We identify the difference in the reconstructed point spread function (PSF) to be the source of this discrepancy. When the same reconstructed PSF was used by both teams, we achieved excellent agreement, within ~0.6σ<, indicating that potential systematics stemming from source reconstruction algorithms and investigator choices are well under control. We recommend that future studies supersample the PSF as needed and marginalize over multiple algorithms or realizations for the PSF reconstruction to mitigate the systematics associated with the PSF. A future study will measure the time delays of the system WGD 2038-4008 and infer the Hubble constant based on our mass models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Sparsifying priors for Bayesian uncertainty quantification in model discovery

We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Constraining physical models at gigabar pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative infor mation. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. As a result, the goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

SU-F-T-132: Variable RBE Models Predict Possible Underestimation of Vaginal Dose for Anal Cancer Patients Treated Using Single-Field Proton Treatments

Purpose: Anal cancer patients treated using a posterior proton beam may be at risk of vaginal wall injury due to the increased linear energy transfer (LET) and relative biological effectiveness (RBE) at the beam distal edge. We investigate the vaginal dose received. Methods: Five patients treated for anal cancer with proton pencil beam scanning were considered, all treated to a prescription dose of 54 Gy(RBE) over 28–30 fractions. Dose and LET distributions were calculated using the Monte Carlo simulation toolkit TOPAS. In addition to the standard assumption of a fixed RBE of 1.1, variable RBE was considered via the application of published models. Dose volume histograms (DVHs) were extracted for the planning treatment volume (PTV) and vagina, the latter being used to calculate the vaginal normal tissue complication probability (NTCP). Results: Compared to the assumption of a fixed RBE of 1.1, the variable RBE model predicts a dose increase of approximately 3.3 ± 1.7 Gy at the end of beam range. NTCP parameters for the vagina are incomplete in the current literature, however, inferring value ranges from the existing data we use D{sub 50} = 50 Gy and LKB model parameters a=1–2 and m=0.2–0.4. We estimate the NTCP for the vagina to be 37–48% and 42–47% for the fixed and variable RBE cases, respectively. Additionally, a difference in the dose distribution was observed between the analytical calculation and Monte Carlo methods. We find that the target dose is overestimated on average by approximately 1–2%. Conclusion: For patients treated with posterior beams, the vaginal wall may coincide with the distal end of the proton beam and may receive a substantial increase in dose if variable RBE models are applied compared to using the current clinical standard of RBE equal to 1.1. This could potentially lead to underestimating toxicities when treating with protons.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Constraining Physical Models at Gigabar Pressures

High-energy-density (HED) experiments in convergent geometry are able to test physical models at pressures beyond hundreds of millions of atmospheres. The measurements from these experiments are generally highly integrated and require unique analysis techniques to procure quantitative information. This work describes a methodology to constrain the physics in convergent HED experiments by adapting the methods common to many other fields of physics. As an example, a mechanical model of an imploding shell is constrained by data from a thin-shelled direct-drive exploding-pusher experiment on the OMEGA Laser System using Bayesian inference, resulting in the reconstruction of the shell dynamics and energy transfer during the implosion. The model is tested by analyzing synthetic data from a 1-D hydrodynamics code and is sampled using a Markov chain Monte Carlo to generate the posterior distributions of the model parameters. The goal of this work is to demonstrate a general methodology that can be used to draw conclusions from a wide variety of HED experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Contributions of Different Cloud Types to Feedbacks and Rapid Adjustments in CMIP5*

When using five climate model simulations of the response to an abrupt quadrupling of CO 2 , the authors perform the first simultaneous model intercomparison of cloud feedbacks and rapid radiative adjustments with cloud masking effects removed, partitioned among changes in cloud types and gross cloud properties. After CO 2 quadrupling, clouds exhibit a rapid reduction in fractional coverage, cloud-top pressure, and optical depth, with each contributing equally to a 1.1 W m -2 net cloud radiative adjustment, primarily from shortwave radiation. Rapid reductions in midlevel clouds and optically thick clouds are important in reducing planetary albedo in every model. As the planet warms, clouds become fewer, higher, and thicker, and global mean net cloud feedback is positive in all but one model and results primarily from increased trapping of longwave radiation. As was true for earlier models, high cloud changes are the largest contributor to intermodel spread in longwave and shortwave cloud feedbacks, but low cloud changes are the largest contributor to the mean and spread in net cloud feedback. The importance of the negative optical depth feedback relative to the amount feedback at high latitudes is even more marked than in earlier models. Furthermore, the authors show that the negative longwave cloud adjustment inferred in previous studies is primarily caused by a 1.3 W m -2 cloud masking of CO 2 forcing. Properly accounting for cloud masking increases net cloud feedback by 0.3 W m -2 K -1 , whereas accounting for rapid adjustments reduces by 0.14 W m -2 K -1 the ensemble mean net cloud feedback through a combination of smaller positive cloud amount and altitude feedbacks and larger negative optical depth feedbacks.

58 GEOSCIENCES↗