Resolution Effects and Enslaved Finite-Difference Schemes for a Double Gyre, Shallow-Water Model

1997 ◽  
Vol 9 (3-4) ◽  
pp. 269-280 ◽  
Author(s):  
Don A. Jones ◽  
Andrew C. Poje ◽  
Len G. Margolin
2019 ◽  
Vol 34 (2) ◽  
pp. 105-117 ◽  
Author(s):  
Gayaz S. Khakimzyanov ◽  
Zinaida I. Fedotova ◽  
Oleg I. Gusev ◽  
Nina Yu. Shokina

Abstract Basic properties of some finite difference schemes for two-dimensional nonlinear dispersive equations for hydrodynamics of surface waves are considered. It is shown that stability conditions for difference schemes of shallow water equations are qualitatively different in the cases the dispersion is taken into account, or not. The difference in the behavior of phase errors in one- and two-dimensional cases is pointed out. Special attention is paid to the numerical algorithm based on the splitting of the original system of equations into a nonlinear hyperbolic system and a scalar linear equation of elliptic type.


2006 ◽  
Vol 36 (7) ◽  
pp. 1265-1286 ◽  
Author(s):  
Tomonori Matsuura ◽  
Mitsutaka Fujita

Abstract A two-layer shallow-water model is used to investigate the transition of wind-driven double-gyre circulation from laminar flow to turbulence as the Reynolds number (Re) is systematically increased. Two distinctly different phases of turbulent double-gyre patterns and energy trajectories are exhibited before and after at Re = 95: deterministic and fully developed turbulent circulations. In the former phase, the inertial subgyres vary between an asymmetric solution and an antisymmetric solution and the double-gyre circulations reach the aperiodic solution mainly due to their barotropic instability. An integrated kinetic energy in the lower layer is slight and the generated mesoscale eddies are confined in the upper layer. The power spectrum of energies integrated over the whole domain at Re = 70 has peaks at the interannual periods (4–7 yr) and the interdecadal period (10–20 yr). The loops of the attractors take on one cycle at those periods and display the blue-sky catastrophe. At Re = 95, the double-gyre circulation reaches a metastable state and the attracters obtained from the three energies form a topological manifold. In the latter, as Re increases, the double-gyre varies from a metastable state to a chaotic state because of the barotropic instability of the eastward jet and the baroclinic instability of recirculation retrograde flow, and the eastward jet meanders significantly with interdecadal variability. The generated eddies cascade to the red side of the power spectrum as expected in the geostrophic turbulence. The main results in the simulation may indicate essential mechanisms for the appearance of multiple states of the Kuroshio and for low-frequency variations in the midlatitude ocean.


Author(s):  
Zinaida I. Fedotova ◽  
Gayaz S. Khakimzyanov

AbstractThe paper contains a description of the most important properties of numerical methods for solving nonlinear dispersive hydrodynamic equations and their distinctions from similar properties of finite difference schemes approximating classic dispersion-free shallow water equations.


2013 ◽  
Vol 3 (1) ◽  
pp. 27-47 ◽  
Author(s):  
Farzana Hussain

AbstractThe Bay of Bengal is surrounded by coastline except to the south, where there is open sea. The coastline bends most sharply along the coast of Bangladesh, and there are many small and large islands in the off shore region there. In order to incorporate the island boundaries and the curved coastline properly, in any numerical scheme it is often necessary to consider a very fine grid resolution along the coastal belts whereas this is unnecessary away from the coasts. However, a very fine resolution involves more memory and more CPU time in the numerical solution process, and invites numerical instability. On the other hand, boundary-fitted curvilinear grids in hydrodynamic models for coastal seas, bays and estuaries not only fit to the coastline but also render the finite difference schemes simpler and more accurate. In this article, the boundary-fitted curvilinear grids for the model represent the complete boundary of the area considered by four curves defined by four functions, and the four boundaries of two of the larger islands are then represented approximately by two general functions. An appropriate independent coordinate transformation maps the curvilinear physical area to a square domain, and each island boundary is transformed to a rectangle within this square domain. The vertically integrated shallow water equations are transformed to the new space domain, and solved by a regular explicit finite difference scheme. The model is applied to compute the water levels due to astronomical tides, and also the water levels due to surges associated with tropical storms that hit the coast of Bangladesh.


Sign in / Sign up

Export Citation Format

Share Document