GenCade:Numerical Stability: Difference between revisions

From CIRPwiki
Jump to navigation Jump to search
mNo edit summary
mNo edit summary
Line 14: Line 14:


where <math> \varepsilon_1 = \frac{2Q_O a_1}{D_B + D_C} </math> and <math> \varepsilon_2 = \frac{Q_O a_2 sin\alpha_b}{D_B + D_C} {\frac{\partial H_b}{\partial x}} </math>
where <math> \varepsilon_1 = \frac{2Q_O a_1}{D_B + D_C} </math> and <math> \varepsilon_2 = \frac{Q_O a_2 sin\alpha_b}{D_B + D_C} {\frac{\partial H_b}{\partial x}} </math>
In the presence of an external current not generated by breaking waves, the second diffusion coefficient in the above partial differential will change to:
<math> \varepsilon_2 = \frac{Q_O sin\alpha_b}{D_B + D_C} (a_2 \frac{\partial H_b}{\partial x} - a_3 \frac{\bar{v}_t}{u_m}) </math>
As the above partial differential is a diffusion-type equation, its stability properties are well known. The numerical stability of the calculation scheme is governed by:
<math> R_s = \frac{\Delta t (\varepsilon_1 + \varepsilon_2)}{(\Delta x)^2} </math>
where the quantity <math> R_s </math> is known as the Courant number in numerical methods; here it is called the stability parameter. The finite difference form of the above governing partial differential equation shows that <math> \Delta y \approx \frac{\Delta t}{(\Delta x)^2} </math> which means that if <math> \Delta x </math> is reduced by a factor of two, the time step will need to be reduced by a factor of four to maintain the same stability of the calculation scheme. If an explicit solution scheme is used to solve the diffusion equation, the following condition must be satisfied (Crank 1975):
<math> R_s \le 0.5 </math>
If the value of <math> R_s </math> exceeds 0.5 at any point on the grid, the calculated shoreline will show an unphysical oscillation that will grow in time if <math> R_s </math> remains above 0.5, alternating in direction at each grid point. The quantities <math> \varepsilon_1 </math> and <math> \varepsilon_2 </math> can change greatly alongshore since they depend on the local wave conditions. Assuming that the grid cell spacing is fixed by engineering requirements, a large wave height would necessitate a small value of <math> \Delta t </math>. The GenCade model will issue a warning if the stability requirement is violated at any point in the scheme. If such a warning is issued, either the time step <math> \Delta t </math> or the spatial resolution <math> \frac{1}{\Delta x} </math> (or both) will need to be reduced. Thus, it is necessary to increase <math> \Delta t </math> and/or decrease <math> \Delta x </math>.

Revision as of 21:24, 3 January 2023

This page is currently in editing. Please check in later to see updates.

The GenCade model utilizes an explicit solution scheme. The main advantages of this scheme are: easy programming, simple (and sometimes the only possible) expressions of boundary conditions, and shorter computer run-time as compared to an implicit scheme (for a single time increment). A major disadvantage is, however, the stability of the solution. This means that smaller time steps are often needed, and thus a larger number, when using an explicit scheme as compared to using an implicit scheme. As a consequence, simulations using explicit schemes often require a longer total computation time. To minimize the computational effort, the longest time step that may be used for a specific calculation must be determined. Under certain idealized conditions, the CERC equation can be reduced to a simpler form to examine the dependence of the solution on the time and space steps. First, rewrite the CERC equation in the form:

where (cubic meters/second). A useful approximate stability criterion can be obtained by linearizing the governing equation with respect to y. The linearization is made by assuming small breaking wave and shoreline angles, which leads to:

and

where is the angle of breaking waves to the local shoreline orientation, is the angle of breaking waves to the x-axis, and is the local shoreline orientation. Assuming that is zero, the governing equation becomes (Kraus and Harikai 1983):

where and

In the presence of an external current not generated by breaking waves, the second diffusion coefficient in the above partial differential will change to:

As the above partial differential is a diffusion-type equation, its stability properties are well known. The numerical stability of the calculation scheme is governed by:

where the quantity is known as the Courant number in numerical methods; here it is called the stability parameter. The finite difference form of the above governing partial differential equation shows that which means that if is reduced by a factor of two, the time step will need to be reduced by a factor of four to maintain the same stability of the calculation scheme. If an explicit solution scheme is used to solve the diffusion equation, the following condition must be satisfied (Crank 1975):

If the value of exceeds 0.5 at any point on the grid, the calculated shoreline will show an unphysical oscillation that will grow in time if remains above 0.5, alternating in direction at each grid point. The quantities and can change greatly alongshore since they depend on the local wave conditions. Assuming that the grid cell spacing is fixed by engineering requirements, a large wave height would necessitate a small value of . The GenCade model will issue a warning if the stability requirement is violated at any point in the scheme. If such a warning is issued, either the time step or the spatial resolution (or both) will need to be reduced. Thus, it is necessary to increase and/or decrease .