International Journal of Computational Engineering Research||Vol, 03||Issue, 10||

||Issn 2250-3005 || ||October||2013|| Page 47

Modelling the rip current flow on a quadratic beach profile

Dr. Evans F. Osaisai
Department of Mathematics, Niger Delta University, PMB 071, Yenagoa, Nigeria
Corresponding address: P.O.Box 1693, Yenagoa, Bayelsa State Tue Oct 15 14:35:24 2013
























I. INTRODUCTION
Wave transformation in the surf zone is the dominant factor in the hydrodynamics of the nearshore
circulation, and consequent sediment transport. A global description in terms of the spatial variation of such
quantities such as wave action, wave action flux and wave radiation stress are the driving entities we have used
to describe the generation by waves of mean currents in the nearshore zone. Studies on the interaction between
waves and currents span nearly half a century. Mainly driven by a combination of engineering, shipping and
coastal interests, there has been much research on shoaling nonlinear waves, on how currents affect waves and
how waves can drive currents. This last aspect is the main concern in this paper.

The basis for this subject was laid down by Longuet-Higgins & Stewart(1960, 1961), who analyzed the
nonlinear interaction between short waves and long waves (or currents), and showed that variations in the
energy of the short waves correspond to work done by the long waves against the radiation stress of the short
waves. In the shoaling zone this radiation stress leads to what are now known as wave setup and wave setdown,
surf beats, the generation of smaller waves by longer waves, and the steepening of waves on adverse currents,
see Longuet-Higgins & Stewart (1962, 1964). The divergence of the radiation stress was shown to generate an
alongshore current by obliquely incident waves on a beach (Longuet-Higgins 1970).

The action of shoaling waves, and wave breaking in the surf zone, in generating a wave-generated
mean sea-level is well-known and has been extensively studied, see for instance the monographs of Mei (1983)
and Svendsen (2006). The simplest model is obtained by averaging the oscillatory wave field over the wave
phase to obtain a set of equations describing the evolution of the mean fields in the shoaling zone based on
small-amplitude wave theory and then combining these with averaged mass and momentum equations in the
surf zone, where empirical formulae are used for the breaking waves. These lead to a prediction of steady set-
down in the shoaling zone, and a set-up in the surf zone. This agrees quite well with experiments and
observations, see Bowen et al (1968) for instance. However, these models assume that the sea bottom is rigid,
and ignore the possible effects of sand transport by the wave currents, and the wave-generated mean currents.
Hydrodynamic flow regimes where the mean currents essentially form one or more circulation cells are known
as rip currents. These form due to forcing by longshore variability in the incident wave field, or the effect of
ABSTRACT
In this paper we develop an analytical theory for the interaction between waves and currents
induced by breaking waves on time-scales longer than the individual waves. We employed the wave-
averaging procedure that is commonly used in the literature. The near-shore zone is often
characterized by the presence of breaking waves. Hence we develop equations to be used outside the
surf zone, based on small-amplitude wave theory, and another set of equations to be used inside the
surf zone, based on an empirical representation of breaking waves. Suitable matching conditions are
applied at the boundary between the offshore shoaling zone and the near-shore surf zone. Essentially
we derive a model for the interaction between waves and currents. Both sets of equation are obtained
by averaging the basic equations over the wave phase. Thus the analytical solution constructed is a
free vortex defined in both shoaling and surf zones. The surf zone solution is perturbed by a
longshore component of the current. Thus the presence of the rip current cell combined with the
longshore modulation in the wave forcing can drive longshore currents along the beach. The
outcome, for our set of typical beach profile, is a description of rip currents.


KEYWORDS: Wave-current interactions, surf zone, shoaling zone, matching conditions, wave-
averaging, rip currents, radiation stress.
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 48

longshore variability in the bottom topography (Kennedy 2003, 2005 ,Yu & Slinn 2003, Yu 2006 and others).
They are often associated with significant bottom sediment transport, and are dangerous features on many surf
beaches (Lascody 1998 & Kennedy 2005).
There is a vast literature on rip current due to wave-current interactions, see the recent works by
(Horikawa 1978 , Damgaard et al. 2002, Ozkan-Haller & Kirby 2003 ,Yu & Slinn 2003 , Yu 2006, Falques,
Calvete & Monototo 1998a and Falques et al 1999b, Zhang et al 2004 and others) and the references therein.
Our purpose in this paper is to extend the wave-current interaction model examined in [5] to the quadratic depth.
In section 2 we record the usual wave-averaged mean field equations that are commonly used in the
literature. In section 3 we introduce a description of the rip current formation and examine the consequences for
both shoaling and surf zones. Then in section 4 we employ section 3 to extend our previous results see
Osaisai, E. F (2013), now to include the rip currents behavior in a quadratic depth profile. We conclude with a
discussion in section 5.

II. MATHEMATICAL FORMULATION
2.1 Wave field
In this section we recall the wave-averaged mean flow and wave action equations that are commonly
used to describe the near-shore circulation (see Mei 1983 or Svendsen 2006 for instance). We suppose that the
depth and the mean flow are slowly-varying compared to the waves. Then we define a wave-phase averaging
operator f f >= < , so that we can express all quantities as a mean field and a wave perturbation, denoted by a
tilde overbar. For instance,
.
~
= , , , + (1)
where , is the free surface elevation above the undisturbed depth ) ( = x h h . Then outside the surf zone, the
representation for slowly-varying, small-amplitude waves is, in standard notation,

. cos ) , (
~
u , a t x : (2)
Here ) . ( = t x a a is the wave amplitude and ) , ( = t x u u is the phase, such that the wavenumber k , frequency
O are given by

. = = ) , ( =
t
l k k u u O V (3)
The local dispersion relation is
H g U k k k e e tanh = , . =
2
+ O (4)


. =
2 2 2
l k where + k
Here ) , ( t x U is the slowly-varying depth-averaged mean current (see below), and ) , ( ) ( = ) , ( t x x h t x H , + .
To leading order, the horizontal and vertical components of the wave velocity field are respectively

. sin
sinh
) ( sinh
~
, cos
sinh
) ( cosh
~
u
k
k
u
k
k
k H
h z
a w
H
h z
a
k
u
+
O
+
O : : (5)

Importantly, note that we have ignored here any reflected wave field, which is assumed to be very
weak when the bottom topography is slowly varying. The basic equations governing the wave field is then the
kinematic equation for conservation of waves

, 0 = e V +
t
k (6)
which is obtained from (3) by cross-differentiation, the local dispersion relation (4), and the wave action
equation for the wave amplitude
. 0 = ) ( A c A
g t
V + (7)
Here e / = E A , where /2 =
2
ga E is the wave energy per unit mass, and
) / = ( , / = = k e k e d d c k c U c
g g k g
+ V is the group velocity. Using the dispersion relation (4) in (6) we
get
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 49

, = O V V +
ex g t
k c k (8)
. =
et g t
c O O V + O (9)
Here the subscript e denotes the explicit derivative of ) , , ( t x k O with respect to either x or t , when the
wavenumber k is held fixed. In this water case these explicit derivatives arise through the dependence of O on
the mean height H and the mean current U .

2.2 Mean fields
The equations governing the mean fields are obtained by averaging the depth-integrated Euler
equations over the wave phase. Thus the averaged equation for conservation of mass is
. 0 = ) ( HU H
t
V + (10)
Here , + h H = where ) ( = x h h is the time-independent undisturbed depth. For the velocity field we
proceed in a slightly different way, that is we define
, =
'
u U u + (11)
where we define U so that the mean momentum density is given by
, > < = = dz u HU M
h
}

,
(12)
But now we need to note that
'
u does not necessarily have zero mean, and that U and u are not necessarily
the same. Indeed, from (11) and (12) we get that
. 0 >= < , > < = dz u and u U u
'
h
'
}

+
,

But ) ( =
2
a O u u
'
+ , so that > <
'
u is ) (
2
a O and it follows that, correct to second order in wave amplitude,
. >= ) ,0, ( >=< < = , =
k
,
k
c
E
t x u u H M where M Hu M
'
w w
+ (13)
The term
w
M in (13) is called the wave momentum.
Next, averaging the depth-integrated horizontal momentum equation yields (Mei 1983)
. > ) = ( < > < . = ) .( ) ( h h z p dz pI u u HUU HU
' '
h
t
V + + V V +
}

,

Next an estimate of the bottom pressure term is made by averaging the vertical momentum equation to get
. > < > < . = ) ( > ) = ( <
t
h h
dz w dz u w h g h z p
} }

+ V +
, ,
,
For slowly-varying small-amplitude waves, the integral terms on the right-hand side may be neglected, and so
) ( > ) = ( < h g h z p + ~ , . Using this in the averaged horizontal momentum equation, and replacing the
pressure p with the dynamic pressure ) ( = , + z p q yields

, V V V + gH S HUU HU
t
. = ) .( (14)

. >
~
2
< > ] [ =<
2
I
g
dz qI uu S where
h
,
,
+
}

(15)

Here S is the radiation stress tensor. In the absence of any background current, so that U is ) (
2
a O ,
we may use the linearized expressions (2, 5) to find that

. ]
2
1
[ I
c
c
E
E
k c S
g
g
+ ~
e
(16)
where the phase speed k e / = c , correct to second order in the wave amplitude.

Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 50

2.3 Shoaling zone
These equations hold in the shoaling zone outside the surf zone (defined below). In summary, the
equations to be solved are that for the conservation of waves (6) combined with the dispersion relation (4), the
wave action equation (7), the averaged equation for conservation of mass (10) and the averaged equation for
conservation of horizontal momentum (14), where the radiation stress tensor is given by (16). In this shoaling
zone, we assume that wave amplitudes are small, and that there is no background current. Then all mean
quantities are ) (
2
a O , and in particular we can systematically replace H with h throughout these equations.
Next we shall suppose that ) ( = x h h depends only on the offshore co-ordinate 0 > x , where the
undisturbed shoreline is at 0 = x defined by 0 = (0) h . Further, in the near-shore region, including all of the
surf zone, we shall assume that 0 >
x
h . Then we seek steady solutions of the equation set in which all variables
are independent of the time t , and are also independent of the transverse variable y . It then follows from the
mean mass equation (10) that HU is constant, and since 0 = H at the shoreline, it follows that we can set
0 = U everywhere. Then in the dispersion relation (4) e = O . From the equation for conservation of waves
(6) we see that the frequency e and the transverse wavenumber l are constants, and the the offshore
wavenumber k is then determined from the dispersion relation (4). As is well-known, it then follows that as
0 H , | | k , that is the waves refract towards the onshore direction, where we assume that the waves
are propagating towards the shoreline so that 0 < k . The wave action equation (7) reduces to
g
Ec is constant.
Near the shore, we can assume that the shallow water approximation holds and then
1/2
) ( gh c
g
~ , so that
,
1/2
0
2
0
1/2 2
h a h a ~ (17)
where
0
a is the wave amplitude at a location offshore where
0
= h h . The surf zone ) ( = < , <
b b b
x h h h x x
can then be defined by the criterion that
b
h is that depth where
cr
A h a = / , defining an empirical breaking
condition. A suitable value is 0.44 =
cr
A , see Mei (1983) or Svendsen (2006).
The last step is to find the wave set-up , from the mean momentum equation (14), which here
becomes
, , 0 =
x x
gH S , + (18)

. )
2
1
(
cos
=
2
E
c
c
E
c
c
S where
g g
+ |
Here | is the angle between the wave direction and the onshore direction, and S is the xx component of
the tensor S . As 0 h , /2 3 , 0, E S c c
g
~ ~ ~ | , and we recover the well-known result of a wave set-
down in the shoaling zone
.
4
=
4
=
3/2
1/2
0
2
0
2
h
h a
h
a
, (19)
Here we have assumed that , is zero far offshore. Note that the first expression for , does not need the use
of the shallow water approximation, as shown by Longuet-Higgins and Stewart (1962).

2.4 Surf zone
In the surf zone
b b
h h x x < < ,0 < < 0 , we make the usual assumption (see Mei (1983) for instance)
that the breaking wave height a 2 is proportional to the total depth H , so that
,
8
= , = 2
2 2
H g
E or H a

(20)
Here the constant is determined empirically, and a typical value is 0.88. = . To determine the mean
height , + h H = , we again use the mean momentum equation (14), but now assume that
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 51

/2 = /2 3 =
2
gH E S I where /8 3 =
2
I , so that
,
) (1
= , 0 = ) (
I +

+ + I
b
b x x
h h
H H that so h H H HH (21)
where the constant
b b b
h H , + = is determined by requiring continuity of the total mean height at
b
x x = .
Note that using (19)
, /4 =
3/2 1/2
0
2
0 b b b
h h a h H
and since
b
H must be positive, there is a restriction on either the deep-water wave amplitude
0
a or on the
breaker depth
b
h ,
. / < /4
5/2
0
5/2 2
0
2
0
h h h a
b
(22)
Note that the expression (21) is valid for any depth ) ( x h , although in the literature it is often derived only for a
linear depth profile x h o = .
We are now in a position to determine the displaced shoreline
s
x x = , defined by the condition that
0 = H . That is, if ) ( =
s s
x h h then I + )/1 ( =
s
h h H , where
, ) (1 =
b b s
h h , I + I (23)
Note that to use the expression (23) it may be necessary to extend the definition ) ( x h into 0 < x . For instance
for a linear beach, x h o = , this is straightforward, but for a quadratic beach profile,
2
= x h | , the extension
for negative x should be
2
= x h | say. Note that from (19),
,
4
=
3/2
1/2
0
2
0
b
b
h
h a
,
and, on combining this with the condition (22), it follows that the shoreline recedes (advances), that is
0) 0(> <
s
h when
. / <
4
<
1
,
1
<
4
5/2
0
5/2
2
0
2
0
5/2
0
5/2
5/2
0
5/2
2
0
2
0
h h
h
a
h
h
or
h
h
h
a
b
b
b
I +
I
I +
I
(24)
This anomalous result does seem only to have been recently noticed [5] [5] and [5]. Since there is an
expectation that the shoreline should advance (see Dean and Dalrymple (2002) for instance), essentially it states
that the present model is only valid for sufficiently small waves far offshore, defined by the first inequality in
(24), which slightly refines the constraint (22). Next we can normalized the generalized bottom profile
. 1 = I +
s
h h H (25)
so that it can be written in the form
). (
1
1
=
b
s
b b
h
h
h
h
h
H

I +
(26)
Thus all the depths profiles we have examined can be plotted as normalized functions.

2.5 Linear depth
Now for the linear depth, this is just, for our parameters ) / / )( 1/(1
b s b
x x x x I + where

5/2
5/2
0
2
0
2
0
4
) (1 =
b b
s
h
h
h
a
x
x
I + + I (27)

and so depends on the wave input
2
0
2
0
/ h a and the ratio . /
0 b
h h The normalized plots are shown in
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 52

figure [1] below. It shows wave height for linear and quadratic beaches. Observe, in linear case the slope is,
) /(1 I + o , and hence is smaller than the undisturbed slope, o . The plot is a function of
b
x x/ from
b s
x x / to
1 where we have evaluated
b s
x x / . Hence the plots depend on these two parameters, and we choose say
0.1 = /
0 0
h a and 1/4 = /
0
h h
b
(gives 0 <
s
x ), which gives 0.2. = /
b s
x x

2.6 Quadratic depth
The general bottom depth profile is given by equation (25). Our normalized parameters are
). ) / ( / | | )( 1/(1
2 2
b s b
x x x x x I + Thus equation (27) is the same as that above for the linear depth case, except
that the left-hand side becomes
.
4
) (1 =
5/2
5/2
0
2
0
2
0
2
2
b b
s
h
h
h
a
x
x
I + + I
Here the depth is
2
x | where the sign is for 0 <
s
x , or 0 >
s
x which can be written as . | | x x | It follows
that for the same parameters the right-side is again 0.2 , and so 0.45. = /
b s
x x Figure [1] shows that they
agree at 0 = x as the expression (25) already shows, but that the linear depth gives a greater setup in 0 > x ,
but is weaker in the region 0 < x .


Figure 1: Plot of normalized wave height given by equation (25) for linear and quadratic profiles which
depends on the ratios 0.1 = /
0 0
h a and 1/4. = /
0
h h
b
The values of x ranges from 1 < / < 0.2
b
x x for the
linear depth and 1 < / < 0.45
b
x x for the quadratic depth profile. Thus the graphs are plotted in the range
1 < / < 0.45
b
x x .

III. THE GENERAL DESCRIPTION FOR THE RIP CURRENT FORMATION
Here we consider a steady-state model driven by an incident wave field which has an imposed
longshore variability. The wave field satisfies equation (7) which in the present steady-state case reduces to
. 0 = ) sin ( ) cos ( u u
g g
Ec y Ec x cc + cc (28)
Here we again assume that ) ( = x h h and that consequently the frequency e and the longshore wavenumber
l are constants, while the onshore wavenumber K is then determined from equation (4). We have the wave
energy E of the form

, ) ( sin ) ( cos ) ( =
0 0 0
x G Ky x F Ky x E E + + (29)
where the longshore period K / 2t is imposed. These equations in the shoaling zone yields
0 = sin ) cos (
0 0
u u
g x g
c KF c E + (30)
0 = sin ) cos (
0 0
u u
g x g
c KE c F (31)
0. = ) cos (
0 x g
c G u (32)

Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 53

on collecting terms in Ky cos , Ky sin and the constant term, which form three equations for
0
E ,
0
F
and
0
G . Equation (32) easily yields that . = cos
0
constant c G
g
u In shallow water, we may approximate by
putting
1/2
gh c
g
~ and 1 cos ~ u , so that then
1/2
0
/ h constant G ~ . For the remaining equations we can use
Snells law,
b b b
c c o u u = / sin = / sin (the constant value, here evaluated at the breaker line), and the shallow-
water approximation to get that

, 0 = } )/ ) {(
2
0
2 2
0
c E K c c E
b x x
o + (33)
while although
0
F satisfies the same equation, once
0
E has been found, then
0
F is given by either (30) or
(31). In practice, 1 << Kc and so approximately we can assume that constant c F E ~ ) , (
0 0
, the usual
shallow-water expressions. Note that here gh c ~ . In the surf zone, the expressions ) ( ), ( ), (
0 0 0
x G x F x E is
determined empirically.
Once the expression (29) has been determined, we may then substitute into the expressions (30,31 &
32) to obtain the radiation stress fields. Our aim here then is to describe how steady-state rip currents are forced
by this longshore modulation of the incident wave field, especially in the surf zone.
The forced two-dimensional shallow water equations that we use here are characteristic of many
nearshore studies (Horikawa 1978 , Damgaard et al. 2002, Ozkan-Haller & Kirby 2003 ,Yu & Slinn 2003 , Yu
2006, Falques, Calvete & Monototo 1998a and Falques et al 1999b, Zhang et al 2004 and others). Then,
omitting the overbars as before, then equations (14) in the present steady-state case reduce to

], [ = ] [
x
x H g y U V x U U H t , c c c c + c c (34)
], [ = ] [
y
y H g y V V x V U H t , c c c c + c c (35)
where the stress terms are defined;

. = =
22 21 12 11
y S x S and y S x S
y x
c c + c c c c + c c t t (36)

Next we observe that equation (10) can be solved using a transport stream function ) , ( y x , that is
,
1
=
1
= x
H
and y
H
U c c c c (37)

Next, eliminating the pressure, we get the mean vorticity equation

x
y
y
x
x y y x
H H H H
] [ ] [ = ) ( ) (
t
t

O

O
(38)
where O is define as
. ) ( ) ( = =
y
y
x
x
y x
H H
U V

+ O (39)

We shall solve this equation (38) in the shoaling zone
b
x x > and in the surf zone
b
x x < , where as
before
b
x x = is the fixed breaker line. It will turn out that the wave forcing occurs only in the surf zone, but
continuity implies that the currents generated in the surf zone must be continued into the shoaling zone.

3.1 Shoaling zone
In
b
x x > we shall assume that h H ~ as , is ) (
2
a O . Then we shall use the expressions [30 ,31]
to evaluate the radiation stress tensor. For simplicity, we shall also use the shallow-water approximation that
gh c c
g
~ ~ , and so we get that
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 54

)
2
1
sin
( = , cos sin = = , )
2
1
cos
( =
2
22 21 12
2
11
+ + u u u u E S E S S E S (40)

These expressions are in principal known at this stage, and so we can proceed to evaluate the forcing
term on the right-hand side of (38). To assist with this we recall Snells law

b b
h h u u sin = sin
where
b
h and
b
u are the water depth and incidence angle at the breaker-line. Now the energy equation (28) has
the approximate form
, 0 = ) sin ( ) cos (
y x
c E c E u u +
and using Snells law, this can be written as
0, = ) cos sin ( )
cos
(
2
c
c
E E E
x
y x
+ + u u u

.
2
1
=
c
c
E E so and
x
x x
t
We can also deduce from (28) that
, 0 = )
sin
( ) cos sin (
2
y x
E E u u u +

.
2
1
=
y y
E so and t
We can now evaluate the right-hand side of (38), and find that its identically zero,
0. = ] [ ] [
x
y
y
x
h h
t
t

Thus in the shoaling zone there is no wave forcing in the mean vorticity equation, although of course there will
be a mean pressure gradient. However, this does not concern us since here our aim is to find only the flow field.
Note that the result that there is no wave forcing in the vorticity equation does not need the specific form (29),
and is based solely on the steady-state wave energy equation (28). The specific form (29) is only used in the surf
zone.
With no forcing term, the vorticity equation (38) can be solved in the compact form, noting that we
again approximate H with h ,
. ) ( = F
h
O
(41)

But here 0 = ) ( F from the boundary conditions in the deep water as x , where the flow field
is zero. Thus our rip current model has zero vorticity in the shoaling zone. It follows that we must solve the
equation

0, = )
1
( )
1
( =
y y x x
h h
+ O (42)

in
b
x x > . Since ) ( = x h h we can seek solutions in the separated form
) ( ) ( = y Y x X (43)

with the outcome that
0. = , 0 = ) (
2
2
Y K Y
h
X K
h
X
yy x
x
+ (44)
We note the separation constant L K / 2 =
2
t must not be zero, and is in fact chosen to be consistent with the
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 55

modulation wavenumber of the wave forcing. Without loss of generality, we can choose
. sin = Ky Y (45)

For each specific choice of ) ( x h we must then solve for ) ( x X in
b
x x > , with the boundary
condition that 0 X as x . We shall give details in the following subsections. Otherwise we complete
the solution by solving the system (38) in the surf zone, and matching the solutions at the breakerline,
b
x x =
where the streamfunction must be continuous, and in order to have a continuous velocity field we must also
have that
x
is continuous.

3.2 Surf zone
To make sense of wave forcing, we assume that the expression (29) holds in this region. The functions
) ( ), ( ), (
0 0 0
x G x F x E are then determined empirically. To determine the wave forcing term in the mean
vorticity equation (38) we shall assume that 1 << =
b
u u so that, on using (??) and (40) we get that
.
2
1
= ,
2
3
=
y y x x
E E t t
Then (38) now becomes, where we again approximate H with ) ( x h ,

,
) (
=
2
=
~ ~
3/2
1/2
2
h
E h
h
h E
h
E
x y y x
x y x y xy
+ c O c c c c O c c c (46)

where here h / =
~
O O is the potential vorticity. Since the wave forcing is given by (29), that is
, ) ( sin ) ( cos ) ( =
0 0 0
x G Ky x F Ky x E E + + (47)

we observe that the unmodulated term ) (
0
x G plays no role here at all, although of course it will
contribute to the wave setup. In order to match at
b
x x = with the expression (45) for the streamfunction in the
shoaling zone, we should try for a solution of (46) of the form
. < , ) ( sin ) ( =
b
x x in x G Ky x F + (48)

The matching conditions for the streamfunction and velocity field at the breakerline
b
x x = require
that
. 0 = ) ( , 0 = ) = ( , ) ( = ) ( , ) ( = ) (
b x b b x b x b b
x G x x G x X x F x X x F
The expression (48) yields
G Ky F
~
sin
~
= + O (49)

where F
~
and G
~
are differential operators where they are defined as;

h
F K
h
F
F
x
x
2
) ( =
~
(50)

h
G
Z Z G
x
x
= , =
~
(51)

From equation (46) we get a set of three equations that are used to determine the rip-current flow field
in the surf zone. These are namely;
, 0 = )
~
(
~
x x
h
F
F
h
F
F (52)
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 56

,
) (
= ) )
~
(
~
(
3/2
0
1/2
h
F h
h
G
F
h
F
G
x
x
(53)
. 0 =
)
(
3/2
0
1/2
h
E h
x
(54)

Equation (54) gives
1/2
0
1/ h E : , which is an unacceptable singularity as 0 h . Hence we must infer
that in the surf zone at least, 0 =
0
E . The first of the three equations, that is (3.2a) suggests that

, =
~
onstant whereCisac CF
h
F
(55)
and the second (3.2b) yields that

3/2
0
1/2
) (
= ) )
~
( (
h
F h
h
G
CG F
x
x
(56)

The boundary conditions at 0 = x where 0 = h are that both mass transport fields V U , should
vanish, that is from (37) constant = and 0 = / h
x
, which implies that
. 0 = , 0 = , = , 0 = = x at
h
G
constant G F F
x
x
(57)
As above there are also the matching conditions for both F and G separately at the breakerline, that is for F
we have that

b
b
b x
b
b x
x x at
x X
x X
x F
x F
= ,
) (
) (
=
) (
) (

where we note that here the right-hand side is a known quantity, depending only on K and
b
x . Next we see
that equation (55) reduces to

ChF
h
F K
h
F
x
x
= ) (
2
(58)

Together with the boundary conditions at
b
x x x = 0, = this is essentially an eigenvalue problem for
) ( x F with eigenvalue C . In general it is solved approximately since we shall assume that 1 <<
b
Kx . Once
) ( x F is known we can solve (56), together with the appropriate boundary conditions to get ) ( x G to complete
the solution.
Note that the amplitude of ) ( x F is an arbitrary constant in this solution, and so we can fix it by
specifying its value at
b
x x = say. Indeed the solution we have constructed is essentially a free vortex defined
by Ky x X sin ) ( in the shoaling zone
b
x x > , and Ky x F sin ) ( in the surf zone
b
x x < , perturbed by a
longshore component ) ( x G in the surf zone. Note that in the presence of the wave forcing, both G F , are
non-zero, see (56). It is significant that unlike the longshore currents considered in the basic state which depend
on an ad hoc frictional parametrization, the presence of the rip current cell combined with the longshore
modulation in the wave forcing can drive a longshore current.

IV. APPLICATION TO A QUADRATIC DEPTH PROFILE
Osaisai, E.F (2013) examined the behavior of the rip-current for the case . = x h o We here extend the
case for which
2
= x h | . In the basic state we show that the linear depth gives a greater setup in 0 > x , but is
weaker in the region 0 < x , where both depths agreed at 0 = x .
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 57


4.1 Shoaling zone
Now we let
2
= x h | and then equation (44) now admits the differential equation of the form
0 = ) ( ) ( 2 ) (
2
x xX K x X x xX
' ' '
(59)

whose solutions are explicitly the exponential functions of the form
1). ( 1) ( = ) (
2 1
+ +

Kx e C Kx e C x X
Kx Kx

From the behavior of the solutions as x we see that 0 =
1
C and so
. 1) ( = ) (
2
+

Kx e C x X
Kx
(60)
Again note that the constant | does not appear in this solution.

4.2 Surf zone
Similarly, let
2
= x h | in (58), so that proceeding as for Osaisai, E.F (2013) we obtain the equation

, =
2
7
0
2
x A F K F
x
F
x xx
(61)

where now
2
= | C , and as 0 x ,
3
0
x A F ~ . The solutions of the homogeneous equations are
known, these are 1) ( =
1
+

Kx e j
Kx
and 1) ( =
2
Kx e j
Kx
where
1
j and
2
j are linearly independent
solutions. The general solution is given by

, ) ( ) ( = ) (
2 2 1 1 2 1
j C j C j x B j x A x F + + + (62)

where the Wronskian of
1
j and
2
j is
2 3
2 = x K W , and
,
12 2
= ) (
3
6
0 5
2
0
3
0
K
x A
dx x j
K
A
x A
x

~
}


.
12 2
= ) (
3
6
0 5
1
0
3
0
K
x A
dx x j
K
A
x B
x

~
}

Thus A and B vanish at 0 = x as
5
x and so
2 1
= C C . Also the normalization of F as 0 x gives
2
0 1
/2 3 = K A C .
Now to apply the boundary condition at
b
x x = which yields , we simplify the calculation by first
approximating A and B as above for small x . The outcome is that
. )
36
(
9
3
0
x
x A F

+ ~
The boundary condition at
b
x x = again yields an explicit equation for which can be simplified by the
assumption that 1. <<
b
Kx Thus scales as
6
b
x which may not be so good an approximation. In spite of
that the RHS may be approximated by .
2 2
b
x K Hence to leading order

. 12 =
6

b
x (63)
This leads to the simple expression
, )
3
(1
6
6
3
0
b
x
x
x A F ~ (64)
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 58

and evidently 0 ~
x
F at
b
x x = .



Figure 2: Plot of
0
)/ ( A x F against
b
x x/ , where
0
A is arbitrary as given by equation (64). Observe there is no
dependence on the slope | but only a weak dependence on K . The value of ) ( x F also reaches a maximum
value of 0.7. )/ (
0
~ A x F This shows that irrespective of ), ( x h the maximum value ) ( x F can admit is
0.7 ~ .

As in [5] we can now add correction terms, letting /12 1 =
6
b
x c + . Expanding A and B we find
that
. )
9
2
36
1
[
72
1
4
1
( = ) (
11 2 9 6 6 4 4 3
0
x K x x K x K x A x F + +
As before, we now find the leading order term for c ,
. /12 1 =
6
b
x c +
Finally we get that
).
3 3 4
1
(1 = ) (
5
6
6
6
4 3
0
b b
x
Kx
x
x
x K x A x F (65)
Next, as before, we need to solve for ) ( x G from (56). As above, we approximate
3
0
= x A F , and also we use
the empirical expression /8 =
2 2
0
h F , see equation (20). Thus we get

0
2 2
4 2
8
5
=
2
A
g
Z x C Z
x
Z
x xx
|
| (66)

Letting
3
= x u we get that
.
8
5
= 9
3
4
0
2 2
u A
g
Z Z
uu
|

As before the dominant balance in the particular solution is between
uu
Z and the right-hand side, so that
2/3
u constant Z
p
~ . We find that
.
16
45
=
0
2 2 2
A
x g
Z
p
|
(67)

Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 59

As in the linear depth case, this have the form as 0 <
.
3
cos
3
sin =
3 2
u
C
u
C Z
h

+
Here we recall that . = The total solution is then
. =
h p
Z Z Z +
The boundary condition at 0 = u gives 0. =
3
C Again imposing the boundary conditions that 0 = Z at the
breakerline
b
x x = gives,

]).
3
2
[ sin
] 3 [2/ sin
1
(
16
45
=
3
3
2
2
0
2 2 2
b
b
b
x
x
x
x
A
x g
Z
|
(68)

Finally we get G from hZ G
x
= and 0 = G at , =
b
x x
. )
] 3 [2/ tan 3 2
1
5
1
]
3
2
[ cos
] 3 [2/ sin 3 2
1
5
(
16
45
=
3
3
5
5
0
5 3 2
+
b
b
b
x
x
x
x
A
x g
G
|
(69)




Figure 3: Depicts the plots of normalized ) ( x Z and ) ( x G given by equations (68) and (69) where each is
normalized by
0
2 2 2
/16 45 A x g
b
| and
0
5 2 2
/16 45 A x g
b
| respectively with
0
A arbitrary. We observe as in
[5], here too as depicted in the figure, there is a small region of reversed flow near the breaker line.


The combined expressions (60, 64, 69) complete the solution, where we recall that the constant C is
given by (63) (since
2
= | C ), or their respective higher-order corrections. Now the amplitude of ) ( x F at
b
x x = is given by
. 1) ( = ) (
2
+

Kx e C x F
Kx
b
(70)
On using the approximation 1 <<
b
Kx , and the approximate expression (64), this reduces to
. =
3
2
= ) (
2
3
0
C
x A
x F
b
b

The rip-current system contains a free parameter
0
A or its equivalent. We choose to define this free parameter
to be the value of ) (
b
x F and normalize the full solution by this value. Thus we get from (43, 45) in
b
x x > ,
and (48, 64, 69) in
b
x x < that the normalized streamfunction
n
is given by

Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 60

, > , ) ( sin
) (
) (
=
b
b
n
x x for Ky
x X
x X
(71)
. < < 0 ,
(0)
) (
sin
) (
) (
=
b
b
n
x x for
G
x G
R Ky
x F
x F
+ (72)
Here again ) ( (0)/ =
b
x F G R is a free parameter. From (64, 69) we find that here
. )
] 3 [2/ tan 3 2
1
5
1
3 [2/ sin 3 2
1
(
32
135
=
2
0
2 2 3

A
x g
R
b
|
(73)

Note that again 0 < R , and that | | R increases as the wave forcing increases, or as the curvature |
increases, or as the depth
2
b
x | at the breaker line increases. In order to estimate typical values for R we again
note that from (64) the longshore velocity field in the ) ( sin Ky -component scales as
b c
x A V | / =
0
, while
the longshore component then scales with
c
RV . Taking account of the actual numerical values in the
expressions given above, we find that a suitable values are 0.1 ~ R . Plots of
n
are shown in figure [4 & 5
] for same values of R as in the linear case, and again 0.2 =
b
Kx .



Figure 4: Plot of the rip current streamlines for a quadratic depth profile, given by equation (72) where ) ( x F
and ) ( x G are equations (64) and (69) respectively for 0.02 = R in the left panel and 0.1 = R in the right
panel.



Figure 5: As for figure 4 but 0.5 = R in the left panel and 2 = R in the right panel.

Overall these plots show the same kind of behaviour as those for the linear depth profile see [5].
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 61

However, the major difference is that the flow in the surf zone is rather weaker, and so the vortex centre is
slightly further offshore.

V. CONCLUSION
We described qualitative solutions for rip currents which are essentially free vortices in both zones.
The free vortex in the surf zone is perturbed by a longshore modulation in the wave forcing. Rip current cell
combining with the longshore modulation in the wave forcing can drive longshore currents along the beach.
Thus the dynamics of the shoaling zone is only dependent on the state-state wave energy equation.The wave
forcing in the surf zone sets the wave activities different from those of the shoaling zone. To determine wave
forcing in the mean vorticity equation we assume that the wave angle becomes smaller. We also note here that
the component of the radiation stress in the y momentum remains unchanged across the entire flow domain.
This shows that it is only the x component of the radiation stress that play a leading role in the wave forcing.
However, wave forcing encountered in the surf zone has an unmodulated term that does not play a role in the
vorticity equation but only contribute to wave setup.To ensure continuity of the streamfunctions in the shoaling
zone we match the solution at the breakerline by a matching condition with appropriate boundary conditions.
Thus the rip currents solution in the surf zone is provided by the terms in the matching condition. The terms in
the matching condition has a cross-shore width and a modulated longshore component. The cross-shore width
was determined by the application of perturbation method and variation of parameter. It would to interesting to
examine the effect of friction on the rip currents.

REFERENCES
[1] Baldock, T.E, ( 2006), Long wave generation by the shoaling and breaking of transient wave groups on a beach, Proceedings of the
Royal society, vol.462, pp 1853-1876.
[2] Baldock, T.E., OHare, T.J. & Huntley, D.A, (2004), Long wave forcing on a barred beach, Journal of Fluid Mechanics, vol 503,
pp. 321-343.
[3] Battjies, J. A., (1988), Surf zone dynamics, Annual Rev. Fluid Mechanics, vol. 20, pp.257-293.
[4] Billigham, J & King, A.C., (2000), Wave Motion- Cambridge Texts in Applied Mathematics, Cambridge University Press.
[5] Brocchini, M., Kennedy, A., Soldini, L & Mancinelli, A., (2004) Topographically controlled, breaking-wave-induced
macrovortices. Part 1. Widely separated breakwaters, J. Fluid Mechanics, vol. 507, pp. 289-307.
[6] Buhler, O., (2000), On the vorticity transport due to dissipating or breaking waves in shallow-water flow, Journal of Fluid
Mechanics, vol 407, pp 235-263.
[7] Buhler, O. & Jacobson, T. E., (2001), Wave-driven currents and vortex dynamics on barred beaches, Journal of Fluid Mechanics,
vol 449, pp 313-339.
[8] Damgaard, J., Dodd, N., Hall, L. & Chesher, T., (2002), Morphodynamic modelling of rip channel growth, Coastal Engineering,
vol.45, pp 199-221.
[9] Davies, A. M. & Lawrence, J. (1995), Modelling the effect of wave-current interaction on the three-dimensional wind driven
circulation of the Eastern Irish Sea, Journal of Physical Oceanography, vol. 25, pp. 29-45.
[10] Falques, A., Calvete, D.& Monototo, A. (1998a), Bed-flow instabilities of coastal currents. In Physics of Estuaries and Coastal
Seas, pp. 417- 424. A. A. Balkema.
[11] Falques, A., Montoto, A. & Vila, D. (1999), A note on hydrodynamic instabilities and horizontal circulations in the surf zone,
Journal of Geophysical Research, vol. 104(C9), pp 20605-20615
[12] Falques, A., Coco, G & Huntley, D.A (2000), A mechanism for the generation of wave-driven rhythmic patterns in the surf zone,
Journal of Geophysical Research, vol. 105, No. C10, pp 24701-20615.
[13] Guza, R. T. & Davis, R. E. (1974), Excitation of edge waves by waves incident on a beach, Journal of Geophysical Research, vol.
79, C9, pp. 1285-1291.
[14] Horikawa, K (1978), An Introduction to Ocean Engineering, University of Tokyo Press.
[15] Kennedy, A.B. (2003), A circulation description of a rip current neck, Journal of Fluid Mechanics, vol. 497, pp. 225-234
[16] Kennedy, A.B. (2005), Fluctuating circulation forced by unsteady multidirectional breaking waves, Journal of Fluid Mechanics,
vol. 538, pp. 189-198
[17] Kennedy, A.B., Brocchini, M., Soldini, L., & Gutierrez. E., (2006), Topographically controlled, breaking-wave-induced
macrovortices. Part 2. Changing geometries, J. Fluid Mechanics, vol. 559, pp. 57-80.
[18] Lane, E.M & Restrepo, J.M (2007), Shoreface-connected ridges under the action of waves and currents, Journal of Fluid
Mechanics, vol. 582, pp. 23-52
[19] Lascody, L. L (1998), East central Florida rip current program, National Weather Digest, vol. 22, pp. 25-30.
[20] List, J.H. (1992), A model for the generation of two dimensional surf beat, Journal of Geophysical Research, vol. 97, pp. 5623-
5635.
[21] Longuet-Higgins, M.S. (1970), Longshore currents generated by obliquely incident sea waves, 1, J. Geophysical Research, vol. 75,
pp 6778-6789.
[22] Longuet-Higgins, M.S. (1970), Longshore currents generated by obliquely incident sea waves, 2, J. Geophysical Research, vol. 75,
pp 6790-6801.
[23] Longuet-Higgins, M.S & Stewart, R.W. (1960), The changes in the form of short gravity waves on long waves and tidal currents,
Journal of Fluid Mechanics, vol. 8, pp 565-583
[24] Longuet-Higgins, M.S & Stewart, R.W. (1961), The changes in amplitude of short gravity waves on on steady non-uniform
currents, Journal of Fluid Mechanics, vol. 10, pp 529-5549
[25] Longuet-Higgins, M.S & Stewart, R.W. (1962), Radiation stress and mass transport in gravity waves, with application ro surf beats,
Journal of Fluid Mechanics, vol. 13, pp 481-504
Modelling the rip current flow-----
||Issn 2250-3005 || ||October||2013|| Page 62

[26] Longuet-Higgins, M.S & Stewart, R.W. (1964), Radiation stress in water waves: a physical discussion with applications, Deep-sea
Research, vol. 11, pp. 529-562.
[27] Longuet-Higgins, M.S. (2005), On wave set-up in shoaling water with a rough sea bed, Journal of Fluid Mechanics, vol. 527, pp.
217-234
[28] Mei, C.C. (1983), Applied Dynamics of Ocean Surface Waves, John Wiley and Sons, Inc.
[29] Osaisai, E. F (2013), An analytical model of the rip current flow International Journal of Computational Engineering Research,
vol. 3, issue 9, pages 1-12
[30] Osaisai, E. F (2008), The interaction between waves and currents in the nearshore zone, PhD thesis, Loughborough University.
[31] Peregrine, D.H. (1998), Surf zone currents, Theoretical and Computational Fluid Dynamics, vol. 10, pp 295-309.
[32] Rahman, M. (2002), Mathematical methods with applications, WITPress, Southampton.
[33] Restrepo, J. M. (2001), Wave-current interactions in shallow waters and shore-connected ridges, Continental Shelf Research, vol.
21, pp. 1331-1360.
[34] Roger Grimshaw and Evans Osaisai (2012), Modeling the effect of bottom sediment transport on beach profiles and wave set -up,
ocean Modelling, 59-60 (2012) 24-30
[35] Shepard, F. P., Emery, K. O. & La Fond, E. C (1941), Rip currents: a process of geological importance, Translation of American
Geophysical Union, vol. 31, pp. 555-565.
[36] Tucker, M.J., (1950), Surf beats: sea waves of 1 to 5 min. period, Proc. R. Soc. London, vol. 202, pp. 565-573.
[37] Yu, J. & Mei, C.C. (2000), Formation of sand bars under surface waves, Journal of Fluid Mechanics, vol. 416, pp. 315 - 348.
[38] Yu, J. & Mei, C.C. (2000), Do longshore bars shelter the shore?, Journal of Fluid Mechanics, vol. 404, pp. 251 - 268.
[39] Yu, J. & Slinn, D.N. (2003), Effects of wave-current interaction on rip currents, Journal of Geophysical Research, vol 108 (C3),
3088.
[40] Yu, J. (2006), On the instability leading to rip currents due to wave-current interaction, Journal of Fluid Mechanics, vol. 549, pp.
403-428.