Results for “Jump conditions” — Page 6

NASA NTRS1992

Acceleration of positrons in supernova shocks

During this project we investigated the acceleration of leptons (electrons and positrons) in collisionless shock waves. In particular, we were interested in how leptons are accelerated in the blast waves existing in the remnants of supernova explosions. Supernova remnants (SNRs) have long been considered as the most likely source of galactic cosmic rays but no definite connection between SNRs and the cosmic rays seen at earth can be made. Only by understanding lepton acceleration in shocks can the rich SNR data base be properly used to understand cosmic ray origins. Our project was directed at the neglected aspects of lepton acceleration. We showed that the efficiency of lepton acceleration depended critically on the lepton injection energy. We showed that, even when infection effects are not important, that proton and lepton distribution functions produced by shocks are quite different in the critical energy range for producing the observed synchrotron emission. We also showed that transrelativistic effects produced proton spectra that were not in agreement with standard results from radio observations, but that the lepton spectra were, in fact, consistent with observations. We performed simulations of relativistic shocks (shocks where the flow speed is a sizable fraction of the speed of light) and discovered some interesting effects. We first demonstrated the power of the Monte Carlo technique by determining the shock jump conditions in relativistic shocks. We then proceeded to determine how relativistic shocks accelerate particles. We found that nonlinear relativistic shocks treat protons and leptons even more differently than nonrelativistic shocks. The transrelativistic effects on the shock structure from the heavy ion component reduces the lepton acceleration to a tiny fraction of the ion acceleration. This effect is dramatic even if high energy leptons (many times thermal energy) are injected, and was totally unexpected. Our results have important consequences for astrophysical environments expected to harbor relativistic flows such as extra-galactic radio sources and accretion onto compact objects.

Ellison, Donald C.↗

NASA NTRS1993

Gas-grain energy transfer in solar nebula shock waves: Implications for the origin of chondrules

Meteoritic chondrules provide evidence for the occurrence of rapid transient heating events in the protoplanetary nebula. Astronomical evidence suggests that gas dynamic shock waves are likely to be excited in protostellar accretion disks by processes such as protosolar mass ejections, nonaxisymmetric structures in an evolving disk, and impact on the nebula surface of infalling 'clumps' of circumstellar gas. Previous detailed calculations of gas-grain energy and momentum transfer have supported the possibility that such shock waves could have melted pre-existing chondrule-sized grains. The main requirement for grains to reach melting temperatures in shock waves with plausibly low Mach numbers is that grains existed in dust-rich zones (optical depth greater than 1) where radiative cooling of a given grain can be nearly balanced by radiation from surrounding grains. Localized dust-rich zones also provide a means of explaining the apparent small spatial scale of heating events. For example, the scale size of at least some optically thick dust-rich zones must have been relatively small (less than 10 kilometers) to be consistent with petrologic evidence for accretion of hot material onto cold chondrules. The implied number density of mm-sized grains for these zones would be greater than 30 m(exp -3). In this paper, we make several improvements of our earlier calculations to include radiation self-consistently in the shock jump conditions, and we include heating of grains due to radiation from the shocked gas. In addition, we estimate the importance of momentum feedback of dust concentrations onto the shocked gas which would tend to reduce the efficiency of gas dynamic heating of grains in the center of the dust cloud.

Hood, L. L.↗

NASA NTRS1995

Entropy jump across an inviscid shock wave

The shock jump conditions for the Euler equations in their primitive form are derived by using generalized functions. The shock profiles for specific volume, speed, and pressure are shown to be the same, however density has a different shock profile. Careful study of the equations that govern the entropy shows that the inviscid entropy profile has a local maximum within the shock layer. We demonstrate that because of this phenomenon, the entropy, propagation equation cannot be used as a conservation law.

Salas, Manuel D.↗

NASA NTRS1994

Motion of the heliospheric termination shock. 2: Energy loss due to energetic particle acceleration

We present a generalization of earlier analysis of the motion of the heliospheric termination shock in response to heliospheric disturbances (Barnes, 1993) (paper 1), to allow jump conditions that include an energy sink at the shock front. The motivation for this study is that acceleration of the anomalous cosmic ray component may in fact represent such a sink. We have idealized the situation by assuming an infinitely thin shock parameterized by a quantity lambda(0 less than or equals lambda less than or equals 1), defined as the fraction of solar wind energy that is lost due to acceleration of the energetic particle component. If the sink is strong (consuming, say, 50% or more of the incident solar wind energy), the model leads to the following principal conclusions: (1) the shocked plasma would be much denser and cooler than in the standard gasdynamic case, thereby leading to more favorable conditions for direct observation of the shocked plasma; (2) the equilibrium shock position would be slightly farther (less than 10%) from the Sun than in the standard model; (3) as in the gasdynamic case, the shock would normally be in motion, so that the instantaneous position of the termination shock is not determined by interstellar conditions but by the recent history of the wind that has passed through the termination shock; and (4) the response of the shock to upstream disturbances would be similar to the response in the gasdynamic case, but the speed of the new termination shock would be somewhat smaller (probably by a f actor of 4 or less). We estimate that this speed is typically, approximately a few tens of kilometers per second corresponding to an inward or outward excursion of order of less than 1 to several AU, rather less than conventional estimates of several tens of AU.

Barnes, Aaron↗

NASA NTRS1994

Relativistic particle acceleration in plerions

We discuss recent research on the structure and particle acceleration properties of relativistic shock waves in which the magnetic field is transverse to the flow direction in the upstream medium, and whose composition is either pure electrons and positrons or primarily electrons and positrons with an admixture of heavy ions. Particle-in-cell simulation techniques as well as analytic theory have been used to show that such shocks in pure pair plasmas are fully thermalized -- the downstream particle spectra are relativistic Maxwellians at the temperature expected from the jump conditions. On the other hand, shocks containing heavy ions which are a minority constituent by number but which carry most of the energy density in the upstream medium do put approximately 20% of the flow energy into a nonthermal population of pairs downstream, whose distribution in energy space is N(E) varies as E(exp -2), where N(E)dE is the number of particles with energy between E and E+dE. The mechanism of thermalization and particle acceleration is found to be synchrotron maser activity in the shock front, stimulated by the quasi-coherent gyration of the whole particle population as the plasma flowing into the shock reflects from the magnetic field in the shock front. The synchrotron maser modes radiated by the heavy ions are absorbed by the pairs at their (relativistic) cyclotron frequencies, allowing the maximum energy achievable by the pairs to be gamma(sub +/-)m(sub +/-)c squared = m(sub i)c squared gamma(sub 1)/Z(sub i), where gamma(sub 1) is the Lorentz factor of the upstream flow and Z(sub i) is the atomic number of the ions. The shock's spatial structure is shown to contain a series of 'overshoots' in the magnetic field, regions where the gyrating heavy ions compress the magnetic field to levels in excess of the eventual downstream value. This shock model is applied to an interpretation of the structure of the inner regions of the Crab Nebula, in particular to the 'wisps,' surface brightness enhancements near the pulsar. We argue that these surface brightness enhancements are the regions of magnetic overshoot, which appear brighter because the small Larmor radius pairs are compressed and radiate more efficiently in the regions of more intense magnetic field. This interpretation suggests that the structure of the shock terminating the pulsar's wind in the Crab Nebula is spatially resolved, and allows one to measure gamma(sub 1) and a number of other properties of the pulsar's wind. We also discuss applications of the shock theory to the termination shocks of the winds from rotation-powered pulsars embedded in compact binaries. We show that this model adequately accounts for (and indeed predicted) the recently discovered X-ray flux from PSR 1957+20, and we discuss several other applications to other examples of these systems.

Arons, Jonathan↗

NASA NTRS1995

Dissipative MHD solutions for resonant Alfven waves in 1-dimensional magnetic flux tubes

The present paper extends the analysis by Sakurai, Goossens, and Hollweg (1991) on resonant Alfven waves in nonuniform magnetic flux tubes. It proves that the fundamental conservation law for resonant Alfven waves found in ideal MHD by Sakurai, Goossens, and Hollweg remains valid in dissipative MHD. This guarantees that the jump conditions of Sakurai, Goossens, and Hollweg, that connect the ideal MHD solutions for xi(sub r), and P' across the dissipative layer, are correct. In addition, the present paper replaces the complicated dissipative MHD solutions obtained by Sakurai, Goossens, and Hollweg for xi(sub r), and P' in terms of double integrals of Hankel functions of complex argument of order 1/3 with compact analytical solutions that allow a straight- forward mathematical and physical interpretation. Finally, it presents an analytical dissipative MHD solution for the component of the Lagrangian displacement in the magnetic surfaces perpen- dicular to the magnetic field lines xi(sub perpendicular) which enables us to determine the dominant dynamics of resonant Alfven waves in dissipative MHD.

Goossens, Marcel↗

NASA NTRS1995

HPTAM, a two-dimensional Heat Pipe Transient Analysis Model, including the startup from a frozen state

A two-dimensional Heat Pipe Transient Analysis Model, 'HPTAM,' was developed to simulate the transient operation of fully-thawed heat pipes and the startup of heat pipes from a frozen state. The model incorporates: (a) sublimation and resolidification of working fluid; (b) melting and freezing of the working fluid in the porous wick; (c) evaporation of thawed working fluid and condensation as a thin liquid film on a frozen substrate; (d) free-molecule, transition, and continuum vapor flow regimes, using the Dusty Gas Model; (e) liquid flow and heat transfer in the porous wick; and (f) thermal and hydrodynamic couplings of phases at their respective interfaces. HPTAM predicts the radius of curvature of the liquid meniscus at the liquid-vapor interface and the radial location of the working fluid level (liquid or solid) in the wick. It also includes the transverse momentum jump condition (capillary relationship of Pascal) at the liquid-vapor interface and geometrically relates the radius of curvature of the liquid meniscus to the volume fraction of vapor in the wick. The present model predicts the capillary limit and partial liquid recess (dryout) in the evaporator wick, and incorporates a liquid pooling submodel, which simulates accumulation of the excess liquid in the vapor core at the condenser end.

Tournier, Jean-Michel↗

NASA NTRS1998

The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws: The Two Dimensional Time Marching Schemes - 1

A new high resolution and genuinely multidimensional numerical method for solving conservation laws is being, developed. It was designed to avoid the limitations of the traditional methods. and was built from round zero with extensive physics considerations. Nevertheless, its foundation is mathmatically simple enough that one can build from it a coherent, robust. efficient and accurate numerical framework. Two basic beliefs that set the new method apart from the established methods are at the core of its development. The first belief is that, in order to capture physics more efficiently and realistically, the modeling, focus should be placed on the original integral form of the physical conservation laws, rather than the differential form. The latter form follows from the integral form under the additional assumption that the physical solution is smooth, an assumption that is difficult to realize numerically in a region of rapid chance. such as a boundary layer or a shock. The second belief is that, with proper modeling of the integral and differential forms themselves, the resulting, numerical solution should automatically be consistent with the properties derived front the integral and differential forms, e.g., the jump conditions across a shock and the properties of characteristics. Therefore a much simpler and more robust method can be developed by not using the above derived properties explicitly.

Chang, Sin-Chung↗

NASA NTRS1999

Condensation of Forced Convection Two-Phase Flow in a Miniature Tube

A physical/mathematical model of annular film condensation at the inlet of a miniature tube has been developed. In the model, the liquid flow is coupled with the vapor flow along the liquid-vapor interface through the interfacial temperature, heat flux, shear stress, and pressure jump conditions due to surface tension effects. The model predicts the shape of the liquid-vapor interface along the condenser and leads to the conclusion that there is complete condensation at a certain distance from the condenser inlet. The numerical results show that complete condensation of the incoming vapor is possible at comparatively low heat loads and that this is a special case of a more general condensation regime with two-phase bubbly flow downstream of the initial annular film condensation region. Observations from the flow visualization experiment confirm the existence and qualitative features of annular film condensation leading to the complete condensation phenomenon in a small diameter (3.25 mm) circular tube condenser.

Begg, E.↗

NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗

NASA NTRS2001

Nature of Fluctuations on Directional Discontinuities Inside a Solar Ejection: Wind and IMP 8 Observations

A solar ejection passed the Wind spacecraft between December 23 and 26, 1996. On closer examination, we find a sequence of ejecta material, as identified by abnormally low proton temperatures, separated by plasmas with typical solar wind temperatures at 1 AU. Large and abrupt changes in field and plasma properties occurred near the separation boundaries of these regions. At the one boundary we examine here, a series of directional discontinuities was observed. We argue that Alfvenic fluctuations in the immediate vicinity of these discontinuities distort minimum variance normals, introducing uncertainty into the identification of the discontinuities as either rotational or tangential. Carrying out a series of tests on plasma and field data including minimum variance, velocity and magnetic field correlations, and jump conditions, we conclude that the discontinuities are tangential. Furthermore, we find waves superposed on these tangential discontinuities (TDs). The presence of discontinuities allows the existence of both surface waves and ducted body waves. Both probably form in the solar atmosphere where many transverse nonuniformities exist and where theoretically they have been expected. We add to prior speculation that waves on discontinuities may in fact be a common occurrence. In the solar wind, these waves can attain large amplitudes and low frequencies. We argue that such waves can generate dynamical changes at TDs through advection or forced reconnection. The dynamics might so extensively alter the internal structure that the discontinuity would no longer be identified as tangential. Such processes could help explain why the occurrence frequency of TDs observed throughout the solar wind falls off with increasing heliocentric distance. The presence of waves may also alter the nature of the interactions of TDs with the Earth's bow shock in so-called hot flow anomalies.

Vasquez, Bernard J.↗

NASA NTRS2007

A Physical Relationship Between Electron-Proton Temperature Equilibration and Mach Number in Fast Collisionless Shocks

The analysis of Balmer-dominated optical spectra from nonradiative (adiabatic) SNRs has shown that the ratio of the electron to proton temperature at the blast wave is close to unity at v(s) approximate or smaller than 400 km s(-1) but declines sharply down to the minimum value of m(e) /m(p) dictated by the jump conditions at shock speeds exceeding 2000 km s(-1). We propose a physical model for the heating of electrons and ions in non-cosmic-ray-dominated, strong shocks (v(s) 400 km s(-1)) wherein the electrons are heated by lower hybrid waves immediately ahead of the shock front. These waves arise naturally from the cosmic ray pressure gradient upstream from the shock. Our model predicts a nearly constant level of electron heating over a wide range of shock speeds, producing a relationship (T(e)/T(p))sub 0 proportional to v(-2/s) (proportional to M(-2)) that is fully consistent with the observations.

SNR(ONRADIATIVE SUPERNOVA REMNANTS)↗

NASA NTRS2015

Method and System for Weakening Shock Wave Strength at Leading Edge Surfaces of Vehicle in Supersonic Atmospheric Flight

A method and system are provided to weaken shock wave strength at leading edge surfaces of a vehicle in atmospheric flight. One or more flight-related attribute sensed along a vehicle's outer mold line are used to control the injection of a non-heated, non-plasma-producing gas into a local external flowfield of the vehicle from at least one leading-edge surface location along the vehicle's outer mold line. Pressure and/or mass flow rate of the gas so-injected is adjusted in order to cause a Rankine-Hugoniot Jump Condition along the vehicle's outer mold line to be violated.

Daso, Endwell O.↗

NASA NTRS2016

A Strong Merger Shock in Abell 665

Deep (103 ks) Chandra observations of Abell 665 have revealed rich structures in this merging galaxy cluster, including a strong shock and two cold fronts. The newly discovered shock has a Mach number of M =3.0 +/- 0.6, propagating in front of a cold disrupted cloud. This makes Abell 665 the second cluster, after the Bullet cluster, where a strong merger shock of M is approximately 3 has been detected. The shock velocity from jump conditions is consistent with (2.7 +/- 0.7) × 10(exp 3) km s(exp −1). The new data also reveal a prominent southern cold front with potentially heated gas ahead of it. Abell 665 also hosts a giant radio halo. There is a hint of diffuse radio emission extending to the shock at the north, which needs to be examined with better radio data. This new strong shock provides a great opportunity to study the reacceleration model with the X-ray and radio data combined.

galaxies: clusters: general- galaxies: clisters; i↗

NASA NTRS2000

A Computational Model for Rotor-Fuselage Interactional Aerodynamics

A novel unsteady rotor-fuselage interactional aerodynamics model has been developed. This model loosely couples a Generalized Dynamic Wake Theory (GDWT) to a thin-layer Navier-Stokes solution procedure. This coupling is achieved using an unsteady pressure jump boundary condition in the Navier-Stokes model. The new unsteady pressure jump boundary condition models each rotor blade as a moving pressure jump which travels around the rotor azimuth and is applied between two adjacent planes in a cylindrical, non-rotating grid. Comparisons are made between measured and predicted time-averaged and time-accurate rotor inflow ratios. Additional comparisons are made between measured and predicted unsteady surface pressures on the top centerline and sides of the fuselage.

D. Douglas Boyd, Jr.↗

NASA NTRS1993

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗

NASA NTRSDate not supplied

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

Cameron J. Hines↗

NASA NTRSDate not supplied

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗