SEARCH · Search NASA

Results for “grid disturbance”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS2015

Aerodynamic Shape Optimization of a Dual-Stream Supersonic Plug Nozzle

Aerodynamic shape optimization was performed on an isolated axisymmetric plug nozzle sized for a supersonic business jet. The dual-stream concept was tailored to attenuate nearfield pressure disturbances without compromising nozzle performance. Adjoint-based anisotropic mesh refinement was applied to resolve nearfield compression and expansion features in the baseline viscous grid. Deformed versions of the adapted grid were used for subsequent adjoint-driven shape optimization. For design, a nonlinear gradient-based optimizer was coupled to the discrete adjoint formulation of the Reynolds-averaged Navier- Stokes equations. All nozzle surfaces were parameterized using 3rd order B-spline interpolants and perturbed axisymmetrically via free-form deformation. Geometry deformations were performed using 20 design variables shared between the outer cowl, shroud and centerbody nozzle surfaces. Interior volume grid deformation during design was accomplished using linear elastic mesh morphing. The nozzle optimization was performed at a design cruise speed of Mach 1.6, assuming core and bypass pressure ratios of 6.19 and 3.24, respectively. Ambient flight conditions at design were commensurate with 45,000-ft standard day atmosphere.

shape optimization↗

NASA NTRS1991

Prediction of unsteady aerodynamic loads in cascades using the linearized Euler equations on deforming grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. Using the linearized Euler technique, one decomposes the flow into a mean (or steady) flow plus an unsteady, harmonically varying, small disturbance flow. Linear variable coefficient equations describe the small disturbance behavior of the flow, and are solved using a pseudo-time marching Lax-Wendroff scheme. For the blade motion problem, a harmonically deforming computational rid that conforms to the motion of vibrating blades eliminates large error producing mean flow gradient terms that would otherwise appear in the unsteady flow tangency boundary condition. The paper also presents a new, numerically exact, nonreflecting far-field boundary condition based on an eigenanalysis of the discretized equations. Computed flow solutions demonstrate the computational accuracy and efficiency of the present method. The solution of the linearized Euler equations requires one to two orders of magnitude less computer time than solution of the nonlinear Euler equations using traditional time-accurate time-marching techniques. In addition, the deformable grid significantly improves the accuracy of the solution.

Hall, Kenneth C.↗

NASA NTRS1976

Application of a multi-level grid method to transonic flow calculations

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse to fine. The coarser grids are used to diminish the magnitude of the smooth part of the residuals. The method was applied to the solution of the transonic small disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is studied with meshes of both constant and variable step size.

South, J. C., Jr.↗

NASA NTRS1981

Use of MAGSAT anomaly data for crustal structure and mineral resources in the US Midcontinent

The analysis and preliminary interpretation of investigator-B MAGSAT data are addressed. The data processing included: (1) removal of spurious data points; (2) statistical smoothing along individual data tracks, to reduce the effect of geomagnetic transient disturbances; (3) comparison of data profiles spatially coincident in track location but acquired at different times; (4) reduction of data by weighted averaging to a grid with 1 deg xl deg latitude/longitude spacing, and with elevations interpolated and weighted to a common datum of 400 km; (5) wavelength filtering; and (6) reduction of the anomaly map to the magnetic pole. Agreement was found between a magnitude data anomaly map and a reduce-to-the-pole map supporting the general assumption that, on a large scale (long wavelength), it is induced crustal magnetization which is responsible for major anamalies. Anomalous features are identified and explanations are suggested with regard to crustal structure, petrologic characteristics, and Curie temperature isotherms.

Carmichael, R. S.↗

NASA NTRS1982

Cloud fluctuation statistics

A space-time statistical analysis of total outgoing infrared radiation is used to determine the gross features of day-to-day cloudiness fluctuations over the Pacific Ocean in summer and winter. IR fluctuations arise from the passage of cloudiness systems through a grid box as well as the creation and destruction of cloudiness in the box. Which process dominates depends upon the size of the box relative to the size, speed and persistence time of a typical cloudiness system. In most regions the statistical analysis yields advection speeds characteristic of 700 mb mean flow with spatial dependence resembling the 300 mb mean flow. Spatial scales less than 2000 km predominate, smaller scales having less persistence. Characteristic time scales are on the order of one or two days, even for a grid box spanning the entire North Pacific storm track. This result is remarkable in view of the much longer time scales commonly associated with atmospheric disturbances. Apparently many cloudiness systems are created and destroyed during the lifetime of a single disturbance.

Cahalan, R. F.↗

NASA NTRS1990

Two-dimensional Euler and Navier-Stokes Time accurate simulations of fan rotor flows

Two numerical methods are presented which describe the unsteady flow field in the blade-to-blade plane of an axial fan rotor. These methods solve the compressible, time-dependent, Euler and the compressible, turbulent, time-dependent, Navier-Stokes conservation equations for mass, momentum, and energy. The Navier-Stokes equations are written in Favre-averaged form and are closed with an approximate two-equation turbulence model with low Reynolds number and compressibility effects included. The unsteady aerodynamic component is obtained by superposing inflow or outflow unsteadiness to the steady conditions through time-dependent boundary conditions. The integration in space is performed by using a finite volume scheme, and the integration in time is performed by using k-stage Runge-Kutta schemes, k = 2,5. The numerical integration algorithm allows the reduction of the computational cost of an unsteady simulation involving high frequency disturbances in both CPU time and memory requirements. Less than 200 sec of CPU time are required to advance the Euler equations in a computational grid made up of about 2000 grid during 10,000 time steps on a CRAY Y-MP computer, with a required memory of less than 0.3 megawords.

Boretti, A. A.↗

NASA NTRS2002

Development of Computational Aeroacoustics Code for Jet Noise and Flow Prediction

Accurate prediction of jet fan and exhaust plume flow and noise generation and propagation is very important in developing advanced aircraft engines that will pass current and future noise regulations. In jet fan flows as well as exhaust plumes, two major sources of noise are present: large-scale, coherent instabilities and small-scale turbulent eddies. In previous work for the NASA Glenn Research Center, three strategies have been explored in an effort to computationally predict the noise radiation from supersonic jet exhaust plumes. In order from the least expensive computationally to the most expensive computationally, these are: 1) Linearized Euler equations (LEE). 2) Very Large Eddy Simulations (VLES). 3) Large Eddy Simulations (LES). The first method solves the linearized Euler equations (LEE). These equations are obtained by linearizing about a given mean flow and the neglecting viscous effects. In this way, the noise from large-scale instabilities can be found for a given mean flow. The linearized Euler equations are computationally inexpensive, and have produced good noise results for supersonic jets where the large-scale instability noise dominates, as well as for the tone noise from a jet engine blade row. However, these linear equations do not predict the absolute magnitude of the noise; instead, only the relative magnitude is predicted. Also, the predicted disturbances do not modify the mean flow, removing a physical mechanism by which the amplitude of the disturbance may be controlled. Recent research for isolated airfoils' indicates that this may not affect the solution greatly at low frequencies. The second method addresses some of the concerns raised by the LEE method. In this approach, called Very Large Eddy Simulation (VLES), the unsteady Reynolds averaged Navier-Stokes equations are solved directly using a high-accuracy computational aeroacoustics numerical scheme. With the addition of a two-equation turbulence model and the use of a relatively coarse grid, the numerical solution is effectively filtered into a directly calculated mean flow with the small-scale turbulence being modeled, and an unsteady large-scale component that is also being directly calculated. In this way, the unsteady disturbances are calculated in a nonlinear way, with a direct effect on the mean flow. This method is not as fast as the LEE approach, but does have many advantages to recommend it; however, like the LEE approach, only the effect of the largest unsteady structures will be captured. An initial calculation was performed on a supersonic jet exhaust plume, with promising results, but the calculation was hampered by the explicit time marching scheme that was employed. This explicit scheme required a very small time step to resolve the nozzle boundary layer, which caused a long run time. Current work is focused on testing a lower-order implicit time marching method to combat this problem.

Keith, Theo G., Jr.↗

NASA NTRS2009

Conical-Domain Model for Estimating GPS Ionospheric Delays

The conical-domain model is a computational model, now undergoing development, for estimating ionospheric delays of Global Positioning System (GPS) signals. Relative to the standard ionospheric delay model described below, the conical-domain model offers improved accuracy. In the absence of selective availability, the ionosphere is the largest source of error for single-frequency users of GPS. Because ionospheric signal delays contribute to errors in GPS position and time measurements, satellite-based augmentation systems (SBASs) have been designed to estimate these delays and broadcast corrections. Several national and international SBASs are currently in various stages of development to enhance the integrity and accuracy of GPS measurements for airline navigation. In the Wide Area Augmentation System (WAAS) of the United States, slant ionospheric delay errors and confidence bounds are derived from estimates of vertical ionospheric delay modeled on a grid at regularly spaced intervals of latitude and longitude. The estimate of vertical delay at each ionospheric grid point (IGP) is calculated from a planar fit of neighboring slant delay measurements, projected to vertical using a standard, thin-shell model of the ionosphere. Interpolation on the WAAS grid enables estimation of the vertical delay at the ionospheric pierce point (IPP) corresponding to any arbitrary measurement of a user. (The IPP of a given user s measurement is the point where the GPS signal ray path intersects a reference ionospheric height.) The product of the interpolated value and the user s thin-shell obliquity factor provides an estimate of the user s ionospheric slant delay. Two types of error that restrict the accuracy of the thin-shell model are absent in the conical domain model: (1) error due to the implicit assumption that the electron density is independent of the azimuthal angle at the IPP and (2) error arising from the slant-to-vertical conversion. At low latitudes or at mid-latitudes under disturbed conditions, the accuracy of SBAS systems based upon the thin-shell model suffers due to the presence of complex ionospheric structure, high delay values, and large electron density gradients. Interpolation on the vertical delay grid serves as an additional source of delay error. The conical-domain model permits direct computation of the user s slant delay estimate without the intervening use of a vertical delay grid. The key is to restrict each fit of GPS measurements to a spatial domain encompassing signals from only one satellite. The conical domain model is so named because each fit involves a group of GPS receivers that all receive signals from the same GPS satellite (see figure); the receiver and satellite positions define a cone, the satellite position being the vertex. A user within a given cone evaluates the delay to the satellite directly, using (1) the IPP coordinates of the line of sight to the satellite and (2) broadcast fit parameters associated with the cone. The conical-domain model partly resembles the thin-shell model in that both models reduce an inherently four-dimensional problem to two dimensions. However, unlike the thin-shell model, the conical domain model does not involve any potentially erroneous simplifying assumptions about the structure of the ionosphere. In the conical domain model, the initially four-dimensional problem becomes truly two-dimensional in the sense that once a satellite location has been specified, any signal path emanating from a satellite can be identified by only two coordinates; for example, the IPP coordinates. As a consequence, a user s slant-delay estimate converges to the correct value in the limit that the receivers converge to the user s location (or, equivalently, in the limit that the measurement IPPs converge to the user s IPP).

Sparks, Lawrence↗

NASA NTRS2018

Vegetation Change Detection in Southern California Solar Energy Developments

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.

Potter, Christopher S.↗

NASA NTRS1989

An assessment of the structural dynamic effects on the microgravity environment of a reference Space Station

An interim 'Permanently Manned Capability' Space Station configuration and one designated 'Assembly Complete' are modeled by FEM techniques in order to select forcing functions for modal transient response analysis and illustrate sample waveforms. In all, 114 applied-force cases have been executed to simulate such expected operational disturbances as crewmembers' treadmill exercises and EVAs. The present discussion of the results obtained give attention to the acceleration-response environment of the U.S. Laboratory's center grid-point. The magnitude of the acceleration responses obtained indicate that attenuation systems must be incorporated, or operational constraints must be instituted.

Del Basso, Steve↗

NASA NTRS1985

Unsteady transonic flow calculations for interfering lifting surface configurations

Unsteady transonic flow calculations are presented for aerodynamically interfering lifting surface configurations. Calculations are performed by extending the XTRAN3S (Version 1.5) unsteady transonic small-disturbance code to allow the treatment of an additional lifting surface. The research was conducted as a first-step toward developing the capability to treat a complete flight vehicle. Grid generation procedures for swept tapered interfering lifting surface applications of XTRAN3S are described. Transonic calculations are presented for wing-tail and canard-wing configurations for several values of mean angle of attack. The effects of aerodynamic interference on transonic steady pressure distributions and steady and oscillatory spanwise lift distributions are demonstrated. Results due to wing, tail, or canard pitching motions are presented and discussed in detail.

Batina, J. T.↗

NASA NTRS1985

Unsteady transonic flow calculations for interfering lifting surface configurations

Batina, J. T.↗

NASA NTRS2008

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave transitionally in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗

NASA NTRS2008

Turbulence Model Behavior in Low Reynolds Number Regions of Aerodynamic Flowfields

The behaviors of the widely-used Spalart-Allmaras (SA) and Menter shear-stress transport (SST) turbulence models at low Reynolds numbers and under conditions conducive to relaminarization are documented. The flows used in the investigation include 2-D zero pressure gradient flow over a flat plate from subsonic to hypersonic Mach numbers, 2-D airfoil flow from subsonic to supersonic Mach numbers, 2-D subsonic sink-flow, and 3-D subsonic flow over an infinite swept wing (particularly its leading-edge region). Both models exhibit a range over which they behave 'transitionally' in the sense that the flow is neither laminar nor fully turbulent, but these behaviors are different: the SST model typically has a well-defined transition location, whereas the SA model does not. Both models are predisposed to delayed activation of turbulence with increasing freestream Mach number. Also, both models can be made to achieve earlier activation of turbulence by increasing their freestream levels, but too high a level can disturb the turbulent solution behavior. The technique of maintaining freestream levels of turbulence without decay in the SST model, introduced elsewhere, is shown here to be useful in reducing grid-dependence of the model's transitional behavior. Both models are demonstrated to be incapable of predicting relaminarization; eddy viscosities remain weakly turbulent in accelerating or laterally-strained boundary layers for which experiment and direct simulations indicate turbulence suppression. The main conclusion is that these models are intended for fully turbulent high Reynolds number computations, and using them for transitional (e.g., low Reynolds number) or relaminarizing flows is not appropriate.

Rumsey, Christopher L.↗

NASA NTRS2011

Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods

A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.

Atkins, Harold L.↗

NASA NTRS1982

Transonic perturbation analysis of wing-fuselage-nacelle-pylon configurations with powered jet exhausts

A transonic small perturbation method has been developed for the analysis of general wing-fuselage-nacelle-pylon configurations with powered jet exhausts. Finite difference successive line relaxation algorithm is used to solve the small disturbance potential equation in conservative form. The nacelle tangency condition and the jet exhaust plume contact conditions are fulfilled in a quasi-cylindrical fashion on a surface fitting the Cartesian grid. The pylon tangency condition is treated in a quasi-planar manner as for the wing. Viscous displacement effects on the wing are modeled by suitable shape changes including the placement of a viscous ramp at the base of the shock. Computed results of a transport configuration show satisfactory correlation with test data.

Wai, J. C.↗

NASA NTRS1993

Flow transition with 2-D roughness elements in a 3-D channel

We develop a new numerical approach to study the spatially evolving instability of the streamwise dominant flow in the presence of roughness elements. The difficulty in handling the flow over the boundary surface with general geometry is removed by using a new conservative form of the governing equations and an analytical mapping. The numerical scheme uses second-order backward Euler in time, fourth-order central differences in all three spatial directions, and boundary-fitted staggered grids. A three-dimensional channel with multiple two-dimensional-type roughness elements is employed as the test case. Fourier analysis is used to decompose different Fourier modes of the disturbance. The results show that surface roughness leads to transition at lower Reynolds number than for smooth channels.

Liu, Zhining↗

NASA NTRS1984

Transonic analysis of canted winglets

A computational method developed to provide a transonic analysis for upper/lower surface wing-tip mounted winglets is described. Winglets with arbitrary planform, cant and toe angle, and airfoil section can be modeled. The embedded grid approach provides high flow field resolution and the required geometric flexibility. In particular, coupled Cartesian/cylindrical grid systems are used to model the complex geometry presented by canted upper/lower surface winglets. A new rotated difference scheme is introduced in order to maintain the stability of the small-disturbance formulation in the presence of large spanwise velocities. Wing and winglet viscous effects are modeled using a two-dimensional 'strip' boundary layer analysis. Correlations with wind tunnel and flight test data for three transport configurations are included.

Rosen, B. S.↗