CFD Modeling of the Multipurpose Hydrogen Test Bed (MHTB) Self-Pressurization and Spray Bar Mixing Experiments in Normal Gravity: Effect of the Accommodation Coefficient on the Tank Pressure

Abstract

A CFD model for simulating the self-pressurization of a large scale liquid hydrogen storage tank is utilized in this paper to model the MHTB self-pressurization experiment. The kinetics-based Schrage equation is used to account for the evaporative and condensi ng interfacial mass flows in this model. The effect of the accommodation coefficient for calculating the interfacial mass transfer rate on the tank pressure during tank selfpressurization is studied. The values of the accommodation coefficient which were considered in this study vary from 1.0e-3 to 1.0e-1 for the explicit VOF model and from 1.0e-4 to 1.0e-3 for the implicit VOF model. The ullage pressure evolutions are compared against experimental data. A CFD model for controlling pressure in cryogenic storage tanks by spraying cold liquid into the ullage is also presented. The Euler-Lagrange approach is utilized for tracking the spray droplets and for modeling the interaction between the droplets and the continuous phase (ullage). The spray model is coupled with the VOF model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. Droplet-ullage heat and mass transfer are modeled. The flow, temperature, and interfacial mass flux, as well as droplets trajectories, size distribution and temperatures predicted by the model are presented. The ul lage pressure and vapor temperature evolutions are compared with experimental data obtained from the MHTB spray bar mixing experiment. The effect of the accommodation coefficient for calculating the interfacial and droplet mass transfer rates on the tank pressure during mixing of the vapor using spray is studied. The values used for the accommodation coefficient at the interface vary from 1.0e-5 to 1.0e-2. The droplet accommodation coefficient values vary from 2.0e-6 to 1.0e-4.