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 847 records · Page 47

Critical effects of downstream boundary conditions on vortex breakdown

The unsteady, compressible, full Navier-Stokes (NS) equations are used to study the critical effects of the downstream boundary conditions on the supersonic vortex breakdown. The present study is applied to two supersonic vortex breakdown cases. In the first case, quasi-axisymmetric supersonic swirling flow is considered in a configured circular duct, and in the second case, quasi-axisymmetric supersonic swirling jet, that is issued from a nozzle into a supersonic jet of lower Mach number, is considered. For the configured duct flow, four different types of downstream boundary conditions are used, and for the swirling jet flow from the nozzle, two types of downstream boundary conditions are used. The solutions are time accurate which are obtained using an implicit, upwind, flux-difference splitting, finite-volume scheme.

Kandil, Osama↗

A comparative study of two codes with an improved two-equation turbulence model for predicting jet plumes

Compressible jet plumes were studied using a two-equation turbulence model. A space marching procedure based on an upwind numerical scheme was used to solve the governing equations and turbulence transport equations. The computed results indicate that extending the space marching procedure for solving supersonic/subsonic mixing problems can be stable, efficient and accurate. Moreover, a newly developed correction for compressible dissipation has been verified in fully expanded and underexpanded jet plumes. For a sonic jet plume, no improvement in results over the standard two-equation model was seen. However for a supersonic jet plume, the correction due to compressible dissipation successfully predicted the reduced spreading rate of the jet compared to the sonic case. The computed results were generally in good agreement with the experimental data.

Balakrishnan, L.↗

A nonlinear relaxation/quasi-Newton algorithm for the compressible Navier-Stokes equations

A highly efficient implicit method for the computation of steady, two-dimensional compressible Navier-Stokes flowfields is presented. The discretization of the governing equations is hybrid in nature, with flux-vector splitting utilized in the streamwise direction and central differences with flux-limited artificial dissipation used for the transverse fluxes. Line Jacobi relaxation is used to provide a suitable initial guess for a new nonlinear iteration strategy based on line Gauss-Seidel sweeps. The applicability of quasi-Newton methods as convergence accelerators for this and other line relaxation algorithms is discussed, and efficient implementations of such techniques are presented. Convergence histories and comparisons with experimental data are presented for supersonic flow over a flat plate and for several high-speed compression corner interactions. Results indicate a marked improvement in computational efficiency over more conventional upwind relaxation strategies, particularly for flowfields containing large pockets of streamwise subsonic flow.

Edwards, Jack R.↗

Spatial and temporal adaptive procedures for the unsteady aerodynamic analysis of airfoils using unstructured meshes

An algorithm which combines spatial and temporal adaption for the time integration of the two-dimensional Euler equations on unstructured meshes of triangles is presented. Spatial adaption involves mesh enrichment to add elements in high gradient regions of the flow and mesh coarsening to remove elements where they are no longer needed. Temporal adaption is a time accurate, local time stepping procedure which integrates the flow equations in each cell according to the local numerical stability constraint. The flow solver utilizes a four-stage Runge-Kutta time integration scheme with an upwind flux-split spatial discretization. Results obtained using spatial and temporal adaption indicate that highly accurate solutions can be obtained with a significant savings of computing time over global time stepping.

Hooker, John R.↗

Numerical simulations using a dynamic solution-adaptive grid algorithm, with applications to unsteady internal flows

An investigation into the numerical simulation of unsteady flows is undertaken using a two-stage Runge-Kutta scheme coupled with the dynamic solution-adaptive grid algorithm developed by the authors. The inviscid fluxes are described by a modified Advective Upwind Split Method to eliminate the need for artificial dissipation. A well-documented numerical example containing moving discontinuities is presented that demonstrates the ability of the coupled grid/solver scheme to accurately capture unsteady flowfield phenomena. Applications are to a typical inlet diffuser configuration at Mach 3.0 with excessive back pressure inducing inlet unstart.

Benson, Rusty A.↗

Viscous flow past a nacelle isolated and in proximity of a flat plate

Linearized-theory design procedures have proven to be useful in preliminary design stages of supersonic aircraft configurations. These procedures are impaired, however, by their inability to account for certain nonlinear effects inherent in complicated flows. The present computations are aimed at providing necessary information for correction and improvement of a particular linearized design method. Three-dimensional, viscous, supersonic flows past nacelle and nacelle-flat plate configurations are investigated. The thin-layer Navier-Stokes equations are solved using an implicit, upwind-biased, finite-volume method. A hybrid domain decomposition technique is utilized to ease the grid generation task. Computations were made for an unit Reynolds number of 2.0 million per foot and Freestream Mach numbers of 1.6, 2.0, and 2.3.

Fouladi, Kamran↗

Solution strategies and heat transfer calculations for three-dimensional configurations at hypersonic speeds

A procedure which reduces the memory requirements for computing the viscous flow over a modified Orbiter geometry at a hypersonic flight condition is presented. The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) code which incorporates a thermochemical nonequilibrium chemistry model, a finite rate catalytic wall boundary condition and wall temperature distribution based on radiation equilibrium is used in this study. In addition, the effect of choice of 'min mod' function, eigenvalue limiter and grid density on surface heating is investigated. The surface heating from a flowfield calculation at Mach number 22, altitude of 230,000 ft and 40 deg angle of attack is compared with flight data from three Orbiter flights.

Weilmuenster, K. J.↗

Characteristics of the Shuttle Orbiter leeside flow during a reentry condition

A study of the leeside flow characteristics of the Shuttle Orbiter is presented for a reentry flight condition. The flow is computed using a point-implicit, finite-volume scheme known as the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA). LAURA is a second-order accurate, laminar Navier-Stokes solver, incorporating finite-rate chemistry with a radiative equilibrium wall temperature distribution and finite-rate wall catalysis. The resulting computational solution is analyzed in terms of salient flow features and the surface quantities are compared with flight data.

Kleb, William L.↗

The Galileo and Pioneer Venus ultraviolet spectrometer experiments - Solar Lyman-alpha latitude variation at solar maximum from interplanetary Lyman-alpha observations

Solar Ly-alpha latitude variation at solar maximum is examined on the basis of interplanetary Ly-alpha observations made during the Galileo and Pioneer Venus UV spectrometer experiments. A comparison is made of the latitude variation of the interplanetary (IP) Ly-alpha signal in 1986 at solar minimum from Pioneer Venus and in 1990 at solar maximum from Galileo. The Galileo EUV spectrometer shows that a large enhancement of the IP Ly-alpha emission occurred over the intervening four years near the solar equator. An IP Ly-alpha model is developed which considers the latitude variation of the solar Ly-alpha flux. The model fit to the data shows a 25-percent decrease of the full disk solar Ly-alpha flux from solar equator to solar pole in 1990. A detailed study of the Galileo IP Ly-alpha observations on day-of-year 190, 193, 197, and 200 in 1990 reveals that large variations occur in response to the 27-d solar variation. Analysis of these data shows that a maximum variation of 20 percent can be expected in the IP Ly-alpha upwind intensity over this 27-d period.

Pryor, W. R.↗

The Osher scheme for non-equilibrium reacting flows

An extension of the Osher upwind scheme to nonequilibrium reacting flows is presented. Owing to the presence of source terms, the Riemann problem is no longer self-similar and therefore its approximate solution becomes tedious. With simplicity in mind, a linearized approach which avoids an iterative solution is used to define the intermediate states and sonic points. The source terms are treated explicitly. Numerical computations are presented to demonstrate the feasibility, efficiency and accuracy of the proposed method. The test problems include a ZND (Zeldovich-Neumann-Doring) detonation problem for which spurious numerical solutions which propagate at mesh speed have been observed on coarse grids. With the present method, a change of limiter causes the solution to change from the physically correct CJ detonation solution to the spurious weak detonation solution.

Suresh, Ambady↗

A comparative study of advanced shock-capturing schemes applied to Burgers' equation

A systematic evaluation is conducted of all extant numerical schemes for nonlinear scalar transport problems, and several advanced shock-capturing schemes are used to solve the nonlinear Burgers' equation in order to characterize their ability to resolve the sharp discontinuity, expansion zone, and propagation and collision features of shocks. For discontinuous functions, the Warming-Beam scheme generates preshock wiggles, while the Lax-Wendroff scheme generates postshock ones. Such limiters as the MUSCL or the superbee are more compressive than minimod or monotonic limiters. The performance of such TVD schemes as the upwind, the symmetric, and the Roe-Sweby, resemble each other.

Yang, H. Q.↗

Navier-Stokes analysis of three-dimensional unsteady flows inside turbine stages

This study presents a numerical method for solving the 3D Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully 3D Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms, while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flowfield of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flowfields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗

Flowfield of a lifting rotor in hover - A Navier-Stokes simulation

The viscous, three-dimensional flowfield of a lifting helicopter rotor in hover is calculated by using an upwind, implicit, finite-difference numerical method for solving the thin layer Navier-Stokes equations. The induced effects of the wake, including the interaction of tip vortices with successive blades, are calculated as part off the overall flowfield solution without using any ad hoc wake models. Comparison of the numerical results for the subsonic and transonic conditions show good agreement with the experimental data and with the previously published Navier-Stokes calculations using a simple wake model. Some comparisons with Euler calculations are also presented, along with some discussions of the grid refinement studies.

Srinivasan, G. R.↗

Active control of asymmetric vortical flows around cones using injection and heating

The effectiveness of certain active-control methods for asymmetric flows around circular cones is investigated by using computational solution of the unsteady, compressible full Navier-Stokes equations. Two main methods of active control which include flow injection and surface heating are used. For the flow-injection-control method, flow injection is used either in the normal direction to the surface or in the tangential direction to the surface. For the surface-heating-control method, the temperature of the cone surface is increased. The effectiveness of a hybrid method of flow control which combines normal injection with surface heating has also been studied. The Navier-Stokes equations, subjected to various surface boundary conditions, are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme for locally-conical flow solutions.

Kandil, Osama A.↗

Computation of vortex wake flows and control of their effects on trailing wings

The near-vortex-wake flow of a large aspect-ratio rectangular wing is accurately computed by using the thin-layer and full Navier-Stokes (NS) equations. The chordwise section of the wing is a NACA-0012 airfoil and its tip is round. The computations have been carried out on a fine C-O grid using an implicit, upwind, flux-difference splitting, finite-volume scheme. The thin-layer NS results have been obtained with and without flux limiters, and the full NS results have been obtained without flux limiters. Flow transition from laminar to turbulent is mimicked by turning-on the Baldwin-Lomax algebraic model at an experimentally prescribed chord-station location of 0.05. Comparison of computed results and experimental data shows that the full NS results give the best resolution of the near-vortex-wake flow. Next, the strength of the wing-tip vortex has been reduced substantially without reducing the lift coefficient by using flow-injection from a slot along a portion of the wing tip. The flow injection is directed in the wing plane at 45 deg with the wing-tip chord.

Wong, Tin-Chee↗

Transonic solutions of a wing/pylon/finned store using hybrid domain decomposition

Transonic Euler calculations about a complex multicomponent configuration are presented. The 3D Euler equations are solved utilizing an upwind-biased, alternating direction implicit, approximately factored, multigrid algorithm. Computational results are compared to experimental data of the finned store in a carriage position.

Newman, James C., III↗

Computing high-speed flows past an oscillating cylinder near a vertical wall

A computational method to simulate unsteady flows involving moving rigid boundaries and interference has been developed. The method is used to solve inviscid equations governing the fluid flow and the dynamic equations governing the motion of rigid bodies. A second-order accurate, upwind-biased, and alterating-direction-implicit method is used to solve the governing equations of the flow. A kinematic domain decomposition (KDD) procedure is extended to treat 3D problems with a high degree of accuracy and generality. The method under consideration is applied to both transonic and supertransonic flows. Both cases involve flow past a cylinder which is forced to pitch sinusoidally near a vertical wall. Benefits of the proposed approach include accurate calculation of the flow around 3D moving multiple bodies with interference; reduction of a numerical error; in particular, the dispersion error which strongly affects wave propagation; and minimization of the phase error which is accumulated according to the time advance procedure.

Yen, Guan-Wei↗

On an origin of numerical diffusion: Violation of invariance under space-time inversion

The invariant properties of the convection equation du/dt + adu/dx = 0 (where d is the partial differential operator) with respect to spatial reflection, time reversal, and space-time inversion are studied. Generally, a finite-difference analog of this equation may possess some or none of these properties. It is shown that, under certain conditions, the von Neumann amplification factor of an analog satisfies a special relation for each invariant property this analog possesses. Particularly, an analog is neutrally stable and thus free of numerical diffusion if it possesses the invariant property related to space-time inversion. It is also explained why generally (1) an upwind scheme possesses neither the invariant property related to spatial reflection nor that related to space-time inversion, and (2) an explicit scheme possesses neither the invariant property related to time reversal nor that related to space-time inversion. Extension to the viscous case and a remarkable connection between the current work and a new numerical framework for solving conservation laws are also discussed.

Chang, Sin-Chung↗