Engineering topics
Publications and source records attributed to Maccormack, R. W..
A method for simulating the flow of a nitrogen plasma in thermodynamic nonequilibrium is described. The method relies on the modeling of the relevant physical and chemical phenomena and on a fully coupled implicit finite volume upwind numerical approximation of the governing equations. The results agree with experimental data.
A computer code for calculating flow about a circulation control airfoil within a wind tunnel test section was developed. This code is being validated for eventual use as an aid to design such airfoils. The concept of code validation being used is explained. The initial stages of the process were accomplished. The present code was applied to a low subsonic, 2-D flow about a circulation control airfoil for which extensive data exist. Two basic turbulence models and variants thereof were successfully introduced into the algorithm, the Baldwin-Lomax algebraic and the Jones-Launder two equation models of turbulence. The variants include adding a history of the jet development for the algebraic model and adding streamwise curvature effects for both models. Numerical difficulties and difficulties in the validation process are discussed. Turbulence model and code improvements to proceed with the validation process are also discussed.
A computer code for calculating flow about a circulation control airfoil within a wind tunnel test section has been developed. This code is being validated for eventual use as an aid to design such airfoils. The concept of code validation being used is explained. The initial stages of the process have been accomplished. The present code has been applied to a low-subsonic, 2-D flow about a circulation control airfoil for which extensive data exist. Two basic turbulence models and variants thereof have been successfully introduced into the algorithm, the Baldwin-Lomax algebraic and the Jones-Launder two-equation models of turbulence. The variants include adding a history of the jet development for the algebraic model and adding streamwise curvature effects for both models. Numerical difficulties and difficulties in the validation process are discussed. Turbulence model and code improvements to proceed with the validation process are also discussed.
The difficulties of constructing efficient algorithms for three-dimensional flow are discussed. Reasonable candidates are analyzed and tested, and most are found to have obvious shortcomings. Yet, there is promise that an efficient class of algorithms exist between the severely time-step sized-limited explicit or approximately factored algorithms and the computationally intensive direct inversion of large sparse matrices by Gaussian elimination.
Explore the source record for details and available documents.
Implicit numerical procedures have been developed to solve the Reynold's Averaged Navier-Stokes equations for the flowfield within a two-dimensional thrust reversing engine nozzle. The procedures have been tested on four different thrust reverser nozzle geometries with both fixed and moving walls. The results have compared well with available experimental data.
The present work extends the recently reported implicit analogue of MacCormack's earlier widely-used explicit method to external axisymmetric laminar flows with strong entropy gradients. The details of the 'numerics' of the implicit part are provided in a body-oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analogue may appeal to researchers for using it in the analysis of three-dimensional flow problems.
An implicit analogue of a widely used explicit method to external axisymmetric laminar flows with strong entropy gradients is extended. The details of the "numerics" of the implicit part are provided in a body oriented coordinate system with a moving outer (shock) boundary during the transient part of the solutions. The limiting values of the Courant number are obtained when the shock boundary is treated explicitly. The solution algorithm outlined includes the treatment of the source term associated with the equations in weak conservation form. From the results obtained for two sample problems, it becomes clear that accuracy of predictions is, indeed, very good at higher values of the Courant number. There is a significant saving in overall computing time, depending on the Courant number used and the flow Reynolds number. These properties combined with the simplicity of programming the implicit analog may appeal to researchers for using it in the analysis of 3-D flow problems.
The explicit-implicit predictor corrector method of MacCormack (1981) is applied to the analysis of flows past airfoils. By comparing results obtained with different methods and meshes, it is shown that the above method provides, after certain modifications, reasonably good predictions of inviscid and viscous flows about an airfoil. Good results are also obtained for the transonic regime if the free-stream conditions are correct and if a suitable mesh is used.
A new numerical method which was used to reduce the computation time required in fluid dynamics to solve the Navier-Stokes equations at flight Reynolds numbers is described. The method is the implicit analogue of the explicit finite different method. It uses this as its first stage, while the second stage removes the restrictive stability condition by recasting the difference equations in an implicit form. The resulting matrix equations to be solved are either upper or lower block bidiagonal equations. The new method makes it possible and practical to calculate many important three dimensional, high Reynolds number flow fields on computers.
Although much progress has already been made in solving problems in aerodynamic design, many new developments are still needed before the equations for unsteady compressible viscous flow can be solved routinely. This paper describes one such development. A method for solving these equations has been devised that: (1) is second-order accurate in space and time; (2) is unconditionally stable; (3) preserves conservation form; (4) requires no block or scalar tridiagonal inversions; (5) is simple and straightforward to program (estimated 10% modification for the update of many existing programs); (6) is more efficient than present methods; and (7) should easily adapt to current and future computer architectures. Computational results for laminar and turbulent flows at Reynolds numbers from 300,000 to 3 x 10 to the 7th and at CFL numbers as high as 1000 are compared with theory and experiment.
Topics discussed include polygon transformations in fluid mechanics, computation of three-dimensional horseshoe vortex flow using the Navier-Stokes equations, an improved surface velocity method for transonic finite-volume solutions, transonic flow calculations with higher order finite elements, the numerical calculation of transonic axial turbomachinery flows, and the simultaneous solutions of inviscid flow and boundary layer at transonic speeds. Also considered are analytical solutions for the reflection of unsteady shock waves and relevant numerical tests, reformulation of the method of characteristics for multidimensional flows, direct numerical simulations of turbulent shear flows, the stability and separation of freely interacting boundary layers, computational models of convective motions at fluid interfaces, viscous transonic flow over airfoils, and mixed spectral/finite difference approximations for slightly viscous flows.
A semi-implicit method is applied to solve the Navier-Stokes equations numerically and to evaluate the features of the free-interaction phenomenon that occurs when a shock wave impinges on a Blasius boundary layer. Comparisons are made with predictions of the triple-deck theory and experiment. Results include pressure and skin-friction distribution in the free-interaction region for various values of Reynolds number.
The presented review is concerned with the problem of calculating compressible viscous flows. Basic numerical considerations and problems associated with calculating viscous flows are examined and current numerical approaches toward the solution of the Navier-Stokes equations are discussed. It is pointed out that the numerical solution of the full time-dependent equations for turbulent flow is not practical with present computers. Therefore, turbulence effects must be accounted for by modeling. Developments related to turbulence modeling are described. In connection with a discussion of numerical methods for solving viscous flow equations, attention is given to numerical domains of dependence of typical explicit and implicit methods, the diffusion problem, the convection-diffusion problem, and the split-hybrid method.
The paper presents a three-dimensional analysis of the nonlinear light matter interaction in a hydrodynamic context. It is reported that the resulting equations are a generalization of the Navier-Stokes equations subjected to an internal potential which depends solely upon the fluid density. In addition, three numerical approaches are presented to solve the governing equations using an extension of McCormack predict-corrector scheme. These are a uniform grid, a dynamic rezoned grid, and a splitting technique. It is concluded that the use of adaptive mapping and splitting techniques with McCormack two-level predictor-corrector scheme results in an efficient and reliable code whose storage requirements are modest compared with other second order methods of equal accuracy.
The calculation of flow fields past aircraft configuration at flight Reynolds numbers is considered. Progress in devising accurate and efficient numerical methods, in understanding and modeling the physics of turbulence, and in developing reliable and powerful computer hardware is discussed. Emphasis is placed on efficient solutions to the Navier-Stokes equations.
A rapid numerical scheme is used to solve the complete mass-averaged Navier-Stokes equations for supersonic turbulent flow over a three-dimensional compression corner. A simple eddy viscosity model is developed, and the interaction of a swept shock wave and a three-dimensional turbulent boundary layer is studied. Good agreement is obtained between the present results and experimental measurements for the case of a wedge with an angle of 6 deg on a flat-plate sidewall. For the case of a 12-deg wedge angle, the computed results do not show the existence of a peak pressure found experimentally. However, the range of interaction, the plateau pressure, and the peak heat transfer are closely predicted for all cases. The high heat transfer near the axial corner is due to the thinning of the boundary layer and inflow of fresh high-momentum fluid. The heat transfer is relieved through pressure reduction and boundary-layer thickening.
A review is presented of implicit and hybrid methods applicable to solving viscous flows at high Reynolds numbers. Flows within axisymmetric channels containing stationary shock waves, past blunt-nosed lifting airfoils, past sharp-nosed symmetric airfoils with buffet, past three-dimensional compression ramps with side walls, and past ogive- and hemisphere-cylinders at angle of attack have been examined. Reynolds numbers as high as 10 to the 9th power have been used.