Stratified layer flow model - a numerical approach to liquid temperature stratification.
Temperature stratification in liquids contained in heated vessels studied in stratified layer flow model, using turbulent layer equations
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Temperature stratification in liquids contained in heated vessels studied in stratified layer flow model, using turbulent layer equations
A reduced dynamical model describing temperature stratification effects driven by natural convection in a liquid hydrogen cryogenic fuel tank has been developed. It accounts for cryogenic propellant loading, storage, and unloading in the conditions of normal, increased, and micro- gravity. The model involves multiple horizontal control volumes in both liquid and ullage spaces. Temperature and velocity boundary layers at the tank walls are taken into account by using correlation relations. Heat exchange involving the tank wall is considered by means of the lumped-parameter method. By employing basic conservation laws, the model takes into consideration the major multi-phase mass and energy exchange processes involved, such as condensation-evaporation of the hydrogen, as well as flows of hydrogen liquid and vapor in the presence of pressurizing helium gas. The model involves a liquid hydrogen feed line and a tank ullage vent valve for pressure control. The temperature stratification effects are investigated, including in the presence of vent valve oscillations. A simulation of temperature stratification effects in a generic cryogenic tank has been implemented in Matlab and results are presented for various tank conditions.
In recent years, the Engine Combustion Network (ECN) has developed as a worldwide reference for understanding and describing engine combustion processes, successfully bringing together experimental and numerical efforts. Since experiments and numerical simulations both target the same boundary conditions, an accurate characterization of the stratified environment that is inevitably present in experimental facilities is required. The difference between the core-, and pressure-derived bulk-temperature of pre-burn combustion vessels has been addressed in various previous publications. Additionally, thermocouple measurements have provided initial data on the boundary layer close to the injector nozzle, showing a transition to reduced ambient temperatures. The conditions at the start of fuel injection influence physicochemical properties of a fuel spray, including near nozzle mixing, heat release computations, and combustion parameters. To address the temperature stratification in more detail, thermocouple measurements at larger distances from the spray axis have been conducted. Both the temperature field prior to the pre-combustion event that preconditions the high-temperature, high-pressure ambient, as well as the stratification at the moment of fuel injection were studied. To reveal the cold boundary layer near the injector with a better spatial resolution, Rayleigh scattering experiments and thermocouple measurements at various distances close to the nozzle have been carried out. The impact of the boundary layers and temperature stratification are illustrated and quantified using numerical simulations at Spray A conditions. Next to a reference simulation with a uniform temperature field, six different stratified temperature distributions have been generated. These distributions were based on the mean experimental temperature superimposed by a randomized variance, again derived from the experiments. The results showed that an asymmetric flame structure arises in the computed results when the temperature stratification input is used. In these predictions, first-stage ignition is advanced by 24 $\mu$s , while second-stage ignition is delayed by 11 $\mu$s. At the same time a lift-off length difference between the top and the bottom of up to 1.1 mm is observed. Furthermore, the lift-off length is less stable over time. Given the shown dependency, the temperature data is made available along with the vessel geometry data as a recommended basis for future numerical simulations.
During heating of a water heater tank, the vertical temperature stratification of the water can be increased or decreased, depending on the method of heating. Methods that increase stratification during heating include (1) removing cold water from the tank bottom, heating it, and re-introducing it to the tank top at relatively low flow rate, (2) using a heat exchanger wrapped around the tank, through which heating fluid (with finite specific heat) flows from top to bottom, and (3) using an immersed heat element that is relatively high in the tank. Using such methods allows for improved heat pump water heater (HPWH) cycle efficiencies when the heat pump can take advantage of the lower temperatures that exist lower in the tank, and accommodate the resulting glide. Transcritical cycles are especially well-suited to capitalize on this opportunity, and other HPWH configurations (that have been proposed elsewhere) may benefit as well. This work provides several stratification categories of heat pump water heater tank configurations relevant to their stratification potential. To illustrate key differences among categories, it also compiles available experimental data for (a) single pass pumped flow, (b) multi-pass pumped flow, and (c) top-down wrapped tank with transcritical refrigerant.
Ion kinetic temperature measured in E region by ion trap aboard Nike-Apache rocket
In recent years, the processes in the atmospheric planetary boundary layer (PBL) over urban areas were intensely investigated, due to ecological problems related to the air, soil, and water pollution. New pollution sources in new residential districts, when in contradiction to the microclimate and topography requirements of that region, create a number of considerable hazards and problems. The present study is a continuation of our preceding investigations and aims at revealing the aerosol structure and stratification during the transition after sunset as measured by two lidars. Such observation of the nocturnal, stable PBL formation over an urban area in Bulgaria has not been reported before. The lidars' high time and spatial resolutions allow the changes of the internal structure of the PBL's part located above the surface layer to be observed.
Temperature and pressure data obtained from the saturated Freon 113 PCA closed-tank stratification tests are presented. The data presented in tabular form are the test conditions, sensible heat values, and Freon 113 PCA liquid and ullage (vapor) properties. Also included, are graphical representations of the liquid bulk temperature and pressure histories, and dimensionless liquid-ullage delta-temperature profiles. Modified Grashof numbers and Fourier number-history data are also presented graphically.
The basic equations describing the energy equilibria, the conductive heat flux, and the temperature stratifications for stellar transition layers and coronae with 'open' field lines are summarized. The temperature dependence of the emission measures for lines originating in different temperature regions of the transition zone is determined. It is found that the stellar transition regions consist of two basically different parts: the lower part where mechnical energy input is balanced by the radiative losses, and the upper part where the mechanical energy input is balanced by the divergence of the conductive flux and radiative losses. In the lower part, the temperature stratification is determined by an equilibrium between mechanical flux input and radiative energy losses. The coronal temperatures increase with increasing mechanical flux and damping length in the upper transition zone.
Abstract This paper presents a sensible-enthalpy-based conditional moment closure (CMC) model for homogeneous charge compression ignition (HCCI) with temperature inhomogeneity. The focus of study is placed on the effects of a conditioning variable on the prediction of heat release rates and developing the mixing models for sensible enthalpy, a reacting scalar. The sensible-enthalpy-based CMC model is validated a priori with 2-D direct numerical simulation (DNS) data for HCCI, and compared with the total-enthalpy-based formulation. Sensible enthalpy is shown to be strongly correlated with temperature, which makes sensible enthalpy well suited for a conditioning variable in CMC for HCCI. Results show that sensible-enthalpy-based CMC performs well for the prediction of heat release rates for all temperature stratification levels investigated in this study, while total-enthalpy-based CMC leads to the overprediction of heat release rates when the temperature stratification level is high. While conditional fluctuations of species mass fractions are substantial, sensible-enthalpy-based CMC is found to accurately estimate heat release rates during ignition, as well as reaction rates of key species, due to substantially suppressed conditional temperature fluctuations. The mapping closure models for the probability density function (PDF) and conditional dissipation rates of sensible enthalpy are proposed. The proposed models are found to capture the characteristics of reactive scalar mixing.
We study the C8.4-class solar flare SOL2016-05-14T11:34 UT using high-resolution spectral imaging in the Ca ii 8542 Å line obtained with the CRISP imaging spectropolarimeter on the Swedish 1 m Solar Telescope. Spectroscopic inversions of the Ca ii 8542 Å line using the non-LTE code NICOLE are used to investigate the evolution of the temperature and velocity structure in the flaring chromosphere. A comparison of the temperature stratification in flaring and non-flaring areas reveals strong footpoint heating during the flare peak in the lower atmosphere. The temperature of the flaring footpoints between log τ{sub 500} ≈−2.5 and −3.5, where τ {sub 500} is the continuum optical depth at 500 nm, is ∼5--6.5 kK close to the flare peak, reducing gradually to ∼5 kK. The temperature in the middle and upper chromosphere, between log τ{sub 500}≈−3.5 and −5.5, is estimated to be ∼6.5–20 kK, decreasing to preflare temperatures, ∼5–10 kK, after approximately 15 minutes. However, the temperature stratification of the non-flaring areas is unchanged. The inverted velocity fields show that the flaring chromosphere is dominated by weak downflowing condensations at the formation height of Ca ii 8542 Å.
This paper presents the development of a two-phase CFD model that is used to study sloshing of a cryogenic fluid (LN2) undergoing phase change generated by lateral acceleration and its effects on the ensuing heat and mass transfer and pressure drop in a small-scale tank. In this model the interface is captured using the volume-of-fluid (VOF) method [1]. Kinetic-theory based Schrage relation [2] is used for calculating phase change mass transfer at the liquid-vapor interface. Computational results are compared to the data provided by a non-isothermal sloshing experiment with phase change conducted by the University of Tokyo and JAXA in 2018 using liquid Nitrogen in a transparent tank. The effect of different thicknesses of a liquid temperature stratification layer created at the interface on the pressure drop during sloshing was studied experimentally. CFD results for interface movement and tank pressure are presented and compared in this paper to experimental data for the cases with the smallest and largest liquid temperature stratification thickness.
As the offshore wind energy industry begins to emerge on the U.S. northeast coast, there is an urgent need for an accurate characterization of the offshore wind resource. The offshore regime features various resource characteristics not seen onshore, including frequent coastal low-level jets, non-steady surface conditions (e.g. waves), complex air-sea and land-sea interactions, and coastal mesoscale circulations (e.g. sea breeze). An accurate characterization of these phenomena in mesoscale models is critical to support accurate energy production estimates and forecasting for future offshore wind farms. However, a robust validation of mesoscale model performance in the offshore U.S. has been limited due to sparse data and difficulty obtaining offshore measurements. Consequently, it is unclear how mesoscale models perform offshore and which suite of parameterizations and input data are best suited for accurate resource characterization. To begin addressing this knowledge gap, the National Renewable Energy Laboratory, in collaboration with the Rutgers Center for Ocean Observing Leadership, conducted a mesoscale model validation campaign in the offshore environment. The Weather Research and Forecasting (WRF) model was assessed and the New Jersey offshore coastal area served as the validation domain. An ensemble of WRF simulations conducted over a one-year period sampled from a range of model setups, including different reanalysis and forecast data inputs, sea surface temperature data inputs, planetary boundary layer parameterizations, and WRF versions. The performance of each ensemble member was validated against a network of observations including offshore floating lidars, coastal met towers, offshore platform measurements, and offshore buoys. Model performance was assessed on annual, seasonal, and diurnal scales with particular emphasis on the vertical wind speed profile and temperature stratification. Preliminary results from this validation work will be presented here. The key WRF inputs and parameters driving changes in the modeled wind resource will be discussed and preliminary recommendations on WRF setups for offshore wind resource modeling will be proposed. The application of this validation exercise to other proposed offshore wind energy areas in the U.S. will also be discussed. mesoscale circulations (e.g. sea breeze). An accurate characterization of these phenomena in mesoscale models is critical to support accurate energy production estimates and forecasting for future offshore wind farms. However, a robust validation of mesoscale model performance in the offshore U.S. has been limited due to sparse data and difficulty obtaining offshore measurements. Consequently, it is unclear how mesoscale models perform offshore and which suite of parameterizations and input data are best suited for accurate resource characterization. To begin addressing this knowledge gap, the National Renewable Energy Laboratory, in collaboration with the Rutgers Center for Ocean Observing Leadership, conducted a mesoscale model validation campaign in the offshore environment. The Weather Research and Forecasting (WRF) model was assessed and the New Jersey offshore coastal area served as the validation domain. An ensemble of WRF simulations conducted over a one-year period sampled from a range of model setups, including different reanalysis and forecast data inputs, sea surface temperature data inputs, planetary boundary layer parameterizations, and WRF versions. The performance of each ensemble member was validated against a network of observations including offshore floating lidars, coastal met towers, offshore platform measurements, and offshore buoys. Model performance was assessed on annual, seasonal, and diurnal scales with particular emphasis on the vertical wind speed profile and temperature stratification. Preliminary results from this validation work will be presented here. The key WRF inputs and parameters driving changes in the modeled wind resource will be discussed and preliminary recommendations on WRF setups for offshore wind resource modeling will be proposed. The application of this validation exercise to other proposed offshore wind energy areas in the U.S. will also be discussed.
Laminar models agree closely with the pressure evolution and vapor phase temperature stratification but under-predict liquid temperatures. Turbulent SST k-w and k-e models under-predict the pressurization rate and extent of stratification in the vapor but represent liquid temperature distributions fairly well. These conclusions seem to equally apply to large cryogenic tank simulations as well as small scale simulant fluid pressurization cases. Appropriate turbulent models that represent both interfacial and bulk vapor phase turbulence with greater fidelity are needed. Application of LES models to the tank pressurization problem can serve as a starting point.
Hydrostatic adjustment due to diabatic heat in two nonisothermal atmospheres is examined. In the first case the temperature stratification is continuous; in the second case the atmosphere is composed of a warm, isothermal troposphere and a colder, isothermal semi-infinitely deep stratosphere.In both cases hydrostatic adjustment, to a good approximation, follows the pattern found in the Lamb problem (semi-infinitely deep. isothermal atmosphere): Initially we have acoustic waves with the kinetic energy increasing or decreasing at the expense of available elastic energy. After this initial period the acoustic waves evolve into acoustic-gravity waves with the kinetic, available potential and available elastic energies interacting with each other. Relaxation to hydrostatic balance occurs within a few oscillations. Stratification in an atmosphere with a continuous temperature profile affects primarily the shape and amplitude of the disturbances. In the two-layer atmosphere, a certain amount of energy is trapped in the tropospheric waveguide as disturbances reflect off the tropopause and back into the troposphere. With each internal reflection a portion of this trapped energy escapes and radiates to infinity.
Radiant heating and cooling systems have significant energy-saving potential and are gaining popularity in commercial buildings. The main aim of the experimental study reported here was to characterize the behavior of radiant cooling systems in a typical office environment, including the effect of ceiling fans on stratification, the variation in comfort conditions from perimeter to core, control on operative temperature vs. air temperature and the effect of carpet on cooling capacity. The goal was to limit both the first cost and the perceived risk associated with such systems. Two types of radiant systems, the radiant ceiling panel (RCP) system and the radiant slab (RS) system, were investigated. The experiments were carried out in one of the test cells that constitute the FLEXLAB test facility at the Lawrence Berkeley National Laboratory in March and April 2016. In total, ten test cases (five for RCP and five for RS) were performed, covering a range of operational conditions. In cooling mode, the air temperature stratification is relatively small in the RCP, with a maximum value of 1.6 K. The observed stratification effect was significantly greater in the RS, twice as much as that in the RCP. The maximum increase in dry bulb temperature in the perimeter zone due to solar radiation was 1.2 K for RCP and 0.9 K for RS – too small to have a significant impact on thermal comfort. The use of ceiling fans was able to reduce any excess stratification and provide better indoor comfort, if required. The use of thin carpet requires a 1 K lower supply chilled water temperature to compensate for the added thermal resistance, somewhat reducing the opportunities for water-side free cooling and increasing the risk of condensation. In both systems, the difference between the room operative temperature and the room air temperature is small when the cooling loads are met by the radiant systems. In conclusion, this makes it possible to use the air temperature to control the radiant systems in lieu of the operative temperature, reducing both first cost and maintenance costs.
The equation of radiative equilibrium for rapidly expanding spherical shells including the special relativistic terms is discussed, and a temperature-correction procedure for the solution of this equation is given. The temperature-correction scheme combines a partial linearization approach with the approximate Lambda-operator iteration method. For some test problems, the temperature structures of atmospheres in radiative equilibrium in the Lagrangian frame are calculated. The effects of the velocity field on the temperature stratification and the luminosity are also discussed. We find that the temperature structure of a moving photosphere in radiative equilibrium in the Lagrangian frame differs significantly from that of a static photosphere and that, at least, all first-order terms, in particular the advection and aberration terms, have to be included consistently both in the radiative transfer and the energy equation in studies of supernovae atmospheres and other high-speed radiating flows.
This paper presents results from one and two-dimensional direct numerical simulations under Reactivity Controlled Compression Ignition (RCCI) conditions of a primary reference fuel (PRF) mixture consisting of n-heptane and iso-octane. RCCI uses in-cylinder blending of two fuels with different autoignition characteristics to control combustion phasing and the rate of heat release. These simulations employ an improved model of compression heating through mass source/sink terms developed in a previous work by Bhagatwala et al. (2014), which incorporates feedback from the flow to follow a predetermined experimental pressure trace. Two-dimensional simulations explored parametric variations with respect to temperature stratification, pressure profiles and n-heptane concentration. Furthermore, statistics derived from analysis of diffusion/reaction balances locally normal to the flame surface were used to elucidate combustion characteristics for the different cases. Both deflagration and spontaneous ignition fronts were observed to co-exist, however it was found that higher n-heptane concentration provided a greater degree of flame propagation, whereas lower n-heptane concentration (higher fraction of iso-octane) resulted in more spontaneous ignition fronts. A significant finding was that simulations initialized with a uniform initial temperature and a stratified n-heptane concentration field, resulted in a large fraction of combustion occurring through flame propagation. The proportion of spontaneous ignition fronts increased at higher pressures due to shorter ignition delay when other factors were held constant. For the same pressure and fuel concentration, the contribution of flame propagation to the overall combustion was found to depend on the level of thermal stratification, with higher initial temperature gradients resulting in more deflagration and lower gradients generating more ignition fronts. Statistics of ignition delay are computed to assess the Zel’dovich (1980) theory for the mode of combustion propagation based on ignition delay gradients.
We present deep spectroscopic (3600–24700 Å ) X-shooter observations of the bright Herbig–Haro object HH 1, one of the best laboratories to study the chemical and physical modifications caused by protostellar shocks on the natal cloud. We observe atomic fine structure lines, H i and He i recombination lines and H{sub 2} ro-vibrational lines (more than 500 detections in total). Line emission was analyzed by means of Non-local Thermal Equilibiurm codes to derive the electron temperature and density, and for the first time we are able to accurately probe different physical regimes behind a dissociative shock. We find a temperature stratification in the range 4000 K ÷ 80,000 K, and a significant correlation between temperature and ionization energy. Two density regimes are identified for the ionized gas, a more tenuous, spatially broad component (density ∼10{sup 3} cm{sup −3}), and a more compact component (density ⩾ 10{sup 5} cm{sup −3}) likely associated with the hottest gas. A further neutral component is also evidenced, having a temperature ≲10,000 K and a density >10{sup 4} cm{sup −3}. The gas fractional ionization was estimated by solving the ionization equilibrium equations of atoms detected in different ionization stages. We find that neutral and fully ionized regions co-exist inside the shock. Also, indications in favor of at least partially dissociative shock as the main mechanism for molecular excitation are derived. Chemical abundances are estimated for the majority of the detected species. On average, abundances of non-refractory/refractory elements are lower than solar of about 0.15/0.5 dex. This indicates the presence of dust inside the medium, with a depletion factor of iron of ∼40%.