The results presented above are limited by the experimental conditions, such as the limited length of wave basin; as the numerical model has been validated above, the corresponding numerical simulations are employed to extend the experimental results in order to obtain a comprehensive picture of the rip current system near a shore-normal vertical structure. This allows the effect of longshore current (for the barred beaches) and other factors to be examined for more general cases. As seen for the planar beach case mentioned above, a direct result of the inclusion of longshore current is the formation of a strong offshore “exit flow”, which is the deflection rip as mentioned earlier. Here, the formation process of this rip is examined in more details and the focus is on the influence of different factors involved, including the wave reflection strength of the structure, wave incident angles and bottom topographies (with or without sandbar and different rip channels). In order to be closer to the real world condition, the simulations are conducted for a prototype case, and are performed by mapping the water depth and wave conditions in the present experiment to the prototype with the scale 1:16. Therefore, the prototype case considered has the water depth , the wave height and the period , the normal case present in the field. This wave condition corresponds to only one set of experimental wave height and period, and . The effects of different wave heights and periods are not considered due to the limited paper scope.
In the computation, the downstream boundary still corresponds to the experiment, the shore-normal structure with the length
= 363 m of the prototype. The upstream boundary is the continuous flow boundary, located at
y = 1200 m far enough from the structure in order to form a fully-developed longshore-uniformed longshore current outside the reflection wave area. The wave making line is set at
x = 406.4 m, according to
Figure 4. The other settings and considerations for the simulations include: the grid sizes in
x- and
y- directions are 0.8 m and 1.6 m, respectively, and the time step is 0.04 s. Therefore, there are 621 grid points in the
x-direction, 751 grid points in the
y-direction and a total of about 500,000 grid points in the computation domain. The simulation duration is taken as 1200 s and the time series of 240–1200 s (160 wave periods) was used to obtain the time-averaged values (for wave height, surface elevation and flow velocity). The results obtained over such a time period correspond to the first development stage of the rip current system, which can generally represent the features of the present rip current system. With these computation costs considered, the terminal times required for the present N2D4-GPU model were found about 1.8 h to complete the simulation duration 1200 s.
5.1. Effects of Wave Reflection
In the present experiment, the wave reflection strength of the shore-normal structure is the total reflection corresponding to
in (23) (
in the numerical simulation). In the simulations here, this boundary condition is extended to the other two wave reflection strengths,
and
; the resulting differences in the computed wave height distribution are shown in
Figure 9, in which the cross-shore distributions of wave heights along the 2nd and 3rd anti-node lines for these three wave reflection strengths are presented. It is seen in the figure that the differences are remarkable: with
reduced from 1000 to 4.4, the wave heights are reduced by about 30% in average and from 1000 to 3.6 by about 50%. Here, these three wave reflection strengths,
, 4.4 and 3.6, are used to represent the three different wave reflection cases: the total reflection (
), the stronger reflection (
) and the weaker reflection (
).
Figure 10 shows the simulated results, wave heights, mean water levels and corresponding mean velocity fields, for the above three wave reflection strengths on the planar beach. The results are given for a fixed incident angle 30°, with the other angles considered in the following subsection. It is seen that for the three wave reflection strengths two different rip current flow patterns appear. With
, the flow pattern forms with only the deflection rip current, the concentrated offshore flow appears through the deflection of incoming longshore current. The incoming longshore current can penetrate into the total wave reflection area. This type of deflection rip is just similar to that studied in the previous researches [
24,
25,
26,
27,
28,
29]. With
increasing to 4.4 and 1000, the rip flows at the node points in the reflection wave area begin to appear and the incoming longshore current is now deflected at the location far away from the structure with the penetration extent into the reflection wave area dependent on the strength of the wave reflection. For the stronger reflection
, the deflection appears in the central region of reflection wave area, with two node points being overlapped by the penetration current; for the total reflection
, it appears near the boundary of reflection wave area, with only one node point being overlapped by the penetration current.
To describe the above deflection rip quantitatively, the offshore flow volume rate of it is examined here, as it is an effective measure of the flow strength of the “exit flow” of the rip current system. One important thing for this examination is that it is not equal to the flow volume rate of incoming longshore current, although this rip flow is formed by the offshore deflection of this longshore current. This means that there are other contributions to the flow volume rate and the sources of these contribution are the concern here, from which the formation mechanism of deflection rip flow can be found. For this examination, the longshore extension of the deflected longshore current needs to be defined first and this is illustrated in
Figure 11, which presents the longshore profiles of the cross-shore velocity along wave breaking line
(determined according to the beginning of wave breaking outside the wave reflection area). The node rip can be identified in the figure by the alternative appearing of offshore (positive value) and onshore (negative value) velocities: the extent and value of offshore velocity indicate the width and magnitude of the node rip, and the deflection rip can be identified as the flow with only a single sign, the positive sign (flowing offshore). According to this, a cross-shore dividing line which going through the point of last zero velocity,
, can be drawn as the boundary between the node rip and the deflection rip, which is called the deflection rip down boundary hereafter. The area between it and the structure is called the inner domain hereafter, as it contains only the node rips. Apart from this boundary, another boundary, the deflection rip upper boundary, also needs to be defined, and this is demonstrated by the line
, which is defined as the position beyond which the reducing of uniform cross-shore volume flux of longshore current is less than 5%. The area between the above two lines is the area of deflection rip and is called the outer domain hereafter (the area beyond
is the domain of uniform incoming longshore current). Based on the above definitions, the amount of offshore flow volume rate of deflection rip can be calculated by the depth integration of offshore (positive) velocity on the segment of line
between above two vertical lines, i.e.,
where
is the offshore flow velocity components (the positive parts of Eulerian velocity
).
Apart from the incoming longshore current, one of the other contributions to the deflection rip flow volume is the Stokes drift. For shore-normal wave incidence, the Stokes drift is balanced by the offshore mean flow (the undertow), but for the obliquely incident wave considered here, it merges into the deflection rip current, forming the horizontal two-dimensional circulation flow in outer domain. Its flow volume along the offshore boundary of outer domain can be calculated by
Another feeding flow to the deflection rip is the flow volume through the lateral boundary
(see
Figure 11 for the strong wave reflection cases
and 4.4), the mean longshore current flowing from the inner domain into the outer domain, called the lateral flow as mentioned in
Section 4. This flow volume rate can be calculated by
With this flow included, the flow volume balance in the outer domain can be realized, that is,
where
(
at
) is the volume rate through the lateral boundary
of the incoming longshore current.
Table 3 presents the values of flow volume rate components contained in (29) for different wave conditions. The values are given non-dimensionally, divided by the incoming longshore current volume rate
, in order to see their magnitudes relative to the incoming longshore current. The corresponding relative net rate
=
is also listed in the table for checking the error with using (29). Its absolute value is less than 5%, meaning this error is negligible. The table shows that the relative value of
is 2.88, meaning that the deflection rip has the volume rate which is more than two times than that of incoming longshore current. Subtracting the longshore current contribution 1.0 gives the contributions of Stokes drift and lateral flow, 1.88, with the former possessing 1.23 and the latter 0.68. This demonstrates that the lateral is as important as the Stokes drift for the formation of deflection rip for the strong wave reflection on the structure. For the stronger reflection
, these two contributions reduce to 1.03 and 0.03, respectively, leading to
reducing to 2.03. For the weaker reflection
, the Stokes drift contribution reduces to 0.51 and the lateral flow contribution is absent (because the inner domain is absent), leading to
equal to 1.16. That is, the deflection rip is almost totally from the longshore current. The above results show that the contribution of Stokes drift and the lateral flow increase with increasing structure wave reflection. The former can be explained by the theoretical formula of this flow flux through unit width,
with
(
is the fluid density and
c the wave celerity). This means that this flow volume increases linearly with squared wave height, so the stronger the wave reflection strength, the larger the wave height and the larger the contribution of Stokes drift. The latter can be explained by the more rapid decrease in longshore direction of wave height of the longshore standing wave for the stronger structure wave reflection. Moreover, more words need to be added about this contribution. As mentioned before, for the case of normal wave incidence (
) the Stokes drift also has the contribution to the offshore flow, the undertow (the flow going offshore just beneath the Stokes drift flow), and leads to a vertical circulation. However, for the oblique wave incidence (
) considered here, the undertow does not form and the Stokes drift flow merges into the longshore current and goes offshore finally as a part of the deflection rip.
Apart from making a contribution to the flow volume balance in outer domain, the lateral flow also has the dynamic effect on the outer domain: it exerts the resistance to the longshore current, preventing it from penetrating further into the wave reflection area. The description of this is given by the values of penetration rate for different wave reflections, defined as the ratio of penetration distance () to shoreline side length () of reflection wave triangle (the area covered by the reflection wave), with () measured from the end of shoreline side of reflection wave triangle to the meeting point of downstream longshore flow and upstream longshore flow. It is seen that for the total reflection , the penetration rate is 0.33, and for the stronger reflection , the rate increases to 0.56, meaning that the penetration is very sensitive to the wave reflection strength of the structure. Actually, this rate is dependent on the magnitude of the lateral flow volume . The smaller this magnitude, the smaller the resistance to the longshore current and the larger the penetration. For the weaker reflection , is smaller than zero, and this corresponds to the near nonresistance produced by the lateral flow, so the longshore current can go directly to the region very close to the structure and is deflected offshore there ( = 1). For the strong reflections and 1000, becomes 0.03 and 0.68, so the lateral flow can produce larger resistance to the longshore current, and the longshore current cannot reach to the structure, being deflected half way to the structure, with = 0.33 and 0.56, respectively.
5.2. Effects of Wave Incident Angle
To examine the effects of wave incident angle, the other two incident angles,
and 45°, are considered with the same wave height and period as before. The comparison is made with the above
case but only for the total wave reflection case, as the interaction between the incoming longshore current and the reflection wave field is strongest for this case. It is known theoretically that with the increasing or decreasing in wave incident angle, the incoming longshore current will be increased or decreased accordingly (being proportional to
[
45]). The reflection wave area will also be increased or decreased accordingly (the wave reflection angle is equal to the wave incident angle).
Figure 12 presents the flow patterns for these two wave angles. Comparing to the case of
, the reflection wave area becomes smaller and the standing wave length becomes larger for the smaller angle
, so there is only one theoretical node and one node rip appearing in the inner domain. More important is that there is a stem wave appearing (see Yoon et al. [
46] for the stem wave), which has the width about half the longshore extent of the reflection wave area. This leads to the 1st node deviating away largely from the theoretical node line by about half the standing wavelength and the flow pattern of inner domain mainly characterized by a single large-scale vortex, with the larger offshore directed velocity of the vortex serving as the node rip current. The deflection rip current is also affected largely by such a flow feature of inner domain, that is, the longshore current is deflected at the location outside the wave reflection area, with no overlapping of deflection rip current with the node rip. The lateral flow (indicated visually by the large upstream directed velocity vectors in
Figure 12a) plays significant role for this phenomenon, as seen from its merging into the deflection rip after leaving the inner domain.
On the contrary, for the larger angle
, the reflection wave area is much larger and contains a larger number of nodes (ten), the incoming longshore current penetrates into the reflection wave area a large distance and leads to three node rips being covered by the deflection rip current, with only seven node rips remained in the inner domain. Another special flow feature is the extra offshore current appearing near the structure, called structure rip hereafter. This flow is apparently not a node rip, as the location of it being near the 1st anti-node along the structure wall (the onshore flow over this anti-node disappears). This flow also leads to the 1st node rip becoming an offshore “exit flow” due to the two flows being very close to each other. As such a flow is not found for the mid large wave angle
(see
Figure 11), its appearance is apparently due to the large wave incident angle. To give the detailed description of this flow, the other larger wave incident angles up to
are considered, the corresponding longshore distributions of cross-shore velocities
at
x = 120 m together with the result of
are given in
Figure 13. The results show that the condition for appearing of structure rip is
, indicated by
being offshore directed for these cases but being directed onshore for
and
(the normal onshore flow at an anti-node). The reason for the appearance of this type of offshore flow is the shifting of 1st node rip towards the structure. The similar shifting also appears for the case
in
Figure 11. This shifting can attribute to the reduced wave reflection on the structure, as for the strong wave reflection
there is not such a shifting (as seen in the same figure). However, the shifting here is due to the stronger longshore current in the outer domain and the weaker lateral flow in the inner domain which appears companying the increasing in wave angle. This can be confirmed in
Table 3 for
: with increasing wave angle, the former increases to 27.1 m
3/s in magnitude and the latter reduces to 0.53.
Table 3 presents the corresponding flow volume rate of deflection rip current for the three wave incident angles and the total wave reflection
, including the flow volume rate of deflection rip
and its constitute components,
and
. The results show that
,
and
all decrease with increasing wave incident angles. For
increasing from
to
, the decrease in
can be explained by the larger wave height caused by the appearance of stem wave at smaller wave angle
and the increase in
by the larger longshore gradient of mean water level caused mainly by the stem wave. However, for
increasing from
to
, the decreasing in
is due to the cross-shore wave orbital velocity
decreasing with increasing wave angle for both the propagation wave and the reflection wave and the resulting longshore standing wave (from the wave theory,
is the projection on cross-shore direction of the total velocity vector and contains the factor
in amplitude). The decreasing in
is due to the smaller longshore gradient of mean water level for larger wave angle caused by the slow longshore decreasing in longshore standing wave height. The different widths of outer domain,
, for different wave angles may also be a reason. However, this is not a determinant because the difference is not sufficiently large to change the above conclusion. All these decreases in constitute components leads to the decrease in
, as the former is the feeding flow to the latter (see (29)).
5.3. Effects of Sand Bar
The discussions above are only for a simplest bottom topography, the planar beach; here, it is extended to the case of barred beach in order to examine the sensitivity of the deflection rip flow to the bottom topography. Here, the discussion is only for the barred beach without rip channel, leaving the more complex topography case, the barred beach with a rip channel, discussed in the next subsection. The barred beach has the similar form to that in the experiment: a bar of Gaussian type section with width 32 m and height 1.28 m which is superimposed on the planar beach. The merit for the discussion without rip channel involved is that the flow pattern does not contain the rip current caused by the rip channel, leaving the flow system consisting of only the deflection rip and the node rip as discussed before for the planar beach. Therefore, the effect of sand bar can be examined in a direct way and this simplifies the problem investigated.
Different from the case of planar beach, the wave motion on the barred beach is sensitive to the water level: the location of breaking point measured from shoreline will be different for different wave levels (this is not the case for planar beach, for which the breaking point is the same for different wave levels observed from shoreline). Therefore, three water levels are chosen here in order to illustrate the effect of water level on rip current flow pattern, which are 7.2 m, 8.0 m and 8.8 m, representing low tide, mid tide and high tide, respectively. The corresponding locations of sandbar crests from the still water shoreline are at
x = 80 m, 112 m and 144 m, respectively. The corresponding breaking points (at the 2nd anti-node) are at the offshore side of the bar (at
x = 90 m), the bar crest (at
x = 112 m) and the planar slope beyond the bar (at
x = 70 m) for the above three water levels, respectively.
Figure 14 presents the mean flow velocity vectors together with corresponding wave height field and mean water level for the three water levels. As the different water levels lead to different extents of surf zone and, thus, the different longshore currents,
Figure 15 shows the cross-shore profiles of the incoming longshore currents corresponding to the three water levels. It is seen in
Figure 15 that, different from the low water level, the longshore currents for the mid and high water levels do not appear over the bar region but appear over the planar beach region close to shoreline. This is because the breaking points shift shoreward beyond the bar region for the latter two levels, as mentioned above. This means that the effect of sand bar on the incoming longshore current is weaker for the two higher water levels, which will be further discussed in the following section. With the above differences in wave breaking points, the offshore boundary of the outer and inner domains is defined a little bit differently for different water levels: the offshore boundary is selected at the bar crest for the low and mid water levels, whereas it is at the wave breaking point (
x = 70 m) for the high water level.
It is seen in
Figure 14 that for different water levels the flow patterns of rip currents have large differences. The obvious difference is that for the low water level, the offshore flowing mean current over the bar region outside the wave reflection region is continuously distributed along the bar crest, but for the mid and high water levels, the flows become concentrated distribution, forming the deflection rips, as in the case without the bar (planar beach). For the latter two rips, the difference also exists, the major one is in the rip flow locations. For the mid water level, the longshore current is deflected outside wave reflection area (a distance from wave reflection area), being unable to penetrate into wave reflection area, whereas for the high water level it is deflected within reflection wave area, being able to penetrate into wave reflection area, similarly to the case for planar beach shown in the upper panel in
Figure 10. This failure of the penetration for the mid water level is because the sand bar effect is stronger than the high water level and this leads to the strength of upstream longshore lateral flow increasing from the high water level to the mid water level. This increase can be seen visually from the longshore extent of lateral flow: the lateral flow extends upstream to
x = 310 m for the high water level but increases to
x = 448 m for the mid water level. This increase is a reflection of the bar effect: For the mid water level, the bar effect is stronger, as the wave breaking is over the bar. For the high water level, the bar effect is weaker, as the wave breaking is not over the bar (beyond the bar), at a topography that is similar to the plane slope, the bottom without the bar. Therefore, the bar strengthened the lateral flow. In fact, the appearance of the continuous offshore flow instead of a concentrated offshore flow for low water levels is also a bar effect. The bar effect is stronger than the mid water level (for the latter, the concentrated offshore flow appears, as mentioned previously). This strong bar effect leads to the longshore extent of lateral flow increasing to
x = 563.2 m, and leads to the offshore deflection of longshore current occurring over large distance over the bar, which forms the continuous deflection flow.
Table 4 lists the relative offshore volume flux of deflection rip
and its components,
,
and
for the three water levels. It is seen that
for the low water level has a smaller value comparing to that of planar beach (about 20 m
3/s), but for the mid water level, this value increases to near 35 m
3/s, whereas for the high water level, it decreases to near 26 m
3/s, which is still larger than that for the planar beach. This variation trend is related to the difference in velocity profile of longshore current for different water levels, as shown in
Figure 15.
Table 4 also demonstrates that although the value of
does not keep increasing with decreasing water level, its relative value and its components
,
and
do so. The reason for this is because the width of deflected current increases with decreasing water level, as indicated by the value of
in the table.
5.4. Effects of Rip Channel
To investigate the channel effect on the rip current flow pattern, the above analysis is extended by cutting a channel through the bar to allow the mean flow pattern to be influenced by the channel. The existence of the channel reduces the mean water level and wave height over the channel area, thus producing the time-mean offshore flow (the channel rip) converging towards the channel. Four channel widths are chosen to examine the channel rip, which are , , and ( is the wavelength of the longshore standing wave), respectively. The central lines of all the channels are located at the 3rd theoretical anti-node point (87.2 m away from the structure). As the anti-node is the location where the onshore flow will occur, the location of the rip channel chosen above will invert the flow by inducing a channel rip. This is because the decrease in the mean water level and wave height induced by the channel was larger than that owing to the location being an anti-node. Only the low water level is considered because the effect of the bar in this case is the strongest among the three water levels, as seen in the previous section.
Figure 16 shows the simulated velocity field, vorticity field, and mean water level for different channel widths. The presence of a rip channel will lead to an incoming longshore current deflecting offshore directly along the channel. Therefore, the boundaries between the inner and outer domains (
) must be selected at the channel edge near the structure.
Compared to the rip currents on the barred beach without a rip channel, an obvious flow feature is the appearance of the channel rip. The volume flow rate through the channel was calculated using the following formula:
where
and
are the locations of the channel’s two lateral boundaries. Now the corresponding deflection rip should include the channel rip, as this “exit flow” now contains two parts, the part flowing through the channel and the part flowing over the bar crest, which is present for the case without the channel. Therefore, the definition of a deflection rip (20) needs to be changed by setting the inner domain boundary to be
(as mentioned above). To see the effect of channel width on the deflection current,
Table 4 presents the flow volume rates of deflection rip
and its components,
,
and
for different channel widths. It is seen that
and
do not increase significantly with increasing the channel width. Especially for the channel rip, its flow volume rate does not increase when the channel width is larger than
. This is because the two main components of the feeding flow for the deflection rip, the incoming longshore current and the Stokes drift, nearly do not vary for different channel widths (the former does not depend on the presence of the channel, and the latter is nearly fixed with the small increase in the length of deflection current area caused by increasing channel width). This means that the appearance of channel rip nearly does not lead to the increase in the total amount of offshore flow volume of the whole flow system (the deflection rip).
The corresponding physical process can be found in
Figure 16. For all these channel widths, only the 1st node rip is remained, with the others disappearing. This disappearance is due to the overlap of the node rips with the channel rip. For
, the channel width is smaller than the wavelength of the standing wave, and the 2nd and 3rd node rips merge into the channel rip because they are close to the edges of the channel. The 4th node rip is directly overlapped by the longshore current that penetrated the inner domain. The large offshore flow volume discharged by the channel deflects the incoming longshore current along the channel. This also drives the formation of a new large-scale vortex centered at the channel edge near the structure, with its offshore flow serving as the rip flow in the channel. For
, the channel width is equal to the wavelength of the standing wave, similar to
, all the node rips except for the 1st one disappear owing to the penetration of the longshore current, and a new large-scale vortex centered at the channel edge near the structure forms, accompanied by the longshore current flowing into the channel. The difference from
is that the flow volume rate of deflection rip
increases from 136.8 to 140.2 m
3/s due to the increase in channel width. For these cases of channel width larger than the standing wavelength
and
, the channel rip is concentrated over a narrow region, which is much smaller than the channel width; thus, the increase in channel width does not lead to an increase in the channel rip width. This is particularly evident for
: the channel rip only has a width similar to that of
. For this largest width, the absolute and relative volumes (
and
) are even smaller than those of the width
(6.58 vs. 6.77), owing to the decrease in the component
(1.20 vs. 0.91). This is because a larger channel width causes the onshore flow of the new large-scale vortex to be within the channel and cancels the offshore flow of the new large-scale vortex during the calculation of the flow volume rate through the offshore boundary; this is particularly true for
.