Primitive variable, strongly implicit calculation procedure for viscous flows at all speeds

A coupled solution procedure is described for solving the compressible form of the time-dependent, twodimensional Navier-Stokes equations in body-fitted curvilinear coordinates. This approach employs the strong conservation form of the governing equations but uses primitive variables (M, v, p, T) rather than the more traditional conservative variables (p, pw, pv, et) as unknowns. A coupled modified strongly implicit procedure (CMSIP) is used to efficiently solve the Newton-linearized algebraic equations. It appears that this procedure is effective for Mach numbers ranging from the incompressible limit (Mx ~ 0.01) to supersonic. Generally, smoothing was not needed to control spatial oscillations in pressure for subsonic flows despite the use of central differences. Dual-time stepping was found to further accelerate convergence for steady flows. Sample calculations, including steady and unsteady low-Mach-number internal and external flows and a steady shock-boundarylayer interaction flow, illustrate the capability of the present solution algorithm.

[1]  A. Roshko.  On the development of turbulent wakes from vortex streets , 1953 .

[2]  Leon Trilling,et al.  The Interaction of an Oblique Shock Wave with a Laminar Boundary Layer , 1959 .

[3]  J. F. Osterle,et al.  Finite difference analysis of plane Poiseuille and Couette flow developments , 1961 .

[4]  O. Burggraf.  Analytical and numerical studies of the structure of steady separated flows , 1966, Journal of Fluid Mechanics.

[5]  I. Flügge-Lotz,et al.  The interaction of a shock wave with a laminar boundary layer , 1968 .

[6]  H. L. Stone.  ITERATIVE SOLUTION OF IMPLICIT APPROXIMATIONS OF MULTIDIMENSIONAL PARTIAL DIFFERENTIAL EQUATIONS , 1968 .

[7]  A. A. Amsden,et al.  Numerical calculation of almost incompressible flow , 1968 .

[8]  Numerical Solutions of the Navier-Stokes Equations in Inlet Regions , 1972 .

[9]  R. T. Cheng,et al.  Numerical solution of the viscous flow in the entrance region of parallel plates , 1973 .

[10]  S. G. Rubin,et al.  A diagonally dominant second-order accurate implicit scheme , 1974 .

[11]  W. L. Hankey,et al.  Study of incompressible Navier-Stokes equations in primitive variables using implicit numerical technique , 1977 .

[12]  R. F. Warming,et al.  An Implicit Factored Scheme for the Compressible Navier-Stokes Equations , 1977 .

[13]  B. Fornberg.  A numerical study of steady viscous flow past a circular cylinder , 1980, Journal of Fluid Mechanics.

[14]  Gerry E. Schneider,et al.  A MODIFIED STRONGLY IMPLICIT PROCEDURE FOR THE NUMERICAL SOLUTION OF FIELD PROBLEMS , 1981 .

[15]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[16]  Richard G. Hindman,et al.  Generalized Coordinate Forms of Governing Fluid Equations and Associated Geometrically Induced Errors , 1982 .

[17]  T. Taylor,et al.  Computational methods for fluid flow , 1982 .

[18]  Gene H. Golub,et al.  Matrix computations , 1983 .

[19]  R. Pletcher,et al.  Computational Fluid Mechanics and Heat Transfer. By D. A ANDERSON, J. C. TANNEHILL and R. H. PLETCHER. Hemisphere, 1984. 599 pp. $39.95. , 1986, Journal of Fluid Mechanics.

[20]  Hassan Hassan,et al.  A STRONGLY IMPLICIT PROCEDURE FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS , 1984 .

[21]  Michele Napolitano,et al.  An incremental block-line-Gauss-Seidel method for the Navier-Stokes equations , 1985 .

[22]  J. Piquet,et al.  On the use of several compact methods for the study of unsteady incompressible viscous flow round a circular cylinder , 1985 .

[23]  C. Merkle.  Computation of compressible flows at very low Mach numbers , 1986 .

[24]  M. Visbal,et al.  Evaluation of an implicit Navier-Stokes solver for some unsteady separated flows , 1986 .

[25]  Robert W. Walters,et al.  Upwind relaxation algorithms for the Navier-Stokes equations , 1987 .

[26]  G. D. Raithby,et al.  The segregated approach to predicting viscous compressible fluid flows , 1986 .

[27]  S. Patankar,et al.  Pressure based calculation procedure for viscous flows at all speeds in arbitrary configurations , 1988 .

[28]  Implicit solutions of the incompressible Navier-Stokes equations in primitive variables , 1988 .

[29]  S. Rogers,et al.  An upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations , 1988 .

[30]  COMPARISON OF Tl!O SCHEMES APPLIED TO THE NAVIER-STOKES EQUATIONS , 1989 .

[31]  Meng-Sing Liou,et al.  A Newton/upwind method and numerical study of shock wave/boundary layer interactions , 1989 .

[32]  W. Y. Soh,et al.  Time-Dependent viscous incompressible Navier-Stokes equations: the finite difference Galerkin formulation and streamfunction algorithms , 1989 .

[33]  Stuart E. Rogers,et al.  Upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations , 1990 .

[34]  Karl Gustafson,et al.  Hopf bifurcation in the driven cavity , 1990 .

[35]  Kuo-Huey Chen.  A primitive variable, strongly implicit calculation procedure for two and three-dimensional unsteady viscous flo ws: applications to compressible and incompressible flo ws including flo ws with free surfaces , 1990 .

[36]  A primitive variable, strongly implicit calculation procedure for viscous flows at all speeds , 1990 .

[37]  Coupled space-marching method for the Navier-Stokes equations for subsonic flows , 1991 .

[38]  戸田 保幸.  The Finite Analytic Method , 1992 .

[39]  Simulation of three-dimensional liquid sloshing flows using a strongly implicit calculation procedure , 1993 .

[40]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.