SEARCH · Search NASA

Results for “Grid Stability”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Euler/Navier-Stokes Flow Computations on Flexible Configurations for Stability Analysis

Abstract Longitudinal dynamic stability derivatives required for design of aircraft are computed by using the state-of-the-art numerical methods for wing-body configurations. The flow is modeled using the Euler/Navier-Stokes equations with turbulence models and solved using an efficient finite-difference scheme suitable for patched structured grids. Computations are made at a flow regime that is beyond the limits of the current linear methods mostly used for computing stability derivatives. Flow conditions include shock-waves and viscous dominated vortical flows. Effect of Mach number and angle-of-attack on stability derivatives are demonstrated for a typical wing-body configuration. For the same configuration the effects of wing flexibility on the magnitude and phase angles of stability derivatives are also demonstrated.

Guruswamy, Guru P.↗

NASA NTRS2000

High Resolution Experiments of Compressible Turbulence Interacting With Shock Waves

Most of the second year of our research program focused on exploring the potentially favorable the effects of expansion waves on homogeneous and isotropic turbulence, which is formed downstream of a grid. Expansion waves are associated with compressible flows and may reduce the drag over airfoils by suppressing turbulence. In the very few previous investigations of interactions of turbulence with expansion waves the effects due to stabilizing streamline curvature substantially masked the effects of turbulence suppression due to flow expansion though the waves. In the present flow configuration planar expansion waves interact with grid generated turbulence in our high-resolution shock tube research facility. This approach will assess directly the effects of the interaction on turbulence. The first objective of our study was to identify the nature of expansion waves present in our shock tube facility. Our time-dependent numerical simulations of the flow in our facility indicated the existence of two regions of traveling expansion waves. The system of expansion waves utilized in this investigation is generated by the exiting shock wave and the induced flow behind it at the end of the driver. Several new measuring techniques are being developed which are capable of providing velocity-gradient-related quantities in compressible flows for the first time.

Agui, Juan↗

NASA NTRS1993

Time-accurate Navier-Stokes calculations with multigrid acceleration

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is implemented by modifying the multigrid-multiblock version of the steady Navier-Stokes equations solver, TLNS3D. The scheme is fully implicit in time and uses TLNS3D to iteratively invert the equations at each physical time step. The design objective of the scheme is unconditional stability (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability which is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. An analysis of the iterative procedure and the implementation of this procedure in TLNS3D are discussed. Numerical results are presented to show both the capabilities of the scheme and its speed up relative to the use of global minimum time stepping. Reductions in computational times of an order of magnitude are demonstrated.

Melson, N. Duane↗

NASA NTRS1996

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗

NASA NTRS1998

Evaluation of Boundary Conditions for the Gust-Cascade Problem

Using a high-order accuracy finite-difference time-domain algorithm, the acoustic scattering from a flat-plate cascade is computed. Keeping the grid and time step fixed, the effect of four different boundary conditions on the accuracy and stability of the computed solution is compared.

Hixon, R.↗

NASA NTRS2017

Stabilized Finite Elements in FUN3D

A Streamlined Upwind Petrov-Galerkin (SUPG) stabilized finite-element discretization has been implemented as a library into the FUN3D unstructured-grid flow solver. Motivation for the selection of this methodology is given, details of the implementation are provided, and the discretization for the interior scheme is verified for linear and quadratic elements by using the method of manufactured solutions. A methodology is also described for capturing shocks, and simulation results are compared to the finite-volume formulation that is currently the primary method employed for routine engineering applications. The finite-element methodology is demonstrated to be more accurate than the finite-volume technology, particularly on tetrahedral meshes where the solutions obtained using the finite-volume scheme can suffer from adverse effects caused by bias in the grid. Although no effort has been made to date to optimize computational efficiency, the finite-element scheme is competitive with the finite-volume scheme in terms of computer time to reach convergence.

Anderson, W. Kyle↗

NASA NTRS2018

Global Patterns of Crop Yield Stability Under Additional Nutrient and Water Inputs

Agricultural production must increase to feed a growing and wealthier population, as well as to satisfy increasing demands for biomaterials and biomass-based energy. At the same time, deforestation and land-use change need to be minimized in order to preserve biodiversity and maintain carbon stores in vegetation and soils. Consequently, agricultural land use needs to be intensified in order to increase food production per unit area of land. Here we use simulations of AgMIP's Global Gridded Crop Model Intercomparison (GGCMI) phase 1 to assess implications of input-driven intensification (water, nutrients) on crop yield and yield stability, which is an important aspect in food security. We find region- and crop-specific responses for the simulated period 1980+/-2009 with broadly increasing yield variability under additional nitrogen inputs and stabilizing yields under additional water inputs (irrigation), reflecting current patterns of water and nutrient limitation. The different models of the GGCMI ensemble show similar response patterns, but model differences warrant further research on management assumptions, such as variety selection and soil management, and inputs as well as on model implementation of different soil and plant processes, such as on heat stress, and parameters. Higher variability in crop productivity under higher fertilizer input will require adequate buffer mechanisms in trade and distribution/storage networks to avoid food price volatility.

Maize↗

NASA NTRS1994

Effects of Artificial Viscosity on the Accuracy of High-reynolds-number Kappa-epsilon Turbulence Model

Wall functions, as used in the typical high Reynolds number k-epsilon turbulence model, can be implemented in various ways. A least disruptive method (to the flow solver) is to directly solve for the flow variables at the grid point next to the wall while prescribing the values of k and epsilon. For the centrally-differenced finite-difference scheme employing artificial viscocity (AV) as a stabilizing mechanism, this methodology proved to be totally useless. This is because the AV gives rise to a large error at the wall due to too steep a velocity gradient resulting from the use of a coarse grid as required by the wall function methodology. This error can be eliminated simply by extrapolating velocities at the wall, instead of using the physical values of the no-slip velocities (i.e. the zero value). The applicability of the technique used in this paper is demonstrated by solving a flow over a flat plate and comparing the results with those of experiments. It was also observed that AV gives rise to a velocity overshoot (about 1 percent) near the edge of the boundary layer. This small velocity error, however, can yield as much as 10 percent error in the momentum thickness. A method which integrates the boundary layer up to only the edge of the boundary (instead of infinity) was proposed and demonstrated to give better results than the standard method.

Chitsomboon, Tawit↗

NASA NTRS1997

On Approximate Factorization Schemes for Solving the Full Potential Equation

An approximate factorization scheme based on the AF2 algorithm is presented for solving the three-dimensional full potential equation for the transonic flow about isolated wings. Two spatial discretization variations are presented, one using a hybrid first-order/second-order-accurate scheme and the second using a fully second-order-accurate scheme. The present algorithm utilizes a C-H grid topology to map the flow field about the wing. One version of the AF2 iteration scheme is used on the upper wing surface and another slightly modified version is used on the lower surface. These two algorithm variations are then connected at the wing leading edge using a local iteration technique. The resulting scheme has improved linear stability characteristics and improved time-like damping characteristics relative to previous implementations of the AF2 algorithm. The presentation is highlighted with a grid refinement study and a number of numerical results.

Holst, Terry L.↗

NASA NTRS1991

Hi-alpha forebody design. Part 2: Determination of body shapes for positive directional stability

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at low speed high angle-of-attack conditions with sideslip. The purpose is to define forebody geometries which provide good directional stability characteristics under these conditions. The flows of the F-5A forebody and Erickson forebody were recomputed with better and refined grids. The results were obtained using a modified version of cfl3d to solve either the Euler equations or the Reynolds equations employing a form of the Baldwin-Lomax turbulence model. Based on those results, we conclude that current CFD methods can be used to investigate the aerodynamic characteristics of forebodies to achieve desirable high angle-of-attack characteristics. An analytically defined generic forebody model is described, and a systematic study of forebody shapes was then conducted to determine which shapes promote a positive contribution to directional stability at high angle-of-attack. A novel way of presenting the results is used to illustrate how the positive contribution arises. Based on the results of this initial parametric study, some guidelines for aerodynamic design to promote positive directional stability are presented.

Ravi, R.↗

NASA NTRS2002

Flutter Analysis of a Transonic Fan

This paper describes the calculation of flutter stability characteristics for a transonic forward swept fan configuration using a viscous aeroelastic analysis program. Unsteady Navier-Stokes equations are solved on a dynamically deforming, body fitted, grid to obtain the aeroelastic characteristics using the energy exchange method. The non-zero inter-blade phase angle is modeled using phase-lagged boundary conditions. Results obtained show good correlation with measurements. It is found that the location of shock and variation of shock strength strongly influenced stability. Also, outboard stations primarily contributed to stability characteristics. Results demonstrate that changes in blade shape impact the calculated aerodynamic damping, indicating importance of using accurate blade operating shape under centrifugal and steady aerodynamic loading for flutter prediction. It was found that the calculated aerodynamic damping was relatively insensitive to variation in natural frequency.

Srivastava, R.↗

NASA NTRS1996

Precision Lower Stratospheric Temperature Monitoring with the MSU: Technique, Validation, and Results, 1979-1991

Microwave Sounding Unit channel 4 data from the TIROS-N (Television and Infrared Operational Satellite)-N series of NOAA (National Oceanic and Atmospheric Administration) satellites are intercalibrated to provide a continuous global record of deep-layer averaged lower stratospheric temperatures during 1979-1991. A 13-year record of temperature anomalies is time averaged into pentads and months on a 2.5 deg. grid. The monthly gridpoint anomalies are validated with ten years of radiosonde data during 1979-88. The calibration stability of each satellite's measurements is evaluated during satellite overlap periods, the longest of which reveal no measurable instrumental drift at the level of 0.01 C yr (exp -1). Intercomparisons between NOAA-6 and NOAA- 7 anomalies indicate monthly gridpoint precision of 0.05 C in the tropics to around 0.10 C in the extratropics, and signal-to-noise ratios precision and stability statistics are much better than have been previously reported by other investigators for MSU channel 4. Pentad precision is about 0.10 C in the tropics to around 0.25 C at high latitudes and signal-to- noise ratios generally over 250 in the tropics and high latitudes, but 100-200 in the middle latitudes. Radiosonde comparisons to the monthly gridpoint anomalies have correlations ranging from 0.90 in the tropics (where the interannual variability is smallest) to as high as 0.99 at high-latitude stations. The corresponding standard error of estimate is generally around 0.3 C. A significant difference in decadal trends is found between the satellite and radiosonde systems, with a step change of 0.21 C (sondes cooler) compared to the satellite measurements. Investigations of the possible sources of the discrepancy lead us to suspect that the gradual transition from on-site calibration of sondes with thermometers to factory calibration of sondes around 1982 might have caused a change in the calibration, although this conclusion must be viewed as tentative. The largest globally averaged temperature variations during 1979-91 occur after the El Chichon (1982) and Pinatubo (1991) volcanic eruptions. These warm events are superimposed upon a net downward trend in temperatures during the period. This cooling trend has more of a step function than linear character, with the step occurring during the El Chichon warm event. It is strongest in polar regions and the Northern Hemisphere middle latitudes. These characteristics are qualitatively consistent with radiative adjustments expected to occur with observed ozone depletions.

Spencer, Roy W.↗

NASA NTRS1988

The effect of perspective displays on altitude and stability control in simulated rotary wing flight

The effect of perspective displays on flight performance is investigated using two simulation experiments. In the first, a perspective grid display was superimposed on computer-generated terrain and subjects attempted to maintain their initial attitude in a simulated hover using terrain and/or one of four grid patterns. Horizontal lines produced the best attitude control performance. In the second experiment, a square grid was studied in combination with various visual display configurations and grid attachment conditions. It was found that performance with the panel-mounted display was significantly worse than with the out-the-window or helmet-mounted displays. The results suggested that the partial grid attachment condition improved hovering performance with the panel-mounted display.

O'Donnell, K. A.↗

NASA NTRS2009

Addition of Improved Shock-Capturing Schemes to OVERFLOW 2.1

Existing approximate Riemann solvers do not perform well when the grid is not aligned with strong shocks in the flow field. Three new approximate Riemann algorithms are investigated to improve solution accuracy and stability in the vicinity of strong shocks. The new algorithms are compared to the existing upwind algorithms in OVERFLOW 2.1. The new algorithms use a multidimensional pressure gradient based switch to transition to a more numerically dissipative algorithm in the vicinity of strong shocks. One new algorithm also attempts to artificially thicken captured shocks in order to alleviate the errors in the solution introduced by "stair-stepping" of the shock resulting from the approximate Riemann solver. This algorithm performed well for all the example cases and produced results that were almost insensitive to the alignment of the grid and the shock.

Burning, Pieter G.↗

NASA NTRS1982

A NASTRAN model of a large flexible swing-wing bomber. Volume 2: NASTRAN model development-horizontal stabilzer, vertical stabilizer and nacelle structures

The NASTRAN model plans for the horizontal stabilizer, vertical stabilizer, and nacelle structure were expanded in detail to generate the NASTRAN model for each of these substructures. The grid point coordinates were coded for each element. The material properties and sizing data for each element were specified. Each substructure model was thoroughly checked out for continuity, connectivity, and constraints. These substructures were processed for structural influence coefficients (SIC) point loadings and the deflections were compared to those computed for the aircraft detail models. Finally, a demonstration and validation processing of these substructures was accomplished using the NASTRAN finite element program installed at NASA/DFRC facility.

Mock, W. D.↗

NASA NTRS1986

Analysis of coaxial spray combustion flames and related numerical issues

An approach to the simulation of strongly coupled multiphase flows in combustion hardware is sketched and its unique requirements highlighted. An example of a successful application to a coaxial injector flame is presented. Furthermore, several numerical issues that tend to interact with the physics of the problem are discussed with special regard to their potential impact on the choices of numerical parameters by the analyst. These include the issues of stability, numerical diffusivity, stiffness, and boundary conditions. The theme of this paper focuses on the intriguing relationships among the grid, the solution algorithm, and the actual physical mechanisms themselves.

Liang, P. Y.↗

NASA NTRS1987

Improving the geological interpretation of magnetic and gravity satellite anomalies

Quantitative analysis of the geologic component of observed satellite magnetic and gravity fields requires accurate isolation of the geologic component of the observations, theoretically sound and viable inversion techniques, and integration of collateral, constraining geologic and geophysical data. A number of significant contributions were made which make quantitative analysis more accurate. These include procedures for: screening and processing orbital data for lithospheric signals based on signal repeatability and wavelength analysis; producing accurate gridded anomaly values at constant elevations from the orbital data by three-dimensional least squares collocation; increasing the stability of equivalent point source inversion and criteria for the selection of the optimum damping parameter; enhancing inversion techniques through an iterative procedure based on the superposition theorem of potential fields; and modeling efficiently regional-scale lithospheric sources of satellite magnetic anomalies. In addition, these techniques were utilized to investigate regional anomaly sources of North and South America and India and to provide constraints to continental reconstruction. Since the inception of this research study, eleven papers were presented with associated published abstracts, three theses were completed, four papers were published or accepted for publication, and an additional manuscript was submitted for publication.

Hinze, William J.↗

NASA NTRS1997

Evaluating and Understanding Parameterized Convective Processes and their Role in the Development of Mesoscale Precipitation Systems

Research efforts during the second year have centered on improving the manner in which convective stabilization is achieved in the Penn State/NCAR mesoscale model MM5. Ways of improving this stabilization have been investigated by (1) refining the partitioning between the Kain-Fritsch convective parameterization scheme and the grid scale by introducing a form of moist convective adjustment; (2) using radar data to define locations of subgrid-scale convection during a dynamic initialization period; and (3) parameterizing deep-convective feedbacks as subgrid-scale sources and sinks of mass. These investigations were conducted by simulating a long-lived convectively-generated mesoscale vortex that occurred during 14-18 Jul. 1982 and the 10-11 Jun. 1985 squall line that occurred over the Kansas-Oklahoma region during the PRE-STORM experiment. The long-lived vortex tracked across the central Plains states and was responsible for multiple convective outbreaks during its lifetime.

Fritsch, J. Michael↗