Hydrodynamic Performance Assessment of Emerged, Alternatively Submerged and Submerged Semicircular Breakwater: An Experimental and Computational Study

<jats:p>Coastal protection structures are essential defenses against wave energy, safeguarding coastal communities. This study aims to refine coastal protection strategies by employing a semicircular breakwater (SBW) model. Through a combination of physical and computational models, the hydrodynamic properties of the SBW under regular wave conditions were thoroughly examined. The primary objectives included delineating the hydrodynamic characteristics of SBWs, developing a computational model to validate experimental findings. Hydrodynamic characteristics of the SBW model were scrutinized across various wave conditions. Experimental testing in a wave flume covered a range of relative water depths (d/h) from 0.667 to 1.667, wave steepness (Hi/L) spanning 0.02 to 0.06 and wave periods ranging from 0.8 to 2.5 s. Notably, analysis of an emerged SBW with d/h = 0.667 revealed superior wave reflection, while an alternative submerged SBW with d/h = 1.000 showed the highest energy loss. These findings are further corroborated by the validation of computational models against experimental outcomes for d/h = 0.667, 1.000, 1.333 and 1.667. Moreover, the investigation of forces revealed an inverse correlation between horizontal forces and wave height, while vertical forces showed nuanced variations, including a slightly decreasing average vertical force with greater relative wave period (B/L) for different immersion scenarios.</jats:p>

section_content: Coastal protection is crucial for ensuring the safety and sustainability of communities and maritime activities. Breakwaters are essential structures designed to mitigate the impact of waves and currents, reduce wave heights and prevent wave-induced erosion [1]. They are built in various shapes and sizes using materials such as concrete, rocks or steel, depending on specific coastal conditions [2]. Three main types of breakwaters are commonly used: vertical, rubble mound and composite breakwaters [3]. Rubble mound breakwaters are the most effective in dissipating wave energy, while vertical breakwaters reflect the maximum wave energy but disperse only a limited amount [4]. Composite breakwaters use a combination of components to reflect and dissipate wave energy through a semi-circular caisson, providing an efficient and stable option for coastal protection [5]. Among the different types of breakwaters, semicircular breakwaters offer several advantages such as efficient wave energy dissipation, stability from the semi-circular caisson design, cost-effectiveness and a visually pleasing appearance [6] [7] [8]. Designing breakwaters is challenging due to harsh wave environments and poor ground conditions [9]. Despite the challenges, breakwaters are crucial for providing a safe and stable environment for maritime activities, while safeguarding coastal communities [10]. 

Mangroves, thriving in saline coastal environments, play a vital role in mitigating coastal erosion, protecting against storm surges and stabilizing shorelines through their intricate root systems [11]. Between 2000 and 2014, mangrove forests worldwide shrank from 137,760 km 2 to 81,484 km 2 due to human activity, despite their importance as green infrastructure [12, 13]. This loss significantly contributes to climate change, given mangroves' high carbon sequestration capacity [14]. Efforts to restore mangrove forests have been undertaken globally, driven by policies from government, non-government and nonprofit organizations [12, 15]. However, restoration projects often fall short of expectations and even fail due to a range of challenges, including physical and biological factors, emphasizing the need for innovative approaches to preserve mangroves while ensuring coastal protection [16] [17] [18]. 

Due to recent advancements in design and construction, several new model designs have been implemented [19, 20]. One concept is the SBW which is focused on by several researchers in China and Japan [6]. The SBW has a semi-circular caisson and bottom slab which is placed on a prepared rubble mound [21]. SBW has been categorized into four types: solid, permeable, front dissipation and rear dissipation [22]. Figure 1 illustrates these distinct SBW types. The SBW offers several advantages, including high stability in the presence of wave action due to its circular shape and the minimization of vertical forces on the foundation due to its hollow structure [23]. 

J. Mar. Sci. Eng. 2024, 12, 1105 2 of 31 semi-circular caisson, providing an efficient and stable option for coastal protection [5]. Among the different types of breakwaters, semicircular breakwaters offer several advantages such as efficient wave energy dissipation, stability from the semi-circular caisson design, cost-effectiveness and a visually pleasing appearance [6] [7] [8]. Designing breakwaters is challenging due to harsh wave environments and poor ground conditions [9]. Despite the challenges, breakwaters are crucial for providing a safe and stable environment for maritime activities, while safeguarding coastal communities [10]. Mangroves, thriving in saline coastal environments, play a vital role in mitigating coastal erosion, protecting against storm surges and stabilizing shorelines through their intricate root systems [11]. Between 2000 and 2014, mangrove forests worldwide shrank from 137,760 km 2 to 81,484 km 2 due to human activity, despite their importance as green infrastructure [12, 13]. This loss significantly contributes to climate change, given mangroves' high carbon sequestration capacity [14]. Efforts to restore mangrove forests have been undertaken globally, driven by policies from government, non-government and nonprofit organizations [12, 15]. However, restoration projects often fall short of expectations and even fail due to a range of challenges, including physical and biological factors, emphasizing the need for innovative approaches to preserve mangroves while ensuring coastal protection [16] [17] [18]. 

Due to recent advancements in design and construction, several new model designs have been implemented [19, 20]. One concept is the SBW which is focused on by several researchers in China and Japan [6]. The SBW has a semi-circular caisson and bottom slab which is placed on a prepared rubble mound [21]. SBW has been categorized into four types: solid, permeable, front dissipation and rear dissipation [22]. Figure 1 illustrates these distinct SBW types. The SBW offers several advantages, including high stability in the presence of wave action due to its circular shape and the minimization of vertical forces on the foundation due to its hollow structure [23]. The study of SBWs began in the 1980s with Tanimoto and his research team [24, 25]. Large-scale reinforced concrete structure SBWs with a semicircular caisson and a bottom slab made of precast were previously built at China's Yangtze Estuary, Japan's Miyazaki Port, Vietnam's Nha Mat Bac Place and China's Tianjin Port [26] [27] [28]. Prior research on coastal structures of this nature primarily focused on assessing hydrodynamic variables, including pressure, reflection, transmission and energy loss coefficients, wave diffraction, The study of SBWs began in the 1980s with Tanimoto and his research team [24, 25]. Large-scale reinforced concrete structure SBWs with a semicircular caisson and a bottom slab made of precast were previously built at China's Yangtze Estuary, Japan's Miyazaki Port, Vietnam's Nha Mat Bac Place and China's Tianjin Port [26] [27] [28]. Prior research on coastal structures of this nature primarily focused on assessing hydrodynamic variables, including pressure, reflection, transmission and energy loss coefficients, wave diffraction, rundown, runup and factors associated with energy dissipation by the SBW. Several semiempirical design formulas have been proposed [29] [30] [31]. 

Field tests have been conducted on SBWs constructed in the port of Miyazaki. The wave pressure data were reported to confirm a decrease in the horizontal force component, enhancing the structure's ability against sliding [32]. Prior research on SBWs featured finite element models, offering detailed analyses of structural behavior and performance [33]. To compute the pressures on the SBW, wave pressures were measured at various heights on the structure. The pressure decreased dynamically and exponentially from the water surface to the bed with increased wavelength and a lower immersion depth. Additionally, it was observed that Goda's modified approach accurately calculated the pressure magnitude, especially near the Stillwater level [34]. Another study found that the modified Goda method tends to underpredict the pressure at relative depth (d/L) < 0. 35 and overpredict the pressure on impermeable structures at d/L > 0. 58, where d is the water depth and L is the wavelength [35]. Investigating how differences in rubble mound height, water depth and perforation percentage affect the ability of semicircular breakwaters to convert non-breaking waves [36]. 

The interaction between the wave and the cylinder was studied using physical and computational modelling [37]. The equivalence of wave height on the porous cylinder to that of incident waves was observed. Investigations were carried out on models of emerged breakwaters with and without perforations. The findings revealed that an increase in the height and inclination of the incident waves resulted in a higher reflection coefficient [38]. 

Another study [39] examined the impact of submerged vertical and SBWs on local wave characteristics, specifically focusing on determining wave reflection coefficients. The proposed semi-empirical parameterizations rely on the dimensionless submergence parameter a/H (a-the depth of submergence of the breakwater; H i -the height of the incident wave at the breakwater). 

A recent study investigated the critical weight required to ensure sliding stability. Various parameters, including water depth and wave characteristics, were considered. A physical model of emerged and perforated (SBW) facing waves was utilized. The study also established a correlation between the hole diameter and spacing. Results showed that the stability parameter decreased with an increase in the slope of wave incident across all depth parameter ranges. Additionally, the study revealed a positive correlation between SBW stability and water depth, where an increase in the depth parameter led to a higher minimum weight required to maintain stability in sliding [40]. 

Analysis of this kind of structure using computational models has been made possible by the emergence of tools for Computational Fluid Dynamics (CFD) modelling, in tandem with the rise in processing capacity. Computational simulations can be utilized to investigate the interaction between breakwaters and fluids, thereby facilitating the assessment of wave loads on intricate geometric breakwaters. Flow-3D has been widely used in various fluid-structure interaction studies related to coastal structures, such as breakwater armour arrangements, concrete armoured blocks and floating breakwaters. According to their findings, the Volume of Fluids (VOF) algorithm-based computational method performed better than other computational codes [41]. Utilizing Flow-3D ® v11. 2 software, an assessment was conducted to evaluate the effectiveness of various shapes of concrete armored blocks in the context of breakwaters [42]. Conducting a study that combines Flow-3D software and a physical model to analyze the performance of a floating breakwater [43], researchers demonstrated the software's capability to accurately simulate hydrodynamic interaction effects with the analyzed structure. 

This research aims to identify the hydrodynamic characteristics of emerged, submerged and alternatively submerged small-scale SBWs and establish a computational model to authenticate the experimental findings. By substituting geotextile tubes with durable SBWs, the study addresses drawbacks such as short lifespan and susceptibility to punctures. SBWs provide a sustainable solution for safeguarding mangrove ecosystems and have wide-ranging applications in coastal protection. The insights obtained contribute to bolstering coastal resilience and promoting biodiversity conservation under regular wave conditions. The comprehensive model takes into account variables such as water depth and wave periods, aligning with prevalent conditions in regions like Malaysia. The paper is organized as follows: Section 1 outlines the objectives and background; Section 2 elaborates on the methodology, including the setup of the computational model for validation; Section 3 discusses an analytical error analysis; Section 4 presents the results and discussion, juxtaposing computational findings with experimental data to validate the model; and Section 5 summarizes the conclusions and key finding.
section_name: Introduction

section_content: The procedure described in this study involves two phases: the experimental phase and the computational phase. In the experimental phase, a physical model of the SBW is constructed and tested in a wave tank to measure its hydrodynamic responses to different wave conditions. The aim is to compute the transmission coefficient (C T ), reflection coefficient (C R ) and energy loss coefficient (C L ). In the computational phase, a computational model of the SBW is created using Flow-3D software and simulated under various wave conditions. The results of the physical tests and computational simulations are analyzed to evaluate the performance of the SBW in reducing the height of the waves.
section_name: Methodology

section_content: The study utilized an SBW model that was built with specific dimensions to evaluate its hydraulic characteristics under different wave conditions. The model had an external radius (R) of 0. 6 m, a thickness of 0. 1 m and a length of 0. 8 m. The width of the SBW (B) was 1. 2 m, while the height (h) was the same as R. To ensure that the model accurately matched the geometrical requirements of the test facility, a Froude scaling ratio of 1:2. 5 was employed. The model was constructed using concrete material with a density of 2400 kg/m 3, as shown in Figure 2. The experiment was conducted at the Offshore Laboratory in Universiti Teknologi PETRONAS. It featured a wave generator integrated into a flume, with a shallow basin with a maximum water depth of 1 m and a length of 20 m. Waves were generated and directed towards the SBW model, positioned 10 m away from the wave generator at one end of the flume. Six wave probes were strategically placed along the flume to measure water profiles. Offshore probes WP1, WP2 and WP3 captured reflected and incident waves, while shoreward probes WP4, WP5 and WP6 measured transmitted wave height. 

contribute to bolstering coastal resilience and promoting biodiversity conservation under regular wave conditions. The comprehensive model takes into account variables such as water depth and wave periods, aligning with prevalent conditions in regions like Malaysia. The paper is organized as follows: Section 1 outlines the objectives and background; Section 2 elaborates on the methodology, including the setup of the computational model for validation; Section 3 discusses an analytical error analysis; Section 4 presents the results and discussion, juxtaposing computational findings with experimental data to validate the model; and Section 5 summarizes the conclusions and key finding.
section_name: Description of Experimental Setup

section_content: The procedure described in this study involves two phases: the experimental phase and the computational phase. In the experimental phase, a physical model of the SBW is constructed and tested in a wave tank to measure its hydrodynamic responses to different wave conditions. The aim is to compute the transmission coefficient (CT), reflection coefficient (CR) and energy loss coefficient (CL). In the computational phase, a computational model of the SBW is created using Flow-3D software and simulated under various wave conditions. The results of the physical tests and computational simulations are analyzed to evaluate the performance of the SBW in reducing the height of the waves.
section_name: Methodology

section_content: The study utilized an SBW model that was built with specific dimensions to evaluate its hydraulic characteristics under different wave conditions. The model had an external radius (R) of 0. 6 m, a thickness of 0. 1 m and a length of 0. 8 m. The width of the SBW (B) was 1. 2 m, while the height (h) was the same as R. To ensure that the model accurately matched the geometrical requirements of the test facility, a Froude scaling ratio of 1:2. 5 was employed. The model was constructed using concrete material with a density of 2400 kg/m 3, as shown in Figure 2. The experiment was conducted at the Offshore Laboratory in Universiti Teknologi PETRONAS. It featured a wave generator integrated into a flume, with a shallow basin with a maximum water depth of 1 m and a length of 20 m. Waves were generated and directed towards the SBW model, positioned 10 m away from the wave generator at one end of the flume. Six wave probes were strategically placed along the flume to measure water profiles. Offshore probes WP1, WP2 and WP3 captured reflected and incident waves, while shoreward probes WP4, WP5 and WP6 measured transmitted wave height. In configuring the experimental setup, the lateral side of the SBW was meticulously sealed using silicon and cement to prevent any gaps. Additional trials confirmed the effectiveness of this sealing, ensuring no waves breached through the sealed gap and that the results remained unaffected. To mitigate reflections, the experimental setup employed a piston-type active absorbing wavemaker system, renowned for its absorption capability exceeding 96. 5% at target wave fields [44]. This system actively absorbs waves, minimizing spurious reflections from the wavemaker. Complementary passive absorption mechanisms prevent wave reflection from fixed boundaries within the flume [45]. Absorbing wavemakers offer dynamic cancellation of re-reflective waves by adjusting paddle motion, with a specially designed controller proving effective for both regular and irregular waves [46]. Hydrodynamic feedback signals are integrated into paddle control to eliminate re-reflection of waves generated by the wavemaker [47]. Furthermore, optimization of parameters for the absorbing beach in the wave tank significantly reduces reflection from the beach, ensuring accurate measurements [48]. 

The SBW model was tested with wave steepness values of 0. 02, 0. 04 and 0. 06, alongside varying relative water depth parameters (d/h) of 0. 667, 1. 000, 1. 333 and 1. 667 under regular wave conditions. Before the experiment, all wave probes underwent meticulous calibration. A wave absorber was strategically placed at the end of the flume to nullify incoming waves effectively. This structure minimized wave reflection from the end wall during experiments by featuring a 1:10 slope and comprised of right-angle triangular-shaped structures, each filled with filtering materials and covered by a steel grating to enhance efficiency. A total of 216 experiments were conducted using the least-squares approach of Mansard and Funke. 

Table 1 presents the experimental and computational parameters, including the wave steepness (H i /L), relative water depth (d/h, where d is water depth and h is SBW height), relative wave period (B/L where B is SBW width and L is wavelength) [49] and the total number of the experiments conducted. Figure 3 illustrates the laboratory setup through both a photograph and a schematic with components. Figure 4 presents photographs of the laboratory setup at various d/h ratios. Figure 5 delves into a comparative time series analysis of wave data at designated measurement points for various d/h ratios (a) d/h = 0. 667, (b) d/h = 1. 000, (c) d/h = 1. 333 and (d) d/h = 1. 667, with wave period 1. 6 and wave steepness 0. 02.
section_name: Description of Experimental Setup

section_content: The research utilized commercial computational fluid dynamics (CFD) software called Flow-3D to perform the computational runs. Flow-3D uses the Volume of Fluid (VOF) algorithm to model free surface flows [50] and wave structure interactions in coastal engineering [51]. VOF tracks the interface between fluids, which is governed by partial differential equations [52]. While CFD allows detailed simulations of coastal dynamics, its accuracy depends on data quality and field validation due to the complex nature of coastal environments.
section_name: Description of Computational Modeling Approach

section_name: Governing Equations

section_content: To simulate fluid flows in three-dimensional form, the mass continuity and Navier-Stokes momentum equations were used for incompressible fluids, within the Cartesian coordinate system [53]. The mass continuity equation is initially presented in Equation (1). However, for simulating incompressible fluid flow, this equation can be simplified, as shown in Equation ( 2 ). Meanwhile, Navier-Stokes equations were employed to solve the fluid velocity components (u, v, w) in the (x, y, z) directions. Equations ( 3 )-( 5 ) illustrate the general forms of the momentum equations utilized in Flow-3D for simulating incompressible fluid dynamics using Cartesian coordinates. 

∂u ∂t

∂v ∂t

∂w ∂t As mentioned earlier, the Flow-3D solver detected the flow free surface using the VOF function. The VOF technique has three main components: defining the VOF function, solving the VOF transport equation (Equations ( 6 ) and ( 7 )) and setting the boundary conditions at the free surface. The VOF method, based on the concept of a fractional volume of fluid, is more flexible and efficient for treating complicated free boundary configurations in numerical simulations [54]. In one-fluid simulations using Flow-3D, the volume-of-fluid (VOF) function (F) is essential for delineating fluid presence within computational cells. Void cells, indicated by (F = 0), represent regions devoid of fluid where variations in pressure, temperature, inertia and friction at the fluid interface are negligible, typically found in scenarios where gas density is significantly lower than that of the fluid and gas velocity matches fluid velocity. This approach efficiently models free surfaces in applications like mold filling with liquid metal, water flow in rivers and micro-fluidic devices, circumventing the computational expense of explicitly modeling gas flow [55]. (F) signifies fluid volume: (F = 1) for fully filled cells and (0 < F < 1) for partially filled cells, crucial for accurate simulation of complex boundaries and free surface phenomena [55, 56]. In this range, the variable (F) quantifies the fraction of the cell volume that is occupied by the fluid, allowing for a detailed representation of the varying degrees of fluid presence within the computational grid, as demonstrated in Figure 6. The fractional volume open to flow (VF) is influenced by fluid density (), veloci components (u, v, w) in specific coordinate directions, fractional area open to flow (Ax, A Az) in corresponding directions, density source term (RSOR), turbulence diffusion ter (RDIF) and coefficients (R and ξ) dependent on the coordinate system. When using cyli drical coordinates, y derivatives are converted to azimuthal derivatives, and ξ is set to For Cartesian coordinates, R is set to unity, and ξ is set to zero.
section_name: Governing Equations

section_content: To model turbulence in this study, the widely recognized two-equation k-ε mod was employed, which is among the most advanced equations commonly used for flu engineering [57] [58] [59]. This turbulence model utilizes two transport equations, one for tu bulent kinetic energy (k) and another for its dissipation (ε). The turbulent kinetic energ represents the energy associated with the turbulent variations in fluid flow, while the d sipation rate quantifies the rate at which this energy is dissipated as turbulence is atten ated. When using the K-ε turbulence model, the boundary condition on the SBW's surfa was set to be no slip, and the maximum turbulent mixing length was chosen to be dynam ically calculated using the algorithm of FLOW-3D software. By incorporating the transport equations, the study aimed to accurately simulate and analyze the wave chara teristics within the fluid flow surrounding the SBW under investigation. The two-equatio k-ε turbulence model is widely regarded as a sophisticated and extensively employe model for solving turbulent flows. Its versatility has been demonstrated in providing re sonable approximations for various types of flows [53].
section_name: ∂F ∂t

section_content: The computational model developed for simulating SBW hydrodynamics usin Flow-3D is illustrated in Figure 7. This figure is made up of five sub-figures that dep specific configurations of the breakwater at varying depths. They are (a) setup and mes (b) d/h = 0. 667 emerged, (c) d/h = 1. 000 alternately submerged, (d) d/h = 1. 333 and (e) d/h 1. 667 fully submerged SBWs. The computational mesh is carefully designed with varie cell sizes of 0. 04 m for the flume and 0. 01 m for the SBW, to align with the experiment parameters and balance resolution and computational efficiency as shown in Table SolidWorks -crafted SBW geometry in STL format is integrated into the model to ensu seamless representation. To optimize computational efficiency while preserving expe mental relevance, a 0. 04 m strip was extracted from the 0. 8 m width. At the left boundar the incident wave conditions were applied, while at the right boundary, an absorptio boundary was used to avoid wave reflection. Symmetry conditions were utilized for The fractional volume open to flow (V F ) is influenced by fluid density (ρ), velocity components (u, v, w) in specific coordinate directions, fractional area open to flow (A x, A y, A z ) in corresponding directions, density source term (R SOR ), turbulence diffusion term (R DIF ) and coefficients (R and ξ) dependent on the coordinate system. When using cylindrical coordinates, y derivatives are converted to azimuthal derivatives, and ξ is set to 1. For Cartesian coordinates, R is set to unity, and ξ is set to zero. 

To model turbulence in this study, the widely recognized two-equation k-ε model was employed, which is among the most advanced equations commonly used for fluid engineering [57] [58] [59]. This turbulence model utilizes two transport equations, one for turbulent kinetic energy (k) and another for its dissipation (ε). The turbulent kinetic energy represents the energy associated with the turbulent variations in fluid flow, while the dissipation rate quantifies the rate at which this energy is dissipated as turbulence is attenuated. When using the K-ε turbulence model, the boundary condition on the SBW's surface was set to be no slip, and the maximum turbulent mixing length was chosen to be dynamically calculated using the algorithm of Flow-3D software. By incorporating these transport equations, the study aimed to accurately simulate and analyze the wave characteristics within the fluid flow surrounding the SBW under investigation. The two-equation k-ε turbulence model is widely regarded as a sophisticated and extensively employed model for solving turbulent flows. Its versatility has been demonstrated in providing reasonable approximations for various types of flows [53].
section_name: Computational Model

section_content: The computational model developed for simulating SBW hydrodynamics using Flow-3D is illustrated in Figure 7. This figure is made up of five sub-figures that depict specific configurations of the breakwater at varying depths. They are (a) setup and mesh, (b) d/h = 0. 667 emerged, (c) d/h = 1. 000 alternately submerged, (d) d/h = 1. 333 and (e) d/h = 1. 667 fully submerged SBWs. The computational mesh is carefully designed with varied cell sizes of 0. 04 m for the flume and 0. 01 m for the SBW, to align with the experimental parameters and balance resolution and computational efficiency as shown in Table 2. SolidWorks-crafted SBW geometry in STL format is integrated into the model to ensure seamless representation. To optimize computational efficiency while preserving experimental relevance, a 0. 04 m strip was extracted from the 0. 8 m width. At the left boundary, the incident wave conditions were applied, while at the right boundary, an absorption boundary was used to avoid wave reflection. Symmetry conditions were utilized for all other open boundaries. Figure 8 shows the time step size used by the Flow-3D solver on selected mesh sizes for SBW simulation. other open boundaries. Figure 8 shows the time step size used by the FLOW3D solver on selected mesh sizes for SBW simulation. other open boundaries. Figure 8 shows the time step size used by the FLOW3D solver on selected mesh sizes for SBW simulation.
section_name: Computational Model

section_content: The assessment of SBW performance commonly involves the use of transmission coefficient (C T ), reflection coefficient (C R ) and energy dissipation coefficients (C L ) [60]. These coefficients are expressed mathematically as:

In the equations, H i, H T and H R represent the mean values of the incident, transmitted and reflected waves, respectively. Estimating energy dissipation at the breakwater poses a challenge due to its complexity in measurement [61] [62] [63] [64], influenced by factors like relative submerge depth, crest width and wave slope [61]. Consequently, the law of conservation of energy, as illustrated by Equation (10), is employed to approximate the amount of energy loss.
section_name: Performance Evaluation Measurements

section_content: The coefficient of determination (R 2 or r-squared) in regression models, calculated using Equation (11), is used to validate experimental results through computational analysis. A higher R 2 value indicates a stronger correlation between the datasets, suggesting a better fit of the two outcomes [65]. Alongside R 2, error functions like Mean Square Error (MSE), Root Mean Square Error (RMSE), Mean Absolute Deviation (MAD) and Mean Absolute Percentage Error (MAPE) (Equations ( 12 )-( 15 )) assess disparities between experimental and computational results [66] [67] [68]. These metrics play a critical role in evaluating computational model performance in validating experimental results, with smaller RMSE and MAPE values indicating closer agreement between experimental and computational values, signifying more accurate representation of experimental data by computational models. Thus, analyzing these error functions is essential for assessing computational model accuracy and reliability in validating experimental results. 

where A t represents experimental data, M t represents computational data, A X signifies the average of experimental data and n stands for the number of experiments conducted.
section_name: Analytical Error Analysis

section_content: The experimental data, analyzed and portrayed through wave transmission, reflection and energy loss graphs, undergoes parallel scrutiny with computational model results for validation. These graphical representations establish correlations between wave characteristics (C T, C R and C L ) and horizontal and vertical forces, incorporating crucial dimensionless parameters such as relative wave period (B/L), wave steepness (H i /L) and relative water depth (d/h). Notably, this investigation spans four relative water depths-d/h = 0. 667, 1. 000, 1. 333 and 1. 667-encompassing emerged, alternatively submerged and fully submerged configurations. In addition to wave transmission, reflection and energy loss graphs, the investigation includes analysis of wave forces on SBW, horizontal force analysis on SBW and vertical force analysis on SBW. The accuracy and reliability of the models are further evaluated through quantitative metrics, including MSE, RMSE, MAD, MAPE and R 2, offering a comprehensive understanding of model performance across diverse conditions. The forthcoming section will delve into the discussion of these mentioned graphs, providing detailed insights into the observed correlations and dynamics, as well as the effects of wave forces on the SBW structure, both horizontally and vertically. Figure 9 presents the free surface elevation comparison between the experimental (Exp) and computational (Comp) models at different d/h ratios: Figure 9a d
section_name: Results and Discussion

section_content: In the examination of wave transmission, both experimental and computational analyses are undertaken to evaluate the SBW's capacity to convey wave energy. This assessment is primarily conducted through the wave transmission coefficient (CT), which is de-
section_name: Wave Transmission

section_content: In the examination of wave transmission, both experimental and computational analyses are undertaken to evaluate the SBW's capacity to convey wave energy. This assessment is primarily conducted through the wave transmission coefficient (C T ), which is defined as the ratio of wave height transmitted (H T ) to the height of the incident wave (H i ) [69] [70] [71] [72]. A lower C T value corresponds to increased wave attenuation by the SBW [69]. Figure 10 displays C T values for d/h = 0. 667, 1. 000,1. 333 and 1. 667. In Figure 10a, the emerged SBW (d/h = 0. 667) shows a decline in the transmission coefficient (C T ) with increasing B/L, regardless of H i /L. This is consistent with previous studies which found that C T decreases as B/L increases in cases with vertical structures [73] and a lower arc-plate breakwater [74]. C T values at maximum and minimum, approximately 0. 32 and 0, indicate enhanced wave attenuation with increasing B/L, particularly effective for intercepting shorter period waves when B/L > 0. 4. The expanded SBW width facilitates more efficient interception and attenuation of waves with shorter wavelengths, carrying less energy and displaying less developed profiles. In Figure 10b, C T for the alternatively submerged SBW (d/h = 1. 000) declines with increasing B/L, irrespective of H i /L. Recorded maximum and minimum C T values are approximately 0. 85 and 0. 32, respectively. This suggests the alternatively submerged SBW's superior capability to intercept shorter-period waves during regular wave actions. Short waves, which convey lower energy flux and are less developed, undergo instability in the vicinity of the alternately submerged SBW. Certain waves have the potential to pass over, getting transmitted to the sheltered side of the structure with diminished wave height. In Figure 10c, for the fully submerged SBW with the relative water depth (d/h) of 1. 333, C T values exhibit relatively high values, ranging from 0. 59 to 0. 99, in comparison to those observed for the SBW with d/h ratios of 0. 667 and 1. 000. This implies that in comparison to the emerging and alternatively submerged SBW configurations, the submerged SBW exhibits less effective wave attenuation. The submerged SBW proves effective in intercepting longer period waves when B/L is less than 0. 6. With shorter period waves, the size of the circular water particle orbits diminishes as the water column extends in the z-direction [75]. The influence of water particle orbits on the SBW is nearly negligible [76, 77], enabling incident waves to pass through the submerged SBW with minimal flow interference. In Figure 10d, for the fully submerged SBW with d/h of 1. 667, C T values exhibit relatively high values ranging from 0. 68 to 0. 98 compared to other d/h values (0. 667, 1. 000 and 1. 333), and C T increases when B/L increases. In Figure 10, wave steepness significantly influences C T for both d/h = 0. 667, 1. 000, 1. 333 and 1. 667 [78]. Similar findings from studies on submerged SBWs and rectangular breakwaters support this observation. The weaker effects of the submerged breakwater on waves result in higher C T values in the numerical scheme, indicating that waves pass over the breakwater with minimal interaction [79]. C T values decrease with increasing H i /L for d/h = 0. 667, 1. 333 and 1. 667. In contrast, an inverse relationship is noted when (B/L) exceeds 0. 5 for (d/h) equal to 1. 000. During this phase, waves with higher steepness are partially intercepted by the SBW, and the remaining waves overtop the crest, forming transmitted waves of considerable sizes. Therefore, the C T of higher steepness waves exhibits larger values. 

Figures 11 12 13 validate the C T for the SBWs across various H i /L conditions. They illustrate C T at low (H i /L < 0. 02), moderate (0. 02 < H i /L < 0. 04) and high (0. 04 < H i /L < 0. 06) wave steepnesses for different d/h ratios (0. 667, 1. 000, 1. 333 and 1. 667). These figures visually confirm the model's accuracy in predicting C T, with computational results closely aligning with experimental data. Additionally, Table 3 presents a Comparative Analysis of Experimental and Computational Error Deviation Data for SBWs' C T, providing essential metrics like-and R 2 across various H i /L ratios and d/h scenarios. 

servation. The weaker effects of the submerged breakwater on waves result in higher CT values in the numerical scheme, indicating that waves pass over the breakwater with minimal interaction [79]. CT values decrease with increasing Hi/L for d/h = 0. 667, 1. 333 and 1. 667. In contrast, an inverse relationship is noted when (B/L) exceeds 0. 5 for (d/h) equal to 1. 000. During this phase, waves with higher steepness are partially intercepted by the SBW, and the remaining waves overtop the crest, forming transmitted waves of considerable sizes. Therefore, the CT of higher steepness waves exhibits larger values.
section_name: Wave Transmission

section_content: The assessment of the SBW's efficacy in reflecting wave energy is pivotal for gauging its effectiveness in coastal protection. This evaluation centers on the wave reflection coefficient, C R, which precisely quantifies the ratio of the reflected wave height to the incident wave height. This section delves into both experimental and computational investigations, providing comprehensive insights into the intricate dynamics governing wave reflection by SBW models. By meticulously examining various parameters such as wave steepness (H i /L), wave period and water depths, these experiments offer valuable insights into the nuanced behavior of the wave reflection coefficient, C R. The evaluation of the SBW's capacity to reflect wave energy revolves around the wave reflection coefficient, C R, which is described as the ratio of the height of the reflected wave (H R ) to the height of the incident wave (H i ) [75, 80]. Figure 14 meticulously scrutinizes the wave reflection characteristics of SBW models, examining the impact of H i /L and d/h ratios. It becomes apparent that the C R variation in relation to H i /L is nearly imperceptible for the fully submerged and alternatively submerged SBW when d/h = 1. 000, 1. 333 and 1. 667, indicating its insignificance in this context. In stark contrast, for d/h = 0. 667, wave steepness emerges as a significant parameter, with C R exhibiting an increase as H i /L decreases. Turning to specific d/h values, the C R for the emerged SBW (d/h = 0. 667) showcases an increase with increasing B/L, featuring relatively high values (ranging from 0. 46 to 0. 95) compared to d/h = 1. 000 and 1. 333. Concurrently, C T values exhibit an increase with increasing H i /L for d/h = 0. 667. The influence of B/L on C R is particularly noteworthy, especially at B/L less than 0. 4, with a sharp C R decline observed for alternatively and fully submerged SBW, reaching a minimum at 0. 3 < B/L < 0. 5. Subsequently, C R exhibits a small peak at B/L = 0. 4 before gradually declining for the larger B/L range, showcasing intriguing behavior consistent for d/h = 1. 000, 1. 333 and 1. 667. The bragging effect of C R, reported in large-scale SBWs [81, 82], underscores the complexity of the relationship. In summary, the proposed SBW demonstrates high reflectivity (0 < C R < 0. 7) against longer waves (B/L < 0. 3) for d/h = 1. 000, 1. 333 and 1. 667). A comparison of C R graphs in Figure 14a-d suggests that wave reflection by the emerged SBW is more dominant than the alternatively submerged SBW, likely due to its greater disruption of incident waves. This is likely due to the greater disruption of incident waves by the emerged structure. Similar findings from studies on submerged SBWs and rectangular breakwaters support this observation, showing that greater submersion results in less wave reflection [79]. The pronounced wave reflection characteristics of the SBW within the lower B/L range in Figure 14b -d result in diminished wave transmission, as illustrated in Figure 10. Ultimately, the proposed SBW emerges as a robust reflector for d/h = 0. 667 for all B/L values and only for B/L less than 0. 3 for d/h = 1. 000-1. 667, with an improvement of approximately 10% to 20%, the alternatively submerged SBW outperforms the submerged SBW in terms of reflective performance. rectangular breakwaters support this observation, showing that greater submersion results in less wave reflection [79]. The pronounced wave reflection characteristics of the SBW within the lower B/L range in Figure 14b -d result in diminished wave transmission, as illustrated in Figure 10. Ultimately, the proposed SBW emerges as a robust reflector for d/h = 0. 667 for all B/L values and only for B/L less than 0. 3 for d/h = 1. 000-1. 667, with an improvement of approximately 10% to 20%, the alternatively submerged SBW outperforms the submerged SBW in terms of reflective performance. The experimental validation of C R using computational analysis is presented, comparing across different Hi/L and various d/h ratios. illustrate this comparison at low, moderate and high wave steepnesses, respectively, showing consistency across d/h ratios. Results from both analyses demonstrate remarkable consistency, affirming the reliability and accuracy of the computational approach in predicting reflection coefficients. Additionally, Table 4 offers a detailed comparison of error deviation data for the C R between experimental and computational results across all H i /L and d/h ratios. The computational model consistently shows effective predictive capabilities, with error analyses like MSE, RMSE, MAD, MAPE and R 2 affirming its accuracy in capturing the reflective properties of SBW models across diverse conditions. 

The experimental validation of CR using computational analysis is presented, c paring across different Hi/L and various d/h ratios. Figures 15-17 illustrate this comp son at low, moderate and high wave steepnesses, respectively, showing consistency ac d/h ratios. Results from both analyses demonstrate remarkable consistency, affirming reliability and accuracy of the computational approach in predicting reflection co cients. Additionally, Table 4 offers a detailed comparison of error deviation data for CR between experimental and computational results across all Hi/L and d/h ratios. computational model consistently shows effective predictive capabilities, with error a yses like MSE, RMSE, MAD, MAPE and R 2 affirming its accuracy in capturing the re tive properties of SBW models across diverse conditions.
section_name: Wave Reflection

section_content: The examination of wave energy loss associated with SBW models constitutes a critical aspect of understanding their efficacy as coastal protection measures. This analysis provides a detailed exploration of both experimental and computational wave energy loss, experimental validation of energy loss coefficients using computational analysis, and error analysis for wave energy loss. By synthesizing insights from these analyses, a comprehensive understanding of the factors influencing wave energy dissipation and the effectiveness of SBWs in mitigating coastal hazards is provided. The effectiveness of the SBW is gauged by its wave reflection, transmission and energy dissipation attributes [5]. The wave dissipation coefficient, C L, serves as an indicator of hydraulic efficiency, derived from the energy conservation law

. The efficacy of SBW in mitigating wave energy is pivotal for assessing their viability as coastal defence solutions. SBWs, distinguished by their curvature, exhibit notably higher energy dissipation compared to conventional types [86]. Figure 18 elucidates the variation of C L across SBWs with different B/L and H i /L under regular wave conditions. Notably, emerged and fully submerged SBWs generally display lower wave energy attenuation compared to alternatively submerged SBWs, a phenomenon attributed to the pronounced curvature, particularly prominent for d/h = 1. 000. Interestingly, the alternatively submerged SBW with d/h = 1. 000 effectively intercepts incoming waves, leading to enhanced energy dissipation as waves traverse the crest. Figure 18b underscores this effect, showcasing a significant increase in C L with B/L, especially for shorter wavelengths, highlighting efficient energy release upon wave impact. Conversely, for d/h = 0. 667, 1. 333 and 1. 667, C L decreases with increasing B/L, demonstrating varied behavior dependent on relative water depth. The influence of H i /L on C L varies across different d/h and B/L values, with distinct trends observed for d/h = 0. 667, 1. 333 and 1. 667, primarily driven by wave breaking over the SBW. In summary, the alternatively submerged SBW (d/h = 1. 000) exhibits superior wave dissipation performance compared to fully submerged SBWs (d/h = 1. 667 and 1. 333), with emerged SBWs at d/h = 0. 667 showing similar efficacy. In fully submerged breakwaters, similar findings from studies on SBWs and rectangular breakwaters indicate that higher submergence leads to a lower dissipation coefficient. Conversely, less submerged breakwaters cause the dissipation coefficient to rise, reflecting increased wave interaction and energy dissipation [79]. 

Figures 19 20 21 provide a comparison between experimental and computational C L values across different H i /L categories and various d/h ratios (d/h = 0. 667, d/h = 1. 000, d/h = 1. 333 and d/h = 1. 667). These figures offer a comprehensive examination of C L under varying wave conditions. The comparison between experimental and computational results within each figure category enables a thorough assessment of the computational model's accuracy. Remarkably consistent results between experimental and computational analyses are observed, indicating the reliability and precision of the computational approach in predicting energy loss coefficients. Additionally, Table 5 offers a detailed comparison of error deviation data for C L between experimental and computational results across all H i /L and d/h ratios. The computational model consistently demonstrates a robust fit with experimental outcomes, showcasing effective predictive capabilities across various H i /L and d/h ratios. Metrics such as Mean Squared Error (MSE), Root Mean Squared Error (RMSE), Mean Absolute Deviation (MAD), Mean Absolute Percentage Error (MAPE) and R-squared (R 2 ) values affirm the accuracy of the model in capturing the energy dissipation properties of SBW models across the entire spectrum of H i /L and d/h ratios. 

uation compared to alternatively submerged SBWs, a phenomenon attributed to the pro nounced curvature, particularly prominent for d/h = 1. 000. Interestingly, the alternativel submerged SBW with d/h = 1. 000 effectively intercepts incoming waves, leading to en hanced energy dissipation as waves traverse the crest. Figure 18b underscores this effec showcasing a significant increase in CL with B/L, especially for shorter wavelengths, high lighting efficient energy release upon wave impact. Conversely, for d/h = 0. 667, 1. 333 an 1. 667, CL decreases with increasing B/L, demonstrating varied behavior dependent on re ative water depth. The influence of Hi/L on CL varies across different d/h and B/L values with distinct trends observed for d/h = 0. 667, 1. 333 and 1. 667, primarily driven by wav breaking over the SBW. In summary, the alternatively submerged SBW (d/h = 1. 000) ex hibits superior wave dissipation performance compared to fully submerged SBWs (d/h 1. 667 and 1. 333), with emerged SBWs at d/h = 0. 667 showing similar efficacy. In fully sub merged breakwaters, similar findings from studies on SBWs and rectangular breakwater indicate that higher submergence leads to a lower dissipation coefficient. Conversely, les submerged breakwaters cause the dissipation coefficient to rise, reflecting increased wav interaction and energy dissipation [79]. The contours depict the pressure distribution on the surface of the SBW, revealing that an increase in d/h results in a proportional escalation of pressure magnitude. This observation infers that larger waves exert more significant forces on the SBW. Notably, the pressure levels at the SBW's surface are comparatively lower than those encountered at greater depths, suggesting a gradual rise in pressure from the surface towards the depths of the SBW. This phenomenon is attributed to the additional hydrostatic pressure exerted by the water column above the SBWs surface. In Figures 23 and 24, the temporal evolution of horizontal (Fh) and vertical force (Fv) variations on the SBW is depicted under different d/h ratios (d/h = 0. 667, d/h = 1. 000, d/h = 1. 333, d/h = 1. 667). Figure 23 illustrates that Fh consistently fall within the specified range of -200 N to +200 N for all four cases on the interval from 20 to 40 s, providing a detailed examination of the SBW's horizontal forces. Simultaneously, Figure 24 presents the evolution of Fv, showcasing distinct ranges for each d/h ratio. For instance, Fv for d/h = 0. 667 (Figure 24b ) range from 360 N to 480 N, while for d/h = 1. 000 (Figure 24c ), the Fv vary from 660 N to 800 N, illustrating the dynamic forces exerted on the SBW as influenced by wave conditions and water depth. When d/h = 1. 333 (Figure 24d ), Fv exhibit a range of 680 N to 880 N, further emphasizing the influence of water depth on the forces exerted. For d/h = 1. 667 (Figure 24d ), Fv exhibit a range of 699 N to 791 N. These outcomes utilize a wave period of 1. 6 s and Hi/L of 0. 02. In Figures 23 and 24, the temporal evolution of horizontal (F h ) and vertical force (F v ) variations on the SBW is depicted under different d/h ratios (d/h = 0. 667, d/h = 1. 000, d/h = 1. 333, d/h = 1. 667). Figure 23 illustrates that F h consistently fall within the specified range of -200 N to +200 N for all four cases on the interval from 20 to 40 s, providing a detailed examination of the SBW's horizontal forces. Simultaneously, Figure 24 presents the evolution of F v, showcasing distinct ranges for each d/h ratio. For instance, F v for d/h = 0. 667 (Figure 24b ) range from 360 N to 480 N, while for d/h = 1. 000 (Figure 24c ), the F v vary from 660 N to 800 N, illustrating the dynamic forces exerted on the SBW as influenced by wave conditions and water depth. When d/h = 1. 333 (Figure 24d ), F v exhibit a range of 680 N to 880 N, further emphasizing the influence of water depth on the forces exerted. For d/h = 1. 667 (Figure 24d ), F v exhibit a range of 699 N to 791 N. These outcomes utilize a wave period of 1. 6 s and H i /L of 0. 02. lution of Fv, showcasing distinct ranges for each d/h ratio. For instance, Fv for d/h = 0. 667 (Figure 24b ) range from 360 N to 480 N, while for d/h = 1. 000 (Figure 24c ), the Fv vary from 660 N to 800 N, illustrating the dynamic forces exerted on the SBW as influenced by wave conditions and water depth. When d/h = 1. 333 (Figure 24d ), Fv exhibit a range of 680 N to 880 N, further emphasizing the influence of water depth on the forces exerted. For d/h = 1. 667 (Figure 24d ), Fv exhibit a range of 699 N to 791 N. These outcomes utilize a wave period of 1. 6 s and Hi/L of 0. 02. showing dynamic interactions for F h. t (ranging from 4 to 194 for d/h = 0. 667, 2 to 454 for d/h = 1. 000, 16 to 308 for d/h = 1. 333 and 14. 3 to 432. 2 for d/h = 1. 667), indicating the complex interplay of forces. Figure 26 illustrates F h. c for d/h = 0. 667, d/h = 1. 000, d/h = 1. 333 and d/h = 1. 667 conditions, revealing distinct relationships with B/L and H i /L. The observed ranges (ranging from 30 to 292 for d/h = 0. 667, 16 to 500 for d/h = 1. 000, 4 to 432 for d/h = 1. 333 and 0. 6 to 371 for d/h = 1. 667) signify varying magnitudes of forces due to wave crests, attributed to the dynamic interplay between wave characteristics and water depth.
section_name: Wave Energy Loss

section_content: In the vertical force analysis of the SBW, Figures 27 and 28 provide the vertical peak force due to the wave trough (F v. t ) and the average vertical peak force due to the wave crest (F v. c ), respectively. The study encompasses various H i /L ranges (H i /L < 0. 02, 0. 02 < H i /L < 0. 04, 0. 04 < H i /L < 0. 06) across different d/h conditions (0. 667, 1. 000, 1. 333 and 1. 667). In the emerged SBW (Figure 27a -d/h = 0. 667), F v. t ranges from 400 to 426, while for d/h = 1. 000 (Figure 27b ), d/h = 1. 333 (Figure 27c ) and d/h = 1. 667 (Figure 27d ) conditions, the ranges are 660 to 736, 648 to 772 and 702. 2 to 755. 2, respectively. These fluctuations highlight the influence of water depth on interaction dynamics, as increased B/L leads to reduced peak forces. In the emerged SBW (Figure 28a-d
section_name: Vertical Force Analysis on SBW

section_content: This study presents the development and evaluation of small SBW designed for diverse coastal protection applications. The hydraulic characteristics of the proposed SBW were investigated through a combined approach involving physical and computational modelling. Notably, the analysis of the CT revealed nuanced effectiveness across SBW, with the emerged SBW (d/h = 0. 667) exhibiting a smaller CT compared to d/h = 1. 000, d/h = 1. 333 and 1. 667. Crucially by correlating CT values with B/L ratios and comparing the transmitted wave height with the permissible wave height criteria for mangrove sampling survival, the success of the SBW at different B/L dimensions have been assessed. The SBW demonstrates a reflection accuracy of nearly 95% when d/h equals 0. 667. The alternatively submerged SBW (d/h = 1. 000) demonstrated superior wave attenuation performance, surpassing the fully submerged SBW (d/h = 1. 667) and achieving a substantial 70% reduction in incident wave height when exposed to regular waves with shorter periods. Moreover, the SBW with (d/h = 1. 000) displayed higher wave reflection., approximately 10% more than the fully submerged SBW when facing longer period waves. This reflection occurred notably at B/L less than 0. 4 for d/h = 1. 000-1. 667. The SBW, when alternatively submerged SBW and subjected to shorter waves, showcased superior energy dissipation performance, reaching up to 90% wave energy dissipation. However, it is crucial to acknowledge that wave energy dissipation tended to deteriorate as the immersion depth increased. Computational results aligned well with experimental findings, demonstrating the computational model's success in simulating t free surface flows in coastal engineering. Analyzing horizontal forces in the SBW, the average horizontal peak force due to wave crest (Fh. c) displayed variations across scenarios, with Fh. c being higher on submerged SBW compared to emerged SBW. Vertical forces in the SBW showed that Fv. t was lower in the emerged condition than in alternative and full submersion, indicating increased forces during wave troughs with deeper immersion. Fv. c was higher during wave crests in the emerged condition and for alternative and full submersion, with limited impact of Hi/L on forces.
section_name: Conclusions

section_content: This study presents the development and evaluation of small SBW designed for diverse coastal protection applications. The hydraulic characteristics of the proposed SBW were investigated through a combined approach involving physical and computational modelling. Notably, the analysis of the C T revealed nuanced effectiveness across SBW, with the emerged SBW (d/h = 0. 667) exhibiting a smaller C T compared to d/h = 1. 000, d/h = 1. 333 and 1. 667. Crucially by correlating C T values with B/L ratios and comparing the transmitted wave height with the permissible wave height criteria for mangrove sampling survival, the success of the SBW at different B/L dimensions have been assessed. The SBW demonstrates a reflection accuracy of nearly 95% when d/h equals 0. 667. The alternatively submerged SBW (d/h = 1. 000) demonstrated superior wave attenuation performance, surpassing the fully submerged SBW (d/h = 1. 667) and achieving a substantial 70% reduction in incident wave height when exposed to regular waves with shorter periods. Moreover, the SBW with (d/h = 1. 000) displayed higher wave reflection, approximately 10% more than the fully submerged SBW when facing longer period waves. This reflection occurred notably at B/L less than 0. 4 for d/h = 1. 000-1. 667. The SBW, when alternatively submerged SBW and subjected to shorter waves, showcased superior energy dissipation performance, reaching up to 90% wave energy dissipation. However, it is crucial to acknowledge that wave energy dissipation tended to deteriorate as the immersion depth increased. Computational results aligned well with experimental findings, demonstrating the computational model's success in simulating t free surface flows in coastal engineering. Analyzing horizontal forces in the SBW, the average horizontal peak force due to wave crest (F h. c ) displayed variations across scenarios, with F h. c being higher on submerged SBW compared to emerged SBW. Vertical forces in the SBW showed that F v. t was lower in the emerged condition than in alternative and full submersion, indicating increased forces during wave troughs with deeper immersion. F v. c was higher during wave crests in the emerged condition and for alternative and full submersion, with limited impact of H i /L on forces.
section_name: Conclusions

section_content: Acknowledgments: Special appreciation to Yayasan Universiti Teknologi PETRONAS (YUTP) and Ministry of Higher Education Malaysia (MOHE) for funding the research through grants 015LC0-352 and FRGS/1/2018/TK10/UTP/02/8, respectively.

section_name: Data Availability Statement:

section_content: The data that support the findings of this study are available on request from the corresponding authors.
section_name: Data Availability Statement: