A00, B00, C00
A00, B00, C00 are the weighting factors (at time levels k+1, k, k-1, respectively) for the free surface and boundary fluxes in the GWCE, and must sum to 1. Most critically, the weighting factors are used in the discretization of the linear gravity wave (pressure gradient) term and are responsible for determining the inherent implicity (impacting solution stability), in addition to order of accuracy and dispersive characteristics of the numerical method.
If the consistent mass-matrix solver is chosen (see IM parameter) then a semi-implicit method is possible and encouraged. In this case the most common choice for the weighting factors are:
If the lumped mass-matrix solver is chosen then only an explicit method is possible (the weighting A00 must be zero as no matrix solve is conducted), and the weighting factors that are typically chosen become simply:
Wave Continuity Equation
Theory is dominated by analysis of the Wave Continuity Equation (WCE), a special case of the Generalized Wave Continuity Equation (GWCE) where the TAU0 parameter is equal to the linear friction coefficient. In what has been determined to be a third-order accurate method centered in time, which was first introduced by Lynch and Gray (1979), the choice of A00, B00, C00 is reduced to depend on a single parameter, :
In other words, k+1 and k-1 weightings are chosen to be equal. It would however not appear that any restriction other than the requirement that A00, B00, C00 must sum to 1 is necessary to obtain second-order accuracy. Unconditional stability is achieved with the prescription of . Likely because of this fact, the typical choice for ADCIRC has become , i.e., as noted above. Different values of may be motivated by the following expression for optimal dispersive accuracy for the consistent mass-matrix solver:
where is the Courant number based on the linear gravity wave speed.
- : lumped mass-matrix solved in 1-D
- : consistent mass-matrix solved in 1-D
- : lumped mass-matrix solved in 2-D
- : consistent mass-matrix solved in 2-D
These conditions are for linear finite-elements (ADCIRC uses these) with even node spacings and constant bathymetry. Other conditions for quadratic finite-elements, uneven node spacings, and non-constant bathymetry are shown in Kinnmark and Gray (1985).
In the explicit method case, stability (and optimal dispersive accuracy) is superior for the lumped mass-matrix solver versus the consistent mass-matrix solver, hence the lumped solver should always be chosen when employing an explicit method (see IM parameter for setting the solver type).
Generalized Wave Continuity Equation
For the theory based on the WCE, stability is shown to be independent of the choice of TAU0 (). However, experience tell us that in the GWCE a larger value of tends to be more unstable than a smaller value. This makes sense since the behavior of the equations will become more and more similar to the Primitive Continuity Equation with greater , which is responsible for 2Δx instabilities - the motive for using the GWCE in the finite-element method.
Kinnmark's 1986 monograph does offer some analysis that backs up this empirical experience. For the consistent mass-matrix solver employing the centered scheme it can be shown that unconditional stability is only possible if,
where is the linear friction coefficient. In the deep ocean when employing a quadratic drag law this could be a very restrictive requirement on (e.g., assuming , m/s, m). In addition, it is in opposition to our desire to choose the recommended value, , for good mass conservation and dispersive properties. An alternative way to view the criteria is to look for a Courant number which eliminates the criteria on the choice of . It can be shown that this is the case in 1-D when,
Assuming the same relationship between 1-D and 2-D as in the WCE, then the model should be stable for any in 2-D when . This means the stability requirement is slightly improved from the lumped explicit method of the WCE (unless one can choose a sufficiently small to obtain greater stability).
One aspect missing from the theory above are different weighting schemes which are possibly non-centered, i.e., . Empirical evidence suggests stability of the GWCE could be improved if,
Moreover, the equations analyzed are always linearized (a requirement of the von Neumann stability analysis), thus stability may be more restricting in real-world problems, especially those with fine grid sizes and where nonlinearities are non-trivial. In such cases where it is not possible to achieve time steps more than twice that possible with an explicit method it becomes preferable to employ the explicit lumped mass-matrix solver since it is computationally twice as fast per time step solve.
- M.G.G. Foreman, An Analysis of the “Wave Equation” Model for Finite Element Tidal Computations, J. Comput. Phys. 52 (1983) 290–312.
- D.R. Lynch, W.G. Gray, A Wave Equation Model for Finite Element Tidal Computations, Computers & Fluids. 7 (1979) 207–228.
- I.P.E. Kinnmark, W.G. Gray, Stability and accuracy of spatial approximations for wave equation tidal models, J. Comput. Phys. 60 (1985) 447–466. doi:10.1016/0021-9991(85)90030-0.
- Kinnmark, I.P.E., Gray, W.G., 1984. A Two-Dimensional Analysis of the Wave Equation Model for Finite Element Tidal Computations. Int. J. Numer. Methods Eng. 20, 369–383.
- Kinnmark, I., 1986. The Shallow Water Wave Equations: Formulation, Analysis and Application, Lecture Notes in Engineering. Springer Berlin Heidelberg, Berlin, Heidelberg. doi:10.1007/978-3-642-82646-7.
- Kolar, R.L., Westerink, J.J., Cantekin, M.E., Blain, C.A., 1994. Aspects of Nonlinear Simulations using Shallow-water Models based on the Wave Continuity Equation. Comput. Fluids 23, 523–538.