SEARCH · Search NASA

Results for “Upwinding”

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 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.↗

NASA NTRS1986

Three-dimensional AOTV flowfields in chemical nonequilibrium

A technique for upwind differencing of the three-dimensional species continuity equations is presented which permits computation of steady flows in chemical equilibrium and nonequilibrium. The capabilities and shortcomings of the present approach for equilibrium and nonequilibrium flows is discussed. Modifications now being investigated to improve computational time are outlined.

Gnoffo, P. A.↗

NASA NTRS1985

Implicit finite-difference methods for the Euler equations

The present paper is concerned with two-dimensional Euler equations and with schemes which are in use of the time of this writing. Most of the development presented carries over directly to three dimensions. The characteristics of the two-dimensional Euler equations in Cartesian coordinates are considered along with generalized curvilinear coordinate transformations, metric relations, invariants of the transformation, flux Jacobian matrices and eigensystems, numerical algorithms, flux split algorithms, implicit and explicit nonlinear control (smoothing), upwind differencing in supersonic regions, unsteady and steady-state computation, the diagonal form of implicit algorithm, metric differencing and invariants, boundary conditions, geometry and mesh generation, and sample solutions.

Pulliam, T. H.↗

NASA NTRS1986

A numerical study of shock wave diffraction by a circular cylinder

The nonstationary shock wave diffraction patterns generated by a blast wave impinging on a circular cylinder are numerically simulated using a second-order hybrid upwind method for solving the two-dimensional inviscid compressible Euler equations of gasdynamics. The complete diffraction patterns, including the transition from regular to Mach reflection, trajectory of the Mach triple point and the complex shock-on-shock interaction at the wake region resulting from the Mach shocks collision behind the cylinder are reported in detail. Pressure-time history and various contour plots are also included. Comparison between the work of Bryson and Gross (1961) which included both experimental schlieren pictures and theoretical calculations using Whitham's ray-shock theory and results of the present finite difference computation indicate good agreement in every aspect except for some nonideal gas and viscous effects which are not accounted for by the Euler equations.

Yang, J.-Y.↗

NASA NTRS1986

Measurements of ocean wave spectra and modulation transfer function with the airborne two-frequency scatterometer

The directional spectrum and the microwave modulation transfer function of ocean waves can be measured with the airborne two frequency scatterometer technique. Similar to tower based observations, the aircraft measurements of the Modulation Transfer Function (MTF) show that it is strongly affected by both wind speed and sea state. Also detected are small differences in the magnitudes of the MTF between downwind and upwind radar look directions, and variations with ocean wavenumber. The MTF inferred from the two frequency radar is larger than that measured using single frequency, wave orbital velocity techniques such as tower based radars or ROWS measurements from low altitude aircraft. Possible reasons for this are discussed. The ability to measure the ocean directional spectrum with the two frequency scatterometer, with supporting MTF data, is demonstrated.

Weissman, D. E.↗

NASA NTRS1986

Discrete models for the numerical analysis of time-dependent multidimensional gas dynamics

A possible technique is explored for extending to multidimensional flows some of the upwind-differencing methods that are highly successful in the one-dimensional case. Emphasis is on the two-dimensional case, and the flow domain is assumed to be divided into polygonal computational elements. Inside each element, the flow is represented by a local superposition of elementary solutions consisting of plane waves not necessarily aligned with the element boundaries.

Roe, P. L.↗

NASA NTRS1986

Implicit, vectorizable schemes for the flux-difference split, three-dimensional Navier-Stokes equations

Two hybrid upwind models are defined for solving the Euler equations. The algorithms both employ approximate factorization (AF) in crossplane and symmetric block Gauss-Seidel relaxation in the third direction. One approach adds an additional factorization step to lower the number of required grid point operations for inversion of the block tridiagonal matrices; however, the move permits only one third of the operations to be vectorized. Finite difference solutions are calculated on a C-H-type grid, in this case enveloping a slender, sharp-edged delta wing. Sample data are provided for the calculated vortex flow for Re of 10,000, at a 20.5 deg angle of attack, represented in a crossflow velocity vector plot and in a spanwise pressure coefficient distribution. The AF scheme, without additional factorization, when used with a grid covering 51 x 51 x 72 points provides a convergent solution with no time step lasting longer than 0.00001 sec.

Liu, C. H.↗

NASA NTRS1985

AOTV bluff body flow - Relaxation algorithm

Preliminary investigations have led to the consideration of two classes of energy efficient aeroassisted orbital transfer vehicle (AOTV) concepts, with high and low L/D. The AOTV represent reentry space vehicles which are to fly hypersonically in the continuum upper reaches of the atmosphere. Problems arise in connection with convective and nonequilibrium radiative heating. As the possibilities for obtaining needed information by experimentation are limited, reliance must be placed on computational techniques to simulate testing of design concepts. A new relaxation algorithm for the compressible Navier-Stokes equations has been developed in the context of recent work with a conservative flux difference eigenvector split implicit upwind scheme CSCM (conservative supracharacteristics method). Attention is given to the CSCM implicit method, and early results for axisymmetric bluff body flow including complete forebody and base flow for a free flight range model drag brake AOTV.

Lombard, C. K.↗

NASA NTRS1986

Meso-beta scale perturbations of the wind field by thunderstorm cells

Data from the high density storm-scale rawinsonde network of the Severe Environmental Storms and Mesoscale Experiment revealed temporal and spatial changes in the divergence fields of the troposphere in response to severe storm evolution on May 2, 1979; these changes were detectable on the meso-beta scale. This unique set of data was subsequently used to study the evolution of the wind, divergence and vertical motion fields in the presence of intense convection. Mid- and upper-tropospheric divergence was superimposed over low-level convergence. The divergence, which has a maximum value of .0004/s, occurred 75 to 100 km upwind as well as over the tornadic cells. To the south of the storm cells, the kinematic pattern was in reverse, upper level convergence was superimposed over low-level divergence. A vertical motion doublet was found to ascend over the squall line and descend about 70 km south of the squall line. It is suggested that the following effects are accountable for the nature of the kinematic fields: (1) blocking of tropospheric environmental flow by the storm cells, (2) anvil outflows, particularly from the tornadic cells, and (3) divergence from the exit region of the jet stream.

Ulanski, S. L.↗

NASA NTRS1986

Navier-Stokes computations of lee-side flows over delta wings

Solutions to the Navier-Stokes equations for the flow over delta wings are computed with emphasis on the separated vortical flows developing on the lee side at high angles of attack. A recently developed implicit algorithm is used which employs upwind differencing for the pressure and convection terms and central differencing for the shear stress and heat transfer terms. Solutions to both the three-dimensional equations and the approximate conical flow equations are compared parametrically with an extensive experimental data base at supersonic speeds. The computations indicate that the conical flow approximation provides results in close agreement with the three-dimensional equations, even to angles of attack as high as 20 degrees. Good agreement with experimentally measured pressures and vapor screen photographs is obtained for the conditions investigated. The method predicts the classical pattern of vortical flow over a delta wing and transition to other flow patterns as the leading edge sweep angle and leading edge normal Mach number are varied.

Thomas, J. L.↗

NASA NTRS1986

Three-dimensional, conservative, Euler computations using patched grid systems and explicit methods

The method of 'zonal approach' (in which the flow field is partitioned into regions with independent grids) for computation of flow over complex geometries, such as aircraft configuration, requires application of a grid-interfacing procedure. A three-dimensional conservative, boundary scheme for patched grids, applicable in generalized coordinates for arbitrary point distributions on a planar zonal surface is presented. The computation technique is derived within the framework of the Osher (1983) upwind scheme, using Euler equations. The three-dimensional interfacing method is applied to the computation of flow about a wing-canard combination using a two-zone, patched grid.

Hessenius, K. A.↗

NASA NTRS1985

A high-precision finite element method for shock-tube calculations

A two-pass explicit scheme is developed in order to exploit some of the capabilities of finite difference modeling (FDM) for finite element modeling (FEM), which offers the opportunity to account for any type of geometry in fluid flow modeling. Features of the first-order upwind and the Lax-Wendroff high precision explicit finite difference algorithms are reviewed. A flux limiter is developed for FEM to serve as an analog for the single limiter function which has been defined for the various FDMs. It is shown that an antidiffusive limiter must be introduced into the weighting function which normally multiplies the time-derivative term in the variational equation. The two-pass scheme which results is demonstrated to be the equivalent of FDMs with five-point support. However, the present scheme is valid only for one-dimensional calculations and linear shape functions for shock tube flow phenomena. Further work is required for its use with nonlinear hyperbolic systems.

Hughes, T. J. R.↗

NASA NTRS1986

Implicit hybrid schemes for the flux-difference split, three-dimensional Navier-Stokes equations

Implicit hybrid algorithms employing symmetric planar Gauss-Seidel (SPGS) relaxation and either block-tridiagonally structured coefficient matrices (AF-SPGS) or block-triangular coefficient matrices (LU-SPGS) are derived to solve the flux-difference-split Navier-Stokes equations for three-dimensional incompressible flow in an upwind scheme. The physical basis of the approach is discussed, and results for problems involving vortex flow around a thin delta wing at Reynolds numbers 900,000 and 10,000 are presented graphically. It is found that AF-SPGS converges faster on vector computers which depend on long vector lengths to achieve optimum performance, whereas LU-SPGS is preferable on sequentially operating machines and vector computers using shorter vector lengths.

Hartwich, P.-M.↗

NASA NTRS1986

On the Nonlinearity of Modern Shock-Capturing Schemes

The development is reviewed of shock capturing methods, paying special attention to the increasing nonlinearity in the design of numerical schemes. The nature is studies of this nonlinearity and its relation to upwind differencing is examined. This nonlinearity of the modern shock capturing methods is essential, in the sense that linear analysis is not justified and may lead to wrong conclusions. Examples to demonstrate this point are given.

Harten, A.↗

NASA NTRS1986

Implicit time-marching solution of the Navier-Stokes equations for thrust reversing and thrust vectoring nozzle flows

An implicit finite volume method is investigated for the solution of the compressible Navier-Stokes equations for flows within thrust reversing and thrust vectoring nozzles. Thrust reversing nozzles typically have sharp corners, and the rapid expansion and large turning angles near these corners are shown to cause unacceptable time step restrictions when conventional approximate factorization methods are used. In this investigation these limitations are overcome by using second-order upwind differencing and line Gauss-Siedel relaxation. This method is implemented with a zonal mesh so that flows through complex nozzle geometries may be efficiently calculated. Results are presented for five nozzle configurations including two with time varying geometries. Three cases are compared with available experimental data and the results are generally acceptable.

Imlay, S. T.↗

NASA NTRS1986

Euler solutions using an implicit multigrid technique

A coarse-grid correction algorithm has been implemented into an implicit upwind Euler solver and tested for transonic airfoil problems. The Euler solver uses split-flux formulation and penta-diagonal scalar equations, respectively, for the explicit and implicit operators. The multigrid sequence starts at the fine grid level, then steps down to each coarse grid level to smooth error components using implicit operators. Estimate of residuals can be obtained by two approaches, which differ in the level where the residuals are collected. Both approaches will lead to a work reduction factor of 12 for a Mach 0.75 flow at 2 degrees incidence on a 65x26 grid. The work reduction factor is found to increase proportional to the number of grid levels.

Li, Chien-Peng↗

NASA NTRS1984

Design of an Advanced Wood Composite Rotor and Development of Wood Composite Blade Technology

In support of a program to advance wood composite wind turbine blade technology, a design was completed for a prototype, 90-foot diameter, two-bladed, one-piece rotor, with all wood/epoxy composite structure. The rotor was sized for compatibility with a generator having a maximum power rating of 4000 kilowatts. Innovative features of the rotor include: a teetering hub to minimize the effects of gust loads, untwisted blades to promote rotor power control through stall, joining of blades to the hub structure via an adhesive bonded structural joint, and a blade structural design which was simplified relative to earlier efforts. The prototype rotor was designed to allow flexibility for configuring the rotor upwind or downwind of the tower, for evaluating various types of teeter dampers and/or elastomeric stops, and with variable delta-three angle settings of the teeter shaft axis. The prototype rotor was also designed with provisions for installing pressure tap and angle of attack instrumentation in one blade. A production version rotor cost analysis was conducted. Included in the program were efforts directed at developing advanced load take-off stud designs for subsequent evaluation testing by NASA, development of aerodynamic tip brake concepts, exploratory testing of a wood/epoxy/graphite concept, and compression testing of wood/epoxy laminate, with scarf-jointed plies.

Stroebel, Thomas↗

NASA NTRS1986

A second-order accurate kinetic-theory-based method for inviscid compressible flows

An upwind method for the numerical solution of the Euler equations is presented. This method, called the kinetic numerical method (KNM), is based on the fact that the Euler equations are moments of the Boltzmann equation of the kinetic theory of gases when the distribution function is Maxwellian. The KNM consists of two phases, the convection phase and the collision phase. The method is unconditionally stable and explicit. It is highly vectorizable and can be easily made total variation diminishing for the distribution function by a suitable choice of the interpolation strategy. The method is applied to a one-dimensional shock-propagation problem and to a two-dimensional shock-reflection problem.

Deshpande, Suresh M.↗