Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase change”

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 1,009 records · Page 56

Zero-Boil-Off Tank (ZBOT) Experiment: Ground-Based Validation of Self-Pressurization and Pressure Control Two-Phase CFD Model

Integral to all phases of NASA's projected space and planetary expeditions is affordable and reliable cryogenic fluid storage for use in propellant or life support systems. Cryogen vaporization due to heat leaks into the tank from its surroundings and support structure can cause self-pressurization relieved through venting. This has led to a desire to develop innovative pressure control designs based on mixing of the bulk tank fluid together with some form of active or passive cooling to allow storage of the cryogenic fluid with zero or reduced boil-off. The Zero-Boil-Off Tank (ZBOT) Experiments are a series of small scale tank pressurization and pressure control experiments aboard the International Space Station (ISS) that use a transparent volatile simulant fluid in a transparent sealed tank to delineate various fundamental fluid flow, heat and mass transport, and phase change phenomena that control storage tank pressurization and pressure control in microgravity. The hardware for ZBOT-1 flew to ISS on the OA-7 flight in April 2017 and operations are planned to begin in September 2017, encompassing more than 90 tests. This paper presents preliminary results from ZBOT's ground-based research delineating both pressurization and pressure reduction trends in the sealed test tank. Tank self-pressurization tests are conducted under three modes: VJ heating, strip heating and simultaneous VJ and strip heating in attempt to simulate heat leaks from the environment, the support structure and both. The jet mixing pressure control studies are performed either from an elevated uniform temperature condition or from thermally stratified conditions following a self-pressurization run. Jet flow rates are varied from 2-25 cm/s spanning a range of jet Re number in laminar, transitional, and turbulent regimes and a range of Weber numbers covering no ullage penetration, partial penetration and complete ullage penetration and break-up (only in microgravity). Numerical prediction of a two-phase CFD model are compared to experimental 1g results to both validate the model and also indicate the effect of the residual non-condensable gas on evolution of pressure and temperature distributions in the tank during pressurization and pressure control.

microgravity evaporation condensatio↗

Novel thermal energy storage component: Development, performance, and phase transition diagnosis

Thermal energy storage (TES) using phase change materials (PCMs) is a promising technology for capturing and storing excess thermal energy for later use. However, challenges such as poor heat transfer efficiency and a lack of modular, scalable designs have limited widespread adoption of TES in real-world applications. This study developed and evaluated modular brick-type and blade-type TES prototypes featuring an aluminum housing, an embedded serpentine coil for active or passive thermal exchange, and a cost-effective metal mesh to enhance PCM thermal conductivity. The blade-type TES achieved notable geometric efficiency, with a thickness-to-length ratio of 0.03 and a thickness-to-width ratio of 0.08, enabling highly compact and modular thermal storage suitable for space-constrained applications. The paper presents a detailed evaluation of the TES prototypes’ performance. The comparative analysis indicated that the TES prototypes provide a highly cost-effective, thermally optimized alternative for compact energy storage and load shifting. A novel diagnostic technique was also introduced: using a portable endoscope to capture real-time visualizations of PCM phase transitions inside the TES. This method provides critical insights into internal heat transfer mechanisms, identifies potential issues, and offers valuable support for optimizing the TES design and developing the control algorithm. Overall, the modular brick-type and blade-type TES designs demonstrated in this work provide a scalable, efficient, and economically viable solution for advancing TES across residential, commercial, and industrial sectors. The designs’ compact structure, enhanced thermal performance, and integrated diagnostic capabilities make them strong candidates for future deployment in energy-efficient systems.

Gao, Zhiming [ORNL] (ORCID:0000000271397995)↗

Convective Melting of Particles in Flow Under Microgravity Conditions

Study of melting of dispersed or packed solid particles in a fluid under gravity and microgravity conditions provides benchmark information for many engineering applications such as material processing, environmental assessment and protection and space fire protection. During such processes, packed or dispersed solid particles are interacting with fluid flow at above-melting temperatures. By unmasking the buoyancy effects in coupled flow, phase change and heat transport phenomena, a better understanding of melting rate in non-thermal equilibrium, convective conditions can be studied. A series of flight experiments were conducted onboard the NASA KC-135 Microgravity Research Airplane. The Particle Melting in Plow (PMF) module was designed to allow flow through the initially packed ice particles at controlled temperature and velocity. To achieve this, a close-loop flow system was designed. Video images were taken to record the visualization of the melting process, from which a time variation of packed particle thickness distribution at different times can be obtained by the image analysis method. The fluid temperature distribution within the melting zone is measured by thermocouples. An infrared camera was mounted from the top of the test section to record the ice-water thermal images at a given location. The results from thermal images yield local temperature variation between melting solid and liquid and local Stephan number. Typical results for a number of cases are presented. The mathematical model, describing mass, energy, and momentum balance equations for the liquid and solid phases, is presented. It is found that melting rate is influenced mainly by the ratio of Reynolds number (based on the initial particle diameter) to the Froud number, and Stephan number. At the absence of gravity, Froud number approaches zero, Reynolds number and Stephan number become dominant factors governing the melting rate. The numerically determined results are compared with the experimental ones. It is found that the discrepancy between the predicted and measured melting rate is largely due to the inaccuracy in the constitutive equations for effective thermal physical properties, such as effective thermal conductivity and diffusivity, and transport properties, such as particle interaction coefficient, and local heat transfer coefficient of particles.

Stewart, C.↗

Improved Damage Resistant Composite Materials Incorporating Shape Memory Alloys

Metallic shape memory alloys (SMA) such as nitinol have unique shape recovery behavior and mechanical properties associated with a material phase change that have been used in a variety of sensing and actuation applications. Recent studies have shown that integrating nitinol-SMA actuators into composite materials increases the composite material's functionality. Hybrid composites of conventional graphite/epoxy or glass/epoxy and nitinol-SMA elements can perform functions in applications where monolithic composites perform inadequately. One such application is the use of hybrid composites to function both in load bearing and armor capacities. While monolithic composites with high strength-to-weight ratios function efficiently as loadbearing structures, because of their brittle nature, impact loading can cause significant catastrophic damage. Initial composite failure modes such as delamination and matrix cracking dissipate some impact energy, but when stress exceeds the composite's ultimate strength, fiber fracture and material perforation become dominant. One of the few methods that has been developed to reduce material perforation is hybridizing polymer matrix composites with tough kevlar or high modulus polyethynylene plies. The tough fibers increase the impact resistance and the stiffer and stronger graphite fibers carry the majority of the load. Similarly, by adding nitinol-SMA elements that absorb impact energy through the stress-induced martensitic phase transformation, the composites' impact perforation resistance can be greatly enhanced. The results of drop-weight and high velocity gas-gun impact testing of various composite materials will be presented. The results demonstrate that hybridizing composites with nitinol-SMA elements significantly increases perforation resistance compared to other traditional toughening elements. Inspection of the composite specimens at various stages of perforation by optical microscope illustrates the mechanisms by which perforation is initiated. Results suggest that the out-of-plane transverse shear properties of the composite and nitinol elements have a significant effect on the perforation resistance. Applications that can utilize the hybrid composites effectively will also be presented with the experimental studies.

Paine, Jeffrey S. N.↗

Modeling of Melt Growth During Carbothermal Processing of Lunar Regolith

The carbothermal processing of lunar regolith has been proposed as a means to produce carbon monoxide and ultimately oxygen to support human exploration of the moon. In this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon monoxide that bubbles out of the melt. Carbon monoxide is further processed in other reactors downstream to ultimately produce oxygen. The amount of oxygen produced crucially depends on the amount of regolith that is molten. In this paper we develop a model of the heat transfer in carbothermal processing. Regolith in a suitable container is heated by a heat flux at its surface such as by continuously shining a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and its temperature rises until it attains the melting point. The energy from the heat flux is then used for the latent heat necessary to change phase from solid to liquid, after which the temperature continues to rise. Thus a small melt pool appears under the heated zone shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat and supplies the energy required to melt more of the regolith, and the size of the molten zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the melt is balanced by radiative losses from the surface. In this paper, we model the melting and the growth of the melt zone with time in a bed of regolith when a portion of its surface is subjected to a constant heat flux. The heat flux is assumed to impinge on a circular area. Our model is based on an axisymmetric three-dimensional variation of the temperature field in the domain. Heat transfer occurs only by conduction, and effects of convective heat transport are assumed negligible. Radiative heat loss from the surface of the melt and the regolith to the surroundings is permitted. We perform numerical computations to determine the shape and the mass of the melt at steady state and its time evolution. We first neglect the volume change upon melting, and subsequently perform calculations including it. Predictions from our model are compared to test data to determine the effective thermal conductivities of the regolith and the melt that are compatible with the data

Balasubramaniam, R.↗

Engineering the fracture resistance of 2H-transition metal dichalcogenides using vacancies: An in-silico investigation based on HRTEM images

Vacancy engineering of 2H-transition metal dichalcogenides (2H-TMDs) has recently attracted great attention due to its potential to fine-tune the phonon and opto-electric properties of these materials. From a mechanical perspective, this symmetry-breaking process typically reduces the overall crack resistance of the material and adversely affects its reliability. However, vacancies can trigger the formation of heterogeneous phases that synergistically improve fracture properties. In this study, using MoSe2 as an example, we characterize the types and density of vacancies that can emerge under electron irradiation and quantify their effect on fracture. Molecular dynamic (MD) simulations, employing a re-parameterized Tersoff potential capable of accurately capturing bond dissociation and structural phase changes, reveal that isolated transition metal monovacancies or chalcogenide divacancies tend to arrest the crack tip and hence enhance the monolayer toughness. In contrast, isolated chalcogenide monovacancies do not significantly affect toughness. The investigation further reveals that selenium vacancy lines, formed by high electron dose rates, alter the crack propagating direction and lead to multiple crack kinking. Using atomic displacements and virial stresses together with a continuum mapping, displacement, strain, and stress fields are computed to extract mechanistic information, e.g., conditions for crack kinking and size effects in fracture events. The study also reveals the potential of specific defect patterns, “vacancy engineering,” to improve the toughness of 2H-TMDs materials.

36 MATERIALS SCIENCE↗

Thermocatalytic Heat Pipes for Geothermal Resource Recovery

Heat pipes are an important technology that allow orders of magnitude faster heat transfer than simple conduction. However, operating principles in heat pipes place fundamental bounds on their performance (critical heat flux and efficiency). Conventional heat pipe functionality is inherently tied to vaporization and condensation of the working fluid charged in the heat pipe. These fluids each have different operating temperature ranges based on the capillary, entrainment, sonic, and boiling limits of the heat pipe design. These limits, typically the capillary limit, dictate the maximum heat flux a heat pipe can carry, and most importantly for geothermal systems, the distance over which the pipes can operate (100 to 200 m maximum under optimum conditions). A thermocatalytic heat pipe breaks the inherent limitations of phase change thermo- and hydrodynamics and can transform heat pipe technology as a potentially more efficient means of extracting heat from a geothermal resource. The thermocatalytic heat pipe uses a working fluid to transport both sensible and chemical heat. An endothermic chemical reaction at depth removes heat from the reservoir and produces reactive intermediates, which are transported to the surface and used to run a reverse exothermic reaction that releases heat for use in power generation or other useful purposes. This technology offers two distinct advantages over conventional geothermal heat recovery technologies: (1) lower heat loss to the rock outside of the geothermal reservoir, and (2) higher heat transfer rates to the well field within the geothermal reservoir. Both advantages offer opportunity to reduce risks and lower costs of geothermal energy recovery. In this report, we discuss an initial effort to assess the efficacy and limitations of this technology for extracting heat from both porous/permeable and nominally impermeable geothermal reservoirs. Numerical simulation capabilities of the STOMP-GT code were enhanced to enable simulations of thermochemical heat pipes traversing geothermal reservoirs. An array of potential thermochemical reaction systems was evaluated and screened. Of these, an ethanol dehydration reaction was most promising in the vapor-liquid reaction set. A solid-phase dehydration reaction (CuSO4·5H2O) showed the highest reaction enthalpy per unit volume but would require development of a nonaqueous carrier fluid to implement it in a heat pipe. Subsurface reservoir simulations predicted long-term performance of the heat pipes for each geothermal reservoir type. The performance of U-shaped wells and coaxial wells was evaluated for a suite of reactions for both hydrothermal and hot dry rock reservoirs and was compared with a baseline case of simply pumping water through the wells. The heat pipe technology was additionally evaluated for an enhanced geothermal system (EGS) with an injection borehole, production borehole, and intervening hydraulically conductive fracture. All reservoir types showed significant improvement in heat recovered over a 20-year operating period ranging from a 1.8X increase for the hot dry rock case to more than 2.5X more energy recovered for the EGS case.

15 GEOTHERMAL ENERGY↗

Predictions of Line Chilldown Boiling Regime Transitions by a Coupled CFD-Sub-Grid Boiling Model Validated against 1G LN2 Experiments

Propellant storage and transfer during future long-duration missions will involve fuel depot operations in which a donor depot tank is used to fill a receiver spacecraft tank with a cryogenic propellant. Prior to the on-Orbit filling operation both the receiver tank and the transfer line must be cooled. The line chilldown process involves transition between boiling regimes in microgravity that will be quite different from their 1g ground-based counterparts. Since the cryogenic propellant itself will be used to perform the chilldown process, the time constants to cool the wall and the amount of fuel that will be used become important design considerations. In this light, the focus of the present work is to capture the flow boiling characteristics of a cryogenic fluid, namely, liquid nitrogen, during the chill-down of a transfer line using CFD modelling and simulations. The cryogenic chill-down process involves different flow boiling regimes: film boiling, transition boiling and nucleate boiling. The prediction of transition between these regimes in a CFD framework is new and challenging. The present work addresses this challenge by employing a volume-of-fluid (VOF) based methodology with Lee phase change model to predict the film boiling regime of the chill-down process in ANSYS Fluent®. The transition and nucleate boiling regimes are predicted by incorporating a sub-grid model that accounts for bubble nucleation, growth, shedding frequency, and departure diameter. The sub-grid model is implemented into Fluent via a user-defined function for wall-fluid heat flux calculations. The sub-grid model is similar in formulation to the well-known Rensselaer Polytechnic Institute (RPI) boiling model. The CFD model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Predicted results shows good agreement of wall temperature, rewetting temperature, and transition between film and nucleate boiling with the experimental measurements for several different LN2 flowrates in the vertical orientation. The strengths and weaknesses of this modeling approach are presented and discussed.

Cryogenic Fluid Management↗

Validation of CFD Model Prediction of Flow Boiling Regime Transitions during LN2 Line Chilldown

Propellant storage and transfer during future long-duration missions will involve fuel depot operations in which a donor depot tank is used to fill a receiver spacecraft tank with a cryogenic propellant. Prior to the on-Orbit filling operation both the receiver tank and the transfer line must be cooled. The line chilldown process involves transition between boiling regimes in microgravity that will be quite different from their 1g ground-based counterparts. Since the cryogenic propellant itself will be used to perform the chilldown process, the time constants to cool the wall and the amount of fuel that will be used become important design considerations. In this light, the focus of the present work is to understand the flow boiling characteristics of a cryogenic fluid, namely, liquid nitrogen, during the chill-down of a transfer line using CFD modelling and simulations. The cryogenic chill-down process involves different flow boiling regimes: film boiling, transition boiling and nucleate boiling. The prediction of transition between these regimes in a CFD framework is new and challenging. The present work addresses this challenge by employing a volume-of-fluid (VOF) based methodology with Lee phase change model to predict the film boiling regime of the chill-down process in ANSYS Fluent®. The transition and nucleate boiling regimes are predicted by incorporating a sub-grid model developed by Craft Tech that accounts for bubble nucleation, growth, shedding frequency, and departure diameter. The sub-grid model is implemented into Fluent via a user-defined function for wall-fluid heat flux calculations. The sub-grid model is similar in formulation to the well-known Rensselaer Polytechnic Institute (RPI) boiling model. The model constants are tested for different operating conditions and designated values are reported. The CFD model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Predicted results showing good agreement of wall temperature, rewetting temperature, and transition between film and nucleate boiling with the experimental measurements are presented and discussed.

Boiling regimes↗

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗

Transient Mesoscopic Immiscibility, Viscosity Anomaly, and High Internal Pressure at the Semiconductor–Metal Transition in Liquid Ga 2 Te 3

Gallium tellurides appear to be promising phase-change materials (PCMs) of the next generation for brain-inspired computing and reconfigurable optical metasurfaces. They are different from the benchmark PCMs because of sp 3 gallium hybridization in both cubic Ga 2 Te 3 and amorphous pulsed laser deposition (PLD) films. Liquid Ga 2 Te 3 also shows a viscosity η(T) anomaly just above melting when η(T) first increases and only then starts decreasing. We used high-energy X-ray diffraction to observe a transient mesoscopic immiscibility that suggested dense metallic liquid droplets in a semiconducting melt. The η(T) shape was consistent with this finding. A vanishing first sharp diffraction peak that also shifts to a higher Q indicates a high internal pressure in the metallic melt, which produces a remarkable asymmetry of the Ga–Te nearest neighbor distances and is reminiscent of high-pressure rhombohedral Ga 2 Te 3 . Finally, the observed phenomena provide a realistic scenario for a fast, multilevel SET–RESET response, which also unravels similar trends in the purported density-driven liquid polyamorphism of water, phosphorus, sulfur, and other materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CFD Predictions of Boiling Regime Transitions during Line Chilldown validated against a 1G LN2 Experiment

Introduction Before filling a propellant tank on the ground or in Space, the transfer line between the donor and receiver tanks must be cooled down preferably by sacrificing a minimum amount of the cryogenic fluid. The cryogenic line chill-down process involves a transition between different flow boiling regimes, namely, film boiling, transition film boiling, and nucleate boiling which are complex and may be quite gravity-dependent. Capturing these boiling phenomena and predicting the transition between them in a CFD framework is new and challenging both for 1g and microgravity applications. Materials & Methods The present work addresses this challenge by employing a two-phase Eulerian approach in the context of a homogeneous fluid mixture together with the Lee phase change model to capture the film boiling regime of the chill-down process using ANSYS Fluent®. The nucleate boiling regime is predicted by incorporating an in-house developed sub-grid model that accounts for bubble nucleation, bubble growth, bubble departure diameter, and their shedding frequency. The sub-grid model is encoded and implemented into Fluent via a user-defined function for the wall-fluid heat flux calculations. The mathematical formulation and numerical implementation of the CFD model are described in detail. The coupled CFD-Subgrid model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Results Numerical simulation results show good agreements between the CFD predictions of the wall temperature evolution, rewetting temperature, and transition between film and nucleate boiling, with the experimental measurements published by Darr et al [2] for several different LN2 flowrates in the vertical pipe orientation. The CFD predictions for the wall temperature distribution indicate a rapid quenching of the wall at two upstream and downstream temperature sensing locations as compared to the experimental measurement. The only tuning parameter in the CFD model is the Lee mass transfer coefficient. The CFD Model predicts the Liedenfrost rewetting temperature in close agreement with the experiment. This marks a transition between stable and transitionary flow boiling regimes. The CFD-predicted boiling curve for the downstream sensor location is also compared against its experimental counterpart and indicates that the model is able to predict all the key temperature and heat flux parameters during the transitions from stable to transitionary film boiling to nucleate boiling in close agreement with the experiment. A sequence of predicted volume fraction, and temperature contours depicting these transitions will be presented.

Evaporation Condensation↗

VALENTInE: A Concept for a New Frontiers Class Long Duration In-Situ Balloon Mission to Venus

Venus and Earth are similar in bulk composition, size, density, and approximate distance from the Sun, yet Venus’s modern-day climate and surface geology is distinctly different [1]. Previous missions to Venus revealed unusual volcanic features, possible continental crust, widespread volcanic plains, a weak magnetic field [2], and insights into the Venusian atmosphere [3,4]. Unfortunately, data from these missions were limited in spatial and temporal resolution and global extent. A future mission to Venus is critical to address fundamental questions surrounding the chemical composition and dynamics of the Venusian atmosphere [5], its geologic history [6,7], its internal structure [8,9], and its habitability throughout time [10,11]. We present the Venus Air and Land Expedition: a Novel Trailblazer for In situ Exploration (VALENTInE) mission to meet this need. VALENTInE is a variable altitude balloon that will passively float in Venus’s atmosphere between 45 and 55 km altitude. VALENTInE will acquire atmospheric data at varying latitudes and longitudes in addition to mapping the surface geomorphology and mineralogy across multiple terrains. Mission Objectives: The VALENTInE mission concept is driven by four main science objectives: 1. Determine whether the driving force of the superrotation of Venus’s atmosphere is caused by horizontal or vertical momentum transport. 2. Determine whether the atmospheric composition and noble gas inventory of the Venusian atmosphere is a product of outgassing from the initial protoplanetary source or if there are significant contributions from exogenic sources. 3. Determine whether the tesserae regions (particularly Aphrodite Terra) are felsic and relatively older than surrounding regions. 4. Determine if there is any evidence of a recent dynamo preserved in the rock record of Venus. Mission Overview: A balloon architecture provides a robust system that can survive long-term in the Venusian environment while taking accurate measurements of the lower cloud deck and surface. Two prior balloon missions to Venus, VEGA 1 and VEGA 2 in 1984, have demonstrated the potential for such planetary exploration; however, these missions were short-lived (46 hr), with limited range (54 km). Our spacecraft concept consists of a battery powered balloon with a gondola and flyby carrier stage (Fig. 1). Our novel design allows for at least 15 days of atmospheric exploration, including multiple ascents and descents in the Venus atmosphere. VALENTInE is baselined to launch in 2032, cruise on solar power for 128 days, enter, descend, and inflate (EDI) into Venus’s atmosphere above Aphrodite Terra, and float in Venus’s atmosphere for a total nominal mission duration of 15 days. The balloon, able to control its altitude by changing its buoyancy, will make one full circumnavigation every 4–8 days and move poleward ~1 o latitude per day; it will be passively directed by the horizontal air currents on Venus. The balloon can only be controlled in the z direction, and the expected latitude range is ±10° from EDI. The balloon itself will be a tracer for atmospheric structures such as zonal winds [5, 12]. During the 15-day operational period, there will be five dips to 45 km for a compositional study of the lower atmosphere and geological and magnetic mapping of the surface. Dips are used to take images and measurements closer to the surface and to obtain in situ vertical atmospheric profiles between 45 and 55 km. Instruments: The VALENTInE instrument payload will allow for extensive study of the geology, atmosphere, and interior of Venus. The instrument suite consists of six instruments. The mission profile for each instrument is shown in Fig. 2. Lower Atmosphere Meteorology Analyzer (LAMA) is an atmospheric structure instrument consisting of a thermometer, barometer, and accelerometer to continuously measure temperature and pressure as a function of altitude, longitude, and latitude. TracE and Noble Gas Investigator (TENGI) is aneutral mass spectrometer used to sample the dense atmosphere. TENGI will operate continuously at 45 km and 55 km and will measure D/H ratios, as well as Ne, Ar, and O isotopic ratios. Kilometer Scale Spectral Imager (KSSI) is amultispectral imager (850-1150 nm) used to resolve surface features on the order of km. Near InfraRed Multispectral Photometer (NIRMP) is an Infrared (IR) photometer used to image the surface through the clouds to characterize mineral assemblages at five different areas of the surface. NIRMP will be operated during dwell. ELevation REConnaissance (ELREC) is a light detection and ranging (LIDAR) instrument used for measuring topography at five different areas of the surface. The ELREC data will be combined with those of NIRMP to determine mineral assemblages and how they correlate to topography. Magnetic Exploration and Interior Detective (MEID) is a magnetometer used to detect any near-surface magnetic anomalies. The magnetometer will continuously operate at all altitudes. Mission Design: The spacecraft’s propulsion system will have a wet mass of 930 kg and will be launched on an Atlas V rocket. The spacecraft will jettison the payload upon arrival in Venus’s atmosphere before decelerating to orbital velocity. Figure 2. Notional power-dip profile while the balloon is on the dayside. Risks and Challenges: Raising and lowering the spacecraft to dip beneath the haze layer require a large amount of energy. This issue was partly compensated by limiting the number of dips over the 15 day prime mission. Power needs for pumping helium are reduced by slowing the descent speed, which consequently increases spatial resolution of in situ measurements between 45 and 55 km. The main limitation for the duration of the mission is the need to carry 15 days’ worth of batteries. The bus experiences external Venusian temperatures ranging from ~27°C at 55 km altitude to ~110°C at 45 km altitude. These high atmospheric temperatures require thermal protection for the bus (instruments, electronics, and flight systems) during its descent, dwell at 45 km, and ascent. The bus is maintained at a mechanically safe temperature range of -10°C to 50°C using white external paint, multi-layer insulation (MLI), mechanical/thermal isolation (e.g. Ti, composites), and ~57 kg of Phase Change Materials (PCM). Images and spectra taken below the haze dominate the available data transmission regardless of the time spent at 45 km. Therefore, less time at the lower altitude mainly reduces the coverage of in situ measurements there. Further, we assumed an unrealistically low-density material for the helium storage tanks on the gondola. However, this issue is partially resolved if we were to jettison 75% of the spent helium storage tanks. After the entry process, the balloon volume remains mostly inflated and storage tanks are required only for reducing balloon volume in dipping to 45 km. This mission was planned against an expected New Frontiers 5 (NF5) cost cap of $1B, as the NF5 call had not yet been released. The mission we describe fits within the predicted cost cap if the flyby carrier stage, responsible for powering the spacecraft during cruise and for relaying in situ measurements back to Earth, can be contributed by another space agency. Acknowledgments: We thank the JPL Planetary Science Summer School, especially A. Nash, J. Scully, K. Mitchell, L. Lowes, and J. Armijo, and our mentors from JPL Team X. Thank you to our review panel for their insightful review. References: [1] Kane, S. et al. (2019) JGR:Ps, 124 , 2015–2028. [2] Phillips, J. L., & Russell, C. T. (1987). JGR: Space Physics , 92 (A3), 2253-2263. [3] Nakamura, M. et al. (2018) Earth Planets Space, 70 (1), 144. [4] Svedhem, H. et al. (2009) JGR:P 114 (E5). [5] Limaye et al. (2009) Decadal Survey White Paper . [6] Ivanov, M. and Head, J. (2011) PSS, 59 (13) , 1559-1600. [7] Smrekar, S. et al. (2018) Space Sci Rev, 214 (5), 88. [8] O’Rourke, J. et al. (2018) EPSL, 502 , 46-56. [9] O’Rourke J. et al. (2019) GRL , (46), 5768–5777 [10] Way, M. et al. (2016) GRL, 43 (16), 8376-8383. [11] Way, M. et al. (2020) JGR:P, 125 (5), e2019JE006276. [12] Preston, R. A. et al. (1986) Science, 231 (4744), 1414-1416.

Mission Concept↗

Application of Ibuprofen Sodium Dihydrate for Thermochemical Energy Storage

Thermochemical energy storage (TCES) offers a transformative approach to address grid instability by harnessing reversible chemical reactions for efficient heat storage and release. Here, we introduce pharmaceutical organic salt hydrates as a class of materials with exceptional performance for low-grade waste heat recovery. We demonstrate ibuprofen sodium dihydrate (ISD) as an example organic hydrate exhibiting a dehydration temperature range of 60-110 °C and a remarkable dehydration enthalpy of up to 59.5 kJ/mol of water, ideally suited for capturing industrial and residential waste heat. Using rigorous multimodal characterization, including thermogravimetric analysis, differential scanning calorimetry, in-situ FTIR, in-situ PXRD, and NMR, we demonstrate ISD's superior thermal, chemical, and structural stability over 150 hydration-dehydration cycles, achieving an unprecedented cycling efficiency of ~99.9%. Compared to conventional inorganic salt hydrates like strontium chloride hexahydrate and calcium oxalate monohydrate, ISD showcases enhanced durability without deliquescence or pulverization, even under high-humidity conditions. In-situ analyses confirm the transition from ISD to ibuprofen sodium anhydrous (ISA) proceeds with structural reorganization, thereby combining the dehydration mechanism with phase transitions, resulting in higher energy storage capacity. Microstructural analyses reveal that repeated water intercalation and structural transitions aid in creating significant porosity that enhances water transport kinetics, further improving the hydration/dehydration performance. By combining the phase change and chemical dehydration mechanisms, ISD paves the way for designing a new class of organic salt hydrates, offering tunable properties to meet diverse thermal energy storage demands and supporting sustainable grid resilience.

Thangaraj, Kavin C.↗

In‐operando crystallization study of zinc nitrate hexahydrate using zinc oxide nucleators

Abstract Supercooling (∆T Sup ) is one of the crucial problems in utilizing a phase change material (PCM), which has been attempted to overcome using nucleators possessing small lattice disregistry without comprehensive understanding of the nucleation phenomenon. Here, this work studies the interactions between nucleators and PCM via in‐operando direct visualization for the first time, to the best of our knowledge, to better understand the crystallization process during freezing of a PCM using zinc nitrate hexahydrate (ZNH) and a zinc oxide (ZnO) nucleator. According to our in‐operando study, freezing was randomly initiated by only a fraction of the nucleator particles during each thermal cycling. However, previously unresponsive nucleators also suddenly crystalized PCM when the propagating crystal encountered them, improving the crystal initiation and thereby ∆T Sup . Consequently, for equal wt% of ZnO, better nucleation behavior was obtained using uniformly distributed small nucleator particles throughout the PCM, as compared to either large or poorly distributed aggregated nucleators. DSC results using ZnO‐needle as nucleator confirmed a 38% and 33.3% lower ∆T Sup upon employing 5 wt% small (individual) particles (~4 μm on average), as compared to using equal wt% of small (aggregated) particle clusters (~200 μm) and large particles (~46 μm), respectively. Crystallization of ZNH caused additional hair‐like ZnO growth preferentially decorated along lateral faces of the original nucleator particles. The in‐operando studies are valuable tools to correlate inherent crystallization phenomena to the practical thermal energy storage properties of the system.

Chakraborty, Anirban↗

On the Pressure Generated by Thermite Reactions Using Stress-Altered Aluminum Particles

Here we examine pressure build-up and decay in thermites upon impact ignition and interprets reactivity based on the holistic pressure history. The thermite is a mixture of aluminum (Al) combined with bismuth trioxide (Bi 2 O 3 ) powder. Four different Al particles sizes were examined that ranged from 100 nm to 18.5 μm mean diameter and for each size, two different Al powder treatments were examined: stress-altered compared to untreated, as-received Al powder. Stress-altered Al powders have been shown to be more reactive, such that the stress-altered Al powder thermites offer a metric for analyzing thermite reactivity in terms of pressure development compared to untreated Al powder. In a binary thermite system, multiple phase changes and interface chemistry influence the transient pressure response during reaction. Results reveal three key pressure metrics that need consideration specifically for thermite combustion: (1) delay time to peak pressure, (2) peak pressure, and (3) decay after peak pressure. Our experiments show that a lower peak pressure corresponds with higher thermite reactivity because aluminum consumption of oxygen generated by decomposing solid oxidizer reduces the peak pressure. Faster rates of reaction consume oxygen at higher rates such that pressure development becomes more limited than less reactive thermites and the result is a lower peak pressure. This conclusion is opposite of traditional studies using metal fuels with a gaseous environment that typically show higher peak pressures correspond with greater reactivity.

42 ENGINEERING↗