SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A linearized version of the Hartree-Fock method is used as a probe to investigate phase transitions in fermionic systems with many-body interactions. An application to a new exactly solvable model which includes two- and three-body forces is shown.
Rossby waves are difficult to detect with in situ methods. However, as we show in this paper, they can be clearly identified in multi-parameters in multi-mission satellite observations of sea surface height (SSH), sea surface temperature (SST) and ocean color observations of chlorophyll-a (chl-a), as well as 1/12-deg global HYbrid Coordinate Ocean Model (HYCOM) simulations of SSH, SST and sea surface salinity (SSS) in the Indian Ocean. While the surface structure of Rossby waves can be elucidated from comparisons of the signal in different sea surface parameters, models are needed to gain direct information about how these waves affect the ocean at depth. The first three baroclinic modes of the Rossby waves are inferred from the Fast Fourier Transform (FFT), and two-dimensional Radon Transform (2D RT). At many latitudes the first and second baroclinic mode Rossby wave phase speeds from satellite observations and model parameters are identified.
The nucleus of comet Halley was observed by the Halley Multicolor Camera (HMC) during the last minutes of the flyby sequence before reaching the point of closest approach. The phase angle change during the observational period was less than 17 deg, so HMC data can only define a two-dimensional contour of the nucleus at a certain time. The overall geometry during the HMC observations is outlined and the findings are related to information available from other sources at different times. The nucleus of the comet is described as a triaxial ellipsoid with major axes of 16, 10, and 9 km. The rotation axis is directed towards delta 1950 = -(40 + or - 5) deg and alpha 1950 = (50 + or - 10) deg. The rotation with a period of P = (54 + or - 1) is in a prograde sense with respect to the orbital motion.
Large Eddy Simulation (LES) is conducted of a three-dimensional temporal mixing layer whose stream is initially laden with liquid drops which may evaporate during the simulation. The gas-phase equations are written in Eulerian frame for two perfect gas species (carrier gas and vapor emanating from the drops), while the liquid-phase equations are written in a Lagrangian frame.
We present formal models and results of formal analysis of two different asynchronous systems. We first examine a mid-value select module that merges the signals coming from three different sensors that are each asynchronously sampling the same input signal. We then consider the phase locking protocol proposed by Daly, Hopkins, and McKenna. This protocol is designed to keep a set of non-faulty (asynchronous) clocks phase locked even in the presence of Byzantine-faulty clocks on the network. All models and verifications have been developed using the SAL model checking tools and the Hybrid SAL abstractor.
Results are reported for VLBI observations of 3C 345 at a wavelength of 3.8 cm with left circular polarization. Sixteen sets of correlated flux densities and three sets of closure phases were obtained from observations made over a three-year period with six independent interferometers formed from up to four antennas. The mathematical models employed to represent the radio brightness distribution of the source are described, and the data are analyzed on the basis of a two-component model. The results show that the compact source in 3C 345 consists of two components whose separation increased by about 0.3 milliarcsec over the observational period, that the smaller and weaker component is to the east and south of the larger and stronger one, and that the apparent expansion velocity would be about 2.5c for a redshift of 0.6, a Hubble constant of 60 km/sec per Mpc, and a deceleration parameter of unity. It is concluded that these findings are difficult to reconcile with the 'Christmas tree' model for the apparent change of separation with time over any given series of observations.
The main objective of the NASA sponsored Aerothermal Modeling Program, Phase 2--Element A, is to develop an improved numerical scheme for predicting combustor flow fields. This effort consists of the following three technical tasks. Task 1 involves the selection and evaluation of various candidate numerical techniques. Task 2 involves an in-depth evaluation of the selected numerical schemes. Task 3 involves the convection-diffusion scheme and the direct solver that will be incorporated in the NASA 3-D elliptic code (COM3S).
A robust feedback algorithm is presented for a near-minimum-fuel ascent of a two-stage launch vehicle operating in the equatorial plane. The development of the algorithm is based on the ideas of neighboring optimal control and can be derived into three phases. In phase 1, the formalism of optimal control is employed to calculate fuel-optimal ascent trajectories for a simple point-mass model. In phase 2, these trajectories are used to numerically calculate gain functions of time for the control(s), the total flight time, and possibly, for other variables of interest. In phase 3, these gains are used to determine feedback expressions for the controls associated with a more realistic model of a launch vehicle. With the Advanced Launch System in mind, all calculations are performed on a two-stage vehicle with fixed thrust history, but this restriction is by no means important for the approach taken. Performance and robustness of the algorithm is found to be excellent.
A three-dimensional model of cirrus cloud formation and evolution, including microphysical, dynamical, and radiative processes, was used to simulate cirrus observed in the FIRE Phase 2 Cirrus field program (13 Nov. - 7 Dec. 1991). Sulfate aerosols, solution drops, ice crystals, and water vapor are all treated as interactive elements in the model. Ice crystal size distributions are fully resolved based on calculations of homogeneous freezing of solution drops, growth by water vapor deposition, evaporation, aggregation, and vertical transport. Visible and infrared radiative fluxes, and radiative heating rates are calculated using the two-stream algorithm described by Toon et al. Wind velocities, diffusion coefficients, and temperatures were taken from the MAPS analyses and the MM4 mesoscale model simulations. Within the model, moisture is transported and converted to liquid or vapor by the microphysical processes. The simulated cloud bulk and microphysical properties are shown in detail for the Nov. 26 and Dec. 5 case studies. Comparisons with lidar, radar, and in situ data are used to determine how well the simulations reproduced the observed cirrus. The roles played by various processes in the model are described in detail. The potential modes of nucleation are evaluated, and the importance of small-scale variations in temperature and humidity are discussed. The importance of competing ice crystal growth mechanisms (water vapor deposition and aggregation) are evaluated based on model simulations. Finally, the importance of ice crystal shape for crystal growth and vertical transport of ice are discussed.
The goal of our research on flame spread across a pool of liquid fuel is the quantitative identification of the mechanisms that control the rate and nature of flame spread when the initial temperature of the liquid pool is below the fuel's flash point temperature. Besides numerous experiments in drop towers and 1 g laboratories, we have flown five microgravity (mu-g) experiments on sounding rockets. As described in earlier papers, the first three flights examined the effect of forced opposed airflow over a 2.5 cm deep x 2 cm wide x 30 cm long pool of 1-butanol in mu-g. It was found that the flame spread is much slower and steadier than in 1 g where flame spread has a pulsating character. It was speculated that the flame spread in mu-g resembled the character of pseudo-uniform spread in 1 g; Ito et al later confirmed this conclusively in 1 g experiments. Much of the mu-g flame is also farther from the surface, dimmer, and with less soot, when compared to the 1 g flame. Three-dimensional liquid-phase flow patterns that control the liquid preheating were discovered in both 1 g and mu-g. Our numerical model, restricted to two dimensions, had predicted faster, pulsating flame spread in mu-g for opposed airflow. In examining the differences in the dimensionality of the model and experiment, it was noted that the experiment allowed gas expansion in the lateral direction (across the width of the pool), for which the model could not account. Such lateral expansion could reduce the expansion in the forward and upward directions. Because only these latter directions could be modeled, it was decided to artificially reduce the gas thermal expansion in the predictions. When this was done, satisfactory agreement could be obtained between the predicted and observed spread rates and the steadiness of the spread in microgravity. In 1 g, however, the predicted flame spread character also changed to pseudo-uniform, which disagreed with our 1 g experiments where the spread is pulsating. It was then speculated that gas-phase buoyant flow might oppose the lateral gas expansion, so that the 1 g experiments retained their pulsating flame spread character. If this speculation was valid, a difference in lateral gas expansion should be observable when comparing 1 g and mu-g experiments. Specifically, it was anticipated that greater flow divergence caused by lateral expansion would be measured in mu-g in the absence of a buoyant flow directed towards the flame.
Flutter-free operation of advanced transonic fan designs continues to be a challenging task for the designers of aircraft engines. In order to meet the demands of increased performance and lighter weight, these modern fan designs usually feature low-aspect ratio shroudless rotor blade designs that make the task of achieving adequate flutter margin even more challenging for the aeroelastician. This is especially true for advanced forward swept designs that encompass an entirely new design space compared to previous experience. Fortunately, advances in unsteady computational fluid dynamic (CFD) techniques over the past decade now provide an analysis capability that can be used to quantitatively assess the aeroelastic characteristics of these next generation fans during the design cycle. For aeroelastic applications, Mississippi State University and NASA Glenn Research Center have developed the CFD code TURBO-AE. This code is a time-accurate three-dimensional Euler/Navier-Stokes unsteady flow solver developed for axial-flow turbomachinery that can model multiple blade rows undergoing harmonic oscillations with arbitrary interblade phase angles, i.e., nodal diameter patterns. Details of the code can be found in Chen et al. (1993, 1994), Bakhle et al. (1997, 1998), and Srivastava et al. (1999). To assess aeroelastic stability, the work-per-cycle from TURBO-AE is converted to the critical damping ratio since this value is more physically meaningful, with both the unsteady normal pressure and viscous shear forces included in the work-per-cycle calculation. If the total damping (aerodynamic plus mechanical) is negative, then the blade is unstable since it extracts energy from the flow field over the vibration cycle. TURBO-AE is an integral part of an aeroelastic design system being developed at Honeywell Engines, Systems & Services for flutter and forced response predictions, with test cases from development rig and engine tests being used to validate its predictive capability. A recent experimental program (Sanders et al., 2002) was aimed at providing the necessary unsteady aerodynamic and vibratory response data needed to validate TURBO-AE for fan flutter predictions. A comparison of numerical TURBO-AE simulations with the benchmark flutter data is given in Sanders et al. (2003), with the data used to guide the validation of the code and define best practices for performing accurate unsteady simulations. The agreement between the analyses and the predictions was quite remarkable, demonstrating the ability of the analysis to accurately model the unsteady flow processes driving stall-side flutter.
Colloidal particles are important in a broad range of applications in technology and materials processing. Due to high sensitivity to thermal fluctuations, a suspension of colloidal particles is considered as an assembly of super-atoms in condensed matter physics. This feature provides an excellent experimental model in statistical mechanics for fundamental studies of microstructural evolution and phase transformation at the individual particle level. However, using colloids to probe three-dimensional phase transitions on Earth remains a challenge as they can be influenced by the gravity force causing the particles to sediment and the solvent to flow. To suppress the undesirable gravity effects, experiments on concentrated colloids were carried out in the International Space Station (ISS) Light Microscopy Module. We will present results of these experiments.
After at least 6 years of quiescence, Anomalous X-ray Pulsar (AXP) 4U 0142+61 entered an active phase in 2006 March that lasted several months and included six X-ray bursts as well as many changes in the persistent X-ray emission. The bursts, the first seen from this AXP in >11 years of Rossi X-ray Timing Explorer monitoring, all occurred in the interval between 2006 April 6 and 2007 February 7. The burst durations ranged from 8-3x10(exp 3)s. The first five burst spectra are well modeled by blackbodies, with temperatures kT approx. 2 - 6 keV. However, the sixth burst had a complicated spectrum that is well characterized by a blackbody plus three emission features whose amplitude varied throughout the burst. The most prominent feature was at 14.0 keV. Upon entry into the active phase the pulsar showed a significant change in pulse morphology and a likely timing glitch. The glitch had a total frequency jump of (1.9+/-0.4)x10(exp -7) Hz, which recovered with a decay time of 17+/-2 days by more than the initial jump, implying a net spin-down of the pulsar. We discuss these events in the context of the magnetar model.
To better characterize and model optical data acquired from ground-based telescopes, the Optical Measurements Center (OMC) at NASA/JSC attempts to emulate illumination conditions seen in space using equipment and techniques that parallel telescopic observations and source-target-sensor orientations. The OMC uses a 75 Watt Xenon arc lamp as a solar simulator, an SBIG CCD camera with standard Johnson/Bessel filters, and a robotic arm to simulate an object's position and rotation. The laboratory uses known shapes, materials suspected to be consistent with the orbital debris population, and three phase angles to best match the lighting conditions of the telescope based data. The fourteen objects studied in the laboratory are fragments or materials acquired through ground-tests of scaled-model satellites/rocket bodies as well as material samples in more/less "flight-ready" condition. All fragments were measured at 10 increments in a full 360 rotation at 6 , 36 , and 60 phase angles. This paper will investigate published color photometric data for a series of orbital debris targets and compare it to the empirical photometric measurements generated in the OMC. Using the data acquired over specific rotational angles through different filters (B, V, R, I), a color index is acquired (B-R, R-I). Using these values and their associated lightcurves, this laboratory data is compared to observational data obtained on the 1 m telescope of the Astronomical Institute of the University of Bern (AUIB), the 0.9 m operated by the Small- and Medium-Aperture Research Telescope System (SMARTS) Consortium and the Curtis-Schmidt 0.6 m Michigan Orbital Debris Space Debris Telescope both located at Cerro Tololo Inter-American Observatory (CTIO). An empirical based optical characterization model will be presented to provide preliminary correlations between laboratory based and telescope-based data in the context of classification of GEO debris objects.
The design of effective new technologies to reduce aircraft propulsion noise is dependent on an understanding of the noise sources and noise generation mechanisms in the modern turbofan engine. In order to more fully understand the physics of noise in a turbofan engine, a comprehensive aeroacoustic wind tunnel test programs was conducted called the 'Source Diagnostic Test.' The text was cooperative effort between NASA and General Electric Aircraft Engines, as part of the NASA Advanced Subsonic Technology Noise Reduction Program. A 1/5-scale model simulator representing the bypass stage of a current technology high bypass ratio turbofan engine was used in the test. The test article consisted of the bypass fan and outlet guide vanes in a flight-type nacelle. The fan used was a medium pressure ratio design with 22 individual, wide chord blades. Three outlet guide vane design configurations were investigated, representing a 54-vane radial Baseline configuration, a 26-vane radial, wide chord Low Count configuration and a 26-vane, wide chord Low Noise configuration with 30 deg of aft sweep. The test was conducted in the NASA Glenn Research Center 9 by 15-Foot Low Speed Wind Tunnel at velocities simulating the takeoff and approach phases of the aircraft flight envelope. The Source Diagnostic Test had several acoustic and aerodynamic technical objectives: (1) establish the performance of a scale model fan selected to represent the current technology turbofan product; (2) assess the performance of the fan stage with each of the three distinct outlet guide vane designs; (3) determine the effect of the outlet guide vane configuration on the fan baseline performance; and (4) conduct detailed flowfield diagnostic surveys, both acoustic and aerodynamic, to characterize and understand the noise generation mechanisms in a turbofan engine. This paper addresses the fan and stage aerodynamic performance results from the Source Diagnostic Test.
We performed deep observations to search for radio pulsations in the directions of 375 unassociated Fermi Large Area Telescope γ-ray sources using the Giant Metrewave Radio Telescope (GMRT) at 322 and 607 MHz. In this paper we report the discovery of three millisecond pulsars (MSPs), PSR J0248+4230, PSR J1207–5050, and PSR J1536–4948. We conducted follow-up timing observations for ∼5 yr with the GMRT and derived phase-coherent timing models for these MSPs. PSR J0248+4230 and J1207–5050 are isolated MSPs having periodicities of 2.60 ms and 4.84 ms. PSR J1536–4948 is a 3.07 ms pulsar in a binary system with an orbital period of ∼62 days about a companion of a minimum mass of 0.32 M⊙. We also present multifrequency pulse profiles of these MSPs from the GMRT observations. PSR J1536–4948 is an MSP with an extremely wide pulse profile having multiple components. Using the radio timing ephemeris we subsequently detected γ-ray pulsations from these three MSPs, confirming them as the sources powering the γ-ray emission. For PSR J1536–4948 we performed combined radio–γ-ray timing using ∼11.6 yr of γ-ray pulse times of arrival (TOAs) along with the radio TOAs. PSR J1536–4948 also shows evidence for pulsed γ-ray emission out to above 25 GeV, confirming earlier associations of this MSP with a ≥10 GeV point source. The multiwavelength pulse profiles of all three MSPs offer challenges to models of radio and γ-ray emission in pulsar magnetospheres.
Considerations leading to the numerical design of the GLAS fourth-order global atmospheric model are discussed, including changes recently introduced into the model. The computation time and memory requirements for the fourth-order model are similar to those of the present second-order GLAS model with the same 4 deg latitude, 5 deg longitude, and 9 vertical-level resolution. However, the fourth-order model forecast skill is significantly better than that of the current GLAS model, and after three days it is comparable to the 2.5 by 3 deg version of the GLAS model in the sea level pressure maps, and has less phase errors in the 500 mb maps.
A nonparametric procedure for determination of the evolution of the luminosity function of extragalactic objects and use of this for prediction of expected redshift and luminosity distribution of objects is described. The relation between this statistical evolution of the population and their physical evolution, such as the variation with cosmological epoch of their luminosity and formation rate is presented. This procedure when applied to a sample of optically selected quasars with redshifts less than two shows that the luminosity function evolves more strongly for higher luminosities, indicating a larger quasar activity at earlier epochs and a more rapid evolution of the objects during their higher luminosity phases. It is also shown that absence of many quasars at redshifts greater than three implies slowing down of this evolution in the conventional cosmological models, perhaps indicating that this is near the epoch of the birth of the quasar (and galaxies).