Next Article in Journal
The Suitability of Snow and Meteorological Conditions of South-Central Slovakia for Ski Slope Operation at Low Elevation—A Case Study of the Košútka Ski Centre
Next Article in Special Issue
Vulnerability Analysis of the Venetian Littoral and Adopted Mitigation Strategy
Previous Article in Journal
Analyzing the Effect of Soil Hydraulic Conductivity Anisotropy on Slope Stability Using a Coupled Hydromechanical Framework
Previous Article in Special Issue
Simulating Future Flows and Salinity Intrusion Using Combined One- and Two-Dimensional Hydrodynamic Modelling—The Case of Hau River, Vietnamese Mekong Delta
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Applying a New Force–Velocity Synchronizing Algorithm to Derive Drag Coefficients of Rigid Vegetation in Oscillatory Flows

1
School of Marine Science, Sun Yat-sen University, Guangzhou 510275, China
2
State Key Laboratory of Hydrology Water Resources and Hydraulic Engineering, Nanjing Hydraulic Research Institute, Nanjing 210029, China
3
Guangdong Province Engineering Research Center for Coasts, Islands and Reefs, Guangzhou 510275, China
4
Guangdong Research Institute of Water Resources and Hydropower, Guangzhou 510630, China
5
School of Marine Engineering and Technology, Sun Yat-sen University, Guangzhou 510275, China
*
Authors to whom correspondence should be addressed.
Water 2018, 10(7), 906; https://doi.org/10.3390/w10070906
Submission received: 24 May 2018 / Revised: 5 July 2018 / Accepted: 6 July 2018 / Published: 9 July 2018

Abstract

:
Coastal vegetation is effective in dissipating incident wave energy during storm conditions, which offers valuable protection to coastal communities. Determining vegetation drag coefficient (CD) is of great importance to the quantification of vegetation-induced wave dissipation. Recently, a direct measuring approach has been developed to derive vegetation drag coefficient more accurately compared to the conventional calibration approach. However, as this approach requires perfectly in-phase force and velocity signals, there are two difficulties associated with it. The first difficulty is the availability of a suitable force sensor to compose synchronized force–velocity measuring systems. The second difficulty is related to realigning the obtained timeseries of force and velocity data. This technical note develops a new synchronized force–velocity measuring system by using standard force sensors and an acoustic doppler velocimeter (ADV). This system is applied together with an automatic realignment algorithm to ensure in-phase data for CD deviation. The algorithm reduces the phase shift between force–velocity signals from ca. 0.26 s to 0.003 s. Both time-varying and period-averaged CD can be obtained using this method. The derived CD can be used to accurately reproduce the measured maximum total acting force on vegetation (R2 = 0.759), which shows the reliability of the automatic alignment algorithm. The newly-developed synchronized force–velocity measuring system and alignment algorithm are expected to be useful in future experiments on vegetation–wave interactions with various hydrodynamic and vegetation settings.

1. Introduction

Mangroves, saltmarshes, and seagrasses are important coastal ecosystems that are widely distributed in world’s coasts [1,2,3]. The wave-damping capacity of these coastal wetlands has been increasingly recognized [1,4,5,6,7,8]. These coastal wetlands can significantly reduce wave energy even under storm or tsunami conditions [9,10], which provides valuable protection to the coastal communities and properties [11,12]. Over the past decades, the wave heights have a clear increase trend in extreme conditions, together with accelerated sea level rise [13,14,15,16,17]. Therefore, there is a demand for better understanding and predictive ability of vegetation–wave interaction process to reduce coastal flooding risks [4,18,19,20].
The main impact of vegetation on incident waves is exerting an additional force on water motion [21,22]. This force can be described by the Morison equation, which is composed by drag force (FD) and inertia force (FM) [23]. For normal field conditions, the drag force is the dominant force, and most relevant for wave energy dissipation. In the Morison equation, FD is proportional to the square of impact velocity on vegetation stems. When the velocity scale is determined, the magnitude of FD varies linearly with vegetation drag coefficients (CD). In oscillatory (wavy) flows, the CD values have a large range of variations (i.e., 0.1 to 100) [24]. The CD values depend on canopy density, hydrodynamic conditions, as well as the morphology of the individual canopy elements. Thus, choosing appropriate CD values are important for accurate simulation of FD, and the resultant wave dampening in many modelling studies [18,25,26,27,28,29,30].
Currently, there are two methods available in determining CD: the calibration method and the direct measurement method. The calibration method is a convectional method developed in the 1990s [31,32], and has been widely used since [33,34,35]. It derives CD by calibrating its values to obtain the best fit between modelled and measured wave height evolution over vegetation fields. The direct measurement method is a new method, which has been developed since the 2010s [36,37,38]. This method directly applies the Morrison equation and measured in-phase force and velocity data to determine CD. The main differences between the calibration and the direct measurement method are: (1) the calibration method can only provide period-averaged CD, but the direct measurement method can derive both period-averaged and time-varying CD; (2) the direct measurement method can eliminate the potential errors often associated with the calibration method, and lead to CD–Re (Reynolds number) relations with better fits, which are desirable for model applications [37].
Since the direct measurement method relies on in-phase force and velocity data, there are two difficulties when applying this method. The first difficulty is the availability of suitable force sensors to assemble synchronized force–velocity measuring systems. The force sensors should be waterproofed and durable in wave flumes, where frequent water logging and splashing occur. Additionally, the sensors should be small enough to fit into wave flumes. The second difficulty associated with this method is the data processing technique required to obtain perfectly aligned force–velocity data for CD derivation [37]. The data alignment is critical for the direct measurement method, as there are time lags between original force data and velocity data signals (ca. 0.2 s), which may lead to large errors in the derived CD. These time lags may originate from small misalignments between force sensors and velocity measurement [37]. They may also be induced by intrinsic time shifts in instrument recordings. The maximum wave energy dissipation occurs at the peak wave orbital velocity in phase with the peak drag force when CD values matter the most. Thus, in order to obtain accurate CD values, it is important to minimize the time lags. Previous studies firstly set an intrinsic time lag between the force and velocity data, and then started iterations to reduce the time lag. Note that this intrinsic time lag varies with different instrument set-ups. This intrinsic time lag needs to be carefully tuned to obtain in-phase data. It is, however, preferable to have an automatic algorithm that can provide generic solutions to the alignment problem.
In this technical note, we have developed (1) a synchronized force–velocity measuring system by using standardized force sensors that can easily fit in wave flumes; (2) and an automatic alignment algorithm to obtain in-phase force–velocity data for CD derivation. The new force–velocity measuring systems are applied in a flume at four locations in a mimicked mangrove canopy, which were tested with various simulated wave conditions. The automatic alignment algorithm was then applied to reduce the time lags between force–velocity signals. The processed data were subsequently used to derive both time-varying and period-averaged CD. To evaluate the accuracy of the derived CD (and also the alignment algorithm), we used the derived CD to reproduce the total acting force on mimicked vegetation, and compared it with the measurements.
The rest of the technical note is organized as follows: Section 2 introduces the automatic alignment algorithm and the set-up of the synchronized force–velocity measuring system in the wave flume at the Fluid Mechanics Laboratory at Sun Yat-sen University, Zhuhai Campus. Section 3 demonstrates the original and processed force–velocity data, as well as the reproduced force data results. Section 4 discusses the current limitations and provides an outlook for future applications of this method. Finally, Section 5 provides conclusions of the current note.

2. Materials and Methods

2.1. The Direct Measuring Method for CD Derivation

The force acting on a single stem can be expressed by Morison equation [23] as
F = F D + F M = 1 2 ρ C D h v b v U | U | + π 4 ρ C M h v b v 2 U t
where F is the total acting force on vegetation which can be obtained by force measurement. FD is drag force, and FM is inertia force. ρ is the density of the fluid. CD and CM are the drag and inertia coefficients, respectively. hv is the height of vegetation in water, and bv is the diameter of circular cylinder. U is the depth-averaged flow velocity. When linear wave theory is applied, the U varies as a function of sine:
U = U w sin ( ω t )
where Uw is the amplitude of horizontal wave orbital velocity. Following linear wave theory, Uw can be expressed as
U w = π H T cos h [ k ( h + Z 0 ) ] sin h ( k h )
where H is wave height, T is wave period, k is wave number, h is water depth, and Z0 is the vertical position of the considered point, which is 0 at the still wave level, and −h at the sea bed. Equation (3) was used to estimate Uw when the velocity measurement is unavailable. CM is often assumed to be equal to 2 for cylinders (e.g., [39]). To derive the time-varying CD, we can apply the following equation:
C D = 2 F D ρ h v b v U | U | = 2 ( F F M ) ρ h v b v U | U |
where FM can be derived based on U t using the timeseries of velocity data, and other parameters (i.e., ρ , bv, hv, CM) in FM are known. Thus, time-varying CD can be obtained readily when in-phase force and velocity data is obtained.
Period-averaged CD is relevant to vegetation-induced wave dissipation. It was not computed as the temporal mean of the time-varying CD. Time-varying CD has great variability over one wave period [37]. Specifically, its value is infinite when the velocity is close to zero. However, those CD values are not relevant for vegetation-induced wave energy dissipation, as dissipation is highest at the velocity peaks. Thus, CD values at high velocity matter the most. To obtain relevant period-averaged CD values, the direct measurement method applies the technique of quantifying the power and work done by the acting force ( ε ) [37]. The time-varying power of FD and FM is evaluated as follows:
P D = F D U
P M = F M U
The work done by the total acting force (F) over a wave period (T) is
W = 0 T F U d t
If we substitute F with Equation (1), then we obtain
W = W D + W M = 1 T 0 T F D U d t + 1 T 0 T F M U d t = 1 2 T 0 T ρ C D h v b v U 2 | U | d t + π 4 T 0 T ρ C M h v b v 2 U t U d t
WD and WM is the work done by FD and FM over a full period, respectively. As U is a sine function (Equation (2)), the work done by F M over a full wave period (i.e., second term on the right) is zero. Thus, the work done by F is equal to the work done by F D :
W = 0 T F U d t = 1 2 0 T ρ C D h v b v U 2 | U | d t
Finally, the period-averaged CD can be derived based on the above equation:
C D = 2 0 T F D U d t 0 T ρ h v b v U 2 | U | d t = 2 0 T F U d t 0 T ρ h v b v U 2 | U | d t
The in-phase time series data of total force (F) and velocity (U) can be used directly in Equation (10) to drive period-averaged CD values. As W is proportional to U 3 (Equation (9)), the integration of FU over a period is largely contributed to by the moments with relatively high velocity, and to a very limited extent, by the moments with low velocity. Thus, deriving period-averaged CD via the technique of quantifying ε can automatically assign large weight to the moments with high velocities in a wave period, resulting in most relevant CD values for wave dissipation analysis. To check the validity of the direct measuring method, we used the derived period-averaged CD values to reproduce the total acting forcing (Frep) using Equation (1), and compare it with the actual measurement. Additionally, another reproduced total force Frep is included by assuming CD = 1. It is used as a reference for the Frep.

2.2. Synchronized Force–Velocity Measuring System

In-phase force–velocity data are critical to the direct measurement method. To obtain in-phase data, a synchronized force–velocity measuring system was developed, which was composed by a force sensor and an acoustic doppler velocimeter (ADV) (Figure 1a,d). Four measuring systems were deployed in the wave flume at Fluid Mechanics Laboratory at Sun Yat-sen University, Zhuhai Campus (Figure 1a). The wave flume is 20 m long, 0.8 m wide, and 0.6 m deep. A series of capacitance-type wave gauges were installed to monitor wave height changes in the wave flume. The mimicked vegetation canopy was 8 m long, and it was constructed by PVC pipes. The pipes were 0.2 m tall and their diameter was 0.02 m. The mimicked vegetation canopy was built following a stagger pattern with a density of 139 stems/m2. It was built on top of a false bottom in order to elevate the canopy so that the force sensors can be mounted underneath it. The mimicked vegetation canopy was submerged in water (water depth = 0.25 m), and it was subjected to various wave conditions. The tested wave height varied from 0.03 to 0.09 m, and the tested wave period varied from 0.6 to 1.2 s. A space-averaged CD can be obtained by taking the mean CD values of the four measuring spots in each test.
The force sensors selected were model M140 built by Utilcell, Spain (Figure 1b). As the size of this sensors is small, they can be installed at multiple locations in one flume test. The reading of the sensor is in “gram”, which can be easily translated into “Newton” by multiplying the acceleration of gravity. The minimum division of the sensor is 3 × 10 3 N, and the maximum measuring load is 30 N. In the current experiment, the measuring frequency of the force sensors is set as 20 Hz. These sensors were chosen also because they are robust and can be easily waterproofed by sealing the cable connection point with glue. Furthermore, the sensor is small (15 cm × 4 cm × 2.5 cm), and can be easily fitted in the flume (Figure 1b). To prevent the sensor being affected by any force acting on itself, an aluminum case was put around the sensor.
For each force sensor, a mimicked vegetation stem (PVC pipe) was firmly screwed to it, so that the force acting on a vegetation stem can be detected (Figure 1b). The PVC pipes that were attached to the force sensor were not different from other pipes in the mimicked vegetation canopy. After the sensor was attached to a PVC pipe, it measured the total acting force on the pipe. For the same force, the reading is constant, regardless of the location of the acting force. This is desirable for our current experiment, in which wave-induced forces acted over the full length of the pipe. Additionally, the sensors can detect the force in both the following and the opposing direction as the wave propagation, which is ideal to measure the force generated by oscillatory flows.
As a first test of the force sensors, we put known weights on one of the sensors attached with PVC pipes. This sensor was held horizontally in this test. The obtained readings were subsequently compared to the known weights. Different weights were put at three different positions on the pipe, i.e., bottom, middle, and tip. The generated forces were in both positive and negative directions. The comparison between the readings and the weights is listed in Table 1. It is clear that for the force measurement in both directions, the relative errors of the sensors are within 1%, compared to the known weights. This first test shows that the force sensors can provide precise measurements.
For the velocity measurement, we used four ADVs. They were deployed at the same cross-section as the force sensors to obtain roughly synchronized signals. The cases wave0312 (3 cm wave height and 1.2 s wave period) and wave0709 (7 cm wave height and 0.9 s wave period) were selected to conduct velocity profile measurements in order to obtain information on velocity structure and depth-averaged in-canopy velocity. The profile was obtained by manually adjusting the vertical measuring location of the ADVs in repeated sequences. For most cases, the velocity measurement was taken at half of the water depth, which was a representative value of the depth-averaged in-canopy velocity. The accuracy of this treatment is acceptable when the submergence ratio (i.e., h/hv = 1.25) is small [24,37], and it can significantly reduce the labor involved. Two of the ADVs were made by Nortek (Vectrino see http://www.nortekusa.com/usa/products/acoustic-doppler-velocimeters/vectrino-1), and the other two were made by SonTek (MicroADV, see https://www.sontek.com/argonaut-adv). These four ADVs are common instruments in fluid mechanics labs. Their basis measurement technology is coherent Doppler processing. They measure 3D water velocity of a small cylinder (i.e., within 1 cm3) that is a few centimeters away from measuring probes in the water. They can measure at frequencies as high as 64 Hz, which are desirable for the direct measuring method. In our experiment, the ADV data acquisition followed their respective user manuals. The measuring frequency was set as 40 Hz in order to accommodate the measuring frequency of the force sensor (i.e., 20 Hz). The obtained data are filtered through a low-pass filter to remove high frequency spikes following a similar method described in Strom and Papanicolaou [40].

2.3. Automatic Alignment Algorithm

Although the force and velocity measurements were deployed at the same cross-section of the wave flume, it did not ensure perfectly in-phase data. In fact, small time lags commonly existed between the obtained original force and velocity time series. These time lags were induced by small misalignments between force and velocity measurement, and/or by intrinsic delays of these electronic devices. To reduce the time-lag between velocity and force measurement, the original force and velocity time series should be realigned.
According to the Morison equation [23], velocity (U) and drag force (FD) should be in phase, which can be used to evaluate the time-lag between velocity and force measurement. A flow chart for the data realignment is shown in Figure 2. The inputs are the timeseries of force (F) and velocity (U). As we can assume that CM = 2 [39], the inertia force can be calculated based on the velocity. Then, the drag force (FD) can be computed by subtracting inertia force (FM) from the total force (F). Subsequently, we can determine the phase shift (∆t) between the velocity and drag force peaks. Lastly, this phase shift (∆t) will be recorded and used to adjust the velocity timeseries, aiming to obtain more in-phase velocity and force data. The obtained new velocity and force data will be used as input in the same loop. This loop continues 30 times, and we chose the minimum phase shift (∆t) and the resultant velocity and force timeseries as outputs for CD derivation. The automatic alignment algorithm is provided in the Appendix A as a MATLAB script. To verify the 30 loop count criterion, a sensitivity analysis is conducted by changing the loop count to 10, 20, 30, and 50. The resulting phase shifts with those loop counts are subsequently compared.

3. Results

3.1. Velocity Profiles in the Vegetation Canopy

Figure 3 shows the velocity profiles were measured at location 3 in case wave0312 (3 cm wave height and 1.2 s wave period) and wave0709 (7 cm wave height and 0.9 s wave period). For the rest of the tested cases, velocity was measured at the half water depth as a proxy of depth-averaged in-canopy velocity. Figure 3a shows that in the case wave0312, the velocity profiles are rather uniform in the vertical direction. The depth-averaged in-canopy velocity amplitude is 0.052 m/s, whereas Uw measured at the half water depth is 0.049 m/s. The difference between these two is small. The velocity profiles in wave0709 have greater vertical gradient, i.e., higher velocity at the top and lower velocity near the bottom. Overall, the difference between Uw measured at the half water depth (0.106 m/s) and the amplitude of depth-averaged in-canopy velocity (0.115 m/s) is small. Therefore, it is acceptable to use Uw measured at the half water depth as a representative value of the depth-averaged in-canopy velocity.

3.2. Wave Height and Wave Orbital Velocity in the Mimicked Vegetation Canopy

Reductions of wave height (H) and magnitude of wave orbital velocity (Uw) through mimicked vegetation canopy can be observed in Figure 4. The wave height reduces continuously from the canopy front to the end. The final wave height reduction rate was 55% (Figure 4a). The shown wave orbital velocity is obtained by ADV measurement at location 1–3. The ADV measurement at location 4 failed during the experiment. The shown Uw is obtained by using Equation (3) based on an average wave height between x = 6–8 m in Figure 4a. With the reduced wave height, the magnitude of wave orbital velocity also reduces from 0.155 m/s to 0.095 m/s from the beginning to the end of the canopy. The reduced wave orbital velocity leads to variations in acting force on vegetation stem as well as in vegetation drag coefficient (CD), which are shown in the following sections.

3.3. Data Alignment and Time-Varying CD

The total acting force (F) and velocity measured at four locations in the wave flume are shown in Figure 5. Waves reach these four locations at different moments. The velocity data at the last measurement point were missing due to the failure of the ADV measurement. It is clear that the acting force (F) and velocity reduce as waves pass through mimicked vegetation canopy. It is also apparent that there are time lags in synchronized force–velocity measurements at all locations. However, it should be noted that these seeming time lags (i.e., time shifts between F and U) are not the real time lags (i.e., time shifts between FD and U), which lead to errors in CD.
Following the previous study [37], we only tracked the first 2–3 full wave periods after the start-up, but before waves reached the back end of the flume, to avoid possible influence of wave reflection. Figure 5e shows the force and velocity data of the chosen first two wave periods after the start-up. Based on the original data and Equation (1), both FD and FM can be estimated as shown in Figure 5e. However, it should be noted that as the total force and velocity data are not in phase. The derived FD and FM are the first estimation. The derived peaks of the FD are even higher than the total force, which is not possible. This result highlights the necessity of obtaining in-phase data. Judging from the peaks between U and estimated FD peaks, the time lag between those two signals is 0.26 s.
In order to eliminate the time lags between these signals, the realignment algorithm was applied to obtain synchronized velocity and force data (Figure 6a–d). It is clear that after the realignment procedure, the time lag between U and FD is largely reduced (Figure 6a,b), but it cannot be completely eliminated, as small shifts in signal peaks and troughs still exist. Based on the phase shifts of the two peaks and two troughs (as indicated by the red double arrow lines), the time shift of the realigned U and FD is significantly reduced to 0.003 s, which is only 1% of the original time shift before the realignment. Note that this optimized time shift is the mean shift of the tested two wave periods (i.e., shifts of two peaks and two troughs). Apart from the reduced time shift, the magnitude of FD is also reduced to be lower than the total force (F), which is in line with the original Morison equation. It is noted that the peak drag force before alignment is about twice as large as the peak drag force after the alignment. Thus, if the original F timeseries were used, the derived drag coefficient would be considerably overestimated.
Based on the synchronized velocity and force data and Equation (4), the time-varying CD can be derived (Figure 6e). The time-varying CD values vary periodically with the changing velocity signals. The values are small when the velocity is large, but they reach to infinity when the velocity is close to zero. The reason for the unrealistically large CD values is because the FD are divided by very small velocity values in Equation (4). These unrealistically large CD values are not useful in modelling wave dissipation by vegetation, since they are associated with period with very low velocity when the energy dissipation is very limited.
To verify the loop count criterion applied in the realignment algorithm, a sensitivity analysis of the loop is conducted (Figure 7). It is clear that with the increase of loops, the phase shifts can be sufficiently reduced. When 30 loop count is applied, the phase shifts can be reduced to minimum, i.e., 0.003 s for the case 1 and almost zero for the case 2. If the number of loops is increased to 50, the time shifts cannot be further reduced. Thus, the results of sensitivity analysis indicate that applying 30 loop count in the realignment algorithm is a reasonable criterion.

3.4. Deriving Period-Averaged CD

Based on the realigned force–velocity data, we can also derive period-averaged CD by quantifying the power and work done, FD and FM. It is clear that the time-varying power of FD is always positive, and its magnitude varies in phase with the velocity magnitude (Figure 8). For ideal sinusoidal velocity signals, PD at the velocity peaks should be equal to the troughs. However, in our test case (and in real field conditions), the wave orbital velocity is asymmetrical: higher in the positive direction and lower in the negative direction. Thus, PD is larger near the wave peaks and smaller near the wave trough. The difference between peaks and troughs are much more apparent in PD compared to the difference in velocity. It is because that PD is to the third power of velocity. Small asymmetry in velocity will be greatly magnified in PD. The variation of PM is different from that of the PD. It is clear that PM varies between positive and negative values. The zero-crossings in PM occur when velocity is zero or when velocity is at its maximum in both directions, i.e., when FM is zero in Equation (5). Note that there are small fluctuations at the peaks of the PD. These fluctuations may be induced by the small phase shifts between the FD and U timeseries.
By integrating PD and PW over one wave period, we can obtain the work done by drag force (WD) and inertia force (WM) as shown in Equation (8). The averaged WD and WM over the two periods in Figure 8 is 2.0   ×   10 3   J and 2.0   ×   10 4   J , respectively. In case of ideal sinusoidal velocity signals, WM should be exactly zero. Due to the asymmetric wave velocity, the WM is not zero, but the magnitude of WM is fairly small, i.e., one tenth of the WD. Since the magnitude of WM is considerably small compared to WD, the assumption that WM can be ignored in the deriving period-averaged CD (in Equation (8)) is still valid, and the period-averaged CD in the shown case is derived as 3.25.
Following the same method, period-averaged CD at three functional measuring locations of all the tested cases are listed in Table 2. It shows that relatively large deviations in CD values exist among different measuring locations. The CD values generally increase from locations from 1 to 3. Previous studies have shown that CD values increase with the reduced velocity (i.e., Reynolds number) [24,37]. The obtained increase of CD values may be related to the reduction of wave orbital velocity from the front to the end of the vegetation canopy, as shown in Figure 4b Thus, the spatial variation in CD values is in-line with previous studies. Additionally, it is noted that the cases with larger wave height and wave period (i.e., higher wave orbital velocity) generally lead to smaller spatially averaged CD, which is also in agreement with previous studies [24,37].

3.5. Assessing the Derived CD by Reproducing Acting Force

In order to test the derived period-averaged CD, we used the derived values to reproduce the total force from the velocity signals using Equation (1). The reproduced total force (Frep) is subsequently compared with the measured actual total force (Figure 9). The reproduced total force (Frep) using CD = 1 is also included as reference. It is clear that Frep is in good agreement with the measured force over the shown two wave periods, although small differences exist between them. Notably, the measured maximum force is well captured in Frep near x = 0.5 π, which is important as the maximum force is critical not only for energy dissipation but also for assessing the stem strength to wave loading. As a comparison, the difference between Frep and the measurement is large, which shows the validity of using period-averaged CD to reproduce the total force.
Figure 10 compares the maximum Frep and maximum measured total force obtained at all three functional measuring locations in all test cases. In general, the reproduced maximum force is well in-line with the measurement, as most of the data points are fairly close to the 1:1 reference line. The R2 value is 0.759 for data of all three functional measuring locations in all test cases. This result indicates the CD deriving procedure is valid, and the intrinsic errors associated with this procedure are acceptable.

4. Discussion

4.1. Advantages of the Current Measuring System and Alignment Algorithm

Our results have shown that large spatial variations can exist in the wave particle velocity and CD (Figure 5 and Table 2). Thus, it is important to have synchronized force–velocity measurement at multiple locations by a number of force–velocity measuring systems. The selected standard force sensors are small enough to be installed at multiple locations in wave flumes. Additionally, these sensors are designed with built-in tapped holes, which facilitate testing various vegetation mimics, e.g., rigid cylinders, flexible stripes, and real vegetation stems.
Our results further show that the realignment process is important to derive both time-varying and period-averaged CD values (Figure 6 and Figure 7). In our experiment, we aligned the instruments as good as possible (please see Figure 1d of the manuscript), but it is inevitable to have small misalignment to cause the delay. The main source causing the delay may be the inherent difference in instruments’ speed of recording and receiving data, as the force and velocity measurements have their separate data acquisition systems. Thus, the automatic synchronizing algorithm is necessary and valuable in the current study. By using this algorithm, the time shift between two signals can be reduced to 0.003 s, which is only 1% of the original time shift before the realignment. The obtained time shifts are believed to be acceptable when comparing to normal wave periods (1–2 s) tested in our lab flume. The time shifts are merely 0.3% to the tested wave period.
The overall good performance between maximum measured force and reproduced force shows the reliability of this alignment algorithm. Importantly, this algorithm can automatically process the force and velocity data. No manual tuning is needed. Hence, it provides a generic solution to the alignment problems in deriving CD. Furthermore, this algorithm can run very efficiently, which is desirable when processing large data sets from multiple measuring locations. Lastly, this alignment algorithm is applied to process the velocity data from ADVs, but it is worth noticing that this algorithm is also applicable for other velocity measuring technologies, e.g., EMF (electromagnetic flow manufacture meter) and PIV (particle image velocimetry) [41,42,43].

4.2. Current Limitations and Future Applications

It is noted that the time shifts after the realignment are non-zero, but they are significantly reduced. To further reduce the time shifts, high frequency force and velocity measurements (e.g., >100 Hz) are required to obtain finite time steps for the realignment algorithm. However, it is perhaps not possible to completely eliminate the time shifts for all the tested cases, especially when multiple wave periods are included in the analysis, as the realignment procedure needs to account for time shifts at multiple peaks.
The velocity measurement in the current experiment was conducted by ADV measurement, which is a conventional method in flume experiments. The main limitation of the ADV measurement is that it is a point measurement. To obtain vertical velocity profile, it is required to manually adjust the ADV-measuring locations and repeat the same test conditions for each measuring location. This process is very time-consuming. Thus, we only conducted the velocity profile measurement for two cases, whereas for other cases, the velocity data is taken at the half water depth, which roughly equaled to the mean in-canopy velocity (Figure 3). The same practice is also done in previous study [37]. However, it is possible that the small deviation between the point velocity and depth-averaged in-canopy velocity can lead to errors in the derived CD values. This may partly explain the difference between the maximum measured force and the reproduced force. In order to improve the velocity measuring accuracy and reduce the labor involved, PIV system can be applied in future experiments. The PIV system can provide detailed velocity information of velocity field [41,42]. By applying such a system, it is also possible to obtain the relative velocity between water motion and the motion of flexible vegetation stems. Thus, the developed technics in the current study can be further applied in flexible vegetation canopies, e.g., saltmarshes and seagrasses, which is interesting to both coastal engineers and ecologists.

5. Conclusions

This technical note provides a practical set-up to derive both time-varying and period-averaged vegetation drag coefficients (CD) following the direct measuring method [36,37,38]. Different from previous studies, standard force sensors are applied to compose four synchronized force–velocity measuring systems in the current experiment. These standard force sensors are robust and suitable for flume applications. The composed force–velocity measuring systems can provide synchronized force–velocity measurement. Although one of the ADV instruments failed, the other three ADVs functioned well during the experiment. Importantly, an automatic algorithm was developed to realign the obtained force and velocity signals for direct CD deviation. This algorithm is expected to be able to accommodate a variety of velocity measuring techniques, providing possibilities to extend current application range. The developed force–velocity measuring systems and the automatic realignment algorithm may assist future experiments on vegetation–wave interactions for better understanding and prediction of vegetation-induced wave dissipation.

Author Contributions

P.Y., Z.H., B.H. and Q.Y. conceived and designed the methodology and experiments. P.Y., H.C., C.T. and L.R. performed the experiments. P.Y., L.R. and B.H. analyzed the data. All the authors contributed to the writing and editing of the manuscript.

Funding

The authors gratefully acknowledge financial support of the National Key R&D Program of China (No. 2016YFC0402607), the National Natural Science Foundation of China (No. 51609269, 51709288), the Joint Research Projects NSFC (No. 51761135022)—NWO (No. ALWSD.2016.026)—EPSRC (No. EP/R024537/1): Sustainable Deltas, and the Open Foundation of State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering (Grant No. 2016491711), and the 2017 Guangdong Oceanic and Fishery Department scientific research and technological development program.

Acknowledgments

We thank the three anonymous reviewers for their comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. MATLAB Code for Force and Velocity Data Realignment

This appendix provides the MATLAB (R2016a, The MathWorks, Natick, U.S.) script that can be applied to realign force and velocity timeseries for CD derivations.
for i=1:3 % number of Synchronized Force–Velocity Measurement systems
 NN=0;
 D=0;
 DD=0;
 while NN<=30 % number of loops
  NN=NN+1
  Delay(i,NN)=round(abs(DD)+D);
%Adjust the velocity bases on the phase difference, while keep the force is not change
   if i==1
    secADV1=AllADV1new(orpo1−Delay(i,NN):orpo1+8*TT*80−Delay(i,NN),:);
    secADV(:,1)=secADV1(:,2);
   end
 
   if i==2
    secADV2=AllADV2new(orpo2−Delay(i,NN):orpo2+8*TT*80−Delay(i,NN),:);
    secADV(:,2)=secADV2(:,2);
   end
 
   if i==3
    secADV3=AllADV3new(orpo3−Delay(i,NN):orpo3+8*TT*80−Delay(i,NN),:);
    secADV(:,3)=secADV3(:,2);
   end
   % Calculation of inertia force
    Fm(:,i)=2*1000*pi*0.02*0.02*hv*(diff(secADV(:,i))/0.0125)/4;
   % Calculation of drag force
    Fd(:,i)=secCLN(2:end,i)−Fm(:,i)
   %Calculation of in−time Cd
    Cd(:,i)=2*Fd(:,i)./(1000*0.02*hv*secADV(2:end,i).*abs(secADV(2:end,i)));
   %Intercept the peak value of velocity
    [maxE,minE]=peakdet(secADV(:,i),0.3*(max(secADV(:,i))−min(secADV(:,i))));
   %Intercept the peak value of drag force
    [maxFd,minFd]=peakdet(Fd(:,i),0.3*(max(Fd(:,i))−min(Fd(:,i))));
 
   if maxE(1)<10
    for kk=2:length(maxE)
     maxE(kk−1)=maxE(kk);
    end
    maxE(kk)=NaN;
   end
 
   if minE(1)<10
    for kk=2:length(minE)
     minE(kk−1)=minE(kk);
    end
    minE(kk)=NaN;
   end
 
   if maxFd(1)<10
    for kk=2:length(maxFd)
     maxFd(kk−1)=maxFd(kk);
    end
    maxFd(kk)=NaN;
   end
 
   if minFd(1)<10
    for kk=2:length(minFd)
     minFd(kk−1)=minFd(kk);
    end
    minFd(kk)=NaN;
   end
    loc_E(:,i)=cat(1,maxE(1:2,1),minE(1:2,1));
    loc_Fd(:,i)=cat(1,maxFd(1:2,1),minFd(1:2,1));
    %Phase difference between drag force and velocity
    delay(i,NN)=mean(loc_E(:,i)−loc_Fd(:,i))
    clear maxFd minFd maxE minE
 
     D=abs(delay(i,NN));
     DD=Delay(i,NN);
  end
end

References

  1. Bouma, T.J.; van Belzen, J.; Balke, T.; Zhu, Z.; Airoldi, L.; Blight, A.J.; Davies, A.J.; Galvan, C.; Hawkins, S.J.; Hoggart, S.P.G.; et al. Identifying knowledge gaps hampering application of intertidal habitats in coastal protection: Opportunities & steps to take. Coast. Eng. 2014, 87, 147–157. [Google Scholar] [CrossRef]
  2. Balke, T.; Friess, D.A. Geomorphic knowledge for mangrove restoration: A pan-tropical categorization. Earth Surf. Process. Landf. 2016, 41, 231–239. [Google Scholar] [CrossRef]
  3. Barbier, E.B. Valuing the storm protection service of estuarine and coastal ecosystems. Ecosyst. Serv. 2015, 11, 32–38. [Google Scholar] [CrossRef]
  4. Temmerman, S.; Meire, P.; Bouma, T.J.; Herman, P.M.J.; Ysebaert, T.; De Vriend, H.J. Ecosystem-based coastal defence in the face of global change. Nature 2013, 504, 79–83. [Google Scholar] [CrossRef] [PubMed]
  5. Horstman, E.M.; Dohmen-Janssen, C.M.; Narra, P.M.F.; van den, B.; Siemerink, M.; Hulscher, S.J. Wave attenuation in mangroves: A quantitative approach to field observations. Coast. Eng. 2014, 94, 47–62. [Google Scholar] [CrossRef]
  6. Zhang, X.; Chua, V.P.; Cheong, H.-F. Hydrodynamics in mangrove prop roots and their physical properties. J. Hydro-Environ. Res. 2015, 9, 281–294. [Google Scholar] [CrossRef]
  7. Zhang, X.; Chua, V.P.; Cheong, H.-F. Geometrical and material properties of Sonneratia alba mangrove roots. Trees-Struct. Funct. 2014, 29, 285–297. [Google Scholar] [CrossRef]
  8. Mazda, Y.; Wolanski, E.; King, B.; Sase, A.; Ohtsuka, D.; Magi, M. Drag force due to vegetation in mangrove swamps. Mangroves Salt Marshes 1997, 1, 193–199. [Google Scholar] [CrossRef]
  9. Alongi, D.M. Mangrove forests: Resilience, protection from tsunamis, and responses to global climate change. Estuar. Coast. Shelf Sci. 2008, 76, 1–13. [Google Scholar] [CrossRef]
  10. Möller, I.; Kudella, M.; Rupprecht, F.; Spencer, T.; Paul, M.; van Wesenbeeck, B.K.; Wolters, G.; Jensen, K.; Bouma, T.J.; Miranda-Lange, M.; Schimmels, S. Wave attenuation over coastal salt marshes under storm surge conditions. Nat. Geosci. 2014, 7, 727–731. [Google Scholar] [CrossRef]
  11. Arkema, K.K.; Guannel, G.; Verutes, G.; Wood, S.A.; Guerry, A.; Ruckelshaus, M.; Kareiva, P.; Lacayo, M.; Silver, J.M. Coastal habitats shield people and property from sea-level rise and storms. Nat. Clim. Chang. 2013, 3, 913–918. [Google Scholar] [CrossRef]
  12. Zhu, Z.; Van, B.; Hong, T.; Kunihiro, T.; Ysebaert, T.; Herman, P.M.J.; Bouma, T.J. Sprouting as a gardening strategy to obtain superior supplementary food: Evidence from a seed-caching marine worm. Ecology 2016, 97, 3278–3284. [Google Scholar] [CrossRef] [PubMed]
  13. Young, I.R.; Zieger, S.; Babanin, A.V. Global Trends in Wind Speed and Wave Height. Science 2011, 332, 451–455. [Google Scholar] [CrossRef] [PubMed]
  14. Cheng, Y.; Ezer, T.; Hamlington, B.D. Sea Level Acceleration in the China Seas. Water 2016, 8, 293. [Google Scholar] [CrossRef]
  15. Lin-Ye, J.; Garcia-Leon, M.; Gracia, V.; Isabel Ortego, M.; Stanica, A.; Sanchez-Arcilla, A. Multivariate Hybrid Modelling of Future Wave-Storms at the Northwestern Black Sea. Water 2018, 10, 221. [Google Scholar] [CrossRef]
  16. Wolf, J.; Lowe, J.; Howard, T. Climate downscaling: Local mean sea level, surge and wave modelling. In Broad Scale Coastal Simulation: New Techniques to Understand and Manage Shorelines in the Third Millennium; Springer: Dordrecht, The Netherlands, 2015; pp. 79–102. ISBN 978-94-007-5258-0. [Google Scholar]
  17. Wolf, J.; Brown, J.M.; Bolaños, R.; Hedges, T.S. Waves in Coastal and Estuarine Waters. In Treatise on Estuarine and Coastal Science; Wolanski, E., McLusky, D., Eds.; Academic Press: Waltham, MA, USA, 2011; pp. 171–212. ISBN 978-0-08-087885-0. [Google Scholar]
  18. Van Loon-Steensma, J.M.; Slim, P.A.; Decuyper, M.; Hu, Z. Salt-marsh erosion and restoration in relation to flood protection on the Wadden Sea barrier island Terschelling. J. Coast. Conserv. 2014, 18, 1–16. [Google Scholar] [CrossRef]
  19. Hu, Z.; Van Belzen, J.; Van Der Wal, D.; Balke, T.; Wang, Z.B.; Stive, M.; Bouma, T.J. Windows of opportunity for salt marsh vegetation establishment on bare tidal flats: The importance of temporal and spatial variability in hydrodynamic forcing. J. Geophys. Res. Biogeosci. 2015, 120, 1450–1469. [Google Scholar] [CrossRef] [Green Version]
  20. Zhang, M.; Zhang, H.; Zhao, K.; Tang, J.; Qin, H. Evolution of wave and tide over vegetation region in nearshore waters. Ocean Dyn. 2017, 67, 973–988. [Google Scholar] [CrossRef]
  21. Dalrymple, R.; Kirby, J.; Hwang, P. Wave Diffraction Due to Areas of Energy Dissipation. J. Waterw. Port Coast. Ocean Eng. 1984, 110, 67–79. [Google Scholar] [CrossRef]
  22. Koch, D.L.; Ladd, A.J.C. Moderate Reynolds number flows through periodic and random arrays of aligned cylinders. J. Fluid Mech. 1997, 349, 31–66. [Google Scholar] [CrossRef]
  23. Morison, J.R.; O’Brien, M.P.; Johnson, J.W.; Schaaf, S.A. The Force Exerted by Surface Waves on Piles. J. Pet. Technol. 1950, 2, 149–154. [Google Scholar] [CrossRef]
  24. Nepf, H.M. Flow over and through Biota. In Treatise on Estuarine and Coastal Science; Wolanski, E., McLusky, D., Eds.; Academic Press: Waltham, MA, USA, 2011; pp. 267–288. ISBN 978-0-08-087885-0. [Google Scholar]
  25. Augustin, L.N.; Irish, J.L.; Lynett, P. Laboratory and numerical studies of wave damping by emergent and near-emergent wetland vegetation. Coast. Eng. 2009, 56, 332–340. [Google Scholar] [CrossRef]
  26. Tang, J.; Shen, S.; Wang, H. Numerical model for coastal wave propagation through mild slope zone in the presence of rigid vegetation. Coast. Eng. 2015, 97, 53–59. [Google Scholar] [CrossRef]
  27. Cao, H.; Feng, W.; Hu, Z.; Suzuki, T.; Stive, M.J.F. Numerical modeling of vegetation-induced dissipation using an extended mild-slope equation. Ocean Eng. 2015, 110, 258–269. [Google Scholar] [CrossRef]
  28. van Loon-Steensma, J.M.; Hu, Z.; Slim, P.A. Modelled Impact of Vegetation Heterogeneity and Salt-Marsh Zonation on Wave Damping. J. Coast. Res. 2016, 32, 241–252. [Google Scholar] [CrossRef]
  29. Peruzzo, P.; De Serio, F.; Defina, A.; Mossa, M. Wave Height Attenuation and Flow Resistance Due to Emergent or Near-Emergent Vegetation. Water 2018, 10, 402. [Google Scholar] [CrossRef]
  30. Tang, J.; Causon, D.; Mingham, C.; Qian, L. Numerical study of vegetation damping effects on solitary wave run-up using the nonlinear shallow water equations. Coast. Eng. 2013, 75, 21–28. [Google Scholar] [CrossRef]
  31. Kobayashi, N.; Raichle, A.W.; Asano, T. Wave attenuation by vegetation. J. Waterw. Port Coast. Ocean Eng.-ASCE 1993, 119, 30–48. [Google Scholar] [CrossRef]
  32. Méndez, F.J.; Losada, I.J.; Losada, M.A. Hydrodynamics induced by wind waves in a vegetation field. J. Geophys. Res. Oceans 1999, 104, 18383–18396. [Google Scholar] [CrossRef] [Green Version]
  33. Ozeren, Y.; Wren, D.G.; Wu, W. Experimental investigation of wave attenuation through model and live vegetation. J. Waterw. Port Coast. Ocean Eng. 2014, 140. [Google Scholar] [CrossRef]
  34. Anderson, M.E.; Smith, J.M. Wave attenuation by flexible, idealized salt marsh vegetation. Coast. Eng. 2014, 83, 82–92. [Google Scholar] [CrossRef]
  35. Losada, I.J.; Maza, M.; Lara, J.L. A new formulation for vegetation-induced damping under combined waves and currents. Coast. Eng. 2016, 107, 1–13. [Google Scholar] [CrossRef]
  36. Infantes, E.; Orfila, A.; Bouma, T.J.; Simarro, G.; Terrados, J. Posidonia oceanica and Cymodocea nodosa seedling tolerance to wave exposure. Limnol. Oceanogr. 2011, 56, 2223–2232. [Google Scholar] [CrossRef] [Green Version]
  37. Hu, Z.; Suzuki, T.; Zitman, T.; Uijttewaal, W.; Stive, M. Laboratory study on wave dissipation by vegetation in combined current-wave flow. Coast. Eng. 2014, 88, 131–142. [Google Scholar] [CrossRef]
  38. Maza, M.; Adler, K.; Ramos, D.; Garcia, A.M.; Nepf, H. Velocity and Drag Evolution From the Leading Edge of a Model Mangrove Forest. J. Geophys. Res. Oceans 2017, 122, 9144–9159. [Google Scholar] [CrossRef]
  39. Dean, R.; Dalrymple, R. Water Wave Mechanics for Engineers and Scientists; Advanced Series on Ocean Engineering; World Scientific: Tokyo, Japan, 1991. [Google Scholar]
  40. Strom, K.B.; Papanicolaou, A.N. ADV measurements around a cluster microform in a shallow mountain stream. J. Hydraul. Eng.-ASCE 2007, 133, 1379–1389. [Google Scholar] [CrossRef]
  41. Persoons, T.; O’Donovan, T.S. High Dynamic Velocity Range Particle Image Velocimetry Using Multiple Pulse Separation Imaging. Sensors 2011, 11, 1–18. [Google Scholar] [CrossRef] [PubMed]
  42. Martino, R.; Paterson, A.; Piva, M. Double-average mean flow and local turbulence intensity profiles from PIV measurements for an open channel flow with rigid vegetation. Environ. Fluid Mech. 2012, 12, 45–62. [Google Scholar] [CrossRef]
  43. You, Z.; Chen, Y. The Use of Tactile Sensors and PIV Analysis for Understanding the Bearing Mechanism of Pile Groups. Sensors 2018, 18, 476. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) Flume experiment set-up. The numbers 1–4 indicate the locations of the in-phase force–velocity measurement; (b) A force sensor connected to a PVC pipe; (c) Mimicked vegetation canopy constructed by PVC pipes; (d) Instrument deployments in the flume without mimicked vegetation. “WG” stands for wave gauge, “ADV” stands for acoustic doppler velocimeter, and “FS” stands for force sensor. The dashed line indicates the ADV and force sensor are placed at the same cross-section to ensure in-phase measurements.
Figure 1. (a) Flume experiment set-up. The numbers 1–4 indicate the locations of the in-phase force–velocity measurement; (b) A force sensor connected to a PVC pipe; (c) Mimicked vegetation canopy constructed by PVC pipes; (d) Instrument deployments in the flume without mimicked vegetation. “WG” stands for wave gauge, “ADV” stands for acoustic doppler velocimeter, and “FS” stands for force sensor. The dashed line indicates the ADV and force sensor are placed at the same cross-section to ensure in-phase measurements.
Water 10 00906 g001
Figure 2. Flow chart to realign velocity and force data signals. The algorithm is provided in the Appendix A as a MATLAB script.
Figure 2. Flow chart to realign velocity and force data signals. The algorithm is provided in the Appendix A as a MATLAB script.
Water 10 00906 g002
Figure 3. (a) Measured velocity profile for case wave0312 with 3 cm wave height and 1.2 s wave period, Umin is the highest wave orbital velocity in negative direction (oppose to wave propagation), Umax is the highest wave orbital velocity in positive direction (same as wave propagation), and Uw is the amplitude of wave orbital velocity; (b) Measured velocity profile for case wave0709 with 7 cm wave height and 0.9 s wave period.
Figure 3. (a) Measured velocity profile for case wave0312 with 3 cm wave height and 1.2 s wave period, Umin is the highest wave orbital velocity in negative direction (oppose to wave propagation), Umax is the highest wave orbital velocity in positive direction (same as wave propagation), and Uw is the amplitude of wave orbital velocity; (b) Measured velocity profile for case wave0709 with 7 cm wave height and 0.9 s wave period.
Water 10 00906 g003
Figure 4. (a) Spatial variations of wave height (H) through the vegetation canopy (i.e., x = 0–8 m) indicated as green bars; (b) Spatial variations of the magnitude of wave orbital velocity (Uw) through the vegetation canopy. The first 3 Uw data points (x = 0–6 m) are obtained from ADV measurement, whereas the last data point at x = 7 m is obtained by using Equation (3) based on an average wave height between x = 6 m and 8 m in panel (a). The shown test case is wave0712 with 7 cm wave height and 1.2 s wave period.
Figure 4. (a) Spatial variations of wave height (H) through the vegetation canopy (i.e., x = 0–8 m) indicated as green bars; (b) Spatial variations of the magnitude of wave orbital velocity (Uw) through the vegetation canopy. The first 3 Uw data points (x = 0–6 m) are obtained from ADV measurement, whereas the last data point at x = 7 m is obtained by using Equation (3) based on an average wave height between x = 6 m and 8 m in panel (a). The shown test case is wave0712 with 7 cm wave height and 1.2 s wave period.
Water 10 00906 g004
Figure 5. (ad) Raw velocity and total force data measured at four locations (1–4) in the order of the wave propagation in the mimicked vegetation canopy; (e) The shaded area in panel (c) is blown up in panel (e) for detailed analysis, where inertia force and drag force were derived based on non-synchronized data. The time shift (∆t) between the FD and U is 0.26 s. The shown test case is wave0712 with 7 cm wave height and 1.2 s wave period. The velocity data in panel (d) is not included, due to the ADV measurement failure at location 4.
Figure 5. (ad) Raw velocity and total force data measured at four locations (1–4) in the order of the wave propagation in the mimicked vegetation canopy; (e) The shaded area in panel (c) is blown up in panel (e) for detailed analysis, where inertia force and drag force were derived based on non-synchronized data. The time shift (∆t) between the FD and U is 0.26 s. The shown test case is wave0712 with 7 cm wave height and 1.2 s wave period. The velocity data in panel (d) is not included, due to the ADV measurement failure at location 4.
Water 10 00906 g005
Figure 6. (ad) Time-varying U, FD, FM, and F data over two wave periods. The gray lines are before realignment, and the red ones are after realignment. The vertical double-arrowed lines in panel (a) and (b) indicate the synchronization status before and after the realignment; (e) Time-varying CD derived based on realigned U and FD data. The shown test case is wave0712 with 7 cm wave height and 1.2 s wave period.
Figure 6. (ad) Time-varying U, FD, FM, and F data over two wave periods. The gray lines are before realignment, and the red ones are after realignment. The vertical double-arrowed lines in panel (a) and (b) indicate the synchronization status before and after the realignment; (e) Time-varying CD derived based on realigned U and FD data. The shown test case is wave0712 with 7 cm wave height and 1.2 s wave period.
Water 10 00906 g006
Figure 7. Sensitivity analysis of the number of loops used in the realignment algorithm. Case 1 is the case wave0712 with 7 cm wave height and 1.2 s wave period, and case 2 is the case wave0512 with 5 cm wave height and 1.2 s wave period.
Figure 7. Sensitivity analysis of the number of loops used in the realignment algorithm. Case 1 is the case wave0712 with 7 cm wave height and 1.2 s wave period, and case 2 is the case wave0512 with 5 cm wave height and 1.2 s wave period.
Water 10 00906 g007
Figure 8. Time-varying U, PD, and PM over two wave periods. The shown test case is wave0712, with 7 cm wave height and 1.2 s wave period. The averaged work done by drag force (WD) and inertia force (WM) over the shown two periods is 2.0   ×   10 3   J and 2.0   ×   10 4   J , respectively. The derived period-averaged CD is 3.25.
Figure 8. Time-varying U, PD, and PM over two wave periods. The shown test case is wave0712, with 7 cm wave height and 1.2 s wave period. The averaged work done by drag force (WD) and inertia force (WM) over the shown two periods is 2.0   ×   10 3   J and 2.0   ×   10 4   J , respectively. The derived period-averaged CD is 3.25.
Water 10 00906 g008
Figure 9. Comparison between reproduced total acting force and measured total force. The red solid line is the measured total force; The black dash line is quantified by using derived period-averaged CD in Equation (1) (i.e., Frep); The blue dash line is quantified by using CD = 1 in Equation (1) for reference (i.e., Frep). The shown test case is wave0712, with 7 cm wave height and 1.2 s wave period.
Figure 9. Comparison between reproduced total acting force and measured total force. The red solid line is the measured total force; The black dash line is quantified by using derived period-averaged CD in Equation (1) (i.e., Frep); The blue dash line is quantified by using CD = 1 in Equation (1) for reference (i.e., Frep). The shown test case is wave0712, with 7 cm wave height and 1.2 s wave period.
Water 10 00906 g009
Figure 10. Comparison between the maximum reproduced total acting force (Frep) and the measured maximum total force at three functional measuring locations in all the tested cases.
Figure 10. Comparison between the maximum reproduced total acting force (Frep) and the measured maximum total force at three functional measuring locations in all the tested cases.
Water 10 00906 g010
Table 1. Comparison between known weights and force sensor reading.
Table 1. Comparison between known weights and force sensor reading.
DirectionKnown Weights (g)1st Reading a (g)2nd Reading b (g)3rd Reading c (g)Mean Reading (g)Absolute Error (g)Relative Error
+5 g5.105.005.005.030.030.60%
10 g10.1010.1010.0010.070.070.70%
20 g20.1020.0020.0020.030.030.15%
50 g49.9049.9049.9049.900.100.02%
5 g−5.00−5.00−5.10−5.03−0.03−0.60%
10 g−10.00−10.10−10.10−10.07−0.07−0.70%
20 g−20.00−20.00−20.00−20.0000%
50 g−50.0−49.9−49.9−49.930.070.14%
a,b,c The 1st, 2nd, and 3rd time readings were taking when the weights were put at the bottom, middle, and tip of the testing pipe, respectively.
Table 2. Period-averaged CD in all the tested cases.
Table 2. Period-averaged CD in all the tested cases.
Test NumberWave Height (m)Wave Period (s)CD at Location 1CD at Location 2CD at Location 3Space-Mean CD Standard Deviation
10.030.65.4110.0210.148.527.28
20.030.93.194.887.095.053.82
30.031.23.602.795.253.881.56
40.050.67.615.835.966.460.98
50.050.93.822.363.623.270.62
60.051.22.843.033.973.280.37
70.070.63.433.017.044.494.90
80.070.91.971.832.942.250.37
90.071.21.772.773.512.680.76
100.090.63.025.766.104.962.86
110.090.91.261.792.641.890.49
120.091.21.442.543.002.330.64

Share and Cite

MDPI and ACS Style

Yao, P.; Chen, H.; Huang, B.; Tan, C.; Hu, Z.; Ren, L.; Yang, Q. Applying a New Force–Velocity Synchronizing Algorithm to Derive Drag Coefficients of Rigid Vegetation in Oscillatory Flows. Water 2018, 10, 906. https://doi.org/10.3390/w10070906

AMA Style

Yao P, Chen H, Huang B, Tan C, Hu Z, Ren L, Yang Q. Applying a New Force–Velocity Synchronizing Algorithm to Derive Drag Coefficients of Rigid Vegetation in Oscillatory Flows. Water. 2018; 10(7):906. https://doi.org/10.3390/w10070906

Chicago/Turabian Style

Yao, Peng, Hui Chen, Bensheng Huang, Chao Tan, Zhan Hu, Lei Ren, and Qingshu Yang. 2018. "Applying a New Force–Velocity Synchronizing Algorithm to Derive Drag Coefficients of Rigid Vegetation in Oscillatory Flows" Water 10, no. 7: 906. https://doi.org/10.3390/w10070906

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