Search NASA⌕ Search

Engineering topics

Hafez, M. M.

Publications and source records attributed to Hafez, M. M..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

3D Structured Grid Adaptation

Grid adaptation for structured meshes is the art of using information from an existing, but poorly resolved, solution to automatically redistribute the grid points in such a way as to improve the resolution in regions of high error, and thus the quality of the solution. This involves: (1) generate a grid vis some standard algorithm, (2) calculate a solution on this grid, (3) adapt the grid to this solution, (4) recalculate the solution on this adapted grid, and (5) repeat steps 3 and 4 to satisfaction. Steps 3 and 4 can be repeated until some 'optimal' grid is converged to but typically this is not worth the effort and just two or three repeat calculations are necessary. They also may be repeated every 5-10 time steps for unsteady calculations.

Banks, D. W.↗
NASA NTRS1996

Theoretical and Computational Studies on Sonic Boom Propagation and Its Submarine Impact

Sonic boom impact on the environment of populated area and habitat has been a major concern for the design, operation, and program planning of super/hypersonic vehicles as well as sE launch. Recent development in sonic boom studies reviewed has made evident need for amending the near-field analysis with nonlinear calculations, but an unambiguous matching procedure to assure waveform-prediction improvement is still lacking. Another problem receiving much attention recently is the renewed issue with 'transition focus booms' particularly the 'superboom' which occurs during a speed change through a threshold Mach number and gives rise to strong wave-focussing effects; however, its intensity and the extent of the impact area not be established from existing methods. A third aspect of a more recent concern is the potential sonic boom impact on pelagic and coastal environment, of which the methodology for defining impact has yet to be adequately developed. The study addresses these three aspects in the frame work of a wave-field analysis for a stratified atmosphere, employing coordinates fixed to the vehicle in steady horizontal motion.

Cheng, H. K.↗
NASA NTRS1994

Hyperbolic Prismatic Grid Generation and Solution of Euler Equations on Prismatic Grids

A hyperbolic grid generation method is used to generate prismatic grids and an approach using prismatic grids to solve the Euler equations is presented. The theory of the stability and feasibility of the hyperbolic grid generation method is presented. The hyperbolic grid generation method of Steger et al for structured grids is applied to a three dimensional triangularized surface definition to generate a grid that is unstructured on each successive layer. The grid, however, retains structure in the body-normal direction and has a computational cell shaped like a triangular prism. In order to take advantage of the structure in the normal direction, a finite-volume scheme that treats the unknowns along the normal direction implicitly is introduced and the flow over a sphere is simulated.

Pandya, S. A.↗
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

Vortex breakdown simulation based on a nonlinear inviscid model

It is shown that the inviscid equations governing steady axisymmetric flow with swirl, admit solutions with closed streamlines. Results are obtained using two different numerical algorithms. The first is based on a multigrid method for nonlinear eigenvalue problems, while the second is based on a least squares formulation.

Hafez, M. M.↗
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

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

Stability analysis of intermediate boundary conditions in approximate factorization schemes

The paper discusses the role of the intermediate boundary condition in the AF2 scheme used by Holst for simulation of the transonic full potential equation. It is shown that the treatment suggested by Holst led to a restriction on the time step and ways to overcome this restriction are suggested. The discussion is based on the theory developed by Gustafsson, Kreiss, and Sundstrom and also on the von Neumann method.

South, J. C., Jr.↗
NASA NTRS1985

A modified least squares formulation for a system of first-order equations

Second order equations in terms of auxiliary variables similar to potential and stream functions are obtained by applying a weighted least squares formulation to a first order system. The additional boundary conditions which are necessary to solve the higher order equations are determined and numerical results are presented for the Cauchy-Riemann equations.

Hafez, M. M.↗
NASA NTRS1985

Some iterative schemes for transonic potential flows

The minimal residual (MR) method for the numerical solution of transonic potential flows is closely related to the conjugate gradient method, which has found widespread use in the solution of large sparse, symmetric, and positive-definite linear equations. The primary advantage of the MR method is its applicability to both symmetric and nonsymmetric matrices.

Wong, Y. S.↗
NASA NTRS1985

Numerical algorithms for transonic inviscid flow calculations

It is pointed out that the supercritical wing is one of the most important features of modern transonic aerodynamics. The design of its shock free airfoil section depends on potential flow calculations. The present paper is concerned with the development of inviscid flow simulation methods based on potential formulations, taking into account also the problem of nonuniqueness of the potential solution. Nonisentropic and nonisoenergetic models are considered, and an alternative approach using the stream function is discussed. Attention is given to transonic small disturbance calculations, calculations based on the full potential equation, iterative methods, wave drag calculations, and an alternative form of Euler equations.

Hafez, M. M.↗
NASA NTRS1984

Calculation of quasi-one-dimensional flows with shocks

The two-point subsonic method of Wornom (1983) is modified to permit calculation of transonic flows with shocks using artificial density and/or artificial pressure. This method, which requires no numerical boundary conditions, is then used to calculate transonic flow for large and small CFL (Courant-Friedrichs-Lewy) numbers. Convergence acceleration is achieved for small CFL numbers by using multi-grid.

Wornom, S. F.↗
NASA NTRS1984

A modified least squares formulation for a system of first order equations

Second order equations in terms of auxiliary variables similar to potential and stream functions are obtained by applying a weighted least squares formulation to a first order system. The additional boundary conditions which are necessary to solve the higher order equations are determined and numerical results are presented for the Cauchy-Riemann equations.

Hafez, M. M.↗
NASA NTRS1984

Entropy condition satisfying approximations for the full potential equations 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

A rule for selecting analytical boundary conditions for the conservative quasi-one-dimensional nozzle flow equations

For the one-dimensional nozzle flow equations, the number of analytical boundary conditions that may be applied corresponds to the number of exterior characteristics reaching the boundary points. Since the slopes of the characteristics are the characteristic speeds u, u+c, and u-c, one concludes that two analytical boundary conditions are required at a subsonic inflow boundary and one at a subsonic outflow boundary. However, using this guideline alone can lead to nonunique solutions in subsonic cases. This paper presents a rule for selecting analytical boundary conditions for the dependent variables which will yield unique solutions if they exist.

Wornom, S. F.↗
NASA NTRS1983

Stability analysis of intermediate boundary conditions in approximate factorization schemes

In many cases, approximate factorization schemes have provided a significant increase in efficiency over previously used solution methods in certain problems. The present investigation is concerned with the importance of intermediate boundary conditions in approximate factorization schemes, taking into account a specific example regarding a boundary-induced stability restriction in a scheme for the transonic full-potential equation. The considered scheme has been discussed by Holst (1979). Holst's scheme is a variation of the AF2 schemes described by Ballhaus and Steger (1975). The application of the AF2 scheme to the two-dimensional Laplace equation in a rectangle is studied giving attention to the stability of the scheme in connection with various boundary conditions for the intermediate variable.

South, J. C., Jr.↗
NASA NTRS1982

A minimal residual method for transonic potential flows

For transonic flow calculations, a combination of the successive line over-relaxation (SLOR) and the preconditioned conjugate gradient (CG) method has been suggested by Wong and Hafez (1981). This paper studies the method of minimal residual (MR) which avoids a combined iteration. This method is closely related to the CG method, may be regarded as a first-order gradient method, and is applicable to symmetric and nonsymmetric matrices. The problem is formulated as a nonlinear mixed elliptic-hyperbolic partial differential equation which includes an artificial viscosity and a switching function which is zero in subsonic regions and nonzero in supersonic regions. Alternatives to the SLOR method which provide faster convergence rates are introduced. The preconditioned MR algorithm is developed, and transonic potential flows around NACA 0012 airfoil are calculated for different Mach numbers and angles of attack. Preliminary results are presented, demonstrating that the MR algorithm requires no parameter estimation and rapidly converges for subsonic flows.

Wong, Y. S.↗