Search NASA⌕ Search

SEARCH · Search NASA

Results for “IDEAL FLUID”

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

High-order cumulants and correlation functions near the critical point from molecular dynamics

We present a systematic investigation of particle-number fluctuations in the crossover region near the critical end point of a first-order phase transition using molecular dynamics simulations of the classical Lennard-Jones fluid. We extend our prior studies to third- and fourth-order cumulants in both coordinate- and momentum-space acceptances and integrated correlation functions (factorial cumulants). We find that, even near the critical point, non-Gaussian cumulants equilibrate on timescales comparable to those of the second-order cumulants, but show stronger finite-size effects. The presence of interactions and of the critical point leads to strong deviations of the cumulants from the ideal-gas baseline in coordinate space; these deviations are expected to persist in momentum space in the presence of collective expansion. In particular, the kurtosis becomes strongly negative, κσ 2 ≪ − 1 , on the crossover side of the critical point. However, this signal is significantly diluted once an efficiency cut used to distinguish protons from baryons is applied, leading to |κσ 2 | ≲ 1 even in the presence of the critical point. We discuss our results in the context of ongoing measurements of proton-number cumulants in heavy-ion collisions in RHIC-BES-II.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Simulating many-engine spacecraft: Exceeding 1 quadrillion degrees of freedom via information geometric regularization

We present an optimized implementation of the recently proposed information geometric regularization (IGR) for unprecedented scale simulation of compressible fluid flows applied to multi-engine spacecraft boosters. We improve upon state-of-the-art computational fluid dynamics (CFD) techniques in terms of computational cost, memory footprint, and energy-to-solution metrics. Unified memory on coupled CPU–GPU or APU platforms increases problem size with negligible overhead. Mixed half/single-precision storage and computation are used on well-conditioned numerics. We simulate flow at 200 trillion grid points and 1 quadrillion degrees of freedom, exceeding the current record by a factor of 20. A factor of 4 wall-time speedup is achieved over optimized baselines. Ideal weak scaling is observed on OLCF Frontier, LLNL El Capitan, and CSCS Alps using the full systems. Strong scaling is near ideal at extreme conditions, including 80% efficiency on CSCS Alps with an 8 node baseline and stretching to the full system.

Wilfong, Benjamin [Georgia Institute of Technology↗
DOE OSTI2025

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression↗
DOE OSTI2024

Generative Physics-Informed Neural Network Solving Multi-Scale and Multi-Phase Plasma Chemical Flow Field

Low-temperature plasmas (LTPs) are non-equilibrium systems with near-room-temperature gas and highly energetic electrons. This makes them ideal for delicate applications in biomedicine and semiconductor manufacturing, enabling processes like wound healing, sterilization, etching, and plasma-enhanced chemical vapor deposition without thermal damage. However, LTPs involve complex chemistries, with hundreds of species and thousands of reactions, complicating their diagnosis, prediction, and control. Conventional diagnostics, such as Fourier-transform infrared spectroscopy (FTIR), laser-induced fluorescence (LIF), and optical emission spectroscopy (OES), offer limited species detection, while mass spectrometry (MS) struggles with low-sensitivity species. Additionally, LTP simulations face multi-scale challenges, as macroscopic fluid dynamics and microscopic particle collisions operate on vastly different timescales. To address these issues, we developed an artificial intelligence (AI) based diagnostic system: a generative physics-informed neural network (PINN-Gen) that can predict spatially resolved species concentrations and temperatures in LTPs by integrating experimental data from planar LIF with microscopic plasma chemical kinetics and macroscopic fluid mechanics, including plasma-liquid interactions at the interface between two phases. PINN-Gen solves no equations but checks the errors of physical laws by substituting the output from neural network, and the comparison with the experimental results. Thus, it naturally avoids the multi-scale difficulty of numerical simulations and predicts the results of conventionally unsolvable multi-scale and multi-phase problems. The real-time prediction will be robust due to the physical information used in the training of such a neural network, and only very limited input of condition required due to its generative feature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

One-dimensional temperature measurement of supersonic jet flow by resonantly ionized photoemission thermometry of molecular nitrogen

As the field of fluid dynamics progresses, the demand for sophisticated diagnostic methods to accurately assess flow conditions rises. In this work, resonantly ionized photoemission thermometry (RIPT) has been used to directly target and ionize diatomic nitrogen (N 2 ) to measure one-dimensional (1D) temperature profiles in a supersonic jet flow. This technique can be considered non-intrusive as the premise uses resonantly enhanced multiphoton ionization (REMPI) to target molecular nitrogen. This resonance excites N 2 into absorption bands of the P, Q, and R rotational branches of N 2 ( b 1 Π u ). The ideal (3 + 1) REMPI scheme excites from the ground state and ionizes N 2 ( b 1 Π u ← X 1 Σ g + ) where de-excitation results in photoemission from the first negative band of ionized N 2 + ( B 2 Σ u + → X 2 Σ g + ) as nitrogen returns to the ground state. The resulting emission can be observed using an intensified camera, thus permitting inference of the rotational temperature of ground-state molecular nitrogen. A linearly regressive Boltzmann distribution is applied based on previous calibration data for this technique to quantify the temperature along the ionized line. This work applies this technique to a pure N 2 supersonic jet in cross-flow and counter-flow orientations to demonstrate N 2 RIPT’s applications in a supersonic flow. Temperature variations are observed at different locations downstream of the exit in cross-flow, and axisymmetric in counter-flow, to generate profiles characterizing the flow dynamics. Due to the collisional effects resulting from the number density of N 2 at higher pressures, a (3 + 2) REMPI scheme is observed throughout this text.

Clark, Aleksander (ORCID:0009000855600042)↗
DOE OSTI2026

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

PowderJet: Spherical metal powder production via multi-orifice droplet-on-demand metal jetting

Leading metal additive manufacturing techniques, such as laser powder bed fusion and directed energy deposition, rely on high-quality spherical metal powders. However, traditional powder production methods like gas atomization face limitations, including low in-spec yield, asphericity, and internal porosity. We introduce PowderJet, a powder production platform that uses electromagnetic pulses to eject liquid metal droplets from a multi-orifice nozzle. Unlike stochastic methods, PowderJet tightly controls powder size, distribution, and purity through a droplet-on-demand approach. We detail the system’s design, operation, and performance using a combined experimental and computational fluid dynamics (CFD) framework. Initial results with Al4008 and Cu110 alloys demonstrate successful production, yielding unsieved aluminum powder batches with a mean diameter of 200 µm and a narrow size distribution (15 µm standard deviation). The produced powders are highly spherical, achieving a roundness > 0.95. PowderJet operates with a small melt volume (3 mL) and supports continuous refilling, enabling production rates between 30 and 140 cm³/hr depending on jetting frequency, number of orifices and particle size. CFD simulations show that future systems could achieve rates exceeding 1000 cm³/hr for particle sizes as small as 40 µm. PowderJet’s high yield of in-spec powder makes it ideal for producing precious or hazardous materials that are inefficient to manufacture using conventional methods. This platform offers a scalable, precise, and efficient solution for producing high-quality powders tailored for advanced manufacturing applications.

Atomization↗
DOE OSTI2024

Studying CPU and memory utilization of applications on Fujitsu A64FX and Nvidia Grace Superchip

ARM-based manycore CPU architectures are well-positioned to provide the rising memory throughput requirements of modern data intensive scientific applications in High Performance Computing (HPC). The Fujitsu A64FX CPU platform is based on the ARM v8.2A architecture, and is the processor of the flagship Japanese supercomputer - "Fugaku", which was previously ranked as the #1 supercomputer in the world according to the Top500 list. The Nvidia Grace superchip features 144 Neoverse V2 cores based on the ARMv9 architecture with 4x128b SVE2, providing exceptional computational power. The chip supports up to 480GB of memory, making it ideal for AI, machine learning, and scientific computing workloads. In this paper, we conduct a thorough performance exploration of a variety of parallel bandwidth-sensitive benchmarks and applications compiled with the native Fujitsu compiler on a Fugaku A64FX compute node and ARM (LLVM) Compiler on an NVIDIA Grace superchip compute node, engaging all the computational cores per cluster using OpenMP multithreading (assuming the cores can drive the available bandwidth). Our ultimate goals are to study the resource utilization of scientific applications and benchmarks on A64FX and Grace superchip, considering graph application scenarios ( GAP Benchmark suite) and eleven appli- cation proxies from the Rodinia heterogeneous benchmark suite (considering domains such as Data Mining, Bioinformatics, Fluid Dynamics, Pattern Recognition, etc.). Through exhaustive performance monitoring, we quantify the resource utilization of diverse OpenMP-based HPC applications on both the Fujitsu A64FX and the Nvidia Grace Superchip platforms.

benchmarking, Performance Analysis, High performan↗
DOE OSTI2025

An alternative electrode design for accurate spectral induced polarization measurements in variably saturated porous media

SUMMARY Applications of spectral induced polarization (SIP) require electrodes that maintain hydrologic contact with surrounding soils to capture small electrical responses, often observed as phase shifts in milliradians. For unsaturated soils, electrodes must overcome increased electrical contact impedance due to reduced pore fluid. Traditional designs use a ceramic membrane electrode (CME) with a water reservoir and metal conductor, requiring periodic maintenance to retain electrolytic solution. For field applications where maintenance is impractical, alternative designs are needed. This study evaluated a new electrode design (silica flour electrode, SFE) alongside a CME. SFEs use packed silica flour to store water via capillary forces against a metal conductor. This study examined both designs in three variably saturated soils at soil suctions up to 700 mbar and soil water contents below 1 percent, with SIP measurements across 0.01 to 10 000 Hz frequencies. SFEs match CMEs at high frequencies and perform better at lower frequencies, without requiring ongoing maintenance, making them ideal for field use. In water-only experiments, CMEs produced errors and high noise below 1.5 Hz, whereas SFEs were more accurate. However, CMEs performed better above 300 Hz. In fine sand, SFEs performed better due to the relatively lower contact impedance as compared to CMEs. Both electrode types performed comparably in silty sand and silt loam soils, although CMEs required ongoing maintenance, suggesting potential for long-term reliability issues.

58 GEOSCIENCES↗
DOE OSTI2025

Deconvoluting thermomechanical effects in X-ray diffraction data using machine learning

X-ray diffraction is ideal for probing the sub-surface state during complex or rapid thermomechanical loading of crystalline materials. However, challenges arise as the size of diffraction volumes increases due to spatial broadening and because of the inability to deconvolute the effects of different lattice deformation mechanisms. Here, we present a novel approach that uses combinations of physics-based modeling and machine learning to deconvolve thermal and mechanical elastic strains for diffraction data analysis. The method builds on a previous effort to extract thermal strain distribution information from diffraction data. The new approach is applied to extract the evolution of the thermomechanical state during laser melting of an Inconel 625 wall specimen which produces significant residual stress upon cooling. A combination of heat transfer and fluid flow, elasto-plasticity and X-ray diffraction simulations is used to generate training data for machine-learning (Gaussian process regression, GPR) models that map diffracted intensity distributions to underlying thermomechanical strain fields. First-principles density functional theory is used to determine accurate temperature-dependent thermal expansion and elastic stiffness used for elasto-plasticity modeling. The trained GPR models are found to be capable of deconvoluting the effects of thermal and mechanical strains, in addition to providing information about underlying strain distributions, even from complex diffraction patterns with irregularly shaped peaks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Multiport high-pressure synchrotron x-ray microscopy cell (CRADA Final Report)

The LBNL team developed a multiport cell allowing for fluid flow, with electrical contacts for sample bias, and with reference electrodes. The cell is currently installed and operating in beam-line 11.0.2 of the ALS, the Berkeley Synchrotron. The cell is closed by a membrane consisting of several stacked layers of graphene (1 to 3). The graphene membrane is supported on a 100 nm thick Si3N4 membrane covered with a 20 nm thin gold film for improved adhesion of the graphene and to ensure electrical conductivity. The Si3N4 film is perforated with a periodic array of holes 1 to 2 micrometers diameter. The cell can withhold an internal pressure (gas or liquid) of up to 10 atmospheres for in situ/operando studies using X-ray Photoelectron emission (XPS), and X-ray Absorption Spectroscopy (XAS). It is also ideal for tip-enhanced nano-Infrared Spectroscopy studies of molecules at the interface between the graphene and the liquid or gas filling the cell. Materials studied in the LBNL laboratory are in the form of nanoparticles, or thin films (few nm) deposited on the side of the graphene membrane facing the solution or gas inside the cell.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Kink Instability of Flux Ropes in Partially Ionized Plasmas

In the solar atmosphere, flux ropes are subject to current-driven instabilities that are crucial in driving plasma eruptions, ejections, and heating. A typical ideal magnetohydrodynamics instability developing in flux ropes is the helical kink, which twists the flux rope axis. The growth of this instability can trigger magnetic reconnection, which can explain the formation of chromospheric jets and spicules, but its development has never been investigated in a partially ionized plasma (PIP). Here, we study the kink instability in PIP to understand how it develops in the solar chromosphere, where it is affected by charge-neutral interactions. Partial ionization speeds up the onset of the nonlinear phase of the instability, as the plasma β of the isolated plasma is smaller than the total plasma β of the bulk. The distribution of the released magnetic energy changes in fully ionized plasma and PIP, with a larger increase in internal energy associated with the PIP cases. The temperature in PIP increases faster also due to heating terms from the two-fluid dynamics. PIP effects trigger kink instability on shorter time scales, which is reflected in more explosive chromospheric flux rope dynamics. These results are crucial to understanding the dynamics of small-scale chromospheric structures—minifilament eruptions—that thus far have been largely neglected but could significantly contribute to chromospheric heating and jet formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Experimental Study of Instabilities in Hydrogen-Air Fueled Rotating Detonation Combustion Presentation

Conventional gas turbine engines rely on an idealized constant pressure combustion process that in reality produces a pressure decrease as a result of viscous and other non-reversible losses. An alternative approach is rotating detonation combustion (RDC) which is a form of pressure gain combustion in which one or more detonation waves propagate an annular channel resulting in an increase in pressure across, subsequently providing greater work availability compared to deflagration ultimately leading to opportunities for greater thermodynamic efficiency when used in gas turbine engines that conventionally relies on constant. Modern gas turbine engines often rely on pre-mixed reactants to limit NOx emissions, although this may result in greater susceptibility to instabilities such as flashback and thermoacoustic oscillation, particularly for applications that utilize hydrogen as the fuel. Research in RDC has focused on non-premixed reactants thus limiting the occurrence of flashback, and high frequency detonation wave propagation (kHz) may interfere with the occurrence of thermoacoustic oscillations. Thermal NOx emissions are controlled through rapid combustion and sudden expansion of the working fluid. Although RDC may not be susceptible to instabilities encountered in conventional state of the art gas turbine engine combustion, there may be other mechanisms occurring that support instabilities that could be detrimental to performance.

Weber, Justin↗
DOE OSTI2026

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗
DOE OSTI2024

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗
DOE OSTI2025

Wave topology in Hall magnetohydrodynamics

Hall magnetohydrodynamics (HMHD) extends ideal MHD by incorporating the Hall effect via the induction equation, making it more accurate for describing plasma behavior at length scales below the ion skin depth. Despite its importance, a comprehensive description of the eigenmodes in HMHD has been lacking. In this work, we derive the complete spectrum and eigenvectors of HMHD waves and identify their underlying topological structure. We prove that the HMHD wave spectrum is homotopic to that of ideal MHD, consisting of three distinct branches: the slow magnetosonic-Hall waves, the shear Alfvén-Hall waves, and the fast magnetosonic-Hall waves, which continuously reduce to their ideal MHD counterparts in the limit of vanishing Hall parameter. Contrary to a recent claim [Mahajan, Sharma, and Lingam, Phys. Plasmas 31, 090701 (2024)], we find that HMHD does not admit any additional wave branches beyond those in ideal MHD. In conclusion, the key qualitative difference lies in the topological nature of the HMHD wave structure: it exhibits nontrivial topology characterized by a Weyl point—an isolated eigenmode degeneracy point—and associated nonzero Chern numbers of the eigenmode bundles over a 2-sphere in 𝐤-space surrounding the Weyl point.

Alfvén waves↗