Analysis of the Bending Height of Flexible Marine Vegetation

: Marine vegetation is increasingly viewed as a living shoreline that protects coastal communities and ecosystems from the damaging effects of wave energy. Many studies have explored the potential of marine vegetation in terms of reducing wave height, but more work is needed. Here, we used particle image velocimetry, fluid–structure interaction simulation, and multiple regression analysis to estimate the bending behaviors of flexible marine vegetation in water flow, and we predicted the wave height reduction in the downstream vegetation meadow. We considered different vegetation types and water flow velocities, constructed a total of 64 cases, and derived a multiple regression equation that simply estimates the vegetation bending height with a tolerance of ~10%. When the bending height rather than the vegetation height was applied, wave height reduction was alleviated by 1.08–9.23%. Thus, flexible vegetation reduced wave height by up to ~10% less than rigid vegetation in our investigation range. This implies that the impact of bending behavior becomes more pronounced with a larger vegetation meadow. The relative % decrease in wave height reduction was greater for fully submerged vegetation compared to partially submerged vegetation.


Introduction
Considerable efforts have been made over the past few decades to investigate the effects of climate change on coastal and terrestrial vegetation.These studies were conducted to better understand how environmental changes (e.g., temperature) affect vegetation covers and distributions [1][2][3][4].Significant efforts have also been made to explore the impacts of marine vegetation (e.g., seagrass, mangroves, salt marshes, kelp, and red algae) on the aquatic environment.Such studies have sought to better understand how marine vegetation regulates fluid flow patterns [5][6][7][8] and the transport of nutrients, pollutants, and sediments [9][10][11][12].Using various tools and techniques, researchers have explored the complex interactions between marine vegetation and fluid flow, yielding insights into the ecological and environmental benefits afforded by such ecosystems.Research in this area is becoming increasingly important in terms of protecting and preserving the health of the oceans and coastal areas.
Waves tend to increase in height as they approach the shore, triggering erosion and other damage to coastal habitats and infrastructure.Marine vegetation meadows serve as natural barriers that slow waves and reduce wave energies and heights [13,14].This not only protects against erosion but also provides habitats for many marine species.The roots of marine vegetation stabilize the sediments; these are then not washed away by strong currents or storms [15,16].Today, marine vegetation is increasingly viewed as a living shoreline that protects coastal communities and ecosystems from the damaging effects of wave energy [14,17].To explore how marine vegetation reduces wave height or energy, laboratory-scale experiments [18][19][20], field experiments [21,22], simulations [23,24], and theoretical calculations [25] have been employed using either living [26] or artificial [20,27] vegetation.Some studies used simple computations that treated vegetation as completely rigid [25,28,29], thus ignoring real-world flexibility.Other studies indeed investigated flexible vegetation [13,27], but only a few specific types.The distinction between rigid and flexible vegetation lies in their motion.Flexible vegetation can bend and sway in response to water flow forces, while rigid vegetation tends to maintain its original posture.Various studies have explored the movement of vegetation, revealing its significant impact on wave attenuation within vegetation meadows by reducing drag [30][31][32][33].A more general method is needed to quantify how flexible vegetation reduces wave height.There is a need for a more widely accepted approach, such as simulation, which has broad applicability, or an empirical equation to address this complexity.
Here, we studied the bending behaviors and reduced heights of flexible marine vegetation under water flow and predicted wave height reductions in the downstream vegetation meadow.This approach may optimally address how vegetation reduces wave height because it considers a reduced vegetation height (i.e., the bending height) rather than the original height (i.e., the vegetation height).We first used particle image velocimetry (PIV) in the laboratory to determine the bending height.We confirmed the results via fluid-structure interaction (FSI) simulation using the ANSYS Fluent software 2023 R1 package (ANSYS, Inc., Canonsburg, PA, USA).Based on the data obtained, we established a multiple regression equation that rather simply yields the bending heights.We used these values to derive the wave height reductions downstream and compared the results to those obtained using the vegetation heights.The contribution of this study can be highlighted as follows.First, we estimated the bending behaviors of flexible vegetation in water flow; the estimated equation was effective in calibrating the vegetation height due to bending behavior.Second, we predicted the effect of bending behavior on wave height; rigid vegetation proved to be more effective in reducing wave height than flexible vegetation.Finally, we observed that the relative % decrease in wave height reduction is greater for fully submerged vegetation compared to partially submerged vegetation.

Artificial Vegetation
Table 1 lists the vegetation species examined in published studies; these varied widely in terms of vegetation height (0.07-0.55 m), width (1.8-38.1 mm), and Young's modulus (0.56-3530 ± 390 MPa).The table also lists the comparable features of alternative (or artificial) materials used.For example, materials used to represent artificial rigid or flexible vegetables included polyethylene [18], polyvinyl chloride [19], and bamboo dowels with silicon sealants [34].It is difficult to perform experiments using live vegetation; real plants are irregular in shape, susceptible to weather conditions, and decompose within a few days after removal from the natural environment.We used thin polystyrene, a synthetic polymer of styrene (an aromatic hydrocarbon) monomers, to simulate the behaviors of real vegetation in water flow; the characteristics are listed in Table 2.We created nine types of artificial vegetation by fixing the thickness (0.001 m) and changing the height (0.250, 0.300, and 0.350 m) and width (0.010, 0.015, and 0.020 m).

Experimental Setup
We performed PIV to quantify the bending heights.PIV is a form of optical measurement that instantaneously measures the velocity field in a flow region [38,39].In a typical PIV experiment, a fluid is illuminated using a pulsed laser light that is scattered by small tracer particles in the fluid, and then the flow field is captured in two images of a highspeed camera taken at two different times during the laser pulse [40].Figure 1 shows the experimental setup comprising a PIV pulse laser unit (DualPower 200-15; Dantec, Skovlunde, Denmark), pulse generator, synchronizer, charge-coupled device (CCD) camera with a 15 Hz sampling rate (EO 4M-32 FlowSense; Dantec; Skovlunde, Denmark), and host computer.The seed particles were hollow and smooth borosilicate glass spheres that both followed the flow and scattered light that was captured by the CCD camera, which converted the light intensities into electrical signals stored in computer memory after analog-to-digital conversion.The synchronizer timed the PIV system's operation; imaging and laser lighting operated simultaneously.The computer featured PIV software (DynamicStudio 7.0) that processed the image pairs to produce velocity flow maps.The micron-sized seed particles scattered light that was imaged using a laser light sheet less than 1 mm in thickness.
The water flume (0.5 m wide, 0.5 m high, and 10.0 m long) was made of steel and tempered glass.When choosing the water flow velocity, we considered previous studies that employed 0.08-0.25 m s −1 [35] and 0.1-0.2m s −1 [41] to replicate real-world flows by location and season.We used velocities of 0.089-0.155m s −1 , thus within the abovementioned ranges for coastal environments.We used an inverter motor attached to a water tank (Figure 1).At rotation frequencies of 30, 40, and 50 Hz, the flow velocities were 0.089, 0.123, and 0.155 m s −1 , respectively.
Using these three water velocities and nine types of artificial vegetation, we evaluated 27 experimental cases (Table 3).Each vegetation was installed in the center of the flume, and five measurements were obtained at each water flow velocity to confirm reproducibility.No waves were generated, and all vegetation was fully submerged in the flume.The dimensions of the velocity profile mapping region were 27 L × 33 L, where L = 10.765mm, and that region was vertically divided into ten local locations, thus from line 0 L to line 27 L at intervals of 3 L. Accordingly, ten local velocities, from  1 to  10 , were obtained (Figure 2a); the window used to measure water velocity was ~290 × 355 mm (Figure 2b).The laser illuminated the vegetation, and the CCD camera captured vegetation movement and velocity.We obtained the bending heights (i.e., the heights after bending by the water flow) by observing the deformations apparent on PIV images when the deformation had stabilized.

Fluid-Structure Interaction Simulation
FSI refers to the interaction between a movable or deformable structure and an internal or surrounding fluid flow.In general, an FSI problem is structured to ensure that the field equations for the solid and fluid affect and modify each other.A solid not only induces fluid motion but also modifies that motion by creating stresses at the solid-fluid interface [42].Two conditions must be met when balancing the motions of a solid and a fluid.First, the solid and fluid movements must match at the interface.Second, the mutual forces exerted at the interface must be equal.In an FSI simulation, the fluid and solid fields are solved independently at each time step and then coupled until the required conditions at the interface are satisfied.Initially, the pressure from the fluid field affects the solid, causing the solid to deform or move; these changes are then back-transferred to the fluid field, changing the properties thereof [43].A calculation cycle is complete when the modified pressure again affects the solid.FSI simulations effectively model the real behaviors of flexible objects and their surrounding environments and have thus found applications in automobile, aircraft, spacecraft, engine, bridge design, and medical fields [44][45][46].
The FSI equation is then combined with the equations governing the fluid and structure domains [47].The fluid domain is described using Navier-Stokes, continuity, volume fraction, and energy equations, as shown in Equations ( 1)-( 4), respectively: where  is time,  the fluid density,  ⃗ the velocity vector, ∇ the gradient,  the pressure, ̿ the stress tensor,  ⃗ gravitational acceleration,  ⃗ an external body force vector,   a mass source term,  the fluid volume fraction,  the total energy,  eff the effective fluid thermal conductivity,  the temperature, ℎ the enthalpy,  ⃗ the diffusion flux vector, ̿ eff the viscous stress tensor, and  ℎ a volumetric energy source.The computational model of the solid domain features three-dimensional strain and nodal displacement equations, and a stress equation [47], as shown in Equations ( 5)-( 7): where  is the strain displacement,  the four-dimensional gradient (time and space),  an element shape function,  the strain,  the nodal displacement,  the stress, and  the modulus of elasticity.

Simulation
We performed an FSI simulation to evaluate the experimental results.To model turbulence, the Reynolds stress model was chosen because it is both reliable and cost-effective [48].We used the "transient structural module" and the "fluent module".The former module determines the displacement or deformation of a solid body caused by fluid pressure within the domain.The result is then transferred to the fluent module, which recalculates the characteristics of the fluid domain.Finally, all results are back-transferred to the solid body; this completes the cycle.The following assumptions were made.First, for the transient structural module, the vegetation was assumed to be uniform in terms of shape and homogeneous; the shape did not change from bottom to top, and the properties thereof were identical throughout the vegetation height.We also assumed that the connection between vegetation and the seabed was rigid; we thus imposed a fixed-support boundary condition.Five surfaces of each vegetation, including the front, back, two sides, and the top, served as the interfaces between the fluid domain and the vegetation.Finally, it was assumed that the seabed was absolutely flat, and friction occurred between the flow and the seabed.
We carefully considered the modeling parameters; we performed tests that independently identified a reasonable domain length, the number of iteration steps, the analysis time, and the number of elements (Figure 3).In this figure, the percentage represents the relative error between the current deformation of the calculation point (  ) and the selected point (  ) and is defined by (  −   )   × 100(% ⁄ ).The top, horizontal vegetation deformation reached convergence values at an early stage with respect to the domain length and number of elements, respectively.However, the deformation was sensitive to the number of iteration steps and analysis time; hence, the selected points were determined to minimize the relative errors of the calculation points.As a result, we chose a domain length of 3.0 m, 200 iteration steps, an analysis time of 5 s, and ~2 million elements.The domain height and width were the actual dimensions of the rectangular flume (50 cm in height and 50 cm in width).Here, the relative % errors are with respect to the selected point; those errors, by definition, must all be positive; however, the "-" sign was intentionally added to highlight the differences with respect to the reference values.
The boundary conditions included an inlet, an outlet, symmetry, and a no-slip wall (Figure 4a).The inlet was the boundary from which the flow entered; the inlet flow velocity was 0.089, 0.123, or 0.155 m s −1 , identical to those of the PIV experiments.The outlet boundary forced the flow to exit through an outlet.The symmetry condition ensured that the physical geometry and the anticipated flow pattern exhibited mirror symmetry [49].The no-slip wall condition was applied to the lateral and bottom sides.We created mesh layers, gradually decreasing the mesh size from 15 mm in the fluid domain to 5 mm in the solid domain (each vegetation) (Figure 4b).A tetrahedral mesh was used because this is optimal in terms of automated meshing, stress visualization, and predictions of contact pressure and shear stress [50].We used three dynamic mesh schemes: smoothing, layering, and re-meshing to explore the movements of individual vegetation [39].Smoothing enables adjustment of a zone mesh by deforming or moving the boundary; the interior mesh nodes are displaced, but the node number is maintained, as are all connections.Dynamic layering facilitates the specification of an ideal layer height for each moving boundary.Re-meshing involves the clustering of cells that fail to meet skewness or size criteria; such cells or faces are subsequently re-meshed [51].We used WebplotDigitizer to extract vegetation bending heights [52].

Experimental and Simulation Results
Figure 5(a1-a10) show ten local velocity distributions at an inlet velocity of 0.089 m s −1 measured at ten vertical lines (Figure 2a).The horizontal axis indicates the normalized velocity     ⁄ , which is the ratio of each local velocity at the  ℎ vertical line to the inlet velocity   .Here,  ranges from 1 to 10 and  from 1 to 3 (e.g.,  1 = 0.089 m s −1 ).The vertical axis indicates the vertical location from 0 to 33 L (Figure 2a).Similarly, Figure 5(b1-b10, c1-c10) show the velocity profiles for  2 = 0.123 m s −1 and  3 = 0.155 m s −1 , respectively.In each figure, the velocity near the boundary (the bottom of the flume) is smaller than elsewhere because of friction between the water flow and the bottom of the flume.Table 4 shows the vegetation bending height estimations for all 27 experimental cases, derived by sequentially combining the nine vegetation with the three velocities.For clarity, the height of vegetation after bending is the "bending height (ℎ  )" and the original height the "vegetation height (ℎ  )".From the PIV measurement, all bending heights were less than the vegetation heights, i.e., 86.86-96.00% of the latter heights (Table 4).

Case
To further assess accuracy, we derived error factors, i.e., the ratio of the FSI and PIV bending heights (Figure 7).Of the 27 cases, only 3 (Exp_12, Exp_22, and Exp_26) exhibited errors exceeding 5%; the remaining 24 had errors less than 5%.The simulation was thus accurate, inexpensive, and applicable to various physical materials and geometries for the current study.This demonstrated the physical effects of water flow on vegetation flow.It should be noted that the error gradually increases as the bending height increases.This is mainly due to the fixed boundary condition in the simulation.The vegetation bottom was ideally modeled as a fixed end; hence, it did not exactly reflect the experiment condition.Moreover, in the simulation, the material was assumed to be ideally isotropic, homogeneous and consistent across all the cross-sections.Such simulation conditions were not fully realized in the experiment.In this sense, the simulation may not adequately handle a large meadow with a high shoot density.

Regression Analysis
We performed 27 additional simulations (Table 5) to create a large dataset for multiple regression analysis; we considered the densities, Young's moduli, geometrical properties, and flow velocities of previous studies (Table 1).Of the 27 additional results, 15 exhibited rather small horizontal deformations attributable to the chosen densities, Young's moduli, vegetation heights, second moments of cross-sectional areas, and flow velocities; the bending heights were nearly identical to the vegetation heights.These data were used for multiple regression analysis, as were the 27 experimental data (Table 4).We considered vegetation height (ℎ  ), vegetation thickness (  ), vegetation width (  ), vegetation density (  ), Young's modulus (  ), and water velocity () as independent variables in the regression analysis because the bending behavior of a long and slender vegetation correlates with cross-sectional area, stiffness, and external force.
We applied the Buckingham Π theorem [53,54] to find the relationship of the bending height with the relevant variables.First, we considered the bending height (ℎ  ) of vegetation as a function of the six independent variables, as shown in Equation ( 8): We carried out the factor analysis to understand the complex relationship among the six independent variables.As a result, we identified the latent factors as follows: (a) vegetation height (ℎ  ) and width (  ) representing the cross-section property of the vegetation stem; (b) water velocity () and vegetation density (  ) representing the external force; and (c) vegetation height (ℎ  ) and Young's modulus (  ) representing the stiffness.
The dimensions of those seven variables, including ℎ  , are the combinations of M (mass), L (length), and T (time), which are represented by the units of kilogram (kg), meter (m), and second (s), respectively.According to the three fundamental dimensions and seven variables, we considered four non-dimensional parameters, i.e.,   ( = 1, 2, 3, and 4).Of those variables, we selected   , , and ℎ  as the repeating variables as these present the material property, flow property, and geometric property, respectively.Then, we took the remaining variables along with these variables one by one to create the non-dimensional groups.The first non-dimensional group is shown in Equation (9a): To make Equation (9a) dimensionless, we selected  =  = 0 and  = −1 ; hence, Equation (9a) becomes Equation (9b): Similarly, we obtained Equations ( 10) and ( 11) for the second and third non-dimensional groups, respectively: The fourth non-dimensional group is shown in Equation (12a), and by selecting  = −1,  = 0, and  = 2, its dimensionless form becomes Equation (12b).
According to the four non-dimensional parameters  1 = ℎ  /ℎ  ,  2 =   /ℎ  ,  3 =   /ℎ  , and  4 =    2 /  , we can establish Equation ( 13): Considering Equation ( 13), we can obtain Equations (14a) and (14b): ) = ln  0 +  1 ln ( We used the experimental and simulation data to derive the coefficients  0 =1.2831,  1 =0.0145,  2 =0.0676, and  3 = −0.0051;hence, the regression equation is determined as Equation ( 15): In Figure 8, the white and gray circles are the bending heights from multiple regression analysis of the 27 experimental data (Table 4) and 27 simulation data (Table 5), respectively.Most circles are within the 10% overestimation and underestimation lines and thus close to the 100% matching line; the multiple regression results were in good agreement with the experimental and simulation results.To further verify the accuracy of the regression analysis, the FSI simulation and multiple regression equation were used to calculate the bending heights for a further ten cases (blue circles in Figure 8).All blue circles are within the 10% lines and thus close to the 100% matching line, further validating the accuracy of the regression analysis.
Figure 9 shows the relative % errors of the bending heights predicted from the multiple regression equation using the 27 experimental results, the 27 simulations, and the 10 additional simulations.Of the 64 cases, 63 had errors < 10%, and 1 had an error of 10.71%.Thus, the multiple regression equation adequately estimated bending heights using the vegetation density, Young's modulus, the vegetation height, width, thickness, and the water flow velocity.⁄ ) of the present study with those of Fonseca and Kenworthy [32], Luhar and Nepf [30], and Zeller et al. [33].The data from Fonseca and Kenworthy [32], denoted by FK87 in the figure, were derived from a lab-scale experiment involving live vegetation exposed to unidirectional flow.The results from Luhar and Nepf [30], denoted by LN11, were obtained from a theoretical model predicting the bending height of a single artificial flexible vegetation stem under unidirectional flow.The data from Zeller et al. [33], denoted by Z14, were obtained from a lab-scale experiment using a blade model exposed to both wave and current conditions.

Effect of Bending Behavior on Wave Height
Most studies focused on wave attenuation solely through wave forces, overlooking the impact of currents.However, the presence of currents significantly diminishes the wave attenuation capacity, particularly in scenarios involving flexible vegetation [36,55].Currents induce vegetation bending, while waves induce swaying [31,56].Currents cause vegetation to move, reducing the front area perpendicular to the flow, lowering drag and reducing wave attenuation [30,31].Thus, to account for the effect of bending behavior on wave height, we estimated the reduction in front area of vegetation and wave attenuation.To quantify the significance of this effect, we considered the experiment conducted by Losada et al. [31].They used live vegetation (Puccinellia maritima), specified by Young's modulus of stems (13 MPa), Young's modulus of leaves (7.8 MPa), stem height (0.473 m), leaf height (0.230 m), number of leaves per stem (5.5), leaf width (0.003 m), and shoot density (2436).The water and wave properties [31] considered are shown in Table 6.To calculate the wave height, Losada et al. [31] developed Equation ( 16) based on linear wave theory and energy conservation: where () is the downstream wave height measured at a distance  (unit: m) away from the forefront of the vegetation meadow, and is  a damping factor.The damping factor is shown in Equation (17), which is a function of wave height at upstream   and two coefficients  0 (Equation ( 18)) and  (Equation ( 19)): Equation (20) shows the dispersion relation for the wave number : where   is the drag coefficient,  the shoot density,  the angular frequency, and ℎ the water depth.By using the damping factor provided by Losada et al. [31], we estimated   from Equations ( 17)- (19).Then, the wave height downstream was calculated by Equation (16).Finally, the effect of bending behavior could be explored by Equation (21): where  is the relative % decrease in wave height reduction considering the wave heights downstream from  _ and  _ due to the bending height ℎ  and vegetation height ℎ  , respectively.
Figure 11 shows the ratios of downstream wave heights and upstream wave heights (    ⁄ ) obtained in this study; the corresponding meadow lengths are 0.79-5.26m for case 1, 0.80-5.32m for case 2, 2.27-5.26m for case 3, and 0.72-5.23 m for case 4. From the figure, we can observe that the     ⁄ ratios calculated by ℎ  appear larger than those calculated by ℎ  for the four cases.This indicates that the wave height is reduced less when bending behavior is taken into account.This is because the bending behavior of vegetation reduces the frontal area, consequently reducing the energy lost in a wave travelling over vegetation [34,57].Therefore, calculating the wave height from the bending height decreases the wave height reduction compared with calculating it from the vegetation height.As the experiment of Losada et al. [31] pointed out, vegetation stems bend in the flow direction; hence, utilizing the bending height for wave height calculation reflects the actual motion of the vegetation stem, unlike when using the vegetation height.

Discussion
Luhar and Nepf [30] considered relatively extensive vegetation properties, including width (1 cm), thickness (0.45 mm), height (10 cm to 30 cm), and elastic modulus (0.4 GPa to 2.4 GPa).Thus, their results (three lines in Figure 10) correspond to the maximum (dotsolid line; upper LN11), minimum (dotted line; lower LN11), and middle values (solid line; LN11) of ℎ  ℎ  ⁄ , respectively.The experimental results of Z14 are near the lower limit of LN11.This is because Zeller et al. [33] used an elastic modulus of 0.3 GPa for a single artificial vegetation blade.In contrast, FK87 reached the upper limit of LN11.This is because the vegetation stems caused a decrease in velocity inside the meadow, leading to increased bending height, and the elastic modulus of live vegetation was relatively high in the work of Fonseca and Kenworthy [32].This point was also discussed by Luhar and Nepf [30].The regression equation, Equation (15), was derived from data with an elastic modulus of 0.5 GPa to 7.5 GPa, which is close to the upper limit data (2.4 GPa) of Luhar and Nepf [30]; hence, our results were also close to the upper limit.
As the vegetation meadow increases, the relative % decrease in wave height reduction () increases linearly, up to ~10% (9.23%) in the selected case studies (Figure 12).This implies that the impact of bending behavior becomes more pronounced with a larger vegetation meadow.In the figure,  is greater when the water depth is 0.60 m compared to 0.40 m.Considering the vegetation height of 0.588 m in our calculation, this indicates that the vegetation is fully submerged at a water depth of 0.60 m but partially submerged at 0.40 m; hence,  is greater in a fully submerged vegetation meadow than in a partially submerged vegetation meadow.In a fully submerged meadow, the current acts along the entire stem height, resulting in greater deformation or a smaller bending height.Calculating the wave height reduction using the vegetation height (assuming no bending occurs) implicitly implies that the vegetation is rigid.In contrast, when bending occurs, the vegetation can be considered flexible.Therefore, our results also indicate that rigid vegetation can reduce wave height ~10% more effectively than flexible vegetation.This argument aligns with the findings of other studies.For example, van Veelen et al. [34] reported that flexible vegetation attenuated waves by up to 70% less than rigid vegetation because flexible plants swayed, reducing the frontal area, the total work performed by the drag force, and thus the total wave energy lost when traveling over vegetation.Liu et al. [58] also noted that rigid vegetation was more effective in reducing wave heights compared to flexible vegetation in their laboratory-scale study.However, they did not specify the extent of such a reduction.Moreover, Mullarney and Henderson [59] found in field experiments that wave dissipation increased with the stiffness of vegetation stems.The overall wave dissipation in the flexible stems was approximately 30% of that observed in the rigid stems.
The work had certain limitations.First, we considered only a single type of vegetation and assumed that all vegetation stems bent identically in water flow.Vegetation on the front of a meadow is more likely to bend than vegetation behind the front; hence, this assumption may yield overestimates of the extents to which vegetation meadows reduce wave height.Second, during the FSI simulation, the water velocity was assumed to be uniform and temporally independent.In other words, the velocity remained constant throughout the vegetation meadow during the simulation.Third, although vegetation both bends and sways in water flow, we only focused on bending.

Conclusions
We considered nine types of marine vegetation and three water velocities, constructed 27 experimental cases, and obtained the bending heights via PIV measurements and FSI simulations.Of the 27 simulation results, 24 had errors < 5% in the experimental results; thus, the simulation was reliable.An additional 27 simulations were performed to provide a large dataset for multi-regression analysis (54 cases).To further validate the regression equation, an additional ten cases were considered.As a result, 63 of the 64 cases had errors < 10%, and 1 had an error of 10.71%.Moreover, a reasonable comparative result was achieved when we compared our regression equation with the existing lab-scale experiments and theoretical solutions.Thus, the multiple regression equation simply estimated the bending height with a tolerance of ~10%, showing the effectiveness of calibrating the vegetation height due to bending behavior.Based on these estimations, we predicted the effect of bending behavior on wave height.First, flexible vegetation was found to be up to ~10% less effective in reducing wave height than rigid vegetation.This is because bending reduced the frontal area of vegetation, the total work performed by the drag force, and the total energy lost by a wave traveling over the vegetation meadow.Second, the relative % decrease in wave height reduction increased as the vegetation meadow size increased, ranging from 1.08% to 9.23% in our investigation, particularly with a vegetation meadow of approximately 5 m.This implies that the impact of bending behavior becomes more pronounced with a larger vegetation meadow.Finally, we observed that the relative % decrease in wave height reduction was greater for fully submerged vegetation compared to partially submerged vegetation.

Figure 2 .
Figure 2. Mapping region for capturing a velocity profile: (a) 10 locations (0 L-27 L) vertically divided and 12 locations (0 L-33 L) horizontally divided, where L = 10.765mm, and (b) window (~290 mm × 355 mm) for measuring water velocity, showing the instantaneous velocity vector field obtained from the PIV measurement.

Figure 3 .
Figure 3. Independence tests for (a) domain length, (b) number of iteration steps, (c) analysis time, and (d) number of elements.Here, the relative % errors are with respect to the selected point; those errors, by definition, must all be positive; however, the "-" sign was intentionally added to highlight the differences with respect to the reference values.

Figure 4 .
Figure 4. Domain, boundary conditions, and mesh configuration: (a) domain and boundary conditions and (b) mesh configuration.

Figure 7 .
Figure 7.Comparison of bending heights obtained from the experiment (i.e., PIV measurement) and the simulation (i.e., FSI simulation).The error factor is defined by the ratio of bending height estimated by the FSI simulation (ℎ − ) to one obtained from the PIV measurement (ℎ − ).

Figure 8 .
Figure 8.Comparison of the bending heights estimated from the multiple regression, Equation (15), and obtained from the experiment and simulation (i.e., 27 experimental results (white circles), 27 simulation results (grey circles), and 10 additional simulation results (blue circles)).

Figure 9 .
Figure 9. Relative % errors of the bending heights estimated from multiple regression analysis.

Figure 10
Figure10compares the ratio (ℎ  ℎ  ⁄ ) of the present study with those of Fonseca and Kenworthy[32], Luhar and Nepf[30], and Zeller et al.[33].The data from Fonseca and Kenworthy[32], denoted by FK87 in the figure, were derived from a lab-scale experiment involving live vegetation exposed to unidirectional flow.The results from Luhar and Nepf[30], denoted by LN11, were obtained from a theoretical model predicting the bending height of a single artificial flexible vegetation stem under unidirectional flow.The data from Zeller et al.[33], denoted by Z14, were obtained from a lab-scale experiment using a blade model exposed to both wave and current conditions.

Figure 12 .
Figure 12.Relative % decrease in wave height reduction when we consider bending height.

Table 2 .
The characteristics of polystyrene used in this study.

Table 3 .
The characteristics of polystyrene used in this study.The density, Young's modulus, and thickness of the polystyrene are fixed as 1040 kg m -3 , 2.73 × 10 9 Pa, and 0.001 mm, respectively.

Table 4 .
Bending height (ℎ  ) of vegetation obtained from the PIV measurement and FSI simulation.

Table 5 .
Bending height of vegetation estimated by the FSI simulation considering 27 additionally designed simulation cases.