SEARCH · Engineering Papers
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.
PHASTRAN analyzes heat-transfer and flow behaviors of materials undergoing phase changes. Many phase changes operate over range of accelerations or effective gravitational fields. To analyze such thermal systems, it is necessary to obtain simultaneous solutions for equations of conservation of energy, momentum, and mass, and for equation of state. Written in APL2.
Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers is typically computationally intensive due to the relative complex time-dependent underlying physics. This paper presents a study investigating a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 37oC nominal melting temperature, leveraging automated parametric CFD simulations. The computational domain is reduced to 2D axisymmetric for a segment of the HX comprising a region in between two fins with constant wall temperature as boundary condition. This approach allows for faster evaluations of the transient thermal characteristics for a given design space of interest. Correlations and reduced order models (ROM’s) are developed. Results indicate that the ROM’s reproduce the CFD simulations with high accuracy for designs from the training data set, however, designs from cross-validation and test sets do not always have good agreement. The limited number of designs in the training set is a possible cause for the lack of robustness of the ROM’s and thus additional training designs will be added in future work. The ROM’s solve the simulations on average 500,000 times faster than CFD.
Natural convection effects in phase change thermal control devices were studied. A mathematical model was developed to evaluate natural convection effects in a phase change test cell undergoing solidification. Although natural convection effects are minimized in flight spacecraft, all phase change devices are ground tested. The mathematical approach to the problem was to first develop a transient two-dimensional conduction heat transfer model for the solidification of a normal paraffin of finite geometry. Next, a transient two-dimensional model was developed for the solidification of the same paraffin by a combined conduction-natural-convection heat transfer model. Throughout the study, n-hexadecane (n-C16H34) was used as the phase-change material in both the theoretical and the experimental work. The models were based on the transient two-dimensional finite difference solutions of the energy, continuity, and momentum equations.
Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.
Phase-change materials play a critical role in industrial energy storage applications to drive efficiency improvements, thermal energy management, and carbon emissions reductions. Recently, it has been shown that rapid solidification of alloys with metastable immiscibility in the liquid phase has the potential to form unique microstructures in which a low-melting phase is uniformly distributed in a high-melting matrix. This feature can be exploited using additive manufacturing to produce components with complex geometries containing such unique phase-change microstructures. Phase-field simulations utilizing high-performance computing were used to provide a detailed description of the evolution of the active phase during service in terms of their morphology and composition in different polycrystalline matrix grain morphologies that are typically produced during additive manufacturing. Phase field simulations were performed using, MEUMAPPS-SL (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid Liquid) code that was developed in-house by the Oak Ridge National Laboratory. The simulations utilized the capabilities of the Kestrel supercomputer at the National Renewable Energy Laboratory. The simulation results were compared with experimental results generated at Siemens Energy, Inc. The results indicate that the kinetics of liquid spreading along grain boundaries is largely determined by the mobility of the triple line along the intersection of the grain boundary liquid and the grain boundary plane.
Energy storage technologies, particularly those utilizing phase change materials (PCMs), have gained attention for their high energy density and efficient thermal management. PCMs, which store energy through solid-liquid phase transitions, can efficiently capture and release thermal energy, but face the challenge of leakage during the phase change process. Inorganic PCMs, such as salt hydrates, offer high energy storage capacity, but are difficult to encapsulate due to their corrosive nature. Conventional encapsulation techniques for inorganic PCMs are limited, particularly for scalable applications. In this work, we present an innovative method for the encapsulation of salt hydrate-based inorganic PCMs (CaCl 2 ·6H 2 O) using co-axial electrospinning. The process involves the creation of co-axial fibers, with salt hydrate as the core and polymer (e.g., PVP) as the outer shell, effectively preventing leakage and improving the stability of the PCM. This approach demonstrates the potential for scalable microencapsulation of inorganic PCMs, marking the first report of using co-axial electrospinning for this purpose. This novel technique could contribute to enhancing the performance and applicability of PCMs in thermal energy storage systems and other energy efficiency applications.
The synthesis and characterization of phase change materials infused into metal–organic framework crystals to demonstrate high latent heats with suppressed phase change temperatures.
A highly accurate and efficient finite-difference method for phase-change problems with multiple moving boundaries of irregular shape is developed by employing a coordinate transformation that immobilizes moving boundaries and preserves the conservative forms of the original governing equations. The numerical method is first presented for one-dimensional phase-change problems (involving large density variation between phases, heat generation, and multiple moving boundaries) and then extended to solve two-dimensional problems (without change of densities between phases). Numerical solutions are obtained non-iteratively using an explicit treatment of the interfacial mass and energy balances and an implicit treatment of the temperature field equations. The accuracy and flexibility of the present numerical method are verified by solving some phase-change problems and comparing the results with existing analytical, semi-analytical and numerical solutions. Results indicate that one- and two-dimensional phase-change problems can be handled easily with excellent accuracies.
There is empirical evidence that freely-suspended triboelectrostatically charged particulate clouds of dielectric materials undergo rapid conversion from (nominally) monodispersed "aerosols" to a system of well-defined grain aggregates after grain motion or fluid turbulence ceases within the cloud. In United States Microgravity Laboratory Space Shuttle experiments USML-1 and USML-2, it was found that ballistically-energized grain dispersions would rapidly convert into populations of filamentary aggregates after natural fluid (air) damping of grain motion. Unless continuously disrupted mechanically, it was impossible to maintain a non-aggregated state for the grain clouds of sand-size materials. Similarly, ground- based experiments with very fine dust-size material produced the same results: rapid, impulsive "collapse" of the dispersed grains into well-defined filamentary structures. In both ground-based and microgravity experiments, the chains or filaments were created by long-range dipole electrostatic forces and dipole-induced dielectric interactions, not by monopole interactions. Maintenance of the structures was assisted by short-range static boundary adhesion forces and van der Waals interactions. When the aggregate containers in the USML experiments were disturbed after aggregate formation, the quiescently disposed filaments would rearrange themselves into fractal bundles and tighter clusters as a result of enforced encounters with one another. The long-range dipole interactions that bring the grains together into aggregates are a product of randomly-distributed monopole charges on the grain surfaces. In computer simulations, it has been shown that when the force vectors of all the random charges (of both sign) on a grain are resolved mathematically by assuming Coulombic interaction between them, the net result is a dipole moment on individual grains, even though the grains are electrically neutral insofar as there is no predominance, on their surface, of one charge sign over another. The random charges of both sign derive from natural grain-to-grain interactions that produce triboelectrification via charge exchange every time grain surfaces make contact with one another. The conversion from a random distribution of grains (upon which there are randomly distributed charges) into an organization of electrostatically-ordered aggregates, can be regarded (within the framework of granular-material science) as an "electrical or Coulombic phase change" of the particulate cloud. It is not totally dissimilar from the more normal phase-change concept in which, for example, a gas with long free-path-molecules suddenly becomes a solid as a result of structural ordering of the molecules (notably, also the result of electronic forces, albeit at a different scale). In both the gas-to-solid case, and the aerosol-to-aggregate case, the same materials and charges are present before and after the phase change, but their arrangement now has a higher degree of order and a lower-energy configuration. An input of energy into the system is required to reverse the situation. The aggregates in the USML experiments were observed to undergo at least two phase changes as noted above. The point about phase changes, and by implication, the "electrostructural" reorganizations in particulate clouds, is the following: (a) they can occur very rapidly, almost spontaneously, above a critical cloud density, (b) in going from a higher energy state to a lower energy state, they convert to a denser system, (c) energy must be required to reverse the situation, implying that energy is released during the high-to-low energy phase change. In applying this information to natural particulate clouds, some inferences can be made (it is stressed that reference is still to dielectric materials attracted by dipole forces). There are several natural settings to which the USML observations apply, and to which the phase-change implications likewise apply. Dense clouds of triboelectrically-charged, kinetically-energized grains are to be found in volcanic eruptions (particularly on earth), aeolian dust storms (particularly on Mars), meteorite impact ejecta curtains (on all planets), in "immature" debris rings around planets (e.g., that from which our own Moon may have condensed), and in gravitationally collapsing protoplanetary dust/planetesimal debris disks where dielectric granules are being increasingly brought into collisional relationships with one another (increasing both charge exchange and physical proximity). It is noted that in many of these cases, the degree of electrical charging on the grains is likely to be much higher than that in the USML experiments where charging was not enhanced above the "normal", naturally encountered level for the particular materials (quartz, glass, and various silicate minerals). Application of the phase-change concept suggests that volcanic, aeolian, and impact debris clouds may, under certain circumstances, undergo rapid, impulsive, or even catastrophic collapse into a denser state that will lead to rapid precipitation or fall-out of suspended particulates. Although this idea has been suggested previously by the author , the phase change concept possibly permits some new insights into cloud-system behavior. For example, in a protoplanetary debris disk, the work of gravity may suddenly be enhanced by electrostatically driven collapse of the system when materials reach a critical intergranular spacing or grain number density. This might reduce the rate of planet formation by orders of magnitude (considering the ratio of g-forces to electrostatic ones for very small grains in close proximity), and indeed, it might drive the collapse into a system configuration that would perhaps not be created by gravity alone. Additional information is contained in the original.
Deploying phase change materials (PCMs) in building envelopes can be efficacious in reducing space heating/cooling loads and providing load shedding and shifting capacity. However, the full potential of PCM-integrated envelopes can only be harnessed if the PCM undergoes phase change using free ambient heating/cooling, and the stored energy is effectively transferred between the exterior and the interior environments. Traditional thermal insulation (with a fixed thermal resistance) limits PCM utilization, which restrains the energy saving potential of a PCM-integrated envelope to a small percentage. Proposed dynamic insulation material and system (DIMS) provides the option of varying its thermal resistance based on the indoor and outdoor conditions. Although it has been established that employing PCM as well as DIMS in building envelopes separately improves buildings’ energy performance, no prior studies that analyzed the combined influence of both technologies were identified. In this study, we examine a novel wall design, comprising a layer of PCM between two layers of DIMS. We note that the PCM-DIMS-integrated wall provides significantly higher energy saving potential than the DIMS-only integrated wall or the PCM-only integrated wall in all the climates and wall orientations analyzed in this study. Depending on the climate, the PCM-DIMS-integrated wall could provide 15–72% reduction in annual heat gain and 7–38% reduction in annual heat loss. The analysis presented in this study supports the need to develop scalable dynamic insulations combined with thermal energy storage systems for buildings.
Phase change material heat sinks have been recognized as an important tool in optimizing thermal control systems for space exploration vehicles and habitats that must deal with widely varying thermal loads and environments. In order to better focus technology investment in this arena, NASA has supported a trade study with the objective of identifying where the best potential pay-off can be found among identified aqueous and paraffin wax phase change materials and phase change material heat sink design approaches. The study used a representative exploration mission with well understood parameters to support the trade. Additional sensitivity studies were performed to ensure the applicability of study results across varying systems and destinations. Results from the study indicate that a water ice PCM heat sink has the potential to decrease the equivalent system mass of the mission s vehicle through a combination of a smaller heat sink and a slight 5% increase in radiator size or the addition of a lightweight heat pump. An evaluation of existing and emerging PCM heat sink technologies indicates that further significant mass savings should be achievable through continued development of those technologies. The largest mass savings may be realized by managing the location of the liquid and the solid in the heat sink to eliminate the melting and freezing pressure of wax and water, respectively, while also accommodating the high structural loads expected on future manned launch vehicles.
Phase change material heat sinks have been recognized as an important tool in optimizing thermal control systems for space exploration vehicles and habitats that must deal with widely varying thermal loads and environments. In order to better focus technology investment in this arena, NASA has supported a trade study with the objective of identifying where the best potential pay-off can be found among identified aqueous and paraffin wax phase change materials and phase change material heat sink design approaches. The study used a representative exploration mission with well understood parameters to support the trade. Additional sensitivity studies were performed to ensure the applicability of study results across varying systems and destinations. Results from the study indicate that replacing a wax PCM heat sink with a water ice PCM heat sink has the potential to decrease the equivalent system mass of the mission s vehicle through a combination of a smaller heat sink and a slight 5% increase in radiator size or the addition of a lightweight heat pump. An evaluation of existing and emerging PCM heat sink technologies indicates that further mass savings should be achievable through continued development of those technologies. The largest mass savings may be realized by eliminating the melting and freezing pressure of wax and water, respectively.
This white paper focuses on the scientific rationale and motivation for deploying phase change processes-based spacecraft thermal management systems for space missions [1,2]. Phase change processes need to be implemented in the design of various life support systems, sub-critical cryogenic and thermal control systems, taking advantage of the latent heat of phase change so as to significantly increase the heat transfer rate per unit fluid mass. Experimental data are critical to understanding gravitational effects on complex phase change processes using novel materials and working fluids. Experimental data can be acquired by either designing new testing platforms or performing experiments with simulant fluids on existing microgravity platforms. With the acquired data, theoretical flow physics models and numerical simulations can then be validated for use under long-duration microgravity conditions. The acquired data and validated models will enable the deployment of game-changing thermal management systems with remarkably improved efficiency and stability in heat acquisition, transport, and rejection for long-duration space missions.
Abstract The primary challenge of salt hydrates as phase change materials (PCMs) is their high degree of supercooling (∆ T Sup ). Past studies have attempted to alleviate ∆ T Sup by incorporating nucleators possessing similar lattice structure, without a thorough analysis of the nucleator/PCM interactions. This work presents a novel in‐operando study to visually observe the crystallization process of zinc nitrate hexahydrate (ZNH) in real‐time using suitable nucleators and a thickener. We have introduced a new dataset where zinc acetate dihydrate (ZnAc.2H) has been employed as an additive to ZNH for thermal energy storage (TES). According to our in‐operando study, ZnAc.2H underwent hydrolysis in molten ZNH to precipitate needle‐like ZnO particles, which acted as nucleators for the PCM. In the absence of a thickener, the crystal propagation of ZNH exhibited a preferential directional tendency, whereas incorporation of thickener made the phase change process uniform along all directions. Self‐hydrolysis of ZnAc.2H generated excess water which caused the undesirable effect of broadening the endothermic peak. The addition of carboxymethyl cellulose (CMC) thickener could restrict the activity of excess water and reduce the broadening of the endothermic peak. Differential scanning calorimetry (DSC) revealed that 5 wt% ZnAc.2H additive and 2 wt% CMC thickener added to ZNH can considerably improve the TES properties of ZNH PCM (∆ H Fusion 125.9 J·g −1 , ∆ T Sup 3.0°C, melting point [M.P.] 29.0°C). Our in‐operando studies can unveil real‐time phase change behaviors to better design PCM systems with desired characteristics.
Spacecraft solar dynamic power systems typically use high-temperature phase-change materials to efficiently store thermal energy for heat engine operation in orbital eclipse periods. Lithium fluoride salts are particularly well suited for this application because of their high heat of fusion, long-term stability, and appropriate melting point. Considerable attention has been focused on the development of thermal energy storage (TES) canisters that employ either pure lithium fluoride (LiF), with a melting point of 1121 K, or eutectic composition lithium-fluoride/calcium-difluoride (LiF-20CaF2), with a 1040 K melting point, as the phase-change material. Primary goals of TES canister development include maximizing the phase-change material melt fraction, minimizing the canister mass per unit of energy storage, and maximizing the phase-change material thermal charge/discharge rates within the limits posed by the container structure.
An important subclass of metal–organic frameworks (MOFs) exhibits cooperative flexibility, wherein individual crystallites undergo global structural phase changes in response to external stimuli. Where cooperative flexibility results in reversible changes between crystalline states of distinct accessible porosity, these frameworks can exhibit rare yet desirable behaviors that cannot be explained by local dynamics alone. Yet, the chemical and structural origins of cooperative flexibility and how frameworks undergo these reversible phase changes at the atomic level remain poorly understood. Deliberate design for specific applications is therefore exceedingly difficult, and there is great impetus to develop a fundamental understanding of this phenomenon. Here, an effective and widely accessible computational approach is developed, which is designed to provide microscopic resolution via direct comparison to experimental data along the desorption-guided pathway. The strategy is applied to explain the desorption-induced phase change in an experimentally well-characterized framework, CdIF-13 (sod-Cd(benzimidazolate)2), where experiment alone was unable to resolve the atomistically detailed phase change landscape. Our findings reveal that the cooperative phase change pathways are adsorbate dependent with thermodynamics of intermediate structural states dictated by a nuanced interplay of ligand orientation, skeletal symmetry, and modes of surface adsorption. The results reveal that this isotropically flexible framework is “chaperoned” through a complex energy landscape by specific adsorbates, revealed by the reported computational approach with atomic-level insight and validated by experimentally determined structures. Thus, this work facilitates both understanding and future design of flexible materials for applications in gas storage, transport, delivery, and separation technologies.
The transient response of a single cell of a typical phase change material type thermal capacitor has been modeled using numerical conductive heat transfer techniques. The cell consists of a base plate, an insulated top, and two vertical walls (fins) forming a two-dimensional cavity filled with a phase change material. Both explicit and implicit numerical formulations are outlined. A mixed explicit-implicit scheme which treats the fin implicity while treating the phase change material explicitly is discussed. A band algorithmic scheme is used to reduce computer storage requirements for the implicit approach while retaining a relatively fine grid. All formulations are presented in dimensionless form thereby enabling application to geometrically similar problems. Typical parametric results are graphically presented for the case of melting with constant heat input to the base of the cell.
The homogenous relaxation model (HRM) is one of the most widely used models to describe the liquid- gas phase transition. However, in its original formulation, it is unable to handle multispecies vapor-liquid equilibrium (VLE), which limits its applicability to single-component fluids. In this work, a unified non-equilibrium phase change model that considers the VLE of multicomponent mixtures is proposed building upon the HRM's structure. A time factor is introduced to mimic the effect of different phase change timescales due to different mechanisms, e.g., cavitation, flash-boiling, and evaporation. Here to assess the model's performance, computational fluid dynamics simulations of the internal and near-nozzle injection flow of the Engine Combustion Network's Spray G injector were performed using the nine-component PACE-20 fuel with both the unified model and the original HRM. The predicted fuel density in the near-nozzle region matched well with X-ray tomography measurements. The simulation results indicated that, whereas the HRM failed to capture the vaporization due to convective mixing between the fuel and ambient gas, the unified model performed well in predicting the mixing-driven vaporization and the corresponding evaporative cooling. Further comparisons using the nine-component fuel formula and a single-component fuel surrogate demonstrated the unified model's ability to predict preferential vaporization, which affects the predictions of local mixture composition and rate of vaporization. Finally, it is shown that the unified model is capable of representing multiple phase change mechanisms, and the relaxation time factor plays an important role in determining the degree of phase change due to the different mechanisms.