SEARCH · Search NASA

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2024

Relations between Shot Noise, Gain Bandwidth, and Saturation of Instabilities

There are numerous instabilities present in charged particle beams that undergo exponential growth and reach saturation. In various applications, such as free-electron lasers or micro-bunching light sources, achieving saturation is desirable. Conversely, there are applications where these instabilities are utilized as linear broad-band amplifiers for signals embedded in the charged beam. In the latter scenario, the saturation of an instability induces non-linear distortions in the imprinted signal, thereby limiting the useful range of such amplifiers. Accurate evaluation of these instabilities necessitates a complete and comprehensive modeling approach that includes shot noise within the beam. Unfortunately, such modeling is not always feasible or practical. In this paper, we introduce a methodology utilizing the frequency and bandwidth of the instability as key parameters. Through this, we derive an estimation for the range of linear instability growth. Our derivation is conducted in a model-independent manner, making it applicable to a broad spectrum of instabilities. To validate our approach, we employ established and thoroughly benchmarked simulations with a free electron laser (FEL) code as well as self-consistent 3-dimensional simulation of plasma-cascade instability using code SPACE.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Impacts of Bulk Microphysics Scheme Structural Choices on Simulations of Rain Initiation Through Drop Coalescence

This study examines how different structural choices in bulk microphysics schemes impact the simulation of warm rain initiation. A single liquid category (SLC) approach prognosing up to four moments of a single drop size distribution (DSD) is compared to the traditional two-category, two-moment approach with separate DSDs for cloud and rain (four total prognostic variables). Different methods for calculating tendencies of the prognostic variables from drop collision-coalescence are also tested: a discretized numerical-integration approach, machine learning via neural networks, lookup tables, and traditional power law fits. Relative to simulations using a bin microphysics model, SLC gives smaller error overall than the two-category approach when numerical integration is used to calculate the collision-coalescence tendencies for both. Replacing the numerical integration with a pre-computed lookup table reduces computational cost with little loss of accuracy. However, using fitted power laws with SLC to represent the collision-coalescence tendencies substantially reduces accuracy and leads to an order of magnitude increase in error. It is also demonstrated that with SLC, reasonably accurate solutions are obtained using only three prognostic moments, while a two-moment SLC scheme leads to substantial error. Overall, both the choice of prognostic moments (e.g., SLC vs. two-category) and method to calculate the collision-coalescence tendencies are important to consider for minimizing errors in bulk schemes. SLC with a sufficiently detailed calculation of the collision-coalescence tendencies provides accurate solutions for a reasonable computational cost, providing a viable alternative to the traditional two-category, two-moment approach for bulk microphysics.

320 (cloud physics and chemistry)↗

DOE OSTI2024

Problem-tailored Simulation of Energy Transport on Noisy Quantum Computers

The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L = 12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Numerical modeling and experimental validation of low velocity impact of woven GFRP/CFRP composites

Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. Further, the comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Hybrid quantum simulations with qubits and qumodes on trapped-ion platforms

We explore the feasibility of gate-based hybrid quantum computing using both discrete (qubit) and continuous (qumode) variables on trapped-ion platforms. Trapped-ion systems have demonstrated record one- and two-qubit gate fidelities and long qubit coherence times, while qumodes, which can be represented by the collective vibrational modes of the ion chain, have remained relatively unex- plored for their use in computing. Using numerical simulations, we show that high-fidelity hybrid gates and measurement operations can be achieved for existing trapped-ion quantum platforms. As an exemplary application, we consider quantum simulations of the Jaynes-Cummings-Hubbard model, which is given by a one-dimensional chain of interacting spin and boson degrees of free- dom. Using classical simulations, we study its real-time evolution and develop a suitable variational quantum algorithm for ground state preparation. Furthermore, our results motivate further studies of hybrid quantum computing in this context, which may lead to direct applications in condensed matter and fundamental particle and nuclear physics.

Lower-dimensional field theories↗

DOE OSTI2024

Chromosome-level genome assemblies and genetic maps reveal heterochiasmy and macrosynteny in endangered Atlantic Acropora

Abstract Background Over their evolutionary history, corals have adapted to sea level rise and increasing ocean temperatures, however, it is unclear how quickly they may respond to rapid change. Genome structure and genetic diversity contained within may highlight their adaptive potential. Results We present chromosome-scale genome assemblies and linkage maps of the critically endangered Atlantic acroporids,Acropora palmataandA. cervicornis. Both assemblies and linkage maps were resolved into 14 chromosomes with their gene content and colinearity. Repeats and chromosome arrangements were largely preserved between the species. The family Acroporidae and the genusAcroporaexhibited many phylogenetically significant gene family expansions. Macrosynteny decreased with phylogenetic distance. Nevertheless, scleractinians shared six of the 21 cnidarian ancestral linkage groups as well as numerous fission and fusion events compared to other distantly related cnidarians. Genetic linkage maps were constructed from oneA. palmatafamily and 16A. cervicornisfamilies using a genotyping array. The consensus maps span 1,013.42 cM and 927.36 cM forA. palmataandA. cervicornis, respectively. Both species exhibited high genome-wide recombination rates (3.04 to 3.53 cM/Mb) and pronounced sex-based differences, known as heterochiasmy, with 2 to 2.5X higher recombination rates estimated in the female maps. Conclusions Together, the chromosome-scale assemblies and genetic maps we present here are the first detailed look at the genomic landscapes of the critically endangered Atlantic acroporids. These data sets revealed that adaptive capacity of Atlantic acroporids is not limited by their recombination rates. The sister species maintain macrosynteny with few genes with high sequence divergence that may act as reproductive barriers between them. In the AtlanticAcropora, hybridization between the two sister species yields an F1 hybrid with limited fertility despite the high levels of macrosynteny and gene colinearity of their genomes. Together, these resources now enable genome-wide association studies and discovery of quantitative trait loci, two tools that can aid in the conservation of these species.

Biotechnology & Applied Microbiology↗

DOE OSTI2024

Optimization of the artificial viscosity in Lagrangian staggered discretization codes. Modeling 1D stand-alone shock - case study

We have developed new measures of errors for numerical shock. The new approach is based on analysis of the structure function, and separation of the errors related to oscillations and shock width, which also include error in the position of the ”center” of the numerical shock. We have demonstrated that those measures correctly characterize the numerical solution. We introduced an objective function in, which both types of errors are weighted, and presented optimal values of the coefficients of the linear and quadratic viscosity for different weights and different Mach numbers.

97 MATHEMATICS AND COMPUTING↗