Search NASASearch

SEARCH · Search NASA

Results for “Charge density waves”

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.

686 records · Page 7

Hydrofluoroether (HFE) Cleaning Fluid Replacement and Qualification

This NASA Engineering and Safety Center (NESC) Technical Assessment Report addresses the urgent need to identify, test, and qualify replacement cleaning solvents and technologies following 3M’s decision to discontinue production of hydrofluoroether (HFE) products, including Novec HFE-7100, by the end of 2025. HFE-7100 is widely used across NASA for precision cleaning of ground and flight hardware, especially in oxygen system components, and as a heat transfer fluid in high-power-density thermal control systems. The assessment evaluates commercially available replacement solvents, alternative cleaning technologies, and their compatibility with NASA’s operational requirements. Key findings include the identification of viable solvent replacements, their performance in contaminant removal, materials compatibility, and oxygen system safety. The report also explores novel cleaning approaches, including aqueous, supercritical CO₂, cryogenic aerosol, plasma, UV-ozone, magnetically optimized fluids, and induced charge active filtration, highlighting their merits, limitations, and readiness for NASA adoption. The NESC provides recommendations for solvent selection, facility adaptation, and continued technology maturation to ensure mission continuity and regulatory compliance.

Alternative Cleaning Technologies

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems

Design and Modeling of the Charge Readout of a SiMOS Quantum Dot with a Single Electron Transistor and CryoCMOS

Single electron spin qubits trapped in SiMOS quantum dots are a promising technology for scaling to thousand- or million-qubit systems due to their compatibility with mature CMOS manufacturing processes. A readout system that combines a single electron transistor with a custom cryogenic CMOS amplification and digitization chain offers key advantages by avoiding the use of bulky RF components or room-temperature interconnects. We present design techniques and simulation results for an optimized qubit-SET cryoCMOS interface, culminating in the design of the QNDR1 ASIC, the first cryogenic readout ASIC designed under the Quandarum project, which targets the development of a many-channel spin qubit based detector for use in high energy physics.

Quinn, Adam [Fermilab] (ORCID:0009000605056371)

Order–Disorder Phase Stabilization by Pressure‐Induced Charge Transfer Enhances the Ferroelectric Photovoltaic Effect in Multiferroic BaFe 4 O 7

Abstract Multiferroic ferroelectric photovoltaic (FPV) materials, combining magnetic and ferroelectric properties, are of paramount importance for optoelectronic and photovoltaic applications. However, optimizing both the remanent polarization and the optical bandgap—key factors for enhanced FPV performance—presents a significant challenge due to their trade‐off. This work shows that pressure‐induced charge transfer between different metal sites can break this trade‐off. Above ≈20 GPa, charge transfer between different trivalent iron (Fe) sites in the multiferroic material BaFe 4 O 7 leads to Fe valence disproportionation, FeO 4 tetrahedra disorder, and Jahn–Teller distortion of FeO 6 octahedra. These changes reduce the bandgap, lower resistivity, and enhance ferroelectric polarization, resulting in a 2.5‐fold increase in photocurrent. Upon decompression, BaFe 4 O 7 retains an order–disorder structure, optimal ferroelectric and optical properties at ambient conditions. This work provides a novel pathway to simultaneously optimizing ferroelectricity and bandgap via pressure‐induced charge transfer, overcoming the traditional trade‐off in FPV materials, and offers a promising approach for developing high polarization performance, narrow‐bandgap FPV materials.

Chemistry

Polyimide/Carbon Nanotube Composite Films for Electrostatic Charge Mitigation

Low color, space environmentally durable polymeric films with sufficient electrical conductivity to mitigate electrostatic charge (ESC) build-up have potential applications on large, deployable, ultra-light weight Gossamer spacecraft as thin film membranes on antennas, solar sails, thermal/optical coatings, multi-layer insulation blankets, etc.. The challenge has been to develop a method to impart robust electrical conductivity into these materials without increasing solar absorptivity (alpha ) or decreasing optical transparency or film flexibility. Since these spacecraft will require significant compaction prior to launch, the film portion of the spacecraft will require folding. The state-of-the-art clear, conductive coating (e.g. indium-tin-oxide, ITO) is brittle and cannot tolerate folding. In this report, doping a polymer with single-walled carbon nanotubes (SWNTs) using two different methods afforded materials with good flexibility and surface conductivities in the range sufficient for ESC mitigation. A coating method afforded materials with minimal effects on the mechanical, optical, and thermo-optical properties as compared to dispersal of SWNTs in the matrix. The chemistry and physical properties of these nanocomposites are discussed.

Smith, Joseph G., Jr.

Numerical Investigation of Fluid Flow and Space Charge in Liquid Argon Time Projection Chamber (LArTPC) Detectors

Overview This project focused on developing a high-fidelity numerical framework to simulate the multiphysics environment within Liquid Argon Time Projection Chamber (LArTPC) detectors. The primary objective was to characterize the complex interplay between ion transport, background fluid dynamics, and electric field distortions—a critical factor for the calibration and sensitivity of next-generation High Energy Physics experiments, such as DUNE. Technical Achievements The research successfully yielded a hybrid numerical space-charge solver utilizing a Cell-Centered Finite Volume Method (FVM) for ion transport coupled with a Finite Element Method (FEM) for electric potential. Key accomplishments include: • Verification & Validation: The 3-D solver was rigorously verified against 1-D analytical solutions, demonstrating high numerical accuracy in predicting space-charge-induced field deviations. • Field Distortion Analysis: 3D simulations revealed that space charge effects introduce significant non-uniformities in the electric field. Critically, the research identified that background LAr flow velocities, when comparable to ion drift velocities, markedly exacerbate these distortions. • Technology Transfer: The resulting source code and comprehensive user manuals were successfully transferred to collaborators at Fermilab, providing a portable computational tool for the broader scientific community. Challenges and Future Directions While the space-charge solver achieved all performance metrics, the integrated fluid dynamics modeling encountered convergence challenges stemming from the extreme 200-fold disparity in length scales between the detector's 37 mm inlet pipes and the 8-meter global domain. To address this, the project has identified a clear technical pivot toward Hierarchical Geometric Adaptive Mesh Refinement (HG-AMR). By implementing an h-type refinement strategy with hanging nodes, future iterations of this solver will be capable of resolving localized high-gradient inlet flows without the prohibitive computational costs of regular grids. This advancement, combined with data-driven uncertainty quantification based on MicroBooNE-style calibration, will enable the precise modeling of detector responses in large-scale cryogenic environments where direct measurement remains difficult. Impact The computational tools developed under this award provide a foundation for enhancing the energy resolution and spatial reconstruction of noble liquid detectors. By bridging the gap between theoretical fluid dynamics and experimental field calibration, this work supports the DOE’s mission to advance the frontiers of neutrino physics and dark matter detection.

42 ENGINEERING

Enhancing Charge Storage of Mo 2 Ti 2 C 3 MXene by Partial Oxidation

Abstract Driving the pseudocapacitive redox intercalation in 2DMXenes with neutral electrolytes is important for safer, more sustainable, and improved electrochemical charge storage. Single transition metal MXenes, such as Ti 3 C 2 , have shown great promise for energy storage, owing to their high conductivity and redox activity. Mixed metallic MXenes, such as out‐of‐plane ordered Mo 2 Ti 2 C 3 , have remained underexplored in energy storage because of the absence of redox activity in most of the electrolytes. Simultaneous structural modifications and instigating intercalation pseudocapacitance in neutral electrolytes could be a viable strategy for enhancing their electrochemical properties. Herein, a facile synthesis of partially oxidized Mo 2 Ti 2 C 3 MXene (PO‐Mo 2 Ti 2 C 3 ) exhibiting improved charge storage capability is demonstrated. Optical, structural, and spectroscopic analyses indicate the formation of oxide nanostructures upon thermal oxidation of Mo 2 Ti 2 C 3 . This leads to an enhanced energy storage capability with remarkably improved cyclability as well as a high Coulombic efficiency in a neutral LiCl electrolyte. This work highlights the importance of structural modifications of MXenes to enhance their charge storage and shows the promise of less explored double transition metal MXenes in energy storage.

Chemistry

Ultra-light antennas via charge programmed deposition additive manufacturing

Abstract The demand for lightweight antennas in 5 G/6 G communication, wearables, and aerospace applications is rapidly growing. However, standard manufacturing techniques are limited in structural complexity and easy integration of multiple material classes. Here we introduce charge programmed multi-material additive manufacturing platform, offering unparalleled flexibility in antenna design and the capability for rapid printing of intricate antenna structures that are unprecedented or necessitate a series of fabrication routes. Demonstrating its potential, we present a transmitarray antenna composed of an interconnected, multi-layered array of dielectric/conductive S-ring unit cells, reducing 94% mass of conventional antenna configurations. A fully printed circular polarized transmitarray system fed by a source and a Risley prism antenna system operating at 19 GHz both show close alignment between testing results and numerical simulations. This printing method establishes a universal platform, propelling discovery of new antenna designs and enabling data-driven design and optimizations where rapid production of antenna designs is crucial.

Science & Technology - Other Topics

Performance of Li-Ion Cells Under Battery Voltage Charge Control

A study consisting of electrochemical characterization and Low-Earth-Orbit (LEO) cycling of Li-Ion cells from three vendors was initiated in 1999 to determine the cycling performance and to infuse the new technology in the future NASA missions. The 8-cell batteries included in this evaluation are prismatic cells manufactured by Mine Safety Appliances Company (MSA), cylindrical cells manufactured by SAFT and prismatic cells manufactured by Yardney Technical Products, Inc. (YTP). The three batteries were cycle tested in the LEO regime at 40% depth of discharge, and under a charge control technique that consists of battery voltage clamp with a current taper. The initial testing was conducted at 20 C; however, the batteries were cycled also intermittently at low temperatures. YTP 20 Ah cells consisted of mixed-oxide (Co and Ni) positive, graphitic carbon negative, LIPF6 salt mixed with organic carbonate solvents. The battery voltage clamp was 32 V. The low temperature cycling tests started after 4575 cycles at 20 C. The cells were not capable of cycling. at low temperature since the charge acceptance at battery level was poor. There was a cell in the battery that showed too high an end-of-charge (EOC) voltage thereby limiting the ability to charge the rest of the cells in the battery. The battery has completed 6714 cycles. SAFT 12 Ah cells consisted of mixed-oxide (Co and NO positive, graphitic carbon negative, LiPF6 salt mixed with organic carbonate solvents. The battery voltage clamp was for 30.8 V. The low temperature cycling tests started after 4594 cycles at 20 C. A cell that showed low end of discharge (EOD) and EOC voltages and three other cells that showed higher EOC voltages limited the charge acceptance at the selected voltage limit during charge. The cells were capable of cycling at 10 C and 0 C but the charge voltage limit had to be increased to 34.3 V (4.3 V per cell). The low temperature cycling may have induced poor chargeability since the voltage had to be increased to achieve the required charge input. The battery has completed 6226 cycles. MSA 10 Ah cells consisted of Co oxide positive, graphitic carbon negative, LiPF6 salt mixed with organic carbonate solvents. The battery voltage clamp was 30.8 V. The low temperature cycling tests were started after 2182 cycles at 20 C. The cells were capable of cycling at 10 C and 0 C. Like SAFT, the voltage limit on charge had to be increased to 36 V (4.5 V per cell). There was a cell (cell S/N 13) in the battery that showed poor performance features such as low EOD voltage and high EOC voltage. The battery has completed 3441 cycles. A reconditioning procedure that consisted of C15 charge to a taper current of C/100 and C/20 discharge improved the voltage behavior of SAFT and MSA cells with no significant effect on YTP cells. We have demonstrated that the charge operation with VT clamp at battery rather than at cell level is feasible for onboard Li-Ion battery operation.

Rao, Gopalakrishna M.

Interfacial Adhesion Mechanism of Anionic Polyelectrolyte Brushes Induced by Oppositely Charged Macromolecular Counterions

Abstract Polyelectrolyte brushes are widely used as model systems for investigating electrostatic interactions at soft interfaces and to achieve exceptional lubrication properties. Previous studies have primarily focused on their behavior in the presence of multivalent counterions, which induce brush collapse, ionic crosslinking, and pronounced changes in interfacial structures. However, interactions between polyelectrolyte brushes and oppositely charged macromolecular counterions, such as polycations, remain poorly understood. Here, we prepare well-defined polystyrene sulfonate (PSS) brushes fabricated via surface-initiated grafting and employ surface forces apparatus measurements to investigate their interactions with oppositely charged polycations. In contrast to multivalent ions, polycations do not induce noticeable brush collapse, but instead generate significant adhesion between symmetric PSS brush layers. This adhesion increases with both contact time and applied load, eventually reaching a steady plateau, indicating that the interaction is governed not simply by electrostatic screening but by the gradual formation of polycation-mediated bridging under confinement. Furthermore, the introduction of monovalent Na+ ions disrupts the adhesive interaction even at very low concentrations, suggesting that the bridging function of the adsorbed polycation is highly sensitive to competitive ionic screening. In comparison, measurements with trivalent counterions reveal the expected brush collapse behavior but minimal dependence of adhesion on contact time, highlighting a clear mechanistic distinction from the polymeric counterion case. Collectively, these results demonstrate that molecular size, configurational restriction, and confinement-induced rearrangement, rather than charge valency alone, govern the interaction behavior of macromolecular counterions at brush interfaces. This work provides new insight into the molecular origins of adhesion and the regulation of interfacial interactions in charged polymer brush systems.

Park, Jinwoo [Argonne National Laboratory , , , ,;

Theory of the thermionic current beyond the traditional space charge limit enabled by trapped ions in the virtual cathode

We show that trapped ions in virtual cathode potential wells can raise the transmitted current of emitted electrons into a plasma much closer to the full emission than is predicted by cathode sheath theories without trapped ions. The transmitted current is controlled by the well barrier voltage, which must adjust to balance the creation of low-energy ions within the well, and their loss. Our model considers the case of a plasma-facing cathode where trapped ions are created passively via charge-exchange collisions and lost passively via thermal leakage over the well. We quantify these rates and estimate the current in terms of system parameters for thermionic emission into a plasma with several cathode geometries. A general prediction is that the current as a function of emitted flux does not saturate at the traditional space charge limit (the onset of a well) but can reach far higher values until the trapped ion balance breaks down, causing instability. The maximum stable current depends on parameters but in principle can be arbitrarily high if active techniques are used to manipulate the trapped ion balance. We conclude that major improvements in plasma technologies with hot cathodes might be achieved by optimizing the current enhancement enabled by trapped ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Flexible and Extensive Platinum Ion Gel Condensers for Programmable Catalysis

Catalytic condensers composed of ion gels separating a metal electrode from a platinum-on-carbon active layer were fabricated and characterized to achieve more powerful, high surface area dynamic heterogeneous catalyst surfaces. Ion gels comprised of poly(vinylidene difluoride)/1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl) imide were spin coated as a 3.8 μm film on a Au surface, after which carbon sputtering of a 1.8 nm carbon film and electron-beam evaporation of 2 nm Pt clusters created an active surface exposed to reactant gases. Electronic characterization indicated that most charge condensed within the Pt nanoclusters upon application of a potential bias, with the condenser device achieving a capacitance of ∼20 μF/cm 2 at applied frequencies of up to 120 Hz. The maximum charge of ∼10 14 |e – | cm –2 was condensed under stable device conditions at 200 °C on catalytic films with ∼10 15 sites cm –2 . Grazing incidence infrared spectroscopy measured carbon monoxide adsorption isobars, indicating a change in the CO* binding energy of ∼19 kJ mol –1 over an applied potential bias of only 1.25 V. Condensers were also fabricated on flexible, large area Kapton substrates allowing stacked or tubular form factors that facilitate high volumetric active site densities, ultimately enabling a fast and powerful catalytic condenser that can be fabricated for programmable catalysis applications.

Catalysis

Fingerprinting Superconductors by Disentangling Andreev and Quasiparticle Currents Across Tunable Tunnel Junctions

Tunneling Andreev reflection (TAR) spectroscopy provides a new approach to identify superconducting pairing symmetry at the atomic scale. Using atomistic superconducting transport simulations, we reveal the mechanism by which TAR can distinguish between pairing symmetries, which is complementary to traditional conductance-based techniques. In particular, owing to the additivity of the excess tunneling decay rate, the TAR spectrum is a weighted average of the contributions from quasiparticle currents, Andreev reflection, and higher-order scattering processes, and their relative weights depend on both the superconducting order parameter and the coupling strength. Within the local tunneling model, TAR dominates mid-gap conductance for s-wave superconductors, is suppressed for d-wave, and coexists with quasiparticle tunneling in sign-changing symmetries if the expectation value for the superconducting gap remains finite. Meanwhile, higher-order processes generally enhance the TAR signal when GN exceeds approximately 0.1G0. As a result, TAR provides a rich spectral fingerprint of the underlying pairing symmetry and electronic structure, enabling atomically resolved identification of unconventional superconducting states.

Maksymovych, Petro [Clemson University]

Modular Battery Charge Controller

A new approach to masterless, distributed, digital-charge control for batteries requiring charge control has been developed and implemented. This approach is required in battery chemistries that need cell-level charge control for safety and is characterized by the use of one controller per cell, resulting in redundant sensors for critical components, such as voltage, temperature, and current. The charge controllers in a given battery interact in a masterless fashion for the purpose of cell balancing, charge control, and state-of-charge estimation. This makes the battery system invariably fault-tolerant. The solution to the single-fault failure, due to the use of a single charge controller (CC), was solved by implementing one CC per cell and linking them via an isolated communication bus [e.g., controller area network (CAN)] in a masterless fashion so that the failure of one or more CCs will not impact the remaining functional CCs. Each micro-controller-based CC digitizes the cell voltage (V(sub cell)), two cell temperatures, and the voltage across the switch (V); the latter variable is used in conjunction with V(sub cell) to estimate the bypass current for a given bypass resistor. Furthermore, CC1 digitizes the battery current (I1) and battery voltage (V(sub batt) and CC5 digitizes a second battery current (I2). As a result, redundant readings are taken for temperature, battery current, and battery voltage through the summation of the individual cell voltages given that each CC knows the voltage of the other cells. For the purpose of cell balancing, each CC periodically and independently transmits its cell voltage and stores the received cell voltage of the other cells in an array. The position in the array depends on the identifier (ID) of the transmitting CC. After eight cell voltage receptions, the array is checked to see if one or more cells did not transmit. If one or more transmissions are missing, the missing cell(s) is (are) eliminated from cell-balancing calculations. The cell-balancing algorithm is based on the error between the cell s voltage and the other cells and is categorized into four zones of operation. The algorithm is executed every second and, if cell balancing is activated, the error variable is set to a negative low value. The largest error between the cell and the other cells is found and the zone of operation determined. If the error is zero or negative, then the cell is at the lowest voltage and no balancing action is needed. If the error is less than a predetermined negative value, a Cell Bad Flag is set. If the error is positive, then cell balancing is needed, but a hysteretic zone is added to prevent the bypass circuit from triggering repeatedly near zero error. This approach keeps the cells within a predetermined voltage range.

Button, Robert

Staged Z-pinch experiment at the double-EAGLE pulsed power facility

The staged Z-pinch is a potential high-energy gain fusion concept where a high atomic number liner implodes on a deuterium target using a pulsed multi-MA current source. Over the past several years this concept has been studied on 0.5–1.0 MA facilities using Ar and Kr gas puffs injected at R = 1.2 cm. Neutron yield up to (2.5 ± 0.34) x 10 10 was measured, Ruskov et al (2023 IEEE Trans. Plasma Sci. 51 3310–6). Here, in this study, we present experimental results from the 4 MA, 110 ns current rise time Double-EAGLE facility at L3Harris (currently, Fisica Inc.) where a larger radius nozzle created gas density profiles peaked at R = 2.5 cm. Modeling with the MACH2 and FLASH codes indicates that the larger radius allows stronger acceleration of the liner plasma and generation of shock waves that preheat a target plasma layer next to the liner to temperatures $T_i >$1 keV. The resulting counter thermal pressure on the liner plasma limits the growth of the Magneto-Rayleigh–Taylor (MRT) instability at the liner-vacuum boundary and the implosion proceeds in a relatively stable manner. Near bang time the enormous target plasma thermal pressure smoothens the MRT perturbations developed during the earlier implosion stages. Time integrated x-ray pinhole images with cutoff energy of 100 eV confirm that a long (∼3 cm), stable and uniform high energy density plasma column is formed in the final implosion stage. Consistent neutron yield in the 10 10 –10 11 range was measured for both Ar and Kr liners imploding on a deuterium target.

Physics - Plasma physics

Broadband and Tunable Microwave Absorption Properties from Large Magnetic Loss in Ni–Zn Ferrite

Highly effective electromagnetic (EM) wave absorber materials with strong reflection loss (RL) and a wide absorption bandwidth (EBW) in gigahertz (GHz) frequencies are crucial for advanced wireless applications and portable electronics. Traditional microwave absorbers lack magnetic loss and struggle with impedance matching, while ferrites are stable, exhibit excellent magnetic and dielectric losses, and offer better impedance matching. However, achieving the desired EBW in ferrites remains a challenge, necessitating further composition design. In this study, impedance matching is successfully enhanced and EBW in Ni–Zn ferrite is broadened by successive doping with Mn and Co , without incorporation of any polymer filler. It is found that Ni 0.4 Co 0.1 Zn 0.5 Fe 1.9 Mn 0.1 O 4 material exhibits exceptional EM wave absorption, with a maximum RL of −48.7 dB. It also featured a significant EBW of 10.8 GHz, maintaining a 90% absorption rate (RL < −10 dB) for a thickness of 4.5 mm. These outstanding properties result from substantial magnetic losses and favorable impedance matching. These findings represent a significant step forward in the development of microwave absorber materials, addressing EM wave pollution concerns within GHz frequencies, including the frequency band used in popular 5G technology.

36 MATERIALS SCIENCE

Uncovering the linear boron environment in Na 3 BP 2 through solid-state 11 B NMR spectroscopy

Boron-based compounds exhibit a wide range of structural diversity, with potential applications spanning organic and inorganic chemistry. Herein, we focus on the characterization of the linear boron-phosphorus unit P═B═P in Na 3 BP 2 using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and density functional theory (DFT) calculations. High-resolution 11 B ssNMR spectra were recorded at two fields, and key parameters such as chemical shift anisotropy (CSA), quadrupolar coupling constants (C Q ), and electric field gradient (EFG) tensors were extracted. The 11 B NMR results revealed a distinct chemical environment for the two-coordinate boron atom, with a CSA span (Ω) of 280 ppm and a C Q of 3.0 MHz. These values were further validated through periodic plane-wave DFT calculations, which showed good agreement with experimental results. The obtained spectral parameters are compared to other linear boron units, such as the BO 2 motif, providing a broader context for understanding boron coordination in inorganic compounds. This work expands the body of NMR knowledge on boron-containing materials, particularly for linear boron motifs. The findings contribute to the growing field of boron chemistry and its potential applications in advanced materials.

Porter, Andrew P. [Ames Laboratory (AMES), Ames, I

Construction of an Exact Pressure-Equilibrium Scheme for the Five-Equation Two-Phase Flow Model With Thermal Relaxation

Numerical simulation of compressible multiphase flows based on the four-equation (homogeneous relaxation) model is known to suffer from two fundamental difficulties with (a) wave propagation and (b) pressure equilibrium preservation. First, the mixture sound speed exhibits non-monotonic dependency with respect to the volume fraction, which leads to robustness issues in the resolution of shocks and acoustic wave propagation across two-phase regions. This difficulty can be mitigated by solving Allaire’s five-equation model augmented with infinitely fast phasic temperature equilibrium, from which solutions of the four-equation model can be recovered. However, when temperature is non-uniform, this augmented five-equation formulation still fails to preserve pressure equilibrium across material interfaces. In this work, we propose a fully conservative numerical scheme that exactly preserves pressure equilibrium at the discrete level for the augmented five-equation model, for arbitrary initial distributions of temperature and volume fraction. Combined with the monotonic sound speed property of the five-equation formulation, the proposed pressure-equilibrium preserving scheme significantly improves robustness in the presence of strong multiphase interactions, including shock–interface interactions and advection of material interfaces.

ESG