We discuss the performance of iterative solution methods implemented in a 3D model of the shallow water equations used to simulate the propagation of tidal waves in stratified ocean domains. The model uses a harmonic approach by means of which the transient problem is transformed into a number of steady-state nonlinear problems (one per each tidal frequency considered). For the vertical interpolation of the model variables, a spectral decomposition in terms of an orthogonal set of basis functions is applied. A Picard iteration is used to solve the nonlinear systems of equations, and a preconditioned conjugate gradient method is used to solve the intermediate linear systems. We report 3D numerical results of the tidal circulation in the Strait of Gibraltar, and we compare the performance of different conjugate gradient methods implemented in the model with the performance obtained with classical direct methods.
|Publication status||Published - 8 Jan 1998|
|Event||10th International Conference on Finite Elements in Fluids - Tucson, United States|
Duration: 5 Jan 1998 → 8 Jan 1998
|Conference||10th International Conference on Finite Elements in Fluids|
|Period||5/01/98 → 8/01/98|
- OCEAN DYNAMICS
- ITERATIVE SOLUTION
- ; THREE DIMENSIONAL MODELS; WAVE PROPAGATION; TIDES; CONJUGATE GRADIENT METHOD; NONLINEAR EQUATIONS
Gonzalez, M., Espino, M., Chumbe, S., Garcia, M., & Sanchez-Arcilla, A. (1998). Performance of iterative algorithms for the solution of the shallow-water equations. 381-386. Paper presented at 10th International Conference on Finite Elements in Fluids, Tucson, United States.