Search NASA⌕ Search

SEARCH · Search NASA

Results for “higher heating value”

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.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Enhanced Elliptic Grid Generation

An enhanced method of elliptic grid generation has been invented. Whereas prior methods require user input of certain grid parameters, this method provides for these parameters to be determined automatically. "Elliptic grid generation" signifies generation of generalized curvilinear coordinate grids through solution of elliptic partial differential equations (PDEs). Usually, such grids are fitted to bounding bodies and used in numerical solution of other PDEs like those of fluid flow, heat flow, and electromagnetics. Such a grid is smooth and has continuous first and second derivatives (and possibly also continuous higher-order derivatives), grid lines are appropriately stretched or clustered, and grid lines are orthogonal or nearly so over most of the grid domain. The source terms in the grid-generating PDEs (hereafter called "defining" PDEs) make it possible for the grid to satisfy requirements for clustering and orthogonality properties in the vicinity of specific surfaces in three dimensions or in the vicinity of specific lines in two dimensions. The grid parameters in question are decay parameters that appear in the source terms of the inhomogeneous defining PDEs. The decay parameters are characteristic lengths in exponential- decay factors that express how the influences of the boundaries decrease with distance from the boundaries. These terms govern the rates at which distance between adjacent grid lines change with distance from nearby boundaries. Heretofore, users have arbitrarily specified decay parameters. However, the characteristic lengths are coupled with the strengths of the source terms, such that arbitrary specification could lead to conflicts among parameter values. Moreover, the manual insertion of decay parameters is cumbersome for static grids and infeasible for dynamically changing grids. In the present method, manual insertion and user specification of decay parameters are neither required nor allowed. Instead, the decay parameters are determined automatically as part of the solution of the defining PDEs. Depending on the shape of the boundary segments and the physical nature of the problem to be solved on the grid, the solution of the defining PDEs may provide for rates of decay to vary along and among the boundary segments and may lend itself to interpretation in terms of one or more physical quantities associated with the problem.

Kaul, Upender K.↗
NASA NTRSDate not supplied

Additional Trials & Tribulations in Synthesizing a Sulfide Standard

Well-characterized sulfide reference materials that can serve as matrix-matched calibrants for in-situ trace element analyses via laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) remain elusive. Here, we describe the creation an in-house sulfide standard at NASA JSC that will be used to measure siderophile and chalcophile trace elements in low pressure experimental products, specifically pentlandite ([FeNi]9S8) and pyrrhotite ([FeNi]1-xS), as well as natural sulfides in terrestrial and meteoritic mineral assemblages. Recent methods in creating a homogenous sulfide standard include pressed pellets, synthesized chips, or fused glass [e.g. 1-3]. Other studies, however, have had success in forming homogenous FeS standards via synthesization methods [e.g. 4-6]. We present a variation of the method described in [4], in which we create a pyrrhotite standard doped with a variety of trace elements (Zn, As, Se, Mo, Ru, Rh, Pd, Sn, Sb, Te, W, Os, Ir, Pt, Au) at ~10-40 ppm and Cu at ~200 ppm. Sulfide compositions were constructed using Fe and Ni metal and elemental sulfur powders. Trace elements were added to elemental sulfur from Atomic Absorption Spectroscopy (AAS) elemental standards as either nitrate or chloride solutions to prevent oxidation of the metal powders. The dried sulfur mixture was mixed with Fe and Ni powders and mechanically mixed before two aliquots were placed in separate SiO2 tubes. Each tube was held under vacuum for ~30 minutes, sealed under vacuum, and then heated at 800C for 48 hours. Like [4], the synthesis products were composed of pourous sulfide crystals. Major element analyses of both experimental aliquots, collected using an Electron Probe Microanalyser (EPMA), yield indistinguishable major element compositions (uncertainties in 2SE & 2RSE[%]), with an average of 57.90  0.09 (0.16 %), 4.95  0.03 (0.55 %), and 39.70  0.14 (0.34 %), for Fe, Ni, and S, respectively. Trace element data were measured using a Photon Machines 193nm laser ablation system coupled to a Thermo-Scientific Element-XR ICP-MS. Spot sizes were limited to 50 μm due to the porous nature of the sulfide target material. Trace element abundances, normalized to Fe as an internal standard, were also found to be homogenous between the two aliquots, with weighted mean 2RSE (%) values of <3.0 for all trace elements. Synthesized products were re-powdered and absolute concentrations measured via solution ICP-MS. Although the sulfide appears to be homogenous, sintering experiments will be performed to more closely match the standard density to natural sulfides and minimize differences in ablation behavior. Improved density also allows for higher sensitivity (i.e. more compact target material) and larger spot sizes or traverses, as void space is eliminated. Additional major and trace element analyses on the products of the sintering experiments will be undertaken.

Jacob B Setera↗
NASA NTRS2010

Regolith Volatile Characterization (RVC) in RESOLVE

Resource investigation in the lunar poles is of importance to the potential impact of in-situ resource utilization (ISRU). The RESOLVE project developed a payload to investigate the permanently shadowed areas of the lunar poles and demonstrate ISRU technology. As a part of the RESOLVE project, the regolith volatile characterization (RVC) subsystem was designed to examine the release of volatiles from sample cores. The test sample was heated in the reactor to release the volatiles where they were analyzed with gas chromatography. Subsequently, the volatile sample was introduced into the lunar water resource demonstration (LWRD) subsystem where the released hydrogen and water were selectively captured. The objective of the Regolith Volatile Characterization (RVC) subsystem was to heat the crushed core sample and determine the desorption of volatile species of interest. The RVC subsystem encompasses the reactor and the system for volatile analysis. The system was designed to analyze H2, He, CO, CO2, N2, 02, CH4, H2S and H2O. The GC chosen for this work is a Siemens MicroSAM process GC with 3 columns and 8 TCD detectors. Neon was chosen as the carrier gas to enhance the analysis of hydrogen and helium.The limit of detection for the gases is approx.1000ppm for H2, CO. CO2 , N2, O2 and H2 S. The limit of detection for CH4 is approx.4000ppm and the water limit of detection is -10000 ppm with a sample analysis time of 2-3 minutes. These values (with the exception of water and H2S) were determined by dilution of a six gas mixture from Scott Gas (5% CO2, CO, O2, N2, 4% CH4 and H2) using mass flow controllers (MFC5). Water was calibrated at low levels using an in house relative humidity (RH) generator. H 2S and high concentrations of H2 were calibrated by diluting a pure stream of gas with MFCs. Higher concentrations of N2 and 02 were calibrated using Air again diluting with MFCs. There were three modification goals for the GC in EBU2 that would allow this process GC to be used in the field demo for RESOLVE. The first modification was to decrease the weight associated with the GC, this included eliminating the explosion proof case (Figure 1) and replacing it with a lightweight case as well as using an on board COPV tank for the neon carrier gas. The next goal was to add a second oven for the molecular sieve column to allow for dual temperature control during GC operation; the separation of hydrogen and helium is optimum at lower temperatures while the water analysis required higher temperatures creating a competing design requirement. The second oven also allows a lower limit of detection for water quantification and avoids the possibility of water condensing in the GC which could ruin the column characteristics. The final goal was to modify the column arrangement to optimize the system for our specific application. Figure 2 shows the internal details of the module optimized optimized for our field application. The modifications and performance of the gas analysis system will be discussed in detail.

Captain, Janine↗
NASA NTRS2014

A Thermal Analysis of a Hot-Wire Probe for Icing Applications

This paper presents a steady-state thermal model of a hot-wire instrument applicable to atmospheric measurement of water content in clouds. In this application, the power required to maintain the wire at a given temperature is used to deduce the water content of the cloud. The model considers electrical resistive heating, axial conduction, convection to the flow, radiation to the surroundings, as well as energy loss due to the heating, melting, and evaporation of impinging liquid and or ice. All of these parameters can be varied axially along the wire. The model further introduces a parameter called the evaporation potential which locally gauges the maximum fraction of incoming water that evaporates. The primary outputs of the model are the steady-state power required to maintain a spatially-average constant temperature as well as the variation of that temperature and other parameters along the wire. The model is used to understand the sensitivity of the hot-wire performance to various flow and boundary conditions including a detailed comparison of dry air and wet (i.e. cloud-on) conditions. The steady-state power values are compared to experimental results from a Science Engineering Associates (SEA) Multi-Element probe, a commonly used water-content measurement instrument. The model results show good agreement with experiment for both dry and cloud-on conditions with liquid water content. For ice, the experimental measurements under read the actual water content due to incomplete evaporation and splashing. Model results, which account for incomplete evaporation, are still higher than experimental results where the discrepancy is attributed to splashing mass-loss which is not accounted in the model.

Thermal Analysis↗
NASA NTRS2004

Global Dayside Ionospheric Uplift and Enhancement Associated with Interplanetary Electric Fields

The interplanetary shock/electric field event of 5-6 November 2001 is analyzed using ACE interplanetary data. The consequential ionospheric effects are studied using GPS receiver data from the CHAMP and SAC-C satellites and altimeter data from the TOPEX/ Poseidon satellite. Data from ~100 ground-based GPS receivers as well as Brazilian Digisonde and Pacific sector magnetometer data are also used. The dawn-to-dusk interplanetary electric field was initially ~33 mV/m just after the forward shock (IMF BZ = -48 nT) and later reached a peak value of ~54 mV/m 1 hour and 40 min later (BZ = -78 nT). The electric field was ~45 mV/m (BZ = -65 nT) 2 hours after the shock. This electric field generated a magnetic storm of intensity DST = -275 nT. The dayside satellite GPS receiver data plus ground-based GPS data indicate that the entire equatorial and midlatitude (up to +/-50(deg) magnetic latitude (MLAT)) dayside ionosphere was uplifted, significantly increasing the electron content (and densities) at altitudes greater than 430 km (CHAMP orbital altitude). This uplift peaked ~2 1/2 hours after the shock passage. The effect of the uplift on the ionospheric total electron content (TEC) lasted for 4 to 5 hours. Our hypothesis is that the interplanetary electric field ''promptly penetrated'' to the ionosphere, and the dayside plasma was convected (by E x B) to higher altitudes. Plasma upward transport/convergence led to a ~55-60% increase in equatorial ionospheric TEC to values above ~430 km (at 1930 LT). This transport/convergence plus photoionization of atmospheric neutrals at lower altitudes caused a 21% TEC increase in equatorial ionospheric TEC at ~1400 LT (from ground-based measurements). During the intense electric field interval, there was a sharp plasma ''shoulder'' detected at midlatitudes by the GPS receiver and altimeter satellites. This shoulder moves equatorward from -54(deg) to -37(deg) MLAT during the development of the main phase of the magnetic storm. We presume this to be an ionospheric signature of the plasmapause and its motion. The total TEC increase of this shoulder is ~80%. Part of this increase may be due to a "superfountain effect." The dayside ionospheric TEC above ~430 km decreased to values ~45% lower than quiet day values 7 to 9 hours after the beginning of the electric field event. The total equatorial ionospheric TEC decrease was ~16%. This decrease occurred both at midlatitudes and at the equator. We presume that thermospheric winds and neutral composition changes produced by the storm-time Joule heating, disturbance dynamo electric fields, and electric fields at auroral and subauroral latitudes are responsible for these decreases.

ionosphere↗
NASA NTRS2014

ADEPT for Secondary Payloads

Deployable entry vehicles provide an enabling capability for large payloads at destinations across the solar system. The primary advantage of such systems derives from the fundamental limitation of available launch vehicle shroads preventing tried-and-true rigid aeroshell technology from being applied. However, little mission-pull exists for deployable entry vehicle technology for missions where the entry vehicle is the primary payload and can fit within available shrouds. The propensity toward applying rigid aeroshell technology is logical given that deployable concepts are higher risk and provide no marginal scientific benefit for missions where a rigid aeroshell is an option. A movement toward flying smaller spacecraft has emerged in recent years and is revealing novel applications for deployable entry vehicles as secondary payloads. Secondary payloads must minimize interference with the primary payload in order to ease accomodation. Advancements in small spacecraft technology have caught the attention of government entities who are now recognizing the high value proposition of small spacecraft for frequent, incremental technology development and science return rather than the traditional high cost and infrequent big-bang approaches. Worldwide, well over 200 small satellites (mass < 180 kg) have launched in recent years for such purposes as eduation, technology development, science, commerce, and defense. The Adaptable Deployable Entry and Placement Technology (ADEPT), a mechanically-deployable entry vehicle technology, has been under development at NASA since 2011. An Ames Center Innovation Fund (CIF) proposal was recently awarded to study ADEPT for EDL of small payloads. The unique capability of ADEPT for small science payloads comes from its ability to stow within a slender volume and deploy passively to achieve a mass-efficient drag surface with a high heat rate capability. The low ballistic coefficient results in entry heating and mechanical loads that can be met by a revolutionary three-dimensionally woven carbon fabric supported with rib structures. This carbon fabric has test-proven capability as both primary structure and payload thermal protection. This presentation will inform the EDL community of the capability of ADEPT for delivery of small secondary payloads as well as provide status of on-going technology development activities. Results of a mission applications study for secondary payloads to Mars and Venus will be presented. A primary goal of these mission design activities is to integrate high-value verification approaches wherever possible. For example, for some mission applications it is possible to test in NASA arc jet facilities at full-scale. This enables endto- end thermostructural requirement verification with a single test. Another important challenge that needs to be solved is the reliability of the decelerator to achieve the desired shape using simple passive mechanical actuators (such as springs) that do not require power from the primary payload for deployment. A configuration of ADEPT is being prototyped and used to evaluate the functionality and reliability of the deployment mechanism through bench-top testing. The prototype will then be tested in a low speed wind tunnel in late 2014 to evaluate aerodynamic load-carrying capability and deflected shape. Results of the deployment testing and status of wind tunnel test plans will be discussed.

Smith, B.↗
NASA NTRS2008

Wide-Band Radar for Measuring Thickness of Sea Ice

A wide-band penetrating radar system for measuring the thickness of sea ice is under development. The need for this or a similar system arises as follows: Spatial and temporal variations in the thickness of sea ice are important indicators of heat fluxes between the ocean and atmosphere and, hence, are important indicators of climate change in polar regions. A remote-sensing system that could directly measure the thickness of sea ice over a wide thickness range from aboard an aircraft or satellite would be of great scientific value. Obtaining thickness measurements over a wide region at weekly or monthly time intervals would contribute significantly to understanding of changes in the spatial distribution and of the mass balance of sea ice. A prototype of the system was designed on the basis of computational simulations directed toward understanding what signal frequencies are needed to satisfy partly competing requirements to detect both bottom and top ice surfaces, obtain adequate penetration despite high attenuation in the lossy sea-ice medium, and obtain adequate resolution, all over a wide thickness range. The prototype of the system is of the frequency-modulation, continuous-wave (FM-CW) type. At a given time, the prototype functions in either of two frequency-band/operational-mode combinations that correspond to two thickness ranges: a lower-frequency (50 to 250 MHz) mode for measuring thickness greater than about 1 m, and a higher frequency (300 to 1,300 MHz) mode for measuring thickness less than about 1 m. The bandwidth in the higher-frequency (lesser-thickness) mode is adequate for a thickness resolution of 15 cm; the bandwidth in the lower-frequency (greater-thickness) mode is adequate for a thickness resolution of 75 cm. Although a thickness resolution of no more than 25 cm is desired for scientific purposes, the 75-cm resolution was deemed acceptable for the purpose of demonstrating feasibility. The prototype was constructed as a modified version of a 500-to-2,000-MHz FM-CW radar system developed previously for mapping near-surface internal layers of the Greenland ice sheet. The prototype included two sets of antennas: one for each frequency-band/mode. For Arctic and Antarctic field tests, the prototype was mounted on a sled that was towed across the ice. The Arctic field test was performed in the lower-frequency mode on ice ranging in thickness from 1 to 4 m. In the analysis of the results of the Arctic field test, a comparison of the radar-determined ice thicknesses with actual ice thicknesses yielded an overall mean difference of 14 cm and standard deviation of 30 cm. The Antarctic field test was performed in the higher-frequency mode; analysis of the results led to the conclusion that this mode is useful for measuring thicknesses between 0.5 and 1 m. Several modifications have been conceived for implementation in further development toward an improved practical system: The system would function in a single frequency-band/mode (100 to 1,200 MHz) that would afford a resolution of about 15 cm. There would be a single antenna system that would be optimized for the entire 100-to-1,000-MHz frequency band. To enable ice-thickness surveys over larger areas, the system would be made capable of operating aboard a low-flying aircraft that could be either piloted or robotic. Data-processing techniques to deconvolve the system response have been developed on the basis of impulse-response measurements over a calm ocean. Implementation of these techniques in the system would enable correction for imperfections of the system and would thereby increase the effective sensitivity of the system.

Gogineni, Prasad↗
NASA NTRS2023

Phyllosilicate Decomposition on Bennu Due to Prolonged Surface Exposure

The Origins, Spectral Interpretation, Resource Identification, and Security–Regolith Explorer (OSIRIS-REx) mission to carbonaceous asteroid (101955) Bennu performed detailed mapping with a suite of instruments to characterize the composition and geology of its surface. Here we use data from the OSIRIS-REx Thermal Emission Spectrometer (OTES) instrument to investigate the relationship of OTES-derived spectral indices to other derived data products from OTES, the OSIRIS-REx Camera Suite (OCAMS), and the OSIRIS-REx Visible and InfraRed Spectrometer (OVIRS) at global and local scales. We quantitatively confirm that high values of the OTES silicate stretching slope (from ∼10 to ∼12 μm) in midday spectra that are indicative of thin and/or patchy dust cover are strongly associated with low thermal inertia (high porosity), low albedo boulders on Bennu. These high porosity boulders have brecciated textures with embedded clasts that likely originated on Bennu's parent body or during its disruption. The high porosity of these boulders is a key factor in the local production of the dust or its entrapment, as some large, brecciated boulders with a lower porosity have little evidence of dust. A second OTES spectral parameter, the silicate bending band depth near 22.7 μm applied to early evening spectra, is not correlated to thermal inertia, but is weakly to strongly correlated to albedo, OVIRS-derived 1.05 μm and 2.74 μm band depths, OVIRS-derived hydrogen abundance, and modeled nanophase troilite abundance. In several regions on Bennu there is a strong spatial relationship between these parameters, whereby areas with shallower silicate bending bands also have shallower 1.05 μm and 2.74 μm bands and lower albedo with higher nanophase troilite abundances. These correlations, combined with analysis of the silicate bending band in laboratory experiments of space weathered and mildly heated carbonaceous chondrites, suggests that decreased silicate bending band depths signify decomposition of phyllosilicates, likely Fe-bearing, due to space weathering or mild heating (<600 °C) via solar radiation during Bennu's time in near-Earth space. There is a strong association of larger silicate bending band depths in areas dominated by small rocks and unresolved material and in areas with small (≤ 25 m) craters identified as the spectrally reddest on Bennu, suggesting that this material has been more recently exposed due to impact and/or mass wasting processes. The shallowest silicate bending depths are associated with larger rocks and boulders that appear to have the longest surface exposure history, although there is band depth variation among them suggesting either initial composition variation that resulted in different responses to space weathering or heating, or varied exposure history of individual boulders themselves. We predict that any grains returned from Bennu with a significant surface exposure history will be characterized by shallower 22.7 μm, 1.05 μm and 2.7 μm band depths and increased sulfide (troilite) abundance, as well as textural and chemical evidence for phyllosilicate dehydration.

Asteroids, surfaces↗
NASA NTRS2001

Investigation of the Polytropic Relationship Between Density and Temperature Within Interplanetary Coronal Mass Ejections Using Numerical Simulations

Single-point spacecraft measurements within coronal mass ejections (CMEs) often exhibit a negative correlation between electron density and temperature. At least two opposing interpretations have been suggested for this relationship. If, on one hand, these single spacecraft observations provide direct measures of the polytropic properties of the plasma, then they imply that the polytropic index for the electrons gamma(sub e) is often < 1. Moreover, since the electrons carry the bulk of the pressure (via their significantly higher temperature), this further implies that the dynamics of CME evolution are dominated by an effective polytropic index gamma(sub e)ff < 1. On the other hand, gamma < 1 implies that as the ejecta propagate away from the Sun and expand, they also heat up; a result clearly at odds with in situ observations. In contrast to these CME intervals, many studies have shown that the quiescent solar wind exhibits a positive correlation between electron density and temperature, suggesting that gamma(sub e) > 1. In this study we simulate the evolution of a variety of CME-like disturbances in the solar wind using a one-dimensional, single-fluid model, to address the interpretation of the relationship between electron density and temperature within CMEs at fixed locations in space. Although we strictly impose a polytropic relationship (with gamma = constant) throughout our simulations, we demonstrate that a variety of correlations can exist between density and temperature at fixed points. Furthermore, we demonstrate that the presence of only local uncorrelated random fluctuations in density and temperature can produce a negative correlation. Consequently, we conclude that these single-point observations of negative correlations between electron density and temperature cannot be used to infer the value of gamma(sub e). Instead, we suggest that entropy variations, together with the plasma's tendency to achieve pressure balance with its surroundings, are responsible for the observed profiles.

Riley, Pete↗
NASA NTRS2010

Metal/Silicate Partitioning at High Pressures and Temperatures

The behavior of siderophile elements during metal-silicate segregation, and their resulting distributions provide insight into core formation processes. Determination of partition coefficients allows the calculation of element distributions that can be compared to established values of element abundances in the silicate (mantle) and metallic (core) portions of the Earth. Moderately siderophile elements, including W, are particularly useful in constraining core formation conditions because they are sensitive to variations in T, P, oxygen fugacity (fO2), and silicate composition. To constrain the effect of pressure on W metal/silicate partitioning, we performed experiments at high pressures and temperatures using a multi anvil press (MAP) at NASA Johnson Space Center and laser-heated diamond anvil cells (LHDAC) at the University of Maryland. Starting materials consisted of natural peridotite mixed with Fe and W metals. Pressure conditions in the MAP experiments ranged from 10 to 16 GPa at 2400 K. Pressures in the LHDAC experiments ranged from 26 to 58 GPa, and peak temperatures ranged up to 5000 K. LHDAC experimental run products were sectioned by focused ion beam (FIB) at NASA JSC. Run products were analyzed by electron microprobe using wavelength dispersive spectroscopy. Liquid metal/liquid silicate partition coefficients for W were calculated from element abundances determined by microprobe analyses, and corrected to a common fO2 condition of IW-2 assuming +4 valence for W. Within analytical uncertainties, W partitioning shows a flat trend with increasing pressure from 10 to 16 GPa. At higher pressures, W becomes more siderophile, with an increase in partition coefficient of approximately 0.5 log units.

Shofner, G.↗
NASA NTRS1993

A model of concurrent flow flame spread over a thin solid fuel

A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.

Ferkul, Paul V.↗
NASA NTRS1989

Nonlinear interactions in mixing layers and compressible heated round jets

The nonlinear interactions between a fundamental instability mode and both its harmonics and the changing mean flow are studied using the weakly nonlinear stability theory of Stuart and Watson, and numerical solutions of coupled nonlinear partial differential equations. The first part focuses on incompressible cold (or isothermal; constant temperature throughout) mixing layers, and for these, the first and second Landau constants are calculated as functions of wavenumber and Reynolds number. It is found that the dominant contribution to the Landau constants arises from the mean flow changes and not from the higher harmonics. In order to establish the range of validity of the weakly nonlinear theory, the weakly nonlinear and numerical solutions are compared and the limitation of each is discussed. At small amplitudes and at low-to-moderate Reynolds numbers, the two results compare well in describing the saturation of the fundamental, the distortion of the mean flow, and the initial stages of vorticity roll-up. At larger amplitudes, the interaction between the fundamental, second harmonic, and the mean flow is strongly nonlinear and the numerical solution predicts flow oscillations, whereas the weakly nonlinear theory yields saturation. In the second part, the weakly nonlinear theory is extended to heated (or nonisothermal; mean temperature distribution) subsonic round jets where quadratic and cubic nonlinear interactions are present, and the Landau constants also depend on jet temperature ratio, Mach number and azimuthal mode number. Under exponential growth and nonlinear saturation, it is found that heating and compressibility suppress the growth of instability waves, that the first azimuthal mode is the dominant instability mode, and that the weakly nonlinear solution describes the early stages of the roll-up of an axisymmetric shear layer. The receptivity of a typical jet flow to pulse type input disturbance is also studied by solving the initial value problem and then examining the behavior of the long-time solution.

Jarrah, Yousef Mohd↗
NASA NTRS2015

Cost Scaling of a Real-World Exhaust Waste Heat Recovery Thermoelectric Generator: A Deeper Dive

Cost is equally important to power density or efficiency for the adoption of waste heat recovery thermoelectric generators (TEG) in many transportation and industrial energy recovery applications. In many cases the system design that minimizes cost (e.g., the $/W value) can be very different than the design that maximizes the system's efficiency or power density, and it is important to understand the relationship between those designs to optimize TEG performance-cost compromises. Expanding on recent cost analysis work and using more detailed system modeling, an enhanced cost scaling analysis of a waste heat recovery thermoelectric generator with more detailed, coupled treatment of the heat exchangers has been performed. In this analysis, the effect of the heat lost to the environment and updated relationships between the hot-side and cold-side conductances that maximize power output are considered. This coupled thermal and thermoelectric treatment of the exhaust waste heat recovery thermoelectric generator yields modified cost scaling and design optimization equations, which are now strongly dependent on the heat leakage fraction, exhaust mass flow rate, and heat exchanger effectiveness. This work shows that heat exchanger costs most often dominate the overall TE system costs, that it is extremely difficult to escape this regime, and in order to achieve TE system costs of $1/W it is necessary to achieve heat exchanger costs of $1/(W/K). Minimum TE system costs per watt generally coincide with maximum power points, but Preferred TE Design Regimes are identified where there is little cost penalty for moving into regions of higher efficiency and slightly lower power outputs. These regimes are closely tied to previously-identified low cost design regimes. This work shows that the optimum fill factor Fopt minimizing system costs decreases as heat losses increase, and increases as exhaust mass flow rate and heat exchanger effectiveness increase. These findings have profound implications on the design and operation of various thermoelectric (TE) waste heat 3 recovery systems. This work highlights the importance of heat exchanger costs on the overall TEG system costs, quantifies the possible TEG performance-cost domain space based on heat exchanger effects, and provides a focus for future system research and development efforts.

cost analysis↗
NASA NTRS2012

Assessment and Enhancement of MERRA Land Surface Hydrology Estimates

The Modern-Era Retrospective analysis for Research and Applications (MERRA) is a state-ofthe-art reanalysis that provides, in addition to atmospheric fields, global estimates of soil moisture, latent heat flux, snow, and runoff for 1979-present. This study introduces a supplemental and improved set of land surface hydrological fields ("MERRA-Land") generated by re-running a revised version of the land component of the MERRA system. Specifically, the MERRA-Land estimates benefit from corrections to the precipitation forcing with the Global Precipitation Climatology Project pentad product (version 2.1) and from revised parameter values in the rainfall interception model, changes that effectively correct for known limitations in the MERRA surface meteorological forcings. The skill (defined as the correlation coefficient of the anomaly time series) in land surface hydrological fields from MERRA and MERRA-Land is assessed here against observations and compared to the skill of the state-of-the-art ERA-Interim (ERA-I) reanalysis. MERRA-Land and ERA-I root zone soil moisture skills (against in situ observations at 85 US stations) are comparable and significantly greater than that of MERRA. Throughout the northern hemisphere, MERRA and MERRA-Land agree reasonably well with in situ snow depth measurements (from 583 stations) and with snow water equivalent from an independent analysis. Runoff skill (against naturalized stream flow observations from 18 US basins) of MERRA and MERRA-Land is typically higher than that of ERA-I. With a few exceptions, the MERRA-Land data appear more accurate than the original MERRA estimates and are thus recommended for those interested in using MERRA output for land surface hydrological studies.

Reichle, Rolf H.↗
NASA NTRS2021

Methodology to Assess the Changing Risk of Yield Failure Due to Heat and Drought Stress Under Climate Change

While the understanding of average impacts of climate change on crop yields is improving, few assessments have quantified expected impacts on yield distributions and the risk of yield failures. Here we present the relative distribution as a method to assess how the risk of yield failure due to heat and drought stress (measured in terms of return period between yields falling 15% below previous five year Olympic average yield) responds to changes of the underlying yield distributions under climate change. Relative distributions are used to capture differences in the entire yield distribution between baseline and climate change scenarios, and to further decompose them into changes in the location and shape of the distribution. The methodology is applied here for the case of rainfed wheat and grain maize across Europe using an ensemble of crop models under three climate change scenarios with simulations conducted at 25 km resolution. Under climate change, maize generally displayed shorter return periods of yield failures (with changes under RCP 4.5 between −0.3 and 0 years compared to the baseline scenario) associated with a shift of the yield distribution towards lower values and changes in shape of the distribution that further reduced the frequency of high yields. This response was prominent in the areas characterized in the baseline scenario by high yields and relatively long return periods of failure. Conversely, for wheat, yield failures were projected to become less frequent under future scenarios (with changes in the return period of −0.1 to +0.4 years under RCP 4.5) and were associated with a shift of the distribution towards higher values and a change in shape increasing the frequency of extreme yields at both ends. Our study offers an approach to quantify the changes in yield distributions that drive crop yield failures. Actual risk assessments additionally require models that capture the variety of drivers determining crop yield variability and scenario climate input data that samples the range of probable climate variation.

Climate risk assessment↗
NASA NTRS2002

Characteristics of Pool Boiling on Graphite-Copper Composite Surfaces

Nucleate pool boiling performance of different liquids on graphite-copper composite (Gr-Cu) surfaces has been experimentally studied and modeled. Both highly wetting fluids, such as freon-113 and pentane, and a moderately wetting fluid (water) were tested on the Gr-Cu surfaces with different graphite-fiber volume fractions to reveal the enhancement effects of the composite surfaces on the nucleate pool boiling. Results of the experiments show that the graphite-fiber volume fraction has an optimum value. The Gr-Cu composite surface with 25 percent graphite-fiber volume (f=0.25) has a maximum enhancement effect on the nucleate boiling heat transfer comparing to the pure copper surface. For the highly wetting fluid, the nucleate boiling heat transfer is generally enhanced on the Gr- Cu composite surfaces by 3 to 6 times shown. In the low heat flux region, the enhancement is over 6 times, but in the high heat flux region, the enhancement is reduced to about 40%. For the moderately wetting fluid (water), stronger enhancement of nucleate boiling heat transfer is achieved on the composite surface. It shown the experimental results in which one observes the nucleate boiling heat transfer enhancement of 5 to 10 times in the low heat flux region and an enhancement of 3 to 5 times in the high heat flux region. Photographs of bubble departure during the initial stage of nucleate boiling indicate that the bubbles detached from the composite surface are much smaller in diameter than those detached from the pure copper surface. Typical photographs are presented.It shows that the bubbles departed from the composite surface have diameters of only O(0.1) mm, while those departed from the pure copper surface have diameters of O(1) mm. It is also found that the bubbles depart from the composite surface at a much higher frequency, thus forming vapor columns. These two phenomena combined with high thermal conductivity of the graphite fiber are considered the mechanisms for such a significant augmentation in nucleate boiling heat transfer on the composite surfaces. A physical model is developed to describe the phenomenon of bubble departure from the composite surface: The preferred site of bubble nucleation is the fiber tip because of higher tip temperature than the surrounding copper base and poor wettability of the graphite tip compared with that of the base material (copper). The high evaporation rate near the contact line produces the vapor cutback due to the vapor recoil pushing the three-phase line outwards from the fiber tip, and so a neck of the bubble is formed near the bubble bottom. Evaporation and surface tension accelerate the necking process and finally result in the bubble departure while a new small bubble is formed at the tip when the surface tension pushes the three-phase line back to the tip. The process is schematically shown. The proposed model is based on and confirmed by experimental results.

Zhang, Nengli↗
NASA NTRS1994

A hydrodynamic approach to cosmology: The mixed dark matter cosmological scenario

We compute the evolution of spatially flat, mixed cold and hot dark matter models containing both baryonic matter and two kinds of dark matter. Hydrodynamics is treated with a highly developed Eulerian hydrodynamic code (see Cen 1992). A standard particle-mesh (PM) code is also used in parallel to calculate the motion of the dark matter components. We adopt the following parameters: h equivalent to (sub 0)/100 km/s Mpc(exp -1) = 0.5, OMEGA(sub C) = 0.3, and OMEGA(sub B) = 0.06, with amplitude of the perturbation spectrum fixed by the Cosmic Background Explorer Satellite (COBE) Dark Matter Radiation (DMR) measurements (Smoot et al. 1992) being sigma (sub 8) = 0.67. Four different boxes are simulated with box sizes of L = (64, 16, 4, 1) h(exp -1) Mpc, respectively, the two small boxes providing good resolution but little valid information due to absence of large-scale power. We use 128(exp 3) approximate 10(exp 6.3) baryonic cells, 128(exp .3) cold dark matter particles, and 2 x 128(exp 3) hot dark matter particles. In addition to the dark matter we follow separately six baryonic species (H, H(+), He, He(+), He(++), e(-)) with allowance for both (nonequilibrium) collisional and radiative ionization in every cell. The background radiation field is also followed in detail with allowance made for both continuum and line processes, to allow nonequilibrium heating and cooling processes to be followed in detail. The mean final Zeldovich-Sunyaev y parameter is estimated to be y Bar = (5.4 + or - 2.7) x 10(exp -7) below currently attainable observations, with a rms fluctuation of approximately delta bar y = (0.6 + or - 3.0) x 10(exp -7) on arcminute scales. The rate of galaxy formation peaks at an even later epoch (z approximate 0.3) than in the standard (OMEGA = 1, sigma sub 8 = 0.67) cold dark matter (CDM) model (z approximate 0.5) and, at a redshift of z = 4, is nearly a factor of 100 lower than for the CDM model with the same value of sigma sub 8. With regard to mass function, the smallest objects are stabilized against collapse by thermal energy: the mass-weighted mass spectrum has a broad peak in the vicinity of M(sub B) = 10(exp 9.5) solar mass with a reasonable fit to the Schechter luminosity function if the ratio of baryon mass to blue light is approximately 4. In addition, one very large PM simulation was made in a box with size (320 h(exp - 1) Mpc) containing 3 x 200(exp 3) = 10(exp 7.4) particles. Utilizing this simulation we find that the model yields a cluster mass function which is about a factor of 4 higher than observed, but a cluster-cluster correlation length marginally lower than observed, but that both are closer to observations than in the (COBE) normalized CDM model. The one-dimensional pairwise velocity dispersion is 605 + or - 8 km/s at 1/h separation, lower than that of the DCM model normalized to COBE, but still significant higher than observations (Davis & Peebles 1983). A plausible velocity bias b(sub v) = 0.8 + or - 0.1 on this scale will reduce but not remove the discrepancy. The velocity auto-correlat ion function has a coherence length of 40/h Mpc, which is somewhat lower than the observed counterpart. In all these respects the model would be improved by decreasing the cold fraction of the dark OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDB). But formation of galaxies and clusters of galaxies is much later in this model than in COBE-normalized CDM, perhaps too late. To improve on these constraints a larger ratio of OMEGA(sub CDM)/ (OMEGA(sub CDM) + OMEGA(sub HDM)) is required than the value of 0.67 adopted here. It does not seem possible to find a value for this ratio which would satisfy all tests. Overall, the model is similar both on large and intermediate scales to the standard CDM model normalized to the same value of sigma(sub B), but the problem with regard to late formation of galaxies is more severe in this model than in that CDM model. Adding hot dark matter, significantly improves the ability of the COBE-normalized CDM scenario to fit existing observations, but the model is in fact not as good as the CDM model with the same sigma(sub 8) and is still probably unsatisfactory with regard to several critical tests.

Cen, Renyue↗
NASA NTRS1994

Soft X-ray spectral fits of Geminga with model neutron star atmospheres

The spectrum of the soft X-ray pulsar Geminga consists of two components, a softer one which can be interpreted as thermal-like radiation from the surface of the neutron star, and a harder one interpreted as radiation from a polar cap heated by relativistic particles. We have fitted the soft spectrum using a detailed magnetized hydrogen atmosphere model. The fitting parameters are the hydrogen column density, the effective temperature T(sub eff), the gravitational redshift z, and the distance to radius ratio, for different values of the magnetic field B. The best fits for this model are obtained when B less than or approximately 1 x 10(exp 12) G and z lies on the upper boundary of the explored range (z = 0.45). The values of T(sub eff) approximately = (2-3) x 10(exp 5) K are a factor of 2-3 times lower than the value of T(sub eff) obtained for blackbody fits with the same z. The lower T(sub eff) increases the compatibility with some proposed schemes for fast neutrino cooling of neutron stars (NSs) by the direct Urca process or by exotic matter, but conventional cooling cannot be excluded. The hydrogen atmosphere fits also imply a smaller distance to Geminga than that inferred from a blackbody fit. An accurate evaluation of the distance would require a better knowledge of the ROSAT Position Sensitive Proportional Counter (PSPC) response to the low-energy region of the incident spectrum. Our modeling of the soft component with a cooler magnetized atmosphere also implies that the hard-component fit requires a characteristic temperature which is higher (by a factor of approximately 2-3) and a surface area which is smaller (by a factor of 10(exp 3), compared to previous blackbody fits.

Meyer, R. D.↗