Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallel”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mn(II)-induced phase transformation of Mn(IV) oxide in seawater

Manganese (Mn) oxides are key components of oceanic and lacustrine Mn nodules and influence metal cycling through oxidation and adsorption processes. Layered Mn oxides (LMOs) are the most common minerals in these nodules and the immediate products of microbially mediated Mn(II) oxidation by O 2 . LMOs can transform into tunneled Mn oxides (TMOs), Mn oxyhydroxides (MnOOH), or Mn(II,III) phase (Mn 3 O 4 ). LMOs often concur with Mn(II) in the environment and the adsorption and oxidation of Mn(II) by LMOs can greatly promote the transformation of LMOs to those phases. However, the Mn(II)-promoted transformation of LMOs in seawater—rich in various cations (300 mM Na + , 10 mM K + , 50 mM Ca 2+ , and 10 mM Mg 2+ ) remains poorly understood. We examined the transformation of δ-MnO 2 in artificial seawater (pH 8.2) under anoxic conditions with the Mn(II)/MnO 2 ratio (r) ranging from 0.08 to 3.83. To assess the effect of ionic strength (IS), parallel experiments were conducted in a mixed 530 mM NaCl and 10 mM KCl solution (having seawater ionic strength but without Ca 2+ and Mg 2+ ) and in 100 mM NaCl solution as a control. At low r (0.08), δ-MnO 2 transformed into triclinic birnessite and a 4 × 4 TMO in 100 mM NaCl solution, which, however, was suppressed in seawater due to strong interactions of Ca 2+ /Mg 2+ with δ-MnO 2 . In the mixed 530 mM NaCl and 10 mM KCl solution (the same ionic strength as of seawater), the transformation occurred extensively but the products had lower crystallinity compared to in 100 mM NaCl solution. At the high Mn(II)/MnO 2 ratios (0.5 ≤ r ≤ 3.83), δ-MnO 2 transformed extensively into MnOOH phases and hausmannite (Mn 3 O 4 ) in 100 mM NaCl solution. The seawater suppressed the transformation, but the suppression became weaker with increasing Mn(II)/MnO 2 ratio. For example, the transformation was completely suppressed at r = 0.5 but essentially negligible at r = 3.83. The suppression at these high Mn(II)/MnO 2 ratios was mainly ascribed to the influence of Ca 2+ and Mg 2+ rather than of the high IS, and the weaker suppression at the higher Mn(II)/MnO 2 ratio suggests stronger competition of Mn(II) with Ca 2+ /Mg 2+ for interacting with δ-MnO 2 . Moreover, the composition and crystallinity of the transformation products (i.e., the relative abundance of MnOOH (α, β, and γ) and Mn 3 O 4 ) were influenced by both the high ionic strength and the presence of Ca 2+ and Mg 2+ . Therefore, even though Ca 2+ and Mg 2+ concentrations are much lower than Na + in seawater, their impacts on the transformation are dominant. Our study explains why MnOOH phases, hausmannite, and TMOs are less common than LMOs in oceanic environments, partially because seawater chemistry suppresses their formation. LMOs are the most reactive for metal adsorption and oxidation among all Mn oxides. Thus, the high stability of LMOs in an oceanic environment confers the high impacts of Mn oxides on metal cycling in the ocean.

Divalent manganese↗
DOE OSTI2025

Molecular insights into Yb(III) speciation in sulfate-bearing hydrothermal fluids from X-ray absorption spectra informed by ab initio molecular dynamics

Rare earth elements (REEs) are critical for advanced technologies, yet in hydrothermal aqueous solutions the molecular level details of their interaction with ligands that control their geochemical transport and deposition remain poorly understood. Here, this study elucidates the coordination behavior of Yb 3+ in sulfate-rich hydrothermal fluids using in situ extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulations. By integrating multi-angle EXAFS with AIMD-derived constraints, we precisely resolve Yb 3+ coordination structures and ligand interactions under hydrothermal conditions. At room temperature, Yb 3+ is coordinated by five water molecules and two sulfate ligands (coordination number, CN = 8), forming a distorted square antiprism geometry. Increasing temperature induces progressive dehydration, reducing the hydration shell and favoring stronger sulfate complexation. At 200°C, sulfate ligands reorganize around Yb 3+ , shifting its geometry to a capped octahedron (CN = 7). At 300 °C, sulfate binding dominates, leading to structural reorganization that parallels the onset of sulfate mineral precipitation, consistent with the retrograde solubility of REE sulfates. These findings provide direct molecular-scale evidence that sulfate acts as both a transport and deposition ligand, critically influencing REE mobility in geochemical environments. Our results can also help to refine thermodynamic models of REE speciation in high-temperature hydrothermal fluids and improve our understanding of REE ore formation processes in nature.

AIMD↗
DOE OSTI2024

Sustainable aviation fuels from biomass and biowaste via bio- and chemo-catalytic conversion: Catalysis, process challenges, and opportunities

Sustainable aviation fuel (SAF) production from biomass and biowaste streams is an attractive option for decarbonizing the aviation sector, one of the most-difficult-to-electrify transportation sectors. Despite ongoing commercialization efforts using ASTM-certified pathways (e.g., lipid conversion, Fischer-Tropsch synthesis), production capacities are still inadequate due to limited feedstock supply and high production costs. New conversion technologies that utilize lignocellulosic feedstocks are needed to meet these challenges and satisfy the rapidly growing market. Combining bio- and chemo-catalytic approaches can leverage advantages from both methods, i.e., high product selectivity via biological conversion, and the capability to build C-C chains more efficiently via chemical catalysis. Herein, conversion routes, catalysis, and processes for such pathways are discussed, while key challenges and meaningful R&D opportunities are identified to guide future research activities in the space. Bio and chemo-catalytic conversion primarily utilize the carbohydrate fraction of lignocellulose, leaving lignin as a waste product. This makes lignin conversion to SAF critical in order to utilize whole biomass, thereby lowering overall production costs while maximizing carbon efficiencies. Thus, lignin valorization strategies are also reviewed herein with vital research areas identified, such as facile lignin depolymerization approaches, highly integrated conversion systems, novel process configurations, and catalysts for the selective cleavage of aryl C–O bonds. The potential efficiency improvements available via integrated conversion steps, such as combined biological and chemo-catalytic routes, along with the use of different parallel pathways, are identified as key to producing all components of a cost-effective, 100% SAF.

09 BIOMASS FUELS↗
DOE OSTI2024

New opportunities in produced water management: A market-based approach to produced water trading

Produced water (PW) is a byproduct of oil and gas (O&G) production. Obtained alongside the more valuable energy products, PW is usually characterized by high levels of salinity and often contains many contaminants (chemicals, soluble and insoluble oil, organics, etc.) making it unsuitable for release without substantial treatment. Couple this with the fact that PW is typically obtained at multiple times the rate of oil or gas, and the added transport, treatment, and disposal costs become a serious challenge for operators. These realities have led to ad-hoc practices including cooperation between industry competitors to recycle, share, or otherwise mitigate PW costs. The National Energy Technology Laboratory (NETL) in partnership with the Ground Water Protection Council (GWPC) is pursuing novel technology solutions to address PW issues that complement or improve ad-hoc practices adopted by operators. In this paper, we observe that well-established market management practices used in electrical power generation have natural analogues in the PW supply chain. These parallels open up a new line of research where we view PW management as a market equilibrium problem, and explore solutions that foster active and data-based collaboration among operators through market structures similar to power markets, with the ultimate objective of improving PW management costs and recycling rates. Here, we make a case for our observations, present a PW market clearing optimization model that shows how such a market system could operate in the O&G space, and provide an illustrative case study for demonstration.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Modeling CO 2 flow through faulted/fractured reservoirs using tEDFM in corner-point grids

The interest in underground CO 2 storage has increased significantly over the last decade because of the rising concern about global warming due to the growing levels of greenhouse gases in the atmosphere. Considering that CO 2 accounts for 80% of these greenhouse gases, carbon capture, utilization, and storage (CCUS) is regarded as one of the most direct approaches to achieving the net zero carbon target. Although CO 2 storage in deep saline aquifers and depleted gas reservoirs has been studied extensively, most studies use commercial simulators that model faults/fractures by simply modifying the transmissibility in the direction perpendicular to the fault surfaces. Here, this work shows that this simplistic approach ignores the accelerated flow in the directions parallel to the fault plane, leading to significantly higher leakage along the fault surface. To accurately model the flow of CO 2 in faulted reservoirs, we present the first transient embedded discrete fracture model for corner-point grids (tEDFM-CPG). By comparing the results of the tEDFM-CPG to high-resolution reference solutions, we show that this approach is accurate and efficient at predicting CO 2 flow in faulted/fractured reservoirs. Finally, this work presents the use of mixed reality (MR) to efficiently observe CO 2 gas migration in the interior of these corner-point grid systems.

25 ENERGY STORAGE↗
DOE OSTI2024

Monitoring spatiotemporal evolution of fractures during hydraulic stimulations at the first EGS collab testbed using anisotropic elastic-waveform inversion

The EGS Collab project acquired continuous active-source seismic monitoring (CASSM) data before, during, and after hydraulic stimulations at the first testbed at the depth of 4850 ft (1478 m) at the Sanford Underground Research Facility in Lead, South Dakota, for monitoring fracture creation and evolution. CASSM acquisition was conducted using 24 hydrophones, 18 accelerometers, and 17 piezoelectric sources within four fracture-parallel wells and two orthogonal wells. 3D anisotropic traveltime tomography and anisotropic elastic-waveform inversion of the campaign cross-borehole seismic data show that the rock within the stimulation region is a heterogeneous horizontal transverse isotropic medium. Here we use these inversion results as the initial models and apply 3D anisotropic first-arrival traveltime tomography and 3D anisotropic elastic-waveform inversion to the CASSM data acquired after each stimulation in May, 2018 and December, 2018. We observe the spatiotemporal evolution of seismic velocities and anisotropic parameters caused by hydraulic fracture stimulations, showing the regions of rock alternation caused by hydraulic fracture stimulation.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

An engineered lactate oxidase based electrochemical sensor for continuous detection of biomarker lactic acid in human sweat and serum

Lactate levels in humans reveal intensity and duration of exertion and provide a critical readout for the severity of life-threatening illnesses such as pediatric sepsis. Using the lactate oxidase enzyme (Lox) from Aerococcus viridians, we demonstrated its functionality for lactate electrochemical sensing in physiological fluids in a lab setting. The structure and dynamics of LOx were validated by crystallography, X-ray scattering, and hydroxyl radical protein footprinting. This provided a validated protein template for understanding and designing an enzyme-based electrochemical sensing elements. Using this template, LOx enzyme variants were generated and compared. Comparison of the variants demonstrates that one exhibits effective lactate sensing at significantly reduced operating voltages. Additionally, we demonstrate that the four hexahistidine-tags on each enzyme tetramer are sufficient for immobilization to create a durable, functional sensor, with no need for a covalent attachment, enabling self-immobilization and eliminating the need for additional immobilization steps. The functionality of the LOx enzyme variants was verified at physiological lactate concentrations in both human serum (0–4 mM) and artificial sweat (0–100 mM) using 3-electrode setups for analysis of the three variants in parallel. Accuracy of measurement in both artificial sweat and human serum were high. Employing a microfluidic flow cell, we successfully monitored varying lactate levels in physiological fluids continuously over a 2h period. Overall, this optimized LOx enzyme, which self-immobilizes onto gold sensing electrodes, facilitates efficient and reliable lactate detection and continuous monitoring at reduced operating voltages suitable for further development towards commercial use.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Performance evaluation of finned tube heat exchanger using curved wavy delta winglet vortex generators with circular perforations

Vortex generation is recognized as an effective passive approach to improve the heat transfer rate in fin and tube heat exchangers (FTHEs). The current study proposed innovative designs of curved wavy delta winglet vortex generators (CWDWVGs), both without and with circular perforations, to improve the heat transfer efficiency of FTHEs. There is potential to increase heat transfer performance further through various CWDWVG designs. Here, this study explores seven unique CWDWVG configurations, from 1-wave to 7-wave. A 3-D computational numerical model is utilized to evaluate the Thermo-hydraulic performance of FTHEs fitted with these different CWDWVG configurations across Reynolds numbers from 400 to 2000. This comparative analysis of the Thermo-hydraulic performance of FTHEs featuring four parallel circular tube layouts assesses configurations both with and without vortex generators (VGs) and various hole configurations. The evaluation of Thermo-hydraulic performance involves different parameters, including the London area goodness factor (LAGF), Colburn factor (j), friction factor (f), pressure drop (?P), and Nusselt number (Nu). Results demonstrate that the various CWDWVG configurations and the number of holes in them substantially affect the efficiency, as evaluated by the dimensionless Performance Evaluation Criteria (PEC). Notably, the 7-wave CWDWVGs surpassed other configurations, and integrating circular punched perforations further improved the thermal-hydraulic performance of FTHE. Specifically, the 7-wave CWDWVGs without holes demonstrated superior performance over other configurations, showing a significant increase in Nusselt number by 75.18% and 85.16% at Reynolds numbers of 2000 and 400, respectively, alongside an increase in pressure drop by 216.38% to 224.96%. Meanwhile, the 7-wave CWDWVGs with eight holes, in comparison to those without holes, exhibited a Nusselt number increase of 0.85%, a pressure drop decrease of 7.31%, and a reduction in the friction factor by 5.82%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

A high efficiency rooftop air conditioning system using multi-speed compressors

This study delineates a meticulous exploration of technologies to enhance the energy efficiency of rooftop air conditioning units, employing the DOE/ORNL heat pump design model for comprehensive engineering design and optimization. A baseline rooftop air conditioning unit, featuring a 13 ton (45.7 kW) cooling capacity and a 17.9 integrated energy efficiency ratio, served as the point of departure for substantive efficiency enhancements. Key modifications included the consolidation of two refrigerant circuits into one, integrating three parallel 2-stage (dual-speed) compressors, fan replacements with high-efficiency substitutes. Notably, a lower global warming potential refrigerant, R452B, was evaluated as a substitute for R-410A, demonstrating better performance in the lab prototype. Further, the achieved measured integrated energy efficiency ratio of 21.4 in the lab prototype surpassed the baseline integrated energy efficiency ratio. Comparative evaluations between R410A and R452B indicated heightened efficiency with the latter, showcasing a lab-demonstrated integrated energy efficiency ratio of 22.4 at the rated capacity of 13.8 ton (48.5 kW) and 23.9 integrated energy efficiency ratio at the rated capacity of 10 ton (35.2 kW). This research underscores the successful development of a rigorous, energy efficient rooftop air conditioning unit prototype with noteworthy environmental and economic implications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

The evolution of coal porosity during pyrolysis

Gasification of coal, municipal waste, or other organic materials is a potential hydrogen source that entails complex thermal decomposition and transport processes. This study provides a multiscale analysis of these processes for sub-bituminous (Usibelli, Healy, Alaska) and lignite (Center, North Dakota) coals and provides data useful for process design. The chemistry, mineralogy, and pore structures of pyrolyzed coal and their evolution with thermal decomposition are discussed. Samples pyrolyzed at 200–1000 °C were analyzed by small-angle neutron scattering; ultra-small, small-, and wide-angle X-ray scattering; and other complementary techniques. Scanning electron microscopy showed new pores in the high-temperature-pyrolyzed material. Upon heating, the coals became progressively denser, and the concentration of hydrogen decreased. Changes in pore volume fell into three temperature ranges: an initial, low-temperature range that, for the Usibelli coal, involved an increase in overall porosity; a mid-temperature range associated with pore volume loss; and a high-temperature range associated with significant porosity increase and char formation. This transformation was paralleled by changes in fractal dimension and correlation length. The higher the pyrolysis temperature the greater the small-pore-volume fraction and overall surface area became. Pyrolysis increased the lateral size of coal crystallites, decreased the amorphous fraction, and increased the aromatics fraction and overall coal rank. Comparisons of neutron and X-ray scattering data and subsequent water uptake studies showed that pre-dried coals can re-hydrate relatively rapidly upon exposure to air, which can significantly affect the porosity calculated from small-angle-scattering data. Fits to the cumulative porosity curves provide a method for modeling the physical and chemical transformation of hydrogen-containing feedstock during gasification.

Anovitz, Lawrence {Larry} [ORNL] (ORCID:0000000226↗
DOE OSTI2025

Cross-aldol condensation on missing linker sites of metal–organic framework UiO-66

Missing linkers induce coordinatively unsaturated Lewis acid-base pair sites in metal–organic frameworks such as UiO-66 (Zr). The Lewis acid-base pair is induced by replacing the benzene dicarboxylic acid linker with formic acid, which subsequently is thermally decomposed. The concentration of coordinatively unsaturated Zr 4+ cations was titrated with acetone. In parallel to their concentration, the catalytic activity increased while the selectivity for aldol condensation between furfural and acetone remained constant. The missing linkers generate a complex ensemble consisting of an (at least partly) accessible Zr 4+ cation surrounded by oxygen, catalyzing the cross-aldol condensation. As a result, density functional theory calculations, matching the experimental findings provide a more detailed understanding of the elementary steps during the condensation reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Physicochemical control of solvation and molecular assembly of charged amphiphilic oligomers at air-aqueous interfaces

Hypothesis : Understanding the rules that control the assembly of nanostructured soft materials at interfaces is central to many applications. We hypothesize that electrolytes can be used to alter the hydration shell of amphiphilic oligomers at the air-aqueous interface of Langmuir films, thereby providing a means to control the formation of emergent nanostructures. Experiments : Three representative salts – (NaF, NaCl, NaSCN) were studied for mediating the self-assembly of oligodimethylsiloxane methylimidazolium (ODMS-MIM + ) amphiphiles in Langmuir films. The effects of the different salts on the nanostructure assembly of these films were probed using vibrational sum frequency generation (SFG) spectroscopy and Langmuir trough techniques. Experimental data were supported by atomistic molecular dynamic simulations. Findings : Langmuir trough surface pressure – area isotherms suggested a surprising effect on oligomer assembly, whereby the presence of anions affects the stability of the interfacial layer irrespective of their surface propensities. In contrast, SFG results implied a strong anion effect that parallels the surface activity of anions. These seemingly contradictory trends are explained by anion driven tail dehydration resulting in increasingly heterogeneous systems with entangled ODMS tails and appreciable anion penetration into the complex interfacial layer comprised of headgroups, tails, and interfacial water molecules. These findings provide physical and chemical insight for tuning a wide range of interfacial assemblies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels

Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

42 ENGINEERING↗
DOE OSTI2024

A fully implicit, asymptotic-preserving, semi-Lagrangian algorithm for the time dependent anisotropic heat transport equation

In this paper, we extend the operator-split asymptotic-preserving, semi-Lagrangian algorithm for time dependent anisotropic heat transport equation proposed in Chacón et al. (2014) [18] to use a fully implicit time integration with backward differentiation formulas. The proposed implicit method can deal with arbitrary heat-transport anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ $\ggg$ 1 (with $\mathcal{X}$∥, $ \mathcal{X}$⟂ the parallel and perpendicular heat diffusivities, respectively) in complicated magnetic field topologies in an accurate and efficient manner. Further, the implicit algorithm is second-order accurate temporally and demonstrates an accurate treatment at boundary layers (e.g., island separatrices), which was not ensured by the operator-split implementation. The condition number of the resulting algebraic system is independent of the anisotropy ratio, and is inverted with preconditioned GMRES. We propose a simple preconditioner that renders the finite-dimensional linear operator compact, resulting in mesh-independent convergence rates for topologically simple magnetic fields, and convergence rates scaling as ~ (NΔt) 1/4 (with N the total mesh size and Δt the timestep) in topologically complex magnetic-field configurations. We demonstrate the accuracy and performance of the approach with test problems of varying complexity, including an analytically tractable boundary-layer problem in a straight magnetic field, and a topologically complex magnetic field featuring magnetic islands with extreme anisotropy ratios $\mathcal{X}$∥ /$ \mathcal{X}$⟂ = 10 10 ) .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Hardware acceleration for HPS algorithms in two and three dimensions

We provide a flexible, open-source framework for hardware acceleration, namely massively-parallel execution on general-purpose graphics processing units (GPUs), applied to the hierarchical Poincaré–Steklov (HPS) family of algorithms for building fast direct solvers for linear elliptic partial differential equations. To take full advantage of the power of hardware acceleration, we propose two variants of HPS algorithms to improve performance on two- and three-dimensional problems. In the two-dimensional setting, we introduce a novel recomputation strategy that minimizes costly data transfers to and from the GPU; in three dimensions, we modify and extend the adaptive discretization technique of Geldermans and Gillman [1] to greatly reduce peak memory usage. We provide an open-source implementation of these methods written in JAX, a high-level accelerated linear algebra package, which allows for the first integration of a high-order fast direct solver with automatic differentiation tools. We conclude with extensive numerical examples showing our methods are fast and accurate on two- and three-dimensional problems.

Fast direct solvers↗