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 379 records · Page 21

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗

Searching for neutrino polarizability at DUNE

We perform an initial study of DUNE’s sensitivity to enhanced neutrino polarizability within models of light scalar mediators. We identify two possible signatures of polarizability due to neutrino scattering on electrons and due to coherent scattering on argon nuclei. These result in either one or two separated electromagnetic showers, respectively, with no associated hadronic activity. For each signature we compute the signal rates and the relevant backgrounds, obtaining the projected reach of the DUNE near detector. We then compare this with the current astrophysical and terrestrial bounds on light scalar models coupling to neutrinos and/or photons.

Carey, Sam [Wayne State U.] (ORCID:000900029607531↗

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗

Improving creep resistance of Al-0.27Zr-0.08Sn (wt%) via cold swaging while maintaining high electrical conductivity

This study investigates the effect of mechanical cold work (0, 50, and 90% cross-sectional area reduction via swaging) on the creep properties of cast aluminum alloys with high electrical conductivity including high-purity Al (HP-Al), Al-0.15Zr (Al-Zr), and Al-0.27Zr-0.08Sn (Al-Zr-Sn, wt%). The Al-Zr alloy is strengthened by Zr in solid solution while the Al-Zr-Sn alloy is additionally strengthened by Al3Zr nanoprecipitates (8 nm diameter). After 90% cold swaging, the alloys exhibit grains elongated along the swaging direction with their widths ranging between 8 and 152 μm perpendicular to the swaging direction. Creep testing at 200 °C shows that the minimum creep rate of all 90% swaged alloys shows high sensitivity to stress, which can be modeled via a threshold stress σth. Increasing swaging magnitude from 0 to 50 to 90% in Al-Zr-Sn improves the creep resistance without reducing electrical conductivity. The formation of subgrains, increased dislocation density, and columnar grain structure in swaged alloys all enhance creep resistance. HP-Al swaged to 90% exhibits much lower creep resistance (σth = 18 MPa) due to subgrain coarsening and the absence of precipitates compared to Al-Zr (σth = 40 MPa) and Al-Zr-Sn (σth = 40 MPa), which maintain stable subgrain structures and benefit from solid solution strengthening and Al3Zr nanoprecipitates, respectively. Although Al3Zr precipitates are known to stabilize the deformed microstructure by pinning grain-boundaries, this work demonstrates that Zr in solid solution alone can effectively stabilize the deformed microstructure resulting in enhanced creep resistance. This effect of Zr solute on stabilizing grain boundaries for creep performance was not previously well-defined in the literature. Despite similar creep performances of Al-Zr and Al-Zr-Sn, the latter shows 70% higher room-temperature microhardness and slightly improved electrical conductivity (56%IACS (International Annealed Copper Standard) vs. 57%IACS, respectively). This work provides new insights into creep mechanisms of cold-worked Al–Zr alloys that will guide the design of heat-resistant Al conductors for high-demand applications.

Coello ramirez, Ismael [Northwestern University,]↗

Conceptual design study of neutron detectors for safeguards measurement of an irradiated pebble

Nuclear material control and accounting (MC&A) of pebble-bed reactors (PBRs) is challenging because a PBR utilizes hundreds of thousands of identical, unmarked pebbles that are continuously recirculated through the core. To develop tools that enable the implementation of international safeguards, especially in the context of MC&A of spent pebbles, we designed and simulated three neutron detection concepts to determine fissile content in individual pebbles: a differential die-away (DDA) detector, a californium interrogation prompt neutron (CIPN) detector, and a passive neutron albedo reactivity (PNAR) detector using Monte Carlo calculations. Burnup calculations were performed on the spent pebbles from the PBMR-400 classic PBR. The varying neutron and gamma source terms, and isotopic compositions in the spent pebbles calculated at various burnup levels were used in the neutron detector models. DDA was found to be sensitive to the number of passes a pebble has had through the core and to the fissile content contained in a spent pebble. Optimization in the DDA design further increased the neutron count rates and thus reduced counting uncertainty. Meanwhile, passive neutron counting using the same detector body could distinguish pebbles with different numbers of passes, but its response was dominated by neutron-emitting actinides and was not sensitive to fissile content. On the other hand, the PNAR technique was not viable for a single pebble but performed reasonably for a 27-pebble array, which suggested potential use for verification measurements of containers filled with 27 or more spent pebbles.

CIPN↗

A Probabilistic Model for Global EMIC Wave Activity Using Van Allen Probes Observations

Electromagnetic ion cyclotron (EMIC) waves play a key role in radiation belt dynamics through resonant interactions. However, their low occurrence probability, high variability, and spatial intermittency pose challenges for accurate modeling. In this study, we present a machine learning (ML)-based global EMIC wave model built on the entire data set from the Van Allen Probes mission. To capture the distinct statistical characteristics of wave occurrence and amplitude, the model is separated into two modules: an occurrence model trained using ML techniques, and a wave amplitude model sampled from observed probability distributions. The input parameters are limited to real-time or predictable variables to ensure practical applicability. Our model shows strong performance across the entire test set and demonstrates improved predictive capability over a baseline random occurrence model, particularly during quiet geomagnetic conditions. Evaluation during both quiet and active periods confirms the model's ability to represent the clustered and intermittent nature of EMIC wave activity. Furthermore, the model provides global estimates of wave power, enabling integration with radiation belt electron data and showing signatures consistent with wave-induced scattering. We found a good correlation between the global wave activity from the model and relativistic electron observation by Van Allen Probes, regardless of the availability of in situ wave observations. The modular structure of the model also allows for straightforward expansion for additional wave properties, such as wave frequency, which can be modeled independently. This flexible, event-sensitive approach offers a promising framework for data-driven radiation belt simulations and space weather applications.

79 ASTRONOMY AND ASTROPHYSICS↗

Seasonal stabilization effects slowed the greening of the Northern Hemisphere over the last two decades

Rising atmospheric CO₂ and warming spring temperatures increase vegetation growth and the terrestrial carbon sink. However, drought, heat stress, phenology, and resource limitations may stabilize or limit theses projected increases. We investigate the balance between these amplifying and stabilizing ecological factors by asking whether enhanced early-season growth leads to continued late-season growth. Using the Moderate Resolution Imaging Spectroradiometer (MODIS) leaf area index (LAI) dataset, we identify three seasonal growth patterns based on early- and peak-season positive LAI anomalies: (1) amplification, where late-season LAI anomalies exceed earlier ones; (2) weak stabilization, where late-season anomalies remain similar or slightly lower; and (3) strong stabilization, where late-season anomalies become negative. Weak and strong stabilization events dominate across 67% and 26% of Northern Hemisphere ecosystems above 30°N, respectively. The absence of any trend in amplifying or stabilizing events suggests stabilizing factors seasonally offset CO₂ and temperature-induced spring greening. Terrestrial biosphere models underestimate strong stabilization and overestimate amplification events. This inconsistency arises from the models’ underestimation late-season LAI sensitivity to precipitation in water-limited regions; overlook negative legacy effects of early enhanced LAI on late-season soil moisture via evapotranspiration losses in energy-limited regions. Our findings suggest water/heat stress and resource limitations limit greening and the land carbon sink.

Zhang, Wen↗

Statistical inference of collision frequencies from x-ray Thomson scattering spectra

Thomson scattering spectra measure the response of plasma particles to incident radiation. In warm dense matter, which is opaque to visible light, x-ray Thomson scattering (XRTS) enables a detailed probe of the electron distribution and has been used as a diagnostic for electron temperature, density, and plasma ionization. In this work, we examine the sensitivities of inelastic XRTS signatures to modeling details, including the dynamic collision frequency and the electronic density of states. Applying verified Monte Carlo inversion methods to dynamic structure factors obtained from time-dependent density functional theory, we assess the utility of XRTS signals as a way to inform the dynamic collision frequency, especially its direct-current limit, which is directly related to the electrical conductivity.

Collision frequency↗

Combined search for electroweak production of winos, binos, higgsinos, and sleptons in proton-proton collisions at s = 13 TeV

A combination of the results of several searches for the electroweak production of the supersymmetric partners of standard model bosons, and of charged leptons, is presented. All searches use proton-proton collision data at s = 13 TeV recorded with the CMS detector at the LHC in 2016–2018. The analyzed data correspond to an integrated luminosity of up to 137 fb − 1 . The results are interpreted in terms of simplified models of supersymmetry. Two new interpretations are added with this combination: a model spectrum with the bino as the lightest supersymmetric particle together with mass-degenerate Higgsinos decaying to the bino and a standard model boson, and the compressed-spectrum region of a previously studied model of slepton pair production. Improved analysis techniques are employed to optimize sensitivity for the compressed spectra in the wino and slepton pair production models. The results are consistent with expectations from the standard model. The combination provides a more comprehensive coverage of the model parameter space than the individual searches, extending the exclusion by up to 125 GeV, and also targets some of the intermediate gaps in the mass coverage. © 2024 CERN, for the CMS Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

IAM-FIRE: a Climate Emulator–Based Framework to Project Wildfire Impacts and Risks for Integrated Assessment Models

Most Integrated Assessment Models (IAMs) underrepresent dynamic feedbacks from climate-driven disturbances such as wildfires, potentially overestimating the permanence of land-based carbon sinks. In particular, representing the impacts of forest fires is becoming increasingly important, as these are expected to intensify in the coming years. We introduce IAM-FIRE (Integrated Assessment Model – Fire Impacts & Risks Emulator), a novel framework that enables the projection of wildfire burned area (BA) and carbon emissions (CE) directly from IAM outputs. IAM-FIRE combines a spatial climate emulator, land-use downscaling, vegetation productivity modelling, and an empirical fire model to generate global annual wildfire impacts for arbitrary socioeconomic and emissions scenarios at 0.5° resolution for the period 2020–2100. Calibrated against GFEDv5 observations and using inputs from the Global Change Analysis Model (GCAM), we report projections BA and CE derived from IAM-FIRE for four scenarios: SSP1-2.6, SSP2-4.5, SSP3-6.6 and SSP5-7.6. The model reproduces historical global trends for total BA, including the observed global decline since the early 2000s, and for forest BA. Projected fire trajectories differ strongly among scenarios: total BA range from declines under SSP1-2.6 (-3.36 Mha yr-1) to increases under SSP3-6.6 (+1.6 Mha yr-1). Corresponding total CE show a similar divergence ranging from -15 to +10.6 TgC yr-1. Socioeconomic development exerts a dominant suppressing effect on wildfire impacts while climate change and CO2-driven increases in vegetation productivity amplify fire risk, particularly under high-emissions pathways. Compared with CMIP6 fire-enabled Earth System Models, IAM-FIRE exhibits greater sensitivity to radiative forcing and a stronger role for human-driven fire suppression, highlighting substantial structural uncertainties in future fire projections. By providing a computationally efficient and internally consistent approach to represent wildfire impacts within IAMs, IAM-FIRE enables systematic exploration of fire–climate–land feedbacks and supports improved assessments of mitigation permanence and climate risks in future integrated scenarios.

Rouhette, Theo↗

Model Scripts for "Old-Aged Groundwater Contributes to Mountain Hillslope Hydrologic Dynamics"

The partitioning of water inputs between deep and shallow groundwater flow paths is a fundamental processes, yet is challenging to observe. Numerical models provide a valuable tool to further develop insights on these groundwater mixing processes. This repository contains scripts to run the ParFlow-CLM and EcoSLIM integrated hydrologic models along the Pumphouse Hillslope in the East River Watershed and associated python scripts to process model outputs. Files includes input decks for the associated models, python scripts, and text and css files to support the model runs and processing. The 2-D hillslope model simulates from 2000 to 2021 using transient forcing conditions from the publicly available NLDAS-2 dataset. We develop an ensemble of models with variable hydrogeologic and soil parameters. Model outputs from the ensemble of runs are compared to to water level observations at the PLM1 and PLM6 wells published at (https://data.ess-dive.lbl.gov/view/doi:10.15485/1866836) and groundwater mean ages published at (https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1960042. The model results are used to evaluate hillslope water mass-balance transience and mixing dynamics between groundwater flow paths with young (<10 year) and old (>10 year) ages.

54 ENVIRONMENTAL SCIENCES↗

FY24: LEMMs: Long-term, Electrochemical Materials degradation Models

FY24 progress for LEMMs: Long-term, Electrochemical Materials degradation Models LDRD project are presented. Significant progress toward creating validated corrosion models was made in this FY paving success in future FY’s for success with battery model development and validation. Specifically, models for various forms of corrosion were probed for sensitivities showing a significant influence of reactive transport parameters for hydroxide species on the run time of the models. Further, cryo-genic time-of-flight secondary ion mass spectroscopy was utilized to map ions in a frozen corrosion droplet showing the precipitates/precipitate species near corroding locations. Additionally, a realistic corrosion droplet was created that accounted for evaporation and condensation coupled with corrosion. This pushes the boundary of corrosion modeling and will be validated in future FY’s.

36 MATERIALS SCIENCE↗

Once-Through Steam Generator Model Analysis Using Python and Advanced Optimization Tools (Summer Internship Report)

This study focuses on the parametric analysis of design parameters for a once-through steam generator (OTSG) model, using python and advanced optimization tools to facilitate applications such as the flowing autoclave steam generator (FASG) test cases. Building on previous research involving another OTSG with a different design, this project aims to enhance our understanding of how steam generators (SGs) behave and how their outputs are influenced by changes in design. The reason for this design change is to allow for more precise modeling and optimization of SG performance, to provide a comparative analysis between the two designs, and to set up the model for integration with the FASG test case. The OTSG python-model is a mathematical representation (including fluid flow and heat transfer equations/models/correlations) of a steam-generating unit in a pressurized water reactor-type small modular reactor system. Design studies involve changing the model’s input design parameters to observe the resulting effects on the output of the system. By using advanced optimization tools, such as the Risk Analysis Virtual Environment (RAVEN) developed at Idaho National Laboratory, detailed design parametric studies and model optimization were performed. Six input parameters—pressure, temperature, and mass flow rate for the inlet of the primary-side (hot fluid) and secondary-side (cold fluid), respectively, of the SG—were randomly perturbed via RAVEN’s Monte Carlo Sampler module, using uniform distributions (i.e., ±1%, ±5% and ±10% relative changes) for 600 samples. The analysis provides valuable insights into SG optimization and can be used for sensor placement optimization to effectively monitor and obtain experimental data in other tests.

20 FOSSIL-FUELED POWER PLANTS↗

Development of a Liquid Bridge Model for Particle Agglomeration and Defluidization in Plastic Pyrolysis

Molten plastic that forms during the pyrolysis of plastic or municipal solid waste feedstock can lead to particle agglomeration. A liquid bridge model is developed in MFiX-DEM and validated against a cold flow experiment with glass beads coated with silicone oil. The liquid bridge model is then extended to support an evolving liquid layer thickness for pyrolysis applications. The extended model is used to study the sensitivity of the pyrolysis reactor to solids holdup, flow conditions, and plastic properties.

Banerjee, Subhodeep↗

Fermilab complex upgrades and CLFV

I review the current status and future prospects of charged lepton flavor violation (CLFV) searches at Fermilab, with emphasis on their sensitivity to physics beyond the Standard Model and their complementarity to the mission-critical LBNF/DUNE neutrino program. In this context, strategic considerations for the evolution of the complex in the Linac-II era will be discussed, including synergies between CLFV initiatives and mission-critical commitments such as reliable high-power beam delivery to LBNF/DUNE.

Hedges, Michael [Fermilab] (ORCID:0000000165041872↗

Implementing belowground controls on nutrient uptake in ELMv2-SPRUCE improves representation of a boreal peatland ecosystem

Boreal peatlands store 13 %–32 % of the global soil carbon (C) stock, a service dependent on plant-mycorrhizal fungi associations. In these nutrient poor systems, ectomycorrhizal and ericoid mycorrhizal fungi supply up to >80 % of the nutrient requirements of their plant hosts, partly with mined nitrogen (N) and phosphorus (P) from soil organic matter that are otherwise inaccessible to plants. Despite the ecological significance, mycorrhizal associations are only represented in a few land surface or ecosystem models. We modify the peatland branch of version 2 of the Energy Exascale Earth System Land Model (ELMv2-SPRUCE) to replace the default photosynthesis-driven inorganic N and P (NP) uptake process with a more realistic representation of the process via three pathways: (1) direct inorganic NP uptake by uncolonized fine roots, (2) indirect inorganic NP acquisition and (3) indirect NP acquisition from organic sources by mycorrhizal roots. We systematically evaluated the performance of the default and modified models with field observations from a whole ecosystem warming and carbon dioxide fertilization experimental site: Spruce and Peatland Responses Under Changing Environment (SPRUCE), in northern Minnesota, USA. The modified model reduces the underestimation of the growth response of shrubs in the default model to warming from 40 %–80 % to 17 %–35 % and reduces the overall relative absolute error on C fluxes from 1.61 to 1.54 in calibration. Improvements on modeled shrub growths and shrub-moss community net ecosystem exchanges are also seen in validation. The improved growth response of shrubs to warming is accompanied by several-fold increase in direct inorganic NP uptake and decrease in fungal colonization rate. The modified model simulates a smaller magnitude of transition of the ecosystem from C sink to C source under warming due to alleviation of plant nutrient limitation. Equifinality analysis shows the newly added parameters in the modified model can be constrained by the observed C fluxes. Sensitivity analysis shows the newly added parameters have stronger statistical interactions than the preexisting parameters in the default model. Overall, the modified model is an improvement over the default ELMv2-SPRUCE and will be a useful tool for understanding boreal peatland change.

Wang, Yaoping [Oak Ridge National Laboratory (ORNL↗

Searching for Neutrino Tridents in the NOvA Near Detector

This dissertation presents a search for neutrino trident production in the NOvA near detector through the coherent ``dimuon" channel: $\nu_\mu +\hspace{1pt}\text{X} \rightarrow \nu_\mu + \mu^- + \mu^+ +\hspace{1pt}\text{X}$. Trident production is a rare, purely electroweak process with sensitivity to physics beyond the Standard Model. The theoretical background, motivation for studying the process, and previous experimental measurements are reviewed. The analysis uses data collected by the NOvA near detector (ND) from Fermilab's Neutrinos at the Main Injector (NuMI) beam between November 2014 and February 2024, corresponding to an exposure of $25.5\times 10^{20}$ protons on target. The ND is a segmented tracking calorimeter located 800~m from the beam target, receiving neutrinos with a mean energy of 2~GeV. A multi-pass background reduction strategy is implemented, including the development of a novel dimuon-specific tracking technique. Trident candidates are identified using a boost ed decision tree classifier trained on simulated signal and background events. Limited background Monte Carlo statistics necessitate the use of functional fits to sideband data, which are extrapolated to estimate backgrounds in the signal region. The unblinded data contain 9 trident-like events, with an estimated background of 5.66 $\pm$ 5.15 events. This yields a best fit estimate of 3.34 tridents compared to the Standard Model prediction of 4.66. A profiled Feldman-Cousins method is used to determine a 90\% confidence interval of [0,9.1] on the number of signal events, corresponding to an upper limit of 1.95$\times$ the Standard Model prediction. This result represents the lowest energy search for trident events to date, and the first experimental contribution to the process in 27 years.

Bowles, Reed Scott [Indiana U.]↗

Agentic Diagrammatica: Towards Autonomous Symbolic Computation in High Energy Physics

We present Diagrammatica, a symbolic computation extension to the HEPTAPOD agentic framework, which enables LLM agents to plan and execute multi-step theoretical calculations. Symbolic computation poses a distinctive reliability challenge for LLM agents, as correctness is governed by implicit mathematical conventions that are not encoded in a form that can be easily checked in the computational backend. We identify two complementary remedies, tool-constrained computation and targeted knowledge grounding, and pursue the first as the primary architecture. Concretely, we concentrate the agent's action distribution onto tool calls with convention-fixing semantics, in which the agent specifies a compact, human-auditable diagram specification and a trusted backend performs the symbolic or numerical manipulations exactly. The toolkit provides two complementary calculation paths consuming a shared diagram specification: Naive Dimensional Analysis (NDA) for order-of-magnitude rate estimates and Exact Diagrammatic Analysis (EDA) for tree-level symbolic calculations via automatic FeynCalc code generation, both supplemented by automatic Feynman diagram enumeration and a navigable theory knowledge base. The architecture is validated on two benchmarks: (1) an exhaustive catalog of all tree-level, single-vertex $1\to 2$ partial decay widths across scalar, fermion, and vector parents, with complete massless and threshold limits and Standard Model validation; and (2) an NDA sensitivity study of the muon decay multiplicity $μ^+ \to ν_μ\barν_e + n(e^+e^-) + e^-$, determining the maximum observable $n$ at current and planned muon experiments.

Menzo, Tony [Alabama U.; Fermilab] (ORCID:00000002↗