This study investigates the capacity of an open inverted filter to stabilize the interface between the sand fill and core fill in a rubble mound breakwater. The strength of this filter relies on the arching mechanism of the base material, which is directly influenced by factors such as stability ratio and effective stresses. These forces arise from porous flow, and its quantification involves both parallel and perpendicular hydraulic gradients. A full-scale physical model, with a 760 mm open inverted filter, is employed to determine the critical hydraulic gradients required to initiate erosion. The model consists of three compartments, with the central compartment containing the open inverted filter. The two side compartments are filled with water, with one of them featuring a plunger mechanism to induce flow. Pressure sensors and water level gauges are used to measure the hydraulic gradients. In addition to the physical model, a numerical model is developed to gain insights into the flow dynamics within the system. The tests reveal that an increase in effective stress results in higher critical parallel and perpendicular gradients. Furthermore, the addition of fines at 5% or 10% seems to enhance the critical parallel and perpendicular gradients. A negative correlation was found between SR and 𝑖⊥,𝑐𝑟𝑖𝑡, the results were fitted to a hyperbolic function. The findings further indicate that, even after erosion occurs, a stable interface can be reestablished. Therefore, designing solely for the single largest wave is not imperative.

On the east coast of Romania, at Eforie, coastal erosion manifests. To strengthen the coastal area a large coastal protection project was setup involving beach nourishment combined with the construction of breakwaters. The breakwaters are designed with the well-known modified Van der Meer formulas. To ensure confidence in the breakwaters stability the designs were tested during a 3D physical model test. The designs were stable, however the measured damage to the breakwater was larger than was expected. Investigating the physical model test resulted in three important aspects, that need to be investigated. First, a broad list of configurations, that could responsible for the damage, could be identified. Secondly, excessive wave breaking on the foreshore during the model test resulted in a shift in wave energy from higher to lower frequency waves. These low frequency waves, or so-called infragravity waves became increasingly dominant near the coast. Their effect on the stability of the structure is however unknown. Thirdly, the incoming waves were very oblique and on a shallow foreshore, resulting in a different failure mechanism than is incorporated in the modified Van der Meer formulas. A new stability method has to be formulated to investigate breakwater stability for this physical model test. Due to the broad list of configurations and the possible influence of infragravity waves, real scale or physical model tests are not feasible and a numerical model needs to be used. A reliable and widely applicable method to link breakwater stability to a numerical model has however not yet been composed. Therefore, the goal of this thesis is to propose a method that can link breakwater stability to a numerical model and to assess whether the identified configurations resulted in the damage along the statically stable rubble mound breakwater, measured during the Eforie physical model test. First, an investigation is performed on the physical model test set-up and observations, resulting in a final list of 5 configurations, that are investigated in this research: The applied offshore transitional slope (1), the assumption of uni-directional waves (2), the slope of the lower foreshore (3), the depth-contour lines inducing wave focusing (4) and the very oblique wave angle on a shallow foreshore (5). Secondly, a method is proposed linking breakwater stability to a velocity signal from the numerical model SWASH. An equation is formulated, based on the theory of Izbash (1935), with a slope factor included, and scaled with the theory of Shields (1936). It requires a velocity signal, that can be obtained from SWASH, to calculate a stone size required for stability. Thirdly, a numerical model is set up in SWASH, with grid dimensions 3m x 2m, resembling the physical model test. The breakwater is modelled as an impermeable core with a permeable porosity layer placed on top. The thickness of the porosity layer is based on the thickness of the outer armour layer of the original breakwater. The numerical model is validated by comparing the wave characteristics, at several locations along the breakwater, to wave data available from the physical model test. The numerical model shows accurate resemblance of the wave characteristics. Since the wave velocity is linked to the wave height, it is assumed that the wave velocity on the breakwater is also correctly modelled. The model is therefore found valid for the modelling study. In the numerical model along the still waterline measurement points are indicated that provide the velocity and waterlevel signal during a simulation. In the numerical model two layers in the vertical are assumed and tested to be sufficient. The velocity of the top layer resembles the velocity that flows just over the stones. Therefore from the velocity signal of the top layer the governing u_0.2% along the waterline at the breakwater is obtained and from the waterlevel signal the wave spectrum is derived. In the study simulations are performed in which the configurations are tested one by one, and all simulations are assessed on two parameters: the u_0.2% and the wave spectral transformation along the breakwater. The results from the different simulations are compared relatively to identify the relative effect the configurations have on the velocity and wave characteristics. The results of this research show that breakwater stability can be predicted reasonably well from a velocity signal obtained from SWASH. The velocity signal, obtained from SWASH, results in reliable stone sizes. The configurations could be investigated with the proposed method and the results provide reliable and useful insights. In addition, the proposed method is able to identify the effect of infragravity wave energy on the stability of a breakwater. The method is also tested by calculating the relative obliqueness factors for different incoming wave angles, which shows promising results. It is important to reproduce the breakwater porosity well in the numerical model as it can significantly influence the velocity signal. A decrease/increase of the porosity thickness with 30% or 0.6m can result in an increase/reduction of 20-26% in velocity respectively. The five discussed configurations provide partial explanations to (in)directly induce the higher breakwater damage in the physical model test. Both the applied offshore transitional slope (1) as the assumption of uni-directional waves (2) result in a slight underestimation of the breakwater stability and therefore a somewhat conservative design along the entire length of the breakwater. The combined effect resulted in a reduction of 0-6% around the head of the breakwater, h/Hs = 2.5-4.8, and a reduction of 16-24% near the shore, h/Hs = 1.1-1.8. Especially near the shore the breakwater is conservatively designed, due to the fact that both the transitional slope as the assumption of unidirectional waves increases the infragravity wave energy in the system. It is found that incoming waves break around h/Hs = 1.7 after which the infragravity waves induce a temporary increase in waterlevel, around h/Hs = 1.1-1.8. This allows the depth-limited short waves to become bigger resulting in higher velocities and more damage on the breakwater. This affects the breakwater stability closer to shore and needs to be taken into account when designing a breakwater in these conditions. The lower foreshore (3) induces the generation of infragravity waves, which affect the velocity closer to shore as described above. The depth-contour lines (4) result in a wave focusing effect increasing the velocity around h/Hs = 1.1-1.8 with 8-11%. Based on the results of this thesis the very oblique wave angle on a shallow foreshore (5) does not induce higher velocities and breakwater instability. It is however assumed that the effect of a breaking plunging wave, inducing acceleration and pressure difference effects on the stones on a slope, is not sufficiently into account, due to the grid dimensions used in the model. As other plausible causes of the increased damage are disproven, it seems likely that the different oblique wave breaking process that is not modelled in detail leads to the increased damage.

This study investigates the use of an open filter in a land reclamation structure. The possibility of replacing the geotextile that secures the core-backfill interface with a granular filter has been investigated. This is easier to install and therefore a potential cost-saving. Experiments have been done to determine the influence of various parameters on the interface stability of a geometrically open inverted granular filter. Relationships have been found which can assist in designing such a filter.