Search NASA⌕ Search

SEARCH · Search NASA

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Vitality models found useful in modeling tag-failure times in acoustic-tag survival studies

Abstract Acoustic telemetry studies often rely on the assumption that premature tag failure does not affect the validity of inferences. However, in some cases this assumption is possibly or likely invalid and it is necessary to apply a correction to estimation procedures. The question of which approaches and specific models are best suited to modeling acoustic tag failures has received little research attention. In this short communication, we present a meta-analysis of 42 acoustic tag-life studies, originally used to correct survival studies involving outmigrating juvenile salmonids in the Columbia/Snake river basin. We compare the performance of nine alternative parametric models including common failure–time/survival models and the vitality models of Li and Anderson Theor Popul Biol 76:118–131, (2009) and Demogr Res 28:341–372, (2013). The tag-life studies used acoustic tags from three different tag manufacturers, had expected lifetimes between 12 and 61 days, and had dry weights ranging from 0.22 to 1.65 g. In 57% of the cases, the vitality models of Li and Anderson Theor Popul Biol 76:118–131, (2009) and Demogr Res 28:341–372, (2013) fit the tag-failure times best. The vitality models were also the second-best choices in 17% of the cases. Together, the vitality models, log-logistic, (19%), and gamma models (14%) accounted for 90% of the models selected. Unlike more traditional failure–time models (e.g., Weibull, Gompertz, gamma, and log-logistic), the vitality models are capable of characterizing both the early onset of tag failure due to manufacturing errors and the anticipated battery life. We provide further guidance on appropriate sample sizes (50–100 tags) and procedures to be considered when applying precise tag-life corrections in release–recapture survival studies.

Skalski, John R. (ORCID:0000000270702505)↗
DOE OSTI2016

How Adverse Outcome Pathways Can Aid the Development and Use of Computational Prediction Models for Regulatory Toxicology

Efforts are underway to transform regulatory toxicology and chemical safety assessment from a largely empirical science based on direct observation of apical toxicity outcomes in whole organism toxicity tests to a predictive one in which outcomes and risk are inferred from accumulated mechanistic understanding. The adverse outcome pathway (AOP) framework has emerged as a systematic approach for organizing knowledge that supports such inference. We argue that this systematic organization of knowledge can inform and help direct the design and development of computational prediction models that can further enhance the utility of mechanistic and in silico data for chemical safety assessment. Examples of AOP-informed model development and its application to the assessment of chemicals for skin sensitization and multiple modes of endocrine disruption are provided. The role of problem formulation, not only as a critical phase of risk assessment, but also as guide for both AOP and complementary model development described. Finally, a proposal for actively engaging the modeling community in AOP-informed computational model development is made. The contents serve as a vision for how AOPs can be leveraged to facilitate development of computational prediction models needed to support the next generation of chemical safety assessment.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Large-scale inference of gene function through phylogenetic annotation of Gene Ontology terms: case study of the apoptosis and autophagy cellular processes

We previously reported a paradigm for large-scale phylogenomic analysis of gene families that takes advantage of the large corpus of experimentally supported Gene Ontology (GO) annotations. This ‘GO Phylogenetic Annotation’ approach integrates GO annotations from evolutionarily related genes across 100 different organisms in the context of a gene family tree, in which curators build an explicit model of the evolution of gene functions. GO Phylogenetic Annotation models the gain and loss of functions in a gene family tree, which is used to infer the functions of uncharacterized (or incompletely characterized) gene products, even for human proteins that are relatively well studied. Here, we report our results from applying this paradigm to two well-characterized cellular processes, apoptosis and autophagy. This revealed several important observations with respect to GO annotations and how they can be used for function inference. Notably, we applied only a small fraction of the experimentally supported GO annotations to infer function in other family members. The majority of other annotations describe indirect effects, phenotypes or results from high throughput experiments. In addition, we show here how feedback from phylogenetic annotation leads to significant improvements in the PANTHER trees, the GO annotations and GO itself. Thus GO phylogenetic annotation both increases the quantity and improves the accuracy of the GO annotations provided to the research community. We expect these phylogenetically based annotations to be of broad use in gene enrichment analysis as well as other applications of GO annotations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Development of a beam ion velocity detector for the heavy ion beam probe

In an axisymmetric plasma, the conservation of canonical angular momentum constrains heavy ion beam probe (HIBP) trajectories such that measurement of the toroidal velocity component of secondary ions provides a localized determination of the poloidal flux at the volume where they originated. We have developed a prototype detector which is designed to determine the beam angle in one dimension through the detection of ion current landing on two parallel planes of detecting elements. A set of apertures creates a pattern of ion current on wires in the first plane and solid metal plates behind them; the relative amounts detected by the wires and plates determine the angle which beam ions enter the detector, which is used to infer the toroidal velocity component. The design evolved from a series of simulations within which we modeled ion beam velocity changes due to equilibrium and fluctuating magnetic fields, along with the ion beam profile and velocity dispersion, and studied how these and characteristics such as the size, cross section, and spacing of the detector elements affect performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Using onset times from frequent seismic surveys to understand fluid flow at the Peace River Field, Canada

With our limited knowledge of the relationship between changes in the state of an aquifer or reservoir and the corresponding changes in the elastic moduli, that is the rock physics model, hampers the effective use of time-lapse seismic observations for estimating flow properties within the Earth. A central problem is the complicated dependence of the magnitude of time-lapse changes on the saturation, pressure, and temperature changes within an aquifer or reservoir. We describe an inversion methodology for reservoir characterization that uses onset times, the calendar time of the change in seismic attributes, rather than the magnitude of the changes. We find that onset times are much less sensitive than magnitudes to the rock physics model used to relate time-lapse observations to changes in saturation, temperature and fluid pressure. We apply the inversion scheme to observations from daily monitoring of enhanced oil recovery at the Peace River field in Canada. An array of 1492 buried hydrophones record seismic signals from 49 buried sources. Time-shifts for elastic waves traversing the reservoir are extracted from the daily time-lapse cubes. In our analysis 175 images of time-shifts are transformed into a single map of onset times, leading to a substantial reduction in the volume of data. These observations are used in conjunction with bottom hole pressure data to infer the initial conditions prior to the injection, and to update the reservoir permeability model. The combination of a global and local inversion scheme produces a collection of reservoir models that are best described by three clusters. The updated model leads to a nearly 70 percent reduction in seismic data misfit. The final set of solutions successfully predict the observed normalized pressure history during the soak and flow-back into the wells between 82 and 175 days into the cyclic steaming operation.

58 GEOSCIENCES↗
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
DOE OSTI2021

Recognition Method of Vehicle Cluster Situation Based on Set Pair Logic considering Driver’s Cognition

The recognition of vehicle cluster situations is one of the critical technologies of advanced driving, such as intelligent driving and automated driving. The accurate recognition of vehicle cluster situations is helpful for behavior decision‐making safe and efficient. In order to accurately and objectively identify the vehicle cluster situation, a vehicle cluster situation model is proposed based on the interval number of set pair logic. The proposed model can express the traffic environment’s knowledge considering each vehicle’s characteristics, grouping relationships, and traffic flow characteristics in the target vehicle’s interest region. A recognition method of vehicle cluster situation is designed to infer the traffic environment and driving conditions based on the connection number of set pair logic. In the proposed model, the uncertainty of the driver’s cognition is fully considered. In the recognition method, the relative uncertainty and relative certainty of driver’s cognition, traffic information, and vehicle cluster situation are fully considered. The verification results show that the proposed recognition method of vehicle cluster situations can realize accurate and objective recognition. The proposed anthropomorphic recognition method could provide a basis for vehicle autonomous behavior decision‐making.

Liu, Shijie↗
DOE OSTI2017

Beam-target double-spin asymmetry in quasielastic electron scattering off the deuteron with CLAS

The deuteron plays a pivotal role in nuclear and hadronic physics, as both the simplest bound multinucleon system and as an effective neutron target. Quasielastic electron scattering on the deuteron is a benchmark reaction to test our understanding of deuteron structure and the properties and interactions of the two nucleons bound in the deuteron. The experimental data presented here can be used to test state-of-the-art models of the deuteron and the two-nucleon interaction in the final state after two-body breakup of the deuteron. Focusing on polarization degrees of freedom, we gain information on spin-momentum correlations in the deuteron ground state (due to the D-state admixture) and on the limits of the impulse approximation (IA) picture as it applies to measurements of spin-dependent observables like spin structure functions for bound nucleons. Information on this reaction can also be used to reduce systematic uncertainties on the determination of neutron form factors or deuteron polarization through quasielastic polarized electron scattering. Furthermore, we measured the beam-target double-spin asymmetry (A || ) for quasielastic electron scattering off the deuteron at several beam energies (1.6–1.7, 2.5, 4.2, and 5.6–5.8GeV), using the CEBAF Large Acceptance Spectrometer (CLAS) at the Thomas Jefferson National Accelerator Facility. The deuterons were polarized along (or opposite to) the beam direction. The double-spin asymmetries were measured as a function of photon virtuality Q 2 (0.13–3.17(GeV/c) 2 ), missing momentum (p m =0.0–0.5GeV/c), and the angle between the (inferred) spectator neutron and the momentum transfer direction (θ nq ). We compare our results with a recent model that includes final-state interactions (FSI) using a complete parametrization of nucleon-nucleon scattering, as well as a simplified model using the plane wave impulse approximation (PWIA). We find overall good agreement with both the PWIA and FSI expectations at low to medium missing momenta (p m ≤0.25GeV/c), including the change of the asymmetry due to the contribution of the deuteron D state at higher momenta. At the highest missing momenta, our data clearly agree better with the calculations including FSI. Final-state interactions seem to play a lesser role for polarization observables in deuteron two-body electrodisintegration than for absolute cross sections. Our data, while limited in statistical power, indicate that PWIA models work reasonably well to understand the asymmetries at lower missing momenta. In turn, this information can be used to extract the product of beam and target polarization (P b P t ) from quasielastic electron-deuteron scattering, which is useful for measurements of spin observables in electron-neutron inelastic scattering. But, at the highest missing (neutron) momenta, FSI effects become important and must be accounted for.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Uncovering New Physics in the Cosmic Microwave Background: Developing Novel Theoretical Models and Machine-Learning-Powered Constraints

The search for evidence of new physics via its signatures in the cosmos is a cornerstone goal of the Office of High Energy Physics in the DOE Office of Science. Indeed, the current concordance cosmological model provides intriguing hints for beyond-the-standard-model (BSM) physics, such as dark matter and dark energy. Recently, a potential breakdown has appeared in this model, which could be initial evidence toward a further important revision in our fundamental theoretical understanding of cosmology. This breakdown is reflected in disagreements between inferences of the current expansion rate of the universe, H 0 (the Hubble constant), based on indirect, cosmological data (e.g., from the early universe) and based on direct, local measurements. Despite significant effort, a compelling new concordance cosmological model has yet to be found; achieving significant progress on this front was the first major focus of the project. Theoretical considerations indicate that if the observational discrepancies are not due to systematic errors, they strongly suggest new physics operating in the redshift range just prior to recombination, when cosmic microwave background (CMB) photons last scattered. Crucially, almost all such models produce unique signatures in the CMB temperature and polarization power spectra, which will be measured with unprecedented precision by ongoing and upcoming experiments, including the DOE-supported CMB-S4 project. However, these subtle hints of new physics must be uncovered from beneath a swath of Galactic and extragalactic foreground contamination. Current CMB analysis methods, although powerful, do not optimally infer the CMB power spectrum in the presence of non-Gaussian foregrounds. There is thus scope for theoretical improvement in this foundational challenge of cosmological inference, which formed the second major focus of the project. The primary objectives of the project were two-fold: (1) to develop new theoretical models in cosmology that can restore concordance amongst the full suite of cosmological data sets, thereby potentially providing evidence of novel BSM physics; (2) to develop new theoretical machinery to enable significant sensitivity improvements in searches for new physics in cosmology, particularly via the CMB power spectrum. The two objectives are intertwined, as the analysis methodology improvements in (2) will enable the tightest possible constraints on the signatures of new physics predicted by the novel scenarios in (1). The theoretical approaches to restore concordance focused on models involving novel scalar field dynamics in the pre-recombination universe (the “early dark energy” scenario and modifications thereof), as well as couplings between this field and other components in the standard cosmological model, such as dark matter. We also studied a model featuring a generalization of the decaying dark matter scenario, in which a sub-component of dark matter converts into dark radiation at late times in cosmic history. While these ideas are mostly driven by phenomenological considerations, this tactic has proven extremely successful in cosmology throughout the past few decades, including in the early history of evidence for dark matter and dark energy. The new theoretical machinery envisioned in (2) is undergirded by developments in signal processing and machine learning, which will enable improvements in CMB power spectrum estimation in the presence of non-Gaussian foreground contaminants. In turn, this will yield optimal sensitivity in searches for new physics in the CMB, by maximizing the cosmological information that is extracted from this observable. The most ambitious outcome of this work would be the construction of a new cosmological model that restores concordance amongst data sets. Although the individual models studied here did not fully achieve that goal, significant progress in narrowing down the model space was made, as described below. Moreover, the outcome of the methodological improvements in (2) will significantly impact a wide range of theoretical cosmology, by enabling the tightest possible constraints on any model that leaves novel signatures in the CMB temperature and polarization power spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Hysteretic temperature sensitivity of wetland CH 4 fluxes explained by substrate availability and microbial activity

Methane (CH 4 ) emissions from wetlands are likely increasing and important in global climate change assessments. However, contemporary terrestrial biogeochemical model predictions of CH 4 emissions are very uncertain, at least in part due to prescribed temperature sensitivity of CH 4 production and emission. While statistically consistent apparent CH 4 emission temperature dependencies have been inferred from meta-analyses across microbial to ecosystem scales, year-round ecosystem-scale observations have contradicted that finding. Here, we show that apparent CH 4 emission temperature dependencies inferred from year-round chamber measurements exhibit substantial intra-seasonal variability, suggesting that using static temperature relations to predict CH 4 emissions is mechanistically flawed. Our model results indicate that such intra-seasonal variability is driven by substrate-mediated microbial and abiotic interactions: seasonal cycles in substrate availability favors CH 4 production later in the season, leading to hysteretic temperature sensitivity of CH 4 production and emission. Our findings demonstrate the uncertainty of inferring CH 4 emission or production rates from temperature alone and highlight the need to represent microbial and abiotic interactions in wetland biogeochemical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Ensemble Kalman Filter for Dynamic State Estimation of Power Grids Stochastically Driven by Time-correlated Mechanical Input Power

State and parameter estimation of power transmission networks is important for monitoring power grid operating conditions and analyzing transient stability. Wind power generation depends on fluctuating input power levels, which are correlated in time and contribute to uncertainty in turbine dynamical models. The ensemble Kalman filter (EnKF), a standard state estimation technique, uses a deterministic forecast and does not explicitly model time-correlated noise in parameters such as mechanical input power. However, this uncertainty affects the probability of fault-induced transient instability and increased prediction bias. Here a novel approach is to model input power noise with time-correlated stochastic fluctuations, and integrate them with the network dynamics during the forecast. While the EnKF has been used to calibrate constant parameters in turbine dynamical models, the calibration of a statistical model for a time-correlated parameter has not been investigated. In this study, twin experiments on a standard transmission network test case are used to validate our time-correlated noise model framework for state estimation of unsteady operating conditions and transient stability analysis, and a methodology is proposed for the inference of the mechanical input power time-correlation length parameter using time-series data from PMUs monitoring power dynamics at generator buses.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Strong Lensing Cosmology with Population-level Calibrated Neural Ratio Estimation

Strong gravitational lensing contains key information about cosmic acceleration. Modern and next-generation galaxy imaging surveys are expected to provide high-quality data on $\mathcal{O}(10^5)$ galaxy-galaxy lensing systems. The plethora and complexity of the data are likely to present computational challenges for parameter inference methods for fitting high-dimensional likelihoods, which are often analytically intractable. Neural Ratio Estimation (NRE) efficiently computes individual likelihood ratios that can be combined into population-level posteriors. We use simulations to study the capacity of NRE to jointly predict the dark energy equation-of-state parameter $w$ and the total matter density $Ω_{m}$ from lensing images and companion spectroscopic information. We also introduce a post hoc posterior coverage calibration procedure that mitigates the model overconfidence that is typically found in neural density estimation applications. Our experiments show that the errors on both parameters decrease with increasing inference population sizes. In particular, for 100 lenses in a standard $Λ$CDM Universe, our calibrated NRE model achieves median fractional uncertainty of $22.8\%$ in $w$ and $2.9\%$ in $Ω_{m}$. This proof of concept demonstrates a potentially scalable approach for efficient cosmological parameter inference with large populations of galaxy-scale lenses observed in future surveys.

Jarugula, Sreevani [Fermilab] (ORCID:0000000253867↗
DOE OSTI2019

A history of pore water oxygen isotope evolution in the Cretaceous Travis Peak Formation in East Texas

Oxygen isotope analyses of diagenetic cements can provide detailed evidence of sedimentary burial processes and conditions, as the δ 18 O values of precipitating minerals reflect contemporaneous local δ 18 O water and temperature conditions. Uncertainties in the timing and rates of pore water δ 18 O evolution in sedimentary basins can complicate interpretation of these records. Fracture-bridging (0.5–1 mm) quartz cements observed in sandstones of the Cretaceous Travis Peak Formation in the East Texas basin show clear growth-zoning by cathodoluminescence and contain detailed fluid inclusion records of temperature that make them excellent candidates for interrogating prolonged histories of basin temperature and the evolution of δ 18 O in basin pore water. New secondary ion mass spectrometer (SIMS) δ 18 O quartz isotopic data from fluid inclusion-rich quartz bridges in Travis Peak sandstones record a steady increase of pore water δ 18 O values from ~5 to 7‰ (VSMOW; Vienna Standard Mean Ocean Water) as the sandstone warms from ~130 to 150 °C. To help evaluate whether this trend could be generated solely from local water-rock interactions in response to burial compaction, a one-dimensional closed system isotopic burial model was created to simulate how δ 18 O water values change in a quartz-dominated sandstone during diagenesis. Using both directly measured and inferred rates of Travis Peak compaction, the magnitude of change in δ 18 O water that we calculate from quartz bridge geochemistry cannot be reasonably modeled solely by local quartz mechanical compaction, pressure solution, and cementation processes, necessitating significant fluxes of silica and high-δ 18 O water from outside of the sandstones prior to maximum burial. In conclusion, this indicates that even units which appear surrounded by significant barriers to fluid flow (i.e., mudrock-bounded channel sandstones) may have been infiltrated and diagenetically modified by large fluxes of fluid on geologic time scales.

58 GEOSCIENCES↗
DOE OSTI2018

Reconciling divergent trends and millennial variations in Holocene temperatures

Cooling during most of the past two millennia has been widely recognized1,2 and has been inferred to be the dominant global temperature trend of the past 11,700 years (the Holocene epoch)3. However, long-term cooling has been difficult to reconcile with global forcing4, and climate models consistently simulate long-term warming4. The divergence between simulations and reconstructions emerges primarily for northern mid-latitudes, for which pronounced cooling has been inferred from marine and coastal records using multiple approaches3. Here we show that temperatures reconstructed from sub-fossil pollen from 642 sites across North America and Europe closely match simulations, and that long-term warming, not cooling, defined the Holocene until around 2,000 years ago. The reconstructions indicate that evidence of long-term cooling was limited to North Atlantic records. Early Holocene temperatures on the continents were more than two degrees Celsius below those of the past two millennia, consistent with the simulated effects of remnant ice sheets in the climate model Community Climate System Model 3 (CCSM3)5. CCSM3 simulates increases in ‘growing degree days’—a measure of the accumulated warmth above five degrees Celsius per year—of more than 300 kelvin days over the Holocene, consistent with inferences from the pollen data. It also simulates a decrease in mean summer temperatures of more than two degrees Celsius, which correlates with reconstructed marine trends and highlights the potential importance of the different subseasonal sensitivities of the records. Despite the differing trends, pollen- and marine-based reconstructions are correlated at millennial-to-centennial scales, probably in response to ice-sheet and meltwater dynamics, and to stochastic dynamics similar to the temperature variations produced by CCSM3. Although our results depend on a single source of palaeoclimatic data (pollen) and a single climate-model simulation, they reinforce the notion that climate models can adequately simulate climates for periods other than the present-day. They also demonstrate that amplified warming in recent decades increased temperatures above the mean of any century during the past 11,000 years.

Science & Technology - Other Topics↗
DOE OSTI2023

Machine learning the spectral function of a hole in a quantum antiferromagnet

Understanding charge motion in a background of interacting quantum spins is a fundamental problem in quantum many-body physics. The most extensively studied model for this problem is the so-called t-t'-t''-J model, where the determination of the parameter t' in the context of cuprate superconductors is challenging. Here we present a theoretical study of the spectral functions of a mobile hole in the t-t'-t''-J model using two machine-learning techniques: K-nearest neighbor regression (KNN) and a feed-forward neural network (FFNN). We employ the self-consistent Born approximation to generate a dataset of about 1.3 x 10 5 spectral functions. Here we show that, for the forward problem, both methods allow for the accurate and efficient prediction of spectral functions, allowing, e.g., rapid searches through parameter space. Furthermore, we find that for the inverse problem (inferring Hamiltonian parameters from spectra), the FFNN can, but the KNN cannot, accurately predict the model parameters using merely the density of states. Our results suggest that it may be possible to use deep-learning methods to predict materials parameters from experimentally measured spectral functions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Vertical structure of magnetized accretion disks around young stars

We model the vertical structure of the magnetized accretion disks that are subject to viscous and resistive heating and irradiation by the central star. We apply our formalism to the radial structure of the magnetized accretion disks that are threaded by the poloidal magnetic field dragged during the process of star formation, which was developed by Shu and coworkers. We consider disks around low-mass protostars, T Tauri, and FU Orionis stars, as well as two levels of disk magnetization: λ{sub sys}=4 (strongly magnetized disks) and λ{sub sys}=12 (weakly magnetized disks). The rotation rates of strongly magnetized disks have large deviations from Keplerian rotation. In these models, resistive heating dominates the thermal structure for the FU Ori disk, and the T Tauri disk is very thin and cold because it is strongly compressed by magnetic pressure; it may be too thin compared with observations. Instead, in the weakly magnetized disks, rotation velocities are close to Keplerian, and resistive heating is always less than 7% of the viscous heating. In these models, the T Tauri disk has a larger aspect ratio, which is consistent with that inferred from observations. All the disks have spatially extended hot atmospheres where the irradiation flux is absorbed, although most of the mass (∼90%–95%) is in the disk midplane. With the advent of ALMA one expects direct measurements of magnetic fields and their morphology at disk scales. It will then be possible to determine the mass-to-flux ratio of magnetized accretion disks around young stars, an essential parameter for their structure and evolution. Our models contribute to the understanding of the vertical structure and emission of these disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

On the inference of viscoelastic constants from stress relaxation experiments

Several constitutive theories have been proposed in the literature to model the viscoelastic response of materials, including widely used rheological constitutive models. These models are characterized by certain parameters ("time constants") that define the time scales over which the material relaxes. These parameters are primarily obtained from stress relaxation experiments using curve-fitting techniques. However, the question of how best to estimate these time constants remains open. As a step towards answering this question, we propose an optimal experimental design approach based on ideas from information geometry, namely Fisher information and Kullback-Leibler divergence. The material is modeled as a spring element in parallel with multiple Maxwell elements and described using a one- or two-term Prony series. Treating the time constants as unknowns, we develop expressions for the Fisher information and Kullback-Leibler divergence that allow us to maximize information gain from experimental data. Based on the results of this study, we propose that the largest time constant estimated from a stress relaxation experiment for a linear viscoelastic material should be at most one-fifth of the total time of the experiment in order to maximize information gain. Our results also provide confirmation that the equilibrium modulus of the material cannot be reliably determined from curve-fitting to data from a stress relaxation experiment.

Fisher information↗
DOE OSTI2022

Bayesian Inference of the Fluctuating Proton Shape in DIS and Hadronic Collisions

We determine the likelihood distribution for the model parameters describing the event-by-event fluctuating proton geometry at small $x$ by performing a Bayesian analysis within the Color Glass Condensate framework. The exclusive J/ψ production data from HERA is found to constrain the model parameters well, and we demonstrate that complementary constraints can be obtained from simulations of Pb+Pb collisions at the LHC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗