This paper presents a nonconforming finite element method for a streamfunction formulation of the stationary quasi-geostrophic equations, which describe the large scale wind-driven ocean circulation. The streamfunction formulation is a fourth order nonlinear PDE and the nonconforming method is based on C0-elements instead of C1-elements. Existence and uniqueness of the approximation are proved and optimal error estimates in several norms of interest are demonstrated under a small data assumption. Two-grid algorithms based on Picard and Newton type linearizations are then presented to efficiently resolve nonlinearities and computational results are given to demonstrate the efficiency of the algorithm. The Mediterranean sea example is tested with real world coastline data, which illustrates the effectiveness of the two-grid approach in the wind-driven ocean circulation simulation.
|Number of pages||19|
|Journal||Computer Methods in Applied Mechanics and Engineering|
|Publication status||Published - 2018 Jan 1|
Bibliographical noteFunding Information:
Third author was supported by NRF - 2015R1A5A1009350 and NRF- 2016R1A2B4014358 .
All Science Journal Classification (ASJC) codes
- Computational Mechanics
- Mechanics of Materials
- Mechanical Engineering
- Physics and Astronomy(all)
- Computer Science Applications