Influence of Wave-Absorbing Chamber Width on the Wave Attenuation Performance of Perforated Caisson Sitting on Rubble-Mound Foundation

A Jarlan-type perforated caisson consisted of a perforated front wall, a solid rear wall, and a wave-absorbing chamber between them. The wave-absorbing chamber was the main feature of the perforated caisson, and its width had a great effect on wave attenuation performance. In this study, a larger range of the wave-absorbing chamber width was observed in model experiments to investigate the effect on wave attenuation performance including the reflection coefficients and the horizontal wave forces of a perforated caisson sitting on a rubble-mound foundation. A resistance-type porosity numerical model based on the volume-averaged Reynolds-averaged Navier–Stokes (VARANS) equations was validated by comparing the present results with those of previously reported and present experiments. The validated numerical model was then used for extended research. It was found that the reflection coefficients, the total horizontal wave force, and its components all tended to oscillate in a decrease ⟶ increase ⟶ decrease manner with increasing the wave-absorbing chamber width. The reflection coefficients and wave forces acting on both sides of the perforated front wall were found to be synchronized regardless of perforation ratio or the rubble-mound foundation height.


Introduction
Perforated caissons were first proposed by Jarlan [1] in the 1960s. Since then, many researchers studied this coastal structure in analytical, physical, and numerical ways [2][3][4][5][6][7][8]. Some other new types of caisson were also proposed and studied in recent years [9,10]. A Jarlan-type perforated caisson consisted of a perforated front wall, a solid rear wall, and a wave-absorbing chamber between them. When the incident wave overlapped with the reflected waves from the perforated front wall and the solid rear wall, the wave energy can be greatly reduced. Compared with the traditional caisson, a perforated caisson reduced reflection coefficients and the wave forces. In addition, a perforated caisson had the advantages of easy construction and a wide range of applications, making it a key structure in coastal and port engineering.
Chen et al. [11], Ma et al. [12], and Jiang et al. [13,14] carried out physical model experiments to study the irregular wave interaction with a perforated caisson breakwater on an impermeable and a permeable seabed. In those experiments, the relative wave-absorbing chamber width was between 0.064 and 0.208. It was pointed out that both the reflection coefficients and the horizontal wave forces acting on a perforated caisson were both correlated quadratically with the wave-absorbing chamber width, but this conclusion had not been confirmed for a perforated caisson sitting on a rubble-mound foundation. Xing [15] studied the irregular wave interaction with a perforated caisson breakwater on a rubble-mound foundation of three different heights, with the wave-absorbing chamber width between 0.064 and 0.208. Within this range, both the reflection coefficients and the horizontal wave forces acting on a perforated caisson were also correlated quadratically with the wave-absorbing chamber width.
So far, it is generally believed that the wave-absorbing chamber width is one of the key factors affecting the wave attenuation performance of a perforated caisson sitting on a rubble-mound foundation. However, when the wave-absorbing chamber width expanded to a larger range, can the changes in reflection coefficients and horizontal wave forces still be described by a quadratic function? Most experimental studies on perforated caissons limit the wave-absorbing width to a small range, so this problem was not solved. Liu et al. [16] attempted to analytically calculate the reflection coefficients of a perforated caisson based on potential flow theory. In their research, the perforated caisson was bedded on an impermeable step, which differs from a rubble-mound foundation used in engineering practice; therefore, the applicability of their analytical solution had great limitations.
Accordingly, in this study, by giving a larger range of the wave-absorbing chamber width, model experiments were carried out to study the effect of the wave-absorbing chamber width on the attenuation performance of a perforated caisson sitting on a rubble-mound foundation. Meanwhile, a resistance-type porosity numerical model based on the volume-averaged Reynolds-averaged Navier-Stokes (VARANS) equations was validated by comparing the present results with those of previously reported and present experiment. e validated numerical model was then used for extended research.

Numerical Model
In the present research, the resistance-type porosity model proposed by Zhao [17] was adopted. is numerical model used the VARANS equations to uniformly describe the flow inside and outside the porous structure, as well as a volumeaveraged k − ε model to simulate the turbulence effect. In addition, the three-step finite element method (FEM) was applied to solve the VARANS equations numerically, and the computational Lagrangian-Eulerian advection remapvolume of fluid (CLEAR-VOF) method was applied to capture the free water surface. In the following sections, the governing equations, calibration of the resistance coefficients, and validation process were explained. Figure 1 shows the definition of a porous medium. e surface S for generating the averaging volume V included a solid phase and a fluid phase. S was defined by a circle with the radius r 0 . e actual volume of the fluid phase V f varied with the porous medium based on the position of the averaging volume, while the total volume V was constant. e macroscopic length scale L and the pore length scale l were also shown in this figure.

VARANS Equations.
A volume averaging procedure was applied to convert Reynolds-averaged Navier-Stokes (RANS) equations to VARANS equations, and the length scale constraint was given as l << r 0 << L. e superficial volume average denoted by 〈 〉 was defined as the average over the entire volume. e intrinsic average denoted by 〈 〉 f was defined as the average over the fluid volume. e superficial and intrinsic volume averages of a scalar, vector, or tensor denoted by B were defined as follows: (1) e relation between these two averages was as follows: where n was the porosity given by In order to obtain VARANS equations, Navier-Stokes equations were ensemble averaged by applying a time-average operator to obtain the Reynolds-averaged Navier-Stokes (RANS) equations. en, a volume-average operator was applied to the RANS equations to obtain volume-averaged Reynolds-averaged Navier-Stokes (VARANS) equations:  Mathematical Problems in Engineering where f i along with the inertia term was the extended Darcy-Forchheimer resistance term F i defined as follows: e first term on the right-hand side of equation (6) was the linear resistance term, and the second term was the quadratic resistance term. a and b were experimental coefficients that need to be determined in advance. e last term was the inertia term. c m � 0.34(1-n)/n was the added mass coefficient representing the transient interaction between grain and water. It was obvious that n � 1 and c m � 0 outside the porous media, and then VARANS equations automatically changed back to RANS equations.
In equation (5), 〈] t 〉 � C μ 〈k〉 2 /〈ε〉 was the volumeaveraged eddy viscosity, and the balance equations of 〈k〉 and 〈ε〉 were obtained by taking the volume averaging over the standard k and ε equations which turn to the following: where In equations (7) and (8), σ k � 1.00, σ ε � 1.30, C ε1 � 1.44, and C ε2 � 1.92 as Nakayama and Kuwahara [18] (1999) suggested. k ∞ and ε ∞ were additional sources of turbulence due to the presence of porous materials. ey represented small-scale turbulence that was smaller than the averaging volume. Nakayama and Kuwahara [18] also proposed closure schemes for small-scale turbulence based on the numerical simulation of a fluid flowing through a series of square rods:

Calibration Methods for Experimental Coefficients in VARANS Equations.
According to Forchheimer's research, the Darcy-Forchheimer resistance coefficients a and b occurred in equation (6) were related to the permeability of the porous medium. ey were expressed as a � μ/ρK and b � C F / �� K √ , respectively, where μ was the coefficient of kinematic viscosity of the fluid, K was the permeability of the porous medium, which could be written as K � n 3 d 2 50 /α(1 − n) 2 , and C F was an unknown coefficient, which could be written as . us, a and b could be expressed as a � αμ(1 − n) 2 /ρn 3 d 2 50 and b � β(1 − n)/n 3 d 50 equivalently, where d 50 was the mean diameter (or grain diameter if the porous medium was homogeneous) of the porous medium and both α and β were coefficients related to the microscopic geometry of the porous medium. When numerically simulating a specific experimental model, the mean diameter d 50 and porosity n could usually be measured. However, α and β were unknown, and they had to be determined through certain calibration procedures.
It should be noted that research studies gave different values depending on how the coefficients were determined or calibrated, how the resistance coefficients a and b were formulated, and even how the numerical model handled turbulence effects. Zhao [17] proposed a calibration method, that was, when using different values of α and β, the wave damping rate of wave motion over porous seabed could be regarded as an index to evaluate the accuracy of numerical simulation. Although the wave motion over the porous seabed was a dynamic process, after a period of time, the entire damping system became stable, and wave heights along the distance remained unchanged. e total error could be expressed as follows: where k exp i was the wave damping rate measured in experiments; k num i was the wave damping rate calculated by numerical modeling; and ε was the relative error. It was seen that the wave damping rate was the only parameter required in the error estimation procedure, which made this calibration method proposed by Zhao [17] easier to implement. In this paper, the same calibration method was adopted and the coefficients α and β were determined to be 100 and 2.0, respectively.

Validation Using the Existed Model
Experiments. e numerical model proposed by Zhao [17] could be successfully used to simulate the porous media-fluid coupled flow. By comparing the simulation results with the existing model experimental results, the applicability and accuracy of the simulation of the irregular wave interaction with a perforated caisson sitting on a rubble-mound foundation were validated.
Xing et al. [15] studied the irregular wave interaction with a perforated caisson breakwater on a rubble-mound foundation. In their model experiments, the widths of the wave chamber b c were set to 0.15, 0.20, and 0.30 m. e front wall of the caisson was partially perforated from 0.2 m below the water surface to the top of the caisson. e porosity μ of the perforated front wall was set to 20% and 40%. e water depth d during the test was maintained at 0.4 m. e heights of the rubble-mound foundation h m were set to 0.10 m, 0.15 m, and 0.20 m. e berm width of the foundation W was 0.25 m. e core of the rubble-mound foundation was the same as the rock-fill core in the wave-absorbing chamber. e grain diameter d 50 was about 1.1 cm, and the porosity of the core n was about 0.4. e significant wave height H s was 0.053 m. Because the water depth d during the test was maintained at 0.4 m, the significant wave lengths L s were only related to the wave periods. ey were 2.35 m, 1.82 m, and 1.44 m, respectively, when the significant wave periods T s were 1.38 s, 1.15 s, and 0.99 s. After combining all these varied factors, a total of 54 cases were carried out in their work.
As shown in Figure 2, a 2-D numerical wave flume was developed in this study. e length of the flume is 15 m. e left boundary represents the wave-maker, and the boundary condition is the velocity inlet. e rubble-mound foundation and the perforated caisson sitting on it were located in the middle section of the wave flume. e perforated front wall was 7.7 m away from the left boundary, which made the left side of the perforated wall the seaward side and the right side the leeward side.
In order to eliminate the effect of re-reflect waves between the structure and the wave generator, a time-domain method for separating the incident waves and the reflected waves (SIRW method) by Frigaard and Christensen [19,20] was adopted in the present research. It was proved that this method could be used in investigating the irregular wave interaction with a perforated caisson by Tang [21]. In addition, a sponge layer was located on the right side of the wave flume, and the flow velocity in it decayed exponentially.
An irregular wave train could be treated as a stationary stochastic process. It was assumed to be a linear superposition of large numbers of monochromatic wave components with different wave heights, wave angular frequencies, and random initial phases. e incident wave spectrums of experiments and numerical modeling were both the modified JONSWAP spectrum (Goda, 1999 [22]). So, by neglecting the interaction of different component of waves, the reflection coefficient of each component could be calculated by estimation of incident and reflected waves [23]. Once the reflection coefficient of each component was determined, the reflected spectral density was obtained by where K rj was the reflection coefficient of the jth component and S r (ω r ) and S r (ω j ) were the spectral densities of the reflected and incident waves, respectively. For irregular waves, the average reflection coefficient Kr can be calculated by where m 0,r and m 0,I are the zero-order moments of the reflected and incident wave spectra, respectively, and can be calculated by For consistency, the average reflection coefficient Kr was considered as the reflection coefficient of irregular waves and just be written as Kr in the following sections.
In order to accurately calculate the reflection coefficient of each component of waves, five wave height gauges (numbered 1 # to 5 # ) were arranged 2.8 m, 3.1 m, 3.6 m, 3.91 m, and 4.17 m from the seaward face of the caisson, respectively. Appropriate combination of wave gauges was to be selected according to different wave length of each component of waves.
e accuracy of the numerical model could be illustrated by comparing the time series data of wave elevation or the frequency spectrum based on numerical modeling and experimental results. For instance, Figure 3 shows a target spectrum of the irregular wave in experiments and numerical modeling for L s � 1.82 m; Figure 4 shows both the numerical and experimental time series data of wave elevation and the corresponding frequency spectrum measured by the 1 # wave height gauge for h m � 0.15 m; b c � 0.2 m; μ � 0.4; and L s � 1.82 m. It was seen that the numerical results were in good agreement with the experimental results. Only in the high frequency region, some deviation occurred.
In Figure 5, the reflection coefficients Kr calculated based on the numerical and experimental results were compared, where the horizontal axis represented the relative waveabsorbing chamber width b c /L s . It was seen that in Xing's model experiments, the relative wave-absorbing chamber width was about 0.064-0.208, and within this range, K r varied nonlinearly. In some cases, as shown in Figure 5(f ), the minimum of K r was found. erefore, a quadratic function was used to describe the correlation between K r and the relative wave-absorbing chamber width. e curves for the three h m values were calculated from the formula proposed by Xing et al. [15].
In general, Figure 5 indicated a reasonable agreement between the numerical and experimental values of K r in all 54 cases. e numerical model used in the present research was applicable and sufficiently accurate to simulate the wave interaction with a perforated caisson sitting on a rubblemound foundation.

Effect of the Wave-Absorbing Chamber Width on the Wave Attenuation Performance
As mentioned in the previous section, since the range of the wave-absorbing chamber width was very small, within this range, it was reasonable to use a quadratic function to describe the variation of reflection coefficients. However, when the wave-absorbing chamber width continued to increase, what would be the behavior of the reflection   Mathematical Problems in Engineering coefficient and even the horizontal wave force acting on the perforated caisson sitting on the rubble-mound foundation?
In this paper, model experiments and numerical simulations were used to answer this question.

Model Experiments and Numerical
Simulations. e model experiments were carried out at the State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology. e wave flume was 56 m in length and 0.7 m in width. A hydraulic servo push-plate type wave generator was installed on one side of the flume. An antislope device was installed on the other side to dissipate wave energy. e setup of the wave height gauges was the same as in Xing et al. [15]. e roofless caisson was made of PMMA with a thickness of 1 cm. Four rows of uniformly distributed rectangular holes were cut out from 0.2 m below the water surface to the top of the caisson, which led to a perforation ratio μ of 40%. e height of the rubble-mound foundation h m was set to 0.20 m. In these experiments, the range of wave-absorbing chamber width b c was very large, measuring 0.15, 0.20, 0.30, 0.40, 0.50, 0.60, 0.70, and 0.80 m. e perforated caisson model is shown in Figure 6.
In order to calculate the wave forces on the caisson, pressure gauges were installed on both sides of the perforated front wall and the solid rear wall, as shown in Figure 7. Once the time series data of pressure were collected, the horizontal wave forces acting on each wall and the total horizontal wave forces acting on the entire structure could be calculated by a simple integration algorithm.
In the whole experiment, the significant wave height H s was 0.053 m, and other geometric parameters and wave experimental conditions are shown in Table 1.
e numerical model validated in the previous section was adopted here to simulate the model experiment cases. Mathematical Problems in Engineering in a decrease ⟶ increase ⟶ decrease manner of the reflection coefficients. is indicated that for incident waves with the same periods, as the wave-absorbing chamber width increased, the phase differences between the incident wave, the wave reflected from the perforated front wall, and the solid rear wall oscillated. is oscillation behavior is most clearly seen in Figure 8(c) in which the significant wave length L s was the shortest, and the range of the relative wave-absorbing chamber width b c /L s was therefore largest. According to b c    Figure 8, for a greater range of the wave-absorbing chamber width changed, a quadratic function could no longer describe the correlation between the relative wave-absorbing chamber width and reflection coefficients; therefore, the formula of Xing [15] did not apply in such cases.

Horizontal Wave Forces.
Another important subject for studying the wave attenuation performance of a perforated caisson sitting on a rubble-mound foundation was the horizontal wave force. is subject included the horizontal wave forces acting on each wall of a perforated caisson and the total horizontal wave force acting on the entire structure. For a typical Jarlan-type perforated caisson, the total horizontal wave force F h could be decomposed into three different components, as shown in Figure 9(a). e total horizontal wave force could be expressed as F h � F h1 − F h2 + F h3 , where F h1 represented the wave force acting on the seaward side of the perforated front wall, F h2 was the force acting on the leeward side of the perforated front wall, and F h3 was the force acting on the solid rear wall.
Under the same conditions, F h0 was the horizontal wave force acting on a traditional caisson, as shown in Figure 9(b).
To study the horizontal wave forces acting on each part of the perforated caisson, dimensionless ratio parameters F h1 /F h0 , F h2 /F h0 , and F h3 /F h0 were introduced. e dimensionless total horizontal wave force could be expressed by F h / F h0 . Since F h2 was opposite in direction to F h0 , F h2 /F h0 was always negative. Figure 10 shows the numerical and experimental horizontal wave forces as the wave-absorbing chamber width varied from 0.15 m to 0.80 m. In Figures 10(a)-10(c), three different groups of cases with the significant wave lengths L s � 2.35 m, 1.82 m, and 1.44 m were shown, respectively. It could be seen that the numerical results were in good agreement with the experimental results. Based on the numerical results, spline curves are used to smooth the curves in Figure 10 and in other graphs in the following sections.
As the wave-absorbing chamber width increased over a large range, results showed that resonant conditions occurred for various b c /L s values introducing an oscillation in a decrease ⟶ increase ⟶ decrease manner of all F h1 / F h0 , F h2 /F h0 , and F h3 /F h0 . e oscillation of curves of F h1 /F h0 and F h2 /F h0 indicated that for incident waves with the same periods, as the wave-absorbing chamber width increased, the phase differences between the incident wave, the wave reflected from the perforated front wall, and the solid rear wall oscillated. e oscillation of F h3 /F h0 curve indicated that for incident waves with the same periods, as the wave-absorbing chamber width increased, the phase differences between the incident wave, the wave reflected from the solid rear wall, and the re-reflected wave from the leeward side of the perforated front wall oscillated. e curve of F h2 /F h0 showed the most obvious oscillation, and the oscillating amplitude of the F h1 /F h0 curve was smaller than that of F h2 /F h0 . is indicated that the wave forces acting on the seaward and the leeward side of the perforated front wall of the caisson varied asymmetrically with the increase in the wave-absorbing chamber width.
It must be noted that synchronization could be found when investigating the curves of F h1 /F h0 , F h2 /F h0 , and K r .
is meant that when K r reached its maximum or minimum value, F h1 /F h0 and F h2 /F h0 also simultaneously reached their maximum or minimum values. On the other hand, the F h3 / F h0 curve showed some hysteresis. Figure 11 shows the numerical and experimental values of the total horizontal wave force as the wave-absorbing chamber width varied from 0.15 m to 0.80 m. In Figures 11(a)-11(c), three different groups of cases with the significant wave lengths L s � 2.35 m, 1.82 m, and 1.44 m were shown, respectively. Some differences can be seen but generally the numerical results were consistent with the experimental results.
Results showed that resonant conditions occurred for various b c /L s values introducing an oscillation in a decrease ⟶ increase ⟶ decrease manner of F h /F h0 . Compared with the curves of F h1 /F h0 and F h2 /F h0 , the F h /F h0 curve was more similar to the F h3 /F h0 curve. is was because, on the one hand, the horizontal wave force acting on    Figure 10: Horizontal wave forces acting on each part of the perforated caisson with irregular waves: black ○, numerical results of F h1 /F h0 ; black •, experimental results of F h1 /F h0 ; black □, numerical results of F h2 /F h0 ; black ■, experimental results of F h2 /F h0 ; black △, numerical results of F h3 /F h0 ; black ▲, experimental results of F h3 /F h0 . the solid rear wall was always larger than that acting on the perforated front wall under experimental conditions. On the other hand, the wave force acting on the seaward side of the perforated front wall was offset by that acting on the leeward side, so the contribution to the total horizontal force was rather small. In general, for all values of reflection coefficient K r , the total horizontal wave force or its components, results showed that resonant conditions occurred for various b c /L s values introducing an oscillation in a decrease ⟶ increase ⟶ decrease manner. It was clear that their relationships could not be described by a quadratic function.

Extended Research Using the Numerical Model.
e experimental cases mentioned in the previous section were limited. In fact, this experiment only included the cases of rubble-mound foundation height h m � 0.2 m and the perforation ratio μ � 0.4. In other cases, how would the reflection coefficient K r , the total horizontal wave force, and its components vary as the wave-absorbing chamber width increased? Would there be any similarities or differences? To answering these questions, the numerical model was used for extended research to extrapolate the above results beyond the parameters of the model experiments.
Different cases with three rubble-mound foundation heights, two perforation ratios, eight wave-absorbing chamber widths, and three significant wave lengths were considered. It totally included 144 cases after combing all these factors. Table 2 showed all the modeling conditions and geometric parameters.

Effect of Wave-Absorbing Chamber Width on Reflection
Coefficients. Figure 12 shows the obtained numerical and experimental reflection coefficients as the wave-absorbing chamber width varied from 0.15 m to 0.80 m. ree different groups of cases with the significant wave lengths L s � 2.35 m, 1.82 m, and 1.44 m were shown, respectively. In Figures 12(a)-12(c), the perforation ratio was 0.4. In Figures 12(d)-12(f ), the perforated ratio was 0.2.
Generally, as the wave-absorbing chamber widened, the reflection coefficients tended to oscillate in a decrease ⟶ increase ⟶ decrease manner.
is was a common feature regardless of the rubble-mound foundation height or the perforation ratio.
From Figures 12(a)-12(c), it was found that when the perforated ratio μ � 0.4, the oscillating amplitude became smaller for higher foundations. However, in Figures 12(d)-12(f ), this phenomenon was not evident.
Comparison between Figures 12(a), 12(b), and 12(c) with 12(d), 12(e), and 12(f ) in that order showed that as the wave-absorbing chamber width increased, different perforation ratios also resulted in different amplitudes of the oscillating curves. Generally, as the wave-absorbing chamber width increased, F h1 /F h0 , F h2 /F h0 , and F h3 /F h0 tended to oscillate in a decrease ⟶ increase ⟶ decrease manner.

Effect of Wave-Absorbing Chamber
e F h2 /F h0 curve oscillated most obviously, and the oscillating amplitude of F h1 /F h0 was smaller than that of F h2 /F h0 . Synchronization occurred between F h1 /F h0 , F h2 /F h0 , and K r to some extent that they reached their maxima and minima simultaneously. However, F h3 /F h0 showed some hysteresis. ese features were seen regardless of rubble-mound foundation height or the perforation ratio.
For different perforation ratios, it could be found that when μ � 0.4, the horizontal wave force acting on the seaward side of the perforated front wall was smaller than that for μ � 0.2. Conversely, the horizontal wave force acting on the solid rear wall was greater. is was mainly because the perforation ratio of the front wall determined how much incident wave energy passed into the waveabsorbing chamber and how much was reflected from the front wall. When the perforation ratio was small, less wave energy entered, which meant that the perforated front wall was subject to greater horizontal wave forces. In addition, compared with the case of μ � 0.4, the curves of F h1 /F h0 and F h3 /F h0 were flatter near their minima and change more sharply when they began to climb or fall. When μ � 0.2, the curves were less symmetrical than for μ � 0.4.
By comparing figures with different rubble-mound foundation heights, it could be found that the rubble-mound foundation height had limited influence on the horizontal wave forces. When h m became higher, F h1 /F h0 gets smaller only on the seaward side of the perforated front wall. is indicated that the horizontal wave forces acting on the area close to the bottom of the caisson were very small. In addition, some deviation was found in the F h2 /F h0 , curves when h m � 0.2 m especially near their maxima. Figure 14 shows the numerical and experimental total horizontal wave forces as the wave-absorbing chamber width varied from 0.15 m to 0.80 m. ree different groups of cases with the significant wave lengths L s � 2.35 m, 1.82 m, and 1.44 m were shown, respectively. In Figures 14(a)-14(c),      In Figure 14, as the wave-absorbing chamber width increased in a larger range, F h /F h0 tended to oscillate in a decrease ⟶ increase ⟶ decrease manner, but the oscillation was weaker than for either F h1 /F h0 or F h2 /F h0 . e F h / F h0 curve showed a certain amount of hysteresis compared with that of K r . ese were the common features regardless of the rubble-mound foundation height or the perforation ratio.
In addition, the variation of perforation ratio had almost no effect on the F h /F h0 curve although the F h /F h0 curve appeared to resemble the F h3 /F h0 curve more closely as the perforation ratio increased.
is was mainly due to the greater force on the solid rear wall of the perforated caisson with greater perforation ratio.
Comparing the figures for different rubble-mound foundation heights, it could be found that the rubble-mound foundation height had limited influence on the total horizontal wave force. As in the figures for K r , when the perforation ratio μ � 0.4, the oscillation amplitudes became smaller with greater foundation heights, but this phenomenon was no longer evident when μ � 0.2.

Conclusions
Model experiments were carried out to study the effects of a large range of the wave-absorbing chamber width on the wave attenuation performance of a perforated caisson sitting on a rubble-mound foundation. Meanwhile, a porosity-type numerical model based on VARANS equations was validated by comparing with the existing and present model experiments. e validated numerical model was then used for extended research. e conclusions are given as follows: (i) e resistance-type porosity model used in the present research took into account the porosity of the rubble-mound foundation and the fluid viscosity. e numerical results were in good agreement with the existing and present model experiments. erefore, the numerical model could be used to further study the wave interaction with a perforated caisson sitting on a rubble-mound foundation.
(ii) Resonant conditions occurred for various relative wave-absorbing chamber width introducing an oscillation in a decrease ⟶ increase ⟶ decrease manner of the reflection coefficients. is was a common feature regardless of the rubble-mound foundation height or the perforated ratio. (iii) Resonant conditions occurred for various relative wave-absorbing chamber width introducing an oscillation in a decrease ⟶ increase ⟶ decrease manner of the horizontal wave forces acting on each part of the caisson. e horizontal wave force acting on the perforated front wall varied synchronously with the reflection coefficients, but the horizontal wave force acting on the solid rear back wall showed a certain amount of hysteresis. (iv) Resonant conditions occurred for various relative wave-absorbing chamber width introducing an oscillation in a decrease ⟶ increase ⟶ decrease manner of the total horizontal wave force acting on the caisson. A certain amount of hysteresis was evident in the variation of the total horizontal wave force compared to that of the reflection coefficient. (v) With the variation of the wave-absorbing chamber width, the reflection coefficient and horizontal wave force acting on a perforated caisson sitting on a rubble-mound foundation had minimum and maximum values. Considering the attenuation performance and wave force conditions, the waveabsorbing chamber should be neither too narrow nor too wide.

Data Availability
All the useful data including the experimental and numerical results used to support the findings of this study are included within the article. Other data are available upon request to Peihong Zhao (571509076@qq.com).

Conflicts of Interest
e authors declare that they have no conflicts of interest.