This article is about the law in fluid dynamics. It is not to be confused with
Green's Theorem
.
Equation describing evolution of waves in shallow water
In
depth contours
parallel to each other (and the coast), it states:
or
where and are the wave heights at two different locations – 1 and 2 respectively – where the wave passes, and and are the mean water depths at the same two locations.
Green's law is often used in
nonlinear effects become important and Green's law no longer applies.[2][3]
in water of mean depth and width (in case of an open channel) satisfy[4][5]
where is the
fourth root
of Consequently, when considering two cross sections of an open channel, labeled 1 and 2, the wave height in section 2 is:
with the subscripts 1 and 2 denoting quantities in the associated cross section. So, when the depth has decreased by a factor sixteen, the waves become twice as high. And the wave height doubles after the channel width has gradually been reduced by a factor four. For wave propagation perpendicular towards a straight coast with depth contours parallel to the coastline, take a constant, say 1 metre or yard.
For refracting long waves in the ocean or near the coast, the width can be interpreted as the distance between wave rays. The rays (and the changes in spacing between them) follow from the geometrical optics approximation to the linear wave propagation.[6] In case of straight parallel depth contours this simplifies to the use of Snell's law.[7]
Green published his results in 1838,
Liouville–Green method – which would evolve into what is now known as the WKB approximation. Green's law also corresponds to constancy of the mean horizontal wave energy flux for long waves:[4][5]
) of shoaling waves does not change, according to Green's linear theory.
Derivation
Green derived his shoaling law for water waves by use of what is now known as the Liouville–Green method, applicable to gradual variations in depth and width along the path of wave propagation.[9]
Wave equation for an open channel
Starting point are the linearized
one-dimensional Saint-Venant equations for an open channel with a rectangular cross section (vertical side walls). These equations describe the evolution of a wave with free surface
elevation and horizontal flow velocity with the horizontal coordinate along the channel axis and the time:
where is the gravity of Earth (taken as a constant), is the mean water depth, is the channel width and and are denoting partial derivatives with respect to space and time. The slow variation of width and depth with distance along the channel axis is brought into account by denoting them as and where is a small parameter: The above two equations can be combined into one wave equation for the surface elevation:
and with the velocity following from
(1)
In the Liouville–Green method, the approach is to convert the above wave equation with non-homogeneous coefficients into a homogeneous one (neglecting some small remainders in terms of ).
Transformation to the wave phase as independent variable
The next step is to apply a
wave phase
) given by
so
and are related through the
celerity
Introducing the slow variable and denoting derivatives of and with respect to with a prime, e.g. the -derivatives in the wave equation, Eq. (1), become:
The next step is transform the equation in such a way that only deviations from homogeneity in the second order of approximation remain, i.e. proportional to
Further transformation towards homogeneity
The homogeneous wave equation (i.e. Eq. (2) when is zero) has solutions for
travelling waves
of permanent form propagating in either the negative or positive -direction. For the inhomogeneous case, considering waves propagating in the positive -direction, Green proposes an approximate solution:
So the proposed solution in Eq. (3) satisfies Eq. (2), and thus also Eq. (1) apart from the above two terms proportional to and , with The error in the solution can be made of order provided
This has the solution:
Using Eq. (3) and the transformation from to , the approximate solution for the surface elevation is
(4)
where the constant has been set to one, without loss of generality. Waves travelling in the negative -direction have the minus sign in the argument of function reversed to a plus sign. Since the theory is linear, solutions can be added because of the superposition principle.
The horizontal flow velocity in the -direction follows directly from substituting the solution for the surface elevation from Eq. (4) into the expression for in Eq. (1):[10]
Satake, K. (2002), "28 – Tsunamis", in Lee, W. H. K.; Kanamori, H.; Jennings, P. C.; Kisslinger, C. (eds.), International Handbook of Earthquake and Engineering Seismology, International Geophysics, vol. 81, Part A,
Synolakis, C. E.; Skjelbreia, J. E. (1993), "Evolution of maximum amplitude of solitary waves on plane beaches", Journal of Waterway, Port, Coastal, and Ocean Engineering, 119 (3): 323–342,