NOAA

Geophysical Fluid
Dynamics Laboratory

Skip to: [content] [navigation]
If you are using Navigator 4.x or Internet Explorer 4.x or Omni Web 4.x , this site will not render correctly!

gfdl's home page > gfdl's bibliography > gareth p. williams > Journal of Fluid Mechanics, 37 (4), 727-750.

Numerical integration of the three-dimensional Navier-Stokes equations for incompressible flow

Williams, G. P., 1969: Numerical integration of the three-dimensional Navier-Stokes equations for incompressible flow. Journal of Fluid Mechanics, 37 (4), 727-750.
Abstract: A method of numerically integrating the Navier-Stokes equations for certain three-dimensional incompressible flows is described. The technique is presented through application to the particular problem of describing thermal convection in a rotating annulus. The equations, in cylindrical polar co-ordinate form, are integrated with respect to time by a marching process, together with the solving of a Poisson equation for the pressure. A suitable form of the finite difference equations gives a computationally-stable long-term integration with reasonably faithful representation of the spatial and temporal characteristics of the flow.
Trigonometric interpolation techniques provide accurate (discretely exact) solutions to the Poisson equation. By using an auxiliary algorithm for rapid evaluation of trigonometric transforms, the proportion of computation needed to solve the Poisson equation can be reduced to less than 25% of the total time needed to advance one time step. Computing on a UNIVAC 1108 machine, the flow can be advanced one time-step in 2 sec for a 14 x 14 x 14 grid upward to 96 sec for a 60 x 34 x 34 grid.
As an example of the method, some features of a solution for steady wave flow in annulus convection are presented. The resemblance of this flow to the classical Eady wave is noted.
smaller bigger reset
last modified: March 24 2004.
this page visited: 1410 times