Engineering topics

Hafez, M.

Publications and source records attributed to Hafez, M..

29 records · Page 2

NASA NTRS1985

Convergence acceleration of iterative solutions of Euler equations for transonic flow computations

In this paper, two acceleration techniques for Euler calculations are investigated. The first technique is an extrapolation procedure based on the Power Method; it is applicable when the iterative matrix has dominant eigenvalues. Both real and complex conjugate roots are allowed. The second technique is a generalization of the Minimal Residual Method, where the extrapolation step consists of a weighted combination of the corrections at different iteration levels and the weights are chosen to minimize the Lz norm of the residual. Numerical results, using Jameson's Runge-Kutta Multigrid Code, are presented. The extra computational work to apply either technique is negligible and the extra storage is not a problem on current supercomputers.

Hafez, M.↗

NASA NTRS1985

Entropy condition satisfying approximations for the full potential equation of transonic flow

A class of conservative difference approximations for the steady full potential equation was presented. They are, in general, easier to program than the usual density biasing algorithms, and in fact, differ only slightly from them. Rigorous proof indicated that these new schemes satisfied a new discrete entropy inequality, which ruled out expansion shocks, and that they have sharp, steady, discrete shocks. A key tool in the analysis is the construction of a new entropy inequality for the full potential equation itself. Results of some numerical experiments using the new schemes are presented.

Osher, S.↗

NASA NTRS1984

Calculations of separation bubbles using boundary-layer-type equations. I. II

Two-dimensional models of separated flows are presented with stream-function boundary equations. A patching procedure is employed to solve for unknowns in viscous-inviscid interactions simultaneously in order to derive strong coupling. The technique features simplified Navier-Stokes equations and is useful for flows around flat plates and more complex geometries such as airfoils and cascades. Sample calculations are provided for separated flows in diffusers and in the neighborhood of trailing edges. Massive bubbles and unclosed bubbles are noted to present difficulties for the models.

Halim, A.↗

NASA NTRS1984

Improved finite difference schemes for transonic potential calculations

Engquist and Osher (1980) have introduced a finite difference scheme for solving the transonic small disturbance equation, taking into account cases in which only compression shocks are admitted. Osher et al. (1983) studied a class of schemes for the full potential equation. It is proved that these schemes satisfy a new discrete 'entropy inequality' which rules out expansion shocks. However, the conducted analysis is restricted to steady two-dimensional flows. The present investigation is concerned with the adoption of a heuristic approach. The full potential equation in conservation form is solved with the aid of a modified artificial density method, based on flux biasing. It is shown that, with the current scheme, expansion shocks are not possible.

Hafez, M.↗

NASA NTRS1983

Improved relaxation schemes for transonic potential calculations

A block relaxation scheme, grouped in a red-black ordering, is applied to transonic airfoil calculations using body fitted coordinates. The scheme is simple and is easily vectorizable. Detailed comparisons with Approximate Factorization Method (AF2) are presented and it is shown that the improved relaxation scheme is competitive in all cases considered. Transonic results, of engineering accuracy, on an 0-type grid of 149 x 30 points, are ususally obtained within two hundred iterations (approximately 40 seconds on Cyber 175).

Hafez, M.↗

NASA NTRS1983

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 formlations are presented together with the details of the treatment of shocks and wakes, and drag and lift calculations. The nonuniqueness problem is studied using different artificial viscosity forms. Numerical results are compared to Euler solutions.

Hafez, M.↗

NASA NTRS1983

Transonic small disturbance calculations including entropy corrections

Murman's fully conservative mixed type finite-difference operators are first modified. A special sonic point operator with an iterative damping term is introduced which helps the convergence and does not affect the spatial conservative differences. Reliable calculations with second order supersonic schemes are obtained using two sonic operators, the regular sonic point operator followed by a first order supersonic scheme. Also, shock point operator is shown to be equivalent to fitting a locally normal shock terminating the supersonic region. The potential calculations are then modified to account for the non-isentropic jump conditions using a simple shock fitting procedure based on Prandtl relation. The entropy increase across the shock is calculated in terms of the Mach number upstream of the shock and the effect of the generated vorticity is estimated via Crocco relation. Different examples are calculated and extensions to the full potential equation are discussed.

Hafez, M.↗

NASA NTRS1982

Transonic wind tunnel wall interference corrections for three-dimensional models

A procedure for the evaluation of wall interference corrections for three-dimensional models is presented. The Mach number and angle-of-attack corrections require the numerical solution of the potential equation about a simplified representation of the experimental model. Pressure measurements are required near the wind tunnel walls. The correction procedure also requires knowledge of the free-stream Mach number, the model angle of attack, and the lift force experienced by the model. The procedure provides an estimate of the accuracy of the correction. For slender configurations at Mach numbers close to one, the Equivalence Rule formulation is adopted to calculate the wall interference effects. Preliminary results are presented for both general and slender-body configurations.

Rizk, M. H.↗

NASA NTRS1981

Perturbation of transonic flow with shocks

A general formulation of the perturbation problem is studied, and a new approach, perturbation sequence expansion, is introduced for handling shock disturbances. The method is applied to unsteady effects, three-dimensional corrections to axisymmetric and two-dimensional flows, and wind tunnel corrections. The perturbation equations are nonlinear and can be solved by shock capturing methods.

Hafez, M.↗

NASA NTRS1981

Numerical solution of transonic stream function equation

The stream function equation, in conservation form, looks similar to the full potential equation and existing methods (e.g. artificial compressibility) can be readily applied. Rotational flows can be calculated once the vorticity (due to shocks or nonuniformity) is evaluated. There are, however, two main difficulties: First, the density is not uniquely determined in terms of the flux (there are two solutions; the subsonic and the supersonic branch with a square root singularity at the sonic point). Methods to overcome this difficulty are studied and results are presented with some remarks on inviscid separation and closed stream lines. Second, the need of two stream functions for three dimensional calculations is briefly discussed.

Hafez, M.↗

NASA NTRS1979

Artificial compressibility methods for numerical solutions of transonic full potential equation

New methods for transonic flow computations based on the full potential equation in conservation form are presented. The idea is to modify slightly the density (due to the artificial viscosity in the supersonic region), and solve the resulting elliptic-like problem iteratively. It is shown that standard discretization techniques (central differencing) as well as some standard iterative procedures (SOR, ADI, and explicit methods) are applicable to the modified transonic mixed-type equation. Calculations of transonic flows around cylinders and airfoils are discussed with special emphasis on the explicit methods that are suitable for vector processing on the STAR 100 computer.

Hafez, M.↗