Search NASA⌕ Search

SEARCH · Search NASA

Results for “Initialization method”

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 253 records · Page 14

Accelerating Multiphase Simulations With Denoising Diffusion Model Driven Initializations

This study introduces a hybrid fluid simulation approach that integrates generative diffusion models with physics‐based simulations, aiming at reducing the computational costs of flow simulations while still honoring all the physical properties of interest. Pore‐scale simulations enhance our understanding of applications such as assessing hydrogen and storage efficiency in underground reservoirs. Nevertheless, they are computationally expensive and the presence of non‐unique solutions can require multiple simulations within a single geometry. To overcome the computational cost hurdle, we propose a method that couples generative diffusion models and physics‐based simulations. While training the data‐driven model, we simultaneously generate initial conditions and perform physics‐based simulations using these. This integrated approach enables us to receive real‐time feedback on a single compute node equipped with both CPUs and GPUs. By efficiently managing these processes within a single compute node, we can continuously monitor performance and halt training once the model meets the specified criteria. To test our model, we generate realizations in a real Berea sandstone fracture which shows that our technique is up to 4.4 times faster than commonly used flow simulation initializations.

36 MATERIALS SCIENCE↗

A Safe Response to Renewable Energy Hazards

The International Association of Fire Fighters (IAFF) and Underwriters Laboratories, LLC (UL) in conjunction with UL Solutions initiated a joint project in 2022 under an agreement with the United States Department of Energy-Office of Energy Efficiency and Renewable Energy (DOE-EERE). This project focused on two separate and important initiatives related to energy efficiency in residential buildings. Initiative 1: Fire Performance on Energy Efficient Exterior Walls. Initiative 2: Firefighting Tactics in Residential Properties with Building Energy Storage Systems (BESS). The project’s first initiative addresses concerns surrounding new technologies with enhanced, energy-efficient exterior walls installed on residential properties. The concerns of fire rapidly traveling vertically up the exterior of these walls were examined. This addressed a growing concern from the first responder community that many times, the fires on the exterior of residential buildings have already evolved into an attic fire by the time of arrival – making it problematic to address the fire scenario. The test plan for Initiative 1 incorporated a modified version of an American Society for Testing and Materials (ASTM) test method, ASTM E2707, Standard Test Method for Determining Fire Penetration of Exterior Wall Assemblies Using a Direct Flame Impingement Exposure, as the foundation of the research. The test method involved a wall structure intended to represent retrofit construction to evaluate how fire would spread vertically or laterally. The second aspect of the UL-IAFF Project focuses on the fire service response to Residential Battery Energy Storage System (RBESS) incidents. These simulation tests were constructed in the large-scale fire test facility at UL Solutions’ Northbrook, IL campus. A baseline test was conducted that involved a test structure with no batteries—shelving units populated with standardized commodities, representing a typical garage with cellulosic and plastic contents. Three additional tests have been conducted to generate data with the contribution of energy storage system (ESS) batteries to compare fire and explosion hazards against the baseline test. Through this work, fire service tactical considerations can be explored. From the data, the team can determine 1) the visual indicators of a residential fire that has involved an RBESS (or, potentially, other large batteries) and 2) the impact of fire service-initiated ventilation of the structure on the fire conditions and explosion risks.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Predicting Initial Trans-Membrane Pressure for Optimized Operations in UF Unit Using Random Forest

With the growing scarcity of freshwater, innovative process design mechanisms like Ultra-filtration(UF) units are increasingly gaining attention among water treatment utilities to address the rising demand. Ensuring reliable water production necessitates efficient resource utilization, minimizing downtime in UF systems. Recent advancements in machine learning (ML) have enabled the development of accurate data-driven models for Model Predictive Control (MPC), often requiring minimal prior knowledge of underlying physical processes. In this study, we present predictive regression models based on Random Forest (RF) and Auto-Regressive (AR) approaches to forecast the initial Trans-Membrane Pressure (TMP) for each filtration cycle in data generated by Direct Potable Reuse (DPR) systems. The proposed RF-based model demonstrates superior performance compared to baseline methods, including historical mean, Last Observation Carried Forward (LOCF), and naïve AR models, across various forecasting horizons in terms of root mean square (RMSE) metric. Accurate prediction of initial TMP is critical for optimizing CCRO operations, as it enables the development of robust modelling frameworks that enhance process efficiency and reliability. The demonstrated efficacy of the RF-based approach highlights its potential as a tool for real-time decision-making in water treatment systems, paving the way for advanced process optimization and sustainable water resource management.

Mukherjee, Subrata [ORNL] (ORCID:0000000309930338)↗

Modeling Analysis of Ball-Milling Process for Battery-Electrode Synthesis

The mechanical alloying process is a promising method for synthesizing electrode materials for batteries owing to its benefits such as the ability to produce nanostructured, high-performing electrode alloys, no adverse effects on the solid electrolyte for solid-state batteries, stable production of thick electrodes, simple processing steps, and low processing costs. It is gaining intensive attention in the battery industry as one of the best methods to replace the conventional wet-slurry-solvent method, and its application is rapidly increasing these days. However, the operation is currently conducted purely based on trial-and-error methods without fully utilizing the features of its functions. Here, this may be attributed to a lack of understanding of the effect of operating parameters on the alloying process and final products. Surprisingly, there is a scarcity of the literature conducting fundamental research to comprehend the underlying physics of the entire mechanical alloying process, resulting in a significant knowledge gap. To address this knowledge gap, extensive research was conducted. The existing literature on mechanical alloying was reviewed to comprehend the current state of understanding and to discuss the direction for future research. Mathematical expressions were developed to create physics-based models capable of capturing the entire mechanical alloying process, including milling kinetics and defect-enhanced phase evolution. These methods were then applied to investigate the impact of operating parameters such as milling frequency, initial mole ratio of the alloyed materials, density of grinding balls, and energy required for the powders to become amorphous (i.e., the amorphization energy threshold). This research aimed not only to comprehend the direct effects of these operating parameters but also to unveil the physics underlying the ball-milling process. The results of our study can serve as crucial information for the battery industry in designing or operating the ball-milling process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A novel high frequency, high bandwidth, three phase Mach–Zehnder optical data link

In inertial confinement fusion, hydrogen isotopes are fused together under high pressures and temperatures. Typically, the duration of these experiments is incredibly short, on the order of around 60–150 ps. Due to the high radiation environment, a detector’s signal is typically data linked far distances to a protected location. The diagnostic challenge for fusion reaction history measurements is to measure signals of interest maximizing dynamic range while also maintaining time resolution on the order of 10 ps. In this work, we present a new experimental optical data link used to efficiently transport diagnostic signals over great distances to the recording system while not restricting the dynamic range. The concept of a three phase Mach–Zehnder modulator system is introduced as well as a description of the physical prototype. The initial results from testing at the OMEGA facility show that this system is a viable method for signal transportation.

47 OTHER INSTRUMENTATION↗

Equilibrium-gated pattern formation: How molecular dissociation thermodynamics drive emergent behavior in dissipative polymeric systems

Emergent patterns in biological systems arise through dissipative processes that balance reaction and transport phenomena, producing highly functional properties from self-regulating mechanisms. Synthetic fabrication, by contrast, often relies on user-controlled, multistep methods that lack the self-organizing capabilities of natural systems. Inspired by nature, we sought chemical systems that integrate strongly coupled reaction and transport phenomena, identifying frontal ring-opening metathesis polymerization (FROMP) as a method capable of creating diverse forms and functions through reactive processing. By employing discrete molecular initiators, FROMP allows precise control of key reaction steps—inhibition, initiation, and propagation. Using an integrated computational and experimental framework, we uncover how near-equilibrium inhibition dynamics, coupled with far-from-equilibrium reaction kinetics, drive pattern formation in frontally polymerized synthetic materials. We propose the concept of equilibrium-gated pattern formation, demonstrating how initiator chemistry can be tuned to achieve programmable macroscale properties. Our study reveals a surprising insight: Emergent behavior in FROMP systems arises from the inhibition-dominated regime of resin composition, expanding prior observations that such behavior is confined to a narrow compositional space near the boundary between front quenching and uniform front propagation. We identify a broader compositional window, far from the quenching regime, where emergent behavior reliably manifests. This expanded design space significantly enhances the operational flexibility of reactive systems and their capacity for self-organization. Furthermore, these insights provide a roadmap for designing bioinspired materials with self-organizing capabilities, unlocking possibilities in synthetic manufacturing.

Architected Polymers↗

Rotational excitation and de-excitation of magnesium mono-sulphide (MgS) by He collisions

ABSTRACT Magnesium mono-sulphide (MgS) plays a crucial role in astrochemical processes, particularly in the interstellar medium (ISM), where metal-sulphur chemistry influences molecular formation and evolution. This work presents a detailed study of the rotational excitation and de-excitation of MgS through collisions with helium (He) atoms, the second most abundant species in the ISM after hydrogen, which includes both atomic (H) and molecular forms (H2). The focus on MgS–He collisions arises from He's high abundance, chemical inertness, and simpler electronic structure, which make it well suited for quantum scattering calculations. These characteristics establish He as an ideal candidate for initial studies, providing fundamental data for future investigations involving H2. The study uses quantum scattering methods to calculate the collisional rate coefficients over a broad temperature range. These rates are critical for interpreting observational data on MgS and predicting its abundance in space. The interaction potential between MgS and He is calculated using the rigid rotor approximation and the Jacobi coordinate system, employing the CCSD(T)-F12a/aug-cc-pVTZ method for accurate two-dimensional potential energy surface. The study explores the anisotropic nature of the MgS–He interaction, which favours odd Δj rotational transitions at low collision energies. The inelastic cross-sections for rotational transitions involving up to 16 rotational levels of MgS were computed up to 1000 cm−1, enabling the calculation of rate coefficients up to 150 K for Δj = ±1, ±2, and ± 3 rotational transitions. The results show that Δj = 1 transitions dominate at low temperatures, while Δj = 2 transitions become more significant at higher temperatures. This study provides valuable data for interpreting future astrophysical observations of MgS. The findings also propose new rotational transitions for MgS detection in space, enhancing our ability to track and study this molecule in various cosmic environments.

Hendaoui, Hamza (ORCID:0000000218641872)↗

Beta-Decay Half-Lives beyond 54 Ca: A Systematic Survey of Decay Properties Approaching the Neutron Dripline

In an experiment performed at the Facility for Rare Isotope Beams (FRIB) using the FRIB Decay Station initiator, 15 new half-lives of isotopes near 54 Ca were measured. A new method of extracting lifetimes from experimental data, taking into account the unknown 𝛽-delayed neutron emission branches of very neutron-rich nuclei, was developed to enable systematic uncertainty analysis. The experiment observed a dramatic change in the half-life systematics for the isotopes with neutron number 𝑁 = 34. Beyond 𝑁 = 34, the decline of nuclear lifetime is much slower, leading to longer than anticipated lifetimes for near-dripline nuclei. State-of-the-art shell-model calculations can explain the experimental results, revealing the imprint of shell effects and the need for modification of single-particle neutron states. In conclusion, the results from a newly developed quasiparticle random phase approximation model with potential for making global predictions were also tested against the experimental results and good agreement was found.

Physics - Nuclear physics and radiation physics↗

Planar Systems for Quantum Information

This project aims to develop two‐dimensional (2D) moiré materials as a quantum simulator to implement model Hamiltonians and their phase diagrams. Progress in quantum information science (QIS) requires the development of advanced quantum materials systems. The rich family of layered van der Waals materials and their heterostructures present opportunities to create previously unrealized types of applications for QIS. Specifically, when two layers of van der Waals materials are overlaid with a small twist angle or/and lattice mismatch, a moiré superlattice with a period of about ten nanometers is formed. This provides a periodic trapping potential for electrons. Electrons can tunnel between the traps and repel each other by their mutual Coulomb interactions. The platform of 2D moiré materials provides many attractive features, including tunability of length and energy scales, charge density, and even lattice symmetry. It presents new possibilities for realizing quantum simulation of the many-body physics in a solid-state platform. This integrated team of six investigators seeks to develop relevant theoretical treatments to link ab-initio studies of 2D moiré materials to model Hamiltonians and to evaluate correlated phases predicted by these model Hamiltonians in the relevant regimes. On the experimental side, the team aims to develop methods to realize a homogeneous and highly controlled potential landscape for the electrons and to initiate, protect, and measure their quantum many-body states.

36 MATERIALS SCIENCE↗

Physics-driven Explosions of Stripped High-mass Stars: Synthetic Light Curves and Spectra of Stripped-envelope Supernovae with Broad Light Curves

Stripped-envelope supernovae (SESNe) represent a significant fraction of core-collapse supernovae, arising from massive stars that have shed their hydrogen and, in some cases, helium envelopes. The origins and explosion mechanisms of SESNe remain a topic of active investigation. In this work, we employ radiative-transfer simulations to model the light curves and spectra of a set of explosions of single, solar-metallicity, massive Wolf–Rayet stars with ejecta masses ranging from 4 to 11 M ⊙ , which were computed from a turbulence-aided and neutrino-driven explosion mechanism. We analyze these synthetic observables to explore the impact of varying ejecta mass and helium content on observable features. We find that the light curve shape of these progenitors with high ejecta masses is consistent with observed SESNe with broad light curves but not the peak luminosities. The commonly used analytic formula based on rising bolometric light curves overestimates the ejecta mass of these high-initial-mass progenitor explosions by a factor of up to 2.6. In contrast, the calibrated method by Haynie et al., which relies on late-time decay tails, reduces uncertainties to an average of 20% within the calibrated ejecta mass range. Spectroscopically, the He I 1.083 μm line remains prominent even in models with as little as 0.02 M ⊙ of helium. However, the strength of the optical He I lines is not directly proportional to the helium mass but instead depends on a complex interplay of factors such as the 56 Ni distribution, composition, and radiation field. Thus, producing realistic helium features requires detailed radiative transfer simulations for each new hydrodynamic model.

79 ASTRONOMY AND ASTROPHYSICS↗

Novel Synthesis and Characterization of Uranium-Zirconium Carbonitride by Direct Casting

UZrCN has exhibited thermophysical properties beneficial to high temperature reactor applications such as nuclear thermal rockets. The present work investigates a novel liquid phase synthesis method involving the admission of nitrogen gas during arc melting of uranium, zirconium, and carbon. Initial microstructural examinations using SEM indicated that the samples remain heterogenous with zirconium-rich cores in a uranium-rich matrix. Heterogeneity resulted from large differences in melting temperature of the major constituents and rapid solidification. Additional EDS, combustion analysis, and inert gas fusion analysis proved that the core regions have both uranium and zirconium and establish that carbon and nitrogen are retained during the fabrication process. XRD clarifies that the light elements did not form compounds with the uranium matrix but rather incorporated into the core region forming a sub-stoichiometric UZrCN. Heat treatment performed on equimolar U-Zr-C resulted in an increase in the homogenous phase present with an XRD pattern reflecting a ternary UZrCN despite small amounts of segregation in the final button.

36 - MATERIALS SCIENCE↗

On Uranium-Zirconium Carbonitride Formation by Direct Casting - A Novel Synthesis Study

UZrCN has exhibited thermophysical properties beneficial to high temperature reactor applications such as nuclear thermal rockets. The present work investigates a novel liquid phase synthesis method involving the admission of nitrogen gas during arc melting of uranium, zirconium, and carbon. Initial microstructural examinations using SEM indicated that the samples remain heterogenous with zirconium-rich cores in a uranium-rich matrix. Heterogeneity resulted from large differences in melting temperature of the major constituents and rapid solidification. Additional EDS, combustion analysis, and inert gas fusion analysis proved that the core regions have both uranium and zirconium and establish that carbon and nitrogen are retained during the fabrication process. XRD clarifies that the light elements did not form compounds with the uranium matrix but rather incorporated into the core region forming a sub-stoichiometric UZrCN. Heat treatment performed on equimolar U-Zr-C resulted in an increase in the homogenous phase present with an XRD pattern reflecting a ternary UZrCN despite small amounts of segregation in the final button.

36 - MATERIALS SCIENCE↗

Research Progress and Perspectives on Pre‐Sodiation Strategies for Sodium‐Ion Batteries

Sodium‐ion batteries (SIBs) with abundant elements have garnered significant attention from researches as a promise compensation to lithium‐ion batteries (LIBs). However, the large‐scale commercial application of SIBs is partially hindered by the limited initial coulombic efficiency (ICE) due to the irreversible formation of solid electrolyte interphase (SEI) and intercalation into the defects in the anode. Similar to pre‐lithiation techniques, pre‐sodiation approaches are considered to be one of the most direct and effective way to compensate for the loss of active sodium at the anode side of SIBs during the initial cycle. In this context, additional sodium ions are pre‐injected to the cathode/anode material by chemical/electrochemical methods, aiming to improve battery span life and energy density. Here, this review delves into the necessity and impact of pre‐sodiation techniques, compiling the latest research progress, for instance, self‐sacrificing cathode additives, over‐sodiated cathode materials, direct contact and solution chemical pre‐sodiation. Notably, the research mechanisms underlying solution chemical pre‐sodiation are highlighted. This comprehensive overview aims to foster a deeper understanding of the pre‐sodiation techniques and expects to provide guidance for realizing the commercial application of high energy density sodium‐ion batteries.

25 ENERGY STORAGE↗

Rapid Prototyping for Nanoparticle-Based Photonic Crystal Fiber Sensors

The advent of nanotechnology has motivated a revolution in the development of miniaturized sensors. Such sensors can be used for radiation detection, temperature sensing, radio-frequency sensing, strain sensing, and more. At the nanoscale, integrating the materials of interest into sensing platforms can be a common issue. One promising platform is photonic crystal fibers, which can draw in optically sensitive nanoparticles or have its optical properties changed by specialized nanomaterials. However, testing these sensors at scale is limited by the the need for specialized equipment to integrate these photonic crystal fibers into optical fiber systems. Having a method to enable rapid prototyping of new nanoparticle-based sensors in photonic crystal fibers would open up the field to a wider range of laboratories that could not have initially studied these materials in such a way before. This manuscript discusses the improved processes for cleaving, drawing, and rapidly integrating nanoparticle-based photonic crystal fibers into optical system setups. The method proposed in this manuscript achieved the following innovations: cleaving at a quality needed for nanoparticle integration could be done more reliably (≈100% acceptable cleaving yield versus ≈50% conventionally), nanoparticles could be drawn at scale through photonic crystal fibers in a safe manner (a method to draw multiple photonic crystal fibers at scale versus one fiber at a time), and the new photonic crystal fiber mount was able to be finely adjusted when increasing the optical coupling before inserting it into an optical system (before, expensive fusion splicing was the only other method).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

3D mesh regularization within an ALE code using a weighted line sweeping method

The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.

42 ENGINEERING↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

Sustained advances in the mathematics of modeling and simulation have resulted in the capability today for routine simulation of a number of large scale complex DOE-relevant systems. As remarkable as this capability for solving the so-called forward problem is, it is typically only the first step-an inner loop within an outer loop that explores the simulation model's parameter space and decision space to characterize uncertainty in the model's predictions, learn unknown model parameters from data, design the most informative experiments, determine optimal control strategies, and create optimal designs. Broadly, what unifies all of these outer loop problems is that they are, in one form or another, optimization problems over parameter/control/design space that are constrained by complex uncertain models. To fully realize the power of scientific simulation as a basis for scientific discovery, technological innovation, and rational decision-making, it is imperative to move beyond simulation to tackle the outer loop of optimization for learning from data, experimental design, and control with complex uncertain models. When the models under consideration are large-scale and complex, and when the optimization variable and uncertain parameter spaces are high (or infinite) dimensional, this constitutes a grand challenge of the highest order, and is intractable with conventional methods. To overcome these challenges, the AEOLUS Center was established to develop a unified mathematical, computational, and statistical framework for (1) Learning predictive models from complex data via Bayesian inference and optimization, and (2) Optimizing experiments, processes, and designs using the resulting uncertain models. These problems are intractable with conventional methods, for several reasons: (1) The simulation problems that govern the inner loops of the optimization problems are expensive to execute (due to severe nonlinearity, heterogeneity, multiphysics/multiscale coupling); (2) The optimization variable and uncertain parameter spaces are high dimensional, often stemming from discretizations of infinite dimensional fields such as initial conditions, sources, or material properties. We argue that the key to overcoming these challenges is to develop new mathematical, computational, and statistical methods that exploit the structure of the Bayesian inference and optimization problems mediated by their underlying complex uncertain models. This structure includes the regularity, sparsity, geometry, low intrinsic dimensionality, and multifidelity nature of the maps from uncertain parameter/optimization variable spaces to the specific objectives targeted: Bayesian inference, optimal experimental design, and optimal control design. Black box methods developed as generic tools are incapable of exploiting this structure. To be successful, we must create, integrate, and cross-fertilize ideas across multiple areas of applied math--including approximation theory, Bayesian inference, data science, experimental design, information theory, machine learning, model reduction, optimal control theory, parallel algorithms, PDE-constrained optimization, randomized algorithms, stochastic optimization, and uncertainty quantification--all while exploiting the structure of the problems at hand. With this goal in mind, we have marshaled a team of leading authorities in these areas. While the methods we develop will be broadly applicable across a wide spectrum of DOE problems in which experiments inform models and the systems those models describe must be optimized under uncertainty, we have chosen a specific area, advanced manufacturing and materials, to drive our work. AMM is characterized by complex models across multiple scales, and is a rich source of challenging problems in inference, experimental design, and optimal control, requiring multifaceted and integrated advances in applied mathematics. As such, AMM serves as an excellent vehicle to motivate and demonstrate the advances in applied mathematics developed by our center.

97 MATHEMATICS AND COMPUTING↗

Accelerating Hamiltonian Monte Carlo for Bayesian inference in neural networks and neural operators

Hamiltonian Monte Carlo (HMC) is a powerful and accurate method to sample from the posterior distribution in Bayesian inference. However, HMC techniques are computationally demanding for Bayesian neural networks due to the high dimensionality of the network’s parameter space and the non-convexity of their posterior distributions. Therefore, various approximation techniques, such as variational inference (VI) or stochastic gradient MCMC, are often employed to infer the posterior distribution of the network parameters. Such approximations introduce inaccuracies in the inferred distributions, resulting in unreliable uncertainty estimates. In this work, we propose a hybrid approach that combines inexpensive VI and accurate HMC methods to efficiently and accurately quantify uncertainties in neural networks and neural operators. The proposed approach leverages an initial VI training on the full network. We examine the influence of individual parameters on the prediction uncertainty, which shows that a large proportion of the parameters do not contribute substantially to uncertainty in the network predictions. This information is then used to significantly reduce the dimension of the parameter space, and HMC is performed only for the subset of network parameters that strongly influence prediction uncertainties. This yields a framework for accelerating the full batch HMC for posterior inference in neural networks. We demonstrate the efficiency and accuracy of the proposed framework on deep neural networks and operator networks, showing that inference can be performed for large networks with tens to hundreds of thousands of parameters. Finally, we show that this method can effectively learn surrogates for complex physical systems by modeling the operator that maps from upstream conditions to wall-pressure data on a cone in hypersonic flow.

Bayesian inference↗

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↗