Search NASA⌕ Search

SEARCH · Search NASA

Results for “upwind schemes”

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 487 records · Page 27

Natural Rolling Responses of a Delta Wing in Transonic and Subsonic Flows

The unsteady, three-dimensional, full Navier-Stokes (NS) equations and the Euler equations of rigid-body dynamics are sequentially solved to simulate the natural rolling response of slender delta wings of zero thickness at moderate to high angles of attack, to transonic and subsonic flows. The governing equations of fluid flow and dynamics of the present multi-disciplinary problem are solved using the time-accurate solution of the NS equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The main focus is to analyze the effect of Mach number and angle of attack on the leading edge vortices and their breakdown, the resultant rolling motion, and overall aerodynamic response of the wing. Three cases demonstrate the natural response of a 65 deg swept, cropped delta wing in a transonic flow with breakdown of the leading edge vortices and an 80 deg swept delta wing in a subsonic flow undergoing either damped or self-excited limit-cycle rolling oscillations as a function of angle of attack. Comparisons with an experimental investigation completes this study, validating the analysis and illustrating the complex details afforded by computational investigations.

Menzies, Margaret A.↗

Unsteady, Transonic Flow Around Delta Wings Undergoing Coupled and Natural Modes Response: A Multidisciplinary Problem

The unsteady, three-dimensional Navier-Stokes equations coupled with the Euler equations of rigid-body dynamics are sequentially solved to simulate and analyze the aerodynamic response of a high angle of attack delta wing undergoing oscillatory motion. The governing equations of fluid flow and dynamics of the multidisciplinary problem are solved using a time-accurate solution of the laminar, unsteady, compressible, full Navier- Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The primary model under consideration consists of a 65 deg swept, sharp-edged, cropped delta wing of zero thickness at 20 deg angle of attack. In a freestream of Mach 0.85 and Reynolds number of 3.23 x 10(exp 6), the flow over the upper surface of the wing develops a complex shock system which interacts with the leading-edge primary vortices producing vortex breakdown. The effect of the oscillatory motion of the wing on the vortex breakdown and overall aerodynamic response is detailed to provide insight to the complicated physics associated with unsteady flows and the phenomenon of wing rock. Forced sinusoidal single and coupled mode rolling and pitching motion is presented for the wing in a transonic freestream. The Reynolds number, frequency of oscillation, and the phase angle are varied. Comparison between the single and coupled mode forced rolling and pitching oscillation cases illustrate the effects of coupling the motion. This investigation shows that even when coupled, forced rolling oscillation at a reduced frequency of 2(pi) eliminates the vortex breakdown which results in an increase in lift. The coupling effect for in phase forced oscillations show that the lift coefficient of the pitching-alone case and the rolling-moment coefficient of the rolling-alone case dominate the resulting response. However, with a phase lead in the pitching motion, the coupled motion results in a non-periodic response of the rolling moment. The second class of problems involve releasing the wing in roll to respond to the flowfield. Two models of sharp-edged delta wings, the previous 65 deg swept model and an 80 deg swept, sharp-edged delta wing, are used to observe the aerodynamic response of a wing free to roll in a transonic and subsonic freestream, respectively. These cases demonstrate damped oscillations, self-sustained limit cycle oscillations, and divergent rolling oscillations. Ultimately, an active control model using a mass injection system was applied on the surface of the wing to suppress the self-sustained limit cycle oscillation known as wing rock. Comparisons with experimental investigations complete this study, validating the analysis and illustrating the complex details afforded by computational investigations.

Menzies, Margaret Anne↗

ULTRA-SHARP solution of the Smith-Hutton problem

Highly convective scalar transport involving near-discontinuities and strong streamline curvature was addressed in a paper by Smith and Hutton in 1982, comparing several different convection schemes applied to a specially devised test problem. First order methods showed significant artificial diffusion, whereas higher order methods gave less smearing but had a tendency to overshoot and oscillate. Perhaps because unphysical oscillations are more obvious than unphysical smearing, the intervening period has seen a rise in popularity of low order artificially diffusive schemes, especially in the numerical heat transfer industry. The present paper describes an alternate strategy of using non-artificially diffusive high order methods, while maintaining strictly monotonic transitions through the use of simple flux limited constraints. Limited third order upwinding is usually found to be the most cost effective basic convection scheme. Tighter resolution of discontinuities can be obtained at little additional cost by using automatic adaptive stencil expansion to higher order in local regions, as needed.

Leonard, B. P.↗

A Navier-Stokes solver for high speed equilibrium flows and application to blunt bodies

This paper presents a finite element method for the solution of Navier-Stokes equations with the assumption of thermodynamic and chemical equilibrium. The method employs an upwind finite element technique with an implicit time-marching scheme for the solution, and uses an adaptively generated unstructured triangular mesh with several layers of quadrilateral elements near solid walls. The complexity associated with the assumption that the flow is in equilibrium is treated consistently, and the inviscid flux Jacobian matrices are derived. Several problems involving inviscid and viscous hypersonic flow past blunt are solved. Results are compared with other numerical results and experimental data, and generally good agreement is observed.

Prabhu, Ramadas K.↗

Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows

A fundamentally multidimensional convection scheme is described based on vector transient interpolation modeling rewritten in conservative control-volume form. Vector third-order upwinding is used as the basis of the algorithm; this automatically introduces important cross-difference terms that are absent from schemes using component-wise one-dimensional formulas. Third-order phase accuracy is good; this is important for coarse-grid large-eddy or full simulation. Potential overshoots or undershoots are avoided by using a recently developed universal limiter. Higher order accuracy is obtained locally, where needed, by the cost-effective strategy of adaptive stencil expansion in a direction normal to each control-volume face; this is controlled by monitoring the absolute normal gradient and curvature across the face. Higher (than third) order cross-terms do not appear to be needed. Since the wider stencil is used only in isolated narrow regions (near discontinuities), extremely high (in this case, seventh) order accuracy can be achieved for little more than the cost of a globally third-order scheme.

Leonard, B. P.↗

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.↗

A three-dimensional upwind parabolized Navier-Stokes code for real gas flows

A real gas, upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the three-dimensional hypersonic flow of equilibrium air around various body shapes. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence et al. (1986). The upwind algorithm is based on Roe's (1981) flux-difference splitting scheme which has been modified to account for real gas effects using the nearly exact approach of Vinokur and Liu (1988). Simplified curve fits are employed to obtain the thermodynamic and transport properties of equilibrium air. The new code has been validated by computing the M-infinity = 25 laminar flow of air over cones at various angles of attack. The results of these computations are compared with the results from a conventional centrally-differenced, real gas PNS code and the previous axisymmetric, upwind, real gas code. The agreement is excellent in all cases.

Tannehill, John C.↗

ULTRA-SHARP nonoscillatory convection schemes for high-speed steady multidimensional flow

For convection-dominated flows, classical second-order methods are notoriously oscillatory and often unstable. For this reason, many computational fluid dynamicists have adopted various forms of (inherently stable) first-order upwinding over the past few decades. Although it is now well known that first-order convection schemes suffer from serious inaccuracies attributable to artificial viscosity or numerical diffusion under high convection conditions, these methods continue to enjoy widespread popularity for numerical heat transfer calculations, apparently due to a perceived lack of viable high accuracy alternatives. But alternatives are available. For example, nonoscillatory methods used in gasdynamics, including currently popular TVD schemes, can be easily adapted to multidimensional incompressible flow and convective transport. This, in itself, would be a major advance for numerical convective heat transfer, for example. But, as is shown, second-order TVD schemes form only a small, overly restrictive, subclass of a much more universal, and extremely simple, nonoscillatory flux-limiting strategy which can be applied to convection schemes of arbitrarily high order accuracy, while requiring only a simple tridiagonal ADI line-solver, as used in the majority of general purpose iterative codes for incompressible flow and numerical heat transfer. The new universal limiter and associated solution procedures form the so-called ULTRA-SHARP alternative for high resolution nonoscillatory multidimensional steady state high speed convective modelling.

Leonard, B. P.↗

An upwind parabolized Navier-Stokes code for real gas flows

A real gas, upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the two-dimensional/axisymmetric hypersonic flow of equilibrium air around various body shapes. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence, Tannehill and Chaussee. The upwind algorithm is based on Roe's flux-difference splitting scheme which has been modified to account for real gas effects. Simplified curve fits are used to obtain the thermodynamic and transport properties of equilibrium air. The new code has been validated by computing the hypersonic laminar flow of air over a flat plate, a wedge, a ramp, and a cone. The results of these computations are compared with the results from a conventional centrally-differenced, real gas, PNS code and the agreement is excellent, except in the vicinity of shock waves where the present code exhibits far superior shock capturing capabilities.

Tannehill, John C.↗

An upwind parabolized Navier-Stokes code for chemically reacting flows

A new upwind, parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic, viscous, chemically reacting flow around two-dimensional or axisymmetric bodies. The new code is an extension of the upwind (perfect gas) PNS code of Lawrence et al. (1986). The upwind algorithm is based on Roe's flux-difference splitting scheme which has been modified to account for real gas effects. The algorithm solves the gas dynamic and species continuity equations in a 'loosely' coupled manner. The new code has been validated by computing the laminar flow (at free stream Mach number 25) of chemically reacting air over a wedge and a cone. The results of these computations are compared with the results from a centrally-differenced, fully coupled, nonequilibrium PNS code. The agreement is excellent, except in the vicinity of the shock wave where the present code exhibits superior shock capturing capabilities.

Tannehill, John C.↗

Supersonic Navier-Stokes simulations of turbulent afterbody flows

Four numerical algorithms and three turbulence models are used to solve the thin layer three-dimensional Reynolds-averaged Navier-Stokes equations for the supersonic flow past a nonaxisymmetric nozzle. The four numerical algorithms evaluated are based on upwind differencing and feature time-dependent, space relaxation, and parabolized schemes. The three different algebraic turbulence models are the Baldwin-Lomax (1978), the nonequilibrium model of Johnson and King (1985), and the Goldberg (1989) modification for the separated flow region. The present investigation was conducted at freestream Mach numbers of 1.2 and 1.3 and at an angle of attack of 0.0 deg. The Reynolds numbers for the investigation ranged from 20.5 million to 21.5 million, based on the model length. The calculations are compared with experimental data.

Abdol-Hamid, Khaled S.↗

On the use of implicit line-relaxation and multi-zonal computations

The use of an implicit Gauss-Seidel line-relaxation scheme to solve unfactored equations resulting from a finite-difference upwind approach in solutions of the governing equations of fluid flow is discussed. It is shown that the use of boundary conditions to pass information between zonal boundaries can dramatically improve convergence in steady-state calculations. This is due to the iterative nature of the current implementation of the line-relaxation procedure.

Rogers, Stuart E.↗

Three-dimensional unstructured grid Euler computations using a fully-implicit, upwind method

A method has been developed to solve the Euler equations on a three-dimensional unstructured grid composed of tetrahedra. The method uses an upwind flow solver with a linearized, backward-Euler time integration scheme. Each time step results in a sparse linear system of equations which is solved by an iterative, sparse matrix solver. Local-time stepping, switched evolution relaxation (SER), preconditioning and reuse of the Jacobian are employed to accelerate the convergence rate. Implicit boundary conditions were found to be extremely important for fast convergence. Numerical experiments have shown that convergence rates comparable to that of a multigrid, central-difference scheme are achievable on the same mesh. Results are presented for several grids about an ONERA M6 wing.

Whitaker, David L.↗

Advanced Computational Aeroacoustics Methods for Fan Noise Prediction

Direct computation of fan noise is presently not possible. One of the major difficulties is the geometrical complexity of the problem. In the case of fan noise, the blade geometry is critical to the loading on the blade and hence the intensity of the radiated noise. The precise geometry must be incorporated into the computation. In computational fluid dynamics (CFD), there are two general ways to handle problems with complex geometry. One way is to use unstructured grids. The other is to use body fitted overset grids. In the overset grid method, accurate data transfer is of utmost importance. For acoustic computation, it is not clear that the currently used data transfer methods are sufficiently accurate as not to contaminate the very small amplitude acoustic disturbances. In CFD, low order schemes are, invariably, used in conjunction with unstructured grids. However, low order schemes are known to be numerically dispersive and dissipative. dissipative errors are extremely undesirable for acoustic wave problems. The objective of this project is to develop a high order unstructured grid Dispersion-Relation-Preserving (DRP) scheme. would minimize numerical dispersion and dissipation errors. contains the results of the funded portion of the project. scheme on an unstructured grid has been developed. constructed in the wave number space. The characteristics of the scheme can be improved by the inclusion of additional constraints. Stability of the scheme has been investigated. Stability can be improved by adopting the upwinding strategy.

Envia, Edmane↗

Evaluation of a Multigrid Scheme for the Incompressible Navier-Stokes Equations

A fast multigrid solver for the steady, incompressible Navier-Stokes equations is presented. The multigrid solver is based upon a factorizable discrete scheme for the velocity-pressure form of the Navier-Stokes equations. This scheme correctly distinguishes between the advection-diffusion and elliptic parts of the operator, allowing efficient smoothers to be constructed. To evaluate the multigrid algorithm, solutions are computed for flow over a flat plate, parabola, and a Karman-Trefftz airfoil. Both nonlifting and lifting airfoil flows are considered, with a Reynolds number range of 200 to 800. Convergence and accuracy of the algorithm are discussed. Using Gauss-Seidel line relaxation in alternating directions, multigrid convergence behavior approaching that of O(N) methods is achieved. The computational efficiency of the numerical scheme is compared with that of Runge-Kutta and implicit upwind based multigrid methods.

Swanson, R. C.↗

A solution adaptive grid procedure for an upwind parabolized flow solver

A solution adaptive grid procedure based on an error equi-distribution scheme is developed and applied to a Parabolized Navier-Stokes solver. An improved method for selecting weighting functions is introduced which involves normalizing a combination of flowfield gradients and curvature of a number of dependent variables and then selecting the largest at each point. The scheme re-distributes grid points line-by-line, with grid point motion controlled by forces analogous to tensional and torsional spring forces with the spring constants set equal to the weighting functions. Torsional terms are functions of the grid point positions along neighboring grid lines and provide grid smoothness and stability. A grid-fitting scheme is introduced for external flows in which the number of grid points in the freestream are reduced to a minimum. Results for several problems are presented to demonstrate the improvements obtainable with the solution adaptive grid procedure.

Harvey, Albert D.↗

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗

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.↗