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.
Data tapes were decoded and 24 tracks over the Indian region were reduced to common elevation. Profiles of raw scalar and vector field data and the residuals along few passes were prepared. An anomaly data set was created from the Investigator-B tape. Data was sampled on an 0.5 deg by 0.5 deg grid.
The three-dimensional boundary-layer equations in the limit as the normal coordinate tends to infinity are called the surface Euler equations. The present paper describes an accurate method for generating edge conditions for three-dimensional boundary-layer codes using these equations. The inviscid pressure distribution is first interpolated to the boundary-layer grid. The surface Euler equations are then solved with this pressure field and a prescribed set of initial and boundary conditions to yield the velocities along the two surface coordinate directions. Results for typical wing and fuselage geometries are presented. The smoothness and accuracy of the edge conditions obtained are found to be superior to the conventional interpolation procedures.
A multiblock/multigrid computation of the inviscid flow over a wing-mounted propfan transport with propwash is presented. An explicit multistage scheme drives the integral Euler equations to a steady state solution, while an actuator disk approximates the slipstream effects of the propfan blades. Practical applications of detailed surface gridding, multiple block field grids and multigrid convergence acceleration are demonstrated.
Three-dimensional numerical investigations on the nozzle start-up side load physics were performed. The objective of this study is to identify the three-dimensional side load physics and to compute the associated aerodynamic side load using an anchored computational methodology. The computational methodology is based on an unstructured-grid, and pressure-based computational fluid dynamics formulation, and a simulated inlet condition based on a system calculation. Finite-rate chemistry was used throughout the study so that combustion effect is always included, and the effect of wall cooling on side load physics is studied. The side load physics captured include the afterburning wave, transition from free- shock to restricted-shock separation, and lip Lambda shock oscillation. With the adiabatic nozzle, free-shock separation reappears after the transition from free-shock separation to restricted-shock separation, and the subsequent flow pattern of the simultaneous free-shock and restricted-shock separations creates a very asymmetric Mach disk flow. With the cooled nozzle, the more symmetric restricted-shock separation persisted throughout the start-up transient after the transition, leading to an overall lower side load than that of the adiabatic nozzle. The tepee structures corresponding to the maximum side load were addressed.
A computational model to study the stability characteristics of an evolving buoyant helium gas jet in ambient air environment is presented. Numerical formulation incorporates a segregated approach to solve for the transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum using a staggered grid method. The operating parameters correspond to the Reynolds number varying from 30 to 300 to demarcate the flow dynamics in oscillating and non-oscillating regimes. Computed velocity and concentration fields were used to analyze the flow structure in the evolving jet. For Re=300 case, results showed that an instability mode that sets in during the evolution process in Earth gravity is absent in zero gravity, signifying the importance of buoyancy. Though buoyancy initiates the instability, below a certain jet exit velocity, diffusion dominates the entrainment process to make the jet non-oscillatory as observed for the Re=30 case. Initiation of the instability was found to be dependent on the interaction of buoyancy and momentum forces along the jet shear layer.
Space communications are a fundamental part of every NASA mission, and reliable space communications are essential. Precise planning, simulation, and analysis of communications abilities are needed to know exactly what data rates, visibility times, and signal strength can be expected during the real missions. The goal of this work was to take a 3D icosphere, where each face was a point to be evaluated, and display the link performance onto a 2D projection of celestial bodies for links to a node location that could be spatially located at any of the icosahedral faces. 2D projections were produced for a wide range of spatial link grid sizes and compressed into low-size MAT files for ease of access and rapid implementation. The final icosahedral projections can be displayed onto various celestial body maps, using ultra-high-resolution maps for stitching. The results have proven satisfactory and are planned to support various simulation efforts.
This report describes modifications and enhancements to XTRAN3S, a computer program for aerodynamic, static aeroelastic, and dynamic aeroelastic analysis of three-dimensional wings in the transonic speed regime. Modifications to the program include incorporation of a viscous boundary layer, modified coefficient and generalized force integration, direct input of airfoil surface coordinates and slopes, variable dimensions of computational arrays and a modified grid mapping transformation to improve computational stability for highly swept and tapered planforms. Results obtained with the modified program are included. Modifications discussed but not incorporated include a state matrix integration method, enhanced vectorization, and use of a cartesian physical mesh.
The effects of data boundaries on objective analysis are examined. Data from limited-area networks are studied using the objective analysis technique of Barnes (1964). The impact of the data boundaries on the amplitudes of long waves and phases of short waves is described. The relation between the weight function shape parameter and the distance that boundary effects intrude into the interior of the grid is investigated. Examples displaying the influence of data field boundaries on the analyses of limited-area data sets are presented.
The interaction between a vortex ring of finite strength and an axisymmetric rotational background flow bounded by a solid wall is studied by a singular perturbation method. The analysis is carried out by combining a composite solution of a vortex ring and an unsteady Navier-Stokes solution for the background rotational flow. Using the method of averaging, numerical scheme is developed to obtain a Navier-Stokes solution in which the grid and time-step sizes depend solely on the length and velocity scales of the background flow. Numerical results are presented to illustrate the separation of the boundary layer on a solid wall and its interaction with the vortex ring.
A portable computer called BLASIM was developed at NASA LeRC to assess ice impact damage on aircraft engine blades. In addition to ice impact analyses, the code also contains static, dynamic, resonance margin, and supersonic flutter analysis capabilities. Solid, hollow, superhybrid, and composite blades are supported. An optional preprocessor (input generator) was also developed to interactively generate input for BLASIM. The blade geometry can be defined using a series of airfoils at discrete input stations or by a finite element grid. The code employs a coarse, fixed finite element mesh containing triangular plate finite elements to minimize program execution time. Ice piece is modeled using an equivalent spherical objective that has a high velocity opposite that of the aircraft and parallel to the engine axis. For local impact damage assessment, the impact load is considered as a distributed force acting over a region around the impact point. The average radial strain of the finite elements along the leading edge is used as a measure of the local damage. To estimate damage at the blade root, the impact is treated as an impulse and a combined stress failure criteria is employed. Parametric studies of local and root ice impact damage, and post-impact dynamics are discussed for solid and composite blades.
Two fourth-order-accurate finite-difference methods for numerically solving hyperbolic systems of conservation equations on smooth curvilinear grids are presented. The first method uses the differential form of the conservation equations; the second method uses the integral form of the conservation equations. Modifications to these schemes, which are required near boundaries to maintain overall high-order accuracy, are discussed. An analysis that demonstrates the stability of the modified schemes is also provided. Modifications to one of the schemes to make it total variation diminishing (TVD) are also discussed. Results that demonstrate the high-order accuracy of both schemes are included in the paper. In particular, a Ringleb-flow computation demonstrates the high-order accuracy and the stability of the boundary and near-boundary procedures. A second computation of supersonic flow over a cylinder demonstrates the shock-capturing capability of the TVD methodology. An important contribution of this paper is the dear demonstration that higher order accuracy leads to increased computational efficiency.
The objectives of the research were as follows: (1) Extend the Representative Elementary Area (RE) concept, first proposed and developed in Wood et al, (1988), to the water balance fluxes of the interstorm period (redistribution, evapotranspiration and baseflow) necessary for the analysis of long-term water balance processes. (2) Derive spatially averaged water balance model equations for spatially variable soil, topography and vegetation, over A RANGE OF CLIMATES. This is a necessary step in our goal to derive consistent hydrologic results up to GCM grid scales necessary for global climate modeling. (3) Apply the above macroscale water balance equations with remotely sensed data and begin to explore the feasibility of parameterizing the water balance constitutive equations at GCM grid scale.
Analysis of four years of earth radiation budget, cloud, and sea surface temperature data confirms that cloud parameters change dramatically when and where sea surface temperatures increase above approximately 300K. These results are based upon monthly mean values within 2.5 deg x 2.5 deg grid points over the 'warm pool' region of the western tropical Pacific. The question of whether sea surface temperatures are influenced, in turn, by the radiative effects of these clouds (Ramanathan and Collins) is less clear. Such a feedback, if it exists, is weak. The reason why clouds might have so little influence, despite large changes in their longwave and shortwave radiative effects, might be that the sea surface responds to both the longwave heating and the shortwave cooling effects of clouds, and the two effects nearly cancel. There are strong correlations between the rate of change of sea surface temperature and any of the radiation budget parameters that are highly correlated with the incident solar flux-implying that season and latitude are the critical factors determining sea surface temperatures. With the seasonal or both seasonal and latitudinal variations removed, the rate of change of sea surface temperature shows no correlation with cloud-related parameters in the western tropical Pacific.
Change detection from satellite sensor vegetation indices (VIs) presents an opportunity to monitor trends and disturbances at the regional scale for southern California's Mojave and Lower Colorado Deserts. Renewable energy sites are being constructed in this region on public lands under the Bureau of Land Management (BLM). We have developed a framework for VI change detection over the past two decades, with initial focus on three sites, Joshua Tree National Park, Mojave National Preserve, and a proximal group of Development Focus Areas (DFAs), for comparison between protected and development-targeted lands. Three Terra MODIS VIs (normalized difference [NDVI], enhanced [EVI], soil-adjusted [SAVI]) were evaluated in the Breaks for Additive Season and Trend (BFAST) setting for the regional MODIS 250-m resolution grid to estimate significant time series shifts (breakpoints) from February 2000 to May 2018. All three VIs tended to detect the maximum number of breakpoints at a grid location, but cross-correlations with precipitation and comparison with timing of wildfire burns near the study sites for breakpoint density (proportion of area with a breakpoint) showed that NDVI had the strongest response to these major disturbances, supporting its use for subsequent analysis. Time series of NDVI breakpoint change densities for individual solar energy sites did not have a consistent vegetation response following construction. Bootstrapping showed that the DFAs had significantly larger kurtosis and variance in the positive NDVI breakpoint distribution than did the protected sites, but there was no significant difference in the negative distribution for all three sites. The inconsistent post-construction NDVI signal and the large number of breakpoints overall suggested that the largest changes in vegetation cover density were tied to seasonal precipitation amounts. The distributional results indicated that existing site-specific conditions were the main control on VI responses, given the history of human disturbances in the DFAs. Although the results do not support persistent VI disturbances resulting from recent solar energy development, continued monitoring and examination of other ecological variables and surface temperatures will be vital to the long-term protection of this desert environment.
The time-dependent three-dimensional Euler equations of gas dynamics were solved numerically to study the steady compressible transonic flow about ducted propfan propulsion systems. Aerodynamic calculations were based on a four-stage Runge-Kutta time-marching finite volume solution technique with added numerical dissipation. An implicit residual smoothing operator was used to aid convergence. Two calculation grids were employed in this study. The first grid utilized an H-type mesh network with a branch cut opening to represent the axisymmetric cowl. The second grid utilized a multiple-block mesh system with a C-type grid about the cowl. The individual blocks were numerically coupled in the Euler solver. Grid systems were generated by a combined algebraic/elliptic algortihm developed specifically for ducted propfans. Numerical calculations were initially performed for unducted propfans to verify the accuracy of the three-dimensional Euler formulation. The Euler analyses were then applied for the calculation of ducted propfan flows, and predicted results were compared with experimental data for two cases. The three-dimensional Euler analyses displayed exceptional accuracy, although certain parameters were observed to be very sensitive to geometric deflections. Both solution schemes were found to be very robust and demonstrated nearly equal efficiency and accuracy, although it was observed that the multi-block C-grid formulation provided somewhat better resolution of the cowl leading edge region.
The interaction between the main passage flow through a high pressure turbine and the secondary flows due to leakage through the wheelspace rim seals is reviewed. Various experimental and analytical studies of these interactions are outlined and a numerical investigation of the interaction between the main passage flow and a wheelspace cavity seal flow is described. The numerical investigation uses a structured grid method to study the overall interaction between the turbine stage components and the wheelspace seal flow, and an unstructured grid method to resolve the detailed flow features within the geometrically complex cavity seal. The numerical results agree with various observations from experimental studies under similar flow conditions. As the flow rate through the rim cavity seal is increased, the ingestion of fluid from the main passage flow into the rim seal area decreases rapidly. A small amount of main flow gas is ingested into the rim seal area, even at high seal flow rates.
The objective of this study was to investigate compressibility effects on a high-lift flowfield by simulating the flow about a three-dimensional multi-element wing. The computations were performed by solving both the incompressible and compressible Navier-Stokes equations (using the INS3D and OVERFLOW codes) on structured, overset grids. Turbulence was modeled via the one-equation, fully turbulent Spalart-Allmaras model. The computational results were validated with surface pressure measurements acquired at the NASA Ames 7- by 10-Foot Wind Tunnel. The geometry used for all computations consisted of an unswept wing in a landing configuration with a half-span flap and a three-quarter-span slat mounted inside a rectangular duct approximating the wind tunnel walls. The solutions were carefully examined to account for effects due to differences in algorithms. Compressibility effects were demonstrated by comparing surface particle traces, sectional pressure coefficient and boundary layer profile plots. It was found that small regions of compressibility near the slat and main-element leading edge can largely impact the flow. Even small compressibility regions can have significant global effects on the circulation and separation of each of the high-lift elements.
This paper describes the spectral analysis of integrated pressures on patches of the scale model of the Space Launch System Block 1 cargo vehicle with the Unsteady Pressure-Sensitive Paint measurements, which were collected in the ascent transient aerodynamics test with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The patches are defined with x station values, indicating the position along length of the vehicle, and azimuth angles of the scale model. For each patch, the polygons are determined from the surface cells of the grid of the model, clipped with the edges of the patch, and each of the polygons is divided into triangles. The integrated pressure of the patch is determined as the ratio of the sum of the forces on the triangles over the sum of the areas of the triangles. For each run of the test, the cross power spectral density and magnitude squared coherence are computed from the time series of the integrated pressures on the patches. The results of pressure integration and spectral analysis are presented, and the data consistency of the test is demonstrated. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding for this research was provided by the NASA Aeroscience Evaluation and Test Capabilities Project.