Numerical Investigation of the Seabed Dynamic Response to a Perforated Semi-Circular Breakwater

: The semi-circular breakwater (SBW) has been implemented at numerous global locations due to its outstanding wave dissipation effectiveness and high structural performance. This study extends prior research by investigating the seabed dynamic response and hydrodynamic response characteristics around perforated SBWs. A coupled numerical model is developed to integrate waves, a semi-circular breakwater, and a sandy seabed. Wave behavior is simulated using Reynolds-averaged Navier–Stokes equations with a k-ε turbulence closure scheme, and the seabed response is numerically simulated using Biot’s full-dynamic (u-w) equations. After verifying computational accuracy, a series of tests is conducted to explore the effects of marine environments and SBW characteristics. Test results reveal a positive correlation between seabed response and wave height, wave period, and perforation number, while showing a negative correlation between seabed response and water depth and perforation rate. The basic perforation type is more effective than front and rear perforation types in maintaining a stable flow field and seabed response. These findings provide insights for designing SBWs for effective wave dissipation and seabed stability in complex marine environments, offering valuable recommendations for future designs.


Introduction
The semi-circular breakwater is a new structure that was initially developed and has been extensively studied in Japan in recent years [1].This structure demonstrates excellent effectiveness in wave dissipation, stability, and cost-effectiveness.It is suitable for deployment in environments characterized by significant waves and unstable seabed foundations [2,3].SBWs are classified into impermeable, front perforation (seaside wall perforated), rear perforation (leeside wall perforated), and basic perforation types (both seaside and leeside walls perforated) [4].In actual marine environments, waves and current flows coexist.The combined loading induces significant fluctuations in pore water pressure and effective stresses within the seabed [5][6][7].Therefore, it is crucial to investigate the seabed response characteristics around SBWs under wave-current interaction.
To enhance the understanding of the seabed response to wave action, extensive experiments and numerical studies have been conducted over the past decades.Mynett et al. [8] investigated the pore water pressure and stress distribution in a saturated seabed beneath vertical caissons using boundary layer theory.Tsai et al. [9] expanded this model to an infinite-thickness scenario and modified the breakwater from a caisson type to a composite type.Mase et al. [10] investigated the pore water pressure and effective stress distribution around perforated breakwaters using Biot's consolidation theory.They revealed the significant influence of the breakwater's perforation rate on pore water pressure and phase lag phenomena.Jeng et al. [11,12] proposed the GFEM-WSSI finite element model to simulate the pore water pressure distribution in the seabed beneath composite breakwaters.They conducted a comparative analysis of pore water pressure distribution and explored the impact of wave and seabed parameters on the seabed response.Zhang et al. [13] conducted numerical simulations to study the seabed response to perforated breakwaters, specifically discussing the influence of wave period and wave height on seabed pore water pressure.Zhang et al. [14] and Zhao et al. [15] extended a two-dimensional wave-structure-seabed model to a three-dimensional model.They numerically simulated the seabed dynamic response and investigated the effects of wave and seabed characteristics.Building upon the model proposed by Mizutani et al. [16], Mostafa et al. [17] developed a coupled (BEM-FEM) model to study wave-breakwater-seabed interaction.They conducted experimental studies on the deformation of caisson breakwaters and analyzed pore pressure distribution.Yan et al. [18] conducted model experiments to investigate the stability of a silty clay seabed around SBWs.They focused on the pore pressure around breakwaters and the seabed shear failure induced by the destruction of the soil fabric.
In addition, wave reflection can significantly impact the seabed dynamic response, potentially leading to severe coastal erosion [19].Xu et al. [20] conducted flume experiments to investigate the Bragg resonance reflection of wave propagation over submerged breakwaters.They discussed the influence of the permeability, relative width, relative height, and cross-sectional shapes of submerged breakwaters on Bragg resonance reflection.Brunone et al. [21] developed a rational definition of the relevant characteristics of wave-structure interaction.They proposed a second-order model to describe the flow field characteristics on a steep slope.
This study uses a coupled numerical model to investigate the dynamic response of a sandy seabed.Reynolds-averaged Navier-Stokes equations (RANS) are used to model wave-current behavior.Meanwhile, 3D full-dynamic (u-w) governing equations based on Biot's poro-elastic theory are adopted for modeling the seabed response.Due to the inadequacy of the existing research on the seabed response surrounding SBWs, this study investigates the influence of marine environments and SBW characteristics on the seabed dynamic response.
The rest of this paper is organized as follows: Section 2 introduces both the flow model and the seabed model.Model validation and numerical setup are presented in Section 3. Section 4 discusses the effects of the coupled wave-seabed-structure model under different wave conditions and breakwater structures.Section 5 provides the conclusion and offers recommendations.

Flow Model
Numerical simulations are based on Reynolds-averaged Navier-Stokes equations (RANS), which describe the movement of a wave-induced fluid.
In this study, solid geometry is represented using area and volume fractions on Cartesian meshes to enhance calculation efficiency and accuracy.The model equations are as follows: where x, y, and z are longitudinal, vertical, and spanwise directions, respectively; v x , v y , and v z are longitudinal, vertical, and spanwise velocities, respectively; V F is the flowable volume fraction in FAVOR grid processing; A x , A y , and A z are the fractions of flowable area in x, y, and z directions, respectively, in grid processing; ρ f is the density of the fluid; p is the pressure; g is the gravitational acceleration; and f x , f y , and f z are the acceleration of viscous forces in the x, y, and z directions, respectively.The RNG k-ε turbulence scheme proposed by Yakhot et al. [22,23] is used in this study.The model is derived from the transient Navier-Stokes equation using a mathematical method known as the renormalization group (RNG).The RNG k-ε model significantly enhances the simulation accuracy of fluid mechanics by incorporating the effects of largescale movements and adjusting viscosity to account for small-scale effects.It systematically eliminates small-scale movements from the control equation, thereby improving the simulation performance for transient flow and flow with bending characteristics.The prognostic equations of turbulent kinetic energy and its dissipation rate can be written as: where C 1 , C 2 , and C 3 are empirical constants; k T is the turbulent kinetic energy; ε T is the turbulence dissipation; P T is the velocity-gradient-induced turbulence kinetic energy generation term; G T is the buoyancy production term; and D k T and D ε are the diffusive terms.This study uses the Euler method for free-surface tracking, incorporating the marker and cell (MAC) and the volume-of-fluid (VOF) methods.Fluid configurations are defined in terms of a VOF function F (x, y, z, t) [24], which represents the volume of fluid occupying a portion of unitary volume space.The fluid fraction F is a variable that depends on position coordinates and time, representing the ratio of fluid volume within the calculation unit to the total volume of the unit.The F value is 1 when the unit is full of fluid, 0 when there is no fluid in the unit, and between 0 and 1 when the unit contains free surfaces or bubbles.The fluid component F satisfies the transport equation: The position of the free surface is identified within units, where the F value ranges between 0 and 1.Since the fluid fraction F is a step function, numerical oscillations can easily occur at discontinuity points when using a differential scheme for a discrete solution, potentially compromising the original definition of the fluid fraction F. To mitigate this issue, Hirt and Nichols [24] proposed the VOF method and the donor-acceptor method for reconstructing the free surface.

Seabed Model
In numerical simulations, the seabed model is primarily solved using 3D full-dynamic (u-w) governing equations based on Biot's poro-elastic theory [25,26].Additionally, the model can be switched to the partly dynamic (u-p) mode or the quasi-static (Biot's consolidation) mode as an optional configuration.The model is capable of simulating a waveinduced seabed response, wave-current interaction, and wave-structure interaction [27,28].The model incorporates the inertia effects of both soil and pore fluid and accounts for the 3D fully non-homogeneous properties of the seabed and cross-anisotropic soil behavior.A notable advantage of the model is its capability to handle arbitrary wave-seabed-structure configurations through the material identification technique.The model can simulate a wave-induced seabed response with detailed physical resolutions, including pore pressure, effective stresses, displacements of soil and pore fluid, and seepage flow.
This study formulates three governing equations in tensor form for a problem involving wave-seabed-structure interaction, assuming the seabed is a porous elastic medium.These equations include the overall equilibrium equation of the soil, the equilibrium equation of pore fluid flow, and the mass balance equation: . .
where σ is the total stress, ρ is the total density of the porous medium, ρ f is the density of the fluid, g is the gravitational acceleration, u i is the displacement of the soil matrix, w i is the average relative displacement of the fluid to the solid skeleton, k i is the permeability of the porous medium, n e is the porosity of the solid phase, K f is the elastic modulus of the pore fluid, and ε is the strain of the soil skeleton.
In a natural non-homogeneous seabed, the permeability (K) and Young's modulus (E) vary spatially.These parameters are spatially dependent on longitudinal, vertical, and spanwise coordinates in the model.This model considers the non-uniform distribution (R(x, y, z)) of eight parameters: the density of the seabed (ρ s ), the density of the fluid (ρ f ), soil permeability (K), the porosity of the solid phase (n), the elastic modulus of the pore fluid (K f ), Young's modulus (E), vertical Poisson's ratio (u hv ), and lateral Poisson's ratio (u hh ).The model is newly formulated to incorporate the spatial gradient of these non-uniform seabed parameters, resulting in the fully dynamic equations (Equations ( 11) and ( 12)): In this study, the flow model is used to analyze the hydrodynamic response characteristics around SBWs.The pressure data obtained from the analysis are integrated into the seabed model to determine the seabed response.Additionally, a comprehensive dynamic response model is used to evaluate the variation in liquefaction depth.

Model Validation and Numerical Setup
In this section, we use experimental studies by Mizutani et al. [29] and Wu and Hsiao [30] to validate numerical simulation results.These experiments specifically focus on permeable submerged breakwaters, making them valuable validations for the numerical model used in this study.

Model Validation
The experiments conducted by Mizutani et al. [29] investigated the interaction between regular waves, a sandy seabed, and a perforated SBW. Figure 1 shows a sketch of the numerical wave flume.Four gauges (WG1: x = −1.93m; WG2: x = 0.42 m; WG3: x = 1.47 m; WG4: x = 3.13 m) are set up to measure the temporal evolution of the wave surface.The test conditions for the model validation are provided in Table 1.The simulated wave height is 0.03 m, and the wave period is 1.4 s.The time ranges from 0 s to 60 s, with intervals of 1.4 s.The calculation domain is 15 m in length and 0.5 m in height, with a mesh resolution of 0.02 m.The mesh resolution in the breakwater zone is 0.01 m.
horizontal and vertical velocities along the water depth at various positions (x = 0.00 m, x = 0.08 m, x = 0.16 m).
In conclusion, the numerical results of velocity agree well with the experimental data, demonstrating that this model accurately predicts fluid velocity around the breakwater.
In addition, the numerical results of the hydrodynamic response and seabed response characteristics agree well with the experimental results, demonstrating the accuracy of the model in predicting wave propagation over a perforated breakwater.Figures 2 and 3 illustrate the variation in dimensionless pressure p/γH and dimensionless wave elevation η/H (η, wave elevation; H, incident wave height; γ, specific of water) with t/T. Figure 2 shows the temporal evolution of pore pressure, comparing simulation and experimental results.Figure 2 demonstrates that the pore pressure simulation is more accurate in P1 and P3, while the mathematical simulation results in P3 exhibit a smoother trend.As shown in Figure 3, the temporal evolution of the wave elevation at WG1/WG2 demonstrates a general agreement with the experiment.However, the correspondence of the temporal evolution at WG3/WG4 is slightly weaker.This occurs because the numerical model uses a sponge layer to absorb transmitted waves, weakening the reflected waves in the leeside.Consequently, the simulation results are relatively smooth.
The experiment conducted by Wu and Hsiao [30] is used to verify the velocity distribution along the wave depth.The rectangular permeable breakwater is 13 cm long and 6.5 cm high, with a porosity value of 0.52.The calculation domain is 15 m in length and 0.25 m in height.The mesh resolution in the breakwater zone is 0.01 m.The origin of the coordinate system is defined at the intersection of the left side of the breakwater and the bottom of the flume.The water depth is 10.6 cm, and the wave height is 4.77 cm.The time ranges from 0 s to 10 s, with intervals of 0.2 s. Figure 4 presents the comparison results of horizontal and vertical velocities along the water depth at various positions (x = 0.00 m, x = 0.08 m, x = 0.16 m).In conclusion, the numerical results of velocity agree well with the experimental data, demonstrating that this model accurately predicts fluid velocity around the breakwater.
In addition, the numerical results of the hydrodynamic response and seabed response characteristics agree well with the experimental results, demonstrating the accuracy of the model in predicting wave propagation over a perforated breakwater.

Numerical Setup
The sketch of the numerical layout is depicted in Figure 5, where the scale is 1:25.Figure 6 illustrates the characteristics of perforated SBWs.Previous studies have shown that dividing the waveform into 20 grids within the wave height range provides a more accurate representation [31,32].To discretize the domain, a total of 800,000 cells are used.The first block comprises 500,000 cells, with dimensions of 0.03 m, while the second block comprises 300,000 cells, with dimensions of 0.01 m.Initial simulations are conducted using various grid resolutions to achieve an adequate solution in terms of capturing the relevant flow details. evant flow details.
The wave and the current are generated at the inlet boundary.The renormalized group (RNG) turbulence model is used to represent turbulent flow phenomena.To reduce wave reflection from the boundary, a wave-absorbing layer is installed ahead of the outflow boundary.The bottom of the domain is designated as a wall boundary type, while the other boundaries are set as symmetry types.The time ranges from 0 s to 20 s, with intervals of 0.0475 s.The wave and the current are generated at the inlet boundary.The renormalized group (RNG) turbulence model is used to represent turbulent flow phenomena.To reduce wave reflection from the boundary, a wave-absorbing layer is installed ahead of the outflow boundary.The bottom of the domain is designated as a wall boundary type, while the other boundaries are set as symmetry types.The time ranges from 0 s to 20 s, with intervals of 0.0475 s.
The type of incident wave in cases of wave-current interaction is determined by the wave height (H), period (T), water depth (d), and flow velocity (c).In the nearshore of West Africa, the interaction between waves and currents is particularly significant, intensifying the dynamic response of the seabed.This study uses hydrographic data from the fishing port in Ghana, West Africa, as presented in Table 2.The liquefaction criteria [33] used in this study are outlined as follows: where p b is the excess pore water pressure (wave dynamic pressure), (p − p b ) is the excess pore water pressure difference, and σ z0 is the initial effective stress of the soil.

Results and Discussion
In this section, the effects of marine environments and breakwater structures on the seabed dynamic response are investigated.

Impact of the Marine Environment
In the ocean, the marine environment is diverse.This section analyzes the impact of the wave height, water depth, and wave period on the seabed response.The perforation rate is 10%, the perforation number is 7, and the perforation type is basic perforation.The liquefaction depths of the seabed in different marine environments are detailed in Table 3. Figure 7 shows the maximum pore pressure and liquefaction depth for different wave heights.As the wave height increases, the maximum pore pressure at the same seabed depth also increases, especially within the range of −0.3 < z/h < 0. This increase is due to the intensified hydrodynamic forces acting on the seabed, resulting in higher pore pressure.Within the range of −1 < z/h < −0.3, there is a noticeable reduction in the rate of pore pressure dissipation, causing a gradual decrease in pore pressure toward 0 Pa.This phenomenon explains why the seabed in shallow regions is prone to experiencing transient liquefaction.
to the intensified hydrodynamic forces acting on the seabed, resulting in higher pore pressure.Within the range of −1 < z/h < −0.3, there is a noticeable reduction in the rate of pore pressure dissipation, causing a gradual decrease in pore pressure toward 0 Pa.This phenomenon explains why the seabed in shallow regions is prone to experiencing transient liquefaction.
Figure 7b illustrates that the range of seabed liquefaction is within 0.05 m and correlates positively with changes in wave height.Due to wave reflection by the SBW, the maximum liquefaction depth typically occurs in the seaside [18].With an increase in wave height from 0.06 m to 0.10 m, the maximum liquefaction depth increases from 1.22 cm to 3.02 cm.  Figure 7b illustrates that the range of seabed liquefaction is within 0.05 m and correlates positively with changes in wave height.Due to wave reflection by the SBW, the maximum liquefaction depth typically occurs in the seaside [18].With an increase in wave height from 0.06 m to 0.10 m, the maximum liquefaction depth increases from 1.22 cm to 3.02 cm.

Wave Period
Figure 9 illustrates the seabed response to different wave periods.As the wave period increases, the wave energy strengthens, resulting in a slight increase in the pore pressure around the SBW.Moreover, a notable increase in liquefaction depth occurs.This phenomenon is caused by the increased wave transmission associated with longer periods, intensifying wave action in the leeside, consistent with previous studies [34].Specifically, leeside liquefaction increases by 0.39 cm, with the average liquefaction depth increasing by 0.4 cm.The analysis demonstrates that pore pressure increases significantly, leading to seabed liquefaction under the influence of long-period waves.This poses a substantial threat to the stability of the seabed.The liquefaction characteristic values of the seabed in different marine environments can be found in Table 3.

Wave Period
Figure 9 illustrates the seabed response to different wave periods.As the wave period increases, the wave energy strengthens, resulting in a slight increase in the pore pressure around the SBW.Moreover, a notable increase in liquefaction depth occurs.This phenomenon is caused by the increased wave transmission associated with longer periods, intensifying wave action in the leeside, consistent with previous studies [34].Specifically, leeside liquefaction increases by 0.39 cm, with the average liquefaction depth increasing by 0.4 cm.The analysis demonstrates that pore pressure increases significantly, leading to seabed liquefaction under the influence of long-period waves.This poses a substantial threat to the stability of the seabed.The liquefaction characteristic values of the seabed in different marine environments can be found in Table 3.
sifying wave action in the leeside, consistent with previous studies [34].Specifically, leeside liquefaction increases by 0.39 cm, with the average liquefaction depth increasing by 0.4 cm.The analysis demonstrates that pore pressure increases significantly, leading to seabed liquefaction under the influence of long-period waves.This poses a substantial threat to the stability of the seabed.The liquefaction characteristic values of the seabed in different marine environments can be found in Table 3.

Impact of Breakwater Structure
In this section, the effects of the perforation rate, perforation number, and perforation type on hydrodynamic response characteristics and the seabed dynamic response are analyzed.The flow velocity is 0.08 m/s, the wave height is 0.088 m, the water depth is 0.228 m, and the wave period is 1.9 s.
Figure 10 depicts the velocity contour in the flow field at different time points: 9.07 s, 9.30 s, 10.17 s, and 10.64 s.When the wave crest is above the SBW (t = 9.07 s), significant wave transformation occurs above the breakwater, leading to the maximum velocity.As the wave traverses the breakwater at 9.30 s, the flow velocity surrounding the structure decreases, resulting in the expansion of the core vortex.An inverse head is observed, with the water level in the leeside surpassing that in the seaside.The development of a reverse flow leads to intensified water mixing and the formation of multiple high-intensity vortices by 10.17 s.By 10.64 s, the approaching wave causes incident and recirculating waters

Impact of Breakwater Structure
In this section, the effects of the perforation rate, perforation number, and perforation type on hydrodynamic response characteristics and the seabed dynamic response are analyzed.The flow velocity is 0.08 m/s, the wave height is 0.088 m, the water depth is 0.228 m, and the wave period is 1.9 s.
Figure 10 depicts the velocity contour in the flow field at different time points: 9.07 s, 9.30 s, 10.17 s, and 10.64 s.When the wave crest is above the SBW (t = 9.07 s), significant wave transformation occurs above the breakwater, leading to the maximum velocity.As the wave traverses the breakwater at 9.30 s, the flow velocity surrounding the structure decreases, resulting in the expansion of the core vortex.An inverse head is observed, with the water level in the leeside surpassing that in the seaside.The development of a reverse flow leads to intensified water mixing and the formation of multiple high-intensity vortices by 10.17 s.By 10.64 s, the approaching wave causes incident and recirculating waters to mix within the arch, contributing to wave energy dissipation.This process illustrates the primary mechanism for energy dissipation in the SBW. to mix within the arch, contributing to wave energy dissipation.This process illustrates the primary mechanism for energy dissipation in the SBW.

Breakwater Perforation Rate
Figure 11 compares the velocity contours at the wave crest and the wave trough at different perforation rates.The perforated SBW creates a stilling chamber due to its structure.As the perforation rate increases, the flow velocity decreases above the breakwater, while increasing inside the breakwater.This occurs because the effective water-blocking area decreases significantly, allowing water to enter and exit more easily through the structure.At a perforation rate of 15%, excessive perforations reduce disturbance levels, leading to a notable reduction in vortices and the turbulence intensity within the arch.As the reverse flow develops, the increasing perforation rate reduces the extent of the reverse head due to expanded water flow paths.This diminishes the range of vortices within the arch, weakening their impact on the seabed.Additionally, perforations with the minimum internal cross-section area consistently exhibit higher flow velocities in the region.

Breakwater Perforation Rate
Figure 11 compares the velocity contours at the wave crest and the wave trough at different perforation rates.The perforated SBW creates a stilling chamber due to its structure.As the perforation rate increases, the flow velocity decreases above the breakwater, while increasing inside the breakwater.This occurs because the effective water-blocking area decreases significantly, allowing water to enter and exit more easily through the structure.At a perforation rate of 15%, excessive perforations reduce disturbance levels, leading to a notable reduction in vortices and the turbulence intensity within the arch.As the reverse flow develops, the increasing perforation rate reduces the extent of the reverse head due to expanded water flow paths.This diminishes the range of vortices within the arch, weakening their impact on the seabed.Additionally, perforations with the minimum internal cross-section area consistently exhibit higher flow velocities in the region.

Breakwater Perforation Rate
Figure 11 compares the velocity contours at the wave crest and the wave trough at different perforation rates.The perforated SBW creates a stilling chamber due to its structure.As the perforation rate increases, the flow velocity decreases above the breakwater, while increasing inside the breakwater.This occurs because the effective water-blocking area decreases significantly, allowing water to enter and exit more easily through the structure.At a perforation rate of 15%, excessive perforations reduce disturbance levels, leading to a notable reduction in vortices and the turbulence intensity within the arch.As the reverse flow develops, the increasing perforation rate reduces the extent of the reverse head due to expanded water flow paths.This diminishes the range of vortices within the arch, weakening their impact on the seabed.Additionally, perforations with the minimum internal cross-section area consistently exhibit higher flow velocities in the region.Figure 12 illustrates that increasing the perforation rate results in a reduced liquefaction depth.Increasing the perforation rate from 5% to 10% reduces the turbulent intensity Figure 12 illustrates that increasing the perforation rate results in a reduced liquefaction depth.Increasing the perforation rate from 5% to 10% reduces the turbulent intensity within the arch, resulting in a 17.3% decrease in the average seabed liquefaction depth.A perforation rate of 15% results in minimal liquefaction around the SBW, with a 23.2% decrease in the liquefaction depth in the seaside.These findings are consistent with previous studies indicating that wave transmission increases with higher perforation rates [4,34,35].This study proposes the existence of an optimal perforation rate that maximizes the wave dissipation capacity, while ensuring seabed stability.
within the arch, resulting in a 17.3% decrease in the average seabed liquefaction depth.A perforation rate of 15% results in minimal liquefaction around the SBW, with a 23.2% decrease in the liquefaction depth in the seaside.These findings are consistent with previous studies indicating that wave transmission increases with higher perforation rates [4,34,35].This study proposes the existence of an optimal perforation rate that maximizes the wave dissipation capacity, while ensuring seabed stability.

Breakwater Perforation Number
Figure 13 depicts the velocity contours of the flow fields around the SBW with different perforation numbers.Initially, as the perforation number increases, there is a decrease in flow velocity, followed by an increase.This phenomenon occurs because the increased number of pathways for water flow initially reduces flow velocity.Simultaneously, a significant velocity gradient within the arch enhances water mixing due to the influence of reverse flow.However, as the perforation number continues to increase, the decrease in the perforation rate limits the water drainage capacity.Overflowing water over the top of the breakwater not only increases flow velocity but also strengthens the inverse head, leading to expanded vortices within the arch.

Breakwater Perforation Number
Figure 13 depicts the velocity contours of the flow fields around the SBW with different perforation numbers.Initially, as the perforation number increases, there is a decrease in flow velocity, followed by an increase.This phenomenon occurs because the increased number of pathways for water flow initially reduces flow velocity.Simultaneously, a significant velocity gradient within the arch enhances water mixing due to the influence of reverse flow.However, as the perforation number continues to increase, the decrease in the perforation rate limits the water drainage capacity.Overflowing water over the top of the breakwater not only increases flow velocity but also strengthens the inverse head, leading to expanded vortices within the arch.
Figure 13 depicts the velocity contours of the flow fields around the SBW with different perforation numbers.Initially, as the perforation number increases, there is a decrease in flow velocity, followed by an increase.This phenomenon occurs because the increased number of pathways for water flow initially reduces flow velocity.Simultaneously, a significant velocity gradient within the arch enhances water mixing due to the influence of reverse flow.However, as the perforation number continues to increase, the decrease in the perforation rate limits the water drainage capacity.Overflowing water over the top of the breakwater not only increases flow velocity but also strengthens the inverse head, leading to expanded vortices within the arch.Figure 14 illustrates the increase in pore pressure with an increasing perforation number.However, the influence of the perforation number on seabed liquefaction is relatively minor.An increase in the perforation number results in an approximate 8% increase in liquefaction in the seaside, along with a decrease in liquefaction in the leeside.When Figure 14 illustrates the increase in pore pressure with an increasing perforation number.However, the influence of the perforation number on seabed liquefaction is relatively minor.An increase in the perforation number results in an approximate 8% increase in liquefaction in the seaside, along with a decrease in liquefaction in the leeside.When the perforation number is 3, the reduced presence of vortices significantly mitigates the seabed liquefaction.
the perforation number is 3, the reduced presence of vortices significantly mitigates the seabed liquefaction.

Breakwater Perforation Type
Figure 15 illustrates the velocity contours of flow fields for different perforation types.On the one hand, basic perforation allows water to pass through more easily, leading to a more stable flow field and improved seabed stability, as demonstrated in Figure 16.Moreover, basic perforation exhibits strong wave transmission properties and is less susceptible to reverse water flow [4].This reduction significantly minimizes flow field disturbances and increases structural stresses.Sasajima et al. [36] noted that the wave dissipation effect of basic perforation is relatively weak due to the low intensity of wave dissipation vortices.
On the other hand, front and rear perforations rely on one-sided drainage, affecting flow uniformity and water pressure gradients.This influences flow field stability and exacerbates local seabed responses (as shown in Figure 16).The seaside liquefaction for the front perforation increases by 22.9% and the leeside liquefaction for the rear perforation increases by 45.9% compared to the basic perforation.The increased local flow velocity leads to increased disturbances on the seabed.Thus, it is important to account for these factors when designing structures to withstand such impacts and disturbances, ensuring the stability and safety of SBWs.The liquefaction characteristic values of different breakwater structures are detailed in Table 4.

Breakwater Perforation Type
Figure 15 illustrates the velocity contours of flow fields for different perforation types.On the one hand, basic perforation allows water to pass through more easily, leading to a more stable flow field and improved seabed stability, as demonstrated in Figure 16.Moreover, basic perforation exhibits strong wave transmission properties and is less susceptible to reverse water flow [4].This reduction significantly minimizes flow field disturbances and increases structural stresses.Sasajima et al. [36] noted that the wave dissipation effect of basic perforation is relatively weak due to the low intensity of wave dissipation vortices.
acerbates local seabed responses (as shown in Figure 16).The seaside liquefaction for the front perforation increases by 22.9% and the leeside liquefaction for the rear perforation increases by 45.9% compared to the basic perforation.The increased local flow velocity leads to increased disturbances on the seabed.Thus, it is important to account for these factors when designing structures to withstand such impacts and disturbances, ensuring the stability and safety of SBWs.The liquefaction characteristic values of different breakwater structures are detailed in Table 4.

Conclusions
This study numerically investigates the effects of marine environments and breakwater structures on seabed dynamic responses.Specifically, it investigates the influence of waves with varying wave heights, water depths, and wave periods.The structural characteristics, including the perforation rate, perforation number, and perforation type, are considered to meet the construction requirements of SBWs.The main conclusions are as follows: 1.The model developed in this study is well suited for investigating the dynamic response of the seabed.The wave model accurately simulates wave generation, propagation, and reflection processes.Additionally, the seabed model effectively captures liquefaction in the seabed foundation.2. Wave characteristics significantly influence the dynamic response of the seabed.Pore pressure and liquefaction show a positive correlation with wave height and wave period, while exhibiting a negative correlation with water depth.3. The perforation rate of the SBW has a minor effect on pore pressure.Increasing the perforation rate from 5% to 10% leads to a 32% decrease in the average liquefaction On the other hand, front and rear perforations rely on one-sided drainage, affecting flow uniformity and water pressure gradients.This influences flow field stability and exacerbates local seabed responses (as shown in Figure 16).The seaside liquefaction for the front perforation increases by 22.9% and the leeside liquefaction for the rear perforation increases by 45.9% compared to the basic perforation.The increased local flow velocity leads to increased disturbances on the seabed.Thus, it is important to account for these factors when designing structures to withstand such impacts and disturbances, ensuring the stability and safety of SBWs.The liquefaction characteristic values of different breakwater structures are detailed in Table 4.

Conclusions
This study numerically investigates the effects of marine environments and breakwater structures on seabed dynamic responses.Specifically, it investigates the influence of waves with varying wave heights, water depths, and wave periods.The structural characteristics, including the perforation rate, perforation number, and perforation type, are considered to meet the construction requirements of SBWs.The main conclusions are as follows: 1.
The model developed in this study is well suited for investigating the dynamic response of the seabed.The wave model accurately simulates wave generation, propagation, and reflection processes.Additionally, the seabed model effectively captures liquefaction in the seabed foundation.2.
Wave characteristics significantly influence the dynamic response of the seabed.Pore pressure and liquefaction show a positive correlation with wave height and wave period, while exhibiting a negative correlation with water depth.

3.
The perforation rate of the SBW has a minor effect on pore pressure.Increasing the perforation rate from 5% to 10% leads to a 32% decrease in the average liquefaction depth.The increasing perforation number slightly enhances pore pressure and deepens liquefaction due to complex wave reflection and transmission.Among the three different perforation types, basic perforation exerts the minimum seabed pressure.
Front perforation increases liquefaction by 22.9% in the seaside, and rear perforation increases liquefaction by 45.9% in the leeside.4.
In the design of SBWs, it is crucial to consider both wave dissipation and the stability of seabed liquefaction comprehensively.Measures such as reducing the permeability of the seabed can be implementing to enhance the stability of the seabed soil.

Figure 1 .
Figure 1.Sketch of the validation simulation setup for wave interaction.

Figure 1 .
Figure 1.Sketch of the validation simulation setup for wave interaction.

JFigure 2 .
Figure 2. Comparison of the pore pressure between experimental and numerical simulations.

Figure 2 .
Figure 2. Comparison of the pore pressure between experimental and numerical simulations.

Figure 2 .
Figure 2. Comparison of the pore pressure between experimental and numerical simulations.

Figure 3 .
Figure 3.Comparison of the waveforms of experimental and numerical water surfaces.

Figure 4 .
Figure 4. Comparison of the horizontal and vertical velocities of experimental and numerical results.

Figure 3 .
Figure 3.Comparison of the waveforms of experimental and numerical water surfaces.

Figure 2 .
Figure 2. Comparison of the pore pressure between experimental and numerical simulations.

Figure 3 .
Figure 3.Comparison of the waveforms of experimental and numerical water surfaces.

Figure 4 .
Figure 4. Comparison of the horizontal and vertical velocities of experimental and numerical results.

Figure 4 .
Figure 4. Comparison of the horizontal and vertical velocities of experimental and numerical results.

Figure 5 .
Figure 5. Sketch of the simulation setup for wave interaction with a perforated SBW (not to scale).

Figure 5 .
Figure 5. Sketch of the simulation setup for wave interaction with a perforated SBW (not to scale).

Figure 5 .
Figure 5. Sketch of the simulation setup for wave interaction with a perforated SBW (not to scale).

Figure 7 .
Figure 7. Distribution of the seabed response for different wave heights.(The gray half circle represents the location of the semi-circular breakwater).4.1.2.Water Depth Figure 8 illustrates the dynamic response of the seabed to different wave depths.With increasing water depth, the energy transmitted to the seabed surface indirectly decreases, leading to a reduced dynamic response of the seabed.The impact of water depth on liquefaction is weak, with only a slight decrease, of 0.27 cm observed in the liquefaction depth in the seaside.No significant changes are observed in the liquefaction depth in the leeside or in the overall average liquefaction depth.

Figure 7 .
Figure 7. Distribution of the seabed response for different wave heights.(The gray half circle represents the location of the semi-circular breakwater).

Figure 8 Figure 8 .
Figure 8 illustrates the dynamic response of the seabed to different wave depths.With increasing water depth, the energy transmitted to the seabed surface indirectly decreases, leading to a reduced dynamic response of the seabed.The impact of water depth on liquefaction is weak, with only a slight decrease, of 0.27 cm observed in the liquefaction depth in the seaside.No significant changes are observed in the liquefaction depth in the leeside or in the overall average liquefaction depth.J. Mar.Sci.Eng.2024, 12, x FOR PEER REVIEW 10 of 16

Figure 8 .
Figure 8. Distribution of the seabed response for different water depths.(The gray half circle represents the location of the semi-circular breakwater).

Figure 9 .
Figure 9. Distribution of the seabed response for different wave periods.(The gray half circle represents the location of the semi-circular breakwater).

Figure 9 .
Figure 9. Distribution of the seabed response for different wave periods.(The gray half circle represents the location of the semi-circular breakwater).

Figure 10 .
Figure 10.Velocity contours around the SBW with different phases.(Arrows indicate the flow direction; The gray half circle represents the location of the semi-circular breakwater).

Figure 10 .
Figure 10.Velocity contours around the SBW with different phases.(Arrows indicate the flow direction; The gray half circle represents the location of the semi-circular breakwater).

Figure 11 .
Figure 11.Velocity contours around the SBW at different perforation rates.(Arrows indicate the flow direction; The gray half circle represents the location of the semi-circular breakwater).

Figure 11 .
Figure 11.Velocity contours around the SBW at different perforation rates.(Arrows indicate the flow direction; The gray half circle represents the location of the semi-circular breakwater).

Figure 12 .
Figure 12.Distribution of the seabed response at different perforation rates.(The gray half circle represents the location of the semi-circular breakwater).

Figure 12 .
Figure 12.Distribution of the seabed response at different perforation rates.(The gray half circle represents the location of the semi-circular breakwater).

Figure 13 .
Figure 13.Velocity contours around the SBW at different perforation numbers.(Arrows indicate the flow direction; The gray half circle represents the location of the semi-circular breakwater).

Figure 13 .
Figure 13.Velocity contours around the SBW at different perforation numbers.(Arrows indicate the flow direction; The gray half circle represents the location of the semi-circular breakwater).

Figure 14 .
Figure 14.Distribution of the seabed response at different perforation numbers.(The gray half circle represents the location of the semi-circular breakwater).

Figure 14 .
Figure 14.Distribution of the seabed response at different perforation numbers.(The gray half circle represents the location of the semi-circular breakwater).

Figure 15 .
Figure 15.Velocity contours around the SBW with different perforation types.(Arrows indicate the flow direction; The gray half circle represents the location of the semi-circular breakwater).

Figure 15 .Figure 16 .
Figure 15.Velocity contours around the SBW with different perforation types.(Arrows indicate the flow direction; The gray half circle represents the location of the semi-circular breakwater).J. Mar.Sci.Eng.2024, 12, x FOR PEER REVIEW 14 of 16

Figure 16 .
Figure 16.Distribution of the seabed response with different perforation types.(The gray half circle represents the location of the semi-circular breakwater).

Table 1 .
Index properties of wave and breakwater structures.

Table 2 .
Index properties of wave and breakwater structures.

Table 3 .
Liquefaction depths in different marine environments.

Table 4 .
Liquefaction depths of different breakwater structures.

Table 4 .
Liquefaction depths of different breakwater structures.