SEARCH · Search NASA

Results for “Artificial viscosity”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Preconditioning for first-order spectral discretization

Efficient solution of the equations from spectral discretizations is essential if the high-order accuracy of these methods is to be realized. Direct solution of these equations is rarely feasible, thus iterative techniques are required. A preconditioning scheme for first-order Chebyshev collocation operators is proposed herein, in which the central finite difference mesh is finer than the collocation mesh. Details of the proper techniques for transferring information between the meshes are given here, and the scheme is analyzed by examination of the eigenvalue spectra of the preconditioned operators. The effect of artificial viscosity required in the inversion of the finite difference operator is examined. A second preconditioning scheme, involving a high-order upwind finite difference operator of the van Leer type is also analyzed to provide a comparison with the present scheme. Finally, the performance of the present scheme is verified by application to several test problems.

Streett, C. L.↗

NASA NTRS1986

Euler solutions for the flow around a hovering helicopter rotor

A method for computing the flow around a hovering helicopter rotor is presented. The solution procedure is divided into two parts. The first part is a finite volume solution of the Euler equations for the near field flow around a rotor blade. The second part is a free wake approach for determining the wake geometry and induced velocities. These two parts are solved in a coupled fashion. The rolled-up vortex wake of the rotor is included in the Euler computational domain through a nonlinear perturbation technique. This eliminates the numerical diffusion of the wake vorticity due to truncation error and artificial viscosity without the need for excessive grid resolution near the vortex core. The method is used to compute both hovering rotor and wing/vortex interaction flows, and the results are compared to experiment.

Roberts, T. W.↗

NASA NTRS1986

Computation of turbulent pipe and duct flow using third order upwind scheme

The fully developed turbulence in a circular pipe and in a square duct is simulated directly without using turbulence models in the Navier-Stokes equations. The utilized method employs a third-order upwind scheme for the approximation to the nonlinear term and the second-order Adams-Bashforth method for the time derivative in the Navier-Stokes equation. The computational results appear to capture the large-scale turbulent structures at least qualitatively. The significance of the artificial viscosity inherent in the present scheme is discussed.

Kawamura, T.↗

NASA NTRS1986

An entropy correction method for unsteady full potential flows with strong shocks

An entropy correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to account for entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Comparisons between the present method and solutions of the Euler equations and between the present method and experimental data are presented. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation.

Whitlow, W., Jr.↗

NASA NTRS1986

Viscous wing theory development. Volume 1: Analysis, method and results

Viscous transonic flows at large Reynolds numbers over 3-D wings were analyzed using a zonal viscid-inviscid interaction approach. A new numerical AFZ scheme was developed in conjunction with the finite volume formulation for the solution of the inviscid full-potential equation. A special far-field asymptotic boundary condition was developed and a second-order artificial viscosity included for an improved inviscid solution methodology. The integral method was used for the laminar/turbulent boundary layer and 3-D viscous wake calculation. The interaction calculation included the coupling conditions of the source flux due to the wing surface boundary layer, the flux jump due to the viscous wake, and the wake curvature effect. A method was also devised incorporating the 2-D trailing edge strong interaction solution for the normal pressure correction near the trailing edge region. A fully automated computer program was developed to perform the proposed method with one scalar version to be used on an IBM-3081 and two vectorized versions on Cray-1 and Cyber-205 computers.

Chow, R. R.↗

NASA NTRS1987

Relaxation algorithms for the Euler equations

The alternating direction implicit central-difference scheme algorithms may be replaced by successive-line-relaxation (SLR) procedures that are stable in three dimensions. Several Beamand Warming-type codes are presently modified on the basis of SLR solution procedures in order to compute the flows over a cylinder, a NACA-0012 airfoil, and a shock wave reflecting off a plate. Since the codes are written in a 'delta form', the implementation of relaxation procedures that rely only on the corrections is relatively easy; the steady state residual computation and any artificial viscosity terms may be left unchanged.

Wornom, S. F.↗

NASA NTRS1987

Free-wake analysis of a rotor in hover

A numerical method based on the axisymmetric, incompressible Navier-Stokes equations is combined with a lifting surface code to predict the vortex wake of hovering rotors. The lifting surface code, AMI Hover, is used to obtain the circulation distribution on the blade. This circulation distribution is fed into the Navier-Stokes code to compute the vortex wake under this specified circulation distribution. An iteration approach is used between these two codes to converge the circulation distribution and the shape of the vortex wake. A relaxation scheme is developed to resolve the instability encountered among the tip vortices. A reconcentration scheme is used to solve the diffusion problem due to the strong artificial viscosity. The results from the present method are compared with experimental data obtained by smoke-flow visualization and hot-wire measurements for several rotor blade configurations. The comparisons show that the present method is able to predict the complex wake system shed by a hovering rotor.

Chen, C. S.↗

NASA NTRS1987

Application of Runge Kutta time marching scheme for the computation of transonic flows in turbomachines

Numerical solutions of the unsteady Euler equations are obtained using the classical fourth order Runge Kutta time marching scheme. This method is fully explicit and is applied to the governing equations in the finite volume, conservation law form. In order to determine the efficiency of this scheme for solving turbomachinery flows, steady blade-to-blade solutions are obtained for compressor and turbine cascades under subsonic and transonic flow conditions. Computed results are compared with other numerical methods and wind tunnel measurements. The study also focuses on other important numerical aspects influencing the performance of the algorithm and the solution accuracy such as grid types, boundary conditions and artificial viscosity. For this purpose, H, O, and C type computational grids as well as characteristic and extrapolation type boundary conditions are included in solution procedures.

Subramanian, S. V.↗

NASA NTRS1986

Computational fluid dynamics overview

The hierarchy of codes; boundary layers and turbulence models; uniqueness and artificial viscosity; efficiency of solution algorithms; data display and analysis; and cost of experimental validation are briefly examined.

Bobbitt, Percy J.↗

NASA NTRS1988

Application of a nonisentropic full potential method to AGARD standard airfoils

An entropy-correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to model the entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Calculated results are presented for the NLR 7301, NACA 0012, and NACA 64A010A airfoils. Comparisons of the present method and solutions of the Euler equations are presented for the NLR 7301 airfoil, and comparisons of the present method and experimental data are presented for all three airfoils. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation. In addition, it is shown that modeling shock-generated entropy extends the range of validity of the full potential method.

Whitlow, Woodrow, Jr.↗

NASA NTRS1988

Euler equation computations for the flow over a hovering helicopter rotor

A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.

Roberts, Thomas Wesley↗

NASA NTRS1988

Application of a nonisentropic full potential method to AGARD standard airfoils

Whitlow, Woodrow, Jr.↗

NASA NTRS1988

Entropy and vorticity corrections for transonic flows

Different models for inviscid transonic flows are examined. The common assumptions that the flow is isentropic and irrotational are critically evaluated. Entropy and vorticity correction procedures for potential and stream function formulations are presented, together with the details of the treatment of shocks and wakes, and drag and lift calculations. The non-uniqueness problem of the potential formulation is studied using different artificial viscosity forms. Numerical results are compared with Euler solutions.

Hafez, M.↗

NASA NTRS1987

A spectral multi-domain technique for analysis of chemically reacting flows

A technique is presented which models the chemical kinetics and flow kinematics of a nonionized air mixture passing through a fully resolved shock wave, thus alleviating the need for artificial viscosity. The quasi-one-dimensional Navier-Stokes equations and the species conversation equation are the governing equations. The validity of the multidomain Navier-Stokes algorithm is demonstrated via comparison with experiment.

Macaraeg, Michele G.↗

NASA NTRS1987

Simulation of separated flow past a bluff body using Navier-Stokes equations

Two-dimensional flow past a bluff body is presently simulated on the basis of an analysis that employs the incompressible, unsteady Navier-Stokes equations in terms of vorticity and stream function. The fully implicit, time-marching, alternating-direction, implicit-block Gaussian elimination used is a direct method with second-order spatial accuracy; this allows it to avoid the introduction of any artificial viscosity. Attention is given to the simulation of flow past a circular cylinder with and without symmetry, requiring the use of either the half or the full cylinder, respectively.

Ghia, K. N.↗

NASA NTRS1988

An adaptive pseudospectral method for discontinuous problems

The accuracy of adaptively chosen, mapped polynomial approximations is studied for functions with steep gradients or discontinuities. It is shown that, for steep gradient functions, one can obtain spectral accuracy in the original coordinate system by using polynomial approximations in a transformed coordinate system with substantially fewer collocation points than are necessary using polynomial expansion directly in the original, physical, coordinate system. It is also shown that one can avoid the usual Gibbs oscillation associated with steep gradient solutions of hyperbolic pde's by approximation in suitably chosen coordinate systems. Continuous, high gradient solutions are computed with spectral accuracy (as measured in the physical coordinate system). Discontinuous solutions associated with nonlinear hyperbolic equations can be accurately computed by using an artificial viscosity chosen to smooth out the solution in the mapped, computational domain. Thus, shocks can be effectively resolved on a scale that is subgrid to the resolution available with collocation only in the physical domain. Examples with Fourier and Chebyshev collocation are given.

Augenbaum, Jeffrey M.↗

NASA NTRS1988

Simulation of 2-dimensional viscous flow through cascades using a semi-elliptic analysis and hybrid C-H grids

A semi-elliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is neglible but which are significantly influenced by upstream interactions. The IPNS equations are obtained from the Navier-Stokes equations by dropping the streamwise viscous-diffusion terms but retaining upstream influence via the streamwise pressure-gradient. A two-step alternating-direction-explicit numerical scheme is developed to solve these equations. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Also, solutions to compressible as well as nearly compressible flows are obtained without any modification either in the analysis or in the solution process. The procedure is applied to constricted channels and cascade passages formed by airfoils of various shapes. These geometries are represented using numerically generated curilinear boundary-oriented coordinates forming an H-grid. A hybrid C-H grid, more appropriate for cascade of airfoils with rounded leading edges, was also developed. Satisfactory results are obtained for flows through cascades of Joukowski airfoils.

Ramamurti, R.↗

NASA NTRS1989

An embedded mesh procedure for leading-edge vortex flows

A cell-vertex scheme is outlined for solving the flow about a delta wing with M (sub infinity) is greater than 1. Embedded regions of mesh refinement allow solutions to be obtained which have much higher resolution than those achieved to date. Effects of mesh refinement and artificial viscosity on the solutions are studied, to determine at what point leading-edge vortex solutions are grid-converged. A macroscale and a microscale for the size of the vortex are defined, and it is shown that the macroscale (which includes the wing surface properties) is converged on a moderately refined grid, while the microscale is very sensitive to grid spacing. The level of numerical diffusion in the core of the vortex is found to be substantial. Comparisons with the experiment are made for two cases which have transonic cross-flow velocities.

Powell, Kenneth G.↗