Search NASA⌕ Search

SEARCH · Search NASA

Results for “Preconditioning”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Comparison of the Effects of Bipolar Membrane Preparation Conditions on the Mechanical Durability and Electrochemical Performance for Electrodialysis Applications

Bipolar membranes (BPMs) are enabling materials for electrochemical conversion technologies such as water electrolysis, fuel cells, CO 2 electrolysis, and electrodialysis (ED) for direct air/ocean capture of CO 2 . However, current BPM durability can suffer from chemical, mechanical, and performance degradation when operated at high current density (ion flux) and physical scale. Therefore, this limits its adoption in a wider applications space. BPMs have several known degradation mechanisms, including chemical breakdown of ion-exchange polymers, loss of junction adhesion, or physical breakdown due to shearing force and pressure swings in an electrodialysis cell. To assess the electrochemical stability and mechanical durability of BPMs under operational conditions, we investigated how fabrication conditions (including preconditioning, hot-pressing temperature and pressure, and catalyst loading) impact the adhesion of custom-made BPMs. T-peel studies were performed ex situ to quantify adhesive forces of BPMs, and bipolar membrane electrodialysis (BPMED) experiments were performed to assess the electrochemical performance of the corresponding BPMs. The results of this systematic comparison indicate that hydration and heated pressing create improved adhesion during the fabrication of BPMs, and BPMED testing shows that these fabrication techniques are not detrimental to the electrochemical performance of the BPMs.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Minimizing Interfacial Resistance between Polymer Electrolytes and Metal Electrodes Using Applied Current

Reducing the interfacial resistance between different phases in electrochemical systems is crucial for enabling practical applications. In this work, we proposed a process for reducing the interfacial resistance between polymer electrolytes and metal electrodes. Thus far in the literature, the lowest interfacial resistance reported in these systems is 15 Ω·cm2. In this study, assembled and preconditioned symmetric cells with lithium–indium alloy electrodes showed similar values. The current through the cell was increased in steps up to the limiting current. This resulted in a permanent decrease of the interfacial resistance to values as low as 1 Ω·cm2, a value that is comparable to that of optimized lithium-ion batteries. The proposed process is general, and it could be applied to any combination of polymer electrolytes and metal electrodes.

Lee, Jaeyong↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2025

Durability Research Is Pivotal for Perovskite Photovoltaics

Metal halide perovskite solar cells have shown promising power conversion efficiencies, but commercialization requires that decent durability is also demonstrated. Under normal operation, solar cells are subject to a complex combination of stressors, such as visible light, ultraviolet light, heat, humidity, mechanical stress and electric potential, which complicates the understanding of failure mechanisms. Existing stress tests do not act as a time machine. In new materials systems such as perovskite photovoltaics, the tests have no known relationship to field service. In this Perspective we recommend following a durability learning cycle that interleaves photovoltaic module engineering with field testing; accelerated testing; and preconditioning and performance engineering. We advocate for field testing to demonstrate real-world performance and identify field-relevant failure modes, and urge the community to develop accelerated and qualification tests that account for device metastability, variations in material composition and different/various processing methods. In conclusion, these practices are more difficult, but more important, than the simple pursuit of higher initial efficiencies.

14 SOLAR ENERGY↗
DOE OSTI2024

Radiation image reconstruction and uncertainty quantification using a Gaussian process prior

We propose a complete framework for Bayesian image reconstruction and uncertainty quantification based on a Gaussian process prior (GPP) to overcome limitations of maximum likelihood expectation maximization (ML-EM) image reconstruction algorithm. The prior distribution is constructed with a zero-mean Gaussian process (GP) with a choice of a covariance function, and a link function is used to map the Gaussian process to an image. Unlike many other maximum a posteriori approaches, our method offers highly interpretable hyperparamters that are selected automatically with the empirical Bayes method. Furthermore, the GP covariance function can be modified to incorporate a priori structural priors, enabling multi-modality imaging or contextual data fusion. Lastly, we illustrate that our approach lends itself to Bayesian uncertainty quantification techniques, such as the preconditioned Crank–Nicolson method and the Laplace approximation. The proposed framework is general and can be employed in most radiation image reconstruction problems, and we demonstrate it with simulated free-moving single detector radiation source imaging scenarios. We compare the reconstruction results from GPP and ML-EM, and show that the proposed method can significantly improve the image quality over ML-EM, all the while providing greater understanding of the source distribution via the uncertainty quantification capability. Furthermore, significant improvement of the image quality by incorporating a structural prior is illustrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

A physically interpretable precursor framework for sub-seasonal prediction of Northern Hemisphere flash flourishing

Flash flourishing describes rapid vegetation increases that can quickly reshape land–atmosphere exchanges and impacts on ecosystem, yet its large-scale precursors, circulation context, and sub-seasonal predictability remain poorly understood. Here, we identified onset-stage circulation regimes across northern extratropical latitudes (NEL; >30°N) using 200 and 1000 hPa geopotential height, and examined their regional expressions over eastern Asia, western North America, and Europe. Flash flourishing onset in East Asian was associated with a baroclinic circulation regime and was preceded by a North Atlantic sea surface temperature (SST) precursor at a four-pentad lead. In contrast, onset in western North American and European preferentially occurred under barotropic regimes, preconditioned by Great Plains soil moisture at three-pentad lead and North Atlantic SST at a four-pentad lead, respectively. Ridge regression forecasts revealed regime-dependent sub-seasonal predictability, with mean out-of-sample R 2 exceeding 0.3 up to lead times of two pentads in East Asia, three pentads in western North America, and four pentads in Europe. Together, these findings established a mechanistic and regionally specific framework for anticipating rapid vegetation greening at sub-seasonal timescales.

Kong, Xiangxu [Nanjing Univ. of Information Scienc↗
DOE OSTI2024

Solving the Hele–Shaw flow using the Harrow–Hassidim–Lloyd algorithm on superconducting devices: A study of efficiency and challenges

The development of quantum processors for practical fluid flow problems is a promising yet distant goal. Recent advances in quantum linear solvers have highlighted their potential for classical fluid dynamics. In this study, we evaluate the Harrow–Hassidim–Lloyd (HHL) quantum linear systems algorithm (QLSA) for solving the idealized Hele–Shaw flow. Our focus is on the accuracy and computational cost of the HHL solver, which we find to be sensitive to the condition number, scaling exponentially with problem size. This emphasizes the need for preconditioning to enhance the practical use of QLSAs in fluid flow applications. Moreover, we perform shots-based simulations on quantum simulators and test the HHL solver on superconducting quantum devices, where noise, large circuit depths, and gate errors limit performance. Error suppression and mitigation techniques improve accuracy, suggesting that such fluid flow problems can benchmark noise mitigation efforts. Finally, our findings provide a foundation for future, more complex application of QLSAs in fluid flow simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Implicit full-F simulations of neoclassical ion transport

The development of implicit time integration capabilities for axisymmetric full-F continuum simulations of ion neoclassical transport is reported. The approach involves the implicit treatment of the gyrokinetic Vlasov equation coupled to the nonlinear Fokker–Planck collision model in the long-wavelength limit approximation. To facilitate implicit simulations, advanced preconditioning of individual physics operators is developed, and a global multi-physics preconditioner is constructed by adopting an operator splitting methodology. The algorithm is implemented in the finite-volume code COGENT and is applied to study neoclassical transport properties for both the main ion species and the lithium impurity species in the closed-field-line region of the LTX- β tokamak. The implicit COGENT simulations elucidate the role of non-local transport effects, while demonstrating substantial speedup over the corresponding explicit approach.

Dorf, Mikhail [Lawrence Livermore National Laborat↗
DOE OSTI2025

Semi-implicit continuum kinetic modeling of weakly collisional parallel transport in a magnetic mirror

We present implicit-explicit (IMEX) kinetic simulations of weakly collisional parallel plasma transport in magnetic mirror configurations using the continuum code COGENT. The numerical scheme employs a Jacobian-free Newton–Krylov method with algebraic multigrid preconditioning to overcome the severe time step limitations imposed by strong mirror forces in fully explicit schemes. Applied to parameters relevant to the Wisconsin HTS Axisymmetric Mirror experiment, the IMEX approach enables time steps up to 2.5×10 4 times larger than those permitted by explicit methods, resulting in a 2500× speedup in 1D–2V simulations of parallel transport with kinetic ions and Boltzmann electrons. Additionally, a reduced bounce-averaged model for a square mirror is implemented to support the computationally intensive fully kinetic simulations. The bounce-averaged formulation is used to evaluate the numerical convergence of the velocity-space discretization algorithms and to assess the role of the collision model by comparing simulations employing the nonlinear Fokker–Planck and the simplified Lenard–Bernstein–Dougherty collision operators.

Collision theories↗
DOE OSTI2025

Increased frequency of planetary wave resonance events over the past half-century

We demonstrate a tripling in the frequency of planetary wave resonance events over the past halfcentury, coinciding with the rise in persistent boreal summer weather extremes. This increase aligns with changes in the underlying climate conditions favoring these events, including amplified Arctic warming and land-sea thermal contrast. We also observe increased prevalence of resonant amplification events following the mature phase of strong El Niño events, suggesting that such events may precondition the mean state conditions in ways that favor large-scale quasi-stationary wave patterns and quasi-resonant wave amplification. Since the impact of anthropogenic warming on quasi-resonant amplification is not well captured by current-generation climate models, it is likely that models are underpredicting the potential increase, indicating even greater risk of persistent extreme summer weather events with ongoing warming.

Arctic amplification↗
DOE OSTI2026

Routes to high-performance operation in Wendelstein 7-X: turbulence suppression with shaping of the density profile

Steep density gradients generally lead to improved plasma performance in the neoclassically optimized stellarator W7-X. This is evident in the global energy confinement time as well as in the ion temperature and can be explained by a strong reduction of the ion temperature gradient turbulence. Such conditions can experimentally be realized by several methods: injection of cryogenic hydrogen pellets, appropriate combination of neutral beam and electron cyclotron resonance heating (ECRH) and, in some cases, with low power of ECRH after preconditioning of the first wall. The duration of the improved phases is determined by the ability to sustain the steep density gradient, by technical limitations of the involved systems and, eventually, by the plasma stability. This paper gives an overview of relevant experimental results and presents example discharges where the improved confinement conditions could be extended to multiple seconds: up to 4 s using neutral beam injection and from 14 to 40 s with steady state pellet injection. In these plasmas the turbulent thermal diffusivity is reduced by a factor of 3 to 4 in a broad radial range, which allows high ion temperatures of up to 3 keV at the densities of about 1.5 • 10 20 m −3 .

high performance plasma↗
DOE OSTI2025

Unprecedented Beaufort Sea ice loss in late summer 2021 and its relationship to an extended period of unusually stormy weather

Previous case studies have linked cyclone-induced atmospheric forcing and/or upper-ocean processes to notable Arctic sea ice loss events in the summers of 2012 and 2016. This study examines a more recent and noteworthy case in late summer 2021 in which substantial sea ice loss followed a period of surface meteorological extremes in the Beaufort Sea region of the Arctic. We focus on the period from mid-August to mid-September 2021 that coincided with the Office of Naval Research THINICE Pilot Field Campaign and investigate stormy and windy conditions with respect to air-sea processes impacting sea ice conditions. We find that during the stormy first half of the campaign, cyclone-induced energy fluxes into the marginal ice zone and surrounding waters preconditioned the ice pack for more rapid melt later in the campaign. The second half of the campaign, in contrast, was marked by non-cyclone wind events that enhanced turbulent (namely sensible) heat fluxes into the ice and upper ocean that increased melt. Moreover, this latter period had enhanced advection of the Beaufort Sea ice pack into above-freezing waters, increasing bottom melt to >1 cm d −1 over the remainder of the campaign. While findings are shown to vary by surface type and at relatively small (i.e. ice-floe) scales, insights are offered on the roles of late summer coupled processes on rapid ice loss events in today’s Arctic environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗
DOE OSTI2026

Artificial Soiling Replication of Field Losses on Commercial Photovoltaic Modules

Here, this study demonstrates the capabilities of an indoor artificial soiling approach developed to closely replicate the natural, cyclic soil accumulation processes in the field-dust suspension, deposition, and sedimentation/cementation-for a subtropical climate. In this work, a near-field environment is replicated in an artificial soiling cubic chamber through controlled regulation of humidity, temperature, dust type, and dust concentration, based on site-specific historical climate data. Two different models (MA and MB) of full-size commercial photovoltaic modules from a single manufacturer, installed side by side in the mid-Atlantic United States, were retrieved and subjected to artificial soiling experiments and various characterization measurements, including short-circuit current, colorimetry, reflectance, X-ray fluorescence, laser diffraction, and optical microscopy. Both in the field and in our improved field-representative artificial soiling tests, the MA model experienced roughly twice the soiling loss as the MB model. To closely replicate the field soiling losses for a site-specific climate, it is critical to include: 1) The use of field-collected dust with identical dust chemistry and particle distribution instead of standardized ISO 12103 Arizona Road dusts, 2) the use of only a small amount of field-collected dust inside the chamber during the deposition process (e.g., 0.15 g), and 3) the preconditioning of the surface coating for the partial/full dose of UV stress as experienced in the field during sunlight exposure and the abrasion as experienced in the field during regular module cleaning activities, if/as needed. The field-representative artificial soiling method developed here could potentially be adopted for rank ordering of various antisoiling coatings developed by researchers and industry stakeholders.

14 SOLAR ENERGY↗
DOE OSTI2025

Surrogate-Based Autotuning for Randomized Sketching Algorithms in Regression Problems

Algorithms from Randomized Numerical Linear Algebra (RandNLA) are known to be effective in handling high-dimensional computational problems, providing high-quality empirical performance as well as strong probabilistic guarantees. However, their practical application is complicated by the fact that the user needs to set various algorithm-specific tuning parameters which are different from those used in traditional NLA. This paper demonstrates how a surrogate-based autotuning approach can be used to address fundamental problems of parameter selection in RandNLA algorithms. In particular, we provide a detailed investigation of surrogate-based autotuning for sketch-and-precondition (SAP)-based randomized least squares methods, which have been one of the great success stories in modern RandNLA. Empirical results show that our surrogate-based autotuning approach can achieve near-optimal performance with much less tuning cost than a random search (up to about 7.6x fewer trials of different parameter configurations). Moreover, while our experiments focus on least squares, our results demonstrate a general-purpose autotuning pipeline applicable to any kind of RandNLA algorithm.

Cho, Younghyun↗
DOE OSTI2025

Parallel-in-Time Solution of Scalar Nonlinear Conservation Laws

Here, we consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially nonoscillatory (WENO) reconstructions, and in time with high-order explicit Runge–Kutta methods. The solution of the global, discretized space-time problem is sought via a nonlinear iteration that uses a novel linearization strategy in cases of nondifferentiable equations. Under certain choices of discretization and algorithmic parameters, the nonlinear iteration coincides with Newton’s method, although, more generally, it is a preconditioned residual correction scheme. At each nonlinear iteration, the linearized problem takes the form of a certain discretization of a linear conservation law over the space-time domain in question. An approximate parallel-in-time solution of the linearized problem is computed with a single multigrid reduction-in-time (MGRIT) iteration; however, any other effective parallel-in-time method could be used in its place. The MGRIT iteration employs a novel coarse-grid operator that is a modified conservative semi-Lagrangian discretization and generalizes those we have developed previously for nonconservative scalar linear hyperbolic problems. Numerical tests are performed for the inviscid Burgers and Buckley–Leverett equations. For many test problems, the solver converges in just a handful of iterations with a convergence rate independent of mesh resolution, including problems with (interacting) shocks and rarefactions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Zero-Order Reaction Kinetics v. 3.6

Zero-RK is a software package that simulates chemically reacting systems using sparse, preconditioned, adaptive matrix methods to achieve orders-of-magnitude reduction in simulation time while maintaining accurate results.

McNenly, MatthewJ [Lawrence Livermore National Lab↗
DOE OSTI2026

NCCS High Performance GMRES Mixed Precision

HPG-MxP is a software package that performs a fixed number of multigrid preconditioned (using a Gauss-Seidel smoother) Generalized minimal residual (PGMRES) iterations in order to solve a possibly nonsymmetric large sparse linear system of equations. It is designed to be a benchmark to measure a computer's performance for sparse linear algebra workloads typical in scientific computing while allowing the use of mixed precision methods. The solution is required to have convergence characteristics and accuracy similar to double precision GMRES. It is based on the High Performance Conjugate Gradient Benchmark (HPCG) which restricts all implementations to use only the IEEE double precision format (FP64). The original implementation (https://github.com/hpg-mxp/hpg-mxp) was written by Ichitaro Yamazaki, Jennifer Loe, Christian Glusa, Sivasankaran Rajamanickam, Piotr Luszczek, and Jack Dongarra. Please refer to that repository for documentation on the original implementation. This version is maintained by the National Center for Computational Sciences at Oak Ridge National Laboratory. It is highly scalable and optimized for Oak Ridge Leadership Computing Facility (OLCF) systems, particularly Frontier.

Kashi, Aditya [Oak Ridge National Laboratory (ORNL↗