Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “upwind”

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.

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

High angle-of-attack inviscid Shuttle Orbiter computation

Results are presented of a comparison between the pressure distributions predicted by the perfect-gas computational fluid dynamics and the Shuttle Orbiter wind-tunnel data for high angles of attack, using the LAURA (for Langley Aerothermodynamic Upwind Relaxation Algorithm) as applied to the wind-tunnel condition to predict the flow over the vehicle. It is shown that the calculated pressures compare well with the wind tunnel data for both the windward and the leeward sides, indicating that the salient inviscid flow features were properly modeled.

Kleb, William L.↗
NASA NTRS1992

Navier-Stokes computations of cavity aeroacoustics with suppression devices

Effectiveness of two devices to suppress the cavity acoustics was computationally investigated. Two dimensional, computational simulations were performed for the transonic, turbulent flows past a cavity, which was first equipped with a rear face ramp and then with a spoiler. The Reynolds-averaged, unsteady, compressible, full Navier-Stokes equations were solved time accurately by a second-order accurate, implicit, upwind, finite-volume method. The effect of turbulence was included through a Baldwin-Lomax model with modifications for the multiple-wall effects and for the highly vortical flow with a shear layer. The results included instantaneous and time-averaged flow properties, and time-series analyses of the pressure inside the cavity, which compared favorably with the available experimental data. These results were also contrasted with the computed aeroacoustics of the same cavity (length-to-depth ratio of 4.5), but without a device, to demonstrate the suppression effectiveness.

Baysal, Oktay↗
NASA NTRS1993

Spatial adaptation procedures on tetrahedral meshes for unsteady aerodynamic flow calculations

Spatial adaptation procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaptation procedures were developed and implemented within a three-dimensional, unstructured-grid, upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. The paper gives a detailed description of the enrichment and coarsening procedures and presents comparisons with experimental data for an ONERA M6 wing and an exact solution for a shock-tube problem to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady results, obtained using spatial adaptation procedures, are shown to be of high spatial accuracy, primarily in that discontinuities such as shock waves are captured very sharply.

Rausch, Russ D.↗
NASA NTRS1992

CFD applications in hypersonic flight

Design studies are underway for a variety of hypersonic flight vehicles. The National Aero-Space Plane will provide a reusable, single-stage-to-orbit capability for routine access to low earth orbit. Flight-capable satellites will dip into the atmosphere to maneuver to new orbits, while planetary probes will decelerate at their destination by atmospheric aerobraking. To supplement limited experimental capabilities in the hypersonic regime, CFD is being used to analyze the flow about these configurations. The governing equations include fluid dynamic as well as chemical species equations, which are solved with robust upwind differencing schemes. Examples of CFD applications to hypersonic vehicles suggest an important role this technology will play in the development of future aerospace systems. The computational resources needed to obtain solutions are large, but various strategies are being exploited to reduce the time required for complete vehicle simulations.

Edwards, T. A.↗
NASA NTRS1993

Heating analysis for a Lunar Transfer Vehicle at near-equilibrium flow conditions

A heating analysis for a 15.2 m diameter Lunar Transfer Vehicle (LTV) at 0 and 10.6 deg angle of attack for a nominal trajectory through the earth's atmosphere is described. The analysis utilizes the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) with thin-layer, Navier-Stokes, thermochemical nonequilibrium options. Radiative heating levels are calculated using the Langley Optimized RAdiative Nonequilibrium (LORAN) and the Non-EQuilibrium AIr Radiation (NEQAIR) codes. At peak heating, the shock layer is substantially in equilibrium. Comprehensive spatial and spectral grid convergence studies have been implemented to quantify grid effects on the convective and radiative heating levels. Axisymmetric tests including the coupled effects of radiative energy transfer show negligible change to the convective heating and a 20 percent reduction in the radiative heating.

Gnoffo, Peter A.↗
NASA NTRS1993

Fourier analysis of numerical algorithms for the Maxwell equations

The Fourier method is used to analyze the dispersive, dissipative, and isotropy errors of various spatial and time discretizations applied to the Maxwell equations on multi-dimensional grids. Both Cartesian grids and non-Cartesian grids based on hexagons and tetradecahedra are studied and compared. The numerical errors are quantitatively determined in terms of phase speed, wave number, propagation direction, gridspacings, and CFL number. The study shows that centered schemes are more efficient than upwind schemes. The non-Cartesian grids yield superior isotropy and higher accuracy than the Cartesian ones. For the centered schemes, the staggered grids produce less errors than the unstaggered ones. A new unstaggered scheme which has all the best properties is introduced. The study also demonstrates that a proper choice of time discretization can reduce the overall numerical errors due to the spatial discretization.

Liu, Yen↗
NASA NTRS1993

Juncture flow improvement for wing/pylon configurations by using CFD methodology

Transonic flow field around a fighter wing/pylon configuration was simulated by using an implicit upwinding Navier-Stokes flow solver (F3D) and overset grid technology (Chimera). Flow separation and local shocks near the wing/pylon junction were observed in flight and predicted by numerical calculations. A new pylon/fairing shape was proposed to improve the flow quality. Based on numerical results, the size of separation area is significantly reduced and the onset of separation is delayed farther downstream. A smoother pressure gradient is also obtained near the junction area. This paper demonstrates that computational fluid dynamics (CFD) methodology can be used as a practical tool for aircraft design.

Gea, Lie-Mine↗
NASA NTRS1993

Three-dimensional supersonic vortex breakdown

Three-dimensional supersonic vortex-breakdown problems in bound and unbound domains are solved. The solutions are obtained using the time-accurate integration of the unsteady, compressible, full Navier-Stokes (NS) equations. The computational scheme is an implicit, upwind, flux-difference splitting, finite-volume scheme. Two vortex-breakdown applications are considered in the present paper. The first is for a supersonic swirling jet which is issued from a nozzle into a supersonic uniform flow at a lower Mach number than that of the swirling jet. The second is for a supersonic swirling flow in a configured circular duct. In the first application, an extensive study of the effects of grid fineness, shape and grid-point distribution on the vortex breakdown is presented. Four grids are used in this study and they show a substantial dependence of the breakdown bubble and shock wave on the grid used. In the second application, the bubble-type and helix-type vortex breakdown have been captured.

Kandil, Osama A.↗
NASA NTRS1993

Three-dimensional unstructured grid Euler method applied to turbine blades

Flow through a turbine annular cascade is calculated using a three-dimensional Euler method based on unstructured tetrahedral meshes. The equations are integrated in time using an explicit Runge-Kutta time-stepping scheme. The inviscid flux terms are discretized using a cell-centered finite-volume formulation with upwind flux-difference splitting. The tetrahedral meshes around the turbine blade are generated using an advancing-front technique with forced geometric periodicity between the blades. Good agreement is obtained between the present calculation and the experiment for both surface pressure distribution and flow behavior in the passage between the blades, demonstrating the capability of the present methodology for turbomachinery flow applications.

Kwon, Oh J.↗
NASA NTRS1993

Effect of leading-edge porosity on blade-vortex interaction noise

The effect of the porous leading-edge of an airfoil on the blade-vortex interaction noise, which dominates far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near-field are calculated directly from the computation. The porosity effect on the surface is modeled. Results show leading-edge transpiration can suppress pressure fluctuations at the leading-edge during BVI, and consequently reduce the amplitude of propagating noise by 30 percent at maximum in the near-field. The effect of porosity factor on the noise level is also investigated.

Lee, Soogab↗
NASA NTRS1993

On the accuracy of limiters and convergence to steady state solutions

This paper addresses the practical problem of obtaining convergence to steady state solutions when limiters are used in conjunction with upwind schemes on unstructured grids. The base scheme forms a gradient and limits it by imposing monotonicity conditions in the reconstruction stage. It is shown by analysis in one dimension that such an approach leads to various schemes meeting TVD requirements in one dimension. It is further shown that these formally second order accurate schemes are less than second order accurate in practice because of the action of the limiter function in smooth regions of the solution. Modifications are proposed to the limiter that restore the second order accuracy. In multiple dimensions these schemes produce steady state solutions that are monotone and devoid of oscillations. However, convergence stalls after a few orders of reduction in the residual. With the modified limiter, on the other hand, it is shown that converged steady state solutions can be obtained.

Venkatakrishnan, V.↗
NASA NTRS1993

Preconditioned Conjugate Gradient methods for low speed flow calculations

An investigation is conducted into the viability of using a generalized Conjugate Gradient-like method as an iterative solver to obtain steady-state solutions of very low-speed fluid flow problems. Low-speed flow at Mach 0.1 over a backward-facing step is chosen as a representative test problem. The unsteady form of the two dimensional, compressible Navier-Stokes equations are integrated in time using discrete time-steps. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. The new iterative solver is used to solve a linear system of equations at each step of the time-integration. Preconditioning techniques are used with the new solver to enhance the stability and the convergence rate of the solver and are found to be critical to the overall success of the solver. A study of various preconditioners reveals that a preconditioner based on the lower-upper (L-U)-successive symmetric over-relaxation iterative scheme is more efficient than a preconditioner based on incomplete L-U factorizations of the iteration matrix. The performance of the new preconditioned solver is compared with a conventional line Gauss-Seidel relaxation (LGSR) solver. Overall speed-up factors of 28 (in terms of global time-steps required to converge to a steady-state solution) and 20 (in terms of total CPU time on one processor of a CRAY-YMP) are found in favor of the new preconditioned solver, when compared with the LGSR solver.

Ajmani, Kumud↗
NASA NTRS1993

Recent developments in high order K-exact reconstruction on unstructured meshes

This paper presents recent improvements in high-order K-exact reconstruction on unstructured meshes. The new reconstruction procedures are incorporated into a basic upwind finite-volume scheme suitable for solving scalar advection-diffusion equations as well as the Euler and Navier-Stokes equations. Numerical calculations are performed comparing the present method with lower order accurate reconstruction procedures (piecewise constant and piecewise linear) and various competing technologies such as the fluctuation splitting method of Roe (1987) and Deconinck et al. (1992) and a system-variant of the streamline diffusion Petrov-Galerkin method developed by Hansbo (1991) and Hansbo and Johnson (1991). Five test problems are used in the numerical comparisons: scalar circular advection, transonic and supersonic Euler flow, laminar boundary-layer flow, and general compressible Navier-Stokes flow.

Barth, Timothy J.↗
NASA NTRS1993

Incompressible Navier-Stokes computations of rotating flows

Flow through pump components, such as an inducer and an impeller, is efficiently simulated by solving the incompressible Navier-Stokes equations. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. The resulting computer code is applied to the flow analysis inside a generic rocket engine pump inducer, a fuel pump impeller, and SSME high-pressure fuel turbopump impeller. Numerical results of inducer flow are compared with experimental measurements. Flow analyses at 80-, 100-, and 120-percent of design conditions are presented.

Kiris, Cetin↗
NASA NTRS1992

Prediction of asymmetric vortical flows around slender bodies using Navier-Stokes equations

Steady and unsteady asymmetric vortical flows around slender bodies at high angles of attack are solved using the unsteady, compressible, thin-layer Navier-Stokes equations. An implicit, upwind-biased, flux-difference splitting, finite-volume scheme is used for the numerical computations. For supersonic flows past point cones, the locally conical flow assumption have been used for efficient computational studies of this phenomenon. Asymmetric flows past a 5-deg semiapex-angle circular cone at different angles of attack, free-stream Mach numbers, and Reynolds numbers have been studied in responses to different sources of disturbances. The effects of grid fineness and computational domain size have also been investigated. Next, the responses of three-dimensional supersonic asymmetric flow around a 5-deg circular cone at different angles of attack and Reynolds numbers to short-duration sideslip disturbances are presented. The results show that flow asymmetry becomes stronger as the Reynolds number and angles of attack are increased. One of the cases of flow over a cone-cylinder configuration is validated fairly well by experimental data.

Liu, C. H.↗
NASA NTRS1991

Euler solutions to nonlinear acoustics of non-lifting rotor blades

For the first time a computational fluid dynamics (CFD) method is used to calculate directly the high-speed impulsive (HSI) noise of a non-lifting hovering rotor blade out to a distance of over three rotor radii. In order to accurately propagate the acoustic wave in a stable and efficient manner, an implicit upwind-biased Euler method is solved on a grid with points clustered along the line of propagation. A detailed validation of the code is performed for a rectangular rotor blade at tip Mach numbers ranging from 0.88 to 0.92. The agreement with experiment is excellent at both the sonic cylinder and at 2.18 rotor radii. The agreement at 3.09 rotor radii is still very good, showing improvements over the results from the best previous method. Grid sensitivity studies indicate that with special attention to the location of the boundaries a grid with approximately 60,000 points is adequate. This results in a computational time of approximately 40 minutes on a Cray-XMP. The practicality of the method to calculate HSI noise is demonstrated by expanding the scope of the investigation to examine the rectangular blade as well as a highly swept and tapered blade over a tip Mach number range of 0.80 to 0.95. Comparisons with experimental data are excellent and the advantages of planform modifications are clearly evident. New insight is gained into the mechanisms of nonlinear propagation and the minimum distance at which a valid comparison of different rotors can be made: approximately two rotor radii from the center of rotation.

Baeder, J. D.↗
NASA NTRS1991

Finite element computation of compressible flows with the SUPG formulation

Finite element computation of compressible Euler equations is presented in the context of the streamline-upwind/Petrov-Galerkin (SUPG) formulation. The SUPG formulation, which is based on adding stabilizing terms to the Galerkin formulation, is further supplemented with a shock capturing operator which addresses the difficulty in maintaining a satisfactory solution near discontinuities in the solution field. The shock capturing operator, which has been derived from work done in entropy variables for a similar operator, is shown to lead to an appropriate level of additional stabilization near shocks, without resulting in excessive numerical diffusion. An implicit treatment of the impermeable wall boundary condition is also presented. This treatment of the no-penetration condition offers increased stability for large Courant numbers, and accelerated convergence of the computations for both implicit and explicit applications. Several examples are presented to demonstrate the ability of this method to solve the equations governing compressible fluid flow.

Le Beau, G. J.↗
NASA NTRS1993

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗