SEARCH · Search NASA
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.
The artificial viscosity method of Kuwahara and Takami (1973) is used to calculate the roll-up of trailing vortices behind a number of practical aerodynamic configurations. Where possible, the results are compared for core location with available experimental data.
Artificial viscosity is added either implicitly or explicitly in practically every numerical scheme for suppressing spurious oscillations in the solution of fluid dynamics equations. In the present central-difference scheme, artificial viscosity is added explicitly for suppressing high frequency oscillations and achieving good convergence properties. The amount of artificial viscosity added is controlled through the use of pre-selected coefficients. In the standard scheme, scalar coefficients based on the spectral radii of the Jacobian of convective fluxes are used. However, this can add too much viscosity to the slower waves. Hence, use the matrix-valued coefficients, which give appropriate viscosity for each wave component is suggested. With the matrix-valued coefficients, the central-difference scheme produces more accurate solutions on a given grid, particularly in the vicinity of shocks, while still maintaining good convergence properties.
Artificial viscosity is added either implicity or explicitly in practically every numerical scheme for suppressing spurious oscillations in the solution of fluid-dynamics equations. In the present central-difference scheme, artificial viscosity is added explicitly for suppressing high-frequency oscillations and achieving good convergence properties. The amount of artificial viscosity added is controlled through the use of preselected coefficients. In the standard scheme, scalar coefficients based on the spectral radii of the Jacobian of the convective fluxes are used. However, this can add too much viscosity to the slower waves. Hence, the use of matrix-valued coefficients, which give appropriate viscosity for each wave component, is suggested. With the matrix-valued coefficients, the central-difference scheme produces more accurate solutions on a given grid, particularly in the vicinity of shocks and boundary layers, while still maintaining good convergence properties.
A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over the airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small perturbation about a steady inviscid flow over a NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 degrees. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time marching scheme are also presented. This procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of terms for artificial viscosity.
A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small time-dependent perturbation about a steady inviscid flow over an NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 deg. The method has been applied to inviscid flows here, but as discussed is also applicable to viscous flows. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundaries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time-marching scheme are also presented. The eigenvalue procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of artificial viscosity.
Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.
Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.
A finite volume method is used to calculate compressible inviscid flows over blunt bodies using, in general, unstructured grids. Artificial viscosity forms are derived based on a simplified least squares procedure. The extra second order terms are consistent with the governing equations, hence a systematic treatment of the numerical boundary conditions can be easily implemented. A special treatment of blunt bodies may be required. The discrete equations are linearized and the resulting system is solved by a relaxation method. Preliminary results indicate that the effect of the numerical dissipation is minimal. For subsonic flows over smooth bodies, the solution is practically vorticity-free and the total pressure loss is of the same order as the truncation error. Finally, some extensions of the present method are briefly discussed.
The total variation diminishing (TVD) finite difference scheme can be interpreted as a Lax-Wendroff scheme plus an upwind weighted artificial dissipation term. If a particular flux limiter is chosen and the requirement for upwind weighting is removed, an artificial dissipation term which is based on the theory of TVD schemes is obtained which does not contain any problem dependent parameters and which can be added to existing MacCormack method codes. Numerical experiments to examine the performance of this new method are discussed.
The computational fluid dynamics code, PARC3D, is tested to see if its use of non-physical artificial dissipation affects the accuracy of its results. This is accomplished by simulating a shock-laminar boundary layer interaction and several hypersonic flight conditions of the Pegasus(TM) launch vehicle using full artificial dissipation, low artificial dissipation, and the Engquist filter. Before the filter is applied to the PARC3D code, it is validated in one-dimensional and two-dimensional form in a MacCormack scheme against the Riemann and convergent duct problem. For this explicit scheme, the filter shows great improvements in accuracy and computational time as opposed to the nonfiltered solutions. However, for the implicit PARC3D code it is found that the best estimate of the Pegasus experimental heat fluxes and surface pressures is the simulation utilizing low artificial dissipation and no filter. The filter does improve accuracy over the artificially dissipative case but at a computational expense greater than that achieved by the low artificial dissipation case which has no computational time penalty and shows better results. For the shock-boundary layer simulation, the filter does well in terms of accuracy for a strong impingement shock but not as well for weaker shock strengths. Furthermore, for the latter problem the filter reduces the required computational time to convergence by 18.7 percent.
Wall functions, as used in the typical high Reynolds number k-epsilon turbulence model, can be implemented in various ways. A least disruptive method (to the flow solver) is to directly solve for the flow variables at the grid point next to the wall while prescribing the values of k and epsilon. For the centrally-differenced finite-difference scheme employing artificial viscocity (AV) as a stabilizing mechanism, this methodology proved to be totally useless. This is because the AV gives rise to a large error at the wall due to too steep a velocity gradient resulting from the use of a coarse grid as required by the wall function methodology. This error can be eliminated simply by extrapolating velocities at the wall, instead of using the physical values of the no-slip velocities (i.e. the zero value). The applicability of the technique used in this paper is demonstrated by solving a flow over a flat plate and comparing the results with those of experiments. It was also observed that AV gives rise to a velocity overshoot (about 1 percent) near the edge of the boundary layer. This small velocity error, however, can yield as much as 10 percent error in the momentum thickness. A method which integrates the boundary layer up to only the edge of the boundary (instead of infinity) was proposed and demonstrated to give better results than the standard method.
Analytical and experimental techniques for modeling the aerothermodynamics of hypersonic flight are assessed, together with the problems which will be encountered in developing reusable hypersonic vehicles. Emphasis is placed on a numerical coupling between a nonequilibrium chemistry model and hypersonic flow kinematics. Finite difference and finite element descriptions of flow fields in which molecules encounter a shock wave and undergo various motion (and thereby energy) transformations are discussed. The effects of artificial smearing of the shock wave are considered in terms of the resulting effects on the distribution of the energies and chemical composition of the transition region. Results are provided from schlieren photographs of shock-tube experiments, Navier-Stokes calculations of axisymmetric flow over a conical body, and calculations using a spectral multidomain approach for chemically reacting flows.
A numerical technique is presented for the design of two-dimensional supercritical wing sections with low wave drag. The method is a design mode of the analysis code H which gives excellent agreement with experimental results and is widely used in the aircraft industry. Topics covered include the partial differential equations of transonic flow, the computational procedure and results; the design procedure; a convergence theorem; and description of the code.
We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.
We use central differences to solve the time dependent Euler equations. The schemes are all advanced using a Runge-Kutta formula in time. Near shocks, a second difference is added as an artificial viscosity. This reduces the scheme to a first order upwind scheme at shocks. The switch that is used guarantees that the scheme is locally total variation diminishing (TVD). For steady state problems it is usually advantageous to relax this condition. Then small oscillations do not activate the switches and the convergence to a steady state is improved. To sharpen the shocks, different coefficients are needed for different equations and so a matrix valued dissipation is introduced and compared with the scalar viscosity. The connection between this artificial viscosity and flux limiters is shown. Any flux limiter can be used as the basis of a shock detector for an artificial viscosity. We compare the use of the van Leer, van Albada, mimmod, superbee, and the 'average' flux limiters for this central difference scheme. For time dependent problems, we need to use a small enough time step so that the CFL was less than one even though the scheme was linearly stable for larger time steps. Using a total variation bounded (TVB) Runge-Kutta scheme yields minor improvements in the accuracy.
The application of multigrid relaxation to transonic potential flow calculation was investigated. Conservative potential flows around aerofoils were taken as test problems. The solution algorithm was based on Newton iteration. It was found that the iteration to the circulation has to be kept outside the multigrid algorithm. To obtain meaningful norms of residuals, difference formulas with asymptotic scaling were introduced. Nonlinear instability problems were solved by upwind differencing using mass flux vector splitting instead of artificial viscosity or artificial density. It is suggested that the algorithms efficiency be increased by improving the iteration on the shock positions even though this is a highly nonlinear process.
The stream function equation is solved for steady two-dimensional (and axisymmetric) rotational flows. Both finite differences and finite volumes discretization techniques are studied, using generalized body fitted coordinates and unstructured staggered grids, respectively. For inviscid transonic flows, a new artificial viscosity scheme which does not produce any artificial vorticity is introduced, for the stability of the mixed flow calculations and for capturing shocks. The solution of Euler equations, in primitive variables, are also considered. The effects of the artificial viscosity and numerical boundary conditions on the total enthalpy and the vorticity distributions are demonstrated.
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.