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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Calculation of AGARD Wing 445.6 Flutter Using Navier-Stokes Aerodynamics

An unsteady, 3D, implicit upwind Euler/Navier-Stokes algorithm is here used to compute the flutter characteristics of Wing 445.6, the AGARD standard aeroelastic configuration for dynamic response, with a view to the discrepancy between Euler characteristics and experimental data. Attention is given to effects of fluid viscosity, structural damping, and number of structural model nodes. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. The V-g analysis indicates that fluid viscosity has a significant effect on the supersonic flutter boundary for this wing.

Lee-Rausch, Elizabeth M.↗
NASA NTRS1993

Active control of asymmetric conical flow using spinning and rotatory oscillations

The effectiveness of active control on asymmetric flows around circular cones is investigated computationally using cone spinning and rotatory oscillation around its axis. The investigation uses the time-accurate solution of the unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, flux-difference splitting, finite-volume scheme. The present solutions are obtained under the locally-conical-flow assumption in order to understand the flow physics using very fine grids for reasonable flow resolution at low computational cost. For all the computational solutions, a grid of 241 x 81 x 2 points in the wrap-around, normal and axial directions, respectively, is used. The grid is spinning or oscillating rigidly with the cone according to its motion and the kinematical and dynamical boundary conditions are modified accordingly. The computational applications include the effects of uniform spinning rates and periodic rotatory oscillations at different amplitudes and frequencies on the flow asymmetry.

Kandil, Osama A.↗
NASA NTRS1993

Shock-vortex interaction over a 65-degree delta wing in transonic flow

Transonic flow over a 65-deg swept-back, sharp-edged, cropped delta wing is investigated computationally using the time-accurate solution of the unsteady, compressible, full Navier-Stokes equations with an implicit, upwind, flux-difference splitting, finite-volume scheme. Coarse and fine O-H grids are used to obtain the solution. The grid consists of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively. The results are presented for an angle of attack of 20 deg Mach number of 0.85 and Reynolds number of 3.23 x 10 exp 6. With the fine grid, the results show that a system of shocks has been captured over the upper wing surface and that the leading-edge vortex core experiences an unsteady supersonic vortex breakdown after passing through a spanwise shock near the wing trailing edge. The computed results at a certain time are in good agreement with the experimental data. Topological aspects of the vortex breakdown flowfield are also presented and discussed.

Kandil, Osama A.↗
NASA NTRS1993

Simulation of tail buffet using delta wing-vertical tail configuration

Computational simulation of the vertical tail buffet problem is accomplished using a delta wing-vertical tail configuration. Flow conditions are selected such that the wing primary-vortex cores experience vortex breakdown and the resulting flow interacts with the vertical tail. This multidisciplinary problem is solved successively using three sets of equations for the fluid flow, aeroelastic deflections and grid displacements. For the fluid dynamics part, the unsteady, compressible, full Navier-Stokes equations are solved accurately in time using an implicit, upwind, flux-difference splitting, finite-volume scheme. For the aeroelastic part, the aeroelastic equation for bending vibrations is solved accurately in time using the Galerkin method and the four-stage Runge-Kutta scheme. The grid for the fluid dynamics computations is updated every few time steps using a third set of interpolation equations. The computational application includes a delta wing of aspect ratio 1 and a rectangular vertical tail of aspect ratio 2, which is placed at 0.5 root-chord length downstream of the wing trailing edge. The wing angle of attack is 35 deg and the flow Mach number and Reynolds number are 0.4 and 10,000, respectively.

Kandil, Osama A.↗
NASA NTRS1992

Point-implicit relaxation strategies for viscous, hypersonic flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's symmetric total variation diminishing scheme. Newton relaxation of the fully coupled equation set is employed on a cell-to-cell basis. Under-relaxation of the inviscid and over-relaxation of the viscous contributions to the residual are implemented. Computational work is easily partitioned among many processors in an asynchronous, dynamic mode for convergence acceleration. An overview of the physical models employed herein for thermochemical nonequilibrium is included. Several test cases and comparisons with experimental data are presented involving hypersonic flow over blunt bodies which illustrate the qualitative and quantitative capabilities of this approach.

Gnoffo, P. A.↗
NASA NTRS1993

CFD validation of subsonic turbulent planar shear layers

The primary objective of the present study is to assess the limitations and capabilities of RPLUS in predicting the entrainment, mixing, and burning characteristics of a high subsonic planar shear layer; these computations are performed in conjunction with an experiment being conducted at NASA Lewis Research Center. Turbulence is modeled by a two-equation k-epsilon closure modified for compressibility effects. Combustion is modeled by two different models: one in which a finite-rate laminar chemistry model is used for H2-O2 oxidation, and the other based on a composition joint pdf approach to account for the turbulence-chemistry interaction effects. Between upwind and centered differencing schemes, minimal differences are observed for the flow properties of a nonreacting shear layer case studied. In comparison with experimental data, computed growth rates are underpredicted in both reacting and nonreacting cases of a shear layer, but more severely for the reacting flows. However, the qualitative agreement between computation and experiment is reasonably good for the other observable characteristics, and indicative of potential for further improvements.

Lai, H. T.↗
NASA NTRS1993

Unsteady aerodynamic flow phenomena in a transonic compressor stage

A three-dimensional unsteady, viscous aerodynamic analysis has been developed for the flow inside a transonic, high-through-flow, single stage compressor. The compressor stage is comprised of a low-aspect-ratio rotor and a closely coupled stator. The analysis is based on a numerical method for solving the three-dimensional Navier-Stokes equation for unsteady viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equation with an implicit scheme. A two-equation turbulence model with a low-Reynolds-number modification is applied for the turbulence closure. 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 the viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The numerical results are compared with available experimental data.

Hah, C.↗
NASA NTRS1993

Computations of spray, fuel-air mixing, and combustion in a lean-premixed-prevaporized combustor

A code was developed for computing the multidimensional flow, spray, combustion, and pollutant formation inside gas turbine combustors. The code developed is based on a Lagrangian-Eulerian formulation and utilizes an implicit finite-volume method. The focus of this paper is on the spray part of the code (both formulation and algorithm), and a number of issues related to the computation of sprays and fuel-air mixing in a lean-premixed-prevaporized combustor. The issues addressed include: (1) how grid spacings affect the diffusion of evaporated fuel, and (2) how spurious modes can arise through modelling of the spray in the Lagrangian computations. An upwind interpolation scheme is proposed to account for some effects of grid spacing on the artificial diffusion of the evaporated fuel. Also, some guidelines are presented to minimize errors associated with the spurious modes.

Dasgupta, A.↗
NASA NTRS1993

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results are shown a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Steinthorsson, E.↗
NASA NTRS1993

Central difference TVD schemes for time dependent and steady state problems

A Runge-Kutta formula in time is presently used to advance schemes in which central differences are used to solve the time-dependent Euler equations; a second difference is added near shocks as an artificial viscosity to reduce the given scheme to a first-order upwind one at shocks. A matrix-valued dissipation is introduced and compared with the scalar viscosity; a connection is shown between this artificial viscosity and flux limiters. The use of various flux limiters for this central difference scheme is compared.

Jorgenson, P.↗
NASA NTRS1993

Three-dimensional flow analysis inside turbomachinery stages with steady and unsteady Navier-Stokes method

This study presents a numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional 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 in the present study. 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 flow field of a subsonic turbine stage and the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid.

Copenhaver, W. W.↗
NASA NTRS1993

Multiple solutions in a rotating annulus flow model

A series of numerical experiments is conducted for rotating annulus flow using Miller et al.'s (1992) Geophysical Flow Simulation (GFS) model; a mixture of 25-percent upwind-differencing and 75-percent center-differencing is employed to approximate the temperature advective terms. Attention is given to the wavenumber selection time and wavenumber regimes, the sensitivity in the wavenumber transition regions, and hysteresis and irregular wavenumber selections.

Lu, Huei-Iin↗
NASA NTRS1993

A laboratory study of friction-velocity estimates from scatterometry - Low and high regimes

Measurements from scatterometers pointing at wind-waves in three large wave tanks are examined to study fetch effects and the correlation with wind friction velocity. Time-series measurements were made at 13, 35, and 95 m with a Ka-band scatterometer aimed upwind at 30 deg incidence angle and vertical polarization. Average normalized radar cross-section (sigma-0) values from all fetches follow a common trend for sigma-0 as a function of wind friction velocity, so the fetch dependence is negligible. An empirical power-law model yields a high correlation between sigma-0 and wind friction velocity, but, because systematic anomalies arise, we reexamine a turbulence approach that delineates low and high regimes with a transition at a wind friction velocity of approximately 25 cm/s. Using this criteria, the data are well represented by a two-section power-law relationship between sigma-0 and wind friction velocity.

Bliven, L. F.↗
NASA NTRS1993

A CFD-based aerodynamic design procedure for hypersonic wind-tunnel nozzles

A new procedure which unifies the best of current classical design practices, computational fluid dynamics (CFD), and optimization procedures is demonstrated for designing the aerodynamic lines of hypersonic wind-tunnel nozzles. The new procedure can be used to design hypersonic wind tunnel nozzles with thick boundary layers where the classical design procedure has been shown to break down. An efficient CFD code, which solves the parabolized Navier-Stokes (PNS) equations using an explicit upwind algorithm, is coupled to a least-squares (LS) optimization procedure. A LS problem is formulated to minimize the difference between the computed flow field and the objective function, consisting of the centerline Mach number distribution and the exit Mach number and flow angle profiles. The aerodynamic lines of the nozzle are defined using a cubic spline, the slopes of which are optimized with the design procedure. The advantages of the new procedure are that it allows full use of powerful CFD codes in the design process, solves an optimization problem to determine the new contour, can be used to design new nozzles or improve sections of existing nozzles, and automatically compensates the nozzle contour for viscous effects as part of the unified design procedure. The new procedure is demonstrated by designing two Mach 15, a Mach 12, and a Mach 18 helium nozzles. The flexibility of the procedure is demonstrated by designing the two Mach 15 nozzles using different constraints, the first nozzle for a fixed length and exit diameter and the second nozzle for a fixed length and throat diameter. The computed flow field for the Mach 15 least squares parabolized Navier-Stokes (LS/PNS) designed nozzle is compared with the classically designed nozzle and demonstrates a significant improvement in the flow expansion process and uniform core region.

Korte, John J.↗
NASA NTRS1993

A preliminary C-band scatterometer model function for the ERS-1 AMI instrument

Backscatter cross section (sigma(sub 0)) measurements from the ERS-1 scatterometer are collocated with surface wind velocity data from operational Numerical Weather Prediction (NWP) analyses to yield a fully empirical C band model function. The empirical model does not obey a power law at speeds less than 5 m/s, and upwind/crosswind ratios are small for all incidence angles at low wind speeds and for small incidence angles at all wind speeds. Downwind sigma(sub 0) exceeds upwind values for incidence angles below approximately 28 degrees. The full ERS-1 scatterometer data set has been reprocessed using the NWP model function.

Freilich, M. H.↗
NASA NTRS1993

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping, and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results shown are a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Steinthorsson, E.↗
NASA NTRS1993

Efficient simulation of incompressible viscous flow over multi-element airfoils

The incompressible, viscous, turbulent flow over single and multi-element airfoils is numerically simulated in an efficient manner by solving the incompressible Navier-Stokes equations. The solution algorithm employs the method of pseudo compressibility and utilizes an upwind differencing scheme for the convective fluxes, and an implicit line-relaxation scheme. The motivation for this work includes interest in studying high-lift take-off and landing configurations of various aircraft. In particular, accurate computation of lift and drag at various angles of attack up to stall is desired. Two different turbulence models are tested in computing the flow over an NACA 4412 airfoil; an accurate prediction of stall is obtained. The approach used for multi-element airfoils involves the use of multiple zones of structured grids fitted to each element. Two different approaches are compared; a patched system of grids, and an overlaid Chimera system of grids. Computational results are presented for two-element, three-element, and four-element airfoil configurations. Excellent agreement with experimental surface pressure coefficients is seen. The code converges in less than 200 iterations, requiring on the order of one minute of CPU time on a CRAY YMP per element in the airfoil configuration.

Rogers, Stuart E.↗
NASA NTRS1993

A comparison of two incompressible Navier-Stokes algorithms for unsteady internal flow

A comparative study of two different incompressible Navier-Stokes algorithms for solving an unsteady, incompressible, internal flow problem is performed. The first algorithm uses an artificial compressibility method coupled with upwind differencing and a line relaxation scheme. The second algorithm uses a fractional step method with a staggered grid, finite volume approach. Unsteady, viscous, incompressible, internal flow through a channel with a constriction is computed using the first algorithm. A grid resolution study and parameter studies on the artificial compressibility coefficient and the maximum allowable residual of the continuity equation are performed. The periodicity of the solution is examined and several periodic data sets are generated using the first algorithm. These computational results are compared with previously published results computed using the second algorithm and experimental data.

Wiltberger, N. Lyn↗