Numerical Analysis of Dynamic Response and Liquefaction Phenomena in Sandy Seabed Foundation around a Semi-Circular Breakwater under Wave Loading

: Understanding the stability of the seabed foundation holds paramount significance in guaranteeing the safety and structural soundness of the breakwater alongside additional offshore structures. This study aimed to investigate the stability of a sandy seabed foundation around a semi-circular breakwater under wave loading in nearshore areas. A coupled numerical model of waves, a semi-circular breakwater, and the seabed was developed based on the OpenFOAM platform. The VARANS equations were used to govern the wave behavior. Meanwhile, the Biot’s partially dynamic model was employed to numerically simulate the seabed response considering both consolidation under self-weight and dynamic response under wave loading. The effects of various wave parameters, seabed properties, and the radius of the structure on the dynamic response of the seabed and the depth of liquefaction were investigated. The numerical results indicate that an increase in wave height, period, and permeability coefficient intensifies the dynamic response of the seabed soil. Furthermore, an increase in water depth weakened the soil’s dynamic response. There was a negative correlation between the radius of the semi-circular breakwater and the dynamic response. The influence of Poisson’s ratio on the dynamic response of the seabed was relatively small. Furthermore, a stronger dynamic pore pressure response was observed at the connection between the semi-circular breakwater and the rubble foundation.


Introduction
With the continuous development of the maritime industry, new requirements are constantly being proposed for the design and construction of breakwaters.The development of new types of breakwaters has become a significant topic in coastal engineering, among which, semi-circular breakwaters have emerged as an innovative coastal protective structure [1].A semi-circular breakwater consists of a reinforced concrete arc-shaped structure and a bottom plate, situated on a rubble foundation.The decision to incorporate wave-dissipating openings, referred to wave-dissipating holes, on the side facing the waves can be made based on specific needs to decrease wave reflection and wave forces.Moreover, openings can be created on both the wave-facing and leeward sides of the breakwater to facilitate the exchange of seawater.
Throughout the service life of the breakwater, it encounters intricate marine environments.Marine structures are subject to not only substantial wind and wave loading, but they also often rest upon recently deposited, loosely compacted soil layers with a low bearing capacity.Under periodic loading, these soil layers are prone to instability, resulting in structural damage and a reduced service life.The previous literature has documented occurrences of breakwater damage due to soil liquefaction [2][3][4].Therefore, gaining a more comprehensive comprehension of the seabed dynamic response and accurate prediction of potential liquefaction zones near the breakwater have evolved into pivotal factors in the design and maintenance of marine structures, especially under cyclic wave loading and in the presence of new types of breakwaters.
In the initial stages of researching such topics, analytical solutions are often employed [5][6][7].However, these analytical solutions assume that breakwaters exist as impermeable straight lines, neglecting the influence of their geometry and characteristics on wave propagation, as well as the effect of the breakwater's self-weight on the initial stress distribution of the seabed foundation.Compared to theoretical analytical methods, physical model tests can effectively address these issues.Mizutani et al. [8] conducted model tests with submerged breakwaters to investigate the interaction between regular waves, breakwaters, and seabed soil.Building on their research, Mostafa et al. [9] investigated the dynamic response of the seabed under composite breakwaters, presenting a BEM-FEM (Boundary Element Method-Finite Element Method) model that considered the influence of caisson deformation.Further research has been conducted through physical model tests by Oumeraci [10] and Sumer et al. [11].Additionally, due to advancements in computer technology, numerical simulations have been widely employed.Mase et al. [12] developed a numerical model based on the Biot consolidation equation to investigate the distribution of effective stress and pore pressure near composite breakwaters and rubble breakwaters under standing wave conditions.Jeng et al. [13] developed a 2D finite element model based on previous research to investigate the pore water pressure near composite breakwaters, with the maximum pore water pressure occurring at the interface between the caisson and the rubble foundation.Elsafti et al. [14] conducted an investigation into the self-rotational displacement of breakwaters and the dynamic response of surrounding seabed soil under wave action.
Numerous scholars worldwide have conducted comprehensive numerical simulation studies on the interactions between waves, breakwaters, and seabed soil.Zhang Jisheng et al. obtained an improved two-dimensional (2D) numerical model by improving the ultrastatic pore pressure accumulation Sassa model [15].Zheng Dongsheng et al. developed an integrated model (PORO-WSSI II) for wave-seabed-structure interactions (WSSI) and found that under the action of wave loads, there is a dense fluid exchange between the water and the seabed, and there is a strong seepage force in the seabed.An excessive upward seepage force leads to seabed liquefaction [16].Cui Lin et al. established an integrated 3D model of wave-seabed-structure interactions, and studied the liquefaction and settlement of breakwater heads at estuaries [17].
Many studies on the interactions between waves, breakwaters, and the seabed mainly focus on vertical breakwaters [18][19][20].As coastal engineering continues to progress, there remains a dearth of research on emerging breakwater designs.In this study, we developed a numerical model for assessing the dynamic response of the seabed using OpenFOAM.By coupling existing fluid models, a comprehensive numerical framework for studying the interaction among waves, a semi-circular breakwater, and the seabed was established.The numerical model was validated, indicating its accuracy and applicable to the wave-breakwater-seabed interaction problem.The shape and size of the breakwater in the numerical model were designed with reference to the semi-circular breakwater in the South Channel Phase I Project of the Yangtze River Estuary Navigation Improvement Project.The effects of different wave heights, periods, water depths, permeability coefficients, Poisson's ratios, and breakwater radius on the pore water pressures, effective stresses, and liquefaction depths in the seabed soils near semi-circular breakwaters were investigated in order to investigate the dynamic response of the seabed foundations near semi-circular breakwaters and the liquefaction of the seabed foundation near them under the action of wave loads.

Numerical Model 2.1. Flow Model
The flow model is based on the olaFlow solver within OpenFOAM, used to simulate the generation and propagation of waves, as well as their interaction with breakwater structures.The fluid in the model is assumed to be incompressible and continuous.Fluid movement and seepage in porous media are governed by the VARANS equation, which follows the laws of conservation of mass and conservation of momentum. where The coefficients A and B were proposed by Engelund [21] and the specific formulas are shown in Equations ( 3) and (4).
where ρ is density, and u i is velocity vector.p * , g j , ϕ, and D 50 are the pseudo-dynamic pressure, gravitational acceleration, porosity of the material, and average equivalent particle size of the material, respectively.The symbol X represents the position vector.Moreover, µ eff is the effective dynamic viscosity, which considers both the molecular motion viscosity and the turbulent effect (µ eff = µ + ρυ t ).Here, the value of υ t is the turbulent motion viscosity, determined by the selected turbulence model; the k-ε turbulence model was used in this study.
The VOF method was employed within the flow model to accurately capture the free surface.Each grid cell in the calculation domain consists of water and air, and α is used to represent the water volume per unit volume in each grid.A value of α = 1 indicates a cell entirely occupied by water, while α = 0 signifies a cell entirely occupied by air.When α assumes a value between 0 and 1, it implies a cell containing a mixture of both water and air.Within the monomer, alterations in fluid properties are comprehensively described through the indicator function α, which is utilized as follows: where Φ water is the nature of water and Φ air is the nature of air.
In numerical simulation and physical flume tank tests, due to limitations in calculation conditions and test site size, waves are affected and reflected as they propagate to the boundary, resulting in the generation of reflected waves [22,23].The wave model for this numerical simulation was solved using the olaFlow solver, and the active wave absorption method was used at the outlet boundary of the wave model to eliminate wave reflection.In addition, taking 3-5 times the wavelength as the length of the flow field, but only taking the result of the middle stable region, in this region, the wavefront diagram is basically stable and the influence of reflected waves can be ignored.

Seabed Model
Assuming that both pore water and soil skeleton in porous media are compressible, and considering the effect of displacement acceleration of soil particles, the seabed soil is assumed to be isotropic.Additionally, and also considering that fluid flow in porous media follows Darcy's law of permeability, a porous elastic seabed model was established, incorporating the soil inertia term based on the u-p formulation of the Biot equation [24].The governing equation for the dynamic response of the seabed soil is presented below: where x, y, and z are the axes of the Cartesian coordinate system.t, k s , and p e are time, the permeability coefficient of the seabed soil, and the oscillating pore water pressure within the seabed soil caused by waves, respectively.γ w denotes the heaviness of water, n s is the porosity of the soil, and ρ is the average density of the porous seabed, which is calculated using the following equation: where ρ s is the seabed soil density, ρ f represents the pore water density, ε s denotes the soil volumetric strain, and β is the pore water compressibility, where ε s and β can, respectively be expressed as follows: where u s , v s , and w s are the displacements of the soil particles in the x, y, and z directions, respectively.
where K ′ is the apparent bulk modulus of water, S r is the saturation of the soil, P w0 is the absolute pressure of water, and K w is the true bulk modulus of water, taken as 1.95 × 10 9 N/m 2 .The relationship between the pore water pressure and the effective stress is: where σ ′ sx , σ ′ sy , and σ ′ sz are the normal stress of the soil in the x, y, and z directions, and τ sxy , τ syz , and τ szx are the shear stress component of the soil in the x, y and z directions, respectively.Normal stress is the stress perpendicular to the same cross-section.The normal stress here is a positive compression.
Assuming that the seabed soil is a homogeneous elastic material, the stress-strain relationship satisfies the generalized Hooke's law according to the theory of linear elasticity: where G is the shear modulus of the soil, and the relationship between G and elastic modulus E is: where µ is the Poisson's ratio.Bringing Equations ( 12)-( 17) into ( 9)- (11), the relationship between displacement and pore pressure in the seabed soil is: where is the Laplacian Operator.

Model Coupling and Boundary Conditions
The fluid model and the seabed model are coupled together using a one-way coupling algorithm, which relies on the continuity of pressure at the common interface between the two models.The fluid model was employed to simulate the generation, propagation, and interaction of waves with the structure.By solving the VARANS equation and the VOF equation, the fluid model provides information such as the velocity field, free surface elevation, and dynamic water pressure at the fluid-structure interface.On the other hand, the seabed model was used to capture the consolidation process under the self-weight of the structure and the dynamic response under wave loading.The dynamic wave pressure at the bottom of the fluid model is applied to the seabed surface.By iteratively solving the governing equations of the seabed model, the dynamic response of the seabed soil is computed, and the corresponding results are obtained [25].To illustrate a more intuitive depiction of the one-way coupling between the fluid model and the seabed model, the computational process of the wave-structure-seabed system at every time step is presented in Figure 1.
seabed soil is computed, and the corresponding results are obtained [25].To illustrate a more intuitive depiction of the one-way coupling between the fluid model and the seabed model, the computational process of the wave-structure-seabed system at every time step is presented in Figure 1.Boundary conditions play a significant role in determining the behavior and performance of structures and soils in marine environments [26].Boundary conditions clearly influence variations in pressures, forces, displacements, and other variables related to structures and soils [27].In general, there are various methods for defining border conditions in order to extract meaningful results based on them [28,29].Based on the specific boundary conditions, the governing equations of the seabed model are solved to obtain the dynamic response of the surrounding seabed soil under wave loading for the semi- Boundary conditions play a significant role in determining the behavior and performance of structures and soils in marine environments [26].Boundary conditions clearly influence variations in pressures, forces, displacements, and other variables related to structures and soils [27].In general, there are various methods for defining border conditions in order to extract meaningful results based on them [28,29].Based on the specific boundary conditions, the governing equations of the seabed model are solved to obtain the dynamic response of the surrounding seabed soil under wave loading for the semi-circular breakwater.The boundary conditions for the seabed computational domain are set as follows in Figure 2: (1) At the sea bed surface (i.e., at z = h), it is commonly assumed that there are no vertical effective normal stresses and shear stresses.The pore water pressure at this location is equal to the dynamic pressure applied on the seabed surface due to wave loading: (2) At the lower boundary of the seabed (i.e., at z = 0), it is generally assumed to be an impermeable rigid bedrock, i.e., no displacement in the horizontal and vertical directions and no vertical seepage:

Mesh Convergence Verification
To ensure the accuracy and efficiency of the model calculations, a convergence test was performed on the grid resolution.Figures 3 and 4 show the temporal variations in the wave surface and the pore water pressure under different grid resolutions.The specific parameters used in the wave model and seabed model are as follows: wave height H = 0.03 m, wave period T = 1.4 s, water depth d = 0.21 m, seabed thickness h = 0.2 m, saturation degree Sr = 0.99, shear modulus G = 1 × 10 8 Pa, Poisson's ratio = 0.33, porosity = 0.3, and permeability coefficient = 2.2 × 10 −3 m/s.The grid size  x = L/150,  z = H/20 was se- lected to balance time and cost constraints while ensuring the accuracy of the calculations.In the calculation process, the grid size required for wave model convergence is significantly smaller than that for the seabed model.However, considering the presence of the breakwater and aiming to enhance calculation precision, this study used grid sizes that are identical in both the wave model and the seabed model.Based on the results of the time step test, the time step in the wave model was set to 0.001 s, while in the seabed model, the time step was set to 0.01 s.

Mesh Convergence Verification
To ensure the accuracy and efficiency of the model calculations, a convergence test was performed on the grid resolution.Figures 3 and 4 show the temporal variations in the wave surface and the pore water pressure under different grid resolutions.The specific parameters used in the wave model and seabed model are as follows: wave height H = 0.03 m, wave period T = 1.4 s, water depth d = 0.21 m, seabed thickness h = 0.2 m, saturation degree Sr = 0.99, shear modulus G = 1 × 10 8 Pa, Poisson's ratio = 0.33, porosity = 0.3, and permeability coefficient = 2.2 × 10 −3 m/s.The grid size ∆x = L/150, ∆z = H/20 was selected to balance time and cost constraints while ensuring the accuracy of the calculations.In the calculation process, the grid size required for wave model convergence is significantly smaller than that for the seabed model.However, considering the presence of the breakwater and aiming to enhance calculation precision, this study used grid sizes that are identical in both the wave model and the seabed model.Based on the results of the time step test, the time step in the wave model was set to 0.001 s, while in the seabed model, the time step was set to 0.01 s.
In the calculation process, the grid size required for wave model convergence is significantly smaller than that for the seabed model.However, considering the presence of the breakwater and aiming to enhance calculation precision, this study used grid sizes that are identical in both the wave model and the seabed model.Based on the results of the time step test, the time step in the wave model was set to 0.001 s, while in the seabed model, the time step was set to 0.01 s.

Model Validation with Wave Flume Test
Mizutani et al. [8] conducted a series of laboratory experiments to study the interactions between waves, the seabed, and the breakwater.The experimental setup is depicted in Figure 5: the left side represents the wave generation area, and the right side represents the wave dissipation area.A submerged breakwater was placed on the seabed, and sensors were installed at locations A, B, C, and D to record the pore pressure.Wave gauges were placed at locations a, b, c, and d to measure the wave height.Table 1 presents the values of the wave, seabed soil, and breakwater parameters used in Mizutani's laboratory experiments.

Model Validation with Wave Flume Test
Mizutani et al. [8] conducted a series of laboratory experiments to study the interactions between waves, the seabed, and the breakwater.The experimental setup is depicted in Figure 5: the left side represents the wave generation area, and the right side represents the wave dissipation area.A submerged breakwater was placed on the seabed, and sensors were installed at locations A, B, C, and D to record the pore pressure.Wave gauges were placed at locations a, b, c, and d to measure the wave height.Table 1 presents the values of the wave, seabed soil, and breakwater parameters used in Mizutani's laboratory experiments.Figures 6 and 7 illustrate the comparison between the numerical simulation results in this study and the experimental results of Mizutani et al. regarding wave profiles and pore water pressure at various seabed points under wave action.In Figure 6, it is evident that the numerical simulation results closely match the experimental results at points a and b in front of the breakwater.Meanwhile, at points c and d located behind the breakwater, the wave experienced attenuation processes.These characteristics are clearly reflected in the experimental results of Mizutani et al., and the numerical simulation results in this study also captured these features.While there exists some numerical discrepancy when compared to the experimental results, the overall trend of the wave profile remains in substantial agreement, indicating the reliability of the fluid model used in this study.Figure 7 shows that the simulated curves of pore water pressure in the sandy seabed under wave loading at points A, B, C, and D correspond closely to the monitoring results of Mizutani et al.'s experiments.This demonstrates that the seabed model employed in this study is suitable for simulating the interaction between waves, the seabed, and breakwaters.

Parameter Setting
The shape and size of the semi-circular breakwater were designed with consideration of the semi-circular breakwater size utilized in the South Channel Phase I Project of the Yangtze River Estuary Navigation Improvement Project.Previous studies have indicated that the presence of openings in the breakwater does not have a significant influence on wave attenuation under high water levels.Given that this study concentrated on submerged breakwaters with high water levels, the design of the semi-circular components in the actual engineering can be simplified.The frontal view of the semi-circular breakwater model is shown in Figure 8.The breakwater consisted of an upper semi-circular component, a lower plate section, and a layer of rubble.The width of the rubble was is 44.8 m, the outer diameter of the semi-circular component was 5 m, the inner diameter was 4.5 m, and the length of the lower straight section was 1.5 m.The coordinate origin was located at the seabed's center, with wave propagation occurring along the positive x-axis.The numerical wave flume's length (Ls) was set at four times the wavelength.

Parameter Setting
The shape and size of the semi-circular breakwater were designed with consideration of the semi-circular breakwater size utilized in the South Channel Phase I Project of the Yangtze River Estuary Navigation Improvement Project.Previous studies have indicated that the presence of openings in the breakwater does not have a significant influence on wave attenuation under high water levels.Given that this study concentrated on submerged breakwaters with high water levels, the design of the semi-circular components in the actual engineering can be simplified.The frontal view of the semi-circular breakwater model is shown in Figure 8.The breakwater consisted of an upper semi-circular component, a lower plate section, and a layer of rubble.The width of the rubble was is 44.8 m, the outer diameter of the semi-circular component was 5 m, the inner diameter was 4.5 m, and the length of the lower straight section was 1.5 m.The coordinate origin was located at the seabed's center, with wave propagation occurring along the positive x-axis.The numerical wave flume's length (Ls) was set at four times the wavelength.According to the research conducted by Jeng et al. [13], it has been shown that relatively high pore water pressures exist at the toe of the rubble foundation, at the interface between the breakwater and the rubble foundation, and below the center of the caisson.Therefore, this study concentrated on investigating the dynamic response of the seabed at points A, C, and D. The parameters for the waves, breakwater, and seabed are listed in Table 2.According to the research conducted by Jeng et al. [13], it has been shown that relatively high pore water pressures exist at the toe of the rubble foundation, at the interface between the breakwater and the rubble foundation, and below the center of the caisson.Therefore, this study concentrated on investigating the dynamic response of the seabed at points A, C, and D. The parameters for the waves, breakwater, and seabed are listed in Table 2.

Effect of Wave Characteristics
In this section, we investigate the influence of wave characteristics on the dynamic response of the seabed.The wavelength is directly affected by the wave period and water depth, and changes in these three variables directly affect the wave loads acting on the seabed surface, consequently influencing the stability of the seabed.
Figures 9 and 10 show the distribution curves of the maximum pore water pressure and maximum effective stress concerning the seabed depth at points A, C, and D for different wave heights (i.e., 0.5 m, 1.25 m, 2 m).In Figure 9, it is evident that the maximum pore water pressure at the same seabed depth significantly increased with increasing wave height.This occurs because larger wave heights result in more substantial wave energy impacting the seabed surface, intensifying the dynamic response of the seabed around the semi-circular breakwater and leading to a significant increase in pore water pressure.This effect was more particularly prominent in shallow seabed areas, specifically within the range of 0.8 < z/h < 1.In regions where z/h < 0.8, the influence was less significant, which is also the reason why liquefaction under wave loading tends to occur in shallow seabed regions.As the seabed depth increases, wave energy propagates and dissipates within the seabed, leading to a gradual decrease in pore water pressure.The maximum pore water pressure decreased most rapidly within a depth range of approximately 4 m below the seabed surface, followed by a slower decrease.Furthermore, among points A, C, and D, it became evident that the maximum pore water pressure was situated at point C, aligning with the interface between the semi-circular breakwater and the rubble foundation.This consistency corresponds to the conclusions drawn by Jeng et al. [13].This phenomenon might be attributed to the modified drainage conditions within the foundation soil, induced by the presence of the semi-circular breakwater, along with the heightened wave reflection and superposition effects near this vicinity.Consequently, this region experienced elevated wave heights and stronger dynamic pore pressure response in the seabed.The wave periods considered in this study were 6 s, 8 s, and 10 s.Figures 11 and 12 illustrate the distribution of the maximum pore water pressure and maximum vertical effective stress at different points around the semi-circular breakwater for different wave periods.From the figures, it is evident that the maximum pore water pressure and maximum vertical effective stress in the seabed at the same depth increased with longer wave periods.This occurs because the wave period is closely related to the wave length, and as the wave period increases, so does the wave length.A larger wave length contains more  Figure 10 reveals that the influence of wave height on the maximum vertical effective stress in the soil around the semi-circular breakwater resembles that of the pore water pressure.The maximum vertical effective stress at the same seabed depth increased with increasing wave height.Among points A, C, and D, point C exhibited the greatest variation in maximum vertical effective stress with wave height, followed by point D, while point A showed the smallest variation.This is due to point A experiencing weaker wave reflection and superposition effects, followed by point D, with point C experiencing the strongest effects.Similar to the maximum pore water pressure, the sudden change in the maximum vertical effective stress occurred at approximately z/h = 0.8.Within the range of 0.8 h < z < h, the maximum vertical effective stress was relatively small.According to the effective stress principle, this range is more susceptible to experience seabed liquefaction according to the criteria for transient liquefaction assessment.
The wave periods considered in this study were 6 s, 8 s, and 10 s.Figures 11 and 12 illustrate the distribution of the maximum pore water pressure and maximum vertical effective stress at different points around the semi-circular breakwater for different wave periods.From the figures, it is evident that the maximum pore water pressure and maximum vertical effective stress in the seabed at the same depth increased with longer wave periods.This occurs because the wave period is closely related to the wave length, and as the wave period increases, so does the wave length.A larger wave length contains more energy, leading to an amplified energy input on the seabed surface and intensifying the dynamic response of the seabed.The water depths considered were 5, 7.5, and 10 m.Figures 13 and 14 show the maximum values of the pore water pressure and vertical effective stress at each measurement point at different water depths around the semi-circular breakwater respectively.Upon analyzing Figures 13 and 14, it can be seen that the maximum pore water pressure and the maximum vertical effective stress decreased as water depth increased while maintaining the same seabed depth.This decline can be attributed to the reduction in wave pressure exerted upon the seabed surface due to heightened water depth, subsequently leading to the diminishment of both pore water pressure and vertical effective stress within the seabed soil.Furthermore, the analysis revealed that in cases of a shallow water depth, the maximum vertical effective stress of the seabed soil underwent more noticeable and rapid changes along the depth distribution of the seabed.When the water depth reached 10 m, the pore water pressure within the range of 0 < z < 0.8 h remained almost constant.Figures 11 and 12 show that the alterations in the wave period impact both the maximum vertical effective stress and maximum pore water pressure at the same seabed depth, with variations occurring around z = 0.8 h.Within the interval of 0 < z < 0.8 h, both the maximum pore water pressure and maximum vertical effective stress exhibited a gradual increase.In the range of 0.8 h < z < h, the rate of increase became more rapid.Notably, the maximum pore water pressure tended to increase significantly in this region.Similarly, the maximum vertical effective stress also showed an increasing trend within the range of 0 < z < 0.8 h, reaching its peak value at a depth of approximately 0.8 h.In the range of 0.8 h < z < h, the maximum vertical effective stress decreased as it approached the seabed surface.This indicates that the seabed is more prone to experience transient liquefaction in this region.
The variation in dynamic response due to changes in wave period follows a similar pattern to that observed for wave height.Point C showed the most significant variation, followed by point D, while point A showed the weakest dynamic response among the three points (A, C, and D).
The water depths considered were 5, 7.5, and 10 m.Figures 13 and 14 show the maximum values of the pore water pressure and vertical effective stress at each measurement point at different water depths around the semi-circular breakwater respectively.Upon analyzing Figures 13 and 14, it can be seen that the maximum pore water pressure and the maximum vertical effective stress decreased as water depth increased while maintaining the same seabed depth.This decline can be attributed to the reduction in wave pressure exerted upon the seabed surface due to heightened water depth, subsequently leading to the diminishment of both pore water pressure and vertical effective stress within the seabed soil.Furthermore, the analysis revealed that in cases of a shallow water depth, the maximum vertical effective stress of the seabed soil underwent more noticeable and rapid changes along the depth distribution of the seabed.When the water depth reached 10 m, the pore water pressure within the range of 0 < z < 0.8 h remained almost constant.

Sea Bed Characteristics
Compared to the effects of wave characteristics, changes in the physical properties of the seabed will directly influence the distribution of the pore water pressure and vertical effective stress in the seabed.This section focuses on the effect of soil permeability and Poisson's ratio variation on the dynamic response of the seabed.
Figures 15 and 16 illustrate the variations in the maximum pore water pressure and maximum vertical effective stress within the seabed at specified points A, C, and D, as a

Sea Bed Characteristics
Compared to the effects of wave characteristics, changes in the physical properties of the seabed will directly influence the distribution of the pore water pressure and vertical effective stress in the seabed.This section focuses on the effect of soil permeability and Poisson's ratio variation on the dynamic response of the seabed.
Figures 15 and 16 illustrate the variations in the maximum pore water pressure and maximum vertical effective stress within the seabed at specified points A, C, and D, as a function of depth, for permeability coefficients of 0.0001 m/s, 0.001 m/s, and 0.01 m/s.As shown in Figures 15 and 16, the maximum pore water pressure within the seabed at the same depth increased as the permeability coefficient increased.As the depth of the seabed increased, the maximum pore water pressure gradually decreased.It tended to stabilize as it approached the bottom of the seabed.Furthermore, as the depth of the seabed increased, the decay rate of the maximum pore water pressure under different permeability coefficients varied at different locations.For a permeability coefficient of k s = 0.0001 m/s, the decay rate diminished around z = 0.9 h.In the case of k s = 0.001 m/s, the attenuation rate dropped at z = 0.8 h.In situations where k s = 0.01 m/s, the attenuation rate did not exhibit a significant inflection point, as the change in the permeability coefficient directly impacted the fluid flow within the seabed pores.The increase in permeability coefficient led to an increase in the fluid velocity in the soil pores.Among the three points A, C, and D, the dynamic response at point C was the strongest.As shown in Figure 16, it is evident that at the same seabed depth, the maximum vertical effective stress decreased with the increase in the permeability coefficient.As the seabed's depth increased, the effective stress experienced a gradual increase, and its peak appeared near the seabed's surface.The smaller the permeability coefficient, the closer its peak was to the surface of the seabed.Across the various measurement points near the semi-circular breakwater, the vertical effective stress reached its maximum at the interface between the breakwater and the riprap foundation.that at the same seabed depth, the maximum vertical effective stress decreased with the increase in the permeability coefficient.As the seabed's depth increased, the effective stress experienced a gradual increase, and its peak appeared near the seabed's surface.
The smaller the permeability coefficient, the closer its peak was to the surface of the seabed.Across the various measurement points near the semi-circular breakwater, the vertical effective stress reached its maximum at the interface between the breakwater and the riprap foundation.Figures 17 and 18 show the variations in the maximum pore water pressure and maximum vertical effective stress with depth for different Poisson's ratios (i.e., 0.2, 0.3, 0.4) in the seabed.As depicted in Figure 17, it is evident that the maximum pore water pressure in the seabed soil at the same depth can be divided into two parts as the Poisson's ratio changes: in the range of 0 < z < 0.8 h, the maximum pore water pressure decreased with an increase in the Poisson's ratio; in the range of 0.8 h < z < h, the variation in the Poisson's ratio minimally impacted the pore water pressure, and the maximum values remained almost the same for all Poisson's ratios at the same depth.Analyzing Figure 18, it becomes apparent that the maximum vertical effective stress in the seabed soil decreased with an increase in the Poisson's ratio at the same depth.It initially increased and subsequently decreased as the seabed depth increased, peaking at z = 0.8 h.In comparison to the influences of wave characteristics and soil permeability coefficients on the dynamic response of the seabed, the effect of the Poisson's ratio on the seabed's dynamic response was relatively small.Figures 17 and 18 show the variations in the maximum pore water pressure and maximum vertical effective stress with depth for different Poisson's ratios (i.e., 0.2, 0.3, 0.4) in the seabed.As depicted in Figure 17, it is evident that the maximum pore water pressure in the seabed soil at the same depth can be divided into two parts as the Poisson's ratio changes: in the range of 0 < z < 0.8 h, the maximum pore water pressure decreased with an increase in the Poisson's ratio; in the range of 0.8 h < z < h, the variation in the Poisson's ratio minimally impacted the pore water pressure, and the maximum values remained almost the same for all Poisson's ratios at the same depth.Analyzing Figure 18, it becomes apparent that the maximum vertical effective stress in the seabed soil decreased with an increase in the Poisson's ratio at the same depth.It initially increased and subsequently decreased as the seabed depth increased, peaking at z = 0.8 h.In comparison to the influences of wave characteristics and soil permeability coefficients on the dynamic response of the seabed, the effect of the Poisson's ratio on the seabed's dynamic response was relatively small.

Geometrical Characteristics of the Semi-Circular Breakwater
The dynamic response of the seabed results from the coupled interaction between the waves, seabed, and structures.Structural dimensional changes can significantly impact the oscillatory response of the surrounding seabed soil [30].Serving as a novel protective structure, the semi-circular breakwater boasts a distinctive design with an arched shape, and variations in its radius directly influence the reflection and diffraction of waves,

Geometrical Characteristics of the Semi-Circular Breakwater
The dynamic response of the seabed results from the coupled interaction between the waves, seabed, and structures.Structural dimensional changes can significantly impact the oscillatory response of the surrounding seabed soil [30].Serving as a novel protective structure, the semi-circular breakwater boasts a distinctive design with an arched shape, and variations in its radius directly influence the reflection and diffraction of waves, thereby affecting the oscillatory response of the seabed.Considering the practical applications, breakwater dimension design is also a crucial consideration.This section focuses on studying the influence of semi-circular breakwater radius variations on the dynamic response of the seabed.
In investigating the influence of the radius of the arched structure, a water depth of 10 m was employed, along with specific radius values of 5 m, 7 m, and 10 m.Additional parameter values are cross-referenced from Table 2. Figures 19 and 20 illustrate the distribution of the maximum pore water pressure and maximum effective stress concerning the seabed depth at different locations around the semi-circular breakwater for different radii.Figure 19 indicates that within the range of 0.8 h < z < h, the maximum pore water pressure at the same location decreased with increasing seabed depth, exhibiting a reduced decay rate at z = 0.8 h.For the range of 0 < z < 0.8 h, the pore water pressure within this interval was not high, and the maximum pore water pressure remained relatively constant.This was due to the attenuation of wave energy within a water depth of 10 m.As it propagates downward from the seabed surface, the fluid flow encounters frictional resistance with the soil particles, leading to further energy attenuation.It was observed that the maximum pore pressure decreased with an increase in the radius, especially for R = 10 m within the top 40% of the seabed depth.This phenomenon arose due to the selected reference points, which are situated below the breakwater and remained relatively unaffected by the wave load from above when the radius increased to 10 m. pore pressure decreased with an increase in the radius, especially for R = 10 m within the top 40% of the seabed depth.This phenomenon arose due to the selected reference points, which are situated below the breakwater and remained relatively unaffected by the wave load from above when the radius increased to 10 m.It is evident from Figure 20 that, at the same seabed depth, the maximum effective stress decreased with an increasing radius.Among the three points A, C, and D, point C exhibited significantly higher maximum effective stress compared to the other two points.The dynamic response of the seabed at point D showed the smallest variation.This can be partly attributed to wave reflection and superposition at point C and also to the larger radius of the breakwater, which provides a stronger hindrance to wave propagation.
The tendency of the maximum effective stress for R = 10 m was different from that

Liquefaction Depth
Liquefaction is the phenomenon in which the soil changes from a solid state to a liquid state when the excess pore pressure is equal to the effective stress of the overlying soil layer.This leads to a loss of soil strength and stiffness.Seabed liquefaction can be divided into two types: instantaneous liquefaction and residual liquefaction.Instantaneous liquefaction is generated by an oscillating pore pressure and exhibits significant periodicity.Residual liquefaction is caused by the accumulation of excess pore pressure after multiple wave cycles.This article mainly studies instantaneous liquefaction.During wave propagation, the oscillating pore pressure varied periodically with the variation in wave peaks and valleys.The variation in pore pressure along the depth of the seabed can be effectively used to capture wave-induced liquefaction phenomena.Near the wave trough, the vertical gradient of oscillating pore pressure generates upward seepage force.When the excess pore pressure exceeds the effective weight of the soil, instantaneous liquefaction may occur.
In the natural environment, the seabed soil undergoes consolidation over time under the influence of gravity, reaching an initial state of equilibrium.However, the construction of a breakwater disrupts this initial equilibrium state and generates disturbances in the surrounding seabed.The breakwater and the adjacent seabed foundation enter a new state of consolidation under the influence of gravity [31].This new equilibrium state is used as the initial state for liquefaction analysis of the seabed in this study.
This section aims to investigate the influence of various parameters (wave characteristics, seabed properties, breakwater characteristics) on the maximum liquefaction depth of the seabed around the semi-circular breakwater.The liquefaction depth is the depth of It is evident from Figure 20 that, at the same seabed depth, the maximum effective stress decreased with an increasing radius.Among the three points A, C, and D, point C exhibited significantly higher maximum effective stress compared to the other two points.The dynamic response of the seabed at point D showed the smallest variation.This can be partly attributed to wave reflection and superposition at point C and also to the larger radius of the breakwater, which provides a stronger hindrance to wave propagation.
The tendency of the maximum effective stress for R = 10 m was different from that for R = 5 m and R = 7 m.The reason is that, when the water depth is 10 m, and the radius of the breakwater is 10 m, a part of the breakwater will be exposed to the water surface during wave propagation.This causes a severe increase in wave reflection on the wave-facing side of the breakwater, an increase in wave height, and a decrease in effective stress.The soil below the seabed is less affected by this, so the changes are not significant.

Liquefaction Depth
Liquefaction is the phenomenon in which the soil changes from a solid state to a liquid state when the excess pore pressure is equal to the effective stress of the overlying soil layer.This leads to a loss of soil strength and stiffness.Seabed liquefaction can be divided into two types: instantaneous liquefaction and residual liquefaction.Instantaneous liquefaction is generated by an oscillating pore pressure and exhibits significant periodicity.Residual liquefaction is caused by the accumulation of excess pore pressure after multiple wave cycles.This article mainly studies instantaneous liquefaction.During wave propagation, the oscillating pore pressure varied periodically with the variation in wave peaks and valleys.The variation in pore pressure along the depth of the seabed can be effectively used to capture wave-induced liquefaction phenomena.Near the wave trough, the vertical gradient of oscillating pore pressure generates upward seepage force.When the excess pore pressure exceeds the effective weight of the soil, instantaneous liquefaction may occur.
In the natural environment, the seabed soil undergoes consolidation over time under the influence of gravity, reaching an initial state of equilibrium.However, the construction of a breakwater disrupts this initial equilibrium state and generates disturbances in the surrounding seabed.The breakwater and the adjacent seabed foundation enter a new state of consolidation under the influence of gravity [31].This new equilibrium state is used as the initial state for liquefaction analysis of the seabed in this study.
This section aims to investigate the influence of various parameters (wave characteristics, seabed properties, breakwater characteristics) on the maximum liquefaction depth of the seabed around the semi-circular breakwater.The liquefaction depth is the depth of liquefaction in the vertical direction (i.e., z-direction) of the seabed.Figure 21 presents the variations in the maximum liquefaction depth of the seabed under different wave characteristics, seabed properties, and semi-circular geometries.The rectangular area in the figure represents the location where the semi-circular breakwater and the stone foundation interact.Given that the breakwaters analyzed in this study are all submerged breakwaters, which are more aesthetically appealing than emerged breakwaters, some waves can pass over the breakwater and reach the area behind it, causing disturbances to the soil.Therefore, liquefaction can also occur behind the breakwater.However, the breakwater reflects the waves on its front side, impeding wave propagation.The dynamic response of the seabed soil behind the breakwater is smaller than that on the front side.

Conclusions
This study systematically analyzed the effects of wave parameters, seabed characteristics, and geometric features of a semicircular breakwater on the dynamic response of the surrounding seabed using the established two-dimensional numerical model of wavebreakwater-seabed system.The following conclusions were drawn: 1. 1The model established in this study can effectively be applied to the study of the The analysis indicated that the maximum liquefaction depth was positively correlated with wave height, wave period, and permeability coefficient, while it was negatively correlated with water depth and Poisson's ratio.Liquefaction almost did not occur when the wave height was 0.5 m.As the wave height increased to 2 m, the maximum liquefaction depth reached 1.0 m.Similarly, as the wave period increased from 6 s to 10 s, the maximum liquefaction depth reached 0.9 m.When the water depth was 0.5 m, the maximum liquefaction depth was 0.8 m.As the water depth increased to 10 m, the maximum liquefaction depth was 0.4 m.This phenomenon occurs because with an increase in water depth, the energy transmitted to the seabed surface from the wave load diminishes, leading to a reduced dynamic response of the soil.The larger the permeability coefficient, the greater the maximum depth of the transient liquefaction.When the permeability coefficients were 0.0001, 0.001, and 0.01 m/s, the corresponding maximum liquefaction depths were 0.4, 0.5, and 0.7 m, respectively.When the Poisson's ratio was 0.2, 0.3, and 0.4, the corresponding maximum liquefaction depths were 0.9, 0.7, and 0.5 m, respectively.Furthermore, as the radius increased to 10 m, the breakwater's ability to obstruct waves strengthened.Only a few waves could pass the breakwater and reach its rear position, resulting in reduced disturbances behind it.Therefore, liquefaction did not occur behind the breakwater.Therefore, when considering enhanced protection in coastal areas, increasing the radius of the breakwater is recommended.

Conclusions
This study systematically analyzed the effects of wave parameters, seabed characteristics, and geometric features of a semicircular breakwater on the dynamic response of the surrounding seabed using the established two-dimensional numerical model of wave-breakwater-seabed system.The following conclusions were drawn: 1.
The model established in this study can effectively be applied to the study of the oscillatory response in the seabed surrounding a semicircular breakwater.The fluid model accurately simulates wave generation, propagation, and the reflection and superposition phenomena in front of the breakwater.The seabed model captures the oscillatory response of the seabed foundation under dynamic wave loading; 2.
The dynamic response of the seabed is influenced by wave characteristics, seabed properties, and the geometric features of the semicircular breakwater.The dynamic response of seabed soil is intensified by increasing wave height, wave period, and permeability coefficient, whereas the dynamic response is weakened by increasing water depth.There is a negative correlation between the radius of the semicircular breakwater and the dynamic response of the seabed, while the influence of the Poisson's ratio is relatively small; 3.
Under wave loading, the maximum pore water pressure in the seabed decreases with increasing seabed depth.The attenuation rate of the maximum pore water pressure shows a sudden change in shallow seabed areas, with the peak value occurring within this region.The maximum effective stress in the seabed shows an increasing trend followed by a decrease as the seabed depth increases, with the peak value occurring in the shallow seabed region; 4.
The presence of the semicircular breakwater causes changes in the distribution of the pore water pressure near its foundation, with the most intense dynamic response observed at the interface between the breakwater and the rubble foundation.

Figure 3 .
Figure 3. Free surface time history curve at a point in the wave model with different grid densities.Figure 3. Free surface time history curve at a point in the wave model with different grid densities.

Figure 3 .
Figure 3. Free surface time history curve at a point in the wave model with different grid densities.Figure 3. Free surface time history curve at a point in the wave model with different grid densities.J. Mar.Sci.Eng.2024, 12, x FOR PEER REVIEW 8 of 29

Figure 4 .
Figure 4. Pore water pressure variation curve with time at a point in the wave model with different grid densities.

Table 1 .
Parameters used in the Mizutani et al. flume test.

Figure 4 .
Figure 4. Pore water pressure variation curve with time at a point in the wave model with different grid densities.

Figures 6
Figures 6 and 7 illustrate the comparison between the numerical simulation results in this study and the experimental results of Mizutani et al. regarding wave profiles and pore water pressure at various seabed points under wave action.In Figure 6, it is evident that the numerical simulation results closely match the experimental results at points a and b in front of the breakwater.Meanwhile, at points c and d located behind the breakwater, the wave experienced attenuation processes.These characteristics are clearly reflected in the experimental results of Mizutani et al., and the numerical simulation results in this study also cap-

Figure 6 .
Figure 6.Comparison of the waveforms between numerical simulation results in this study and the experimental results of Mizutani et al. [8].

Figure 6 .
Figure 6.Comparison of the waveforms between numerical simulation results in this study and the experimental results of Mizutani et al. [8].

Figure 7 .
Figure 7.Comparison of the pore water pressure between numerical simulation results and the experimental results.

Figure 7 .
Figure 7.Comparison of the pore water pressure between numerical simulation results and the experimental results.

Figure 8 .
Figure 8.The sketch of the numerical set-up.

Figure 8 .
Figure 8.The sketch of the numerical set-up.

Figure 9 .Figure 9 .
Figure 9. Variation curve of the maximum pore water pressure with depth at each measurement point around the semi-circular breakwater at different wave heights.

Figure 10 .
Figure 10.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater at different wave heights.

Figure 10 .
Figure 10.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater at different wave heights.

J
. Mar. Sci.Eng.2024, 12, x FOR PEER REVIEW 16 of 29 followed by point D, while point A showed the weakest dynamic response among the three points (A, C, and D).

Figure 11 .Figure 11 .
Figure 11.Variation curve of the maximum pore water pressure with depth at each measurement point around the semi-circular breakwater at different periods.

Figure 12 .
Figure 12.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater at different periods.

Figure 12 .
Figure 12.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater at different periods.

Figure 13 .Figure 13 .Figure 14 .
Figure 13.Variation curve of the maximum pore water pressure with depth at each measurement point around the semi-circular breakwater at different water depths.

Figure 14 .
Figure 14.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater at different water depths.

( a )Figure 15 .Figure 15 .
Figure 15.Variation curve of the maximum pore water pressure with depth at each measurement point around the semi-circular breakwater with different levels of soil permeability.

Figure 16 .
Figure 16.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater with different levels of soil permeability.

Figure 16 .
Figure 16.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater with different levels of soil permeability.

Figure 17 .
Figure 17.Variation curve of the maximum pore water pressure with depth at each measurement point around the semi-circular breakwater with different Poisson's ratios.

Figure 17 .
Figure 17.Variation curve of the maximum pore water pressure with depth at each measurement point around the semi-circular breakwater with different Poisson's ratios.

Figure 18 .
Figure 18.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater with different Poisson's ratios.

Figure 18 .
Figure 18.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater with different Poisson's ratios.

Figure 19 .
Figure 19.Variation curve of the maximum pore water pressure with depth at each measurement point around the semi-circular breakwater with different radii.

Figure 19 .
Figure 19.Variation curve of the maximum pore water pressure with depth at each measurement point around the semi-circular breakwater with different radii.

Figure 20 .
Figure 20.Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater with different radii.

JFigure 21 .
Figure 21.Effect of different parameters on the maximum liquefaction depth of the seabed around a semi-circular breakwater.

Figure 21 .
Figure 21.Effect of different parameters on the maximum liquefaction depth of the seabed around a semi-circular breakwater.

Table 1 .
Parameters used in the Mizutani et al. flume test.

Table 2 .
Parameters used in the numerical studies.
Variation curve of the maximum value of effective stress with depth at each measurement point around the semi-circular breakwater with different radii.