& Tornado Alley
The equations for this case are:
where u' and h' are the deviations of fluid velocity and height from the mean background values of u and H, respectively. The parameter g is the acceleration due to gravity.
Let U respresent the column vector of (u', h')T. The above equations in matrix form are then
where A is the matrix
When solutions of the form U(x,t) = U0eik(x-ct) are sought it is found that c has to be such that
Thus c must be an eigenvalue of A. The determinant equation reduces to
Consider now the simplified case in which u=0. In this case the shallow water wave equations reduce to:
Take c to be (gH)1/2 and consider the transformation of the dependent variables
When we compute ∂w/∂t - c(∂w/∂x) we find it is
The term in the brackets is zero because of the governing equation for u'.
The coefficient of the last term is (-c), which can be represented as -(g/c)(c2/g). Therefore the common factor of -(g/c) can be factored from the last terms leaving
which is also zero because of the second governing equation.
A similar computation shows that
Consider now the transformation of independent variables applied to the equation for w:
The differential operators for the old variables are defined in terms of the new variables as:
When these operators are applied in the eqution for q the result is:
This result, ∂q/∂s = 0, is very significant. It means that the variable q, which was created as u'+(g/c)h', is a function of z=(x-ct) alone.
A similar change of variables
establishes that w=u'-(g/c)h' is a function of z=x+ct alone. Therefore (q+w)/2, which is equal to u', is a function of (x-ct) and (x+ct). Likewise h', which is the same as (q-w)/2, is a function of (x-ct) and (x+ct). These are wave functions with one wave traveling to the right with a velocity of c and the other traveling to the left with a velocity of -c.
The solutions for u' and h' are then of the forms
Now the question is how to find the functions U() and G(). The answer is easier than might be supposed. The unknown functions must satisfy the initial conditions, the profile of u' and h' as functions of x at time equal to zero. Let u0(x) and h0(x) be the spatial profiles of u' and h' at t=0. Then putting t=0 in the solutions gives
Thus the solutions for this case of the one dimensional shallow water wave equations are
The solutions to the one dimensional shallow water wave equation is then a wave spreading right and left and maintaining the profile that existed at time zero.
An interesting question is what profile the finite difference scheme approximates in that the full profile is not available to the discrete model.
HOME PAGE of Thayer Watkins
HOME PAGE OF applet-magic.com