Research study for determination of liquid surface profile in a cryogenic tank during gas injection Quarterly progress report no. 9, Jun. 18 - Sep. 17, 1966
Determining liquid surface profiles in cryogenic tank during gas injection
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.
Determining liquid surface profiles in cryogenic tank during gas injection
Static tests were made on the compatibility of liquid fluorine with several nonmetallic materials at -3200 F and at pressures of 0 and 1500 pounds per square inch gage. The results are compared with those from previous work with gaseous fluorine at the same pressures, but at atmospheric temperature. In general, although environmental effects were not always consistent, reactivity was least with the low-temperature, low-pressure liquid fluorine. Reactivity was greatest with the warm, high-pressure gaseous fluorine. None of the liquids and greases tested was found to be entirely suitable for use in fluorine systems. Polytrifluorochloroethylene and N-43, the formula for which is (C4F9)3N, did not react with liquid fluorine at atmospheric pressure or 1500 pounds per square inch gage under static conditions, but they did react when injected into liquid fluorine at 1500 pounds per square inch gage; they also reacted with gaseous fluorine at 1500 pounds per square inch gage. While water did not react with liquid fluorine at 1500 pounds per square inch gage, it is known to react violently with fluorine under other conditions. The pipe-thread lubricant Q-Seal did not react with liquid fluorine, but did react with gaseous fluorine at 1500 pounds per square inch gage. Of the solids, ruby (Al2O3) and Teflon did not react under the test conditions. The results show that the compatibility of fluorine with nonmetals depends on the state of the fluorine and the system design.
A computer program was developed to predict the trajectory, ground deposition, and drift of liquid sprays injected into the wake of an agricultural aircraft in ground effect. The program uses a horseshoe vortex wake model and includes the effects of liquid droplet evaporation, crosswind, the propeller slipstream, ground effect, and tunnel walls on small scale models. This user's guide includes several case examples demonstrating user options. A complete listing of the FORTRAN program is provided.
A computer model has been developed to analyze the three-dimensional two-phase reactive flows in liquid fueled rocket combustors. The model is designed to study the influence of liquid propellant injection nonuniformities on the flow pattern, combustion and heat transfer within the combustor. The Eulerian-Lagrangian approach for simulating polidisperse spray flow, evaporation and combustion has been used. Full coupling between the phases is accounted for. A nonorthogonal, body fitted coordinate system along with a conservative control volume formulation is employed. The physical models built into the model include a kappa-epsilon turbulence model, a two-step chemical reaction, and the six-flux radiation model. Semiempirical models are used to describe all interphase coupling terms as well as chemical reaction rates. The purpose of this study was to demonstrate an analytical capability to predict the effects of reactant injection nonuniformities (injection anomalies) on combustion and heat transfer within the rocket combustion chamber. The results show promising application of the model to comprehensive modeling of liquid propellant rocket engines.
Compressors can be described as the heart of a vapor compression refrigeration system, circulating the working fluid in the system. Significant benefits can be gained from improving its efficiency since energy consumption is responsible for almost the entire system's power input. An effort has been made to achieve isothermal compression through liquid refrigerant injection or inter-stage cooling in refrigeration systems. In recent years, much effort has been invested into isothermal compression technology for air compressors or compressed air energy storage systems with the rise of renewable energy. Here, this work has explored the advancements made in this area and categorized them to discuss the benefits, limitations, and tradeoffs. Evidently, increasing heat transfer area is a significant aid in improving heat transfer, resulting in better isothermal compression. However, liquid pistons open more avenues for this to be pursued, and therefore, are limited to air compression applications. In addition, water injection proved to be the most effective among the methods. Finally, there is a need for more experimental research on refrigeration applications and the field is still in its infancy.
A computational fluid dynamics (CFD) model is developed to simulate pressure control of an ellipsoidal-shaped liquid hydrogen tank under external heating in low gravity. Pressure control is provided by an axial jet thermodynamic vent system (TVS) centered within the vessel that injects cooler liquid into the tank, mixing the contents and reducing tank pressure. The two-phase cryogenic tank model considers liquid hydrogen in its own vapor with liquid density varying with temperature only and a fully compressible ullage. The axisymmetric model is developed using a custom version of the commercially available FLOW-3D software and simulates low gravity extrapolations of engineering checkout tests performed at Marshall Space Flight Center in 1999 in support of the Solar Thermal Upper Stage Technology Demonstrator (STUSTD) program. Model results illustrate that stable low gravity liquid-gas interfaces are maintained during all phases of the pressure control cycle. Steady and relatively smooth ullage pressurization rates are predicted. This work advances current low gravity CFD modeling capabilities for cryogenic pressure control and aids the development of a low cost CFD-based design process for space hardware.
A computational fluid dynamics (CFD) model is developed to simulate pressure control of an ellipsoidal-shaped liquid hydrogen tank under external heating in normal gravity. Pressure control is provided by an axial jet thermodynamic vent system (TVS) centered within the vessel that injects cooler liquid into the tank, mixing the contents and reducing tank pressure. The two-phase cryogenic tank model considers liquid hydrogen in its own vapor with liquid density varying with temperature only and a fully compressible ullage. The axisymmetric model is developed using a custom version of the commercially available FLOW-31) software. Quantitative model validation is ,provided by engineering checkout tests performed at Marshall Space Flight Center in 1999 in support of the Solar Thermal Upper Stage_ Technology Demonstrator (STUSTD) program. The engineering checkout tests provide cryogenic tank self-pressurization test data at various heat leaks and tank fill levels. The predicted self-pressurization rates, ullage and liquid temperatures at discrete locations within the STUSTD tank are in good agreement with test data. The work presented here advances current CFD modeling capabilities for cryogenic pressure control and helps develop a low cost CFD-based design process for space hardware.
Combustion instability, and secondary injection studies in liquid propulsion technology
This paper highlights results from the first of a series of hierarchical simulations aimed at assessing the modeling requirements for application of the large eddy simulation technique to cryogenic injection and combustion processes in liquid rocket engines. The focus is on liquid-oxygen-hydrogen coaxial injectors at a condition where the liquid-oxygen is injected at a subcritical temperature into a supercritical environment. For this situation a diffusion dominated mode of combustion occurs in the presence of exceedingly large thermophysical property gradients. Though continuous, these gradients approach the behavior of a contact discontinuity. Significant real gas effects and transport anomalies coexist locally in colder regions of the flow, with ideal gas and transport characteristics occurring within the flame zone. The current focal point is on the interfacial region between the liquid-oxygen core and the coaxial hydrogen jet where the flame anchors itself.
This paper highlights results from the first of a series of hierarchical simulations aimed at assessing the modeling requirements for application of the large eddy simulation technique to cryogenic injection and combustion processes in liquid rocket engines. The focus is on liquid-oxygen-hydrogen coaxial injectors at a condition where the liquid-oxygen is injected at a subcritical temperature into a supercritical environment. For this situation a diffusion dominated mode of combustion occurs in the presence of exceedingly large thermophysical property gradients. Though continuous, these gradients approach the behavior of a contact discontinuity. Significant real gas effects and transport anomalies coexist locally in colder regions of the flow, with ideal gas and transport characteristics occurring within the flame zone. The current focal point is on the interfacial region between the liquid-oxygen core and the coaxial hydrogen jet where the flame anchors itself.
A theoretical study of the breakup of a circular liquid jet injected transversely into a supersonic air stream is conducted. Two different criteria for breakup are explored in the context of a previously developed model for the behavior of liquid jets in compressible crossflow (Heister et al., 1989). The local sonic point criterion first proposed by Schetz, et al. (1980) is explored, in addition to an auxiliary criterion put forth by Clark (1964) based on surface-tension stability. It is found that the local sonic point appears to provide a more reasonable approximation to the actual location of jet breakup, based on comparisons with limited experimental data.
Group IVA is a large magmatic group of iron meteorites. The mean DELTA O-17 (= delta O-17 - 0.52(raised dot) delta O-18) of the silicates is approx. plus or minus 1.2%o, similar to the highest values in L chondrites and the lowest values in LL chondrites; delta O-18 values are also in the L/LL range. This strongly suggests that IVA irons formed by melting L-LL parental material, but the mean Ni content of IVA irons (83 mg/g) is much lower than that of a presumed L-LL parent (approx. 170 mg/g) and the low-Ca pyroxene present in two IVA meteorites is Fs13, much lower than the Fs20-29 values in L and LL chondrites. Thus, formation from L-LL precursors requires extensive addition of metallic Fe, probably produced by reduction of FeS and FeO. Group IVA also has S/Ni, Ga/Ni, and Ge/Ni ratios that are much lower than those in L-LL chondrites or any chondrite group that preserves nebular compositions, implying loss of these volatile elements during asteroidal processing. We suggest that these reduction and loss processes occurred near the surface of the asteroid during impact heating, and resulted partly from reduction by C, and partly from the thermal dissociation of FeS and FeO with loss of O and S. The hot (approx. 1770 K) low-viscosity melt quickly moved through channels in the porous asteroid to form a core. Two members of the IVA group, Sao Joao Nepomuceno (hereafter, SJN) and Steinbach, contain moderate amounts of orthopyroxene and silica, and minor amounts of low-Ca clinopyroxene. Even though SJN formed after approx. 26% crystallization and Steinbach formed after approx. 77% Crystallization of the IVA core, both could have originated within several tens of meters of the core-mantle interface if 99% of the crystallization occurred from the center outwards. Two other members of the group (Gibeon and Bishop Canyon) contain tabular tridymite, which we infer to have initially formed as veins deposited from a cooling SiO-rich vapor. The silicates were clearly introduced into IVA irons after the initial magma crystallized. Because the y-iron crystals in SJN are typically about 5 cm across, an order of magnitude smaller than in IVA irons that do not contain massive silicates, we infer that the metal was in the gamma-iron field when the silicates were injected. The SJN and Steinbach silicate compositions are near the low-Ca-pyroxene/silica eutectic compositions. We suggest that a tectonic event produced a eutectic-like liquid and injected it together with unmelted pyroxene grains into fissures in the solid metal core. Published estimates of IVA metallographic cooling rates range from 20 to 3000 K/Ma, leading to a hypothesized breakup of the core during a major impact followed by scrambling of the core and mantle debris [Haack, H., Scott, E.R.D., Love, S.G., Brearley, A. 1996. Thermal histories of IVA stony-iron and iron meteorites: evidence for asteroid fragmentation and reaccretion. Geochim. Cosmochim. Acta 60, 3103-3113]. This scrambling model is physically implausible and cannot explain the strong correlation of estimated cooling rates with metal composition. Previous workers concluded that the low-Ca clinopyroxene in SJN and Steinbach formed from protopyroxene by quenching at a cooling rate of 10(sup 12) K/Ma, and suggested that this also supported an impact-scrambling model. This implausible spike in cooling rate by a factor of 10(sup 10) can be avoided if the low-Ca clinopyroxene were formed by a late shock event that converted orthopyroxene to clinopyroxene followed by minimal growth in the clinopyroxene field, probably because melt was also produced. We suggest that metallographic cooling-rate estimates (e.g., based on island taenite) giving similar values throughout the metal compositional range are more plausible, and that the IVA parent asteroid can be modeled by monotonic cooling followed by a high-temperature impact event that introduced silicates into the metal and a low-temperature impact event that partially converted thopyroxene into low-Ca clinopyroxene.
Atomization of single liquid jets injected downstream in high pressure and high velocity airflow was investigated to determine the effect of airstream pressure on mean drop size as measured with a scanning radiometer. For aerodynamic - wave breakup of liquid jets, the ratio of orifice diameter D sub o to measured mean drop diameter D sub m which is assumed equal to D sub 32 or Sauter mean diameter, was correlated with the product of the Weber and Reynolds numbers WeRe and the dimensionless group G1/square root of c, where G is the gravitational acceleration, 1 the mean free molecular path, and square root of C the root mean square velocity, as follows; D sub o/D sub 32 = 1.2 (WeRe) to the 0.4 (G1/square root of c) to the 0.15 for values of WeRe 1 million and an airstream pressure range of 0.10 to 2.10 MPa.
Atomization of single liquid jets injected downstream in high pressure and high velocity airflow was investigated to determine the effect of airstream pressure on mean drop size as measured with a scanning radiometer. For aerodynamic - wave breakup of liquid jets, the ratio of orifice diameter D sub o to measured mean drop diameter D sub m which is assumed equal to D sub 32 or Sauter mean diameter, was correlated with the product of the Weber and Reynolds numbers WeRe and the dimensionless group G1/square root of c, where G is the gravitational acceleration, 1 the mean free molecular path, and square root of C the root mean square velocity, as follows; D sub o/D sub 32 = 1.2 (WeRe) to the 0.4 (G1/square root of c) to the 0.15 for values of WeRe 1 million and an airstream pressure range of 0.10 to 2.10 MPa.
The co-optimization of in-cylinder combustion and after-treatment technology has become a major aspect in engine design and development, with the goal of meeting the increasingly restrictive emission regulations in the transportation industry. Selective Catalytic Reduction is a robust technology to control the emission of NO x , and the injection of urea in water solution is the exhaust tailpipe is a key aspect of its operation. The proposed work uses high-fidelity Computational Fluid Dynamics to characterize the atomization dynamics of the liquid jet in relevant cross-flow conditions. The study focuses on a commercial low-pressure (9 bar) pressure-swirl injector which is characterized in its internal geometry through high-resolution X-ray micro-computational tomography. The internal two-phase flow has been modeled according to the volume-of-fluid approach in a large eddy simulation framework and validated against near-nozzle X-ray radiography measurement. Moreover, characterizing the breakup dynamics for the swirling hollow cone formation, and assessing the influence of the cross-flow in the breakup dynamics was completed. The results have been reported proposing Re-Oh maps and probability density functions of the spray kinematics. Higher cross-flow momentum generates an increase in the jet intact length and a reduction of the liquid droplet diameters. The axial momentum of the jet is affected by the cross-flow already in the near-nozzle region, determining a relevant deviation of the spray velocities. In conclusion, this work aims to inform the initialization of Eulerian-Lagrangian spray models through the assignment of droplet kinematics and static one-way coupling between volume-of-fluid results and Lagrangian spray parcels, to be used for system-size domain simulations.
The Advanced Rocket Injector/Combustor Code (ARICC) that has been developed to model the complete chemical/fluid/thermal processes occurring inside rocket combustion chambers are highlighted. The code, derived from the CONCHAS-SPRAY code originally developed at Los Alamos National Laboratory incorporates powerful features such as the ability to model complex injector combustion chamber geometries, Lagrangian tracking of droplets, full chemical equilibrium and kinetic reactions for multiple species, a fractional volume of fluid (VOF) description of liquid jet injection in addition to the gaseous phase fluid dynamics, and turbulent mass, energy, and momentum transport. Atomization and droplet dynamic models from earlier generation codes are transplated into the present code. Currently, ARICC is specialized for liquid oxygen/hydrogen propellants, although other fuel/oxidizer pairs can be easily substituted.
Mixed mode propulsion concepts are studied for advanced, single stage earth orbital transportation systems (SSTO) for use in the post-1990 time period. These propulsion concepts are based on the sequential and/or parallel use of high density impulse and high specific impulse propellants in a single stage to increase vehicle performance and reduce dry weight. Specifically, the mixed mode concept utilizes two propulsion systems with two different fuels (mode 1 and mode 2) with liquid oxygen as a common oxidizer. Mode 1 engines would burn a high bulk density fuel for lift-off and early ascent to minimize performance penalties associated with carrying fuel tankage to orbit. Mode 2 engines will complete orbital injection utilizing liquid hydrogen as the fuel.
A series of multispecies, multiphase computational fluid dynamics (CFD) analyses of the 24-inch diameter joint government industry industrial research and development (JIRAD) hybrid rocket motor is described. The 24-inch JIRAD hybrid motor operates by injection of liquid oxygen (LOX) into a vaporization plenum chamber upstream of ports in the hydroxyl-terminated polybutadiene (HTPB) solid fuel. The injector spray pattern had a strong influence on combustion stability of the JIRAD motor so a CFD study was initiated to define the injector end flow field under different oxidizer spray patterns and operating conditions. By using CFD to gain a clear picture of the flow field and temperature distribution within the JIRAD motor, it is hoped that the fundamental mechanisms of hybrid combustion instability may be identified and then suppressed by simple alterations to the oxidizer injection parameters such as injection angle and velocity. The simulations in this study were carried out using the General Algorithm for Analysis of Combustion SYstems (GALACSY) multiphase combustion codes. GALACSY consists of a comprehensive set of droplet dynamic submodels (atomization, evaporation, etc.) and a computationally efficient hydrocarbon chemistry package built around a robust Navier-Stokes solver optimized for low Mach number flows. Lagrangian tracking of dispersed particles describes a closely coupled spray phase. The CFD cases described in this paper represent various levels of simplification of the problem. They include: (A) gaseous oxygen with combusting fuel vapor blowing off the walls at various oxidizer injection angles and velocities, (B) gaseous oxygen with combusting fuel vapor blowing off the walls, and (C) liquid oxygen with combusting fuel vapor blowing off the walls. The study used an axisymmetric model and the results indicate that the injector design significantly effects the flow field in the injector end of the motor. Markedly different recirculation patterns are observed in the vaporization chamber as the oxygen velocity and/or spray pattern is varied. The ability of these recirculation patterns to stabilize the diffusion flame above the surface of the solid fuel gives a plausible explanation for the experimentally determined combustion stability characteristics of the JIRAD motor, and suggests how combustion stability can be assured by modifications to the injector design.