Living Breakwaters have evolved from traditional breakwaters to create multi-purpose structures that provide environmental and social benefits. These breakwaters promote a healthy habitat for flora and fauna, while achieving the same structural capabilities as conventional structures. Reefy has developed a modular living breakwater that can serve both as an artificial reef and a stable breakwater. This study aims to model the impact of a Reefy breakwater on the transmitted wave height, using a process-based numerical model, namely XBeach. This study builds upon a previous one conducted by van den Brekel [2021], where scaled experiments were performed in a physical wave flume, to investigate the hydrodynamic and ecological functionalities of the Reefy breakwater. In the wave flume experiments a total of 15 structure configurations together with 35 wave conditions, both regular and irregular, were tested. The first 7 of the configurations were 2D structures with a relatively simple shape and a porosity of 20%, and the rest were complex 3D structures with a porosity above 20%. The wave flume was recreated within the numerical model. The influence of the breakwaters on the wave heights are examined through the comparison between the calculated transmission coefficient of the wave flume and the numerical model. However, XBeach, and more specifically XBeach non hydrostatic+ mode which was used, does not resolve the complex interaction between waves, friction and porosity of the structures. As a result, a simple method had to be found to compensate for the loss of these phenomena, while successfully calculating the transmission coefficient. The findings of this analysis demonstrate that treating the breakwater as an impermeable change in bathymetry, coupled with a 15% reduction in structure height, produces the desired outcomes. From the 35 wave conditions only the 8 irregular wave conditions were chosen to be of interest due to time constraints. Among these 8 wave conditions that were used as an input in XBeach, only 5 of them were satisfying the criterion that describes the range for which the numerical model should produce accurate results (kd<2). The transmission coefficient was defined as the ratio of the transmitted wave height behind the breakwater, and incident wave height in front of the structure. The wave decomposition method used was a modified Guza split method. In pursuit of finding a valid solution, the breakwaters were modeled firstly as impermeable structures with a decreased height, secondly as impermeable structures with a decreased width, thirdly as vegetation with the help of the vegetation module, and lastly the maximum wave steepness criterion as applied in the model was increased (maxbrsteep up to 1.4), without modifying the structure height. Among the suggested solutions and after being validated for 13 experiments with both simple and complex structures, decreasing the width did not achieve the desired results while the vegetation module failed to produce consistent outcomes. As far as transmission coefficient is of interest, a 15% decrease in the structure height had the smallest mean absolute percentage error of almost 10% and a root mean square error of 8.2%. On the other hand, choosing to increase the maximum wave steepness (to a value of 1.4) is not a valid option for the complex 3D structures, but illustrates even better behavior than a change in the structure height when simple structures with 20% porosity are concerned. To summarize, this thesis provides alternatives to replicate in a simplified way the impact of porous structures, such as Reefy breakwaters, without the need of an extensive and time expensive numerical model.