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.
Artificial dissipation in modern computational fluid dynamics is something of an art. In spite of the dramatic strides made over the last decade, dissipation continues to lack a quantitative theory with predictive power. Such a theory, based on entropy production, is presented. An important side-effect is nonlinear stability for schemes which satisfy certain conditions.
Capillary Pumped Loop (CPL) Modeler computer program is amalgamation of software that mathematically models performance of CPL system and its environment. Two-phase heat-transport device capable of transferring heat loads efficiently over large distance with little temperature differential. Utilizes surface-tension forces established in fine-pore capillary wick to circulate working fluid, requiring no external pumpling power. Predicts steady-state or quasi-steady behavior of CPL embedded in spacecraft or other thermal environment. Also predicts location of liquid/vapor interface in each condenser. Written in VAX/VMS FORTAN 77.
An analysis was performed to predict the thermal distortion of the solar dynamic concentrator for Space Station Freedom in low earth orbit and to evaluate the effects of that thermal distortion on concentrator on-orbit performance. The analysis required substructural finite element modeling of critical concentrator structural subsystems, structural finite element modeling of the concentrator, mapping of thermal loading onto the structural finite element model, and the creation of specialized postprocessors to assist in interpreting results. Concentrator temperature distributions and thermally induced displacements and slope errors and the resulting receiver flux distribution profiles are discussed. Results determined for a typical orbit indicate that concentrator facet rotations are less than 0.2 mrad and that the change in facet radius due to thermal flattening is less than 5 percent. The predicted power loss due to thermal distortion effects is less than 0.3 percent. As a consequence the thermal distortions of the solar dynamic concentrator in low earth orbit will have a negligible effect on the flux distribution profiles within the receiver.
A commercial fluid mechanics program, FLUENT, has been applied to the modeling of a hot-filament diamond deposition reactor. Streamlines and contours of constant temperature and species concentrations are obtained for practical reactor geometries and conditions. The modeling is presently restricted to two-dimensional simulations and to a chemical mechanism of ten independent homogeneous and surface reactions. Comparisons are made between predicted power consumption, substrate temperature, and concentrations of atomic hydrogen and methyl-radical with values taken from the literature. The results to date indicate that the modeling can aid in the rational design and analysis of practical reactor configurations.
The dynamic model is a method for computing the coefficient C in Smagorinsky's model for the subgrid-scale stress tensor as a function of position from the information already contained in the resolved velocity field rather than treating it as an adjustable parameter. A variational formulation of the dynamic model is described that removes the inconsistency associated with taking C out of the filtering operation. This model, however, is still unstable due to the negative eddy-viscosity. Next, three models are presented that are mathematically consistent as well as numerically stable. The first two are applicable to homogeneous flows and flows with at least one homogeneous direction, respectively, and are, in fact, a rigorous derivation of the ad hoc expressions used by previous authors. The third model in this set can be applied to arbitrary flows, and it is stable because the C it predicts is always positive. Finally, a model involving the subgrid-scale kinetic energy is presented which attempts to model backscatter. This last model has some desirable theoretical features. However, even though it gives results in LES that are qualitatively correct, it is outperformed by the simpler constrained variational models. It is suggested that one of the constrained variational models should be used for actual LES while theoretical investigation of the kinetic energy approach should be continued in an effort to improve its predictive power and to understand more about backscatter.
A statistical study of diffuse energetic ion events and their related waves upstream of the Earth's bow shock was performed using data from the Active Magnetospheric Particle Tracer Explorers/Ion Release Module (AMPTE/IRM) satellite over two 5-month periods in 1984 and 1985. The data set was used to test the assumption in the self-consistent model of the upstream wave and particle populations by Lee (1982) that the particle acceleration through hydromagnetic waves and the wave generation are directly coupled. The comparison between the observed wave power and the wave power predicted on the observed energetic particle energy density and solar wind parameters results in a high correlation coefficient of about 0.89. The intensity of diffuse ions falls off approximately exponentially with the distance upstream from the bow shock parallel to the magnetic field with e-folding distances which vary from approximately 3.3 R(sub E) to approximately 11.7 R(sub E) over the energy range from 10 keV/e to 67.3 keV/e for both protons and alpha particles. After normalizing the upstream particle densities to zero bow shock distance by using these exponential variations, a good correlation (0.7) of the density of the diffuse ions with the solar wind density was found. This supports the suggestion that the solar wind is the source of the diffuse ions. Furthermore, the spectral slope of the diffuse ions correlates well with the solar wind velocity component in the direction of the interplanetary magnetic field (0.68 and 0.66 for protons and alpha particles) which concurs with the notion that the solar wind plays an important role in the acceleration of the upstream particles.
Nonbuoyant flames offer themselves as an attractive and promising platform to gain a better understanding of soot mechanisms. The effects of buoyancy can be eliminated temporarily in drop towers which sustain brief intervals of reduced gravity-typically lower than 10(exp -3)g- extending up to several seconds at a time. Microgravity facilities have been employed to show that nonbuoyant flames are longer, wider and sootier than their normal-gravity counterparts. Sunderland et al. recently verified the existence of smoke point in laminar nonbuoyant flames. As reported, microgravity flames operating above their smoke point displayed a blunt tip and much broader soot-containing regimes in comparison to their buoyant counterparts. Mortazavi et al. established that residence times in microgravity laminar jet diffusion flames with Re=0(100) tend to be proportional to burner diameter and inversely proportional to burner exit velocity. This offers the capability to alter residence times in nonbuoyant laminar jet diffusion flames when varying the burner exit diameters and velocities. Megaridis et al. presented a quantitative definition of the soot-field structure within laminar microgravity jet diffusion flames which operated well above their smoke point. The experimental methodology involved a full-field laser-light extinction technique and jet diffusion flames of nitrogen-diluted (50% vol.) acetylene fuel burning in quiescent air at atmospheric pressure. The work was conducted at the 2.2s drop tower of the NASA Lewis Research Center (NASA-LeRC). Parallel work on 1-g flames was also presented in (6) to facilitate comparisons on the effect of gravity on the soot fields. As reported, the soot spatial distributions in 0-g flames did not change in a detectable manner after 1s within a typical 2.2s experiment. During that period, the soot field was shown to sustain a pronounced annular structure throughout the luminous nonbuoyant-flame zone. The maximum soot volume fraction measured at 0-g was nearly a factor of two higher than that at 1-g, thus confirming the enhanced sooting tendency of nonbuoyant flames. Greenberg and Ku presented a similar study and reported trends that matched those of for the 50% (vol.) nitrogen-diluted acetylene fuel. Furthermore, they examined pure acetylene flames and reported similar trends with respect to the influence of gravity on maximum soot volume fractions and flame cross-section-averaged soot loadings. Both studies clearly demonstrated the improved spatial resolution of microgravity flames compared to their normal-gravity counterparts. The current study evaluates the influence of moderate fuel preheat on soot formation within 0-g laminar gas jet diffusion flames. While fuel temperature variations have little influence on residence times in 1-g, they have a much more significant effect in 0-g. The primary objective of this program is to quantify this effect and its consequences on sooting by comparing soot volume fraction distributions under preheated and unpreheated-fuel conditions. Furthermore, the current work aims at expanding the limited soot database available for nonbuoyant flames. Soot fields in such flames can be used to perform additional tests of recently developed soot sub-models which have the potential to become powerful predicting tools in combustion design.
While the response surface method is an effective method in engineering optimization, its accuracy is often affected by the use of limited amount of data points for model construction. In this chapter, the issues related to the accuracy of the RS approximations and possible ways of improving the RS model using appropriate treatments, including the iteratively re-weighted least square (IRLS) technique and the radial-basis neural networks, are investigated. A main interest is to identify ways to offer added capabilities for the RS method to be able to at least selectively improve the accuracy in regions of importance. An example is to target the high efficiency region of a fluid machinery design space so that the predictive power of the RS can be maximized when it matters most. Analytical models based on polynomials, with controlled level of noise, are used to assess the performance of these techniques.
Space radiation shield applications and studies of cosmic ray propagation in the Galaxy require reliable cross sections to calculate spectra of secondary particles and yields of the isotopes produced in nuclear reactions induced both by particles and nuclei at energies from threshold to hundreds of GeV per nucleon. Since the data often exist in a very limited energy range or sometimes not at all, the only way to obtain an estimate of the production cross sections is to use theoretical models and codes. Recently, we have developed improved versions of the Cascade-Exciton Model (CEM) of nuclear reactions: the codes CEM97 and CEM2k for description of particle-nucleus reactions at energies up to about 5 GeV. In addition, we have developed a LANL version of the Quark-Gluon String Model (LAQGSM) to describe reactions induced both by particles and nuclei at energies up to hundreds of GeVhucleon. We have tested and benchmarked the CEM and LAQGSM codes against a large variety of experimental data and have compared their results with predictions by other currently available models and codes. Our benchmarks show that CEM and LAQGSM codes have predictive powers no worse than other currently used codes and describe many reactions better than other codes; therefore both our codes can be used as reliable event-generators for space radiation shield and cosmic ray propagation applications. The CEM2k code is being incorporated into the transport code MCNPX (and several other transport codes), and we plan to incorporate LAQGSM into MCNPX in the near future. Here, we present the current status of the CEM2k and LAQGSM codes, and show results and applications to studies of cosmic ray propagation in the Galaxy.
Electromagnetic propagation models for signal strength prediction within aircraft cabins are essential for evaluating and designing a wireless communication system to be implemented onboard aircraft. A model was developed using Wireless Valley's SitePlanner; which is commercial grade software intended for predictions within office buildings. The performance of the model was evaluated through a comparison with test data measurements taken on several aircraft. The comparison concluded that the model can accurately predict power propagation within the cabin. This model can enhance researchers understanding of power propagation within aircraft cabins and will aid in future research.
The density profile of a critical fluid near a solid surface is expected to show an universal shape. This is known as critical adsorption. The measurement of this effect, especially close to the critical point, is often obscured by gravity. We were able to separate the gravitational effect from critical adsorption by using two capacitors, one with a large gap and one with a small gap of approximately 2 m. Within the uncertainty in the measurement, our data, which ranges between 10(exp -3) to 2 x 10(exp -6) in reduced temperatures, is consistent with the predicted power law dependence. This work is carried out in collaboration with Rafael Garcia, Sarah Scheidemantel and Klaus Knorr. It is funded by NASA's office of Biological and Physical Researchunder.
We report further results from our ongoing assessment of magnetogram-based measures of active-region nonpotentiality and size as predictors of coronal mass ejections (CMEs). We have devised improved generalized measures of active-region nonpotentiality that apply to active regions of any degree of magnetic complexity, rather than being limited to bipolar active regions as our initial measures were. From a set of approx.50 active-regions, we have found that measures of total nonpotentiality have a 75-80% success rate n predicting whether an active region will produce a CME in 2 days after the magnetogram. This makes measures of total nonpotentiality a better predictor than either active-region size, or active region twist (size-normalized nonpotentiality), which have a approx.65% success rates. We have also found that we can measure from the line-of-sight magnetograms an active region's total nonpotentiality and the size, which allows use to use MDI to evaluate these quantities for 4-5 consecutive days for each active region, and to investigate if there is some combination of size and total nonpotentiality that have a stronger predictive power than does total nonpotentiality. This work was funded by NASA through its LWS TR&T Program and its Solar and Heliospheric Physics SR&T Program, and by NSF through its Solar Terrestrial Research and SHINE programs.
Statistical Methods outlined in [Ladbury, TNS20071 can be generalized for Monte Carlo Rate Calculation Methods Two Monte Carlo Approaches: a) Rate based on vendor-supplied (or reverse-engineered) model SEE testing and statistical analysis performed to validate model; b) Rate calculated based on model fit to SEE data Statistical analysis very similar to case for CREME96. Information Theory allows simultaneous consideration of multiple models with different complexities: a) Model with lowest AIC usually has greatest predictive power; b) Model averaging using AIC weights may give better performance if several models have similar good performance; and c) Rates can be bounded for a given confidence level over multiple models, as well as over the parameter space of a model.
The goal of this paper is to derive the equation for the turbulence dissipation rate epsilon for a shear-driven flow. In 1961, Davydov used a one-point closure model to derive the epsilon equation from first principles but the final result contained undetermined terms and thus lacked predictive power. Both in 1987 and in 2001, attempts were made to derive the epsilon equation from first principles using a two-point closure, but their methods relied on a phenomenological assumption. The standard practice has thus been to employ a heuristic form of the equation that contains three empirical ingredients: two constants, c(sub 1 epsilon), and c(sub 2 epsilon), and a diffusion term D(sub epsilon) In this work, a two-point closure is employed, yielding the following results: 1) the empirical constants get replaced by c(sub 1), c(sub 2), which are now functions of Kappa and epsilon; 2) c(sub 1) and c(sub 2) are not independent because a general relation between the two that are valid for any Kappa and epsilon are derived; 3) c(sub 1), c(sub 2) become constant with values close to the empirical values c(sub 1 epsilon), c(sub epsilon 2), (i.e., homogenous flows); and 4) the empirical form of the diffusion term D(sub epsilon) is no longer needed because it gets substituted by the Kappa-epsilon dependence of c(sub 1), c(sub 2), which plays the role of the diffusion, together with the diffusion of the turbulent kinetic energy D(sub Kappa), which now enters the new equation (i.e., inhomogeneous flows). Thus, the three empirical ingredients c(sub 1 epsilon), c(sub epsilon 2), D (sub epsilon)are replaced by a single function c(sub 1)(Kappa, epsilon ) or c(sub 2)(Kappa, epsilon ), plus a D(sub Kappa)term. Three tests of the new equation for epsilon are presented: one concerning channel flow and two concerning the shear-driven planetary boundary layer (PBL).
Retrievals of aerosol optical depth (AOD) and related parameters from satellite measurements typically involve prescribed models of aerosol size and composition, and are therefore dependent on how well these models are able to represent the radiative behaviour of real aerosols, This study uses aerosol volume size distributions retrieved from Sun-photometer measurements at 11 Aerosol Robotic Network (AERONET) island sites, spread throughout the world's oceans, as a basis to define such a model for unpolluted maritime aerosols. Size distributions are observed to be bimodal and approximately lognormal, although the coarse mode is skewed with a long tail on the low-radius end, The relationship of AOD and size distribution parameters to meteorological conditions is also examined, As wind speed increases, so do coarse-mode volume and radius, The AOD and Angstrom exponent (alpha) show linear relationships with wind speed, although there is considerable scatter in all these relationships, limiting their predictive power. Links between aerosol properties and near-surface relative humidity, columnar water vapor, and sea surface temperature are also explored. A recommended bimodal maritime model, which is able to reconstruct the AERONET AOD with accuracy of order 0.01-0.02, is presented for use in aerosol remote sensing applications. This accuracy holds at most sites and for wavelengths between 340 nm and 1020 nm. Calculated lidar ratios are also provided, and differ significantly from those currently used in Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) processing.
To facilitate locating archaeological sites before they are compromised or destroyed, we are developing approaches for generating maps of probable archaeological sites, through detecting subtle anomalies in vegetative cover, soil chemistry, and soil moisture by analyzing remotely sensed data from multiple sources. We previously reported some success in this effort with a statistical analysis of slope, radar, and Ikonos data (including tasseled cap and NDVI transforms) with Student's t-test. We report here on new developments in our work, performing an analysis of 8-band multispectral Worldview-2 data. The Worldview-2 analysis begins by computing medians and median absolute deviations for the pixels in various annuli around each site of interest on the 28 band difference ratios. We then use principle components analysis followed by linear discriminant analysis to train a classifier which assigns a posterior probability that a location is an archaeological site. We tested the procedure using leave-one-out cross validation with a second leave-one-out step to choose parameters on a 9,859x23,000 subset of the WorldView-2 data over the western portion of Ft. Irwin, CA, USA. We used 100 known non-sites and trained one classifier for lithic sites (n=33) and one classifier for habitation sites (n=16). We then analyzed convex combinations of scores from the Archaeological Predictive Model (APM) and our scores. We found that that the combined scores had a higher area under the ROC curve than either individual method, indicating that including WorldView-2 data in analysis improved the predictive power of the provided APM.
Astronauts experience sensorimotor disturbances during the initial exposure to microgravity and during the readapation phase following a return to a gravitational environment. These alterations may lead to disruption in the ability to perform mission critical functional tasks during and after these gravitational transitions. Astronauts show significant inter-subject variation in adaptive capability following gravitational transitions. The ability to predict the manner and degree to which each individual astronaut will be affected would improve the effectiveness of a countermeasure comprised of a training program designed to enhance sensorimotor adaptability. Due to this inherent individual variability we need to develop predictive measures of sensorimotor adaptability that will allow us to predict, before actual space flight, which crewmember will experience challenges in adaptive capacity. Thus, obtaining this information will allow us to design and implement better sensorimotor adaptability training countermeasures that will be customized for each crewmember's unique adaptive capabilities. Therefore the goals of this project are to: 1) develop a set of predictive measures capable of identifying individual differences in sensorimotor adaptability, and 2) use this information to design sensorimotor adaptability training countermeasures that are customized for each crewmember's individual sensorimotor adaptive characteristics. To achieve these goals we are currently pursuing the following specific aims: Aim 1: Determine whether behavioral metrics of individual sensory bias predict sensorimotor adaptability. For this aim, subjects perform tests that delineate individual sensory biases in tests of visual, vestibular, and proprioceptive function. Aim 2: Determine if individual capability for strategic and plastic-adaptive responses predicts sensorimotor adaptability. For this aim, each subject's strategic and plastic-adaptive motor learning abilities are assessed using a test of locomotor function designed specifically to delineate both mechanisms. Aim 3: Develop predictors of sensorimotor adaptability using brain structural and functional metrics. We will measure individual differences in regional brain volumes (structural MRI), white matter integrity (diffusion tensor imaging, or DTI), functional network integrity (resting state functional connectivity MRI), and sensorimotor adaptation task-related functional brain activation (functional MRI). We decided to complete the data collection for Specific Aims 1, 2 and 3 simultaneously on the same subjects to increase data capture. By having the same subjects perform all three specific aims we can enhance our ability to detect how a wider range of factors can predict adaptability in a specific individual. This provides a much richer database and potentially a better understanding of the predictive power of the selected factors. In this presentation I will discuss preliminary data obtained to date.
Emergency medicine physicians recently adopted the use of ultrasonography to quickly measure optic nerve sheath diameter (ONSD) as concomitant with increased intracranial pressure. NASA Space and Clinical Operations Division has been using ground and on-orbit ultrasound capabilities since 2009 to consider this anatomical measure as a proxy for intracranial pressure in the microgravity environment. In the terrestrial emergency room population, an ONSD greater than 0.59 cm is considered highly predictive of elevated intracranial pressure. However, this cut-off limit is not applicable to the spaceflight setting since over 50% of US Operating Segment (USOS) astronauts have an ONSD greater than 0.60 cm even before missions. Crew Surgeon clinical decision-making is complicated by the fact that many astronauts have history of previous spaceflights. Data will be presented characterizing the distribution of baseline ONSD in the astronaut corps, longitudinal trends in-flight, and the predictive power of this measure related to increased intracranial pressure outcomes.