Hostname: page-component-77c89778f8-m42fx Total loading time: 0 Render date: 2024-07-19T05:51:35.221Z Has data issue: false hasContentIssue false
  • English
  • Français

Numerical simulation of incompressible flows within simple boundaries: accuracy

Published online by Cambridge University Press:  29 March 2006

Steven A. Orszag
Steven A. Orszag
Affiliation:
National Center for Atmospheric Research, Boulder, Colorado Permanent address: Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139.
Get access Rights & Permissions [Opens in a new window]

Abstract

Galerkin (spectral) methods for numerical simulation of incompressible flows within simple boundaries are shown to possess many advantages over existing finite-difference methods. In this paper, the accuracy of Galerkin approximations obtained from truncated Fourier expansions is explored. Accuracy of simulation is tested empirically using a simple scalar-convection test problem and the Taylor–Green vortex-decay problem. It is demonstrated empirically that the Galerkin (Fourier) equations involving Np degrees of freedom, where p is the number of space dimensions, give simulations at least as accurate as finite-difference simulations involving (2N)p degrees of freedom. The theoretical basis for the improved accuracy of the Galerkin (Fourier) method is explained. In particular, the nature of aliasing errors is examined in detail. It is shown that ‘aliasing’ errors need not be errors at all, but that aliasing should be avoided in flow simulations. An eigenvalue analysis of schemes for simulation of passive scalar convection supplies the mathematical basis for the improved accuracy of the Galerkin (Fourier) method. A comparison is made of the computational efficiency of Galerkin and finite-difference simulations, and a survey is given of those problems where Galerkin methods are likely to be applied most usefully. We conclude that numerical simulation of many of the flows of current interest is done most efficiently and accurately using the spectral methods advocated here.

Type
Research Article
Information
Journal of Fluid Mechanics , Volume 49 , Issue 1 , 13 September 1971 , pp. 75 - 112
Copyright
© 1971 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Arakawa, A. 1966 Computational design for long-term numerical integration of the equations of motion : two-dimensional incompressible flow. Part 1. J. Comp. Phys., 1, 119–143.Google Scholar
Araeawa, A. 1970 Numerical simulation of large-scale atmospheric motions. Numerical Solution of Field Problems in Continuum Physics Proc. SIAM-AMS, vol. 2, 24–40. Providence : American Mathematical Society.
Burstein, S. Z. & Mirin, A. A. 1970 Third-order difference methods for hyperbolic equations. J. Comp. Phys., 5, 547–571.Google Scholar
Crank, J. & Nicolson, P. 1947 A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type. Proc. Camb. Phil. Soc., 43, 50–67.Google Scholar
Crowley, W. P. 1968 Numerical advection experiments. Mon. Weath. Rev., 96, 1–11.Google Scholar
Deardorff, J. W. 1970 A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers. J. Fluid Mech., 41, 453–480.Google Scholar
Deardorff, J. W. 1971 On the magnitude of the subgrid scale eddy coefficient. J. Comp. Phys., 7, 120–133.Google Scholar
Ellsaesser, H. W. 1966 Evaluation of spectral versus grid methods of hemispheric numerical weather prediction. J. Appl. Meteor., 5, 246–262.Google Scholar
Fromm, J. E. 1963 A method for computing nonsteady incompressible, viscous fluid flows. Los Alamos Sci. Lab. Rep. LA-2910.Google Scholar
Fromm, J. E. 1969 Practical investigation of convective difference approximations. Phys. Fluids (suppl. 2) 12, 3–12.Google Scholar
Fromm, J. E. & Harlow, F. H. 1963 Numerical solution of the problem of vortex street development. Phys. Fluids, 6, 975–982.Google Scholar
Gaunt, J. A. 1927 The deferred approach to the limit. Part 2. Interpenetrating lattices. Phil. Trans. A 226, 350–361.Google Scholar
Goldstein, S. 1940 Three-dimensional vortex motion in a viscous fluid. Phil. Mag., 30, 85–102.Google Scholar
Grammeltvedt, A. 1969 A survey of finite-difference schemes for the primitive equations for a barotropic fluid. Mon. Weath. Rev., 97, 384–404.Google Scholar
Harlow, F. H. & Welch, J. E. 1965 Numerical calculation of time-dependent viscous incompressible flow of fluid with free surface. Phys. Fluids, 8, 2182–2189.Google Scholar
Lilly, D. K. 1964 Numerical solutions for the shape-preserving two-dimensional thermal convection element. J. Atmos. Sci., 21, 83–98.Google Scholar
Lilly, D.K. 1965 On the computational stability of numerical solutions of time-dependent non-linear geophysical fluid dynamics problems. Mon. Weath. Rev., 93, 11–26.Google Scholar
Lilly, D. K. 1967 The representation of small-scale turbulence in numerical simulation experiments. Proc. IBM Scientific Computing Symposium on Environmental Sciences, 195–210. White Plains, N.Y.: IBM.
Moleneamp, C. R. 1968 Accuracy of fmite-difference methods applied to the advection equation. J. Appl. Meteor., 7, 160–167.Google Scholar
Orszag, S. A. 1969 Numerical methods for the simulation of turbulence. Phys. Fluids (suppl. 2) 12, 250–257.Google Scholar
Orszag, S. A. 1970 Transform method for the calculation of vector-coupled sums: application to the spectral form of the vorticity equation. J. Atmos. Sci., 27, 890–895.Google Scholar
Orszag, S. A. 1971a Numerical simulation of incompressible flows within simple boundaries. 1. Galerkin (spectral) representations. Stud. in Appl. Math. To be published.Google Scholar
Orszag, S. A. 1971b Galerkin approximations to flows within slabs, spheres, and cylinders. Phys. Rev. Letters, 26, 1100–1133.Google Scholar
Orszag, S. A. 1971c On the resolution requirements of kite-difference schemes. Stud. in Appl. Math. To be published.Google Scholar
Patterson, G. S. & Orszag, S. A. 1971 Spectral calculations of isotropic turbulence: efficient removal of aliasing interactions. Phys. Fluids. To be published.Google Scholar
Phillips, N. A. 1959 An example of non-linear computational instability. The Atmosphere and the Sea in Motion, 501–504. New York: Rockefeller Institute Press.
Piacsek, S. A. & Williams, G. P. 1970 Conservation properties of convection difference schemes. J. Comp. Phys., 6, 392–405.Google Scholar
Platzman, G. W. 1961 An approximation to the product of discrete functions. J. Meteor., 18, 31–37.Google Scholar
Price, H. S. & Varga, R. S. 1970 Error bounds for semidiscrete Galerkin approximations of parabolic problems with applications to petroleum reservoir mechanics. Numerical Solution of Field Problems in Continuum Physics, Proc. SSIAM-AMS vol. 2, 74–94. Providence : American Mathematical Society.
Richardson, L. F. 1910 The approximate arithmetical solution by finite differences of physical problems involving differential equations, with an application to the stresses in a masonry dam. Phil. Trans. A 210, 307–357. (Also Proc. Roy. Soc. A 83, 335–336.)Google Scholar
Richardson, L. F. 1927 The deferred approach to the limit. Part 1. Single lattice. Phil. Trans. A 226, 299–349.Google Scholar
Roberts, K. V. & Weiss, N. O. 1966 Convective difference schemes. Math. Comput., 20, 272–349.Google Scholar
Shuman, F. G. & Stackpole, J. D. 1969 The currently operational NMC model, and results of a recent simple numerical experiment. Proc. WMO/IUGG Symp. on Numerical Weather Prediction, vol. 2, 85–98. Tokyo : Japan Meteorological Agency.
Smauorinsky, J. 1963 General circulation experiments with the primitive equations. Part 1. The basic experiment. Mon. Weath. Rev., 91, 99–164.Google Scholar
Taylor, G. I. & Green, A. E. 1937 Mechanism of the production of small eddies from large ones. Proc. Roy. Soc. A 158, 499–521.Google Scholar
Thompson, P. D. 1961 Numerical Weather Analysis and Prediction. Macmillan.
Wilkinson, J. H. 1965 The Algebraic Eigenvalue Problem. Oxford University Press.
Williams, G. P. 1969 Numerical integration of the three-dimensional Navier-Stokes equations for incompressible flow. J. Fluid Mech., 37, 727–750.Google Scholar