This work concerns the development of a finite-element algorithm for the stationary quasi-geostrophic equations to treat the large scale wind-driven ocean circulation. The algorithm is developed based on the streamfunction formulation involving fourth-order gradients of the streamfunction. Here, we examine the adaptation of a relatively inexpensive, nonconforming method based on C0 basis functions. We develop the variational form of the method and establish consistency. The method weakly enforces continuity of the normal flux across interelement boundaries and stabilization is achieved via Nitsche's method. Explicit expression for the choice of the stabilization parameter is derived. Moreover, the theoretical error estimate for the linear Stommel-Munk model is performed. Numerical results from several benchmark problems on rectangular and curved domains are provided to demonstrate the accuracy and robustness of the method. We also provide the Mediterranean sea example to demonstrate the capability of the approach to the wind-driven ocean circulation simulation.
|Number of pages||20|
|Journal||Computer Methods in Applied Mechanics and Engineering|
|Publication status||Published - 2016 Mar 1|
Bibliographical noteFunding Information:
The research of EJP was supported in part by the National Research Foundation of Korea ( NRF2015 R1A5A1009350 ).
All Science Journal Classification (ASJC) codes
- Computational Mechanics
- Mechanics of Materials
- Mechanical Engineering
- Physics and Astronomy(all)
- Computer Science Applications