Search NASA⌕ Search

SEARCH · Search NASA

Results for “parameter sensitivity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Determination of noise equivalent reflectance for a multispectral scanner: A scanner sensitivity study

The methods used to calculate the sensitivity parameter noise equivalent reflectance of a remote-sensing scanner are explored, and the results are compared with values measured over calibrated test sites. Data were acquired on four occasions covering a span of 4 years and providing various atmospheric conditions. One of the calculated values was based on assumed atmospheric conditions, whereas two others were based on atmospheric models. Results indicate that the assumed atmospheric conditions provide useful answers adequate for many purposes. A nomograph was developed to indicate sensitivity variations due to geographic location, time of day, and season.

Gibbons, D. E.↗
NASA NTRS1974

Influence of propulsion system size, shape, and location on supersonic aircraft design

The effects of various propulsion system parameters on the characteristics of a supersonic transport were investigated. The effects of arbitrarily scaling engine size on wave drag, friction drag, drag-due-to-lift, wing sizing, airplane balance, and airplane weight were studied. These evaluations were made for two families of nacelle shapes, resulting from typical turbojet and turbofan installations. Also examined were effects of nacelle location, and the wing camber plane deformations required to cancel the nacelle interference pressure field at cruise Mach number (2.7 M) were determined. The most drag-sensitive parameter is found to be nacelle shape. Similarly, wing deformation requirements are found to be primarily affected by nacelle shape. Effects of engine size variations are noted primarily in airplane gross weight.

Bonner, E.↗
DOE OSTI2026

Reweighting Underlying Event and Colour Reconnection parameter variations in Sherpa

We propose and validate a new method to trace the impact of parameter variations in the simulation of multi-parton interactions and colour reconnections in the Sherpa event generator. They are reflected, at an event-by-event basis, through relative weights with respect to the central production parameters that give rise to the generated events and distributions. Our method facilitates the tuning of the Monte Carlo event generator at a dramatically reduced computational cost, alleviates parameter sensitivity studies, and enables robust quantification of parametric uncertainties on-the-fly, one of the missing ingredients for future simulations of high-energy particle collisions. The method can easily be adapted to and implemented in other event generators. To illustrate its potential, we here consider combined tunes of the multi-parton-interaction and colour-reconnection models in Sherpa using LHC proton-proton collision data at $\sqrt{s}=7\,\text{TeV}$. We furthermore calibrate the energy-scaling behaviour of dimensionful model parameters based on $\sqrt{s}=13\,\text{TeV}$ LHC data and Tevatron data taken at $\sqrt{s}=1.96\,\text{TeV}$.

Pabst, Moritz [Gottingen U.]↗
DOE OSTI2018

Linking big models to big data: efficient ecosystem model calibration through Bayesian model emulation

Data-model integration plays a critical role in assessing and improving our capacity to predict ecosystem dynamics. Similarly, the ability to attach quantitative statements of uncertainty around model forecasts is crucial for model assessment and interpretation and for setting field research priorities. Bayesian methods provide a rigorous data assimilation framework for these applications, especially for problems with multiple data constraints. However, the Markov chain Monte Carlo (MCMC) techniques underlying most Bayesian calibration can be prohibitive for computationally demanding models and large datasets. We employ an alternative method, Bayesian model emulation of sufficient statistics, that can approximate the full joint posterior density, is more amenable to parallelization, and provides an estimate of parameter sensitivity. Analysis involved informative priors constructed from a meta-analysis of the primary literature and specification of both model and data uncertainties, and it introduced novel approaches to autocorrelation corrections on multiple data streams and emulating the sufficient statistics surface. We report the integration of this method within an ecological workflow management software, Predictive Ecosystem Analyzer (PEcAn), and its application and validation with two process-based terrestrial ecosystem models: SIPNET and ED2. In a test against a synthetic dataset, the emulator was able to retrieve the true parameter values. A comparison of the emulator approach to standard brute-force MCMC involving multiple data constraints showed that the emulator method was able to constrain the faster and simpler SIPNET model's parameters with comparable performance to the brute-force approach but reduced computation time by more than 2 orders of magnitude. The emulator was then applied to calibration of the ED2 model, whose complexity precludes standard (brute-force) Bayesian data assimilation techniques. Both models are constrained after assimilation of the observational data with the emulator method, reducing the uncertainty around their predictions. Performance metrics showed increased agreement between model predictions and data. Our study furthers efforts toward reducing model uncertainties, showing that the emulator method makes it possible to efficiently calibrate complex models.

58 GEOSCIENCES↗
DOE OSTI2021

Sensitivity of a Subject-specific Ankle Sprain Simulation to Extrinsic Versus Intrinsic Biomechanical Factors

Ankle sprains are the most common musculoskeletal injury in sport and military activity, despite existing prophylactic strategies. The purpose of this report was to develop a probabilistic simulation of lateral ankle sprains during single-limb drop landing, towards accelerating innovation in ankle sprain prevention. A deterministic, subject-specific musculoskeletal model was extended with automation and probabilistic distributions on sprain-related biomechanical factors. Probabilistic simulations were generated using traditional Monte Carlo techniques and the advanced mean value method, a more computationally-efficient approach. Predicted distributions of peak ankle joint rotations, velocities, and moments borne by supporting passive structures agreed favorably with the deterministic model and with reports of real sprain biomechanics. Parameter sensitivities identified that predictions were most strongly influenced by drop height, subtalar joint posture at contact, invertor/evertor co-activation, and passive ankle stiffness. The advanced mean value method predicted confidence bounds comparable to a 1000-trial Monte Carlo simulation, and required only 14 model evaluations and 4-min processing time. The extended probabilistic simulation may be useful to virtually test new prophylactic strategies for ankle sprains, and is made available for open-source use (https://simtk.org/projects/sprain-sim).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Entanglement Requirements for Coherent Enhancement in Detectors

Coherent enhancement is a powerful mechanism for improving the sensitivity of a wide range of detectors, but its practical use is often limited by the difficulty of preparing the required quantum states. We show that this difficulty has a fundamental origin: coherent enhancement of a signal interacting with a detector is quantitatively constrained by entanglement. We prove general bounds on how the strength of coherent effects can scale with system size, as a function of the single-mode entanglement entropy of the detector. These bounds smoothly interpolate between the incoherent and fully coherent regimes, and apply both to parameter-estimation problems and to scattering processes. We discuss these results from two complementary perspectives: First, they appear as bounds on the quantum Fisher information of many-body states, which translate directly into limits on parameter sensitivity via the quantum Cramér-Rao bound. Second, they can be interpreted as limits on a class of scattering cross sections, leading to predictions for how minimum detectable interaction strengths scale with target size. Together, these results provide a unified view of coherent enhancement in metrology and scattering experiments, and motivate the development of new techniques for generating entangled detector states.

Bogorad, Zachary [Fermilab] (ORCID:000000019913647↗
DOE OSTI2022

Patient-Specific Mathematical Model of the Clear Cell Renal Cell Carcinoma Microenvironment

The interactions between cells and molecules in the tumor microenvironment can give insight into the initiation and progression of tumors and their optimal treatment options. In this paper, we developed an ordinary differential equation (ODE) mathematical model of the interaction network of key players in the clear cell renal cell carcinoma (ccRCC) microenvironment. We then performed a global gradient-based sensitivity analysis to investigate the effects of the most sensitive parameters of the model on the number of cancer cells. The results indicate that parameters related to IL-6 have high a impact on cancer cell growth, such that decreasing the level of IL-6 can remarkably slow the tumor’s growth.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells’ location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells’ diffusion and immune cells’ influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes’ diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

An examination of anticipated g-jitter on Space Station and its effects on materials processes

This study is concerned with the effects of g-jitter, the residual acceleration aboard spacecraft, on selected classes of materials processes. In particular, the anticipated acceleration environment aboard Space Station Freedom (SSF) and its potential effects are analyzed, but the topic is covered with a sufficient level of generality as to apply to other processes and to other vehicles as well. Some of the key findings of this study include: The present acceleration specifications for SSF are inadequate to assure a quality level low-g environment. The local g vector orientation is an extremely sensitive parameter for certain key processes, but can not be controlled to within the desired tolerance. Therefore, less emphasis should be placed upon achieving a tight control of SSF attitude, but more emphasis should be focused on reducing the overall level of the g-jitter magnitude. Melt-based crystal growth may not be successfully processed in the relatively noisy environment of a large inhabited space structure. Growth from vapor or from solution appears more favorable. A smaller space structure and/or a free flyer can provide better alternatives in terms of g-jitter considerations. A high priority (including budgetary) should be given to coordinated efforts among researchers, SSF designers, and equipment contractors, to develop practical experiment-specific sensitivity requirements. Combined focused numerical simulations and experiments with well-resolved acceleration measurements should be vigorously pursued for developing reliable experiment-specific sensitivity data. Appendices provide an extensive cross-referenced bibliography, a discussion of the merits offered by g-jitter analysis techniques, as well as definitions of relevant nondimensional quantities and a brief description of available accelerometry hardware.

Nelson, Emily S.↗
NASA NTRS2000

A Method for Assessing the Quality of Model-Based Estimates of Ground Temperature and Atmospheric Moisture Using Satellite Data

A method is developed for validating model-based estimates of atmospheric moisture and ground temperature using satellite data. The approach relates errors in estimates of clear-sky longwave fluxes at the top of the Earth-atmosphere system to errors in geophysical parameters. The fluxes include clear-sky outgoing longwave radiation (CLR) and radiative flux in the window region between 8 and 12 microns (RadWn). The approach capitalizes on the availability of satellite estimates of CLR and RadWn and other auxiliary satellite data, and multiple global four-dimensional data assimilation (4-DDA) products. The basic methodology employs off-line forward radiative transfer calculations to generate synthetic clear-sky longwave fluxes from two different 4-DDA data sets. Simple linear regression is used to relate the clear-sky longwave flux discrepancies to discrepancies in ground temperature ((delta)T(sub g)) and broad-layer integrated atmospheric precipitable water ((delta)pw). The slopes of the regression lines define sensitivity parameters which can be exploited to help interpret mismatches between satellite observations and model-based estimates of clear-sky longwave fluxes. For illustration we analyze the discrepancies in the clear-sky longwave fluxes between an early implementation of the Goddard Earth Observing System Data Assimilation System (GEOS2) and a recent operational version of the European Centre for Medium-Range Weather Forecasts data assimilation system. The analysis of the synthetic clear-sky flux data shows that simple linear regression employing (delta)T(sub g)) and broad layer (delta)pw provides a good approximation to the full radiative transfer calculations, typically explaining more thin 90% of the 6 hourly variance in the flux differences. These simple regression relations can be inverted to "retrieve" the errors in the geophysical parameters, Uncertainties (normalized by standard deviation) in the monthly mean retrieved parameters range from 7% for (delta)T(sub g) to approx. 20% for the lower tropospheric moisture between 500 hPa and surface. The regression relationships developed from the synthetic flux data, together with CLR and RadWn observed with the Clouds and Earth Radiant Energy System instrument, ire used to assess the quality of the GEOS2 T(sub g) and pw. Results showed that the GEOS2 T(sub g) is too cold over land, and pw in upper layers is too high over the tropical oceans and too low in the lower atmosphere.

Wu, Man Li C.↗
DOE OSTI2024

A macro-micro approach for identifying crystal plasticity parameters for necking and failure in nickel-based alloy haynes 282

Here, this work develops a two-scales macro-micro approach to address the challenge in calibrating crystal plasticity microstructural models when samples undergo necking prior to fracture. The crystal plasticity models are crucial for predicting the materials’ plastic deformation and failure at the microstructure level, identifying the materials’ intrinsic properties as well as investigating the microstructure-properties relationships. However, after necking occurs, the experimentally measured stress-strain curves fail to reflect the materials ‘true’ stress-strain behavior and cannot be directly fitted into crystal plasticity models. The proposed macro-micro approach employs a top-down strategy to address this challenge, which has been studied with experimental tests on precipitation-strengthened Ni-based superalloy Haynes® 282®. In this approach, a macro rate-dependent anisotropic plasticity model with Voce-type hardening and Rice-Tracey damage law is first utilized to model the deformation and failure of the tensile bar, and calibrated by matching the stress-strain curves, necking strain, and reduction of area. Especially, to match the testing results under different applied strain rates, the rate-sensitivity parameter m and saturation stress in the elasticity model are modified to incorporate dependence on the local strain rate. Then, the ‘true’ stress-strain behaviors are extracted from the necking zone of the macro-model, which are used to calibrate a micro-model with explicit microstructures and governed by an extended crystal plasticity law. The consistency between the micro-model and macro-model are enforced during calibration. The calibration outcomes from the crystal plasticity model elucidate the materials intrinsic properties for slip, hardening, and failure, which is vital for further investigations into the microstructure-properties relationship and for accurate prediction of the material behavior under various test and service conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantitative study of the effect of grain boundary parameters on the slip system level Hall-Petch slope for basal slip system in Mg-4Al

Several theoretical studies have reported that the geometry and structure of grain boundaries in polycrystalline materials could impose a significant effect on the Hall-Petch slope. However, experimental observations are primarily limited by the ability of the techniques to accurately quantify the grain boundary strength and validate these theoretical models. Using high-resolution electron backscatter diffraction (HR-EBSD), the local stress tensor ahead of a slip band blocked by a grain boundary was quantified and coupled with a continuum dislocation pile-up model to assess the barrier strength of specific grain boundaries to specific slip systems, referred to as micro-Hall-Petch coefficient. For basal slip system in a deformed Mg-4Al alloy, the micro-Hall-Petch coefficient ($k^{basal}_{μ}$) varied significantly, from 0.054 to 0.184 MPa - m 1/2 for nine different grain boundaries. These results were correlated with geometric descriptors of the respective grain boundaries, with three-dimensional GB profile additionally measured via focused ion beam milling. It was found that the angle between the two slip plane traces on the grain boundary plane was the most sensitive parameter affecting $k^{basal}_{μ}$, followed by the angle between the slip directions. A functional form for calculation of $k^{basal}_{μ}$ depending on these two angles is proposed to augment crystal plasticity constitutive models with slip resistance dependent on some measure of the grain size. The method allows a new pathway to calibrate grain size strengthening parameters in crystal plasticity models, allowing further computational investigations of the interrelationship between texture, grain morphology, and the Hall Petch effect.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An Integrative Model for Soil Biogeochemistry and Methane Processes: I. Model Structure and Sensitivity Analysis

Abstract Environmental changes are anticipated to generate substantial impacts on carbon cycling in peatlands, affecting terrestrial‐climate feedbacks. Understanding how peatland methane (CH 4 ) fluxes respond to these changing environments is critical for predicting the magnitude of feedbacks from peatlands to global climate change. To improve predictions of CH 4 fluxes in response to changes such as elevated atmospheric CO 2 concentrations and warming, it is essential for Earth system models to include increased realism to simulate CH 4 processes in a more mechanistic way. To address this need, we incorporated a new microbial‐functional group‐based CH 4 module into the Energy Exascale Earth System land model (ELM) and tested it with multiple observational data sets at an ombrotrophic peatland bog in northern Minnesota. The model is able to simulate observed land surface CH 4 fluxes and fundamental mechanisms contributing to these throughout the soil profile. The model reproduced the observed vertical distributions of dissolved organic carbon and acetate concentrations. The seasonality of acetoclastic and hydrogenotrophic methanogenesis—two key processes for CH 4 production—and CH 4 concentration along the soil profile were accurately simulated. Meanwhile, the model estimated that plant‐mediated transport, diffusion, and ebullition contributed to ∼23.5%, 15.0%, and 61.5% of CH 4 transport, respectively. A parameter sensitivity analysis showed that CH 4 substrate and CH 4 production were the most critical mechanisms regulating temporal patterns of surface CH 4 fluxes both under ambient conditions and warming treatments. This knowledge will be used to improve Earth system model predictions of these high‐carbon ecosystems from plot to regional scales.

58 GEOSCIENCES↗
DOE OSTI2019

Evaluating the Importance of Barometric Pumping for Subsurface Gas Transport Near an Underground Nuclear Test Site

Core Ideas Gas transport modeling can estimate ground surface breakthrough after an underground nuclear explosion. Comparison tracer‐injection experimental and model results reveals parameter sensitivity. Barometric pumping may explain gas breakthrough following an underground nuclear explosion. Amplitude and period of barometric pressure signal are key in controlling breakthrough. An underground nuclear explosion (UNE) generates and distributes radioactive gases that can be transported to the ground surface though preexisting and explosion‐induced fractures over timescales of hours to months. If detected, the presence of short‐lived radionuclides in gas is evidence of a recent UNE. Numerical modeling can provide estimates of surface arrival times that can help inform gas monitoring strategies at suspected foreign test sites. Efforts are underway at historic US UNE sites to better understand subsurface gas‐transport processes following a UNE by geologic characterization of the near‐field damage structures, field‐scale tracer experiments, and subsurface air pressure monitoring. The development of numerical models using historical and experimental datasets from former UNE sites can improve predictions by testing conceptual models, highlighting key processes, and constraining parameter ranges. The models developed in this study represent the U20az site at the Nevada National Security Site where the Barnwell device was expended in December 1989. A two‐phase (water and air), dual‐permeability flow and transport model of the U20az site was built to investigate gas transport processes under recent experimental conditions and following the Barnwell nuclear event. Results indicate that the model can explain both the lack of arrival of radioactive gas tracers in a 2013 field experiment as well as the observed arrival of radioactive gases following the 1989 Barnwell event using barometric pressure records from the respective periods, even when additional advective processes associated with the Barnwell detonation are ignored. The results demonstrate that the character of the barometric records may be a key factor in explaining the differences in transport behavior.

58 GEOSCIENCES↗
DOE OSTI2021

Assessing the predictive impact of factor fixing with an adaptive uncertainty-based approach .

Despite widespread use of factor fixing in environmental modeling, its effect on model predictions has received little attention and is instead commonly presumed to be negligible. We propose a proof-of-concept adaptive method for systematically investigating the impact of factor fixing. The method uses Global Sensitivity Analysis methods to identify groups of sensitive parameters, then quantifies which groups can be safely fixed at nominal values without exceeding a maximum acceptable error, demonstrated using the 21-dimensional Sobol’ G-function. Furthermore, three error measures are considered for quantities of interest, namely Relative Mean Absolute Error, Pearson Product-Moment Correlation and Relative Variance. Results demonstrate that factor fixing may cause large errors in the model results unexpectedly, when preliminary analysis suggests otherwise, and that the default value selected affects the number of factors to fix. To improve the applicability and methodological development of factor fixing, a new research agenda encompassing five opportunities is discussed for further attention.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Thermal Transport and Fracture Behavior of Sintered Fuel Pellets: Experimental Validation of NEAMS Tool MARMOT

The project targets at experimental validation of NEAMS tool MARMOT for predicting thermal transport and fuel fracture. The key components of the project include: (1) sintering UO 2 samples with well-controlled microstructures (grain size, pore, fission bubble, porosity and pore distribution either in intra- or inter-granular pores); (2) thermal conductivity measurements of these wellcharacterized polycrystalline samples; (3) indentation testing of sintered UO 2 with various grain size, porosity and stoichiometry to obtain fracture and crack propagation mechanisms (inter or intragranular fracture); and (4) validation and uncertainty quantification of the MARMOT thermal transport and fracture models. The high quality experimental data will be used to validate the predictions of MARMOT for thermal transport and fracture, using average data and local microstructure information. In this work, UO 2 fuels with controlled microstructure (grain size, pore size/structure and fuel chemistry) were successfully prepared by high energy ball milling of the hyper-stoichiometric UO 2+x powders purchased commercially, and consolidated by spark plasma sintering. Dense nano-sized UO 2 fuel pellets with controlled stoichiometry were sintered by SPS. The effects of different SPS sintering parameters on the pore size and distribution were investigated, and the correlation among the SPS parameter – microstructure control – stoichiometry/defect structure was established. To further control the pore structure and fuel chemistry, systematic thermal annealing of the sintered fuels with the controlled pore structure was performed in a reduced. The controlled microstructure and fuel chemistry enable the investigation of the separate effects of pore structure and fuel chemistry on thermal transport behavior. The mechanical properties of the sintered UO 2 fuels were tested by nanoindentation and microindentation testing at different temperatures. The fracture behavior of the sintered fuels as functions of grain structure and temperature was investigated. Nanocrystalline UO 2 display higher hardness than microcrystalline counterpart, consistent with the Hall-Petch strengthening mechanism. Greater Young’s modulus and fracture toughness are also identified for the nanocrystalline UO 2 , and hardness and Young’s modulus decrease with temperature, suggesting better ductility of oxide fuels at greater temperature at small length scale. Hyper-stoichiometric UO 2 specimen displays higher hardness and fracture toughness than stoichiometric UO 2 , due to the impediment of the crack propagation by the oxygen interstitial atoms. These results are useful in understanding the mechanical properties of the high burn-up structure (HBS) formed in nuclear fuels during reactor operation, and also provide critical experimental data as the input for the development and validation of the MARMOT fracture model of nuclear fuels. The thermal properties of the sintered UO 2 fuel pellets with controlled grain structure, fuel chemistry and pore structures are measured by laser flash and the impact of the microstructure characteristics on the thermal transport properties is elucidated. Through a close collaboration with the MARMOT team at Idaho National Laboratory, the thermal properties of the sintered oxide fuel pellets with controlled microstructures are measured experimentally, and modeled using the developed MARMOT thermal transport models for a joint effort of demonstrating high spatial resolution measurements at the microstructural level on as-fabricated nuclear fuels and materials and obtain critical experimental data for model validation and properties predictions. In this work, we also conducted sensitivity analysis and uncertainty quantification (UQ) on a mesoscale simulation that has been applied to understand the effective thermal conductivity of UO 2 reactor fuel using the MARMOT tool. We also conduced UQ on mechanistic macroscale models that have been developed to predict the effective thermal conductivity. The most sensitive parameter is the thermal conductivity of bulk UO 2 for the thermal conductivity models. In the future, the predicted simulation distributions need to be compared to experimental data for validation of the models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Modelling the Sensitivity of Yukon River Biogeochemical Dynamics to Environmental and Chemical Drivers: Implications for Dissolved Organic Carbon

Riverine dissolved organic carbon (DOC) is a critical biogeochemical component that transmits information from Arctic soils to the Arctic Ocean, significantly influencing carbon dynamics in this unique ecosystem. As DOC travels downstream, it undergoes transformations that alter its composition and fate. The Yukon River serves as an effective testbed for modelling these dynamics, offering sufficient scale to capture key biogeochemical processes while having a simpler hydrology than other major Arctic rivers, as well as long-term DOC observational data for model validation. To investigate DOC transformations during transit in the Yukon River, we adapted our Arctic Riverine Organic Macromolecular Model by applying regional-specific parameterisations. Our model simulates the transport and transformation of 15 organic macromolecules, including CDOM (coloured dissolved organic matter), proteins, polysaccharides, lipids, lignin phenols, and humic substances. Initial DOC concentrations were derived from observed soil organic carbon stocks in the surrounding watershed, while chemical transformations and hydrological dynamics were modelled along the river's course. Sensitivity and uncertainty analyses were conducted using a Monte Carlo approach under two experimental setups. Results revealed that variability in DOC and CDOM concentrations at the river mouth were predominantly driven by initial DOC concentration (~70% of variability explained) and dilution at confluence points (~10%). The refractory fraction of DOC explained 21%–88% of the variability in 14 macromolecular concentrations and ranked in the top five sensitive parameters for all outputs when a uniform parameter distribution was assumed. However, when a more likely variability was applied to this parameter, its influence on DOC and CDOM decreased. Given that refractory DOC accounts for ~80% of total DOC in Arctic Rivers, this suggests that most DOC resists degradation and retains its chemical composition during transport to the coastal environment. River velocity, which determines residence time, explained 8%–47% of the variability in protein, polysaccharide, lipid, pigments, and lignin phenols at the river mouth. In contrast, chemical turnover times contributed only 1%–5% to output variability. Our findings underscore the need for improved land-specific headwater observations, including seasonal soil moisture and lateral transport dynamics that control the initial tributary-specific DOC inputs. With accelerated permafrost thaw and increasing river discharge, extending our model to other Arctic River systems and seasons will enhance understanding of Arctic riverine carbon fluxes and their contributions to the Arctic Ocean.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗