Search NASA⌕ Search

SEARCH · Search NASA

Results for “approximation”

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 361 records · Page 20

Approach to startup inventory for viable commercial fusion power plant

With the increasing efforts to commercialize fusion power, private and government organizations are investing heavily in the development of technology to support a viable fusion power plant. Deuterium-Tritium (DT) fueled reactors are more prevalent than other proposed designs, requiring tritium processing and handling technology for safe operations and self-sufficiency. Further, each fusion power plant will need a specific-to-design startup inventory of tritium to begin operations. This startup inventory is required prior to breeding and is the minimum tritium inventory required to fill each processing component in the fuel cycle, to offset radioactive decay losses, and to avoid a zero-fuel situation for continuous operation. We present an approach to calculate the startup tritium inventory for a 500 MW th reactor, with considerations for reserve inventory for maintenance and commissioning. A baseline startup inventory was calculated to be approximately 327 gs. This value was obtained using modest assumptions about the technology and operating parameters of a fusion power plant. The required operating reserve inventory or the inventory necessary to keep a fusion power plant operational using only direct internal recycling for 24 h for the same plant design is approximately 642 gs. The approach and findings of this paper will enable fusion energy stakeholders to better utilize the existing scarce global tritium supply.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled feldspar dissolution and secondary mineral precipitation in batch systems: 6. Labradorite dissolution, calcite growth, and clay precipitation at 60 °C and pH 8.2–8.4

Here, we conducted experiments on concurrent labradorite dissolution, calcite precipitation, and clay precipitation in batch reactor systems and tracked reaction processes using multiple isotope tracers. Labradorite was chosen for its role as a major and reactive component in basalt; the experiments thus directly impact our understanding of CO 2 storage in basalt aquifers and enhanced rock weathering. We doped initial solutions with 29 Si, 43 Ca, and Ca 13 CO 3 (s). Experiments were conducted at 60 °C and pH ~ 8.3 for up to 840 h, with isotope ratios in the experimental aqueous solutions measured using MC-ICP-MS. Unidirectional rates of labradorite dissolution near equilibrium were approximately two orders of magnitude slower than far-from-equilibrium rates reported in the literature. Calcite growth occurred near equilibrium and the rates were limited by the labradorite dissolution rates. In the steady state phase, the interplay of these three heterogeneous reactions—labradorite dissolution, calcite growth, and clay precipitation—results in a coupled system that approaches a near-equilibrium state. The system does not reach true equilibrium because labradorite continues to dissolve, albeit at a much slower rate near equilibrium. The overall reaction can be approximated as, Na 0.4 Ca 0.6 Al 1.6 Si 2.4 O 8 + 0.6HCO 3 - + 1·.7H 2 O + 0.4H + → 0.4Na + + 0.6CaCO 3(s) + 0.5Al 2 Si 2 O 5 (OH) 4(s) + 0.6Al(OH) 4 - + 1.4SiO 2 o (aq). The experimental results show that using short-term far-from-equilibrium rate constants would lead to an overestimation of feldspar weathering rates at the Earth’s surface (e.g., basalt weathering and enhanced rock weathering) and CO 2 mineralization in basalt aquifers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗

Superconductivity in WBe 2

WBe 2 , which occurs in space group 194, with hexagonal symmetry P6 3 /mmc, is prepared by arc-melting at temperatures above 2200 C, where Be vapor loss is significant. This study is motivated by recent work on MoB 2 and WB 2 , both superconductors (T c =32 and 17 K respectively) under high (∼70 GPa) pressure. In order to avoid the known Be-rich superconducting phases (WBe 13 and WBe 22 ) in the complex phase diagram, both known to be superconducting at 4.1 K, the sample was prepared with a slight (∼5 %) excess of W. The resultant sample, prepared using high purity (99.999 %) Be, is essentially single phase WBe 2 , with some spread in its superconducting properties due to the known homogeneity range. (WBe 2 forms in space group 194 between approximately W 1.02 Be 1.98 and W 0.88 Be 2.12 .) Characterization was carried out with x-ray diffraction, electrical resistivity, ρ, in zero and applied magnetic fields, and specific heat. The resistivity in zero and applied fields and specific heat data indicate that our sample of WBe 2 is a bulk superconductor at ambient pressure with a T c onset in ρ at 1.05 K, and T c (ρ→0) at ∼0.86 K. Furthermore, there is no signature of superconductivity in ρ at 4.1 K, indicating successful avoidance of WBe 13 and WBe 22 . The ρ data in field indicate an upper critical field of approximately 400 gauss.

Electrical transport↗

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗

RandONets: Shallow networks with random projections for learning linear and nonlinear operators

Deep neural networks have been extensively used for the solution of both the forward and the inverse problem for dynamical systems. However, their implementation necessitates optimizing a high-dimensional space of parameters and hyperparameters. This fact, along with the requirement of substantial computational resources, pose a barrier to achieving high numerical accuracy, but also interpretability. Here, to address the above challenges, we present Random Projection-based Operator Networks (RandONets): shallow networks with random projections and tailor-made numerical analysis methods that learn accurately and fast linear and nonlinear operators. Building on previous works, we prove that RandOnets are universal approximators of linear and nonlinear operators. Due to their simplicity, RandONets provide a one-step transformation of the input space, facilitating interpretability. For the evaluation of their performance, we focus on operators of PDEs. We show, that RandONets outperform by several orders of magnitude, both in terms of numerical approximation accuracy and computational cost, the “vanilla” DeepONets. Hence, we believe that our method will trigger further developments in the field of scientific machine learning, for the development of new ‘’light”schemes that will provide high accuracy while reducing dramatically the computational cost. A MATLAB toolbox for RandONets, including demos, is available on GitHub at https://github.com/GianlucaFabiani/RandONets.

Interpretable machine learning↗

Multi-plane moment-of-fluid interface reconstruction in 3D

Moment-of-fluid (MOF) methods for interface reconstruction approximate the region occupied by material in each mesh element only through reference to its geometric moments. Here, we present a 3D MOF method that represents the material (POM) in each cell as the convex intersection of the cell and multiple half-spaces, each selected to minimize the least-squares error between computed moments of the approximated material and provided reference moments. This optimization problem is highly non-linear and non-convex, making the numerical result very sensitive to the initial guess. To create an effective initial guess in each cell, we construct an ellipsoid from 0th–2nd order reference moments such that its shape corresponds with that of the POM. Within this ellipsoid we inscribe a polyhedron, and initialize the minimization problem with the half-spaces defined by each of its faces. The inscribed polyhedron has minimally 4 faces, and using up to 3rd order moments permits optimization over up to 20 unknown values. We therefore define MOF methods that utilize 4, 5, or 6 half-spaces, correspondingly initialized with the faces of a single inscribed tetrahedron, triangular prism, or hexahedron. Stability of the non-linear optimization is further improved with a prepossessing step that normalizes the reference moments according to the axes of the reference ellipsoid. Using this approach, the non-linear least-squares solver reliably converges to a near-global minimum from a single initial guess. We demonstrate accuracy and robustness using single-cell and multi-cell examples over a wide spectrum of geometry. In particular, we demonstrate our ability to exactly reproduce several important and complex features defined by up to four half-spaces, such as corners, filaments, filament tips, and embedded material in the cell.

3D interface reconstruction↗

Solving high-dimensional partial integral differential equations: The finite expression method

Partial integro-differential equations (PIDEs) have broad applications in the sciences, from electro-magnetism to options pricing. Here, in this paper, we introduce a new finite expression method (FEX) to solve PIDEs. This approach builds upon the original FEX and its inherent advantages with new advances: 1) A novel method of parameter grouping is proposed to reduce the number of coefficients in high-dimensional function approximation; 2) A Taylor series approximation method is implemented to significantly improve the computational efficiency and accuracy of the evaluation of the integral terms of PIDEs. The new FEX based method, denoted FEX-PG to indicate the addition of the parameter grouping (PG) step to the algorithm, provides both high accuracy and interpretable numerical solutions, with the outcome being an explicit equation that facilitates intuitive understanding of the underlying solution structures. These features are often absent in traditional methods, such as finite element methods (FEM) and finite difference methods, as well as in deep learning-based approaches. To benchmark our method against recent advances, we apply the new FEX-PG to solve benchmark PIDEs in the literature. In high-dimensional settings, FEX-PG exhibits strong and robust performance, achieving relative errors on the order of single precision machine epsilon, significantly outperforming existing approaches based on neural networks.

Combinatorial optimization↗

A method for bounding high-order finite element functions: Applications to mesh validity and bounds-preserving limiters

We introduce a novel method for bounding high-order multi-dimensional polynomials in finite element approximations. The method involves precomputing optimal piecewise-linear bounding boxes for polynomial basis functions, which can then be used to locally bound any combination of these basis functions. This approach can be applied to any element/basis type at any approximation order, can provide local (i.e., subcell) extremum bounds to a desired level of accuracy, and can be evaluated efficiently on-the-fly in simulations. Furthermore, we show that this approach generally yields more accurate bounds in comparison to traditional methods based on convex hull properties (e.g., Bernstein polynomials). Furthermore, the efficacy of this technique is shown in applications such as mesh validity checks and optimization for high-order curved meshes, where positivity of the element Jacobian determinant can be ensured throughout the entire element, and continuously bounds-preserving limiters for hyperbolic systems, which can enforce maximum principle bounds across the entire solution polynomial.

Bounding box↗

Non -degenerate marginal-likelihood calibration with application to quantum characterization

Here, we propose a marginal likelihood strategy within the Kennedy-O’Hagan (KOH) Bayesian framework, where a Gaussian process (GP) models the discrepancy between a physical system and its simulator. Our approach introduces a novel marginalized likelihood by integrating out the degenerate eigenspace of the covariance matrix, rather than approximating the original likelihood. Unlike approximation methods that compromise accuracy for computational efficiency, our method defines an exact likelihood—distinct from the original but preserving all relevant information. This formulation achieves computational efficiency and stability, even for large datasets where the covariance matrix nears degeneracy. Applied to the characterization of a superconducting quantum device at Lawrence Livermore National Laboratory, the approach enhances the predictive accuracy of the Lindblad master equations for modeling Ramsey measurement data by effectively quantifying uncertainties consistent with the quantum data.

general physics↗

Valence instability and collapse of ferromagnetism in EuB 6 at high pressures

Despite the simplicity of their cubic crystal lattice, rare-earth hexaborides display complex physical properties including a (long debated) onset of metallization via magnetic polaron formation at T c1 approximate to 15 K preceding ferromagnetic ordering at T c2 approximate to 12 K. Here, in this work, we used applied pressure to tune the interplay between electronic structure and magnetism in EuB6. We probed the magnetism, valence, and structure of EuB 6 under quasi-hydrostatic pressures up to 30 GPa using X-ray techniques. Our findings show evidence for collapse of ferromagnetism above 20 GPa following a monotonic increase of mean Eu valence. While X-ray diffraction measurements in the paramagnetic state at room temperature show that the lattice retains cubic symmetry, a measurable quadrupole interaction seen by time-domain synchrotron Mössbauer spectroscopy suggests a lowering of symmetry associated with magnetic ordering, becoming more prominent across the magnetic transition. The interplay between conduction band electron count and magnetism observed under applied pressure in EuB 6 opens possibilities for fine-tuning metallization and magnetic properties of similar Eu-based semi-metal systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Osmotic control of the spacing of parallel shear cracks in shale growing subcritically in geologic past

The geological genesis of natural cracks in sedimentary rocks such as shale is a problem that needs to be understood to improve the technology of hydraulic fracturing as well as deep sequestration of harmful fluids. Why are the vertical natural cracks roughly parallel and equidistant, and why is the spacing roughly 10 cm rather than 1 cm or 100 cm? Fracture mechanics of critical cracks cannot answer this question. Neither can the material heterogeneity. The growth of critical parallel cracks is impossible because the relative crack face displacements would immediately localize into one crack, leading to an earthquake. The cracks must have formed, on the tectonic time scale, by a slow growth of subcritical shear cracks governed by the Charles-Evans law. The idea advanced here is that what controls the crack spacing is the balance between the reduction, due to shear dilatancy, of the concentration of ions such as Na + and Cl - in each fracture process zone (PFZ), which decelerates the cracks, and the restoration of ion concentration by diffusion of ions from the space between the cracks into the FPZ. This diffusion of water is driven mainly by the osmotic pressure gradient, which offsets the deceleration and depends strongly on the crack spacing. A simple analytical solution of the steady state is rendered possible by approximating the ion concentration profiles between adjacent cracks by parabolic arcs. Applying this theory to Woodford shale yields the approximate crack spacing of 10 cm, which is realistic. Furthermore, the stability of unlimited parallel mode II frictional crack growth is proven by examining the second variation of the free energy. Water concentration drop in the FPZ due to shear dilatancy and its restoration by water diffusion from the inter-crack space have similar effect, although probably much weaker.

42 ENGINEERING↗

Martensitic transformation induced strength-ductility synergy in additively manufactured maraging 250 steel by thermal history engineering

Maraging steels are known for their exceptional strength but suffer from limited work hardening and ductility. Here, in this study, we report an intermittent printing strategy to tailor the microstructure and mechanical properties of maraging 250 steel via tuning the thermal history during wire-arc directed energy deposition. By introducing a dwell time between adjacent layers, the maraging 250 steel is cooled below the martensite start temperature, triggering thermally-driven martensitic transformation during the printing process. Thermal cycling during subsequent layer deposition results in the formation of reverted austenite which shows a refined microstructure and induces elemental segregation between martensite and reverted austenite. The Ni enrichment in the austenite promotes stabilization of the reverted austenite upon cooling to room temperature. The reverted austenite is metastable during deformation, leading to strain-induced martensitic transformation under loading. Specifically, a 3 min interlayer dwell time produces a maraging 250 steel with approximately 8% reverted austenite, resulting in improved work hardening via martensitic transformation induced plasticity during deformation. Meanwhile, the higher cooling rate and refined prior austenite grains lead to substantially refined martensitic grains (by approximately fivefold) together with an increased dislocation density. With 3 min interlayer dwell time, the yield strength of the printed maraging 250 steel increases from 836 MPa to 990 MPa, and the uniform elongation is doubled from 3.2% to 6.5%. This intermittent deposition strategy demonstrates the potential to tune the microstructure of maraging steels for achieving strength-ductility synergy by engineering the thermal history during additive manufacturing.

Additive manufacturing↗

O'Hare Airport roadway traffic prediction via data fusion and Gaussian process regression

This study proposes an approach of leveraging information gathered from multiple traffic data sources at different resolutions to obtain approximate inference on the traffic distribution of Chicago's O'Hare Airport area. Specifically, it proposes the ingestion of traffic datasets at different resolutions to build spatiotemporal models for predicting the distribution of traffic volume on the road network. Due to its good adaptability and flexibility for spatiotemporal data, the Gaussian process (GP) regression was employed to provide short-term forecasts using data collected by loop detectors (sensors) and supplemented by telematics data. The GP regression is used to make predictions of the distribution of the proportion of sensor data traffic volume represented by the telematics data for each location of the sensors. Consequently, the fitted GP model can be used to determine the approximate traffic distribution for a testing location outside of the training points. Policymakers in the transportation sector can find the results of this work helpful for making informed decisions relating to current and future transportation conditions in the area.

42 ENGINEERING↗

FIB-ToF-SIMS characterization of irradiated U-10Zr

Post-irradiation examination (PIE) is critical for the performance assessment and qualification of nuclear fuels. Secondary ion mass spectrometry (SIMS) is a powerful materials characterization technique that allows for elemental and isotopic mapping with a depth resolution greater than EDS and EPMA. However, it has not yet been applied to PIE of metallic nuclear fuel. Here, in this work, we characterize an fast neutron spectrum irradiated U-10Zr fuel sample using a time-of-flight SIMS (ToF-SIMS) system connected to a FIB/SEM system, which allows for flexible sample analysis compared to a dedicated ToF-SIMS instrument. Analysis of the resulting hyperspectral micrograph data was aided by the development of an unsupervised machine learning (ML) algorithm that iterates on existing methods to segment the 3D micrographic datasets based on the similarity of mass spectra. The results showed that the FIB-ToF-SIMS instrument was potentially capable of spatially resolving closed fission gas bubbles in 3D by continued ion sputtering of the analyzed volume. Additionally, the ML algorithm proved useful in revealing the chemical segregation of light fission products (those with an atomic mass between approximately 85–105 amu, such as ruthenium and rhodium) plus matrix zirconium, heavy fission products (those with an atomic mass between approximately 135–150 amu, such as the lanthanides) and uranium. Future studies are planned to conduct FIB-ToF-SIMS analysis on more irradiated U-Zr samples to study the constituent redistribution.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Effect of stoichiometry and hydration level on water domain size and transport in poly(aryl piperidinium) alkaline anion-exchange membranes

Alkaline water electrolysis holds promise in decarbonizing the global economy by enabling renewable hydrogen production with non-precious group metal catalysts. Anion exchange membranes are an important component of alkaline water electrolyzers and would ideally be durable while allowing for high hydroxide conductivity. The poly(aryl piperidinium) (PAP) class of polymers has attracted recent interest due to their good mechanical robustness and high ionic conductivity. Here, in this work, we perform atomistic molecular dynamics (MD) simulations of several PAP polymers at experimentally relevant hydration levels and polymer ion exchange capacities (IECs) to gain nanoscale insight into their properties and to help elucidate the trade-offs that result from tuning the IECs through the polymer stoichiometry. Our MD-predicted macroscopic polymer properties were found to be in good agreement with experimentally available polymer swelling ratios, water-occupied volumes, X-ray scattering, and ionic conductivities. The models show that for hydration levels greater than 8H 2 O per cation a single water cluster will form that percolates through the system. The growth in water cluster size results in large polymer swelling, the creation of larger channels with widths of 7 Å or larger, and nanophase separation between the hydrophilic domains and the polymer with characteristic length scales of approximately 20–30 Å. The experimentally observed lack of a strong X-ray scattering peak at low wavevectors can be explained by a cancellation between the polymer-polymer/water-water and polymer-water correlations and not a loss in nanophase separation. The overlap in coordination environments of the hydroxide oxygen and polymer nitrogen atoms implies that vehicular diffusion between cationic groups could play a role in hydroxide transport. The polymers' hydroxide and water diffusion constants increase by approximately an order of magnitude between hydration levels of 8 and 20H 2 O per cation. However, there are diminishing returns in hydroxide diffusion constant once the IEC exceeds 2.4 meq/g.

08 HYDROGEN↗

Development of a high-throughput method for processing sponge-stick samples to detect viable Bacillus anthracis spores

Since the national validation of the sponge-stick based method for detection of Bacillus anthracis spores in environmental samples, there have not been focused efforts to address the low throughput nature of the method, which processes only one sample at one time. Sample processing remains a serious bottleneck for rapidly analyzing large numbers of samples expected from a biological warfare attack. Therefore, we developed a high-throughput method to simultaneously process multiple sponge-stick samples to be better prepared for rapid response and recovery after wide area anthrax incidents. In this method, sponges are placed in 50 mL tubes containing 25 mL extraction buffer and shaken to release spores, after which the suspension is recovered for analysis. Here, we determined that an additional extraction step, conducted in the same tubes with 10 mL buffer, further increased spore recovery from sponge-stick by approximately 10 %. We determined that orbital shaking and multi-tube vortexing were both more effective than reciprocating shaking for recovering spores. We conducted simultaneous processing of up to 12 sponge-stick samples and demonstrated comparable spore recovery efficiencies to the traditional low-throughput stomacher-based method (approximately 60 % recovery at 10 2 -spore level and 75 % recovery at 10 4 -spore level for both methods in three replicate experiments, P > 0.05 for two-tailed t-tests for each experiment and spore level). We also demonstrated that our high-throughput method could be integrated with Rapid Viability-Polymerase Chain Reaction (RV-PCR) analysis and could detect levels as low as 40 spores per sponge even when challenged by a PCR particulate contaminant.

Anthrax↗

Uncertainty guided online ensemble for non-stationary data streams in fusion science

Machine Learning (ML) is poised to play a pivotal role in the development and operation of next-generation fusion devices. Fusion data shows non-stationary behavior with distribution drifts, resulted by both experimental evolution and machine wear-and-tear. ML models assume stationary distribution and fail to maintain performance when encountered with such non-stationary data streams. Online learning techniques have been leveraged in other domains, however it has been largely unexplored for fusion applications. In this paper, we investigate online learning for continuous adaptation to drifting data streams in the prediction of Toroidal Field (TF) coils deflection at the DIII-D fusion facility. We further address the short-term performance degradation inherent to standard online learning, which arises because ground truth is unavailable at prediction time. To mitigate this issue, we propose an uncertainty-guided online ensemble framework. The method leverages the Deep Gaussian Process Approximation (DGPA) for calibrated uncertainty estimation and uses these uncertainty measures to guide a meta-algorithm that aggregates predictions from learners trained over different historical horizons. Our results show that online learning reduces prediction error by 80% compared to a static model. The online ensemble and the proposed uncertainty-guided ensemble further reduce error by approximately 6%, and 10% respectively, relative to standard single-model online learning, while also providing calibrated uncertainty estimates to support operational decision-making.

AI↗