3.6 Exercise 2: Wave Interference 25
3.5.10 Superposition of Waves
The superposition of two or more waves of different period and/or wavelength can
lead to various interference patterns such as a standing wave, being a wave of vir-
tually zero phase speed. Interfering wave patterns travel with a certain speed, called
group speed that can be different to the phase speeds of the contributing individual
waves. Interference of storm-generated waves in the ocean can result in waves of
gigantic wave heights (wave height is twice the wave amplitude) of >20 m, known
as freak waves.
3.6 Exercise 2: Wave Interference
3.6.1 Aim
The aim of this exercise is to explore interferences that result from the superpo-
sition of two linear waves. To this end, SciLab will be used to calculate possible
interferences pattern and to produce animations thereof.
3.6.2 Task Description
Consider the interference of two waves of the same amplitude of A
o
=1m.The
resultant wave can be described by:
A(x, t) = A
o

sin

2π

x
λ
1
−
t

instance, leads to a standing wave, being the result of two identical waves travelling
in opposite directions. This is achieved by prescribing a negative value of the wave
period for wave 2. Is this surprising how many different interference patterns can be
created by superposition of just two waves. The reader is encouraged to experiment
with other scenarios!
26 3 Basics of Geophysical Fluid Dynamics
3.6.3 Sample Script
The SciLab script for this exercise, called “interference.sce” can be found in the
folder “Exercise 2” on the CD-ROM.
3.6.4 A Glimpse of Results
Figure 3.4 shows a snapshot result for the last scenario. It took me a while to achieve
what I wanted in terms of graphical display and, perhaps, the reader will come up
with a better solution. It can be seen that, at some locations, wave disturbances add
up to create a signal of doubled amplitude, whereas, in other regions, the waves
compensate each other such that the resultant wave signal almost vanishes.
Fig. 3.4 Snapshot results of Scenario 6
3.6.5 A Rule of Thumb
The wave period should be resolved by at least 10 time steps. Otherwise details of
the wave evolution can get lost. In the following SciLab script, I used 20 time steps
in a period. Adequate choices of time steps and grid spacings play an important role
in the modelling of dynamical processes in fluids
3.7 Forces 27
3.7 Forces
3.7.1 What Forces Do
A non-zero force operates to change the speed and/or the direction of motion of a
flui parcel. In geophysical flui dynamics, forces are conventionally expressed as
forces per unit mass, being directly proportional to acceleration or deceleration of
a flui parcel. So, whenever the term “force” appears in the following, this actually
should mean “force per unit mass”.
In component form, a temporal change of the velocity vector can be formally

+ F
x
2
dv
dt
= F
y
1
+ F
y
2
dw
dt
= F
z
1
+ F
z
2
Acceleration or deceleration results if any component of the resultant net force is
different from zero. In the general case, the sum symbol “

” can be used to write:
du
dt
=
m

i=1
F

motion.
3.7.3 Apparent Forces
Apparent forces come into play in a rotating coordinate system such as our Earth.
One of these apparent forces is the Coriolis force that gives rise to circular weather
patterns in the atmosphere, eddies in the ocean, or Jupiter’s Red Spot.
3.7.4 Lagrangian Trajectories
Imagine that you sit on a flui parcel of a certain temperature moving with the f ow.
The path along which you move is called a Lagrangian trajectory, based on work by
Lagrange (1788). Without any heat exchange with the ambient fluid the temperature
of your flui parcel remains constant and this feature can be formulated as:
dT
dt
= 0 (3.6)
where “T” is temperature and the “d” symbol now refers to a change of temperature
along the pathway of motion.
3.7.5 Eulerian Frame of Reference and Advection
Instead of moving with the f ow, you could stand still at a f xed location and measure
changes in temperatures as the flui moves past. This perspective is called the Eule-
rian system, based on work by Euler (1736). In this case, you would notice a change
in temperature if a fl w exists that carries differences (gradients) in temperature
towards you. This process is called advection. In Cartesian coordinates, the effect
of temperature advection can be expressed as:
∂T
∂t
=−u
∂T
∂x
− v
∂T
∂y

∂t
+ u
∂u
∂x
+ v
∂u
∂y
+ w
∂u
∂z
=
m

i=1
F
x
i
∂v
∂t
+ u
∂v
∂x
+ v
∂v
∂y
+ w
∂v
∂z
=
m

with a spoon much more efficien than just waiting until the soup has mixed itself.
The reader can also blame these terms for the unreliability of weather forecasts for
longer than 5 days ahead.
3.8 Fundamental Conservation Principles
3.8.1 A List of Principles
There are several conservation principles that need to be considered when studying
flui motions. These are:
1. Conservation of momentum (Newton’s laws of motion)
2. Conservation of mass
3. Conservation of interal energy (heat)
4. Conservation of salt.
In addition to this comes the so-called equation of state that links the fiel vari-
ables such as temperature and salinity to the density of the fluid All these equations
are coupled with each other, which makes the equations describing flui motions a
coupled system of partial differential equations.
3.8.2 Conservation of Momentum
Conservation of momentum is an expression of changes in fl w speed and/or direc-
tion as a result of forces. The frictional stress imposed by winds along the sea surface
acts as a boundary source term in the momentum equations. Friction at the sea f ow
acts as a sink term in these equations. Forces of relevance to fluid are explored in
the next sections.
3.8.3 Conservation of Volume – The Continuity Equation
Water is largely incompressible, so that the mass of a given water volume cannot
change much under compression. Conservation of mass can therefore be expressed
in terms of conservation of volume. To understand this important concept, consider
a virtual volume element (Fig. 3.5). For simplicity, we orientate this element in such
a way that its face normals are parallel to the directions of the Cartesian coordinate
system. The side-lengths of this box are δx, δy and δz, and the volume is δV =
δx · δy · δz.
3.8 Fundamental Conservation Principles 31

δu
δx
+
δv
δy
+
δw
δz
= 0
The equation is valid for any finit volume and, accordingly, for a vanishingly
small volume, which can be expressed by the partial differential equation
∂u
∂x
+
∂v
∂y
+
∂w
∂z
= 0 (3.10)
being called the continuity equation. This equation constitutes the local form of
volume conservation. One shortcoming when assuming an incompressible flui is
that acoustic waves in the flui can no longer be described.