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.
Rainfall modification of scatterometer response from the sea surface was simulated in wind-wave tank experiments. Data show that for a given wind speed, radar cross section increases as rainfall rate increases, but this effect decreases as wind speed increases. An empirical model accounts for these observations.
A complex two-dimensional, unsteady, viscous hypersonic shock wave interaction is numerically simulated by a high-resolution, second-order fully implicit shock-capturing scheme. The physical model consists of a nonstationary oblique shock impinging on the bow shock of a blunt body. Studies indicate that the unsteady flow patterns are slightly different from their steady counterparts. However, for the sample cases investigated the peak surface pressures for the unsteady flows seem to occur at very different impingement locations than for the steady flow cases.
Electron distribution functions (EDFs) with a peak oblique to the magnetic field, adjacent to but distinct from loss-cone features, have been observed by the DE 1, Viking, and S3-3 spacecraft in passes through the nightside auroral zone, polar cap, dayside cusp, and extended dayside auroral oval. Using particle simulations, two types of wave excitation and particle acceleration mechanisms which may contribute to producing these electron conic distributions are investigated. The first involves excitation of upper-hybrid waves by the electron loss cone, and the subsequent perpendicular heating of the background and thermal electrons. The second involves excitation of downward propagating parallel modes by an auroral electron beam which frequently accompanies the upflowing electron conics. These modes provide parallel acceleration, which modifies the GDF. Those electrons which are not lost to the atmosphere and mirror back up the magnetic field line give rise to enhancements in the GDF at the edge of the loss cone.
The problem of nonlinear development of Goertler vortices and interaction with Tollmein-Schlichting (TS) waves is considered within the framework of incompressible Navier-Stokes equations which are solved by a Fourier-Chebyshev spectral method. It is shown that two-dimensional waves can be excited in the flow modulated by Goertler vortices. Due to nonlinear effects, this interaction further leads to the development of oblique waves with spanwise wavelength equal to the Goertler vortex wavelength. Interaction is also considered of oblique waves with spanwise wavelength twice that of Goertler vortices.
Rainfall modification of directional scatterometer response from the sea surface was simulated in wind-wave tank experiments. Data show that for the range of conditions in laboratory experiments, rain enhances radar cross section for all azimuthal angles relative to wind direction. This result broadens previous measurements, which showed that scatterometer response increases with increasing rainfall for radars pointing upwind. But more to the point, the data also show that the directional dynamic-range of scatterometry diminishes rapidly as rainfall rate increases. Thus, while it may be possible to determine wind speed and direction during rain, it will require adequate system sensitivity.
Recent observations indicate that a region of enhanced plasma pressure and decreased magnetic field intensity frequently occurs in front of the plasma depletion layer at the dayside magnetopause (Song et al., 1990). This inverse relationship is characteristic of a slow-mode wave. This phenomenon was simulated with a two-dimensional incompressible MHD simulation code. When a normal component of the interplanetary magnetic field is present (Bx not equal to 0), the total magnetic field intensity tends to decrease in front of the depletion layer due to the bending of the magnetic field lines, and the plasma pressure is enhanced in this region. On the other hand, when Bx = 0, this slow-mode structure is not present in the simulation and only the plasma depletion layer is observed.
A new procedure is presented for the simultaneous coarsening and refinement of three-dimensional unstructured tetrahedral meshes. This algorithm allows for localized grid adaption that is used to capture aerodynamic flow features such as vortices and shock waves in helicopter flowfield simulations. The mesh-adaption algorithm is implemented in the C programming language and uses a data structure consisting of a series of dynamically-allocated linked lists. These lists allow the mesh connectivity to be rapidly reconstructed when individual mesh points are added and/or deleted. The algorithm allows the mesh to change in an anisotropic manner in order to efficiently resolve directional flow features. The procedure has been successfully implemented on a single processor of a Cray Y-MP computer. Two sample cases are presented involving three-dimensional transonic flow. Computed results show good agreement with conventional structured-grid solutions for the Euler equations.
Over recent years, Adaptive Mesh Refinement (AMR) algorithms which dynamically match the local resolution of the computational grid to the numerical solution being sought have emerged as powerful tools for solving problems that contain disparate length and time scales. In particular, several workers have demonstrated the effectiveness of employing an adaptive, block-structured hierarchical grid system for simulations of complex shock wave phenomena. Unfortunately, from the parallel algorithm developer's viewpoint, this class of scheme is quite involved; these schemes cannot be distilled down to a small kernel upon which various parallelizing strategies may be tested. However, because of their block-structured nature such schemes are inherently parallel, so all is not lost. In this paper we describe the method by which Quirk's AMR algorithm has been parallelized. This method is built upon just a few simple message passing routines and so it may be implemented across a broad class of MIMD machines. Moreover, the method of parallelization is such that the original serial code is left virtually intact, and so we are left with just a single product to support. The importance of this fact should not be underestimated given the size and complexity of the original algorithm.
A three-dimensional cloud model, radiative transfer model-based simulation system is tested and validated against the aircraft-based radiance observations of an intense convective system in southeastern Virginia on 29 June 1986 during the Cooperative Huntsville Meteorological Experiment. NASA's ER-2, a high-altitude research aircraft with a complement of radiometers operating at 11-micrometer infrared channel and 18-, 37-, 92-, and 183-GHz microwave channels provided data for this study. The cloud model successfully simulated the cloud system with regard to aircraft- and radar-observed cloud-top heights and diameters and with regard to radar-observed reflectivity structure. For the simulation time found to correspond best with the aircraft- and radar-observed structure, brightness temperatures T(sub b) are simulated and compared with observations for all the microwave frequencies along with the 11-micrometer infrared channel. Radiance calculations at the various frequencies correspond well with the aircraft observations in the areas of deep convection. The clustering of 37-147-GHz T(sub b) observations and the isolation of the 18-GHz values over the convective cores are well simulated by the model. The radiative transfer model, in general, is able to simulate the observations reasonably well from 18 GHz through 174 GHz within all convective areas of the cloud system. When the aircraft-observed 18- and 37-GHz, and 90- and 174-GHz T(sub b) are plotted against each other, the relationships have a gradual difference in the slope due to the differences in the ice particle size in the convective and more stratiform areas of the cloud. The model is able to capture these differences observed by the aircraft. Brightness temperature-rain rate relationships compare reasonably well with the aircraft observations in terms of the slope of the relationship. The model calculations are also extended to select high-frequency channels at 220, 340, and 400 GHz to simulate the Millimeter-wave Imaging Radiometer aircraft instrument to be flown in the near future. All three of these frequencies are able to discriminate the convective and anvil portions of the system, providing useful information similar to that from the frequencies below 183 GHz but with potentially enhanced spatial resolution from a satellite platform. In thin clouds, the dominant effect of water vapor is seen at 174, 340, and 400 GHz. In thick cloudy areas, the scattering effect is dominant at 90 and 220 GHz, while the overlaying water vapor can attenuate at 174, 340, and 400 GHz. All frequencies (90-400 GHz) show strong signatures in the core.
An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Olsson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.
An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Oisson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.
This document reports the conclusion and findings of our research activities for this grant. The goal of the project is the development and application of the method of Space-Time Conservation Element and Solution Element, or the CE/SE method, to simulate chemically reacting flows. The product of this project will be a high-fidelity, time-accurate flow solver analyzing unsteady flow fields advanced propulsion concepts, including the low-emission turbojet engine combustion and flow fields of the Pulse Detonation Engines (PDE). Based on the documents and computer software of the CE/SE method that we have received from the CE/SE working group at NASA Lewis, we have focused our research effort on addressing outstanding technical issues related to the extension of the CE/SE method for unsteady, chemically reactive flows. In particular, we have made progresses in the following three aspects: (1) Derivation of the governing equations for reacting flows; (2) Numerical treatments of stiff source terms; and (3) Detailed simulations of ZND detonation waves.
A process and product providing a High Efficiency Microwave Power Amplifier (HEMPA) which propagates S-Band microwave frequency square waves; utilizing a program simulating FETs at high DC-to-RF efficiencies, which analyses linear elements of selected FETs in a frequency domain, and non-linear elements of the FETs in a time domain, and converts the time domain values into the frequency domain, and performs DC and S-parameter simulated measurements based on predefined data for each FET. Individual FET parameters are extracted and isolated by converting the S-parameters to admittance or impedance parameters to derive FET models for each FET, which the program uses to provide a final output of a HEMPA circuit based on iterative simulations of an amplification circuit utilizing microwave topology and frequencies. Iterative simulations of the amplification circuit analyze output values of a plurality of cascaded stages of the FETs, which are arranged in a push-pull configuration.
The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented with a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, under the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions that motivate this research are: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance. This paper presents an initial attempt to the flow behavior in terms (perturbation) using simulation outputs from a customer-deck model obtained from Pratt&Whitney, (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine. Flow behavior is analyzed for the high pressure compressor and combustion chamber employing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000-sample simulation. Furthermore, the characteristic waves and eave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures; which are clearly the case in the combustion chamber and high-pressure turbine. Finally a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.
The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented using a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, to determine the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions motivate this research: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance? This paper presents an initial attempt to model the flow behavior in terms of their eigenfrequencies subject to the assumption of the uncertainty or variation (perturbation). The flow behavior is explored using simulation outputs from a customer-deck model obtained from Pratt & Whitney. Variations of the main variables (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine are modeled. Flow behavior is analyzed for the high-pressure compressor and combustion chamber utilizing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high-pressure and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000 sample simulation. Furthermore, the characteristic waves and wave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures in the case of the combustion chamber and high-pressure turbine. Finally, a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.
This slide presentation reviews the development of Multiple Waveband Temperature Sensor (MWTS). The MWTS project will result in a highly stable, monolithically integrated, high resolution infrared detector array sensor that records registered thermal imagery in four infrared wavebands to infer dynamic temperature profiles on a laser-irradiated ground target. An accurate surface temperature measurement of a target in extreme environments in a non-intrusive manner is required. The development challenge is to: determine optimum wavebands (suitable for target temperatures, nature of the targets and environments) to measure accurate target surface temperature independent of the emissivity, integrate simultaneously readable multiband Quantum Well Infrared Photodetectors (QWIPs) in a single monolithic focal plane array (FPA) sensor and to integrate the hardware/software and system calibration for remote temperature measurements. The charge was therefore to develop and demonstrate a multiband infrared imaging camera with the detectors simultaneously sensitive to multiple distinct color bands for front surface temperature measurements Wavelength ( m) measurements. Amongst the requirements are: that the measurement system will not affect target dynamics or response to the laser irradiation and that the simplest criterion for spectral band selection is to choose those practically feasible spectral bands that create the most contrast between the objects or scenes of interest in the expected environmental conditions. There is in the presentation a review of the modeling and simulation of multi-wave infrared temperature measurement and also a review of the detector development and QWIP capacities.
Phased array beamforming results of the F31/A31 historical baseline counter-rotating open rotor blade set were investigated for measurement data taken on the NASA Counter-Rotating Open Rotor Propulsion Rig in the 9- by 15-Foot Low-Speed Wind Tunnel of NASA Glenn Research Center as well as data produced using the LINPROP open rotor tone noise code. The planar microphone array was positioned broadside and parallel to the axis of the open rotor, roughly 2.3 rotor diameters away. The results provide insight as to why the apparent noise sources of the blade passing frequency tones and interaction tones appear at their nominal Mach radii instead of at the actual noise sources, even if those locations are not on the blades. Contour maps corresponding to the sound fields produced by the radiating sound waves, taken from the simulations, are used to illustrate how the interaction patterns of circumferential spinning modes of rotating coherent noise sources interact with the phased array, often giving misleading results, as the apparent sources do not always show where the actual noise sources are located. This suggests that a more sophisticated source model would be required to accurately locate the sources of each tone. The results of this study also have implications with regard to the shielding of open rotor sources by airframe empennages.
Atmospheric boundary layer (ABL) turbulence causes variability of the sonic boom waveform at the ground. Recent numerical investigations of sonic boom propagation through kinematic velocity fluctuations indicate that loudness metric distributions are positively skewed relative to a normal distribution. This skewness depends on the propagation distance and turbulence intensity. Propagation simulations of N-waves and shaped booms through inhomogeneous ABL turbulence are presented. Meteorological conditions are varied to examine different daytime ABL conditions and their effect on sonic boom loudness distributions. Two outcomes are observed: 1) the loudness metric distributions become increasingly positively skewed as the propagation distance through the ABL increases, and 2) the distributions become increasingly positively skewed at the same lateral distance from the flight path as the convection level of the daytime ABL is increased. Thus, results indicate that ground level measurements of sonic boom loudness from flight tests performed at large lateral distances from the flight path may not be normally distributed, due to turbulence present in the ABL. (This research is supported by the Commercial Supersonic Technology Project of the National Aeronautics and Space Administration under Grant No. 80NSSC19K1685.)