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.
The numerical operations involved in a currently used optimization technique are discussed and analyzed with special attention to the numerical accuracy. Alternative methods for deriving linear system transfer functions, finding the relationships between the transfer function coefficients and the design parameters, and solving a matrix equation are presented for more accurate and cost effective solutions.
Ocean surface evaporation and the latent heat it carries are the major components of the hydrologic and thermal forcing on the global oceans. However, there is practically no direct in situ measurements. Evaporation estimated from bulk parameterization methods depends on the quality and distribution of volunteer-ship reports which are far less than satisfactory. The only way to monitor evaporation with sufficient temporal and spatial resolutions to study global environment changes is by spaceborne sensors. The estimation of seasonal-to-interannual variation of ocean evaporation, using spacebased measurements of wind speed, sea surface temperature (SST), and integrated water vapor, through bulk parameterization method,s was achieved with reasonable success over most of the global ocean, in the past decade. Because all the three geophysical parameters can be retrieved from the radiance at the frequencies measured by the Scanning Multichannel Microwave Radiometer (SMMR) on Nimbus-7, the feasibility of retrieving evaporation directly from the measured radiance was suggested and demonstrated using coincident brightness temperatures observed by SMMR and latent heat flux computed from ship data, in the monthly time scale. However, the operational microwave radiometers that followed SMMR, the Special Sensor Microwave/Imager (SSM/I), lack the low frequency channels which are sensitive to SST. This low frequency channels are again included in the microwave imager (TMI) of the recently launched Tropical Rain Measuring Mission (TRMM). The radiance at the frequencies observed by both TMI and SSM/I were simulated through an atmospheric radiative transfer model using ocean surface parameters and atmospheric temperature and humidity profiles produced by the reanalysis of the European Center for Medium Range Weather Forecast (ECMWF). From the same ECMWF data set, coincident evaporation is computed using a surface layer turbulent transfer model. The sensitivity of the radiance to evaporation over various seasons and geographic locations are examined. The microwave frequencies with radiance that are significant correlated with evaporation are identify and capability of estimating evaporation directly from TMI will be discussed.
A measurement-based radiance estimation approach for vicarious radiometric calibration of spaceborne multispectral remote sensing systems has been developed. This simplified process eliminates the use of radiative transfer codes and reduces the number of atmospheric assumptions required to perform sensor calibrations. Like prior approaches, the simplified method involves the collection of ground truth data coincident with the overpass of the remote sensing system being calibrated, but this approach differs from the prior techniques in both the nature of the data collected and the manner in which the data are processed. In traditional vicarious radiometric calibration, ground truth data are gathered using ground-viewing spectroradiometers and one or more sun photometer( s), among other instruments, located at a ground target area. The measured data from the ground-based instruments are used in radiative transfer models to estimate the top-of-atmosphere (TOA) target radiances at the time of satellite overpass. These TOA radiances are compared with the satellite sensor readings to radiometrically calibrate the sensor. Traditional vicarious radiometric calibration methods require that an atmospheric model be defined such that the ground-based observations of solar transmission and diffuse-to-global ratios are in close agreement with the radiative transfer code estimation of these parameters. This process is labor-intensive and complex, and can be prone to errors. The errors can be compounded because of approximations in the model and inaccurate assumptions about the radiative coupling between the atmosphere and the terrain. The errors can increase the uncertainty of the TOA radiance estimates used to perform the radiometric calibration. In comparison, the simplified approach does not use atmospheric radiative transfer models and involves fewer assumptions concerning the radiative transfer properties of the atmosphere. This new technique uses two neighboring uniform ground target areas having different reflectance values. The target areas can be natural or artificial and must be large enough to minimize adjacent-pixel contamination effects. The radiative coupling between the atmosphere and the terrain needs to be approximately the same for the two targets. This condition can be met for relatively uniform backgrounds when the distance between the targets is within a few hundred meters. For each target area, the radiance leaving the ground in the direction of the satellite is measured with a radiometrically calibrated spectroradiometer. Using the radiance measurements from the two targets, atmospheric adjacency and atmospheric scattering effects can be subtracted, thereby eliminating many assumptions about the atmosphere and the radiative interaction between the atmosphere and the terrain. In addition, the radiometrically calibrated spectroradiometer can be used with a known reflectance target to estimate atmospheric transmission and diffuse- to-global ratios without the need for ancillary sun photometers. Several comparisons between the simplified method and traditional techniques were found to agree within a few percent. Hence, the simplified method reduces the overall complexity of performing vicarious calibrations and can serve as a method for validating traditional radiative transfer models
Since early 2000, NASA's Multi-angle Imaging SpectroRadiometer (MISR) instrument has been performing remote sensing retrievals of aerosol optical properties from the polar-orbiting Terra spacecraft. A noteworthy aspect of MISR observations over the ocean is that, for much of the Earth, some of the multi-angle views have contributions from solar reflection by the ocean surface (glint, or glitter), while others do not. Aerosol retrieval algorithms often discard these glint-influenced observations because they can overwhelm the signal and are difficult to predict without knowledge of the (wind-speed-driven) ocean surface roughness. However, theoretical studies have shown that multi-angle observations of a location at geometries with and without reflected sun glint can be a rich source of information, sufficient to support simultaneous retrieval of both the aerosol state and the wind speed at the ocean surface. We are in the early stages of creating such an algorithm. In this paper, we describe our assessment of the appropriate level of parameterization for simultaneous aerosol and ocean surface property retrievals using sun glint. For this purpose, we use generalized nonlinear retrieval analysis (GENRA), an information content assessment (ICA) technique employing Bayesian inference, and simulations from the Ahmad–Fraser iterative radiative transfer code. We find that four parameters are suitable: aerosol optical depth (τ), particle size distribution (expressed as the fine mode fraction f of small particles in a bimodal size distribution), surface wind speed (w), and relative humidity (r, to define the aerosol water content and complex refractive index). None of these parameters define ocean optical properties, as we found that the aerosol state could be retrieved with the nine MISR near-infrared views alone, where the ocean body is strongly absorbing in the open ocean. We also found that retrieval capability varies with observation geometry and that as τ increases so does the ability to determine aerosol intensive optical properties (r and f, while it decreases for w). Increases in w decrease the ability to determine the true value of that parameter but have minimal impact on retrieval of aerosol properties. We explored the benefit of excluding the two most extreme MISR view angles for which radiative transfer with the plane-parallel approximation is less certain, but we found no advantage in doing so. Finally, the impact of treating wind speed as a scalar parameter, rather than as a two-parameter directional wind, was tested. While the simpler scalar model does contribute to overall aerosol uncertainty, it is not sufficiently large to justify the addition of another dimension to parameter space. An algorithm designed upon these principles is in development. It will be used to perform an atmospheric correction with MISR for coincident ocean color (OC) observations by the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument, also on the NASA Terra spacecraft. Unlike MISR, MODIS is a single-view-angle instrument, but it has a more complete set of spectral channels ideal for determination of optical ocean properties. The atmospheric correction of MODIS OC data can therefore benefit from MISR aerosol retrievals. Furthermore, higher-spatial-resolution data from coincident MISR observations may also improve glint screening.
A physics-based Nuclear Thermal Propulsion (NTP) Testing Reference Design (TRD) power balance model was coded in Simulink to investigate engine performance for various design and parameter modifications. Since the primary mode of heat transfer in NTP engines is convective, the convective heat transfer coefficient (HTC) is a key parameter that requires accurate representation. The industry standard Westinghouse correlation has an uncertainty of ±20% which was investigated in this study. The results showed that a 20% decrease in the HTC led to a 4.14% increase in maximum fuel temperature while a 20% decrease in the HTC led to a 1.81% decrease in maximum fuel temperature suggesting that narrowing the uncertainty of this correlation through experimental work would be a critical step in the development of NTP engines. Furthermore, a maximum fuel temperature relationship with specific impulse was developed for the TRD engine which showed potential engine operation between specific impulse values of 715 and 900 seconds with minimal changes to the engine design. This graph could be useful for high level vehicle performance estimations for fuel types with different maximum operating temperatures.
Results are shown for a three-dimensional Navier-Stokes analysis of both the flow and the surface heat transfer for turbine applications. Heat transfer comparisons are made with the experimental shock-tunnel data of Dunn and Kim, and with the data of Blair for the rotor of the large scale rotating turbine. The analysis was done using the steady-state, three-dimensional, thin-layer Navier-Stokes code developed by Chima, which uses a multistage Runge-Kutta scheme with implicit residual smoothing. An algebraic mixing length turbulence model is used to calculate turbulent eddy viscosity. The variation in heat transfer due to variations in grid parameters is examined. The effects of rotation, tip clearance, and inlet boundary layer thickness variation on the predicted blade and endwall heat transfer are examined.
Results are shown for a three-dimensional Navier-Stokes analysis of both the flow and the surface heat transfer for turbine applications. Heat transfer comparisons are made with the experimental shock-tunnel data of Dunn and Kim, and with the data of Blair for the rotor of the large scale rotating turbine. The analysis was done using the steady-state, three-dimensional, thin-layer Navier-Stokes code developed by Chima, which uses a multistage Runge-Kutta scheme with implicit residual smoothing. An algebraic mixing length turbulence model is used to calculate turbulent eddy viscosity. The variation in heat transfer due to variations in grid parameters is examined. The effects of rotation, tip clearance, and inlet boundary layer thickness variation on the predicted blade and endwall heat transfer are examined.
Motivated by fire safety concerns and the advent of long-term micro-gravity facilities, a cooperative program has been developed to study the mechanisms and material properties that control flow assisted (co-current) flame spread. This program has used as a common fire scenario a reacting steady-state boundary layer. Preliminary studies explored the aerodynamics of a reacting boundary layer by simulating a condensed fuel by means of a gas burner. Stability curves for ethane air flames were obtained and different burning regimes were identified. An important feature of this study was the independent identification of the different mechanisms leading to the instability of the flow. It was observed that fuel injection velocity and thermal expansion independently contributed to the separation of the flow at the leading edge of the burner. The occurrence of separation resulted in complex three-dimensional flow patterns that have a dominant effect on critical fire safety parameters such as the stand-off distance and flame length. This work was extended to a solid fuel (PMMA) leading to a Sounding Rocket experiment (Mini-Texus-6). The solid phase showed similar flow patterns, mostly present at low flow velocities (<100 mm/s) but the results clearly demonstrated that the thermal balance at the pyrolyzing fuel surface is the dominant mechanism that controls both stand-off distance and flame length. This thermal balance could be described in a global manner by means of a total mass transfer or "B" number. This "B" number incorporates surface re-radiation, radiative feedback and in-depth heat conduction as first prescribed by Emmons. The mass transfer number becomes the single parameter that determines the evolution of these fire safety variables (flame length, stand-off distance) and therefore can be used as a ranking criterion to assess the flammability of materials. The particular configuration is representative of the NASA upward flame spread test (Test 1) therefore this approach can be used in the interpretation of the results obtained from this test. Nevertheless, complete validation of this approach has not been fully achieved due, mainly because all the measurements necessary to compare with the theoretical predictions have not been obtained. Following these studies two different directions have been taken. The first attempts to elucidate the details of the gas phase combustion reaction and the associated flow field by means of quantitative and qualitative measurements. The second approach, a more practical one, is to apply this methodology to the assessment of material flammability. The former is currently being conducted with a gas burner because it allows for easier control and longer experimentation time. The results obtained so far will be presented in more detail. The latter is a new program therefore only a brief summary of the objectives will be presented.
Two-phase thermal management subsystems that take advantage of both the sensible and latent heat of a working fluid can potentially yield significant enhancements in overall performance by adopting heat transfer processes that are based on phase transition like boiling and condensation. Performance of terrestrial two-phase flow systems may be predictable because the hydrodynamic and body forces are understood, however, in microgravity, which is predominant during planetary space travel, forces that are masked by the strong body force on Earth (gravitational or buoyancy force) reappear with different magnitude and influence. The need arose for a facility that provides for two-phase flow with phase transition testing in microgravity. The Flow Boiling and Condensation Experiment (FBCE) is a facility that was launched to the International Space Station in August of 2021 and is in operation since February of 2022. This facility enables investigators to perform two-phase flow and phase transition research in flow boiling and condensation. Along with the test module that is experiment specific, the FBCE system consists of the fluid, avionics, and software subsystems. Currently two test modules, namely, the Flow Boiling Module (FBM) and the Condensation Module for Heat Transfer (CM-HT) are available. A third module, the Transfer Line test Module (TL) is being developed. The fluid subsystem conditions and delivers the fluid at the desired thermodynamic state to the test module. It consists of two fluid modules and a heater module that are connected by flex hoses for fluid circulation and by data and electrical cables for control and data acquisition. Two avionics modules acquire pressure and temperature data from various sensors in the flow loop. For FBM, a high-speed camera is available to acquire images of the boiling process. Experiments are operated autonomously by software and are based on an Experiment Parameters Master Table (EPMT) that is uploaded to ISS and is executed by the FBCE flight software. This presentation briefly introduces the objectives of FBCE and provides a system description of the experiment onboard of the ISS/Fluid Integrated Rack (FIR). Results of the test campaign carried out using the FBM are presented. Specifically, microgravity flow boiling of n-perfluorohexane (test fluid) is discussed with subcooled inlet conditions in a single-side-heated rectangular channel of dimensions 114.6-mm heated length, 2.5-mm heated width, and 5.0-mm height. Key operating parameters investigated are mass velocity (199.90 – 3200.13 kg/m2s), inlet subcooling (0.10 – 45.76°C), and inlet pressure (113.30 – 164.29 kPa). Image sequences acquired via high-speed-video are shown to elucidate the interfacial flow physics. The effects of various parameters on flow boiling heat transfer in microgravity, from the onset of boiling to the critical heat flux are discussed. Heat transfer results are presented in terms of flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and parametric trends of local and averaged heat transfer coefficient, and the critical heat flux.
Correlation of second-order rate constants for many reactions involving electron transfer between organic molecules, solvated electron reactions, iodine diffusion coefficients, and triplet state electron transfer reactions has been made with the heat of mixing parameter (HMP) for the aqueous binary solvent systems. The aqueous binary solvents studied are those containing methanol or ethanol (type I solvent); 1-propanol or tert-butyl alcohol (type II solvent); or sucrose or glycerol (type III solvent). A plot of the HMP vs. the diffusion parameter for each reaction yields superimposable curves for these reactions in a particular solvent mixture over the entire solvent mixture range, irrespective of the value of the reaction's rate constant or diffusion coefficient in water.
Impulsive plasma penetration has been proposed as a means by which solar wind plasma may cross the magnetopause and enter the magnetosphere. It has been predicted that plasma entry by this means will be correlated with a magnetopause penetration parameter. One of the most important signatures of plasma entry into the magnetosphere are flux transfer events (FTEs). If plasma penetration is an important mechanism in plasma entry, then it would be expected that the occurrence of flux transfer events will correlate with the penetration parameter. On the other hand, if reconnection is, as widely believed, the dominant mechanism for the entry of solar-wind plasma into the magnetosphere, then plasma entry will be correlated with the north/south IMF component. AMPTE UKS data show that the occurrence of FTE signatures is better correlated with the direction of the IMF than with the penetration parameter. Thus, impulsive plasma penetration is unlikely to be the mechanism responsible for FTEs.
A fundamental problem of interest for crystal growth in micro-gravity applications involves the mixed convection heat transfer from a sphere in a uniform flow of fluid at a differing temperature. Under the combined influence of the imposed free stream as well as an induced buoyancy force due to thermal expansion of the fluid, the heat transfer from the sphere will be different from that of either the pure forced convection flow or the pure free convection flow. For the present study, the method of matched asymptotic expansions is applied to the laminar flow problem of an impulsively heated, impulsively started sphere in an originally quiescent fluid. Time series expansions are developed for the dependent variables by acknowledging the existence of two distinct regions: one, an inner region, near the sphere, in which viscous effects are significant; and two, an outer region in which the fluid may be treated as inviscid. The time series expansions are developed in terms of the Reynolds number and Richardson number (Buoyancy Parameter), and the relevant heat transfer and drag coefficients are calculated and plotted.
The systems engineering description of a wideband communications channel is provided which is based upon the fundamental propagation aspects of the problem. In particular, the well known time variant description of a channel is formulated from the basic multiple scattering processes that occur in a random propagation medium. Such a connection is required if optimal processing methods are to be applied to mitigate the deleterious random fading and multipathing of the channel. An example is given which demonstrates how the effective bandwidth of the channel is diminished due to atmospheric propagation impairments.
A physics-based Nuclear Thermal Propulsion (NTP) Testing Reference Design (TRD) power balance model was coded in Simulink to investigate engine performance for various design and parameter modifications. Since the primary mode of heat transfer in NTP engines is convective, the convective heat transfer coefficient (HTC) is a key parameter that requires accurate representation. The industry standard Westinghouse correlation has an uncertainty of ±20% which was investigated in this study. The results showed that a 20% decrease in the HTC led to a 4.14% increase in maximum fuel temperature while a 20% decrease in the HTC led to a 1.81% decrease in maximum fuel temperature suggesting that narrowing the uncertainty of this correlation through experimental work would be a critical step in the development of NTP engines. Furthermore, a maximum fuel temperature relationship with specific impulse was developed for the TRD engine which showed potential engine operation between specific impulse values of 715 and 900 seconds with minimal changes to the engine design. This graph could be useful for high level vehicle performance estimations for fuel types with different maximum operating temperatures.
The effects of cylindrical leading edge sweep on surface pressure and heat transfer rate for swept shock wave interference were investigated. Experimental tests were conducted in the Calspan 48-inch Hypersonic Shock Tunnel at a nominal Mach number of 8, nominal unit Reynolds number of 1.5 x 10 to the 6th power per foot, leading edge and incident shock generator sweep angles of 0, 15, and 30 deg, and incident shock generator angle-of-attack fixed at 12.5 deg. Detailed surface pressure and heat transfer rate on the cylindircal leading edge of a swept shock wave interference model were measured at the region of the maximum surface pressure and heat transfer rate. Results show that pressure and heat transfer rate on the cylindrical leading edge of the shock wave interference model were reduced as the sweep was increased over the range of tested parameters. Peak surface pressure and heat transfer rate on the cylinder were about 10 and 30 times the undisturbed flow stagnation point value, respectively, for the 0 deg sweep test. A comparison of the 15 and 30 deg swept results with the 0 deg swept results showed that peak pressure was reduced about 13 percent and 44 percent, respectively, and peak heat transfer rate was reduced about 7 percent and 27 percent, respectively.
The trajectory data are presented chronologically and are organized by holding the arrival date constant while varying the Earth departure date in increments of 10 days. Upon completion of the specified range of Earth departure dates, the arrival date is incremented and the range of departure dates is repeated. For long trip times, where the variation of the trajectory parameters is relatively small, the size of the increment of the arrival date is increased. The range of departure and arrival dates and their corresponding increments are given in Table 5-1 for each launch opportunity. The criterion for the selection of these dates is that they encompass the region in which the Earth departure hyperbolic excess speed is less than or equal to 0.65 EMOS. There are two lines of print for each trajectory (departure-date/arrival date pair). In the first line the two left most columns contain the dates of departure and arrival. The next 18 columns of the first line can be divided into three groups: six columns of data related to departure, six columns pertinent to the heliocentric phase of the mission, and six columns related to arrival at the target planet. The second line of print contains, respectively, the Delta V requirements for departure and arrival, the total Delta V requirement, the heliocentric transfer trajectory type, and four parameters defining conditions at arrival. The value computed for the arrival Delta V is for entry into a circular orbit. The radius selected for this orbit, while necessarily somewhat arbitrary, is representative of the broad range of orbit radii which tend to minimize the arrival Delta V for the range of excess speeds between 0.1 and 0.8 EMOS. The value selected for Uranus and Neptune is 3 planet radii. Significant reductions in the computed Delta V can be realized by assuming entry into an elliptical orbit having a periapsis radius equal to the selected circular-orbit radius. The magnitude of the reduction can be determined from Figures 2-6 and 2-8 i n Section 2.
In this work we report the results of a preliminary investigation focused on evaluating the thermal-hydraulic performance of R455A (CO2, R32, R1234yf – 3/21.5/75.5) in a commercial brazed plate heat exchanger. Three key parameters were measured which were refrigerant flow maldistribution, overall refrigerant side heat transfer coefficient, and total refrigerant side pressure drop. These same parameters were also measured for R134a which served as a base-line refrigerant for comparison. The experimental test matrix was set up to investigate the influence of refrigerant mass flux (3 ≤ Gref ≤ 5 kg m-2 s-1), refrigerant inlet vapor quality (0.1 ≤ xin ≤ 0.4), and refrigerant exit superheat (5 ≤ Tsh ≤ 10 K) on the above mentioned thermal-hydraulic parameters. For these experimental conditions, on average, the overall refrigerant side heat transfer coefficients for R455A were 45.2% lower than those of R134a. Refrigerant flow maldistribution for R455A, on average, was observed to be 62.6 % higher than that measured for R134a, and the total pressure drop across the heat exchanger (minor and frictional losses) for R455A was on average 26.5% higher than that measured for R134a. These initial results set the stage for a comprehensive investigation focused on further investigating the thermal-hydraulic performance of zeotropic refrigerants in brazed plate heat exchangers.
A new diagnostic tool is developed for examining relationships between the synoptic scale circulation and regional temperature distributions in GCMs. The 4 x 5 deg GISS GCM is shown to produce accurate simulations of the variance in the synoptic scale sea level pressure distribution over the U.S. An analysis of the observational data set from the National Meteorological Center (NMC) also shows a strong relationship between the synoptic circulation and grid point temperatures. This relationship is demonstrated by deriving transfer functions between a time-series of circulation parameters and temperatures at individual grid points. The circulation parameters are derived using rotated principal components analysis, and the temperature transfer functions are based on multivariate polynomial regression models. The application of these transfer functions to the GCM circulation indicates that there is considerable spatial bias present in the GCM temperature distributions. The transfer functions are also used to indicate the possible changes in U.S. regional temperatures that could result from differences in synoptic scale circulation between a 1XCO2 and a 2xCO2 climate, using a doubled CO2 version of the same GISS GCM.