Bivariate spline method for numerical solution of time evolution Navier-Stokes equations over polygons in stream function formulation

✍ Scribed by Ming-Jun Lai; Chun Liu; Paul Wenston

• Publisher: John Wiley and Sons
• Year: 2003
• Tongue: English
• Weight: 410 KB
• Volume: 19
• Category: Article
• ISSN: 0749-159X
• DOI: 10.1002/num.10072

✦ Synopsis

Abstract

We use a bivariate spline method to solve the time evolution Navier‐Stokes equations numerically. The bivariate splines we use in this article are in the spline space of smoothness r and degree 3__r__ over triangulated quadrangulations. The stream function formulation for the Navier‐Stokes equations is employed. Galerkin's method is applied to discretize the space variables of the nonlinear fourth‐order equation, Crank‐Nicholson's method is applied to discretize the time variable, and Newton's iterative method is then used to solve the resulting nonlinear system. We show the existence and uniqueness of the weak solution in L~2~(0, T; H^2^(Ω)) ∩ L~∞~(0, T; H^1^(Ω)) of the 2D nonlinear fourth‐order problem and give an estimate of how fast the numerical solution converges to the weak solution. The C^1^ cubic splines are implemented in MATLAB for solving the Navier‐Stokes equations numerically. Our numerical experiments show that the method is effective and efficient. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 776–827, 2003.