Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗
DOE OSTI2025

Corrigendum to “The Emissions Model Intercomparison Project (Emissions-MIP): quantifying model sensitivity to emission characteristics” published in Atmos. Chem. Phys., 23, 14779–14799, 2023

Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO 2 injection, seasonality of SO 2 and black carbon (BC) particulate emissions, and the assumed fraction of SO 2 emissions that is injected into the atmosphere as particulate phase sulfate (SO 4 ) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO 2 , SO 4 , and BC, with a particularly large relative variation for SO 2 , which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO 2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of -0.35?W?m-2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8?d), surface SO 2 concentration (up to 59?% decrease), and surface sulfate concentration (up to 23?% increase). Emitting SO 2 at height consistently increased SO 2 and SO 4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO 4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO 4 emission fraction are accurately and consistently represented in global models.

Ahsan, Hamza [Pacific Northwest National Laborator↗
DOE OSTI2026

Model Sensitivities of Biomass-Burning Aerosol Chemical Aging, Sulfate Formation, and Cloud Droplet Activation in the Southeastern Atlantic Using CESM and E3SM

Biomass-burning smoke drives large uncertainty in climate projections of the Earth's radiative balance. This is due to the chemical and physical evolution of smoke and its impact on clouds and radiation. Here we focus on the southeastern Atlantic region and its inflow of African biomass-burning smoke during August 2017. We evaluate smoke properties and processes in two coupled earth-system models, the Energy Exascale Earth System Model (E3SM) and Community Earth System Model (CESM). These are compared against in situ aircraft observations from two field campaigns, ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) and CLoud–Aerosol–Radiation Interaction and Forcing: Year 2017 (CLARIFY-2017). Observations reveal an increase and subsequent decrease in smoke mean diameter, and a steady decrease in the mass ratio of organic aerosol (OA) to black carbon aerosol (BC) (OA:BC) over 4–12 days of aging, neither captured by the base models. Implementation of a photolytic loss scheme for secondary organic aerosol (SOA)—as a proxy for other heterogeneous volatilization chemistry—and a ∼1-day conversion for primary OA to SOA significantly improves the representation of this loss. In the boundary layer, both models show dimethyl sulfide driving a large increase in the sulfate aerosol mass fraction from the free troposphere, which is consistent with observations. Finally, models tend to underpredict cloud droplet number concentration partially due to weak modeled turbulent updraft strength, and model performance improves when the parameterized turbulent updraft strength is increased substantially. These results are expected to provide insights into future model development to reduce climate model uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modeled sensitivity of multi-MA accelerator performance to electrode contaminant inventory

Significant particle-in-cell code development has enabled simulations of power flow in multi-MA accelerators to include the desorption of surface contaminants, their ionization into surface plasmas, and the impact of these plasmas on efficiency. The simulations base desorption on an Arrhenius equation, whose most significant unknown is the surface contaminant inventory. The sensitivity of power-flow simulations to this inventory is studied here using Sandia National Laboratories' Z accelerator with a 7-nH MagLIF load [Phys. Plasmas 17, 056303 (2010)]. Simulations are conducted in 3D cylindrical coordinates for the current-adder, or “convolute,” region of Z and in 2D for the final feed only. Simulated contaminant inventories are varied from 1 to 32 monolayers (MLs) in 2D, and 2 to 4 ML in 3D. The results reveal sensitivities to the local ratio of E/B⁠. The high B-field, low E-field region near the short-circuit load is insensitive to the contaminant inventory, where assumed values of 4–32 ML change the load current by ≤ 2%, and agree with experiment to within 2% at peak current. A 1-ML value is the outlier, increasing the load current by 5%, but still within measurement uncertainty. In contrast, the relatively higher E-field, lower B-field convolute region has slower contaminant desorption and higher-magnitude E-field penetration of the surface plasmas. The current loss in the convolute region does increase with contaminant inventory. The loss assuming 4 ML is 12% larger than for 2 ML, with 4 ML being the better match to experiment.

Arrhenius equation↗
DOE OSTI2024

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

FFTF HT9 Cladding Microstructure Characterization

The sodium-cooled fast reactor (SFR) is a promising candidate for next generation nuclear reactors, operating at extreme conditions which include high temperatures (>500?C core outlet temperature) and significant neutron damage. High-Cr martensitic HT9 steel is an excellent candidate for SFR cladding and duct material due to its compatibility with liquid sodium, good thermal conductivity, resistance to void swelling, and strong creep rupture strength [1-4].However, the harsh in-core environment of SFRs can cause complex microstructural changes and mechanical property degradation in HT-9. Ensuring the safe use of HT9 cladding for metallic fuel requires both a thorough understanding of its mechanical response to microstructure evolution as well as reliable microstructure-sensitive modeling predictions. Microstructure-sensitive modeling of high temperature creep behavior in HT9 cladding for SFR applications currently lack experimental data to model the phenomena accurately. To fill this need, methods to perform microstructural characterization have been developed and performed on HT9.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Toward real-time optimization through model reduction and model discrepancy sensitivities

Optimization problems arise in a range of scenarios, from optimal control to model parameter estimation. In many applications, such as the development of digital twins, it is essential to solve these optimization problems within wall-clock-time limitations. However, this is often unattainable for complex systems, such as those modeled by nonlinear partial differential equations. One strategy for mitigating this issue is to construct a reduced-order model (ROM) that enables more rapid optimization. In particular, the use of nonintrusive ROMs—those that do not require access to the full-order model at evaluation time—is popular because they facilitate the computation of optimization solutions within the wall-clock time requirements. However, the optimization solution will be unreliable if the iterates move outside the ROM training data. This article proposes the use of hyper-differential sensitivity analysis with respect to model discrepancy (HDSA-MD) as a computationally efficient tool to augment ROM-constrained optimization and improve its reliability. The proposed approach consists of two phases: (i) an offline phase where several full-order model evaluations are computed to train the ROM, and (ii) an online phase where a ROM-constrained optimization problem is solved, a limited number of full-order model evaluations are computed, and HDSA-MD is used to enhance the optimization solution. Numerical results are demonstrated for two examples, atmospheric contaminant control and wildfire ignition location estimation, in which a ROM is trained offline using inaccurate atmospheric data. In conclusion, the HDSA-MD update yields a significant improvement in the ROM-constrained optimization solution using only one full-order model evaluation online with corrected atmospheric data.

PDE-constrained optimization↗
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↗
DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

PSA 2025 Presentation: "Modeling and Sensitivity Analysis of a Generation IV Pebble Bed Reactor Using MELCOR 2.2"

Accompanying the advancement of reactor technologies is the need for computational modeling and simulation to predict their behavior under normal operating conditions and accident scenarios. New Generation IV reactor designs which employ non-conventional fuel have a particular need for modeling the behavior and release of radionuclides and other material from the fuel. In this work, MELCOR version 2.2, a system-level safety and accident scenario code developed by Sandia National Laboratories, was used to model a 200-MWth pebble bed modular reactor and calculate the inventories of circulating and deposited graphite, metal dust, and elemental components released from the fuel elements. A base case modeling the reactor under standard operating conditions was calculated using MELCOR and the inventories were extrapolated to 30 years of operation time using a logarithmic regression fit. A sensitivity analysis was also performed in which several key parameters for the base case model were modified to explore the effect of these changes on the inventories calculated by MELCOR. A set of transient scenario simulations for a depressurized loss of forced cooling (DLOFC) accident were also performed. The results of the sensitivity analysis and transient simulations are reported and discussed in relation to the modeling techniques used for this study.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Automated ICRF heating surrogate modeling via machine learning

This work introduces automated machine learning workflows that address critical bottlenecks in surrogate model development for Ion Cyclotron Range of Frequencies (ICRF) heating applications. The automated framework includes data analysis tools that transform raw datasets into actionable insights in seconds, replacing weeks of manual exploratory effort and ensuring consistent, reproducible dataset characterization. By integrating advanced hyperparameter optimization (HPO) methods including Bayesian optimization via BoTorch and Tree-structured Parzen Estimators (TPE), the framework significantly reduces model development time from weeks to hours, decreasing computational cost and required expertise, while enabling high-accuracy surrogate models. Compared to traditional hyperparameter scanning (HPS) techniques such as methodical, randomized, and grid searches, HPO methods achieve superior convergence and predictive performance, even when compared to already well-tuned reference models. On NSTX High Harmonic Fast Wave (HHFW) heating datasets, both Random Forest Regressor (RFR) and neural network surrogates demonstrate improved accuracy, achieving R 2 values beyond 0.97 and 0.98, respectively. The results show that while HPO gains are modest for robust architectures like RFR, they become essential for more sensitive models such as neural networks, highlighting the trade-offs across optimization strategies. Through automated workflows that eliminate manual hyperparameter tuning and require minimal ML expertise, this work enables widespread adoption of high-fidelity surrogate models across the fusion community for real-time plasma control, uncertainty quantification, rapid experimental scenario development, and integrated system optimization.

Sanchez-Villar, Alvaro [Princeton Plasma Physics L↗
DOE OSTI2025

Atmospheric Feedbacks Reverse the Sensitivity of Modeled Photosynthesis to Stomatal Function

Stomata mediate fluxes of carbon and water between terrestrial plants and the atmosphere. These fluxes are governed by stomatal function and can be modulated in many Earth system models by an empirical parameter within the calculation of stomatal conductance, the stomatal slope (𝑔 1⁢𝑀 ). Intuitively, 𝑔 1⁢𝑀 represents the marginal water cost of carbon, relating it to the emergent plant property of water use efficiency. Observations show that 𝑔 1⁢𝑀 can range widely across and within plant types in varying environments, and this distribution of 𝑔 1⁢𝑀 is not captured within Earth system models which represent each plant type with a single 𝑔 1⁢𝑀 value. Here we examine how 𝑔 1⁢𝑀 influences photosynthesis using coupled Earth system model simulations by perturbing 𝑔 1⁢𝑀 to observed 5⁢t⁢h and 95⁢t⁢h percentiles for each plant type. We find that high 𝑔 1⁢𝑀 reduces photosynthesis nearly everywhere, while low 𝑔 1⁢𝑀 has regionally dependent responses. Under fixed atmospheric conditions, low 𝑔 1⁢𝑀 increases photosynthesis in the Amazon and central North America but decreases photosynthesis in boreal Canada. These responses reverse when the atmosphere responds interactively due to spatially differing sensitivity to increases in temperature and vapor pressure deficit. Choice of 𝑔 1⁢𝑀 also influences photosynthetic response to changes in atmospheric carbon dioxide (CO 2 ), with lower and higher 𝑔1⁢𝑀 modifying total global response to elevated 2x preindustrial CO 2 by +6.4% and −9.6%, respectively. Our work demonstrates that atmospheric feedbacks are critical for determining the photosynthetic response to 𝑔 1⁢𝑀 assumptions and some regions are particularly sensitive to choice of 𝑔 1⁢𝑀 .

Liu, Amy X. [University of Washington, Seattle, WA↗
DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗
DOE OSTI2026

Informing Robust Functional Relationship Benchmarks: An Evaluation of the Temperature Sensitivity of Ecosystem Respiration Across the Arctic-Boreal Region

During land model development, simulated carbon dynamics are often benchmarked against observational data sets to evaluate model performance. Functional relationship benchmarks are the relationship between a driving variable (e.g., temperature) and a response variable (e.g., ecosystem respiration) and are a promising tool for assessing model performance by evaluating modeled sensitivities to changing environmental conditions. However, observed functional relationships can be influenced by choices made during data collection and throughout the benchmarking process, impacting the inferred skill of land models. To avoid misrepresenting a model's true performance, it is necessary to systematically evaluate best practices when constructing functional relationship benchmarks. We developed a set of guidelines for constructing functional relationship benchmarks, considering the choice of data set, number of daily observations, temporal extent, and temporal resolution across Alaska and Canada over a 20-year period from 2001 to 2020. The temperature sensitivity of ecosystem respiration from observations, evaluated through an apparent Q 10 , is highly variable both spatially and as a result of the data processing approach applied in the benchmark formation. When benchmarking 13 models from the Warming Permafrost Model Intercomparison Project (WrPMIP), the range in inferred model skill is substantially impacted by the choices applied in constructing functional relationship benchmarks. The inferred performance of a given model is most sensitive to the number of daily observations and temporal extent, followed by choice of benchmark data set and temporal averaging. Results from this analysis can guide the development of consistent and robust functional relationships for future model evaluation studies.

Poe, Jeralyn [Northern Arizona University, Flagsta↗
DOE OSTI2025

Using ARM Observations to Evaluate Process-Interactions in MCS Simulations Across Scales (Final Progress Report)

This project, funded by DOE Atmospheric System Research (DE-SC0020050), focused on improving the representation of mesoscale convective systems (MCSs) in numerical weather and climate models by leveraging high-resolution observations from the DOE Atmospheric Radiation Measurement (ARM) program. The research aimed to evaluate model sensitivities to grid spacing, microphysics, and planetary boundary layer (PBL) schemes, with a particular emphasis on improving convection parameterization for high-resolution modeling. Findings from this work highlight several key advancements. Model validation against ARM radar wind profiler data from the Southern Great Plains (SGP) and Manaus (MAO) sites revealed systematic biases in simulated convective mass flux profiles, leading to the development of an observationally constrained evaluation framework for diagnosing and improving model performance. Sensitivity analyses demonstrated that the representation of Amazonian MCSs was highly dependent on PBL scheme selection, while mid-latitude MCSs were more strongly influenced by microphysics parameterizations. A series of high-resolution WRF simulations, ranging from 4 km to 125 m grid spacing, provided insight into the behavior of convective drafts across scales. While updraft properties converged at sub-kilometer resolutions, biases in downdraft intensity persisted even at the finest resolution tested, emphasizing the need for further refinements in model physics. Additionally, comparisons of MCS vertical structures between mid-latitude and tropical environments revealed stronger updrafts and larger mass flux in mid-latitude MCSs, providing critical insights for improving climate model representations of storm-scale dynamics. The project’s findings have already contributed to advancing numerical modeling capabilities, particularly in WRF, MPAS, ICON, and DOE’s SCREAM model, by refining how convective processes are represented in high-resolution climate simulations. Results were disseminated through peer-reviewed publications, conference presentations, and ARM/ASR Research Highlights, engaging the broader scientific community. The project also provided valuable training opportunities for two postdoctoral researchers, who played central roles in model development, analysis, and dissemination of results. Their work contributed to several publications and conference presentations, helping prepare them for careers in atmospheric modeling. By improving the simulation of MCSs, this research directly supports the development of next-generation climate models capable of more accurately representing extreme precipitation and convective processes. The insights gained will inform future improvements in convective parameterization and guide the design of high-resolution weather and climate simulations, ultimately enhancing the reliability of climate projections and weather forecasts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Validation and moisture content sensitivity analysis of cross-laminated timber wall assemblies in EnergyPlus

Cross-laminated timber buildings are becoming more common in North America, with many numerical studies showing potential energy savings. However, no studies have validated any EnergyPlus heat transfer algorithms or quantified their accuracy in simulating CLT in building envelopes. This study empirically validates the heat flux predictions for each of EnergyPlus's heat transfer algorithms (Conduction Transfer Functions (CTF), Effective Moisture Penetration Depth (EMPD), Conduction Finite Difference (CondFD), and Heat and Moisture Transfer (HAMT)) for two different CLT ply thicknesses with both summer and winter boundary conditions measured in controlled lab experiments. It also evaluates the model sensitivity of heat flux and heating and cooling loads to moisture content. The 1D validation shows that the HAMT model is the most accurate among all algorithms. All EnergyPlus's heat flux predictions are accurate independent of CLT plate thickness for summer conditions. However, the three constant property algorithms (CTF, EMPD, and CondFD) underpredict heat flux throughout the whole day during winter conditions. The 1D sensitivity analysis indicates that elevated moisture content can increase peak heat fluxes through the material by up to 20 %. Finally, the whole building model sensitivity analysis shows increased heating load and slight cooling load variation due to increased moisture content when using constant property models. The analysis shows significantly lower peak thermal demand (7 % lower heating and 6 % lower cooling) and monthly thermal load (8 % less cooling and 6 % less heating) predictions when using HAMT vs a constant property model.

42 ENGINEERING↗