Future operations of many fluid, thermal and power systems and their ability to store, transfer, and manage a variety of single or multiphase fluids in reduced gravity environment are of great importance. For many of these systems, cryogenic conditions will play an important role. Cryogenic vaporization, caused by heat leakage into the tank from the surrounding environment, is one of the main causes of mass loss and leads to self-pressurization of the storage tanks. Available publications on self-pressurization and stratification of cryogenic tanks mainly focus on the convection and surface evaporation influences. Because large superheats increase the likelihood of evaporation in the liquid, the evaporation and its effect on vapor pressure under microgravity is studied numerically in this paper. The effects of reduced gravity, contact angle of the vapor bubble, and surface tension are investigated. The computations are carried out by using the CFD software package, Ansys Fluent, and an in-house developed code to calculate the source term associated with the phase change. A coupled level set and the volume-of-fluid method (CLSVOF) are used to solve a single set of conservation equations for the whole domain and the interface between the two phases is tracked or captured. A heat and mass transfer model is implemented into the Fluent code for solving problems involving evaporation or condensation. Results show that small tiny vapor regions caused by the evaporation process change the pressure rise. Vortices are observed due to the vapor dynamics.
Nomenclaturec p Specific heat at constant pressure J/(kg·K) D Diameter of the tank [m] d bw Departure diameter of a vapor bubble (m) E Energy [J] F Body force [N/m 3 ] f Vapor bubble departure frequency (1/s) g Gravitational acceleration [m/s 2 ] H Height of the tank [m] k Thermal conductivity [W/(m·K)] L Heated length [m] L H Latent heat [J/kg] n Normal direction [m], also number of active nucleation sites P Pressure [Pa] q w Wall heat flux [W/m 2 ] Ra Rayleigh number, 2 4 2 p w Ra g c q L k βρ µ = r Radial coordinate [m] M A N U S C R I P T A C C E P T E D ACCEPTED MANUSCRIPT 2 S h Source term in energy equation [W/m 3 ] s Width of the rib [m] T Temperature [K] t Time [s] t g Growth time [s] t w Waiting time [s] V r Velocity vector [m/s] x, y, z Coordinates [m]Greek symbols α Volume of fluid fraction [-] α Heat transfer coefficient [W/m 2 K]