Next Article in Journal
Evaluation of TRMM Precipitation Dataset over Himalayan Catchment: The Upper Ganga Basin, India
Previous Article in Journal
A Comparison of Continuous and Event-Based Rainfall–Runoff (RR) Modelling Using EPA-SWMM
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Investigation on Characteristics of Sand Waves with Fine Sand under Waves and Currents

1
College of Engineering, Ocean University of China, Qingdao 266100, China
2
Shandong Province Key Laboratory of Ocean Engineering, Ocean University of China, Qingdao 266100, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(3), 612; https://doi.org/10.3390/w11030612
Submission received: 22 February 2019 / Revised: 20 March 2019 / Accepted: 21 March 2019 / Published: 24 March 2019
(This article belongs to the Section Hydraulics and Hydrodynamics)

Abstract

:
A series of physical experiments was conducted to study the geometry characteristics and evolution of sand waves under waves and currents. Large scale bedforms denoted as sand waves and small bedforms represented by ripples were both formed under the experimental hydrodynamic conditions. Combining the experimental data with those from previous research, the characteristics of waves and currents and measured sand waves were listed. Small amplitude wave theory and Cnoidal wave theory were used to calculate the wave characteristics depending on different Ursell numbers, respectively. The results show good agreement between the dimensionless characteristics of sand waves and the dimensionless wave characteristics with a smaller wave steepness. When the wave steepness is large, the results seem rather scattered which may be affected by the wave nonlinearity. Sand wave steepness hardly changed with bed shear stress. A simple linear relationship can be found between sand wave length and wave steepness. It is easy to evaluate the sand wave characteristics from the measured wave data.

1. Introduction

Bed forms with different scales are widely observed on many shallow continental shelves. Small scale bedforms, known as ripples, are often superimposed on the large scale bedforms such as sand waves, sand ridges etc. They coexist under interaction between tides, winds, waves and other complicated hydrodynamic conditions [1,2]. Among these bed forms, sand waves are of particular concern.
Sand waves are dynamic rhythmic bed forms with wave lengths of hundreds of meters and heights of several meters [3,4,5]. They can migrate up to tens of meters per year [6]. From an engineering point of view, migrating sand waves can affect submarine constructions and the safety of pipelines and cables [7,8]. Moreover, the generation and growth of sand waves can also decrease least navigable depths and increase the cost of dredging [7,9,10,11].
Investigating the relationship between sand wave characteristics and hydrodynamics is very important for understanding sand wave generation and evolution. The morphology of submarine barchans and sand waves off the Dongfang coast in Beibu Gulf of South China Sea has been previously studied [12,13]. Many observations indicate that the morphological patterns strongly depend on the hydrodynamical conditions [14]. Belderson [15] indicated that sand waves only occur in areas where depth-averaged tidal currents are stronger than 0.5 m/s. The ripples superimposed on the sand waves can affect the bottom current dramatically. However, there are few studies and detailed descriptions of superimposed ripples.
Stability analysis, which considers the hydrodynamics and bed forms as a coupled system, is the most prevalent method to study the evolution of sand waves. Hulscher [16] described the formation mechanism of sand waves as the free instability of the sandy seabed subject to tidal motions. In Hulscher [16], the interactions between oscillatory tidal flows and sinusoidal bed perturbations gave rise to tidally-averaged residual currents, in the form of vertical recirculating cells directing from the trough to the crest of the sand waves [16,17]. More processes were later implemented into the model, such as residual currents [18], multiple tidal components [19], wind waves [5,20,21], biology [22,23], and graded sediment [24,25,26,27]. In addition to stability models, process-based morphodynamic models, such as Delft3D, were widely applied in the study of sand waves recently. More elaborate processes such as turbulence [17], and suspended sediment transport [28,29] can be handled by these kinds of models. Some Computational Fluid Dynamics (CFD) models were also used to describe sediment transport [30]. To date, few efforts have been dedicated to study the sand wave and its migration by physical experiments, except those of Cataño-Lopera et al. [1], Williams et al. [31], and Zhu et al. [32]. Cataño-Lopera et al. [1] described ripple configurations superimposed upon larger bottom features. Equations to predict the characteristics of sand waves using measured wave data were provided. However, it is still unknown how the wave nonlinearity affects the geometry characteristics and migration rules of sand waves, and therefore this is specially investigated in this manuscript. The aims of this paper are twofold: (1) to describe the experiment results of bedforms formed under waves and currents, (2) to study the relationship between sand wave characteristics and wave configuration using different wave theories. The remainder of the paper is organized as follows: Section 2 describes the experimental set-up and calculation of different wave theories based on the Ursell number. In Section 3, the dynamic and geometry characteristics of sand waves are analyzed. Finally, Section 4 and Section 5 contain the discussion and conclusions, respectively.

2. Experimental Installation and Methods

2.1. Experimental Set-Up

A series of experiments was conducted in the wave-current flume in the Hydraulic Lab at the Ocean University of China. Figure 1 shows the schematic layout of the physical model experiment. The length of the wave flume is 60 m; the height and the width are 1.5 m and 0.8 m, respectively. The waves are generated by a piston type wave maker at one side of the flume. A porous medium is used at the other end to absorb the waves. Uniform currents with positive or negative directions are generated by two circulating pumps. The flume is capable of generating waves of heights between 0.03 m and 0.3 m and wave period between 0.5 s and 3 s. In the middle-rear of the wave flume, the bottom was covered by a sand bed 15 m in length, 0.8 m in width and 0.4 m in thickness with non-cohesive sands. Two 1:10 concrete slopes were set at the both sides of the sand bed to avoid a sudden change of hydrodynamic conditions.
The wave height and wave period are measured with two wave gauges. During the experiments, 50 Hz of sampling rate is used which is sufficient to measure the variation of water surface elevation for all the scenarios in this paper. The depth-averaged velocity is measured by ADV (Acoustic Doppler Velocimetry). The 0.6-depth method (One-point method) is used in the experiments to determine the depth-averaged velocity. A Trasonic Terrain automatic Measurement and analysis System (TTMS) is used to measure evolution of the sand bed. The resolution of the system is ±1 mm in the vertical direction.
The initial bed surface is flat at the beginning of each experiment. Under the waves and currents action, the bed elevation is measured every half an hour in the first hour, and then measured every hour. The experiment is complete when the bed achieves a dynamic equilibrium. In this paper, we assume that the bed has achieved a dynamic equilibrium when the change of bed is less than 5% of the sand wave elevation. Table 1 shows the experimental scenarios. The medium diameter in the first eight tests is 0.17 mm while it is 0.18 mm in the last two tests.

2.2. Water Characteristics Calculation

The Ursell numbers (defined as U r = H L 2 / d 3 , here L denotes wave length) computed for all 10 tests in this paper were less than 26. Therefore, small amplitude wave theory can be adopted to obtain the wave configurations. An Ursell number larger than 26 corresponds to shallow water in almost half the cases of Cataño-Lopera et al. [1], so Cnoidal wave theory must be used. Consequently, different wave theories were used according to an Ursell number of 26.
The bed shear stress is an important factor to induce the sediment incipient motion and transport. The dimensionless bed shear stress T * can be expressed as:
T * = τ c w τ c r τ c r
where τ c w is mean bed shear stress due to current and waves; τ c r is the critical bed shear stress, which is defined as:
τ c r = ( ρ s ρ w ) g D 50 θ c r
where ρ s and ρ w are the densities of sediment and water, respectively; g is gravitational acceleration; θ c r is the critical Shields number, which can be calculated according to the classical Shields curve. Here, we adopt a function of the dimensionless particle diameter D * which is proposed by van Rijn [33]:
θ c r = { 0.24 D * 1         1 < D *   4 0.14 D * 0.64         4 < D *   10 0.04 D * 0.1         10 < D *   20 0.013 D * 0.29         20 < D *   150 0.055         150 < D *  
where the dimensionless particle diameter is,
D * = D 50 [ ( ρ s / ρ w 1 ) g v 2 ] 1 / 3
where ν is the viscosity coefficient of water.
Figure 2 shows the calculation flow chart of total bed shear stress τ c w . With the measured wave period T and water depth d , the wave length is calculated with the dispersion relation,
L = g T 2 2 π tan h ( 2 π L d )
With the measured wave height H , the Ursell number can be obtained. When the Ursell number is less than 26, the amplitude of the near-bottom orbital velocity is calculated according to small amplitude wave theory,
u m = π H T 1 sin h ( k d )
where k is the wave number.
The wave-averaged bed shear stress due to waves is calculated as,
τ w = 0.25 f w ρ w u m 2
where f w is the wave-related friction factor, which is related to the amplitude of the wave particle and the ripple existing on the bed. We adopt the formulation of Swart [34],
f w = { exp ( 5.213 ( a m Δ ) 0.194 5.977 )       a m Δ > 1.57     0.30                                                                                                             a m Δ 1.57
where Δ is the wave related roughness, which is related to the estimated ripple height, Δ = 4 η s , here the ripple height η s is set as 1.5 cm; a m is the near-bed peak orbital excursion, which can be solved as,
a m = T 2 π u m
When the Ursell number is larger than 26, Cnoidal wave theory must be used to evaluate the wave configurations. Based on the Cnoidal wave theory, wave period can be calculated as [35,36]:
T = 4 d 3 g H { κ K ( κ ) 1 + H d [ 1 + 1 κ 2 ( 2 3 E ( κ ) K ( κ ) ) ] }
where κ is the modulus of the elliptic integral; K ( κ ) and E ( κ ) are the first and second types of complete elliptic integrals, respectively. With the measured data of wave period, the parameter of κ can be solved using the Newton iterative method. Then the wave length with Cnoidal wave theory is:
L = 16 d 3 3 H [ κ K ( κ ) ]
The time series of the wave particle velocity calculated with Cnoidal wave theory in shallow water is quite different from the small amplitude wave theory. The near-bed instantaneous velocity due to waves can be calculated by:
u g d = { 5 4 + 3 z t 2 d z t 2 4 d 2 + ( 3 H 2 d z t H 2 d 2 ) cn 2 H 2 4 d 2 cn 4 8 H K 2 ( κ ) L 2 ( d 3 z 2 2 d ) ( κ 2 sn 2 cn 2 + cn 2 dn 2 sn 2 dn 2 ) }
where sn ( ) , cn ( ) , dn ( ) are Jacobian elliptical sinusoid function and Jacobian elliptical cosine function and Jacobian elliptical delta function, respectively; z t is the distance from the wave trough to the seabed, the distance from the wave crest to the seabed is defined as z c .
z c = d + 16 d 3 3 L 2 K ( κ ) [ K ( κ ) E ( κ ) ]
z t = z c H = d + 16 d 3 3 L 2 K ( κ ) [ K ( κ ) E ( κ ) ] H
cn 2 = cn 2 [ 2 K ( κ ) ( x L t T ) , κ   ]
sn 2 = 1 cn 2 [ 2 K ( κ ) ( x L t T ) , κ   ]
dn 2 = 1 κ 2 sn 2 [ 2 K ( κ ) ( x L t T ) , κ   ]
Two scenarios were used to evaluate the difference of velocity between small amplitude wave theory and Cnoidal wave theory. Scenario 3 (with water depth of 0.6 m, wave height of 0.14 m and wave period of 2 s) in Table 1 was used for the first test, while an experimental scenario (with water depth of 0.56 m, wave height of 0.104 m, and wave period of 4.1 s) from Cataño-Lopera et al. [1] was used for the second test. The Ursell numbers for the two tests were 12.33 and 55.5, respectively. Therefore, the small amplitude wave theory was used for the first test and the Cnoidal wave theory was used for the second test. The results can be found in Figure 3.
With Cnoidal wave theory, the wave period-averaged bed shear stress can be expressed as:
τ w = 1 T 0 T 1 2 ρ w f w u 2 d t
It is difficult to get the analytical solution from Equation (18), so a numerical integration method was used. With Δ T = T 100 ,
τ w 1 T i = 0 100 ( 1 2 f w ρ w ( u ( i Δ T ) ) 2 ) Δ T
When the wave is combined with current, the current induced bed shear stress is:
τ c = 1 8 ρ w f c u 0 | u 0 |
where f c is the current friction factor. When 10 < 2 d / Δ < 400 , f c can be expressed as [37]:
f c = 0.113 ( Δ d ) 1 / 3
With all the equations above, the dimensionless bed shear stress can be obtained. Hydrodynamic conditions and characteristics of the measured sand waves are listed in Table 2.

3. Results

3.1. Dynamic Growth of Sand Waves

Take Scenario 4 for example, Figure 4 shows the process of bed evolution during the experiment. First, the bed is flat at the initial time with the help of some calibration lines posted outside the glass of the flume. Second, fill the water into the wave flume until the depth reaches up to 0.6 m upon the sand bed. After the wave action has lasted for several minutes, ripples are formed alternatively on the bed. Third, after about 1-hour’s wave action, the ripples can be found all over the sand bed. Sand waves start to form. Finally, as Figure 4d shows, sand waves are generated. The height of the sand waves is of several centimeters and the length of about several meters. Ripples superimposed upon the sand waves present 2D or 3D patterns depending on their location over the large bedforms. 2D patterns of ripples can often be found at the crest or trough of sand waves, while 3D patterns are often found between the crest and trough. The ripples on the crest are usually larger than those at the trough. It is quite important to investigate the characteristics of ripples at different locations. They can affect the near-bed flow field and bed shear stress dramatically by changing the roughness coefficient.
Using the TTMS, the bed level change was measured with an accuracy of 1 cm. Figure 5 shows the bottom evolution over time. With the FFT (Fast Fourier Transformation) method, the wave height can be obtained. In this case, the height of the sand wave is about 3.79 cm, and the length is about 1.78 m after 5 h. Note in Figure 5 the sand wave reaches a quasi-steady condition after 5 h when no obvious bed changes can be observed.
The growth rate of sand waves is calculated with the FFT method of bed level at different times. From Figure 6 the first half hour shows the growth rate is about 1.04 cm/h which is not very large. At this stage ripples start to appear along the bed and large bedforms are formed gradually. The growth rate reaches the largest value which is about 1.42 cm/h in the second half hour and the main geometric characteristics of sand waves begin to be formed. Then from the second to the fifth hour, the growth rate decreases gradually at about 1.01 cm/h, 0.78 cm/h, 0.50 cm/h, and 0.27 cm/h, respectively. The change of bed is less than 5% of the sand wave elevation at the fifth hour, so with the first 5 h wave action the sand waves have reached a balanced state. At an average of 10 h, the growth rate is about 0.4 cm/h.

3.2. Geometry Characteristics of Sand Waves

Dimensionless bed shear stress was used to evaluate the measured sand wave characteristics. The hydrodynamic conditions and measured data of sand waves are listed in Table 2. The measured data of the bed forms denote the long type sand waves, not the ripples superposed on the waves. Notice that Cnoidal wave theory is used when the Ursell number is larger than 26 which make the results more accurate. Figure 7 shows the relationship between dimensionless sand wave height and dimensionless shear stress. The data is rather scattered. A simple linear relationship was used to evaluate change regulation with all the data. The black line is the fitting curve and the formula is as follows.
h s w d = 0.162 0.000226 T *
Dimensionless sand wave height decreases slightly with increase of dimensionless bed shear stress. The correlation coefficient is about 0.02 which shows poor correlation between dimensionless sand wave height and bed shear stress. The data is rather scattered when the dimensionless shear stress is larger than 30, so the hydrodynamic conditions were analyzed. Among these scenarios, the result of Scenario 18 is quite different which is much larger than the others. In addition, the wave steepness of seven scenarios (9, 10, 11, 15, 23, 26, 31) is quite large (larger than 0.04) which is marked by “×” in the figure. Without considering these data, the nonlinear fitting curve is the blue line. The fitting formula can be expressed as,
h s w d = exp ( 1.583 0.028 T * )
The correlation coefficient is about 0.49 which was obviously improved. A declining trend can also be found from this formula.
Figure 8 shows the negative relationship between dimensionless sand wave length and dimensionless bed shear stress. The results rather show scatter when all the data are considered. The black line shows the fitting curve with linear relationship and the formula is as follows,
l s w d = 9.68 0.121 T *
The dimensionless sand wave length decreases with increasing bed shear stress. The correlation coefficient is about 0.47. When the eight tests are not considered, there is a better nonlinear relation (i.e., the blue line in Figure 8). The fitting equation is,
l s w d = exp ( 2.6 0.0427 T * )
This equation also shows the decline trend of the sand wave length with dimensionless bed shear stress. The correlation coefficient is improved to 0.87. In general, the dimensionless sand wave length decreases with increasing bed shear stress.
The steepness of sand waves is defined as σ s w = h s w / l s w . As Figure 9 shows, the linear fitting formula can be expressed as follows.
σ s w = 0.0156 + 0.000255 T *
Equation (26) shows that the steepness increases slightly with increase of the dimensionless shear stress. The data seems still rather scattered when the eight larger wave steepnesses of the wave conditions are excluded. However, the sand wave steepness seems to be a constant number which is about 0.02.
The dimensionless sand wave characteristics above are rather scattered, and here another method was used. The relationship of the dimensionless sand wave length denoted as l s w / L and the dimensionless shear stress can be found from Figure 10. The black line is the linear fitting line which is,
l s w L w = 0.315 + 0.0115 T *
The correlation coefficient is about 0.44. Excluding the eight tests, a good trend can be found as the blue line. The ratio of sand wave length and wave length is a constant with curve fitting, the value is about 0.44 which corresponds with Cataño-Lopera et al. [1].
Because of the complication of the calculation of dimensionless shear stress with Cnoidal wave theory in shallow water, the relationship between sand wave length and wave steepness is represented in Figure 11. Considering the large wave steepness, the fitting equation reads,
l s w g T 2 = 0.0032 + 2.227 H L
Without considering the large wave steepness, the fitting equation is (the blue one in Figure 11):
l s w g T 2 = 0.124 + 1.151 H L
The correlation coefficients are 0.758 and 0.767, respectively. The wave characteristics are easy to obtain, then the sand wave length can be evaluated directly.

4. Discussion

Tidal currents are the most critical environmental factors that control the formation and evolution of sand waves, as confirmed by many field observations and linear stability analyses. In the present study, we argue that waves and currents both have an important effect on sand waves.

4.1. Comparison with the Results of Cataño-Lopera et al. [1]

Cataño-Lopera et al. [1] chose the representative parameter Reynolds wave number (Rew) to express measured values. Twenty one data of tests with waves alone and combined currents were used to analyze the relationship between dimensionless geometric configuration of sand waves and Rew. The wave characteristics were obtained using small amplitude wave theory for all cases though Ursell number was larger than 26 in about half of the experiments. In this paper, Cnoidal wave theory was used to calculate the wave characteristics for the experiments with a large Ursell number. Additionally, two different sediment diameters were used in this paper, and the Reynolds wave number was no longer suitable as the dimensionless parameter because sediment configuration needs to be considered. Therefore, the dimensionless bed shear stress was used as the hydrodynamic parameter in this paper.
For the cases of Cataño-Lopera, most of the dimensionless shear stress was in the range of 7 to 20, with some data scattered in the range of 30 to 40. The results of our experiments supplement the range of shear stress data in the range from 20 to 30. For dimensionless height and length of sand waves (as Figure 7, Figure 8 and Figure 10 show), they are highly correlated when the dimensionless shear stress is less than 30. The results had poor correlation when the shear stress was larger than 30. According to the wave steepness, the cases larger than 0.04 represent a different trend. Without considering these cases, Equations (23), (25), (29) show good agreement with the measured data.
The measured data is scattered while the wave steepness is large, which may be caused by the wave nonlinearity. When the wave steepness is larger than 0.04, small amplitude wave theory is less accurate in calculating the orbital velocities and bed shear stress. Our equations gave a better relationship for cases when wave steepness was less than 0.04.
Sand wave steepness is 0.02 in our study which is in accordance with Cataño-Lopera et al. [1]. The ratio of sand wave length and wave length is about 0.44 which is also in accordance with Cataño-Lopera et al. [1].

4.2. Bedforms with Dimensionless Shear Stress

Sand waves and sand ripples coexist under the waves and currents. Dimensionless shear stress is a key parameter in this paper. According to van Rijn [38], the flow conditions in an alluvial channel are usually classified into: (1) Lower flow regime with plane beds, mini-ripples, mega-ripples, and dunes. (2) Transition flow regime with washed-out dunes and sand waves. (3) Upper flow regime with plane bed or anti-dunes. The upper regime which is defined as occurring for T * 25 , is characterized by a dominating suspended load transport. It was found from van Rijn [38] that the sand bed would be out of plane and start anti-dunes when the parameter was beyond 25. However, this concept has been challenged by Yalin [39] and Julien [40]. According to their work, lower-regime bedforms can be observed at values of T * well beyond 25 in very large rivers. Measurements from Table 2 include values of T * 25 . Sand waves can be found under the upper flow regime of van Rijn. Our experimental results confirmed the findings of Yalin [39] and Julien [40].
Wave energy dissipation and bottom roughness are important parameters in many hydrodynamic and morphodynamic models for coastal oceans, lakes and estuaries [41]. It is important to describe the topography of ripple pattern properly since it determines the net shear stress over the flow field. Xiong et al. [42] calculated the wave height of ripples based on the formulas of Boyd et al. [43] and Allen [44]. The scales are comparable to the results of Zhao et al. [45] at Dongtai. It was found that the ripple-related bottom shear stress is much larger than the skin-friction shear stress. In our experiments, the ripple pattern can be found all over the sand waves. It is reasonable to use Equation (8) to evaluate the roughness.

5. Conclusions

In this paper, a series of experiments was conducted to study sand bed configurations under waves and currents in the wave flume. Both the experimental results of this paper and previous research show that ripples and sand waves coexist together by waves alone and combined currents. Several conclusions can be drawn from the observations and analyses:
(1) Sand ripples can represent 2D patterns or 3D patterns depending on the location of the large-scale sand waves. The 2D ripples existing on the crest are much larger than those on the trough.
(2) Based on the Ursell numbers of different hydrodynamic conditions, the dimensionless bed shear stresses were calculated with different wave theories. Geometry characteristics of sand waves such as height, length, and steepness were very scattered for all cases. Without the cases of large wave steepness, the sand wave height and length both decreased with increase of the dimensionless bed shear stress. The sand wave steepness hardly changed with bed shear stress.
(3) A simple linear relationship can be found between the sand wave length and wave steepness. It is easy to evaluate the sand wave characteristics from the measured wave data.

Author Contributions

Conceptualization, B.L. and Z.W.; methodology, B.L.; formal analysis, Z.W. and B.L.; investigation, Z.W. and G.W.; writing—original draft preparation, Z.W.; writing—review and editing, Z.W. and G.W.; visualization, Z.W.

Funding

This research was funded by National Natural Science Foundation of China, (Grant No. 51679223, 51739010, 51709243).

Acknowledgments

The authors appreciate the comments from the anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

Notation

The following symbols are used in this paper:
H = wave height
T = wave period
u 0 = depth-averaged current velocity
d = water depth
D 50 = median diameter of sediment
U r = Ursell number
L = wave length
T * = dimensionless bed-shear stress
τ c w = mean bed shear stress due to current and waves
τ c r = critical bed shear stress
ρ s = density of sediment
ρ w = density of water
g = gravitational acceleration
θ c r = critical Shields number
D * = dimensionless particle diameter
ν = viscosity coefficient of water
u m = amplitude of the near-bottom wave orbital velocity
f w = wave friction factor
f c = current friction factor
τ w = wave-averaged bed shear stress due to waves only
τ c = bed shear stress due to current only
k = wave number
a m = near-bed peak orbital excursion
Δ = wave related roughness
η s = ripple height
κ = modulus of elliptic integral
K ( κ ) = the first kind complete elliptic integral
E ( κ ) = the second kind complete elliptic integral
u = near-bed instantaneous velocity due to waves
sn ( ) = Jacobian elliptic sinusoid function
cn ( ) = Jacobian elliptic cosine function
dn ( ) = Jacobian elliptic delta function
z t = distance from the wave trough to the sea bed
z c = distance from the wave crest to the sea bed
h s w = sand wave height
l s w = sand wave length
σ s w = sand wave steepness

References

  1. Cataño-Lopera, Y.A.; García, M.H. Geometry and migration characteristics of bedforms under waves and currents. Part 1: Sandwave morphodynamics. Coast. Eng. 2006, 53, 767–780. [Google Scholar] [CrossRef]
  2. Cataño-Lopera, Y.A.; García, M.H. Geometry and migration characteristics of bedforms under waves and currents: Part 2: Ripples superimposed on sandwaves. Coast. Eng. 2006, 53, 781–792. [Google Scholar] [CrossRef]
  3. Terwindt, J.H.J. Sand waves in the southern bight of the North Sea. Mar. Geol. 1971, 10, 51–67. [Google Scholar] [CrossRef]
  4. Mccave, I.N. Sand waves in the North Sea off the coast of Holland. Mar. Geol. 1971, 10, 199–225. [Google Scholar] [CrossRef]
  5. Campmans, G.H.P.; Roos, P.C.; de Vriend, H.J.; Hulshcher, S.J.M.H. Modeling the influence of storms on sand wave formation: A linear stability approach. Cont. Shelf Res. 2017, 137, 103–116. [Google Scholar] [CrossRef]
  6. Van Dijk, T.A.G.P.; Kleinhans, M.G. Processes controlling the dynamics of compound sand waves in the North Sea, Netherlands. J. Geophys. Res. Earth Surf. 2005, 110, 1–15. [Google Scholar] [CrossRef]
  7. Besio, G.; Blondeaux, P.; Brocchini, M.; Hulshcer, S.J.M.H.; Idier, D.; Knaapen, M.A.F.; Németh, A.; Roos, P.C.; Vittori, G. The morphodynamics of tidal sand waves: A model overview. Coast. Eng. 2017, 55, 657–670. [Google Scholar] [CrossRef]
  8. Berg, J.V.D.; Sterlini, F.; Hulscher, S.J.M.H.; Damme, R.V. Non-linear process based modelling of offshore sand waves. Cont. Shelf Res. 2012, 37, 26–35. [Google Scholar] [CrossRef]
  9. Knaapen, M.A.F.; Hulscher, S.J.M.H.; Tiessen, M.C.H.; Berg, J.V.D.; Parker, G.; García, M.H. Using a sand wave model for optimal monitoring of a navigation depth. In Proceedings of the Forth IAHR Symposium on River, Coastal and Estuarine Morphodynamics, Urbana, IL, USA, 4–7 October 2005; pp. 999–1007. [Google Scholar]
  10. Németh, A.; Hulscher, S.J.M.H.; de Vriend, H.J. Offshore sand wave dynamics, engineering problems and future solutions. Pipeline Gas J. 2003, 230, 67–69. [Google Scholar]
  11. Roos, P.C.; Blondeaux, P.; Hulscher, S.J.M.H.; Vittor, G. Linear evolution of sandwave packets. J. Geophys. Res. Earth Surf. 2005, 110, 1–10. [Google Scholar] [CrossRef]
  12. Ma, X.; Yan, J.; Fan, F. Morphology of submarine barchans and sediment transport in barchans fields off the Dongfang coast in Beibu Gulf. Geomorphology 2014, 213, 213–224. [Google Scholar] [CrossRef]
  13. Li, Y.; Lin, M.; Jiang, W.; Fan, F. Process control of the sand wave migration in Beibu Gulf of the South China Sea. J. Hydrodyn. 2011, 23, 439–446. [Google Scholar] [CrossRef]
  14. Dodd, N.; Blondeaux, P.; Calvete, D.; Swart, H.E.D.; Falqués, A.; Hulscher, S.J.M.H.; Różyński, G.; Vittori, G. Understanding coastal morphodynamics using stability methods. J. Coast. Res. 2003, 19, 849–865. [Google Scholar]
  15. Belderson, R.H. Offshore tidal and non-tidal sand ridges and sheets: Differences in morphology and hydrodynamic setting. In Shelf Sands and Sandstones—Memoir 11; CSPG Special Publications: Calgary, AB, Canada, 1986. [Google Scholar]
  16. Hulscher, S.J.M.H. Tidal-induced large-scale regular bed form patterns in a three-dimensional shallow water model. J. Geophys. Res. Oceans 1996, 101, 20727–20744. [Google Scholar] [CrossRef]
  17. Borsje, B.W.; Roos, P.C.; Kranenburg, W.M.; Hulscher, S.J.M.H. Modeling tidal sand wave formation in a numerical shallow water model: The role of turbulence formulation. Cont. Shelf Res. 2013, 60, 17–27. [Google Scholar] [CrossRef]
  18. Németh, A.A.; Hulscher, S.J.M.H.; de Vriend, H.J. Modelling sand wave migration in shallow shelf seas. Cont. Shelf Res. 2002, 22, 2795–2806. [Google Scholar] [CrossRef]
  19. Besio, G.; Blondeaux, P.; Brocchini, M.; Vittori, G. On the modeling of sand wave migration. J. Geophys. Res. Oceans 2004, 109, 1–13. [Google Scholar] [CrossRef]
  20. Campmans, G.H.P.; Roos, P.C.; de Vriend, H.J.; Hulscher, S.J.M.H. The influence of storms on sand wave evolution: A nonlinear idealized modeling approach. J. Geophys. Res. Earth Surf. 2018, 123, 2070–2086. [Google Scholar] [CrossRef]
  21. Campmans, G.H.P.; Roos, P.C.; Schrijen, E.; Hulscher, S.J.M.H. Modeling wave and wind climate effects on tidal sand wave dynamics: A North Sea case study. Estuar. Coast. Shelf Sci. 2018, 213, 137–147. [Google Scholar] [CrossRef]
  22. Borsje, B.W.; de Vries, M.B.; Bouma, T.J.; Besio, G.; Hulscher, S.J.M.H.; Herman, P.M.J. Modelling biogeomorphological influences for offshore sandwaves. Cont. Shelf Res. 2009, 29, 1289–1301. [Google Scholar] [CrossRef]
  23. Borsje, B.W.; Hulscher, S.J.M.H.; Herman, P.M.J.; de Vries, M.B. On the parameterization of biological influences on offshore sandwave dynamics. Ocean Dyn. 2009, 59, 659–670. [Google Scholar] [CrossRef]
  24. Roos, P.C.; Wemmenhove, R.; Hulscher, S.J.M.H.; Hoeijmakers, H.W.M.; Kruyt, N.P. Modeling the effect of nonuniform sediment on the dynamics of offshore tidal sandbanks. J. Geophys. Res. Earth Surf. 2007, 112, 1–11. [Google Scholar] [CrossRef]
  25. Van Oyen, T.; Blondeaux, P. Grain sorting effects on the formation of tidal sand waves. J. Fluid Mech. 2009, 629, 311–342. [Google Scholar] [CrossRef]
  26. Van Oyen, T.; Blondeaux, P. Tidal sand wave formation: Influence of graded suspended sediment transport. J. Geophys. Res. Oceans 2009, 114, 1–18. [Google Scholar] [CrossRef]
  27. Blondeaux, P. Sediment mixtures, coastal bedforms and grain sorting phenomena: An overview of the theoretical analyses. Adv. Water Resour. 2012, 48, 113–124. [Google Scholar] [CrossRef]
  28. Borsje, B.W.; Kranenburg, W.M.; Roos, P.C.; Matthieu, J.; Hulscher, S.J.M.H. The role of suspended load transport in the occurrence of tidal sand waves. J. Geophys. Res. Earth Surf. 2014, 119, 701–716. [Google Scholar] [CrossRef]
  29. Van Gerwen, W.; Borsje, B.W.; Damveld, J.H.; Hulscher, S.J.M.H. Modelling the effect of suspended load transport and tidal asymmetry on the equilibrium tidal sand wave height. Coast. Eng. 2018, 136, 56–64. [Google Scholar] [CrossRef]
  30. Chou, Y.J.; Shao, Y.C.; Sheng, Y.H.; Cheng, C.J. Stabilized formulation for modeling the erosion/deposition flux of sediment in circulation/CFD models. Water 2019, 11, 197. [Google Scholar] [CrossRef]
  31. Williams, J.J.; Bell, P.S.; Thorne, P.D. Unifying large and small wave-generated ripples. J. Geophys. Res. Oceans 2005, 110, 1–18. [Google Scholar] [CrossRef]
  32. Zhu, X.; Wang, Z.; Wu, Z.; Liu, F.; Liang, B.; He, Y.; Huang, F. Experiment research on geometry and evolution characteristics of sand wave bedforms generated by waves and currents. In 5th International Conference on Coastal and Ocean Engineering, Shanghai, China, 2018; IOP Conference Series: Earth and Environmental Science; IOP Publishing: Bristol, UK, 2018; pp. 9–15. [Google Scholar]
  33. Van Rijn, L.C. Principles of Sediment Transport in Rivers, Estuaries and Coastal Seas; Aqua Publications: Blokzijl, The Netherlands, 1993. [Google Scholar]
  34. Swart, D.H. Offshore Sediment Transport and Equilibrium Beach Profiles. Ph.D. Thesis, Civil Engineering and Geosciences, Delft, The Netherlands, December 1974. [Google Scholar]
  35. Korteweg, D.J.; de Vries, G. On the change of form of long waves advancing in a rectangular canal, and on a new type of stationary waves. Lond. Edinb. Dublin Philos. Magazine J. Sci. 1895, 39, 422–443. [Google Scholar] [CrossRef]
  36. Keller, J.B. The solitary wave and periodic waves in shallow water. Commun. Appl. Math. 1984, 1, 323–339. [Google Scholar] [CrossRef]
  37. Zou, Z.L. Coastal Dynamics; China Communications Press: Beijing, China, 2009. (In Chinese) [Google Scholar]
  38. Van Rijn, L.C. Sediment transport, part III: Bed forms and alluvial roughness. J. Hydraul. Eng. ASCE 1984, 110, 1733–1754. [Google Scholar] [CrossRef]
  39. Yalin, M.S. River Mechanics; Pergamon Press: Oxford, UK, 1992. [Google Scholar]
  40. Julien, P.Y. Erosion and Sedimentation; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  41. Wiberg, P.L.; Sherwood, C.R. Calculating wave-generated bottom orbital velocities from surface-wave parameters. Comput. Geosci. 2008, 34, 1243–1262. [Google Scholar] [CrossRef]
  42. Xiong, J.; Wang, Y.P.; Gao, S.; Du, J.; Yang, Y.; Tang, J.; Gao, J. On estimation of coastal wave parameters and wave-induced shear stresses. Limnol. Oceanogr. Meth. 2018, 16, 594–606. [Google Scholar] [CrossRef]
  43. Boyd, R.; Forbes, D.L.; Heffler, D.E. Time-sequence observations of wave-formed sand ripples on an ocean shoreface. Sedimentology 1988, 35, 449–464. [Google Scholar] [CrossRef]
  44. Allen, J.R.L. Physical Processes of Sedimentation; American Elsevier Pub. Co.: New York, NY, USA, 1970. [Google Scholar]
  45. Zhao, N.; Liu, X.J.; Lu, G. Characteristics and evolution process and their modeling of sand ripples on Dongtai tidal flats, Jiangsu province China. J. Sediment. Res. 2012, 3, 4, (In Chinese with English abstract). [Google Scholar]
Figure 1. Schematic diagram of model experimental layout.
Figure 1. Schematic diagram of model experimental layout.
Water 11 00612 g001
Figure 2. Flow chart to calculate the total bed shear stress.
Figure 2. Flow chart to calculate the total bed shear stress.
Water 11 00612 g002
Figure 3. Time series of the wave particle horizontal velocity with different wave theories. The red dash line represents the first test and the wave velocity is symmetrical in a wave period. The blue solid line represents the second test and the wave velocity is larger following the wave direction but the duration is shorter.
Figure 3. Time series of the wave particle horizontal velocity with different wave theories. The red dash line represents the first test and the wave velocity is symmetrical in a wave period. The blue solid line represents the second test and the wave velocity is larger following the wave direction but the duration is shorter.
Water 11 00612 g003
Figure 4. Bed evolution over time during the experiment. (a) initial time; (b) several minutes later; (c) 1 h later; (d) 10 h later.
Figure 4. Bed evolution over time during the experiment. (a) initial time; (b) several minutes later; (c) 1 h later; (d) 10 h later.
Water 11 00612 g004
Figure 5. Bottom evolution over time: the generation of a sand wave.
Figure 5. Bottom evolution over time: the generation of a sand wave.
Water 11 00612 g005
Figure 6. The growth rate of a sand wave as a function of time.
Figure 6. The growth rate of a sand wave as a function of time.
Water 11 00612 g006
Figure 7. The relationship between the measured dimensionless sand wave height and dimensionless shear stress. The blue circles represent the results of experiments of this paper with a medium diameter of 0.17 mm and the dark ones of 0.18 mm. The red circles represent the results of experiments conducted by Cataño-Lopera et al. [1]. The range of the dimensionless shear stress is about 7–49.
Figure 7. The relationship between the measured dimensionless sand wave height and dimensionless shear stress. The blue circles represent the results of experiments of this paper with a medium diameter of 0.17 mm and the dark ones of 0.18 mm. The red circles represent the results of experiments conducted by Cataño-Lopera et al. [1]. The range of the dimensionless shear stress is about 7–49.
Water 11 00612 g007
Figure 8. The relationship between the measured dimensionless sand wave length and dimensionless shear stress.
Figure 8. The relationship between the measured dimensionless sand wave length and dimensionless shear stress.
Water 11 00612 g008
Figure 9. The relationship between the measured sand wave steepness and dimensionless shear stress. The black line is the fitting curve for all the data.
Figure 9. The relationship between the measured sand wave steepness and dimensionless shear stress. The black line is the fitting curve for all the data.
Water 11 00612 g009
Figure 10. The relationship between the measured dimensionless sand wave length and dimensionless shear stress.
Figure 10. The relationship between the measured dimensionless sand wave length and dimensionless shear stress.
Water 11 00612 g010
Figure 11. The relationship between the measured dimensionless sand wave height and dimensionless wave length.
Figure 11. The relationship between the measured dimensionless sand wave height and dimensionless wave length.
Water 11 00612 g011
Table 1. Experimental scenarios.
Table 1. Experimental scenarios.
ScenariosH (m)T (s)u0 (m/s)d (m)D50 (mm)
10.10200.60.17
20.12200.60.17
30.14200.60.17
40.16200.60.17
50.142−0.20.60.17
60.1420.20.60.17
70.1420.250.60.17
80.1420.350.60.17
90.20200.60.18
100.22200.60.18
Where H is wave height; T is wave period; u0 is depth-averaged current velocity; d is the water depth, and D50 is the median diameter of sediment.
Table 2. Hydrodynamic conditions and characteristics of measured sand waves.
Table 2. Hydrodynamic conditions and characteristics of measured sand waves.
Scenarios aH
(cm)
T
(s)
u0
(m/s)
d
(m)
D50
(mm)
hsw
(cm)
lsw
(m)
Ur τ w τ c τ c r T *
110.0200.60.171.092.288.813.490.000.1522.10
212.0200.60.172.792.210.574.190.000.1526.72
314.0200.60.173.181.9612.334.920.000.1531.49
416.0200.60.173.791.7814.095.660.000.1536.41
514.02−0.20.60.172.651.8412.334.92−0.270.1529.71
614.020.20.60.173.672.0412.334.920.270.1533.26
714.020.250.60.175.642.4112.334.920.420.1534.26
814.020.350.60.173.812.4112.334.920.820.1536.93
920.0200.60.186.23.917.627.220.000.1644.05
1022.0200.60.1812.54.5319.388.030.000.1649.10
1119.6200.560.2526.95.820.127.470.000.2232.57
12 b12.85.200.560.253.35.1117.191.930.000.227.65
1319.26.900.560.2511.96.8347.561.800.000.227.07
1419.74.200.560.259.83.4123.083.730.000.2215.77
1517.71.600.560.2521.510.327.080.000.2230.79
1620.92.900.560.255.92.656.766.070.000.2226.27
1710.73.400.560.257.73.437.882.630.000.2210.81
1817.42.300.560.2539.98.424.836.160.000.2226.66
1915.15.900.560.2511.55.1186.261.860.000.227.37
2010.44.100.560.2564.255.462.120.000.228.51
2114.74.100.560.258.64.182.122.950.000.2212.24
2212.35.100.560.257.45.1107.381.910.000.227.57
2319.61.90.170.560.2513.55.117.757.600.200.2234.04
2419.25.80.170.560.2512.46.8240.362.320.200.2210.34
2519.73.70.170.560.258.24.192.734.380.200.2219.58
2617.71.50.170.560.2553.48.657.100.200.2231.77
2711.33.60.170.560.253.54.145.682.630.200.2211.71
2810.720.170.560.257.52.110.983.930.200.2217.57
2912.15.10.170.560.2514.56.3105.381.880.200.228.33
30134.40.170.560.2511.65.883.152.410.200.2210.74
312320.170.560.2512.52.323.618.950.200.2240.08
a The first 10 scenarios were conducted in this paper, the last 21 scenarios were conducted by Cataño-Lopera et al. [1]. b The rows with brown color denote the shallow water, with Cnoidal wave theory.

Share and Cite

MDPI and ACS Style

Wang, Z.; Liang, B.; Wu, G. Experimental Investigation on Characteristics of Sand Waves with Fine Sand under Waves and Currents. Water 2019, 11, 612. https://doi.org/10.3390/w11030612

AMA Style

Wang Z, Liang B, Wu G. Experimental Investigation on Characteristics of Sand Waves with Fine Sand under Waves and Currents. Water. 2019; 11(3):612. https://doi.org/10.3390/w11030612

Chicago/Turabian Style

Wang, Zhenlu, Bingchen Liang, and Guoxiang Wu. 2019. "Experimental Investigation on Characteristics of Sand Waves with Fine Sand under Waves and Currents" Water 11, no. 3: 612. https://doi.org/10.3390/w11030612

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop