Search NASA⌕ Search

SEARCH · Search NASA

Results for “encounter”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Frhodo: A program for simulating chemical kinetic measurements and optimizing kinetic mechanisms

Chemical kinetic simulations are frequently used to extract kinetic and mechanistic information from experimental data by fitting simulations to the data. For simple systems of consecutive single channel reactions this is a relatively straightforward task. However, as the complexity of the system increases and particularly for systems involving competitive multi-channel reactions manually optimizing a mechanism to obtain the best fits to a range of experimental data quickly becomes a time-intensive and challenging process, especially when rate coefficients have pressure and temperature dependencies. There is considerable potential for exploiting automated optimization methods to rapidly screen reaction mechanisms and optimize the fit by adjusting rate coefficients within well-defined constraints. A new chemical kinetics simulation program, Frhodo, has been developed for this purpose. The program allows either manual simulation of individual experiments or fully automated optimization against a range of experimental data by adjusting user selected reactions within user defined constraints. Frhodo incorporates machine learning-based optimization routines that work by either minimizing standardized residuals or through Bayesian parameter estimation. Frhodo's optimization capabilities are demonstrated by reexamining two previously studied systems of relevance to combustion. Dissociation of diacetyl followed by recombination of methyl radicals is an example of a sequence of single channel, consecutive reactions. Pyrolysis of 2-methyl furan exemplifies multi-channel, unimolecular reactions that are often encountered and have temperature and pressure dependencies in the rate coefficients of the channels and branching ratios between them. For both systems the optimization routines resulted in solutions similar to the original studies. Finally, comments are made on the strengths and limitations of the approach.

0D simulation↗
DOE OSTI2024

Flame kinetics at scramjet-engine-relevant conditions: Role of prompt dissociation of weakly-bound radicals

Combustion in high-speed ram-based propulsion engines occurs under distinct thermodynamic conditions of high reactant temperatures (greater than 1000 K) and relatively low pressures (<5 atm). There is a lack of fundamental flame measurements at such conditions that result in adiabatic flame temperatures (T ad ) exceeding 2500 K. In this work, we have measured laminar flame speeds of oxygen-enriched CH 4 /oxidizer mixtures at sub-atmospheric conditions to probe kinetics at high T ad using the isobaric spherically expanding flame approach. Simulations with recent kinetic models revealed increasing differences between data and model predictions with increasing T ad , reaching up to 25 %. Kinetic analyses reveal that at the thermodynamic conditions in these O 2 -enriched flames, i.e., lower pressures and higher T ad , the effects of HCO prompt dissociation are accentuated. In addition to HCO, the prompt dissociations of CH 2 OH and C 2 H 5 are also considered. Here, the prompt dissociations of all three radicals were evaluated and their effects considered in flame speed simulations. Reaction path analysis for the present flames revealed that approximately half of the reaction flux for HCO formation undergoes prompt dissociation to H + CO. Furthermore, these analyses also revealed that the pathways and sensitive reactions are similar between oxygen-enriched fuel/oxidizer mixtures and preheated fuel/air mixtures, if both have similar T ad . Thus, flames of oxygen-enriched mixtures could be a surrogate to probe the flame chemistry of highly preheated mixtures at relatively low pressures that are often encountered in ram-based propulsion engine combustors.

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

On the numerical sensitivity of cellular automata grain structure predictions to large thermal gradients and cooling rates

Cellular automata (CA) models of as-solidified grain structure, originally developed and applied to casting, have become a common means of predicting grain structure resulting from Additive Manufacturing (AM) processes. The majority of these models are based on the decentered octahedron approach, which attempts to correct for the effect of grid anisotropy on the prediction of competitive solidification of dendritic grains. However, AM solidification occurs under cooling rates ($\dot{T}$) and thermal gradients (G) that are orders of magnitude larger than those encountered in casting, and no systematic investigation on the effect of the CA model cell size (Δx) and time step (Δt) on AM microstructure predictions has been performed. Here, in this study, such an investigation is first performed via simulation of individual grains of various crystallographic orientations with a fixed, unidirectional G, showing that CA prediction of the steady-state undercooling matched the expected values based on the interfacial response function at small G and deviated from the expected values at large G. Simulation of competitive growth of multiple grains showed a weakening of the predicted texture as G and Δx became large. Simulation of solidification under AM conditions, where G and $\dot{T}$ vary spatially across the melt pools, showed that not only does grain selection weaken and deviate from expectations at large Δx, but grains with crystallographic $\langle$100$\rangle$ aligned with the grid directions are more adversely affected by the temperature field discontinuities than grains with other crystallographic orientations. Despite the fact that the exact grain competition results depended on Δt, the overall texture development was notably less sensitive to Δt than Δx, provided that a reasonable value of Δt is selected based on the ratio of Δx to the maximum local solidification velocity in the simulation domain. Finally, from the directional solidification and AM simulation results, an analysis of computational cost compared to simulation resolution is performed based on an equation derived to quantify the relatively inaccuracy in grain selection based on the model and temperature field inputs. From this analysis, it is concluded that there is a need for algorithmic improvements to improve CA grain competition accuracy for large G processing conditions as sufficiently small Δx to resolve the necessary competition is intractable for many AM processing conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

ExaCA v2.0: A versatile, scalable, and performance portable cellular automata application for additive manufacturing solidification

The previously established ExaCA software for performance portable alloy grain structure simulation has been updated to better represent the solidification behavior during complex alloy processing conditions, such as those encountered during metal additive manufacturing (AM), and for improved performance and scalability. Here, an extension to the time–temperature history input data format and the core ExaCA algorithm to include an arbitrary number of melting and solidification events yielded improved prediction of texture for various melt pool geometries, expanding the range of AM-relevant conditions that can be accurately simulated. Improved heat transport process simulation coupling, including the creation of large raster datasets from single track time–temperature history data and in-memory coupling with the new, performance portable finite difference code Finch, were also demonstrated in example studies on the effect of multilayer AM microstructure predictions on hatch spacing and cell size, respectively. Additional new features are detailed and demonstrated, including the ability to perform simulations using various interfacial response function forms, execute simulations on state-of-the-art hardware, improved usability through post-processing versatility, and improved strong and weak scaling performance. The performance, physics, and versatility improvements demonstrated here will further enable large-scale studies on AM process–microstructure relationships that were not previously possible. Furthermore, the usability improvements and ability to run coupled AM process–microstructure simulations using the Finch-ExaCA workflow will facilitate broader use of this open-source software by the computational materials community.

36 MATERIALS SCIENCE↗
DOE OSTI2024

PlasmoData.jl — A Julia framework for modeling and analyzing complex data as graphs

Datasets encountered in scientific and engineering applications appear in complex formats (e.g., images, multivariate time series, molecules, video, text strings, networks). Graph theory provides a unifying framework to model such datasets and enables the use of powerful tools that can help analyze, visualize, and extract value from data. In this work, we present PlasmoData.jl, an open-source, Julia framework that uses concepts of graph theory to facilitate the modeling and analysis of complex datasets. The core of our framework is a general data modeling abstraction, which we call a DataGraph. We show how the abstraction and software implementation can be used to represent diverse data objects as graphs and to enable the use of tools from topology, graph theory, and machine learning (e.g., graph neural networks) to conduct a variety of tasks. We illustrate the versatility of the framework by using real datasets: (i) an image classification problem using topological data analysis to extract features from the graph model to train machine learning models; (ii) a disease outbreak problem where we model multivariate time series as graphs to detect abnormal events; and (iii) a technology pathway analysis problem where we highlight how we can use graphs to navigate connectivity. Further, our discussion also highlights how PlasmoData.jl leverages native Julia capabilities to enable compact syntax, scalable computations, and interfaces with diverse packages. Overall, we show that the DataGraph abstraction and PlasmoData.jl Julia package are able to model data within graphs and enable useful analysis.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2020

Simulation of axial tensile well deformation during reservoir compaction in offshore unconsolidated methane hydrate-bearing formation

Sand production encountered in the 2013 offshore field gas production tests at the Nankai Trough, Japan, could be attributed to well failure during reservoir compaction. In this study, well integrity under various reservoir compaction patterns for the Nankai Trough case is examined using a well-formation finite element model. Here, the modelling details include the inclusion of a cement sheath as well as the modelling of construction processes (such as cement shrinkage). Well elongation in the overburden layer becomes significant when the reservoir subsidence is localized near the wellbore under large depressurization. Results show that the maximum plastic deviatoric strain level in the cement could reach 0.7% when the maximum reservoir subsidence reaches 0.85 m and cement shrinkage is limited. When cement shrinkage rises to 0.75%, the maximum plastic deviatoric strain increases to 2.4% as the cement accumulates additional plastic strain during shrinkage due to its deformation being constrained by the casing. In order to prevent the cement from failure, it might be effective to hold the pressure drawdown at a low level (e.g., several MPa) until the hydrate dissociation front advances to a certain radius from the well (e.g., a couple of tens of metres).

58 GEOSCIENCES↗
DOE OSTI2025

RLGBS: Reinforcement Learning-Guided Beam Search for process optimization in a paper machine dryer section

Paper drying is responsible for over two-thirds of energy consumption in the U.S. pulp and paper industry, presenting significant potential for energy savings through optimization of process parameters. Current approaches often assume fixed operating conditions, neglecting dynamic ambient and process variations that limit achievable savings and real-world applicability. To this end, we develop a physics-based simulation environment for a paper machine dryer section and propose a reinforcement learning (RL) framework to minimize overall energy consumption by optimizing drying process parameters under diverse operating conditions. To mitigate overdrying and numerical instabilities caused by suboptimal local RL actions, we introduce Reinforcement Learning-Guided Beam Search (RLGBS), which explores multiple action sequences in parallel using beam search. Instead of making step-by-step decisions, RLGBS prioritizes solutions based on cumulative probability, reducing the impact of individual suboptimal actions. Experiments demonstrate that RLGBS achieves consistent energy savings under unseen operating conditions not encountered during training, outperforming conventional RL methods. While validated in drying optimization, this framework is broadly applicable to other RL-based industrial process control problems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

LieART 2.0 – A Mathematica application for Lie Algebras and Representation Theory

In this paper, we present LieART 2.0 which contains substantial extensions to the Mathematica application LieART (Lie Algebras and Representation Theory) for computations frequently encountered in Lie algebras and representation theory, such as tensor product decomposition and subalgebra branching of irreducible representations. The basic procedure is unchanged—it computes root systems of Lie algebras, weight systems and several other properties of irreducible representations, but new features and procedures have been included to allow the extensions to be seamless. The new version of LieART continues to be user friendly. New extended tables of properties, tensor products and branching rules of irreducible representations are included in the supplementary material for use without Mathematica software. LieART 2.0 now includes the branching rules to special subalgebras for all classical and exceptional Lie algebras up to and including rank 15.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Implicit highly-coupled single-ion Hall-MHD formulation for hybrid particle-in-cell codes

The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A multi-center quadrature scheme for the molecular continuum

We report a common way to evaluate electronic integrals for polyatomic molecules is to use Becke’s partitioning scheme (Becke and Chem, 1988) in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron–molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Becke’s partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A red algal polysaccharide influences the multicellular development of the choanoflagellate Salpingoeca rosetta

The choanoflagellate Salpingoeca rosetta exemplifies the capacity of marine microeukaryotes to integrate environmental cues into their life histories. Some of the cues that S. rosetta detects are lipids from their bacterial prey that induce the development of multicellular colonies called rosettes. In the aquatic environments where S. rosetta would encounter bacterial cues, microbial communities gain refuge and food from algae. For example, microbes in coastal environments benefit from macroalgae that annually release ∼10 14 g of carbon, primarily in the form of polysaccharides that comprise ∼80% of dissolved organic carbon in surface waters. Despite the ecological impact of algae and historical descriptions of choanoflagellates attached to algae, the influence of algae on choanoflagellate life history transitions has remained unexplored. Here, we report that porphyran, a polysaccharide produced by the red macroalga Porphyra umbilicalis, induces multicellular development in S. rosetta. We first noticed this response when S. rosetta formed multicellular rosette colonies while growing in media prepared from P. umbilicalis, even though rosette-inducing bacteria were absent. By biochemically purifying extracts of P. umbilicalis, we identified porphyran as the rosette-inducing factor. This response provides a biochemical insight that indicates algal polysaccharides may serve as cues to mediate associations between choanoflagellates and algae. Moreover, this observation complements evidence from environmental and geochemical studies that show the impact that algae have exerted on the ecology and evolution of eukaryotes, including a rise in algal productivity during the origin of animals, the closest living relatives of choanoflagellates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The polymeric conformational effect on capacitive deionization performance of graphene oxide/polypyrrole composite electrode

Exploitation of novel faradic materials is an alternative implementation for solving the problem of poor specific electrosorption capacity that conventional carbon materials are encountered in capacitive deionization. Particularly, composite electrode is just a suitable choice because of its potentially high ion-storage ability. Herein, a cyclic voltammetric treatment method with different low limit of potential window was used to manipulate the polymeric conformation and doping level of graphene oxide/polypyrrole (GO/PPy) composite electrode. Based on it, the effect of polymeric structure on the electrosorption performance was systematically studied. When the low limit of potential window is shifted negatively enough, the irreversible polymeric conformational shrinks of GO/PPy are promoted, which not only hinders the insertion process of ions, but also decreases the doping level of polymer due to the intensive interchain-action produced by more entangled polymeric chain. Thus, the number of intercalated ions should decrease, which is expressed by electrochemical impedance spectroscopy (EIS) results and is proportional to the electrosorption capacity of GO/PPy composite electrode in membrane capacitive deionization (MCDI) process. Finally, our work suggests that the less packing density, higher doping level and more charge delocalization on PPy backbone in electrode are beneficial to enhance its capacitive deionization performance.

42 ENGINEERING↗
DOE OSTI2023

Pore-scale simulation of multiphase flow and reactive transport processes involved in geologic carbon sequestration

Multiphase flow and reactive transport are two essential physicochemical processes that govern the effectiveness of geological carbon sequestration (GCS). The interaction and feedback among different phases and components during intricate physicochemical processes hold great significance in understanding CO 2 sequestration. Pore-scale simulations can account for multiphase flow and reactive transport processes in porous media and obtain spatial distributions of parameters (density, velocity, concentration, etc.) in the pore space as well as their temporal evolutions. This proves especially valuable considering that experiments can be hindered by constraints in spatial and temporal resolution. The comprehensive insights garnered from pore-scale research can be leveraged for continuum modeling using the representative elementary volume (REV) concept. In this contribution, four sequential mechanisms of CO 2 -brine-rock interaction in three zones delineated by CO 2 saturation are elaborated to elucidate complicated physicochemical processes involved in GCS, which are followed by general descriptions of mathematical equations and pore-scale numerical methods. In addition, as interested and commonly encountered processes, leakage risks during GCS and CO 2 -enhanced oil recovery (CO 2 -EOR) processes are presented. The existing challenges and future directions are discussed for both the performance of the pore-scale models and the current gaps in the field of GCS. Importantly, we expect that this review will prove beneficial for researchers interested in pore-scale simulations, GCS, and related disciplines.

58 GEOSCIENCES↗
DOE OSTI2025

Uniform secondary current distribution at disk electrodes under Tafel kinetics enabled by concentric current-shielding rings

Non-uniform secondary current distribution at rotating disk electrodes (RDE) is a common problem when using resistive electrolyte media or large applied currents. In a recent publication, we have shown that auxiliary electrodes such as the ring of a rotating ring disk electrode (RRDE) can help suppress current non-uniformities at the disk enabling reliable electroanalytical measurements. However, this previous work considered linear kinetics where current distribution non-uniformities were moderate. In the present contribution, we consider current distribution non-uniformities encountered under Tafel kinetics. We show, for the case of 2,5-dihydroxy-1,4-benzoquinone (DHBQ) reduction, that optimally chosen ring conditions serve to provide effective shielding at the disk edge rendering the overall disk current distribution to be uniform. Numerical modeling and scaling analysis (using the Wagner number) are presented to aid a user in determining the optimal ring current density for achieving uniform disk current distribution under Tafel kinetics. This approach is especially useful when studying soluble-soluble redox transitions for which, unlike deposit distribution in electrodeposition, the current distribution non-uniformity is not visually apparent.

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

Self-supported MoO 2 /MoS 2 nano-sheets embedded in a carbon cloth as a binder-free substrate for high-energy lithium–sulfur batteries

As one of the most prospective candidates for next-generation rechargeable batteries, lithium-sulfur (Li-S) batteries currently still encounter great challenges associated with the low conductivity, severe shuttle effects, and sluggish redox kinetics. Herein, a self-supported sulfur host is fabricated by an in-situ growth of MoO 2 /MoS 2 nano-sheets on a carbon cloth (CC). First, carbon cloth serves as a good template for the morphology-controlled synthesis of nanostructured materials, which not only alleviates the agglomeration of MoO2/MoS2, but also enhances the flexibility and mechanical strengthen of the hybrid architecture as a free-standing host. Second, the polysulfide-trapping ability can be greatly enhanced by both physical and chemical adsorption from the MoO 2 /MoS 2 -anchored carbon cloth. Moreover, the partially sulfurized MoO 2 /MoS 2 nano-sheets integrate the benefits of conductive MoO 2 and sulfiphilic MoS 2 , thus facilitating fast charge transfer and redox kinetics of polysulfide conversion. As a result, these attributes enable the host to hold a high sulfur loading (up to 7.6 mg cm -2 ), which exceeds most of the reported carbon cloth-related cathode work in the literature. Furthermore, the Li-S cells can achieve a high peak capacity of 1350 mA h g -1 , excellent rate capability (C/20 - 2C rate), impressive areal capacity (up to 6.3 mA h cm -2 ), and a high capacity retention of 85% after 100 cycles.

25 ENERGY STORAGE↗
DOE OSTI2024

Gaseous contaminant transfer in membrane-based air-to-air energy exchangers

Membrane-based air-to-air energy exchangers (M–AAEEs) transfer heat and moisture between building exhaust air and fresh ventilation air streams through a membrane, thereby reducing the energy required for conditioning the fresh ventilation air. Energy exchangers are typically used in buildings with relatively clean building exhaust air, such as office buildings and schools. Recently interest in using energy exchangers in a wider range of buildings has grown, to reduce the energy consumption associated with the heating, ventilating, and air-conditioning (HVAC) systems in these buildings. However, if the building exhaust air is not clean, as would be the case for laboratories or factories, new risks are encountered when using energy exchangers. It is possible that gaseous contaminants in the exhaust air may also transfer along with the moisture through the membranes, contaminating the incoming fresh ventilation air. Current test standards provide a test procedure to determine the contamination of the fresh ventilation air by measuring the transfer of an inert tracer gas in an energy exchanger. However, the tracer gas test may not represent the transfer of common indoor air contaminants due to differences in their transport properties. Therefore, in this study, an experimental facility is developed to determine the transfer of seven different contaminants through two membranes (porous and dense) at different flow rates. Contaminant transfer is quantified using a parameter called the exhaust contaminant transfer ratio (ECTR), which gives the fraction of the contaminants transferred from the exhaust air to the fresh ventilation air. A theoretical model based on the effectiveness-number of transfer units (ε-NTU) correlation and moisture transfer resistance of the membrane is presented to determine transfer through porous membranes and validated with experimental results. The major contribution of this paper is that it presents a simple method to predict the transfer of different contaminants through a porous membrane based on the moisture transfer resistance of the membrane at different operating conditions. It is found that as contaminant diffusivity decreases, ECTR generally also decreases, and as the flow rate increases, ECTR decreases, which is consistent with the predictions from the correlation.

42 ENGINEERING↗
DOE OSTI2024

Techno-economic feasibility of geothermal energy production using inactive oil and gas wells for district heating and cooling systems in Tuttle, Oklahoma

Oil and gas wells have been drilled in the United States and often encounter hot co-produced water possibly suitable for geothermal direct-use applications. Here this study evaluated geothermal resources, heating and cooling demand, and techno-economic potential of four oil and gas wells repurposed for geothermal energy production to serve two public schools and 250 nearby houses in Tuttle, Oklahoma, USA. Subsurface geology in the study area primarily consists of Permian to Mississippian sandstone and limestone formations approximately from 2 km to 3.5 km depth. With a geothermal gradient of 25.3 °C/km, the temperature of geothermal energy production was estimated as 68.2 °C at 2.1 km. The potential of calcite scaling within boreholes and pipes was minimal. Using the characterized reservoir properties, techno-economic analyses were conducted for evaluating levelized costs of geothermal energy production for direct-use heating and cooling and district heating systems with five different production scenarios. Heating and cooling demand in two schools and 250 houses and additional 0.3 MW and 0.6 MW base heating loads for potential geothermal direct-use applications in the study area were also incorporated in the modeling. Results indicated that the levelized cost of heating of geothermal district heating systems utilizing existing wells were significantly lower than those of conventional geothermal energy systems that include well drillings. The geothermal energy production with doublet and quartet configurations was generally sufficient to supply total annual heating demand in the district, while a peaking boiler was used to supply peak loads. Doublet configurations showed higher production temperature with a lower thermal drawdown than the quartet configurations. The doublet system at 2.1 km production depth with 9.3 kg/s flow rate represented the optimal levelized cost of $71/MWh and 91.2% utilization factor. These results imply that the Tuttle geothermal district heating system is techno-economically feasible and comparable to geothermal district heating systems in the United States and the regional natural gas rate for residentials by recycling inactive oil and gas wells with no drillings.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Advanced defrosting techniques in air source heat pumps: A review of vapor injection, thermal energy storage, and experimental frost accumulation data

Electrification is a critical step for reducing greenhouse gas emissions from heating. Air source heat pumps (ASHPs) are a promising alternative to fossil fuel-based systems due to their high coefficients of performance (COP), dual heating and cooling capability, and lower carbon footprint. However, for ASHPs to achieve widespread adoption, they must operate reliably across all climates, including cold regions. Additionally, defrosting techniques should be energy efficient and minimally disruptive to indoor comfort. Vapor injection (VI) technology can address the high-pressure and high-temperature lift challenges encountered in low ambient conditions. More recently, in addition to enhancing heating performance, VI has also been shown to improve the speed and efficiency of reverse cycle defrosting. Likewise, thermal energy storage (TES) has steadily gained attention for its ability to serve as an auxiliary heat source during both normal operation and defrosting. This review analyzes the benefits and limitations of VI- and TES-assisted defrosting approaches. While both technologies show strong potential individually, no studies to date have explored their combined use in ASHP systems. Additionally, to support continued development of defrosting strategies, both in modeling and experimental work, it is critical to establish frost accumulation data under a range of operating conditions. By compiling the available data from the literature, this paper also highlights the limited availability of such experimental data and the wide variation in frosting and defrosting durations and termination criteria, which are often influenced by system design and test setups.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗