SEARCH · Engineering Papers

Results for “forward modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2021

The Broadband Counterpart of the Short GRB 200522A at z = 0.5536: A Luminous Kilonova or a Collimated Outflow with a Reverse Shock?

In this work, we present the discovery of the radio afterglow and near-infrared (NIR) counterpart of the Swift short gamma-ray burst (GRB) GRB 200522A, located at a small projected offset of ≈1 kpc from the center of a young, star-forming host galaxy at z = 0.5536. The radio and X-ray luminosities of the afterglow are consistent with those of on-axis cosmological short GRBs. The NIR counterpart, revealed by our Hubble Space Telescope observations at a rest-frame time of ≈2.3 days, has a luminosity of ≈(1.3–1.7) × 10 42 erg s -1 . This is substantially lower than on-axis short GRB afterglow detections but is a factor of ≈8–17 more luminous than the kilonova of GW170817 and significantly more luminous than any kilonova candidate for which comparable observations exist. The combination of the counterpart's color (i - y = -0.08 ± 0.21; rest frame) and luminosity cannot be explained by standard radioactive heating alone. Additionally, we present two scenarios to interpret the broadband behavior of GRB 200522A: a synchrotron forward shock with a luminous kilonova (potentially boosted by magnetar energy deposition), or forward and reverse shocks from a ≈14°, relativistic (Γ0 ≳ 10) jet. Models that include a combination of enhanced radioactive heating rates, low-lanthanide mass fractions, or additional sources of heating from late-time central engine activity may provide viable alternate explanations. If a stable magnetar was indeed produced in GRB 200522A, we predict that late-time radio emission will be detectable starting ≈0.3–6 yr after the burst for a deposited energy of ≈10 53 erg. Counterparts of similar luminosity to GRB 200522A associated with gravitational wave events will be detectable with current optical searches to ≈250 Mpc.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Aligning Standards Communities for Omics Biodiversity Data: Sustainable Darwin Core-MIxS Interoperability

The standardization of data, encompassing both primary and contextual information (metadata), plays a pivotal role in facilitating data (re-)use, integration, and knowledge generation. However, the biodiversity and omics communities, converging on omics biodiversity data, have historically developed and adopted their own distinct standards, hindering effective (meta)data integration and collaboration. In response to this challenge, the Task Group (TG) for Sustainable DwC-MIxS Interoperability was established. Convening experts from the Biodiversity Information Standards (TDWG) and the Genomic Standards Consortium (GSC) alongside external stakeholders, the TG aimed to promote sustainable interoperability between the Minimum Information about any (x) Sequence (MIxS) and Darwin Core (DwC) specifications. To achieve this goal, the TG utilized the Simple Standard for Sharing Ontology Mappings (SSSOM) to create a comprehensive mapping of DwC keys to MIxS keys. This mapping, combined with the development of the MIxS-DwC extension, enables the incorporation of MIxS core terms into DwC-compliant metadata records, facilitating seamless data exchange between MIxS and DwC user communities. Through the implementation of this translation layer, data produced in either MIxS- or DwC-compliant formats can now be efficiently brokered, breaking down silos and fostering closer collaboration between the biodiversity and omics communities. To ensure its sustainability and lasting impact, TDWG and GSC have both signed a Memorandum of Understanding (MoU) on creating a continuous model to synchronize their standards. These achievements mark a significant step forward in enhancing data sharing and utilization across domains, thereby unlocking new opportunities for scientific discovery and advancement.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Causal explicit algorithm for heat conduction in a plasma

Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Accelerating Nuclear Fuels and Materials Qualification by Multi-Level Irradiation Experiment Campaign

The advanced reactor technologies feature fuels, coolants, and materials that promise safer operating conditions under normal and accident scenarios. However, the nuclear fuels and materials qualifications require several decades (for example, new reactor fuel qualification from conceptualization requires about 20 years). Therefore, accelerating nuclear fuels and materials qualification is essential, and it can be achieved by combining high through-put materials irradiation and testing, advanced post-irradiation examinations, and Multiphysics modeling. This paper addressed the associated challenges in accelerating nuclear fuels and material qualification for new and advanced reactor designs, which differ based on fuel, coolant, operating conditions, and structural materials. These challenges vary for radiation level, operating conditions (e.g., temperature and pressure), and coolant type (e.g., corrosion environment). In addition, the challenges and limitations in modeling tools, experimental facilities, and licensing guidelines are also discussed, and a general solution path forward is recommended.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2022

Investigation of the background in coherent $J/\psi$ production at the EIC

Understanding fundamental properties of nucleons and nuclei are among the most important scientific goals of the next-generation machine, the Electron-Ion Collider (EIC). With the unprecedented versatility provided by the EIC, it will provide answers to many standing puzzles and open questions in modern nuclear physics. One of the golden measurements at the EIC is coherent vector meson production in electron-nucleus (eA) scattering in order to obtain the spatial gluon density distribution in heavy nuclei. This requires the experiment to overcome an overwhelmingly large background arising from the incoherent diffractive production, where the nucleus mostly breaks up into fragments of particles in the far-forward direction close to the hadron beam rapidity. In this report, we systematically study the rejection of incoherent J/\psi J / ψ production by vetoing products from the nuclear breakup - protons, neutrons, and photons, which is modeled with the BeAGLE event generator and the most up-to-date EIC Far-forward Interaction Region design.

Chang, Wan↗

DOE OSTI2025

Terra-Populus v0.1: A Python Library for LandScan High-Definition Population Analysis and Modeling

The terra-populus library is designed for use by the LandScan HD technical team, offering a streamlined set of tools for generating and updating LandScan HD datasets from foundational building-level data, referred to as 'molecules,' provided by the building-level attribution team. This document serves as the primary technical documentation for terra-populus. Version 0.1 of the library includes the core modeling components necessary for LandScan HD production. It enables the generation of the LandScan HD Baseline dataset as well as corresponding confidence measures for the occupancy rates used. Parameters have been included for incorporating damaged building indicators and changes in population, to faciliate the creation of rapid updates for LandScan HD. Future iterations of terra-populus will introduce tools for creating a confidence index, and quantifying and propagating uncertainty, facilitating the creation of probabilistic LandScan HD outputs. This report provides an overview of the tools available in the library and the corresponding code implementations. One of the key advancements implemented in terra-populus is a redefinition of the atomic modeling unit for LandScan HD. Traditionally, the LandScan HD vector analytical framework has generated population estimates at the building sub-component (molecule) level. However, terra-populus adopts a building-level modeling approach. This shift is an operational decision aimed at aligning LandScan HD outputs with confidence measures, which are computed and validated at the building level (confidence measures are not included in this version of terra-populus, aside from those associated with the occupancy rates). Additional advancements to the LandScan HD modeling, as implemented by terra-populus, include a minimum population value parameter and an auto assignment of building floor counts. The population minimum value was implemented to prevent buildings and subsequent LandScan HD pixels that contained small values that may not rasterize in production. An 'auto' value has been included as a method for dealing with buildings lacking floor count information, where it is the average floor count of all other buildings with a residential building use type tag. The logic behind this is to remain consistent with the current logic employed for dealing with building use type null instances, where a null use type is defaulted to residential since it is the most common building type. The auto logic is intended to apply the most common building floor count of the most common type of buildings. The tools provided in terra-populus represent a significant step forward in improving the efficiency, reproducibility, and transparency of the LandScan HD modeling process. As the library evolves, it will continue to serve as a foundational resource for high-resolution population modeling.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

DOE OSTI2023

Presentation: Accelerating Nuclear Fuels and Materials Qualification by Multi-Level Irradiation Experiment Campaign

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2024

CO2 response screen in grass Brachypodium reveals the key role of a MAP kinase in CO2-triggered stomatal closure

Abstract Plants respond to increased CO2 concentrations through stomatal closure, which can contribute to increased water use efficiency. Grasses display faster stomatal responses than eudicots due to dumbbell-shaped guard cells flanked by subsidiary cells working in opposition. However, forward genetic screening for stomatal CO2 signal transduction mutants in grasses has yet to be reported. The grass model Brachypodium distachyon is closely related to agronomically important cereal crops, sharing largely collinear genomes. To gain insights into CO2 control mechanisms of stomatal movements in grasses, we developed an unbiased forward genetic screen with an EMS-mutagenized B. distachyon M5 generation population using infrared imaging to identify plants with altered leaf temperatures at elevated CO2. Among isolated mutants, a “chill1” mutant exhibited cooler leaf temperatures than wild-type Bd21-3 parent control plants after exposure to increased CO2. chill1 plants showed strongly impaired high CO2-induced stomatal closure despite retaining a robust abscisic acid-induced stomatal closing response. Through bulked segregant whole-genome sequencing analyses followed by analyses of further backcrossed F4 generation plants and generation and characterization of sodium azide and CRISPR-cas9 mutants, chill1 was mapped to a protein kinase, Mitogen-Activated Protein Kinase 5 (BdMPK5). The chill1 mutation impaired BdMPK5 protein-mediated CO2/HCO3− sensing together with the High Temperature 1 (HT1) Raf-like kinase in vitro. Furthermore, AlphaFold2-directed structural modeling predicted that the identified BdMPK5-D90N chill1 mutant residue is located at the interface of BdMPK5 with the BdHT1 Raf-like kinase. BdMPK5 is a key signaling component that mediates CO2-induced stomatal movements and is proposed to function as a component of the primary CO2 sensor in grasses.

Lopez, Bryn N. K. (ORCID:0009000937288216)↗

DOE OSTI2025

MS25: Materials Science-Focused Benchmark Data Set for Machine Learning Interatomic Potentials

Here, we present MS25, a benchmark data set for evaluating machine learning interatomic potentials (MLIPs) across diverse materials-relevant systems including MgO surfaces, liquid water, zeolites, a catalytic Pt surface reaction, high-entropy alloys (HEAs), and disordered Zr-oxides. Five MLIP architectures (MACE, NequIP, Allegro, MTP, and Torch-ANI) are trained and tested, focusing not only on traditional metrics (energies, forces, and stresses) but also explicitly validating derived physical observables such as lattice constants, volumes, and reaction barriers. We find that most models reach comparable accuracy on standard error metrics across the simple systems, although equivariant MLIPs offer 1.5–2× improvements over nonequivariant MLIPs in energy and force error for structurally complex or compositionally disordered environments such as HEAs and Zr–O systems. Our analysis highlights that low errors in energy and force predictions do not guarantee reliable observables, emphasizing the necessity of explicit validation. We demonstrate limitations in cross-framework transferability, as models trained on one zeolite framework (CHA) fail to reliably generalize to predictions of structurally distinct frameworks (e.g., MFI). Size-extensive tests show some dependence on system size for MgO, resulting from forced periodicity. The HEA and Zr–O data sets are identified as challenging tests for future benchmarks and MLIP model architecture developments as they show significant differentiation in error between MLIP architectures and are still relatively difficult at 1000 training images. Moving forward, we recommend that benchmarking efforts shift their focus from marginal accuracy improvements in energy and force errors toward identifying and understanding model failure modes, rigorously assessing transferability, and evaluating how their errors affect observable predictions. For researchers looking to choose an MLIP architecture, we suggest selecting equivariant MLIP architectures if the complexity of the system is a challenge. For simple materials problems, auxiliary features such as integration with molecular dynamics engines, trade-offs between computational data set generation cost vs MLIP inference speed, and framework integration may play a more important decision factor than small differences in error metrics that are unlikely to matter for production-level research.

chemical structure↗

DOE OSTI2025

Hydrological connectivity: a review and emerging strategies for integrating measurement, modeling, and management

This review synthesizes methods for measuring, modeling, and managing hydrologic connectivity, offering pathways to improve practices and address environmental challenges (e.g., climate change) and sustainability. As a key driver of water movement and nutrient cycling, hydrologic connectivity influences flood mitigation, water quality regulation, and biodiversity conservation. However, traditional field-based methods (e.g., dye tracing), indirect measurements (e.g., runoff analysis), and remote sensing techniques (e.g., InSAR) often struggle to capture the complexity of catchment-scale interactions. Similarly, modeling approaches—including process-based and percolation theory-based models, graph theory, and entropy-based metrics—face limitations in fully representing these interconnected processes. Both modeling and measurement techniques are constrained by inadequate spatial and temporal coverage, high data demands, computational complexity, and difficulties in representing subsurface connectivity. Subsequently, we critique current management practices that prioritize isolated variables (e.g., streamflow, sediment transport) over system-wide strategies and emphasize the need for adaptive, connectivity-based approaches in water resource planning and restoration. Moving forward, we highlight the importance of interdisciplinary collaboration, technological innovations (e.g., AI-driven modeling, real-time monitoring), and integrated frameworks to improve connectivity measurement, modeling, and adaptive management to restore fragmented hydrologic networks. This integrated approach sets the stage for transformative water resource management, fostering proactive policy development and stakeholder engagement.

Dwivedi, Dipankar↗

DOE OSTI2023

Indirect searches at CMS

In these proceedings, we present several new measurements of Standard Model (SM) processes, in the Higgs sector and beyond, that push the precision frontier forward at CMS. Results are presented in the context of a framework parameterizing deviations in Higgs boson couplings, as well as in the context of SM Effective Field Theory, where new analyses targeting Higgs, top, and multi-boson processes probe an increasingly diverse set of operators. Within these frameworks, CMS is efficiently exploring a large space of new physics models. No significant deviation from SM expectations is found.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Analytical Sensitivity Analysis of a Spent Nuclear Fuel Cask

Here, we report nuclear science and engineering is a field increasingly dominated by computational studies resulting from increasingly powerful computational tools. As a result, analytical studies, which previously pioneered nuclear engineering, are increasingly viewed as secondary or unnecessary. However, analytical solutions to reduced-fidelity models can provide important information concerning the underlying physics of a problem and aid in guiding computational studies. Similarly, there is increased interest in sensitivity analysis studies. These studies commonly use computational tools. However, providing a complementary sensitivity study of relevant analytical models can lead to a deeper analysis of a problem. This work provides the analytical sensitivity analysis of the one-dimensional (1D) cylindrical mono-energetic neutron diffusion equation using the forward sensitivity analysis procedure (FSAP) developed by Cacuci. Further, these results are applied to a reduced-fidelity model of a spent nuclear fuel cask, demonstrating how computational analysis might be improved with a complementary analytic sensitivity analysis.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

A Statistician’s Overview of Physics-Informed Neural Networks for Spatio-Temporal Data

The recent success of deep neural network models with physical constraints (so-called, Physics-Informed Neural Networks, PINNs) has led to renewed interest in the incorporation of mechanistic information in predictive models. Statisticians and others have long been interested in this problem, which has led to several practical and innovative solutions dating back decades. In this overview, we focus on the problem of data-driven prediction and inference of dynamic spatio-temporal processes that include mechanistic information, such as would be available from partial differential equations, with a strong focus on the quantification of uncertainty associated with data, process, and parameters. Here, we give a brief review of several paradigms and focus our attention on Bayesian implementations given they naturally accommodate uncertainty quantification. We then show that it is straight-forward to include the Bayesian PINN (B-PINN) within the Bayesian hierarchical model (BHM) framework that has long been considered for modeling dynamic spatio-temporal processes. Such a BHM-PINN is illustrated via a simulation study in which a latent nonlinear Burgers’ equation PDE governs the dynamics of Poisson distributed spatio-temporal data. Supplementary materials for this article are available online, including a standardized description of the materials available for reproducing the work.

Bayesian↗

DOE OSTI2021

Invited Review Paper in Commemoration of Over 50 Years of Oxidation of Metals: Addressing the Role of Water Vapor on Long-Term Stainless Steel Oxidation Behavior

In the past 20 years, it has become clear that the long-term exposure of chromia-forming austenitic stainless steels to air containing ~ 10 vol% water vapor at 650–800 °C significantly increases the Cr consumption rate compared to laboratory air. However, a similar acceleration does not occur in 100% steam because O 2 is not available to form volatile CrO 2 (OH) 2 . As the research has shown, this accelerated Cr loss is particularly important for thin-walled components like gas turbine recuperators. Typically, conventional austenitic stainless steels exhibit accelerated attack (i.e., rapid formation of Fe-rich oxide nodules), while higher alloyed steels resist this type of degradation until the alloy surface becomes Cr depleted. The development of strategies and solutions for this issue are reviewed. Model Fe–Cr–Ni alloys have been used to study composition effects on this behavior. For example, increasing the Mn content up to ~ 4 wt% was not beneficial, but increasing the Ni content improved oxidation resistance in this environment. However, at 650–700 °C even highly alloyed steels showed surface Cr depletion at alloy grain boundaries resulting in Fe-rich oxide nodule formation. The path forward for this issue will require more refined mechanistic understanding and increased used of modeling to develop better application-specific lifetime models to identify the most cost-effective alloy solutions.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Lattice Boltzmann simulation of the dissolution of slag in alkaline solution using real-shape particles

Highlights: • A dissolution numerical model was proposed to capture the real dissolution kinetics of slag in alkaline solution. • The log forward dissolution rate of Si was described as a function of NBO/T and solution pH. • A threshold solid volume fraction of 0.688 was found for a voxel in 3D, 63.8% larger than that for a pixel in 2D. • The proposed dissolution numerical model provides a reliable alternative to study the dissolution kinetics of slag. A dissolution numerical model was proposed in this study to capture the real dissolution kinetics of slag in alkaline solution. It consists of three modules, i.e. (i) simulation of the initial particle parking structure of slag in alkaline solution using real-shape particles of slag, (ii) simulation of the chemical reactions between slag and solution based on the transition state theory, and (iii) simulation of the physical transport of aqueous ions using the lattice Boltzmann method. This dissolution numerical model was verified using experimental results, showing reasonable accuracy. After verification, the dissolution numerical model was implemented to study the influences of temperature and particle shape using a proper recipe of slag in alkaline solution. This recipe was designed to avoid solid phase precipitation or gel formation via thermodynamic analysis. The simulation results showed faster dissolution kinetics of slag when using higher temperatures and more irregular particle shapes.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Single spin asymmetry in forward p A collisions from the Pomeron-odderon interference

Working in the hybrid framework of the high energy p A collisions we identify a new contribution to transverse single spin asymmetry (SSA). The phase necessary for the SSA is provided by the Pomeron-odderon interference in the dense nuclear target. The complete formula for the p A → h X polarized cross section also contains the transversity distribution for the polarized projectile as well as the real part of the twist-three fragmentation function. We numerically estimate the asymmetry A N and its nuclear dependence. Based on a model computation we find that A N can be a percent level in the forward and low- P h ⊥ region. For large nuclei we find significant suppression, with A N ∝ A − 7 / 6 parametrically. As a notable feature we find a node of A N as a function of the P h ⊥ around the values of the initial saturation scale that could be used to test this mechanism experimentally. Published by the American Physical Society 2025

Benić, Sanjin↗

DOE OSTI2026

torch-einshard v1.0

torch-einshard is a Python library for describing local and distributed PyTorch tensor computations with compact, einsum-like notation. Its expressions name logical axes, specify how they are sharded across a PyTorch DeviceMesh, and represent partial reductions. The library automatically performs contractions, permutations, reshaping, splitting, gathering, reduction, reduce-scatter, and repartitioning while preserving autograd. Additional features include sharding-aware FFTs, tensor rolls, halo exchange, sliding windows, 1D–3D convolutions, uneven-shard handling, parameter initialization and gradient management, and cost-based execution planning. It is designed for scientific machine learning and large-model workloads, including tensor-, sequence-, and spatial-parallel MLPs, attention, convolutions, and spectral operations. Compared with manually combining torch.einsum and distributed collectives, torch-einshard expresses both the mathematical operation and data placement in one readable formula. This reduces boilerplate and synchronization errors, keeps forward and backward communication consistent, and allows the library to select optimized collective strategies without changing model code.

Morozov, Dmitriy [Lawrence Berkeley National Labor↗

DOE OSTI2022

Statistics of Green's functions on a disordered Cayley tree and the validity of forward scattering approximation

The accuracy of the forward scattering approximation for two-point Green's functions of the Anderson localization model on the Cayley tree is studied. A relationship between the moments of the Green's function and the largest eigenvalue of the linearized transfer-matrix equation is proved in the framework of the supersymmetric functional-integral method. The new large-disorder approximation for this eigenvalue is derived and its accuracy is established. Using this approximation the probability distribution of the two-point Green's function is found and compared with that in the forward scattering approximation (FSA). It is shown that FSA overestimates the role of resonances and thus the probability for the Green's function to be significantly larger than its typical value. The error of FSA increases with increasing the distance between points in a two-point Green's function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗