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.

481 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Revisiting the Causes of the US Midwest Great Flood of 1993

The floods that occurred in the Midwest during the summer of 1993 remain one of the greatest flooding events in US history. While a considerable amount of work has already been done addressing its causes there still remain uncertainties as to the relative roles of forcing from the tropics and the extra-tropics as well as the extended nature of the excessive rainfall that lasted throughout much of the late spring and summer of that year. Here we use a novel “replay” technique in which the GEOS AGCM is constrained to remain close to the MERRA-2 reanalysis over specified regions of the globe to isolate the causes of the excessive rain. We show the main rainfall event in the upper Midwest in July (extending into August) was the result of the development of a large amplitude Rossby wave (and associated pronounced moisture influx from the Gulf of Mexico) apparently forced by anomalous transient eddies in the western North Pacific. Rainfall that fell in the Midwest during the two months prior to that, while not as intense, acted to precondition the soil, with forcing from the tropics contributing to the May precipitation, while the excessive Midwest rain that fell during June was part of a record hemispheric-wide precipitation signal that appears to be linked to an unusually cold north polar stratosphere. The nature of the stratospheric impact and the possible delayed influence of the Pinatubo eruption is currently being assessed.

Midwest↗
NASA NTRS2021

Lila: Optimal Dispatching in Probabilistic Temporal Networks using Monte Carlo Tree Search

Executing a Probabilistic Simple Temporal Network (PSTN) amounts at scheduling, i.e. \textit{dispatch}, a set of events under time uncertainty. This constitutes a NP-hard online optimization problem. The right execution time must be dynamically assigned to each event of the PSTN such that the temporal constraints are met, whereas activity durations are progressively observed as the execution unfolds. We propose a dispatching algorithm based on Monte Carlo Tree Search, called Lila, with the following characteristics: (i) it is an anytime algorithm, both offline and online, proven asymptotically optimal; (ii) it returns the current probability of success, either before or at any moment during operations; (iii) it handles any possible continuous or discrete, even non-parametric, probability distributions, as well as inter-dependencies between random variables, exogenous and endogenous uncertainty; and (iv) can be easily extended to handle probabilistic external events, PSTNs with resources, PSTNs with cutoff times and precondition chains, etc. Lila is universal in the sense that it can handle any dispatching protocol, simply by specifying it to the algorithm. It has the unlimited flexibility offered by the simulation paradigm, whilst it asymptotically converges to optimal decisions and/or robustness approximations.

Chien, Steve A.↗
NASA NTRS2022

Understanding Structure in Line-Driven Stellar Winds Using Ultraviolet Spectropolarimetry in the Time Domain

The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, and preconditions the surrounding medium for all following events, with consequences for ISM energization, chemical enrichment, and dust formation. Understanding these processes requires accurate observational constraints on the mass-loss rates of the most luminous stars, which can also be used to test theories of stellar wind driving. In the past, mass-loss rates have been characterized via collisional emission processes such as optical Hα and free-free radio emission, but these so-called “density squared” diagnostics require correction in the presence of widespread clumping. Recent observational and theoretical evidence points to the likelihood of a ubiquitously high level of such clumping in hot-star winds, but quantifying its effects requires a deeper understanding of the complex dynamics of radiatively driven winds and their stochastic instabilities. Furthermore, large-scale structures initiating in surface anisotropies and propagating throughout the wind can also affect wind driving and alter mass-loss diagnostics. Time series spectroscopy of high resonance-line opacity in the UV, capable of high resolution and high signal-to-noise, are required to better understand these complex dynamics, and more accurately determine mass-loss rates. The proposed Polstar mission (Scowen et al. 2022, this volume) provides the necessary resolution at the Sobolev (∼10 km s −1 ) or sound-speed (∼20 km s −1 ) scale, for over three dozen bright galactic massive stars with signal-to noise an order of magnitude above that of the celebrated MEGA campaign (Massa et al. 1995) of the International Ultraviolet Explorer (IUE), via continuous observations that track propagating structures through the winds in real time. Supporting geometric constraints are provided by the polarimetric capabilities present in all the datasets of such a mission.

Polstar – NASA MIDEX↗
NASA NTRS2022

The Polstar High Resolution Spectropolarimetry MIDEX Mission

The Polstar mission will provide a space-borne 60 cm spectropolarimeter operating at ultraviolet (UV) wavelengths, capturing all four Stokes parameters (intensity, two linear polarization components, and circular polarization). Polstar’s capabilities are designed to meet its goal of determining how circumstellar gas flows alter and inform massive star evolution, affect the stellar remnant population, and stir and enrich the interstellar medium (ISM). These will be achieved by investigating the dynamical geometries in the winds and disks of hot stars, the composition and magnetic alignment of interstellar dust, and the star-forming accretion disks of UV-bright stars at an important transition boundary. Together these areas map out a kind of two-way interface between massive stars and their effect on our galaxy, wherein the stellar winds enrich the ISM with metals and kinetic energy, preconditioning their environment and the stellar endpoints prior to undergoing supernova. The ISM dust in turn reveals the composition and magnetic environment leading to new star formation, and the accretion disks of Herbig Ae/Be stars reveal how the ISM gas returns to make new massive stars. Polstar will combine high-resolution spectroscopy in the time domain with high-precision UV polarimetry. Doppler-shifted UV resonance line opacity will provide information about circumstellar kinematics, while polarization gives complementary geometric information about unseen structures. The composition and magnetic alignment of the smallest interstellar dust grains provides a probe of the ISM utilizing radiative alignment theory (RAT). Polstar will operate in the far-UV (FUV) at 122–200 nm at high spectral resolution of around R ∼ 30k, and at FUV and near-UV (NUV) wavelengths of 122–320 nm at lower spectral resolutions of 0.1 - 1k. Detection of polarization levels as weak as 0.1% are expected, with a temporal cadence ranging from 5–10 minutes for most wind variability studies, to hours or days for sampling rotation, to days or weeks for sampling binary orbits, to months to a year for sampling substructure in the inner regions of protoplanetary disks. Sub-meter-class aperture is well suited to access this wide array of time domain science, made possible by restricting to a few hundred bright, massive stars, necessarily extincted by a small to moderate column of interstellar dust, informing both the attributes of the stars and the ISM through which they are seen. As such, the focus is on our own galaxy and its evolutionary drivers, but a few targets in the Magellanic clouds offer the potential to extend this understanding to low-metallicity environments.

Paul A Scowen↗
NASA NTRS2023

Recurrent 3he-Rich Solar Energetic Particle Injections Observed By Solar Orbiter at ∼0.5 Au

We report Solar Orbiter observations of six recurrent solar energetic particle injections in 2022 March 3−6 at ∼0.5 au. All but one were associated with jets emanating from a plage near a large sunspot in active region 12 957. We saw large jets in injections with high 3He and Fe enrichments and minor jets in injections with no or lower enrichments. Furthermore, the event with the highest enrichment showed a more compact configuration of the underlying photospheric magnetic field. The higher fluences as well as harder spectra were seen in the event with a simultaneous jet and wider eruption. However, in this case, the energy buildup time in the source might be required to produce such spectra. Extreme-ultraviolet images from Solar Orbiter revealed a number of intersecting loops at the base of jets not seen from 1 au that might be a precondition for the recurrent events.

Sun: particle emission↗
NASA NTRSDate not supplied

On the Prediction of Surface Melt Over the Greenland Ice Sheet at NWP and S2S Timescales

Greenland surface melt “events” – encompassing large swaths of the ice sheet over approximately 4-10 day periods – have occurred with increasing frequency in the 21st Century. Surface processes are significant to the ice sheet’s overall contribution to sea level. In addition, the sudden onset of widespread surface melt produces hazardous conditions for researchers and for communities adjacent to the ice sheet. The mechanisms producing melt events, their associated large-scale forcing, and timescales of interest are highly variable and poorly understood. Here we examine the predictability of melt events in a numerical weather prediction system (NWP) and a subseasonal-to-seasonal system (S2S). The NASA Global Modeling and Assimilation Office forward processing (GMAO-FP) is an NWP system that incorporates an evolving, global hybrid 4D ensemble–variational (4DEnVar) data assimilation system and a state-of-the-art model that is run twice-daily with forecasts out to five and ten days. The GMAO S2S system is a hindcast-anomaly, coupled-model prediction system initialized with an atmospheric analysis, and in situ and satellite ocean data assimilation. Both systems represent ice sheet surface hydrology and snowpack processes. We examine the performance of the systems over the period 2018-2023, a span that includes approximately nine significant events, against analyses and satellite-derived melt extent. In addition to extent, we examine related conditions and their role in the melt forecast, including the onset of both omega and Rex atmospheric blocking patterns, sea surface temperatures, sea ice extent, water vapor transport, precipitation, and surface preconditioning. We further examine the tendency for false-alarm events in each system and assess reasons for the forecast performance. North Atlantic atmospheric blocking patterns leading to surface melt are well forecast by the NWP system out to around eight days, while the S2S system generally underestimates the magnitude of such events at monthly lead times. These systems suggest that while short term forecasts of ice sheet melt and melt events are relatively robust, S2S predictions of surface conditions are highly dependent on the ability of a coupled system to generate the short-term atmospheric conditions associated with surface melt events.

Richard Cullather↗
NASA NTRS2024

Predictability of the Minimum Sea Ice Extent from Winter Fram Strait Ice Area Export: Model vs Observations

Observations show predictive skill of the minimum sea ice extent (Min SIE) from late winter anomalous off-shore ice drift along the Eurasian coastline – leading to local ice thickness anomalies at the onset of the melt season – a signal then amplified by the ice-albedo feedback. We assess whether the observed seasonal predictability of September sea ice extent (Sept SIE) from Fram Strait Ice Area Export (FSIAE, a proxy for Eurasian coastal divergence) is present in Global Climate Model (GCM) large ensembles, namely the CESM2-LE, GISS-E2.1-G, FLOR-LE, CNRM-CM6-1, and CanESM5. All models show distinct periods where winter FSIAE anomalies are negatively correlated with the May sea ice thickness (May SIT) anomalies along the Eurasian coastline, and the following Sept Arctic SIE, as in observations. Counter-intuitively, several models show occasional periods where winter FSIAE anomalies are positively correlated with the following Sept SIE anomalies when the mean ice thickness is large or late in the simulation when the sea ice is thin, and/or when internal variability increases. More importantly, periods with weak correlation between FSIAE and the Sept SIE dominate, suggesting that summer melt processes dominate over late winter preconditioning and May SIT anomalies. In general, we find that the coupling between the FSIAE and ice thickness anomalies along the Eurasian coastline is an ubiquitous feature of GCMs and that the relationship with the following Sept SIE is dependent on the mean Arctic sea ice thickness.

off-shore ice drift↗
DOE OSTI2025

Two-level overlapping additive Schwarz preconditioner for training scientific machine learning applications

In this work we introduce a novel two-level overlapping additive Schwarz preconditioner for accelerating the training of scientific machine learning applications. The design of the proposed preconditioner is motivated by the nonlinear two-level overlapping additive Schwarz preconditioner. The neural network parameters are decomposed into groups (subdomains) with overlapping regions. In addition, the network’s feed-forward structure is indirectly imposed through a novel subdomain-wise synchronization strategy and a coarse-level training step. Through a series of numerical experiments, which consider physicsinformed neural networks and operator learning approaches, we demonstrate that the proposed two-level preconditioner significantly speeds up the convergence of the standard (LBFGS) optimizer while also yielding more accurate machine learning models. Moreover, the devised preconditioner is designed to take advantage of model-parallel computations, which can further reduce the training time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Performance and Energy Study of GPU-Resident Preconditioners for Conjugate Gradient Solvers: In the Context of Existing and Novel Approaches

Optimizing a particular subprogram out of the set of Basic (sparse) Linear Algebra Subprograms (BLAS) for a given architecture is a common topic of research. In applications, however, these BLAS functions rarely appear in isolation; usually, many of them are used together, in various combinations and with varying inputs. As the need to solve a large, sparse linear system is ubiquitous throughout HPC applications, linear solvers constitute a realistic, sufficiently complex and well-defined representative use case for composite BLAS routines. To this end, based on a representative set of matrices drawn from a diverse set of fields, we present a framework to study, from the performance and energy perspective, the efficacy of GPU- resident parallel Conjugate Gradient (CG) linear solver with different preconditioner options, including Gauss-Seidel, Jacobi, and incomplete Cholesky. We also propose a novel GPU-based preconditioner, in which the triangular solves are approximated by an iterative process. The development of this preconditioner was motivated by solving large graph Laplacian linear systems, for which the existing preconditioners either perform slow on GPU-based platforms or are not applicable. We compare the performance of these preconditioners on different hardware accelerator architectures, i.e., AMD MI250X, MI100, Nvidia A100, V100, and Jetson. Our experiments reveal performance trade-offs and provide information on how to select the best strategy for the given linear system, dictated by its properties, and the platform of interest. We demonstrate the application of our novel preconditioner for solving CG and graph Laplacian systems. Overall, the framework can be utilized as a benchmark to guide informed decisions in choosing a specific preconditioner, i.e., whether it is better to rely on the performance of a triangular solver or on the performance of sparse matrix-vector product. Finally, by considering power consumption to solve the linear systems, we report the energy footprint for the solvers.

Preconditioned Conjugate Gradient, GPUs, iterative↗
DOE OSTI2024

Constrained Local Approximate Ideal Restriction for Advection-Diffusion Problems

Herein this paper focuses on developing a reduction-based algebraic multigrid (AMG) method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based AMG approach, $\ell \text{AIR}$ (local approximate ideal restriction), that was developed for solving advection-dominated problems. Though this new solver is very effective in the advection-dominated regime, its performance degrades in cases where diffusion becomes dominant. This is consistent with the fact that in general, reduction-based AMG methods tend to suffer from growth in complexity and/or convergence rates as the problem size is increased, especially for diffusion-dominated problems in two or three dimensions. Motivated by the success of $\ell \text{AIR}$ in the advective regime, our aim in this paper is to generalize the AIR framework with the goal of improving the performance of the solver in diffusion-dominated regimes. To do so, we propose a novel way to combine mode constraints as used commonly in energy-minimization AMG methods with the local approximation of ideal operators used in $\ell \text{AIR}$. The resulting constrained $\ell \text{AIR}$ algorithm is able to achieve fast scalable convergence on advective and diffusive problems. In addition, it is able to achieve standard low complexity hierarchies in the diffusive regime through aggressive coarsening, something that was previously difficult for reduction-based methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Algebraic Multigrid with Filtering: An Efficient Preconditioner for Interior Point Methods in Large-Scale Contact Mechanics Optimization

Large-scale contact mechanics simulations are crucial in many engineering fields such as structural design and manufacturing. In the frictionless case, contact can be modeled by minimizing an energy functional; however, these problems are often nonlinear, nonconvex, and increasingly difficult to solve as mesh resolution increases. In this work, we employ a Newton-based interior-point (IP) filter line-search method, an effective approach for large-scale constrained optimization. While this method converges rapidly, each iteration requires solving a large saddle-point linear system that becomes ill-conditioned as the optimization process converges, largely due to IP treatment of the contact constraints. Such ill-conditioning can hinder solver scalability and increase iteration counts with mesh refinement. Here, to address this, we introduce a novel preconditioner, algebraic multigrid with filtering (AMGF), tailored to the Schur complement of the saddle-point system. Building on the classical AMG solver, commonly used for elasticity, we augment it with a specialized subspace correction that filters near null space components introduced by contact interface constraints. Through theoretical analysis and numerical experiments on a range of linear and nonlinear contact problems, we demonstrate that the proposed solver achieves mesh independent convergence and maintains robustness against the ill-conditioning that notoriously plagues IP methods. These results indicate that AMGF makes contact mechanics simulations more tractable and broadens the applicability of Newton-based IP methods in challenging engineering scenarios. More broadly, AMGF is well suited for problems, optimization or otherwise, where solver performance is limited by a low-dimensional subspace, such as those arising from localized constraints, interface conditions, or model heterogeneities. This makes the method widely applicable beyond contact mechanics and constrained optimization.

Mathematics and Computing↗
DOE OSTI2025

Station Impact Analysis 2025

As part of the U.S. DOE EVs@Scale consortium, the NextGen Profiles (NGP) project presents analysis and results from the study of High Power Charging Electric Vehicles and Battery Charging Infrastructure. High Power Charging equipment is capable of recharging electric vehicle traction batteries at power levels of 200KW and above. The intent of the project is to further understand the most recent technological capabilities of the electric mobility industry related to charging performance. The project aims to develop EV, EVSE, and Fleet characterization testing practices and comprehensive analysis with inputs from key industry stakeholders. The results published in this NextGen Profiles project report provide data and insight for use by numerous entities including modeling and simulation organizations, policy makers, fleet planners, industry stakeholders and the general public involved with the development, deployment and operation of electrified transportation technologies. The factors influencing Electric Vehicle (EV) Direct Current Fast Charging (DCFC), including EV battery specifications, temperature effects on lithium-ion battery and power electronics performance, lithium-ion battery SOC bounding and charging station design considerations are specifically investigated to analyze their impacts on charging station operation and recommendations are made to minimize charge station dwell time, reduce charging costs and mitigate electric grid and charge station congestion. Additional high-power charging results are anticipated in future publications in support of the U.S. DOE EVs@Scale consortium NextGen Profiles project.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Hybridized Discontinuous Galerkin Methods for Computational Fluid Dynamics

Hybridizable Discontinuous Galerkin (HDG) methods hold promise for any applications with significant advection character, including thermal hydraulics in light water reactors and advanced reactor concepts and fluid models of plasmas in magnetic confinement fusion. Its features include natural upwinding, local element conservation, and extensibility to arbitrarily high order accuracy. In the last fiscal year we have implemented HDG in the Multiphysics Object-Oriented Simulation Environment (MOOSE). We developed a first-of-its-kind automatic static condensation system in MOOSE’s underlying finite element library libMesh which can condense out arbitrarily many internal variables. Finally, we developed the first preconditioner for HDG discretizations of the Navier-Stokes equations which shows robust performance across a wide range of problem sizes and Reynolds numbers. This preconditioner yields solution times that are equivalent to the fastest developed for industry standard finite volume methods. Moreover, the arbitrarily high-order nature of HDG makes it a prime candidate for acceleration via graphical processing units (GPUs). We believe these developments will hold significant importance for future DOE Nuclear Energy (NE) and Fusion Energy Science (FES) programs.

97 MATHEMATICS AND COMPUTING↗