Premium
Bivariate spline method for numerical solution of steady state Navier–Stokes equations over polygons in stream function formulation
Author(s) -
Lai MingJun,
Wenston Paul
Publication year - 2000
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/(sici)1098-2426(200003)16:2<147::aid-num2>3.0.co;2-9
Subject(s) - mathematics , spline (mechanical) , nonlinear system , uniqueness , galerkin method , finite element method , mathematical analysis , partial differential equation , thin plate spline , function space , navier–stokes equations , spline interpolation , statistics , physics , structural engineering , quantum mechanics , aerospace engineering , compressibility , engineering , bilinear interpolation , thermodynamics
We use the bivariate spline finite elements to numerically solve the steady state Navier–Stokes equations. The bivariate spline finite element space we use in this article is the space of splines of smoothness r and degree 3 r over triangulated quadrangulations. The stream function formulation for the steady state Navier–Stokes equations is employed. Galerkin's method is applied to the resulting nonlinear fourth‐order equation, 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 H 2 (Ω) of the nonlinear fourth‐order problem and give an estimate of how fast the numerical solution converges to the weak solution. The Galerkin method with C 1 cubic splines is implemented in MATLAB. Our numerical experiments show that the method is effective and efficient. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 147–183, 2000