Dec 2, 1998

PHY 711 - Notes on Hydrodynamics - ("Solitary Waves"[])

Basic assumptions

We assume that we have in incompressible fluid ( r = constant) a velocity potential of the form F(x,z,t). The surface of the fluid is described by h+z(x,t) = z. The fluid is contained in a tank with a structureless bottom (defined by the plane z = 0) and is filled to a vertical height h at equilibrium. These functions satisfy the following conditions.
Poisson equation: 2 F(x,z,t)
x2
 + 2 F(x,z,t)
z2
 = 0
(1)
Zero vertical velocity at bottom of the tank: F(x,0,t)
z
 = 0
(2)
Bernoulli¢s equation: - F(x,z,t)
t
 + 1
2
 æ
ç
è 		F(x,z,t)
x
 ö
÷
ø 		2

 
 +gz(x,t) ú
ú
û

z = h+z 
 = 0
(3)
Surface equation: - F(x,z,t)
z
 + F(x,z,t)
x
 z(x,t)
x
 - z(x,t)
t
 ú
ú
û

z = h+z 
 = 0
(4)

In this treatment, we assume seek the form of surface waves traveling along the x- direction and assume that the effective wavelength is much larger than the height of the surface h. This allows us to approximate the z- dependence of F(x,z,t) by means of a Taylor series expansion:

F(x,z,t) » F(x,0,t) + z F
z
 (x,0,t)+ z2
2
 2 F
z2
 (x,0,t) + z3
3!
 3 F
z3
 (x,0,t) ¼
(5)

and

F(x,z,t)
z
 » + F
z
 (x,0,t)+ z 2 F
z2
 (x,0,t) + z2
2
 3 F
z3
 (x,0,t) + z3
3!
 4 F
z4
 (x,0,t) ¼.
(6)

These expansions can be simplified because of the bottom boundary condition (2) and the Poisson equation (1). In particular, (2) and (1) ensure that all odd derivatives [(n F)/( zn)](x,0,t) vanish from the Taylor expansions (5) and (6). Taking the time derivative of 5 and evaluating it at the surface, z = h+z(x,t), and using the Poisson equation (1) to convert all derivatives with respect to z to derivatives with respect to x, we find:

F(x,z,t)
t
 » F(x,0,t)
t
 -
t
 æ
ç
è 		(h+z(x,t))2
2
 2 F
x2
 (x,0,t) ö
÷
ø 		¼,
(7)

and simplifying (6), we find:

F(x,z,t)
z
 » - (h+z(x,t)) 2 F
x2
 (x,0,t) + (h+z(x,t))3
3!
 4 F
x4
 (x,0,t) ¼.
(8)

We can now use these results to approximate the boundary conditions implied by Bernoulli's equation (3 and the time derivative the of the surface equation (4). Keeping all terms up to leading order in non-linearity and up to fourth order derivatives in the linear terms we find that the surface definition equation becomes:

x
 æ
ç
è 		(h+z(x,t)) F(x,0,t)
x
 ö
÷
ø 		- h3
3!
 4 F(x,0,t)
x4
 - z(x,t)
t
 = 0,
(9)
and the Bernoulli's equation becomes:
- F(x,0,t)
t
 + h2
2
 3 F(x,0,t)
t x2
 + 1
2
 æ
ç
è 		F(x,0,t)
t
 ö
÷
ø 		2

 
 +gz(x,t) = 0.
(10)

We would like to solve Eqs. (9-10) for a traveling wave of the form:

F(x,0,t) = f(x-ct) and z(x,t) = y(x-ct),
(11)
where the speed of the wave c will be determined. Letting u º x-ct, Eqs. (9 and 10) become:
u
 æ
ç
è 		(h + y(u)) f(u)
u
 ö
÷
ø 		- h3
6
 4 f(u)
u4
 +c y(u)
u
 = 0,
(12)
and
c f(u)
u
 - ch2
2
 3 f(u)
u3
 + c2
2
 æ
ç
è 		f(u)
u
 ö
÷
ø 		2

 
 +g y(u) = 0.
(13)

The modified surface equation (12) can be integrated once with respect to u, choosing the constant of integration to be zero and giving the new form for the surface condition:

(h+y(u)) f(u)
u
 - h3
6
 3 f(u)
u3
 +cy(u) = 0.
(14)

Using the modified Bernoulli equation (13) to eliminate the function f(u) from Eq. (14), discarding some small terms as described in the text, we determine the following non-linear equation for y(u):

æ
ç
è 		1 - hg
c2
 ö
÷
ø 		y(u) - 3
2h
 ( y(u) )2 - h2
3
 2 y(u)
u2
 = 0.
(15)

Your text shows that a solution to Eq. (15) is the solitary wave form:

z(x,t) = y(x-ct) = y0
cosh2 é
ê
ë 		æ
ç
è 		  æ
Ö
[(3y0)/ h]
 
 x-ct
2h
 ö
÷
ø 		ù
ú
û
,
(16)
with
c2 = gh
1-y0/h
 .
(17)

References

[]
Alexander L. Fetter and John Dirk Walecka, Theoretical Mechanics of Particles and Continua, (McGraw Hill, 1980), Chapt. 10.

File translated from TEX by TTH, version 1.92.
On 2 Dec 1998, 22:28.