Abstract

We consider the finite element solution of the stream function-vorticity formulation for a large-scale ocean circulation model. First, we study existence and uniqueness of solution for the continuous and discrete problems. Under appropriate regularity assumptions we prove that the stream function can be computed with an error of order h in H1-seminorm. Second, we introduce and analyze an h-adaptive mesh refinement strategy to reduce the spurious oscillations and poor resolution which arise when convective terms are dominant. We propose an a posteriori anisotropic error indicator based on the recovery of the Hessian from the finite element solution, which allows us to obtain well adapted meshes. The numerical experiments show an optimal order of convergence of the adaptive scheme. Furthermore, this strategy is efficient to eliminate the oscillations around the boundary layer. © 2003 Elsevier B.V. All rights reserved.