SEARCH · Engineering Papers

Results for “Numerical relativity”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Ensemble Simulation Techniques and Fast Randomized Algorithms

The major goals of the project were to develop and analyze new ensemble simulation techniques, including trajectory stratification and preconditioned MCMC techniques, as well as develop fast numerical linear algebra techniques closely related to ensemble simulation ideas. The trajectory stratification techniques involve simulating in parallel short trajectory fragments of a Markov process confined to a specific region of space‐time and then patching together the statistics gathered to assemble estimates of very general dynamical properties. We have also developed this approach for rare event simulation and extended the techniques to applications requiring a more general framework (such as electronic structure calculations). The preconditioned MCMC techniques involve simulating multiple Markov chains in parallel and then using information from the ensemble to speed the mixing of each individual chain. The fast randomized linear algebra methods are motivated by the diffusion Monte Carlo technique, but are applicable to finding the dominant eigenvalue of (almost) general matrices. For most non‐negative matrices, the schemes result in an error (compared to the power method) that is constant in the dimension of the problem. For more general matrices, we see a very clear sublinear cost trend in computational tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Finite-volume quantization condition from the N/D representation

We propose a new model-independent method for determining hadronic resonances from lattice QCD. The formalism is derived from the general principles of unitarity and analyticity, as encoded in the N/D representation of a partial-wave two-body amplitude. The associated quantization condition relates the finite-volume spectrum to the infinite-volume numerator $\mathscr{N}$, used to reconstruct the scattering amplitude from dispersive relations. Unlike the original Lüscher condition, this new formalism is valid for energies coinciding with the left-hand cuts from arbitrary one- and multiparticle exchanges.

Dawid, Sebastian M. [Univ. of Washington, Seattle,↗

DOE OSTI2026

Structure-aware Initialization via Numerical Continuation and Informed Priors

Scientific machine learning (SciML) often operates in ill-conditioned, weakly identifiable regimes due to limited data or indirect observations. In such settings, optimization and inference are highly sensitive to the starting point, making initialization--often under-reported--a consequential degree of freedom. Random initialization is not a neutral default as it induces an implicit prior over candidate solutions and can systematically bias the result, producing large run-to-run variability. Here, we formalize this view by treating initialization as a hidden confounder in SciML and develop a unifying theory for structure-aware initialization via numerical continuation, constructing warm starts from related problem instances. Across representative tasks, including physics-informed neural networks, maximum likelihood estimation, and variational inference, warm starts have been shown to consistently reduce optimization effort and improve reliability.

Data integrity↗

DOE OSTI2025

Data from TropiRoot 1.0 database: tropical root characteristics across environments

TropiRoot 1.0 is a new tropical root database with root characteristics across environment gradients. It has data extracted from 104 new sources, resulting in more than 8000 rows of data (either species or community data). Most of the data in TropiRoot 1.0 includes root characteristics such as root biomass, morphology, root dynamics, mass fraction, architecture, anatomy, physiology and root chemistry. This initiative represents an approximately 30% increase in the currently available data for tropical roots in the Fine Root Ecology Database (FRED). TropiRoot 1.0, contains root characteristics from 25 different countries where seven are located in Asia, six in South America, five in Central America and the Caribbean, four in Africa, two in North America, and 1 in Oceania. Due to the volume of data, when ancillary data was available, including soil data, these data was either extracted and included in the database or their availability was recorded in an additional column. Multiple contributors checked the entries for outliers during the collation process to ensure data quality. For text-based observations, we examined all cells to ensure that their content relates to their specific categories. For numerical observations, we ordered each numerical value from least to greatest and plotted the values, checking apparent outliers against the data in their respective sources and correcting or removing incorrect or impossible values. Some data (soil and aboveground) have different columns for the same variable presented in different units, including originally published units, but root characteristics data had units converted to match the ones reported in FRED. By filling a gap from global databases, TropiRoot 1.0 expands our knowledge of otherwise so far underrepresented regions, and our ability to assess global trends. This advancement can be used to improve tropical forest representation in vegetation models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

TropiRoot 1.0: Database of tropical root characteristics across environments

Tropical ecosystems contain the world's largest biodiversity of vascular plants. Yet, our understanding of tropical functional diversity and its contribution to global diversity patterns is constrained by data availability. This discrepancy underscores an urgent need to bridge data gaps by incorporating comprehensive tropical root data into global datasets. Here, we provide a database of tropical root characteristics. This new database, TropiRoot 1.0, will be instrumental in evaluating an array of hypotheses pertaining to root functional ecology and plant biogeography, both within the tropics and relative to other global biomes. The data compilation was conducted by the TropiRoot Initiative, in partnership with the Fine-Root Ecology Database (FRED) and the Global Root Trait (GRooT) database, Colorado State University (CSU) and the Smithsonian Tropical Research Institute (STRI). Literature search and data extraction were conducted between 2020 and 2024. Literature was identified using Web of Science, Scopus, and complemented using the expert knowledge of members of TropiRoot. To provide broad environmental and geographical distributions, literature searches included root characteristics (traits) across global change drivers, natural gradients, and from different continents. We adopted FRED standardized data columns and streamlined the format to enhance accessibility for data extraction across various user groups. This optimized framework resulted in a smaller, yet comprehensive datasheet. To make the database compatible with other global root trait initiatives, column identification was standardized following the codes provided by FRED. These efforts culminated in data extracted from 104 new sources, resulting in more than 8000 rows of data (either species or community data). Most of the data in TropiRoot 1.0 include root characteristics such as root biomass, morphology, root dynamics, mass fraction, architecture, anatomy, physiology, and root chemistry. This initiative represents a 30% increase in the currently available data for tropical roots in FRED. TropiRoot 1.0 contains root characteristics from 25 different countries, where seven are located in Asia, six in South America, five in Central America and the Caribbean, four in Africa, two in North America, and 1 in Oceania. Due to the volume of data, when ancillary data were available, including soil data, these data were either extracted and included in the database or its availability was recorded in an additional column. Multiple contributors checked the entries for outliers during the collation process to ensure data quality. For text-based observations, we examined all cells to ensure that their content relates to their specific categories. For numerical observations, we ordered each numerical value from least to greatest and plotted the values, checking apparent outliers against the data in their respective sources and correcting or removing incorrect or impossible values. Some data (soil and aboveground) have different columns for the same variable presented in different units, including originally published units, but root characteristics data had units converted to match those reported in FRED. By filling a gap from global databases, TropiRoot 1.0 expands our knowledge of otherwise so far underrepresented regions and our ability to assess global trends. This advancement can be used to improve tropical forest representation in vegetation models. The data are freely available and should be cited when used.

FRED↗

DOE OSTI2024

Beyond the Debye–Hückel limit: Toward a general theory for concentrated electrolytes

The phenomenon of underscreening in concentrated electrolyte solutions leads to a larger decay length of the charge–charge correlation than the prediction of Debye–Hückel (DH) theory and has found a resurgence of both theoretical and experimental interest in the chemical physics community. To systematically understand and investigate this phenomenon in electrolytes requires a theory of concentrated electrolytes to describe charge–charge correlations beyond the DH theory. We review the theories of electrolytes that can transition from the DH limit to concentrations where charge correlations dominate, giving rise to underscreening and the associated Kirkwood Transitions (KTs). In this perspective, we provide a conceptual approach to a theoretical formulation of electrolyte solutions that exploits the competition between molecular-informed short-range (SR) and long-range interactions. We demonstrate that all deviations from the DH limit for real electrolyte solutions can be expressed through a single function ΣQ that can be determined both theoretically and numerically. Importantly, ΣQ can be directly related to the details of SR interactions and, therefore, can be used as a tool to understand how differences in representations of interaction can influence collective effects. The precise function form of ΣQ can be inferred through a Gaussian field theory of both the number and charge densities. The resulting formulation is validated by experiment and can accurately describe the collective phenomenon of screening in concentrated bulk electrolytes. Importantly, the Gaussian field theory predictions of the screening lengths appear to be less than ∼1 nm at concentrations above KTs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Local adaptation of both plant and pathogen: an arms‐race compromise in switchgrass rust

In coevolving species, parasites locally adapt to host populations as hosts locally adapt to resist parasites. Parasites often outpace host local adaptation since they have rapid life cycles, but host diversity, the strength of selection, and external environmental influence can result in complex outcomes. To better understand local adaptation in host–parasite systems, we examined locally adapted switchgrass (Panicum virgatum), and its leaf rust pathogen (Puccinia novopanici) across a latitudinal range in North America. We grew switchgrass genotypes in 10 replicated multiyear common gardens, measuring rust severity from natural infection in a ‘host reciprocal transplant’ framework for testing local adaptation. We conducted genome-wide association mapping to identify genetic loci associated with rust severity. Genetically differentiated rust populations were locally adapted to northern and southern switchgrass, despite host local adaptation to environmental conditions in the same regions. Rust resistance was polygenic, and distinct loci were associated with rust severity in the north and south. We narrowed a previously identified large-effect quantitative trait locus for rust severity to a candidate YELLOW STRIPE-LIKE gene and linked numerous other loci to defense-related genes. Overall, our results suggest that both hosts and parasites can be simultaneously locally adapted, especially when parasites impose less selection than other environmental factors.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Uptake and elimination of 137Cs in aquatic biota inhabiting a contaminated effluent canal

Radiocesium (137Cs) is among the most prevalent and concerning radionuclide released from anthropogenic activities due to its relatively slow decay time and high bioavailability. Numerous studies have reported activity concentrations of 137Cs in exposed biota, but few have assessed detailed toxicokinetics. The objective of our study was to quantify and compare the uptake and elimination of 137Cs in aquatic taxa from multiple classes with different ecological roles inhabiting a canal contaminated from nuclear production activities on the Savannah River Site in South Carolina, USA. We tracked accumulation and elimination of 137Cs in four understudied aquatic species: one aquatic macrophyte [American white water lily (Nymphae odorata), two benthic associated omnivorous species (American bullfrog (Lithobates catesbeianus) and red swamp crayfish (Procambarus clarkii), and one pelagic carnivorous fish (Eastern mosquitofish (Gambusia holbrooki)]. We monitored accumulation of 137Cs from initial exposure until equilibrium for biota restricted within the canal and subsequently tracked elimination of 137Cs over time by moving organisms from the canal to a nearby reference pond. Results indicated bullfrog tadpoles reached 95 % equilibrium in approximately 9 days, the fastest among the species evaluated. After bullfrog tadpoles, the next fastest time to equilibrium was achieved by red swamp crayfish (51 days) followed by mosquitofish (86 days). At equilibrium, bullfrog tadpoles were found to have similar 137Cs activity concentrations to mosquitofish (3.16 and 3.57 Bq/g, respectively), and both were significantly higher than red swamp crayfish (1.77 Bq/g). For elimination, we found biological half-lives to range between 8 and 43 days with bullfrog tadpoles eliminating 137Cs the fastest followed by water lilies, mosquitofish, and red swamp crayfish. Overall, we found the uptake and elimination of 137Cs to vary widely among our studied species, but benthic organisms are not at an increased risk to accumulate 137Cs compared to organisms in less frequent contact with the sediment.

Benthic organism Bioaccumulation Bioindicator Elim↗

DOE OSTI2024

ZFP: A compressed array representation for numerical computations

HPC trends favor algorithms and implementations that reduce data motion relative to FLOPS. We investigate the use of lossy compressed data arrays in place of traditional IEEE floating point arrays to store the primary data of calculations. Simulation is fundamentally an exercise in controlled approximation, and error introduced by finite-precision arithmetic (or lossy compression) is just one of several sources of error that need to be managed to ensure sufficient accuracy in a computed result. We describe ZFP, a compressed numerical format designed for in-memory storage of multidimensional arrays, and summarize theoretical results that demonstrate that the error of repeated lossy compression can be bounded and controlled. Furthermore, we establish a relationship between grid resolution and compression-induced errors and show that, contrary to conventional floating point, ZFP reduces finite-difference errors with finer grids. We present example calculations that demonstrate data reduction by 4x or more with negligible impact on solution accuracy. Our results further demonstrate several orders-of-magnitude increase in accuracy using ZFP over IEEE floating point and Posits for the same storage budget.

Lindstrom, Peter↗

DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗

DOE OSTI2024

Elemental profiling and genome-wide association studies reveal genomic variants modulating ionomic composition in Populus trichocarpa leaves

The ionome represents elemental composition in plant tissues and can be an indicator of nutrient status as well as overall plant performance. Thus, identifying genetic determinants governing elemental uptake and storage is an important goal for breeding and engineering biomass feedstocks with improved performance. In this study, we coupled high-throughput ionome characterization of leaf tissues with high-resolution genome-wide association studies (GWAS) to uncover genetic loci that modulate ionomic composition in leaves of poplar ( Populus trichocarpa ). Significant agreement was observed across the three ionomic profiling platforms tested: inductively coupled plasma-mass spectrometry (ICP-MS), neutron activation analysis (NAA) and laser-induced breakdown spectroscopy (LIBS). Relative quantification of 20 elements using ICP-MS across a population of 584 genotypes, revealed larger variation in micro-nutrients and trace elements content than for macro-nutrients across genotypes. The GWAS performed using a set of high-density (>8.2 million) single nucleotide polymorphisms, identified over 600 loci significantly associated with variations in these mineral elements, pointing to numerous uncharacterized candidate genes. A significant enrichment for genes related to ion homeostasis and transport was observed, including several members of the cation-proton antiporters (CPA) family and MATE efflux transporters, previously reported to be critical for plant growth and fitness in other species. Our results also included a polymorphic copy of the high-affinity molybdenum transporter MOT1 found directly associated to molybdenum content. For the first time in a perennial plant, our results provide evidence of genetic control of mineral content in a model tree species.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

High pressure ammonia/methanol oxidation up to 100 atm

Here, high pressure ammonia/methanol oxidation and NOx formations were investigated using a recently developed supercritical pressure jet-stirred reactor (SP-JSR) at 20 and 100 atm with temperatures between 550 and 950 K and equivalence ratios of 0.138 and 1.15. The experimental results show that NH 3 oxidation at high pressure is significantly accelerated by the active OH radicals produced from CH 3 OH oxidation. Furthermore, the kinetic interactions between NH 3 and CH 3 OH are governed mainly by the reactions CH 3 OH + NH 2 = CH 2 OH + NH 3 , CH 3 OH + NH 2 = CH 3 O + NH 3 , and CH 2 O + NH 2 = HCO + NH 3 . A HP-Mech model for high-pressure NH 3 /CH 3 OH oxidation was developed in this study. It consists of the most recent NH 3 and CH 3 OH models including some new reactions and updated rate constants from the literature as well as NH 3 -CH 3 OH interactions where rate constants of CH 3 OH + NH 2 = CH 2 OH + NH 3 , CH 3 OH + NH 2 = CH 3 O + NH 3 , NH 2 + CH 2 O = NH 3 + HCO, and NH 2 + CH 2 O = NH 2 CHO + H were theoretically calculated in this study. Our model with these updates improves the prediction for the measured N 2 O/NO x temperature dependence at 100 atm. In addition, the reaction pathway and sensitivity analysis show that N 2 O/NO x /HONO interactions with HO 2 are very important, especially for a fuel-lean mixture at 100 atm. The HONO mole fraction for the fuel-lean mixture at 100 atm was then measured by off-axis integrated cavity output spectroscopy (ICOS) at wavenumber of 6638.26 cm -1 . The experimental data show a significant HONO formation at intermediate temperature that is strongly underpredicted by numerical simulation at 100 atm. Therefore, the HONO related reactions with notable uncertainty at high pressure such as NO + OH (+M) = HONO (+M) and H 2 NO + NO 2 = HONO + HNO need deeper exploration in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

PETSc/TAO Users Manual Revision 3.22

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations (PDEs) and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for implementing large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Nuclear Criticality Safety Margin and Handbook Data: Concepts and Applications [Slides]

At the end of this briefing, personnel should be able to: • Define concepts associated with criticality safety margin • Provide and understand practical applications that illustrate safety margin concepts • Describe how safety margin is addressed in criticality safety evaluations • Improve documentation and communication of safety margin • Connect these concepts to related statements in the ANS-8 Standards • Find numerous applications for handbook data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

PETSc/TAO Users Manual Revision 3.23

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

PETSc/TAO Users Manual Revision 3.24

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

PETSc/TAO Users Manual Revision 3.25

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

DOE OSTI2026

Current-Phase Relations of Superconducting PtSi Constriction Josephson Junctions

Here, we study the current-phase relations of superconductor–constriction–superconductor Josephson junctions made from platinum silicide thin films by fabricating dc-superconducting quantum interference devices incorporating pairs of junctions and measuring their magnetic-field-dependent electrical transport. By comparing the supercurrent interference patterns with numerical Ginzburg–Landau simulations, we extract the current-phase relations of individual junctions and quantify the degree of nonlinearity. These measurements show that while the constrictions themselves sustain substantial nonlinearity, the effect is moderated by the kinetic inductance of the device leads, which is an important consideration for practical superconducting circuit applications.

77 NANOSCIENCE AND NANOTECHNOLOGY↗