Results for “material flows” — Page 86

NASA NTRS1996

Directional Solidification of Mercury Cadmium Telluride During the Second United States Microgravity Payload Mission (USMP-2)

As a solid solution semiconductor having, a large separation between liquidus and solidus, mercury cadmium telluride (MCT) presents a formidable challenge to crystal growers desiring an alloy of high compositional uniformity. To avoid constitutional supercooling during Bridgman crystal growth it is necessary to solidify slowly in a high temperature gradient region. The necessary translation rate of less than 1 mm/hr results in a situation where fluid flow induced by gravity on earth is a significant factor in material transport. The Advanced Automated Directional Solidification Furnace (AADSF) is equipped to provide the stable thermal environment with a high gradient, and the required slow translation rate needed. Ground based experiments in AADSF show clearly the dominance of flow driven transport. The first flight of AADSF in low gravity on USMP-2 provided an opportunity to test theories of fluid flow in MCT and showed several solidification regimes which are very different from those observed on earth. Residual acceleration vectors in the orbiter during the mission were measured by the Orbital Acceleration Research Experiment (OARE), and correlated well with observed compositional differences in the samples.

Gillies, D. C.↗

NASA NTRS1995

Geology and landscape evolution of the Hellas region of Mars

Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S 60 deg E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (>1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topography, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approx. 1- to 2-km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approx. 100-200 m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approx. 100 km across; however, we disagree with a hypothesis that suggests that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Tanaka, Kenneth L.↗

DOE OSTI2023

Flow Induced Vibration Studies in Pressurized Helium Gas Cooling Channels

Production of metastable Technetium-99 (Tc-99m), a radioactive tracer that emits gamma rays, is vital to the medical imaging community. Tc-99m is extracted from the decay of Molybdenum-99 (Mo-99) which has a half-life of about 2-3 days. The work presented in this report is part of the NNSA’s mission to produce Mo-99 commercially, within the US, without the use of highly enriched uranium (HEU) in support of nonproliferation and global security. Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 through the irradiation of Mo-100 targets using an electron beam. The NMR target comprises a stack of approximately 60-70 Mo-100 disks with diameter 24 mm and thickness 0.74 mm held in stainless steel laminations, each separated using 0.25 mm thick stainless steel spacers. The symmetric target stack is housed in an Inconel vessel with two Inconel windows on either side. Two electron accelerators are used to produce 40 MeV, 3.16 µA electron beams each that penetrate the Inconel windows and irradiate the Mo-100 disks. Approximately 90% of the total 250 kW beam power is deposited in the NMR target during the irradiation process. During irradiation, pressurized helium gas flows through 0.25 mm thin gaps between the disks cooling the beam window, target disks, disk laminations and spacers. Both NMR and LANL have found during cold testing of the target system (no heat deposition) that the Mo-100 disks undergo significant mass loss and disk breakage due to vibrations induced by the flowing helium gas. The mass loss is not only undesirable due to monetary loss from reduced final quantities of Mo-99, but also due to the hazards associated with radioactive material trapped in the cooling lines and particle filters. The effect of flow rate and target geometry on the flow induced vibrations need to be quantified, and recommendations provided to minimize this mass loss. This work describes LANL’s experimental characterization of the flow induced vibrations and disk mass loss in a reduced scale set-up containing 10 Mo-100 disks. We use high speed imaging, displacement measurements and microphone measurements combined with signal processing to estimate the vibration frequency of each disk. The effect of disk thickness, target fit and duration of testing on the mass loss is described. We find that in the current configuration of NMR targets, the vibrations and mass loss on the first disk are minimized, while those in the adjacent disks are highest. The microphone and high-speed image data show that increased flow rates and increased duration of testing increases vibration frequency and mass loss. The mass loss is due to both disk rotation and back and forth motion. There are visible wear marks on the disks with the highest mass loss. We also note that the current NMR window gap reduces flow induced vibrations compared to the previous smaller gaps. Longer duration testing will provide more data and verification for the findings presented in this report. The work will be continued in FY 24.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Window Cooling Studies and Disk Vibration Testing on a Subset of Mo-100 Disks

Production of metastable Technetium-99 (Tc-99m), a radioactive tracer that emits gamma rays, is vital to the medical imaging community. Tc-99m is extracted from the decay of Molybdenum-99 (Mo-99) which has a half-life of about 2-3 days. The work presented in this report is part of the NNSA’s mission to produce Mo-99 commercially, within the US, without the use of highly enriched uranium (HEU) in support of nonproliferation and global security. Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 through the irradiation of Mo-100 targets using an electron beam. The NMR target comprises a stack of approximately 60-70 Mo-100 disks with diameter 24 mm and thickness 0.74 mm held in stainless steel laminations, each separated using 0.25 mm thick stainless-steel spacers. The symmetric target stack is housed in an Inconel vessel with two Inconel windows on either side. Two electron accelerators are used to produce 40 MeV, 3.16 µA electron beams each that penetrate the Inconel windows and irradiate the Mo-100 disks. Approximately 90% of the total 250 kW beam power is deposited in the NMR target during the irradiation process, with a smaller percentage adding up to 2.2 kW of heat deposited on the Inconel window. During irradiation, pressurized helium gas flows through thin gaps between the disks cooling the beam window, target disks, disk laminations and spacers. Both NMR and LANL have found during cold testing of the target system (no heat deposition) that the Mo 100 disks undergo significant mass loss and disk breakage due to vibrations induced by the flowing helium gas. The mass loss is not only undesirable due to monetary loss from reduced final quantities of Mo-99, but also due to the hazards associated with radioactive material trapped in the cooling lines and particle filters. The effect of flow rate and target geometry on the flow induced vibrations need to be quantified, and recommendations provided to minimize this mass loss. LANL has previously also tested NMR’s Inconel beam window by heating the window, while flowing pressurized helium, using the average heat deposited on the window. However, the NMR beam is pulsed with a duty cycle of 12.5%, which introduces oscillation in temperature around the nominal 600 °C steady state value with each pulse. Available fatigue curves for Inconel are few, established for room temperature, and they are based on mechanical strain cycles not thermally induced strain as in the NMR target. The effect of pulsed beam heating on the Inconel window therefore needs to be quantified. This report details the experiments conducted to assess the factors that lead to mass loss in the NMR target disks as well as to understand the effect of a pulsed beam on NMR’s Inconel window. This work describes LANL’s experimental characterization of the flow induced vibrations and disk mass loss in a reduced scale set-up containing 5 to 10 Mo-100 disks. We use high speed imaging, displacement measurements and microphone measurements combined with signal processing to estimate the vibration frequency of each disk. The effect of disk thickness, target fit and duration of testing on the mass loss is described. We find that in the current configuration of NMR targets, the vibrations and mass loss on the first disk are minimized, while those in the adjacent disks are highest. The microphone and high-speed image data show that increased flow rates and increased duration of testing increases vibration frequency and mass loss. The mass loss is due to both disk rotation and back and forth motion. There are visible wear marks on the disks with the highest mass loss. We also note that the current NMR window gap reduces flow induced vibrations compared to the previous smaller gaps. Improved target holders significantly reduce disk mass loss to almost negligible quantities. This work finds that the larger window to first disk gap and improved target holder geometry should allow NMR to successfully conduct irradiations with minimal mass loss. The window tests were conducted to understand the effect of a pulsed beam on both the window longevity and to estimate the window temperature and displacement during pulsing. The experiments presented here were performed at significantly low power, due to the limitations of the induction heating system. The window temperature rose to approximately 73 °C with a significantly reduced power of 45 W without beam pulsing. With a 5 Hz pulse rate, 12.5% duty cycle, the window temperature remained constant at 26 °C. These experiments will be repeated with improved coil geometry and reported in upcoming journal papers.

42 ENGINEERING↗

NASA NTRS1994

Shear stabilization of solidification fronts

A linear-stability analysis is performed on the interface formed during the directional solidification of a dilute binary alloy in the presence of a time-periodic flow. In general, the flow is generated by the elliptical motion of the crystal parallel to the interface. The presence of the flow can either stabilize or destabilize the system relative to the case without flow with the result depending on the frequency and amplitude of the oscillations as well as the properties of the material. In the particular, we find that the morphological instability present in the absence of flow can be entirely suppressed with respect to disturbances of the same frequency as the oscillation.

Davis, Stephen H.↗

NASA NTRS1999

Near Surface Vapor Bubble Layers in Buoyant Low Stretch Burning of Polymethylmethacrylate

Large-scale buoyant low stretch stagnation point diffusion flames over solid fuel (polymethylmethacrylate) were studied for a range of aerodynamic stretch rates of 2-12/ sec which are of the same order as spacecraft ventilation-induced stretch in a microgravity environment. An extensive layer of polymer material above the glass transition temperature is observed. Unique phenomena associated with this extensive glass layer included substantial swelling of the burning surface, in-depth bubble formation, and migration and/or elongation of the bubbles normal to the hot surface. The bubble layer acted to insulate the polymer surface by reducing the effective conductivity of the solid. The reduced in-depth conduction stabilized the flame for longer than expected from theory neglecting the bubble layer. While buoyancy acts to move the bubbles deeper into the molten polymer, thermocapillary forces and surface regression both act to bring the bubbles to the burning surface. Bubble layers may thus be very important in low gravity (low stretch) burning of materials. As bubbles reached the burning surface, monomer fuel vapors jetted from the surface, enhancing burning by entraining ambient air flow. Popping of these bubbles at the surface can expel burning droplets of the molten material, which may increase the fire propagation hazards at low stretch rates.

Olson, Sandra L.↗

DOE OSTI2024

Flow and Heat Transfer Experimental Study for 3D-Printed Solar Receiving Tubes With Helical Fins at Internal Surface

3D-printing technology was applied to fabricate novel solar thermal collection tubes that have internal heat transfer enhancement fins and external surfaces with high solar absorptivity and low emissivity due to the ability to use different materials in one tube. Helical fins were selected to introduce circumferential flow and thus minimize the circumferential temperature difference of the tube that receives sunlight on one side. The structures of the helical fins were previously optimized from computational fluid dynamics (CFD) analysis with the objective of low entropy production rate by looking for high heat transfer coefficient and relatively lower pressure loss. High-temperature alloy, Inconel-718, was used to 3D print the tubes, which can resist corrosion for the potential application of molten chloride salts as heat transfer fluid. Experimental tests were carried out using water as the heat transfer fluid with the high heat flux provided by a tubular furnace heater. The tested Reynolds number ranges from 3.9 × 10 3 to 6.1 × 10 4 . Heat transfer coefficients of up to 2.8 times that of the smooth tube could be obtained with the expense of increased pressure loss compared to that of the smooth tube. The total system entropy generation can be significantly reduced due to the benefit of heat transfer enhancement that is greater than the expenses of the increased pressure loss. The experimental results of the 3D-printed heat transfer tubes confirmed the CFD-based results of fin optimization. Furthermore, the novel heat transfer tube is recommended for application in concentrating solar power systems.

14 SOLAR ENERGY↗

NASA NTRS1976

Accretion by the Galaxy - Effects of radiative cooling on the flow structure and infall rate

The extreme effects of radiative cooling on the structure of an inflow of material into the Galaxy are explored. Significantly, it is found that cooling cannot prevent the formation of a very hot (about 2 million K) corona in the flow close to the Galaxy. When it is important, however, cooling in the flow considerably increases the infall rate over the adiabatic maximum value for the same set of boundary conditions. Limits set on the flow rate and coronal gas parameters by observations of the Galaxy are discussed.

Cox, D. P.↗

DOE OSTI2025

MPACT Safeguards Modeling: FY25 Update

Sandia National Laboratories develops and maintains several open-source software packages to support material accountancy analyses. This includes the Material Accountancy Performance Indicator Toolkit (MAPIT), the Fissile Facility Flow Modeler (F3M) and the Separation and Safeguards Performance Model Library (SSPM-L). MAPIT is responsible for performing statistical safeguards analyses on bulk and itemized data from nuclear fuel cycle facilities and can operate on real or synthetic data. MAPIT is the only open-source software for such analyses. F3M is a library of modules, built in MATLAB Simulink, that contain pre made blocks to represent different generic fuel cycle processes. These blocks can be used together in a modular fashion to represent and simulate nuclear fuel cycle processes with the goal of improving facility-level accountancy during the design phase. F3M is also an open-source library. Finally, the SSPM-L library is a series of completed models built from F3M. The library includes facility models such as a generic PUREX facility and a fuel fabrication facility. The SSPM-L library is not open source, but is available to collaborators with a relevant use case. These tools include modeling and simulation pipelines to simulate nuclear fuel cycle facilities and the underlying software needed to simulate measurement uncertainty and perform statistical analyses. Together, these tools can perform end-to-end nuclear material accountancy analyses. This report documents the various improvements made to these tools in FY25. Specifically, we added new statistical test, new statistical modeling capabilities, new fuel cycle facility models, and launched a new open-source model component library.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

The Impact of Oxidation-Induced Degradation On Materials Used in Hydrogen-Fired Microturbines

Hydrogen-fueled microturbines are being considered as part of the future green microgrid. However, the use of hydrogen as a fuel presents new challenges for selection and development of suitable high temperature materials for hydrogen combustion. The burning of hydrogen is expected to result in higher operating temperatures and higher than typically observed water vapor contents in exhaust gases versus burning natural gas. In the present work, foil specimens of various Fe- and Ni-based alloys were oxidized in air + 10 % H 2 O and air + 60% H 2 O for up to 5,000 h at 700 °C to simulate the exhaust atmosphere of natural gas and hydrogen-fueled microturbines. Here, the impact of alloy composition and water vapor content on the oxidation/ volatilization induced loss of wall thickness was experimentally evaluated. Enhanced external oxidation and volatilization of Cr 2 O 3 and Ti-doped Cr 2 O 3 scales was observed in air + 60% H 2 O compared to air + 10% H 2 O. No significant impact of the higher water vapor content was observed on Al 2 O 3 scales formed on Fe-based alumina forming alloys. Lifetime modeling was employed to predict the combined effects of water vapor content, gas flow rates, temperature and alloy composition on the oxidation-induced lifetime of the investigated materials.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Mechanical behavior and texture evolution of WE43 magnesium-rare earth alloy in Split-Hopkinson Pressure Bar and Taylor Impact Cylinder Testing

Mechanical behavior and texture evolution of Mg rare-earth alloy WE43 is investigated for strain-rates 10 -3 /s to upwards of 10 5 /s for the two material conditions - as-cast (AC) and T6 age hardened, rolled plate (RT6). The high strain-rate behavior is tested using both Taylor cylinder impact tests (TC) and split Hopkinson pressure bar tests (SHB) and bulk textures are obtained using neutron diffraction. Unlike the quasi-static strained material, AC and RT6 SHB retained high hardening rates throughout the test, even up to 30% true strain. Moreover, the high strain-rate data revealed that the RT6 material has a much higher strength than the AC material, but similar hardening rates despite significantly different initial texture. The flow stress near yield increased up to I 0% for RT6 and up to 30% for AC as the strain-rate increased six orders of magnitude from quasi-static rates 10 -3 /s to 10 3 /s. Neither material exhibited significant plastic anisotropy over the broad range of strain rates, despite the fact that the RT6 material had a moderately strong initial texture. In the TC tests, the geometric cross-sectional changes and texture along the cylinder from the cylindrical sample foot to head are measured and from the neutron diffraction texture analysis, upper-bound estimates of twin volume fraction are obtained as well as dislocation density from analyzing diffraction peak broadening. Recorded geometrical changes along several loading directions show that the material has deformed remarkably homogeneously under impact. Analysis of deformed textures indicates that {101¯2} extension deformation twinning occurred in the RT6 condition over the range of strain rates, with an upper bound estimate of 40% twin volume fraction for approximately 0.10-0.25 true strain. The peak texture components after the impact have their c-axes closely aligned with the impact direction. As a result, these observations are presented and rationalized in the paper.

42 ENGINEERING↗

NASA NTRS1994

Thermonuclear runaways in nova outbursts. 2: Effect of strong, instantaneous, local fluctuations

In an attempt to understand the manner in which nova outbursts are initiated on the surface of a white dwarf, we investigate the effects fluctuations have on the evolution of a thermonuclear runaway. Fluctuations in temperature density, or the composition of material in the burning shell may arise due to the chaotic flow field generated by convection when it occurs, or by the accretion process itself. With the aid of two-dimensional reactive flow calculations, we consider cases where a strong fluctutation in temperature arises during the early, quiescent accretion phase or during the later, more dynamic, explosion phase. In all cases we find that an instantaneous, local temperature fluctuation causes the affected material to become Rayleigh-Taylor unstable. The rapid rise and subsequent expansion of matter immediately cools the hot blob, which prevents the lateral propagation of burning. This suggests that local temperature fluctuations do not play a significant role in directly initiating the runaway, especially during the early stages. However, they may provide an efficient mechanism of mixing core material into the envelope (thereby pre-enriching the fuel for subsequent episodes of explosive hydrogen burning) and of mixing substantial amounts of the radioactive nucleus N-13 into the surface layers, making novae potential gamma-ray sources. This suggests that it is the global not the local, evolution of the core-envelope interface to high temperatures which dominates the development of the runaway. We also present a possible new scenario for the initiation of nova outbursts based on our results.

Shankar, Anurag↗

DOE OSTI2024

Hydrogen embrittlement of Zircaloy-4 fabricated by ultrasonic additive manufacturing

Ultrasonic additive manufacturing (UAM) was successfully applied to the zirconium material system to create a planar geometry. Following fabrication, SS-J3 type tensile specimens of the UAM and wrought Zircaloy-4 with a nominal gage section of 5×1.2×0.75 mm were machined for hydriding studies and mechanical testing. The SS-J3 type tensile specimens were gas-charged with hydrogen using a custom system that precisely controls hydrogen gas flow rate, hydrogen partial pressure, and temperature. To avoid altering the UAM material, the maximum process temperature was limited to 550°C. Using different initial hydrogen gas pressures and flow rates, various hydrogen contents (70–1755 wppm) were achieved for Zircaloy-4 specimens. Tensile testing shows that, regardless of hydrogen content, all UAM specimens measured yield strengths in the range of 557±16 MPa and ultimate tensile strengths of 660±4MPa. However, total elongation clearly decreased as a function of increasing hydrogen content. Further, at the lowest hydrogen content, the total plastic elongation measured 21.5%, and this value decreased to less than 2% when the hydrogen content was increased to 1000 wppm. Cross-sectional optical microscopy images revealed that the hydride distributions are randomly oriented. For the low hydrogen content specimen, these randomly distributed hydrides are isolated from each other. A sandwiched structure, consisting of high-density and low-density layers, was also observed for the specimens with hydrogen content between 200 and 400 wppm. As the hydrogen content increases, the hydrides diffuse into the low-density hydride layers to form a network across the whole specimen. Although the tensile properties of UAM and wrought Zircaloy-4 exhibit the same behavior with increasing hydrogen content, the distinctions in grain orientations led to differences in the hydride orientations at low and intermediate hydrogen concentrations.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Incorporating dynamic recrystallization into a crystal plasticity model for high-temperature deformation of Ti-6Al-4V

During hot deformation of (α + β) titanium alloys, the simultaneous action of strain and temperature in the (α + β) regime facilitates dynamic recovery, dynamic recrystallization (DRX), and phase transformations via non-equilibrium paths. DRX is manifested in the form of fine recrystallized α or β grains. Here, in the present study, a two-phase crystal plasticity finite element framework (CP-DRX) has been developed which incorporates DRX kinetics into the crystal plasticity (CP) model to predict the flow characteristics of Ti-6Al-4V alloys during thermo-mechanical processing. The CP slip system parameters, as well as elastic properties from both α and β phases of Ti, are calibrated for different strain rate conditions. An EBSD-informed two-phase microstructure representation has been utilized in the CP-DRX framework to explore the time-dependent evolution of DRX microstructure and crystal orientation for different strain rate conditions. The proposed CP-DRX can capture the evolution of crystal orientation and plastic flow stress–strain response of polycrystalline Ti-6Al-4V during the deformation process. Furthermore, the proposed model is able to capture the softening behavior, observed in average stress–strain response from experiments performed using a Gleeble thermomechanical simulator and predict the recrystallization texture.

(α+ β)-Titanium (Ti) Alloys↗

NASA NTRS1984

Surface acoustic admittance and absorption of highly porous, layered, fibrous materials

Some acoustic properties of Kevlar-29 - a fine fibered, layered material is investigated. Kevlar is characterized by very high strength, uniform filaments arranged in a parallel batt where most filaments are random in the x-y plane but ordered as planes in the z direction. For experimental purposes, volume porosity, static flow resistance and mean filament diameter are used to identify the material. To determine the acoustic surface admittance of Kevlar, batts of the material are cut into small pads and placed into a standing wave tube terminated by a rigid brass plug. The attenuation and relative phase shift are recorded at each frequency in the range of 50 to 6000 Hz. Normalized conductance and susceptance are combined to form the acoustic absorption coefficient. The data are compared with theory by plotting the normalized admittance and normal incident absorption coefficient versus cyclic frequency.

Tesar, J. S.↗

NASA NTRS1973

Cyclic arc plasma tests of RSI materials using a preheater

The results of a test program are reported in which a preheater was used with an arc plasma stream to study the thermal response of samples of candidate reusable surface insulation materials for the space shuttle. The preheater simulated the shuttle temperature history during the first and last portions of the test cycle, which could not be simulated by the air arc plasma flow. Pre- and post-test data taken for each of the materials included magnified views, optical properties, and chemical analyses. The test results indicate that the mullite base samples experience higher surface temperatures than the other materials at heating rates greater than 225 kw/sq m. The ceramic fibrous mullite and silica coatings show noncatalytic wall behavior. Internal temperature response data for the materials are compared and correlated with analytical predictions.

Stewart, D. A.↗

NASA NTRSDate not supplied

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗

NASA NTRSDate not supplied

Simulations of Cryogenic Line Chilldown with Advanced Sub-Grid Wall Boiling Models

A meso-scale model developed at MIT [2][11] for boiling processes in water was adapted for cryogens and demonstrated for chilldown of propellant transfer lines in both liquid nitrogen and hydrogen. The sub-grid boiling model accurately captures the contributions to heat transfer from the generation of bubble nuclei, growth, and interaction of the bubbles in the microlayer as well as quenching of the boiling surface following bubble departure. It was adapted for cryogenic fluids using thermodynamic scaling concepts taking into account non-dimensional pressures and temperatures that are scaled by the corresponding critical values for the fluid. The boiling model was demonstrated for line chilldown in liquid nitrogen (Darr et al. [6]). The predicted wall temperature at which quenching occurs was close to the test data while the slope of the temperature curve after quenching is initiated, showing a steeper variation than the test data. Simulations were also performed for liquid hydrogen by simulating experiments of Hartwig et al. [7]. Chilldown times in liquid hydrogen are much more rapid due to higher heat transfer in the film boiling regime and accounting for the higher turbulence levels were found to be important. Furthermore, at the much lower fluid temperatures in liquid hydrogen flows, accounting for the variable thermal properties of the solid material is critical and has a dramatic impact on the quench times. The need for additional studies to better understand the evolution of the quench front in liquid hydrogen was noted.

Computational Fluid Dynamics↗