SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A formulation representing multicomponent-fuel (MC-fuel) composition as a Probability Distribution Function (PDF) depending on the molar weight is used to construct a model of a large number of MC-fuel drops evaporating in a gas flow, so as to assess the extent of fuel specificity on the vapor composition.
LES of a circular jet laden with evaporating liquid drops are conducted to assess computational-drop modeling and three different SGS-flux models: the Scale Similarity model (SSC), using a constant coefficient calibrated on a temporal mixing layer DNS database, and dynamic-coefficient Gradient and Smagorinsky models.
Cassini ISS observed multiple widespread changes in surface brightness in Titan's equatorial regions over the past three years. These brightness variations are attributed to rainfall from cloud systems that appear to form seasonally. Determining the composition of this rainfall is an important step in understanding the "methanological" cycle on Titan. I use data from Cassini VIMS to complete a spectroscopic investigation of multiple rain-wetted areas. I compute "before-and-after" spectral ratios of any areas that show either deposition or evaporation of rain. By comparing these spectral ratios to a model of liquid ethane, I find that the rain is most likely composed of liquid ethane. The spectrum of liquid ethane contains multiple absorption features that fall within the 2-micron and 5-micron spectral windows in Titan's atmosphere. I show that these features are visible in the spectra taken of Titan's surface and that they are characteristically different than those in the spectrum of liquid methane. Furthermore, just as ISS saw the surface brightness reverting to its original state after a period of time, I show that VIMS observations of later flybys show the surface composition in different stages of returning to its initial form.
Evaporation and condensation at a liquidvapor interface is important for long-term, in-space cryogenic propellant storage. Yet the current understanding of interfacial physics does not predict behavior or evaporation condensation rates. The proposed paper will present a physical model, based on the 1-D Heat equation and Schrages equation which demonstrates thin thermal layers at the fluidvapor interface.
The Gravity Recovery and Climate Experiment (GRACE) mission and its Follow-On (GRACE-FO) mission provide unprecedented observations of terrestrial water storage (TWS) dynamics at basin to continental scales. Established GRACE data assimilation techniques directly adjust the simulated water storage components to improve the estimation of groundwater, streamflow, and snow water equivalent. Such techniques artificially add/subtract water to/from prognostic variables, thus upsetting the simulated water balance. To overcome this limitation, we propose and test an alternative assimilation scheme in which precipitation fluxes are adjusted to achieve the desired changes in simulated TWS. Using a synthetic data assimilation experiment, we show that the scheme improves performance skill in precipitation estimates in general, but that it is more robust for snowfall than for rainfall, and it fails in certain regions with strong horizontal gradients in precipitation. The results demonstrate that assimilation of TWS observations can help correct (adjust) the model’s precipitation forcing and, in turn, enhance model estimates of TWS, snow mass, soil moisture, runoff, and evaporation. A key limitation of the approach is the assumption that all errors in TWS originate from errors in precipitation. Nevertheless, the proposed approach produces more consistent improvements in simulated runoff than the established GRACE data assimilation techniques.
Thermal equipment in buildings is a primary contributor to peak loads on the electrical grid. Thermal energy storage is a cost-effective strategy to decouple electric use from thermal loads, thus reducing grid peak costs for building owners. One method for storing thermal energy in a building is to integrate a phase change material (PCM) directly into the heating, ventilation, and air conditioning system. These systems often require additional glycol loops, pumps, valves, and heat exchangers to couple the storage to the cooling system and building space, which increases the complexity and cost. This work will discuss an alternate approach where the storage is added directly into the heat pump evaporator. A detailed two-dimensional finite difference heat transfer model of a PCM-refrigerant-glycol heat exchanger was developed to simulate the performance of this component. The fluid stream was discretized along the flow direction to capture changes in the fluid properties and local heat transfer rates, and the phase change material was discretized in both the x and y directions to capture the movement of the melt front. The model was used to understand the impact of different material and geometric properties on the charge and discharge characteristics of the device. Finally, a Ragone framework analogous to that used for electrochemical batteries was used to maximize the energy density and round-trip efficiency of the device while supplying loads appropriate for space cooling in buildings.
The Hanford Tank Waste Treatment and Immobilization Plant (WTP) pilot studies found that semi-volatile elements started to leave the molten glass pool and enter the offgas stream when the melter became idle. The rate at which the semi-volatile elements leave the molten glass pool follows a typical half-life decay curve. The WTP Direct Feed Low-Activity Waste (DFLAW) program is using existing LAW offgas treatment systems to capture semi-volatile elements leaving the melters. A new building, Effluent Management Facility (EMF), will be built to process effluents generated by the offgas treatment equipment. The effluents are concentrated by evaporation and are recycled to the melter feed streams. Computer modeling studies show that increased amounts of semi-volatile elements in the recycles due to melter downtime adversely impact the DFLAW operations. The extent to which glass productions and process efficiencies are hampered depends upon the frequency, duration, and pattern of the downtime events. Various scenarios were investigated using randomly generated and operator-anticipated downtime schedules. This paper discusses 1) derivation of the volatility half-life curves, 2) downtime schedules, 3) increased amount of semi-volatile elements in recycled material as a result of melter shut down, 4) impacts to glass production and process efficiency, and 5) comparison of scenario results. This study concludes that plant operators can mitigate much of the negative effects of semi-volatile element recycles by managing melter downtime and idle time when possible. (authors)
An analytic model for pressurization and cryogenic propellant conditions during all mission phases of any liquid rocket based vehicle has been developed and validated. The model assumes the propellant tanks to be divided into five nodes and also implements an empirical correlation for liquid stratification if desired. The five nodes include a tank wall node exposed to ullage gas, an ullage gas node, a saturated propellant vapor node at the liquid-vapor interface, a liquid node, and a tank wall node exposed to liquid. The conservation equations of mass and energy are then applied across all the node boundaries and, with the use of perfect gas assumptions, explicit solutions for ullage and liquid conditions are derived. All fluid properties are updated real time using NIST Refprop.1 Further, mass transfer at the liquid-vapor interface is included in the form of evaporation, bulk boiling of liquid propellant, and condensation given the appropriate conditions for each. Model validation has proven highly successful against previous analytic models and various Saturn era test data and reasonably successful against more recent LH2 tank self pressurization ground test data. Finally, this model has been applied to numerous design iterations for the Altair Lunar Lander, Ares V Core Stage, and Ares V Earth Departure Stage in order to characterize Helium and autogenous pressurant requirements, propellant lost to evaporation and thermodynamic venting to maintain propellant conditions, and non-uniform tank draining in configurations utilizing multiple LH2 or LO2 propellant tanks. In conclusion, this model provides an accurate and efficient means of analyzing multiple design configurations for any cryogenic propellant tank in launch, low-acceleration coast, or in-space maneuvering and supplies the user with pressurization requirements, unusable propellants from evaporation and liquid stratification, and general ullage gas, liquid, and tank wall conditions as functions of time.
A two-dimensional Heat Pipe Transient Analysis Model, 'HPTAM,' was developed to simulate the transient operation of fully-thawed heat pipes and the startup of heat pipes from a frozen state. The model incorporates: (a) sublimation and resolidification of working fluid; (b) melting and freezing of the working fluid in the porous wick; (c) evaporation of thawed working fluid and condensation as a thin liquid film on a frozen substrate; (d) free-molecule, transition, and continuum vapor flow regimes, using the Dusty Gas Model; (e) liquid flow and heat transfer in the porous wick; and (f) thermal and hydrodynamic couplings of phases at their respective interfaces. HPTAM predicts the radius of curvature of the liquid meniscus at the liquid-vapor interface and the radial location of the working fluid level (liquid or solid) in the wick. It also includes the transverse momentum jump condition (capillary relationship of Pascal) at the liquid-vapor interface and geometrically relates the radius of curvature of the liquid meniscus to the volume fraction of vapor in the wick. The present model predicts the capillary limit and partial liquid recess (dryout) in the evaporator wick, and incorporates a liquid pooling submodel, which simulates accumulation of the excess liquid in the vapor core at the condenser end.
This paper examines a novel concept of flow assisted latent heat driven two-phase evaporative cooling (EC) confined in-between slot liner and active-winding of electric motor. Wicking micro-structure enhanced PDMS liner axially sucks coolant in the form of thin film between the PDMS liner and active-winding and eventually enables thin film evaporation on the outer surface of the active-winding. Therefore, EC based thermal management eliminates contact resistance between the winding and slot-liner and, enhances the heat extraction from the winding without compromising the electro-magnetic performance. Two-way coupled electro-magnetic (EM) – computational fluid dynamics/heat transfer (CFD/HT) and EM - lumped parameter thermal network (LPTN) models have been developed to assess the electro-thermal performance of the EC under steady and transient conditions. Taking a case study of a 125 kW jacket cooled BMW i3 motor and dielectric coolant FC-84, EC is shown to be capable of handling a maximum steady state rms current density of 26 A/mm2 at a evaporative heat transfer coefficient of 5,000 W/m2.K, which is about 78.7% higher compared to the traditional jacket cooling (JC). In case of EC, a maximum steady state and transient rms current density of 30 A/mm2 (106.2% higher compared to the JC), and 40.8 A/mm2 have been realized by using high thermal conductivity epoxy (1.9 W/m.K) impregnation material. Thermal performance of the EC is also assessed and compared with JC over a dynamic drive cycle. Lastly, a motorette testing has been performed to demonstrate the applicability of the proposed EC method and to validate the developed modeling framework.
This paper addresses the behavior of a candle flame in a long-duration, quiescent microgravity environment both on the space Shuttle and the Mir Orbiting Station. On the Shuttle, the flames became dim blue after an initial transient where there was significant yellow (presumably soot) in the flame. The flame lifetimes were typically less than 60 seconds. The safety-mandated candlebox that contained the candle flame inhibited oxygen transport to the flame and thus limited the flame lifetime. The flames on the Mir were similar, except that the yellow luminosity persisted longer into the flame lifetime because of a higher initial oxygen concentration, The Mir flames burned for as long as 45 minutes. The difference in the flame lifetime between the Shuttle and Mir flames was primarily the redesigned candlebox that did not inhibit oxygen transport to the flame. In both environments, the flame intensity and the height-to-width ratio gradually decreased as the ambient oxygen content in the sealed chamber slowly decreased. Both sets of experiments showed spontaneous, axisymmetric flame oscillations just prior to extinction. The paper also presents a numerical model of a candle flame. The formulation is two-dimensional and time-dependent in the gas phase with constant specific heats, thermal conductivity and Lewis number (although different species can have different Lewis numbers), one-step finite-rate kinetics, and gas-phase radiative losses from CO2 and H2O. The treatment of the liquid/wick phase assumes that the, fuel evaporates from a constant diameter sphere connected to an inert cone. The model predicts a steady flame with a shape and size quantitatively similar to the Shuttle and Mir flames. The computation predicts that the flame size will increase slightly with increasing ambient oxygen mole fraction. The model also predicts pre-extinction flame oscillations if the rate of decrease in ambient oxygen is small enough, such as that which would occur for a flame burning in a sealed ambient.
The GLAS General Circulation Model (GCM) was applied to the four-month simulation of the thermodynamic part of the Parkinson-Washington sea ice model using atmospheric boundary conditions. The sea ice thickness and distribution were predicted for the Jan. 1-Apr. 30 period using the GCM-fields of solar and infrared radiation, specific humidity and air temperature at the surface, and snow accumulation; the sensible heat and evaporative surface fluxes were consistent with the ground temperatures produced by the ice model and the air temperatures determined by the atmospheric concept. It was concluded that the Parkinson-Washington sea ice model results in acceptable ice concentrations and thicknesses when used with GLAS GCM for the Jan.-Apr. period suggesting the feasibility of fully coupled ice-atmosphere simulations with these two approaches.
The dynamic portable life support system (PLSS) simulation software Virtual Space Suit (V-SUIT) has been under development at the Technische Universitat Munchen since 2011 as a spin-off from the Virtual Habitat (V-HAB) project. The MATLAB(trademark)-based V-SUIT simulates space suit portable life support systems and their interaction with a detailed and also dynamic human model, as well as the dynamic external environment of a space suit moving on a planetary surface. To demonstrate the feasibility of a large, system level simulation like V-SUIT, a model of NASA's PLSS 1.0 prototype was created. This prototype was run through an extensive series of tests in 2011. Since the test setup was heavily instrumented, it produced a wealth of data making it ideal for model validation. The implemented model includes all components of the PLSS in both the ventilation and thermal loops. The major components are modeled in greater detail, while smaller and ancillary components are low fidelity black box models. The major components include the Rapid Cycle Amine (RCA) CO2 removal system, the Primary and Secondary Oxygen Assembly (POS/SOA), the Pressure Garment System Volume Simulator (PGSVS), the Human Metabolic Simulator (HMS), the heat exchanger between the ventilation and thermal loops, the Space Suit Water Membrane Evaporator (SWME) and finally the Liquid Cooling Garment Simulator (LCGS). Using the created model, dynamic simulations were performed using same test points also used during PLSS 1.0 testing. The results of the simulation were then compared to the test data with special focus on absolute values during the steady state phases and dynamic behavior during the transition between test points. Quantified simulation results are presented that demonstrate which areas of the V-SUIT model are in need of further refinement and those that are sufficiently close to the test results. Finally, lessons learned from the modelling and validation process are given in combination with implications for the future development of other PLSS models in V-SUIT.
Production models were developed and the confirmation of each one had significant astrophysical impact. These include radical modifications of propagation models, cosmic ray antiprotons injection from neighboring domains of antimatter, p production by evaporating primordial black holes, and cosmic ray p's as annihilation products of supersymmetry particles that might make up the dark dynamical mass of the Galaxy. It is that p's originating from supersymmetric parents might have distinct spectral features that would survive solar modulation; in one model, higgsino annihilation proceeds through the bb quark-antiquark channel, producing a spectral bump at approx. 0.3 GeV in the p spectrum.
Purpose : Dry-eye disease, an increasingly prevalent ocular-surface disorder, significantly alters tear physiology. Knowing the essential physics of tear dynamics in healthy and dry eyes benefits both diagnosis and treatment of dry eye. We introduce a physiological-based model to describe tear dynamics during blinking. Materials and methods : Tears are compartmentalized over the ocular surface; the blink cycle is divided into three repeating phases. Conservation laws quantify the tear volume and tear osmolarity of each compartment during each blink phase. Lacrimal-supply and tear-evaporation rates are varied to reveal the dependence of tear dynamics on dry-eye conditions, specifically tear osmolarity, tear volume, tear-turnover rate (TTR), and osmotic water flow. Results : Predicted periodic-steady tear-meniscus osmolarity is 309 and 321 mOsM in normal and dry eyes, respectively. Tear osmolarity, volume, and TTR all match available clinical measurements. Osmotic water flow through the cornea and conjunctiva contribute 10 and 50% to the total tear supply in healthy and dry-eye conditions, respectively. TTR in aqueous-deficient dry eye (ADDE) is only half that in evaporative dry eye (EDE). Conclusions : The compartmental periodic-steady tear-dynamics model accurately predicts tear behavior in normal and dry eyes. Inclusion of osmotic water flow is crucial to match measured tear osmolarity. Tear-dynamics predictions corroborate the use of TTR as a clinical discriminator between ADDE and EDE. The proposed model is readily extended to predict the dynamics of aqueous solutes such as drugs or fluorescent tags.
Atmospheric water-vapor isotopes have been proposed as a potentially powerful constraint on convection, which plays a critical role in Earth's present and future climate. It is shown here, however, that the mean tropical profile of HDO in the free troposphere does not usefully constrain the mean convective entrainment rate or precipitation efficiency. This is demonstrated using a single-column analytical model of atmospheric water isotopes. The model has three parameters: the entrainment rate, the precipitation efficiency, and the distance that evaporating condensates fall. At a given relative humidity, the possible range of HDO is small: its range is comparable to both the measurement uncertainty in the mean tropical profile and the structural uncertainty of a single-column model. Therefore, the mean tropical HDO profile is unlikely to add information about convective processes in a bulk-plume framework that cannot already be learned from relative humidity alone.
Spectral lines formed at lower atmospheric layers show peculiar profiles at the "leading edge" of ribbons during solar flares. In particular, increased absorption of the BBSO/GST He I λ10830 line, as well as broad and centrally reversed profiles in the spectra of the Mg II and C II lines observed by the IRIS satellite, has been reported. In this work, we aim to understand the physical origin of such peculiar IRIS profiles, which seem to be common of many, if not all, flares. To achieve this, we quantify the spectral properties of the IRIS Mg II profiles at the ribbon leading edge during four large flares and perform a detailed comparison with a grid of radiative hydrodynamic models using the RADYN+FP code. We also studied their transition region (TR) counterparts, finding that these ribbon front locations are regions where TR emission and chromospheric evaporation are considerably weaker compared to other parts of the ribbons. Based on our comparison between the IRIS observations and modeling, our interpretation is that there are different heating regimes at play in the leading edge and the main bright part of the ribbons. More specifically, we suggest that bombardment of the chromosphere by more gradual and modest nonthermal electron energy fluxes can qualitatively explain the IRIS observations at the ribbon leading front, while stronger and more impulsive energy fluxes are required to drive chromospheric evaporation and more intense TR emission in the bright ribbon. Our results provide a possible physical origin for the peculiar behavior of the IRIS chromospheric lines in the ribbon leading edge and new constraints for the flare models.
Whilst the physics of both classical evaporation and supercritical fluid mixing are reasonably well characterized and understood in isolation, little is known about the transition from one to the other in the context of liquid fuel systems. The lack of experimental data for microscopic droplets at realistic operating conditions impedes the development of phenomenological and numerical models. To address this issue we performed systematic measurements using high-speed long-distance microscopy, for three single-component fuels (n-heptane, n-dodecane, n-hexadecane), into gas at elevated temperatures (700–1200 K) and pressures (2–11 MPa). We describe these high-speed visualizations and the time evolution of the transition from liquid droplet to fuel vapour at the microscopic level. The measurements show that the classical atomization and vaporisation processes do shift to one where surface tension forces diminish with increasing pressure and temperature, but the transition to diffusive mixing does not occur instantaneously when the fuel enters the chamber. Rather, subcritical liquid structures exhibit surface tension in the near-nozzle region and then, after time surrounded by the hot ambient gas and fuel vapour, undergo a transition to a dense miscible fluid. Although there was clear evidence of surface tension and primary atomization for n-dodecane and n-hexadecane for a period of time at all the above conditions, n-heptane appeared to produce a supercritical fluid from the nozzle outlet when injected at the most elevated conditions (1200 K, 10 MPa). This demonstrates that the time taken by a droplet to transition to diffusive mixing depends on the pressure and temperature of the gas surrounding the droplet as well as the fuel properties. We summarise our observations into a phenomenological model which describes the morphological evolution and transition of microscopic droplets from classical evaporation through a transitional mixing regime and towards diffusive mixing, as a function of operating conditions. We provide criteria for these regime transitions as reduced pressure–temperature correlations, revealing the conditions where transcritical mixing is important to diesel fuel spray mixing.