SEARCH · Engineering Papers

Results for “BOILING”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Liquid Nitrogen Zero Boiloff Testing

Cryogenic propellants such as liquid hydrogen (LH2) and liquid oxygen (LO2) are a part of NASAs future space exploration due to their high specific impulse for rocket motors of upper stages suitable for transporting 10s to 100s of metric tons of payload mass to destinations outside of low earth orbit and for their return. However, the low storage temperatures of LH2 and LO2 cause substantial boil-off losses for missions with durations greater than several months. These losses can be eliminated by incorporating high performance cryocooler technology to intercept heat load to the propellant tanks and modulating the cryocooler to control tank pressure. The active thermal control technology being developed by NASA is the reverse turbo-Brayton cycle cryocooler and its integration to the propellant tank through a distributed cooling tubing network coupled to the tank wall. This configuration was recently tested at NASA Glenn Research Center, in a vacuum chamber and cryo-shroud that simulated the essential thermal aspects of low Earth orbit, its vacuum and temperature. Testing consisted of three passive tests with the active cryo-cooler system off, and 7 active tests, with the cryocooler powered up. The test matrix included zero boil-off tests performed at 90 full and 25 full, and several demonstrations at excess cooling capacity and reduced cooling capacity. From this, the tank pressure response with varied cryocooler power inputs was determined. This test series established that the active cooling system integrated with the propellant tank eliminated boil-off and robustly controlled tank pressure.

Cryogenic fluid management↗

NASA NTRSDate not supplied

Modeling and Simulation of Tank Pressure Control using Zero-Boiloff 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↗

NASA NTRSDate not supplied

Improved Models for Simulation of Cryogenic Spray Cooling for Propellant Tank Chilldown

This paper presents advances in spray boiling models for cryogenic fluids that exhibit strong variation in thermodynamic properties; it accounts for the interaction of the spray droplets with walls and the complex heat transfer and boiling processes that play a vital role in chilling solid surfaces. The model has been developed by corroborating and validating it with test data in liquid nitrogen obtained at the University of Connecticut. The test data includes detailed visualization and plate temperature measurements for liquid nitrogen sprays impinging on a warm test plate and subsequent transient chilldown of the plate. Simulations of transient plate chilldown in liquid nitrogen were carried out and temperature histories at various locations on the plate were compared with test measurements to demonstrate the validity of the spray boiling models that have been developed for cryogens.

Boiling↗

NASA NTRSDate not supplied

Improved Models for Simulation of Cryogenic Spray Cooling for Propellant Tank Chilldown

Boiling↗

DOE OSTI2016

Clad to coolant heat transfer during RIA in PWR conditions

The power pulse during an insertion of reactivity accident (RIA) leads to a transient heating of the fuel. Heat transfer toward the coolant through the clad and associated temperature transient of the clad is of safety significance since it may largely impact the risk of clad failure. For sufficiently large enthalpy deposited within the fuel, clad temperature increases rapidly and this can lead to coolant boiling. This paper concerns the critical review of the present understanding and of the models for the clad to coolant heat transfer in RIA for PWR conditions. The results of past experiments are reviewed and the reliability of current models dedicated to these conditions is analyzed. It appears that the use of classical models for wall to fluid heat transfer does not allow predicting the clad temperature transient as measured in RIA-related experiments. The deviation is all the more large when film boiling conditions are reached. To reproduce the experimental data, some correlations of heat transfer can be fitted but the reliability of such model out of the range of conditions of the experimental data is low. It is mainly due to a rather low level of understanding of the heat transfer in very transient conditions (i.e. when the time scale of variation of the main heat transfer parameters, like a wall temperature, is less than the typical internal time scale of the heat transfer process itself). It results in a very large scattering in the numerical simulation of clad temperature when several fuel codes are considered. To improve the reliability of these calculations, it is believed that the study of transient boiling heat transfer thanks to devoted experimental facilities is required.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2016

Uncertainty Quantification with MCMC method and BFBT Data

The reliability of predictions of the system codes is closely related to the validation of their physical models. For example, the accuracy of void fraction prediction in a Boiling Water Reactor (BWR) using TRACE code depends on the uncertainties of closure relations developed for the two-phase flow models. This work focus on quantifying the uncertainties of two physical models, interfacial drag model and subcooled boiling model, using Markov Chain Monte Carlo (MCMC) method and BFBT experimental benchmark data. We applied the MCMC and MLE algorithms to BFBT benchmark data and estimated the uncertainty distribution of two model parameters: subcooled boiling heat transfer coefficient and interfacial drag (bubbly/slug rod bundle) coefficient. These algorithms tend to give various estimation results. Despite the differences in these algorithms, we find that the estimation results are consistent. After obtaining the uncertainty estimates, we validated the uncertainty estimates by applying the estimated distributions of model parameters to TRACE prediction. The new predictions succeed in reducing the void fraction prediction error, but it is not possible to eliminate all prediction error. We are only adjusting two model parameters, and there exists uncertainties in all model parameters. Application of the algorithms to other physical models and experimental data will be helpful for future work. (authors)

42 ENGINEERING↗

DOE OSTI2016

Contaminant Entrainment in a Liquid Fuel Fire with Multi-Component Evaporation Droplet Model

A number of fire scenarios can result in the release of hazardous chemical or radiological contaminants as sub-micron solid particulates, posing a significant danger of inhalation. These particulates may be entrained in the air around the fire through a variety of mechanisms including liquid droplet formation above the boiling surface of the contaminated fuel. The formation of both film drops from the rupture of a liquid film as a boiling bubble bursts and jet drops from the collapse of the cavity of a ruptured bubble may be involved. A 2015 study by Brown and Louie focused on such a fire scenario to investigate the physics of contaminant release. This paper extends that study by supplementing the existing physics models with a multicomponent evaporation mechanism, where liquid droplets formed near the boiling fuel surface experience evaporation of their various components at different rates based on distinct material properties (density, specific heat, heat of vaporization, etc.). As in the Brown and Louie study, the Sierra low-Mach module Fuego is used as the Lagrangian-Eulerian code for simulation of this scenario. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2016

CTF Parallel Performance Improvements

CTF is a modernized version of COBRA-TF (Coolant Boiling in Rod Arrays-Two Fluid), which is a thermal hydraulic (T/H) sub-channel code based on a two-phase, two-fluid model that formulates the conservation equations of mass, energy, and momentum for three fields of vapor, continuous liquid, and entrained liquid droplets. Under the Consortium for Advanced Simulation of Light Water Reactors (CASL) program sponsored by the U.S. Department of Energy, CTF is being further developed and improved for modeling and simulation of light water reactors. The code improvements for predicting T/H responses of pressurized water reactor (PWR) cores are focused on its applications to specific challenge problems deemed significant to the industry, such as departure from nucleate boiling (DNB) and crud-induced power shift. The code is also being used for T/H feedback in normal operating condition depletions in PWR and boiling water reactor designs in the reactor core simulator, Virtual Environment for Reactor Applications Core Simulator (VERA-CS), being developed by CASL. CTF modifications under the CASL program include software optimization, new closure models, and parallelization for modeling full reactor core T/H responses. The purpose of this work was to improve the parallel performance of CTF. CTF was originally parallelized via a domain-decomposition strategy using the Message Passing Interface (MPI), in which the basis of a solution domain was a physical assembly in the model. In other words, a user who wanted to perform a parallel solution with CTF had to provide the code one processor per assembly in the model (e.g., 56 cores for a 56 assembly quarter-core model, or 193 cores for a 193 assembly full-core model). This design was rigid, demanding a very specific number of processors for a model; furthermore, when CTF was run in VERA-CS-which typically uses thousands of processors- many processors sat idle during the CTF solve because CTF was unable to use any more processors than there are physical assemblies. This was not necessarily a problem if the CTF execution speed was very fast, but the average time spent in a CTF solve for prototypical PWR operating conditions was significant (i.e., over 20 % of the simulation wall time was being spent in CTF for a standard core depletion simulation). Larger, more computationally intensive simulations performed in support of the CASL DNB challenge problem demonstrated a need for improvements in this original parallelization approach. The specific goal of this work was to allow CTF to use more processors than there are physical assemblies in the model, which would not only speed up CTF and improve the overall efficiency of VERA-CS but also give users more-flexible options for executing their simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2017

Newton-Krylov method in solving drift-flux two-phase flow problems for both pre-CHF and post-CHF conditions at steady state

Furthermore, this paper presents numerical simulations of two-phase flow heat transfer for both pre-CHF and post-CHF conditions typical in light water nuclear reactors. The conjugate heat transfer problem is modeled with one-dimensional drift-flux two-phase flow model for coolant flow and two-dimensional heat conduction equation in solid walls. Fully implicit time integration schemes were employed, and the resulted non-linear equations were solved with a Newton-Krylov method. For the drift-flux two-phase flow model, the EPRI drift-flux closure correlation is used to model the relative motion between the two phases. Additional closure correlations were adopted from REALP5-3D to determine two-phase flow regime, wall boiling heat and mass transfer, interfacial heat and mass transfer, and two-phase flow frictional pressure drop. Three sets of experiments were used for code validation, including Bartolomei subcooled flow boiling in vertical pipes, FRIGG boiling tests in vertical bundles, and Bennett heated tube tests. RELAP5-3D simulation results were also provided to perform code-to-code benchmark. Overall, numerical results of this work showed good agreements with both experimental data and RELAP5-3D simulation results. Discrepancy between them were also identified in post-CHF region, which suggests necessary further improvement.

42 ENGINEERING↗

DOE OSTI2022

A remarkable CHF of 345W/cm 2 is achieved in a wicked-microchannel using HFE-7100

Boiling heat transfer of dielectric fluids is a promising cooling technique for thermal management of microelectronic systems. However, the critical heat flux achievable is generally low, because of the poor thermophysical properties of these fluids. Here, to address this dilemma, we propose a new cooling concept to substantially enhance global liquid supply during phase-change process using enhanced capillary-driven force. Additionally, dedicated vapor pathways are designed among a bank of micro-pillars to facilitate vapor removal. In this work, new wicks comprised of silicon micro-pinfin arrays are explored to significantly enhance the flow boiling heat transfer performance. To examine the functionalities of this wick, experiments on HFE-7100 were carried out with mass velocities varying from 247 to 3,465 kg/m 2 s. To explore the enhancement mechanisms and to analyze the capillary-assisted flow boiling process, visualization studies were conducted. The results indicate that sustainable evaporation induced by wick microstructures and efficient liquid supply are the enhancement mechanisms compared to parallel microchannels with solid walls. It is found that the overall heat transfer coefficient is substantially increased up to 75%. Remarkably, a high critical heat flux (CHF) of approximately 345 W/cm 2 is recorded at G = 3,465 kg/m 2 s at coolant inlet temperature of ~20 °C. Equally importantly, this noticeable enhancement of CHF value is associated with drastically decreased pressure drops compared to microchannels decorated with μ-pinfin fences.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Liquid hydrogen pump performance and durability testing through repeated cryogenic vessel filling to 700 bar

Here, this paper reports the results of a comprehensive test of liquid hydrogen (LH 2 ) pump performance and durability conducted while cycle testing a prototype thin-lined cryogenic pressure vessel 456 times to 700 bar. This extensive LH2 pump experimental data set provides a wealth of information vital for a complete evaluation of the future potential of this promising technology for ambient temperature and cryo-compressed vessel refueling. The experiment was conducted at Lawrence Livermore National Laboratory (Livermore, CA)'s hydrogen test facility, specifically built for this experiment and including a containment vessel for safe testing of the prototype vessel and a control room for efficient monitoring. Original pump and storage instrumentation was complemented with an electric power analyzer and a boil-off mass flow meter for more complete pump characterization. The results of the experiment confirm most of the expected virtues of the LH 2 pump: rapid (3 min) refueling of the 65-liter prototype vessel at high flow rate (1.55 kgH 2 per minute on average), unlimited back to back refueling, low electricity consumption (1.1 kWh/kg H 2 ), no measurable degradation, and low maintenance. High cryogenic vessel fill density is another key performance metric that was demonstrated in an earlier publication. These virtues derive from the high density of LH 2 enabling pressurization to high density with relatively little energy consumption and high throughput from a small displacement pump. Boil-off losses as high as 27.7% of dispensed hydrogen were measured, at experimental conditions not representative of operation at a hydrogen refueling station. These losses drop to 15.4% for operation that may be representative of a small station (332 kg/day), while we anticipate less than 6% boil-off with the existing pump and Dewar with improved LH 2 delivery truck operations and a more favorable arrangement of the LH 2 pump relative to the Dewar.

08 HYDROGEN↗

DOE OSTI2022

Formation of a super-dense hydrogen monolayer on mesoporous silica

Adsorption on various adsorbents of hydrogen and helium at temperatures close to their boiling points shows, in some cases, unusually high monolayer capacities. The microscopic nature of these adsorbate phases at low temperatures has, however, remained challenging to characterize. Here, using high-resolution cryo-adsorption studies together with characterization by inelastic neutron scattering vibration spectroscopy, we show that, near its boiling point (~20 K), H 2 adsorbed on a well-ordered mesoporous silica forms a two-dimensional monolayer with a density more than twice that of bulk-solid H 2 , rather than a bilayer. Theoretical studies, based on thorough first-principles calculations, rationalize the formation of such a super-dense phase. The strong compression of the hydrogen surface layer is due to the excess of surface–hydrogen attraction over intermolecular hydrogen repulsion. Use of this super-dense hydrogen monolayer on an adsorbent might be a feasible option for the storage of hydrogen near its boiling point, compared with adsorption at 77 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

High-speed imaging of transition from fluid breakup to phase explosion in electric explosion of tungsten wires in air

High-speed visible imaging of sub-microsecond electric explosion of wires at the low specific energy deposition threshold reveals three distinct modes of wire failure as capacitor charge voltage and energy deposition are increased. For 100 micron diameter gold-plated tungsten wires of 2 cm length, deposited energies of 1.9 eV/atom produces a liquid column that undergoes hydrodynamic breakup into droplets with radii of order of wire diameter on timescales of 200 μs. Instability growth, column breakup, and droplet coalescence follow classic Rayleigh-Plateau predictions for instability of viscous fluid column. Above 3.2 eV/atom of deposited energy, wires are seen to abruptly transition to an expanding mixture of micron scale liquid-droplets and vapor within one frame (less than 3.33 μs), which has been termed ‘phase explosion’ in literature. Between these two limits, at 2.5 eV/atom of deposited energy, wire radius is unchanged for the first 10 μs before the onset of a rapid expansion and disintegration that resembles homogenous nucleation of mechanically unstable bubbles. Here, thermodynamic calculations are presented that separate cases by temperature obtained during heating: below boiling point, near boiling point, and exceeding boiling point.

42 ENGINEERING↗

DOE OSTI2024

Stochastic machine learning via sigma profiles to build a digital chemical space

This work establishes a different paradigm on digital molecular spaces and their efficient navigation by exploiting sigma profiles. To do so, the remarkable capability of Gaussian processes (GPs), a type of stochastic machine learning model, to correlate and predict physicochemical properties from sigma profiles is demonstrated, outperforming state-of-the-art neural networks previously published. The amount of chemical information encoded in sigma profiles eases the learning burden of machine learning models, permitting the training of GPs on small datasets which, due to their negligible computational cost and ease of implementation, are ideal models to be combined with optimization tools such as gradient search or Bayesian optimization (BO). Gradient search is used to efficiently navigate the sigma profile digital space, quickly converging to local extrema of target physicochemical properties. While this requires the availability of pretrained GP models on existing datasets, such limitations are eliminated with the implementation of BO, which can find global extrema with a limited number of iterations. A remarkable example of this is that of BO toward boiling temperature optimization. Holding no knowledge of chemistry except for the sigma profile and boiling temperature of carbon monoxide (the worst possible initial guess), BO finds the global maximum of the available boiling temperature dataset (over 1,000 molecules encompassing more than 40 families of organic and inorganic compounds) in just 15 iterations (i.e., 15 property measurements), cementing sigma profiles as a powerful digital chemical space for molecular optimization and discovery, particularly when little to no experimental data is initially available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Characterization of Fuel-to-Coolant Heat Transfer During Reactivity-Initiated Accidents Using Tightly Coupled Thermal Hydraulics and Fuel Thermomechanics

The reactivity-initiated accident (RIA) is a complex scenario with several tightly interacting physical phenomena. Accurately predicting fuel behavior during these transients is difficult due to limitations in the modeling of fuel-to-coolant heat transfer. Common approaches to simulate RIAs involve standalone calculations using either a fuel performance code or a thermal-hydraulic code. The complex interdependencies of thermal-hydraulic and fuel mechanical behavior suggest that a tight coupling between these codes may provide more accurate predictions of fuel-to-coolant heat transfer and cladding mechanical response. Here, RELAP5-3D and BISON are coupled in this paper to simulate RIAs, and a sensitivity analysis is performed to rank key thermal properties and two-phase heat transfer parameters relevant for fuel-to-coolant heat transfer and cladding failure mechanisms in UO 2 –Zircaloy-4 systems. Gas gap conductance, film boiling heat transfer uncertainty, pulse width, fuel-specific heat capacity, and cladding-specific heat capacity were identified as important parameters. Variations in figures of merit resulting from changes to pulse width and the material thermal properties indicate that time-dependent heat transfer rates are significant for safety-relevant mechanical parameters due to the time dependence of cladding ductility and pellet-cladding mechanical interaction loading. The results suggest that the thermal-hydraulic factors have a nonnegligible influence on the thermomechanical solution and vice versa. Tight coupling of both sets of physics is recommended to improve prediction of fuel behavior during RIAs. Highlights include the following: 1. The RELAP5-3D thermal-hydraulic code and the BISON fuel performance code are tightly coupled for simulation of RIA transients with energy depositions at the Zircaloy-4 cladding failure threshold. 2. Departure from nucleate boiling occurred for all simulated cases. Due to the ductility of fresh fuel, substantial ballooning occurred in most cases. 3. Gas gap conductance, fuel-specific heat capacity, cladding-specific heat capacity, transient pulse width, and film boiling heat transfer were the dominant thermal factors impacting the safety figures of merit at energy depositions.

Critical Heat Flux (CHF)↗

DOE OSTI2018

Small-scale cook-off experiments and models of ammonium nitrate

Here, we have completed a series of small-scale cook-off experiments of ammonium nitrate (AN) prills in our Sandia Instrumented Thermal Ignition test at nominal packing densities of about 0.8 g/cm 3 . We increased the boundary temperature of our aluminum confinement cylinder from room temperature to a prescribed set-point temperature in 10 min. Our set-point temperature ranged from 508 to 538 K. The external temperature of the confining cylinder was held at the set-point temperature until ignition. We used type K thermocouples to measure temperatures associated with several polymorphic phase changes as well as melting and boiling. As the AN boiled, our thermocouples were destroyed by corrosion, which may have been caused by reaction of hot nitric acid (HNO 3 ) with nickel to form nickel nitrate, Ni(NO 3 ) 2 . Videos of the corroding thermocouples showed a green solution that was similar to the color of Ni(NO 3 ) 2 . We found that ignition was imminent as the AN boiling point was exceeded. Ignition of the AN prills was modeled by solving the energy equation with an energy source due to desorption of moisture and decomposition of AN to form equilibrium products. A Boussinesq approximation was used in conjunction with the momentum equation to model flow of the liquid AN. We found that the prediction of ignition was not sensitive to small perturbations in the latent enthalpies.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2015

DNS of turbulent flow with hemispherical wall roughness

The present study of the effect of roughness density on the mean flow turbulence parameters is motivated by the need for new generation of boundary conditions for multiphase computational multiphase fluid dynamics (CMFD) models applied to boiling flows. Effect of roughness element density on the turbulent flow in a channel is quantified through direct numerical simulations (DNSs). The Navier--Stokes equations are solved using finite element method and bubbles are approximated as rigid near-hemispherical obstacles at the wall. Six different cases were analysed including channel flow with smooth wall and channel flow with rough wall for five different bubble nucleation site densities. Friction factor and the law of the wall was calculated and compared with the previously published results. Existing correlations for nucleating bubble site density dependency on a wall heat flux were used to obtain a relation between the heat flux and the friction factor, leading to the law of the wall dependency on the heat flux. This separate effect study provides new guidelines on how the heat flux in subcooled boiling regime affects the turbulence behaviour near the wall and guides the computational fluid dynamics model development for boiling two-phase flows.

Mechanics↗

DOE OSTI2020

effect of ambient pressure on the behavior of single-component fuels in a gasoline multi-hole injector

The injection characteristics of neat ethanol and pure isooctane are studied under different ambient pressure and temperature conditions. Injection underflash-boiling conditions can enhance liquid atomization and evaporation, providing the possibility of improvement in the fuel/air mixing. These superheated conditions often introduce phenomena that are not taken into account in the standard modeling of sprays for engine applications. The present work proposes a numerical investigation of the behavior of Engine Combustion Network's 8-hole spray-G injector, starting at the subcooled nominal condition and reducing the ambient pressure at constant low temperature to reach the flare flash-boiling condition. To initialize the properties of the injected fuel, the flow in the nozzle is simulated with a Eulerian approach, handling the two phases with a mixture model and the phase change, due to cavitation and flash boiling, with the Homogenous Relaxation Model. A map of the mixture's kinematic and thermal behavior is obtained at the interface between the injector and the chamber to initialize the Lagrangian simulations. A literature-based vaporization model is implemented to obtain the proper description of the characteristic features of a multi-hole spray under super-heated conditions, like plume-plume interaction. The numerical representation of the spray is validated in terms of penetration and radial spreading on DBI images, reproducing the light attenuation profiles caused by the presence of the liquid spray. Simulations show that coupled nozzle flow and spray calculations capture the spray morphology and shape better compared to calculations performed without considering the nozzle flow simulation details, especially under flare flash conditions.

Nocivelli, Lorenzo↗