1D-2D Numerical Model for Wave Attenuation by Mangroves as a Porous Structure

In this paper, we investigate wave attenuation caused by mangroves as a porous media. A 1-D mathematical model is derived by modifying the shallow water equations (SWEs). Two approaches are used to involve the existing of mangrove: friction term and diffusion term. The model will be solved analytically using the separation of variables method and numerically using a staggered finite volume method. From both methods, wave transmission coefficient will be obtained and used to observe the damping effect induced by the porous media. Several comparisons are shown to examine the accuracy and robustness of the derived numerical scheme. The results show that the friction coefficient, diffusion coefficient and vegetation’s length have a significant effect on the transmission coefficient. Moreover, numerical observation is extended to a 2-D SWEs, where we conduct a numerical simulation over a real bathymetry profile. The results from the 2-D numerical scheme will be validated using the data obtained from the field measurement which took place in Demak, Central Java, Indonesia. The results from this research will be beneficial to determine the characteristics of porous structures used for coastal protection.


Introduction
The ocean covers 71% of the Earth's surface and 40% of the world population, approximately 2.4 billion people, live within 100 km of the coast. On the other hand, more than 600 million people, live in low-elevations coastal areas-less than 10 m above sea level-(Ocean Conference, United Nations, 2017). This causes a severe vulnerability to the sea level rise caused by global warming. It is predicted that the global sea-level rise will reach 20-30 cm by 2050 [1][2][3][4][5][6][7].
Consequently, coastal protection is essential to overcome coastal problems and minimize risks. There are many approaches to protect coastal such as using breakwaters, sea walls, tidal barriers, artificial headlands, and many more [8][9][10]12,13,33]. However, these kinds of coastal protections can be harmful to the ecosystems around them. On that account, this research will observe the phenomenon using porous media, particularly mangroves. Mangroves have many ecological, social, and economic advantages. Mangroves provide a good ecosystem, enhance coastal accretion, cause a considerable wave damping, and decrease flow velocities regarding high tides or flooding [14]. The wave attenuation phenomenon by the porous media can be described as two events: friction and diffusion. The interaction between the wave and the porous media will cause friction. On the other hand, as the wave reaches and passes through the porous media, its energy will be spread, which in this case is called diffusion. simulations will be compared with the analytical solutions and experimental data. Lastly, section six will conclude this paper and give some recommendations for future researches.

Mathematical Model
Here, we will formulate the 1-D and 2-D mathematical model based on the shallow water equations (SWEs). The SWEs are hydrostatic model meaning the wave velocity along the vertical direction is homogenous. Therefore, we assume that the effects of vertical shear of the horizontal velocity are negligible. The "no-shear" assumption is reasonable if the fluid is "shallow", and this partly accounts for the name "shallow water equations" [53]. SWEs are derived from the principles of two partial differential equations. The first is the mass conservation equation, and the other is the momentum equation. In order to capture the wave attenuation phenomenon, the SWEs will be modified.

One-Dimensional (1-D) Mathematical Model
For the 1-D model, we divide the domain of observation into five areas as depicted in Figure 1. The incoming wave enters the observation domain from the left side with the initial amplitude A i moving towards the right-side with the transmitted amplitude A t . In this 1-D model observation, we consider the bottom profile to be flat, with water depth d. Consider the wave profile for the 1-D model is a function of space x and time t. The SWEs are described as follows. In Figure 1, we can see that d(x) represents the depth of the water and z = 0 is the still water level. In this observation, the wave enters the domain by a certain period and frequency with the velocity of u(x, t). It will oscillate as far as η(x, t) which makes the total water thickness is h(x, t) = η(x, t) + d(x).
Next, we modify the momentum equation to capture the phenomenon. As the wave enters domain Ω 2 and Ω 4 , it will interact with the vegetation, which will cause friction. Therefore, the first modification is the addition of a friction factor, c f ωu. This particular friction vector is introduced in [54] to allow evaluation of the wave attenuation due to vegetation friction to be performed from knowledge of periodic wave characteristics. The second modification the addition of a diffusion factor. As we know, the diffusion term itself represents damping phenomenon. Therefore, adding a diffusion factor will allow us to capture the attenuation effect caused by the vegetation. Implementing the modifications in the momentum Equation (2) for each domain, we have with piece-wise constant functions of c f and c D are described as follow and We use the notation c f to denote the friction coefficient and c D to denote the diffusion coefficient, where ω is the wave frequency and g is the gravitational acceleration.

Two-Dimensional (2-D) Mathematical Model
In this subsection, we will modify our model even further by constructing the 2-D model. The 2-D model will be beneficial when observing the phenomenon over a real bathymetry profile. For the 2-D model, the observation domain will be divided into three domains Ω 1 , Ω 2 , and Ω 3 , with Ω 2 is the porous media domain and Ω 1 , Ω 3 are the nonporous media domains. The 2-D wave profile is a function of the longitudinal distance x, span-wise distance y, and time t. The wave velocity along x-direction is denoted by u(x, y, t) and the wave velocity along y-direction is denoted by v(x, y, t). The initial condition of the velocities are u(x, y, 0) = 0 and v(x, y, 0) = 0, with the initial still water condition is η(x, y, 0) = 0. As mentioned before, SWEs are derived from two partial differential equations. Below, Equation (7) represents the mass conservation whereas Equations (8) and (9) are the momentum equations with respect to x and y.
As in the case of the 1-D model, Equations (8) and (9) must be modified to capture the phenomenon. Here, we add the same friction and diffusion factor as previously discussed in this paper. Therefore, the governing modified momentum equations for all the three domains are as follows 1.
Friction term with

Diffusion term
with

Analytical Solution
In this section, we will derive the analytical solution for the 1-D SWEs. The analytical solution is in the form of the transmission coefficient K t which measures the reduction of wave amplitude, hence Here, a i refers to the initial amplitude and a t is the transmitted amplitude. Before we derive the solution for K t , we will construct the solution for η and u.

Solution for η(x, t) and u(x, t)
First, we derive the solution for η and u without the porous media. Recall that we observe the phenomenon in a shallow water environment. Hence, we assume η(x, t) ≤ d(x) which makes h(x, t) ≈ d(x). In this observation, we also assume d(x) is a constant through all domains. By implementing these assumptions to Equation (1), we have Let η, u be the following functions Differentiate Equations (18) and (19) with respect to x and t, then substitute the results to the SWEs represented by Equations (2) and (17), we have iωB 1 (x)e iωt + gA 1x e iωt = 0.
Differentiating Equation (21) over x and substituting it to Equation (20) will give us By using the method of characteristics, we obtain the solution for Equation (22), where i = (−1) , C 1 and C 2 are undetermined coefficients, and Next, substituting the derivative of Equation (23) over x to Equation (21) will result in Substitute A 1 (x) and B 1 (x), respectively, to Equations (18) and (19), we have the solutions for η and u for the domain Ω 1 as follow Let a i = C 1 and a r = C 2 , where a r is the reflected amplitude. We have Note that domains Ω 3 and Ω 5 have the same characteristics as domain Ω 1 . Therefore, by implementing the exact derivations, we have the solution for η and u for domains Ω 3 and Ω 5 are as follows and where C 5 , C 6 , C 9 , and C 10 are undetermined coefficients. Notice that domain Ω 5 is directly connected to the shore and we assume that the wave is absorbed entirely by the shore. This implies that there is no amplitude being reflected from the shore to domain Ω 5 which means a r = 0. Therefore, if we consider C 9 = a t , which is the transmitted amplitude, and C 10 = a r , then we have

Solutions for η(x, t) and u(x, t) with Friction Coefficient
In this subsection, we will construct the solutions for η and u with porous media represented by friction term. Similarly, by substituting Equations (18) and (19) to the SWEs, we have iωA 2 (x)e iωt + dB 2x e iωt = 0, then we substitute B 2 (x) from Equation (37) to Equation (36), so that By using the same characteristics method and procedures as applied before, the solutions for A 2 (x) and B 2 (x) are where Therefore, the solutions for domain Ω 2 are described as follow Furthermore, in a similar way, we have the solutions for domain Ω 4 ,

Solutions for η(x, t) and u(x, t) with Diffusion Coefficient
Similarly, here, we have the solutions for η and u with diffusion factor for domain Ω 2 , while for domain Ω 4 , where In general, the solutions for both η and u for each domain are as follow where . Now, we have successfully derived the solutions for η and u for each domain. Using this knowledge, we can construct the solution for K t . Notice that in reality, wave elevation and horizontal flux are continuous through all domains. As we can see, both η and u have different formula in each domain which is based on whether or not the vegetation is present. Therefore, we must apply η| , where x i is the discontinue points 0, L 1 , L 2 , L 3 , for each domain, to maintain wave continuity.
Let α = 1 + k/k j and β = 1 − k/k j . By implementing the wave continuity in domain Ω 1 and Ω 2 , we have which implies that Next, for domain Ω 2 and Ω 3 , the results are By manipulating Equations (57) and (58), we can express C 5 and C 6 as follow Moving on to domains Ω 3 and Ω 4 , in a similar way, we yield Lastly, for domains Ω 4 and Ω 5 , we have By using all the equations above, we have the final result, where We can see that K t depends on several parameters which are length, friction coefficient c f , and diffusion coefficient c D . The relation between K t and the parameters will be discussed further in Section 5.

Numerical Scheme
In this section, we will discretize our model using the finite volume method on a staggered grid. We will also use an upwind approximation and centered discretization to compute the numerical fluxes. The constructed numerical scheme will be then used for the numerical simulations in Section 5.

One-Dimension (1-D) Scheme
Consider our observation domain to be Ω = [0, L] and the time interval t = [0, T]. The spatial domain is partitioned into N = L ∆x + 1 cells of the length ∆x in a stag-gered way which is x 1/2 = 0, x 1 , ..., x j−1/2 , x j , x j+1/2 , ..., x Nx , x Nx+1/2 = L, as illustrated in Figure 2. Meanwhile, the time interval is divided into Nt = T ∆t + 1 time steps by the length of ∆t. In this numerical scheme, we will evaluate the mass conservation Equation (1) in a cell centered at x j and momentum Equation (2) in cell centered at x j+1/2 . From this setting, we have the horizontal flow η and the water thickness h evaluated in the full grid x j , whereas the wave velocity u is evaluated in the half grid x j+1/2 . Note that u at the boundary x = 0 and x = L are always zero. Therefore, by applying the Forward in Time and Center in Space (FTCS) method, we have the following approximations of the 1-D SWEs: with q n j+1/2 = * d n j+1/2 u n j+1/2 , j = 1, 2, ..., N x , and n = 1, 2, ..., N t . Here, η n j denotes the approximation of η at point x j and time t n where u n j+1/2 is the approximation of u at point x j+1/2 and time t n . Notice that we evaluate η in Equation (67) at t n+1 rather than t n . We do this in order to avoid any instability in the scheme. Additionally, according to the Von Neumann stability analysis, the stability condition of numerical scheme (66)-(67) is 0 ≤ gd 0 ∆t ∆x ≤ 1, with d 0 is the flat bottom depth [55,56]. Note that in calculating q n j+1/2 , we need the information of d n j+1/2 . However, as d shares the same characteristics as h, it is evaluated on the full grid. Thus, we need to approximate the value of d on the half grid, * d n j+1/2 . This is where the upwind approximation will take place. The upwind method will utilize the wave direction to approximate the value of d n j+1/2 as follows * d n j+1 The upwind approximation works by approximating the value of d j+1/2 with the value of d around it. If the wave is moving from left to right, in this case u j+1/2 ≥ 0, then * d n j+1/2 = d n j , which is the d on its left. This implies that the flux in x j+1/2 equals d j u j+1/2 . The same interpretation applies when the wave is moving towards the left side, u j+1/2 < 0, then * d n j+1/2 = d n j+1 , which is the d on its right. Hence, the flux in x j+1/2 equals d j+1 u j+1/2 . Now, as for the modified SWEs, the numerical approximations are as the followings:

Two-Dimension (2-D) Scheme
First, we consider the 2-D shallow water equations with the observation domain Ω = [0, L x ]¯[0, L y ] and the time interval t = [0, T]. The spatial domain is divided using a rectilinear grid of N x¯Ny and the time interval t is divided into N t time steps. With the same approach as in the 1-D scheme, we will evaluate η and h in the mass conservation equations on the full grid x i,j . In contrast, the velocity in the x-direction u and the velocity in y-direction v will be evaluated on the half grid x i−1/2,j and x i,j−1/2 , respectively, as the Figure 3 shown below. Here, the indices are i = 1, 2, ..., N x and j = 1, 2, ..., N y . The approximation of η and h at point x i,j and time t n are represented, respectively, by η n i,j and h n i,j . Now, for u, the approximation at point x i+1/2,j and time t n is denoted by u n i+1/2,j . Meanwhile, at x i,j+1/2 and time t n , the approximation of v is denoted by v n i,j+1/2 , with n = 1, 2, ..., N t . Consider a scheme at time t n+1 , the discretization of the 2-D linear SWEs will yield the following equations where q n u(i+1/2,j) = * d n i+1/2,j u n i+1/2,j and q n v(i,j+1/2) = * d n i,j+1/2 v n i,i+1/2 . The Von Neumann stability analysis for the 2-D linear SWEs provides the stability condition for Equations (71)-(73), which is 0 ≤ ∆t gd 0 ( 1 ∆x 2 + 1 ∆y 2 ) ≤ 1 [57]. Next, we need to approximate the value of d on the half grid x i+1/2,j and x i,j+1/2 . Similarly, using the upwind approximation for all (i, j) ∈ {1, ..., N x } × {1, ..., N y }, we have * d n i+1/2,j = d n i,j , u i+1/2,j ≥ 0, d n i+1,j , u i+1/2,j < 0, For the vertical components, when the flow is moving up, v n i,j+1/2 ≥ 0, we use the information from below, which is d i,j . However, if the flow is moving towards the bottom, v n i,j+1/2 < 0, then the information of the upper flow, d i,j+1 , is used. For the horizontal components, we treat them as described in the 1-D scheme.
As for the modified momentum equations, the discretizations are as follows:

Numerical Simulation and Discussion
Here, we will implement our numerical schemes to simulate the wave propagation phenomena for both the 1-D and 2-D models. The simulations will allow us to observe how porous media affects wave amplitude. Furthermore, we will compare the results with those from an experiment and those of the analytical solutions. For the simulations in this research, all parameters, including the axes in the figures, are in SI units.

One-Dimension (1-D) Numerical Simulations
For the 1-D simulation, we choose 0 < x < L as the spatial domain, with L = 30 and water thickness d = 10, along all x. The porous media will be placed in domains Ω 2 and Ω 4 , where Ω 2 = {x | 5 < x < 10} and Ω 4 = {x | 15 < x < 20}. We consider an incoming monochromatic sinusoidal wave, which enters the observation domain from the left side moving towards the right side. We choose the initial conditions to be the following: u(x, 0) = 0.
For the left boundary condition, we take a monochromatic sinusoidal wave with the initial amplitude A = 0.1 m and wave frequency ω = 2π, and as for the right boundary, we apply the absorbing boundary condition, with g = 9.81, as follows First, we will simulate the wave propagation phenomenon with the friction coefficient. We choose the friction coefficient c f = 0.32, the spatial partition ∆x = 0.2, and the time partition ∆t = ∆x √ gd , in order to maintain stability. The result of the 1-D numerical simulation is presented below. The green areas represent the vegetation (porous media) and the blue line represents the surface wave elevation. From Figure 4, we can see that the wave amplitude is reduced after passing through the vegetation. The reduction of the wave amplitude can be seen more clearly with the numerical transmission coefficient K t computed from the simulation. The numerical simulation gives us the result of K t = 0.3767. This tells us that the porous media reduces the wave amplitude by more than 60% of the initial amplitude. Now, for the diffusion factor, we set the diffusion coefficient, the spatial partition, and the time partition as c D = 0.6, ∆x = 0.23, and ∆t = 0.01, respectively. In Figure 5, we can see the outcome of the 1-D numerical simulation. Notice that after passing through the porous media, there is a slight decrease in the wave amplitude. The numerical simulation gives us the result of the transmission coefficient K t = 0.8166. From both results, we can conclude that both models, with friction and diffusion coefficient, have successfully captured the phenomenon of wave reduction caused by the porous media. However, the reduction caused by friction is more significant than diffusion.
In addition, for this 1-D model, we have performed a simulation to evaluate the computational cost of our model and its comparison against an existing model, which is a Boussinesq-type model. Since we are not studying the Boussinesq-type model for wave-vegetation interaction cases, in this comparison, we only simulate the more simple wave phenomenon, which is standing wave. In the simulations, we use a spatial domain of [0, 250] with observation time of T = 15. All the parameters, including ∆x, ∆t, gravity acceleration, and water depth are set to be exactly the same. After simulating both models, it was found that the computational cost for Boussinesq-type model is about 0.651996 s. Meanwhile, our SWEs give a computational cost of around 0.075074 s, which is almost 10 times smaller than Boussinesq-type model's. This finding confirms that our model is much less expensive than another existing model, in this case, Boussinesq-type model.

Comparing the One-Dimensional (1-D) Numerical and Analytical Result
In this subsection, we will compare the numerical and analytical values of transmission coefficient K t . The numerical K t is computed in the same way as the analytical K t . We compute the numerical K t by comparing the maximum value of η in domain Ω 5 with the maximum value of η in domain Ω 1 . We will compute K t with different parameters of the porous media. This will help us to understand more about the relationship between the porous media and wave amplitude reduction.
The results from the numerical simulations approximate our analytical solutions well when describing the relation between K t with both c f and c D , as depicted in Figure 6. After several simulations, we obtain that the RMSE of the numerical approximations towards the analytical solutions is 3.7 × 10 −5 for friction and 0.874 for diffusion. We also compare the numerical and analytical solutions of the relation between transmission coefficient K t and the vegetation length. We conduct multiple simulations using different vegetation lengths. Similar to the analytical solutions, the numerical solutions give us the result that the relation between the transmission coefficient K t and the length of vegetation is significant. From Figure 7, we can see that the discrepancy between the numerical and analytical transmission coefficient K t is hardly noticeable. This implies that our numerical solutions successfully approximate the analytical solutions with the RMSE of 1.8 × 10 −4 for phenomenon with friction factor and 0.852 with diffusion. From Figure 6 and Figure 7, we can conclude that K t has a strong relation with c f , c D , and the length of vegetation. The transmission coefficient will gradually decrease along with the increase of these parameters. Next, from Figure 8, we can see that K t remains almost constant throughout the changes of vegetation distance. If we take a closer look to the Figure 8, we will see that transmission coefficient K t does oscillate over different vegetation distances. The numerical simulation using friction coefficient c f gives us a clear representation of this phenomenon, as depicted in Figure 9. From these results, we can say that the distance between vegetation patches do not have a significant effect upon K t . The numerical solutions approximate the analytical solutions with RMSE of 2.1 × 10 −5 for friction, whereas for diffusion, the RMSE is 0.765.  We can conclude that the transmission coefficient and the distance of the vegetation do not have a particular relation. In other words, the distance between vegetation is insignificant towards the reduction of wave amplitude by the porous media. The results of the numerical approximations towards the analytical solutions are very satisfactory with the RMSE value close to zero. This implies that our numerical scheme is validated, since it successfully approximates our analytical solutions.

Two-Dimensional (2-D) Numerical Simulations
The 2-D numerical simulation will help us investigate and understand more about the effectiveness of the porous media on the reduction of wave amplitude over an uneven bottom. We have conducted an experiment in Demak and now will compare the experiment results with the numerical results.
Similar to the 1-D numerical simulations, for the left boundaries, we use In this 2-D numerical simulation, we are going to stop the simulation before the wave reaches the dry area to avoid any instability. Therefore, any right boundary for this simulation is applicable.
First, we simulate our 2-D scheme over a flat bottom. This is to make sure that our scheme is stable for the entire 2-D simulation. We will simulate a monochromatic sinusoidal wave entering the observation domain. The observation domain is set to be Ω = {(x, y, t) | 0 < x < 40, 0 < y < 50, t < 7}, ∆x = 0.2, ∆y = 0.2, ∆t = ∆x √ gd , with g = 9.81 and d = 5.
We consider an incoming wave η(0, y, t) = Asin(2.5t) and A = 1 is the wave amplitude. Figure 10 gives us the result of the simulation. Notice that the waveform and amplitude remain the same throughout the entire simulation. Therefore, we can conclude that our 2-D model works perfectly in this setting. Next, we will simulate the wave propagation phenomenon over a flat bottom profile by adding porous media to our simulation. We are going to implement Equations (76)-(79) with the settings of Ω = {(x, y, t) | 0 < x < 20, 0 < y < 20, t < 5}, ∆x = 0.2, ∆y = 0.2, ∆t = ∆x √ gd , g = 9.81, and d = 5. The porous media will take place in Ω 2 = {(x, y) | 5 < x < 10, 0 < y < 50}, with c f = 0.32 and c D = 0.15. Consider a monochromatic sinusoidal wave η(0, y, t) = Asin(ωt)t, with A = 1 and ω = 2π, which enters the observation domain. We have the results of the numerical simulations as follows. From Figure 11, we can see that the addition of the porous media, which represented by the green area, causes the wave amplitude to decrease. The result of the simulation with friction captures this phenomenon perfectly, whereas for diffusion, it is slightly harder to observe in plain view. Another way to observe the reduction is by calculating the transmission coefficient. The numerical transmission coefficient K t computed by our simulation gives us the result of K t = 0.4966 for friction and K t = 0.9574 for diffusion.
Next, we continue to conduct numerical simulation using a real bathymetry profile. In this simulation, the wave will move from the deeper water level towards the shallower water area. The change in water depth will cause some of the wave properties such as length, height, and velocity to change. This phenomenon is called wave shoaling. The interpretation of this phenomenon can be seen in Figure 12, using an artificial bottom profile. As we can see in Figure 12, the wave amplitude increases as the wave approaches the shallow water level. Therefore, in order to measure the wave amplitude reduction by the porous media, we cannot use transmission coefficient K t alone. We need to compare the K t with the shoaling coefficient K s . K s measures the ratio of wave heights in two different water levels. The difference between K s and K t is that K s is calculated without any interference of the porous media. We are going to use the analytical shoaling coefficient as stated in [57], which is described as the following, where H 1 represents the initial wave amplitude at water depth d 1 , whereas H 2 represents the wave amplitude in shallow water level d 2 .
We place the porous media in Ω 2 = {(x, y) | 5 < x < 15, 0 < y < 50}, the initial wave is η(0, y, t) = Asin(15t)t, and the initial amplitude is A = 10. Figure 13 gives us the results of the 2D numerical simulations of wave propagation phenomenon with friction and diffusion coefficient. The wave amplitude reduction was not captured clearly in Figure 13. We can still observe how the porous media reduces the wave amplitude through transmission coefficient K t . The numerical simulation gives us the result of shoaling coefficient K s = 1.0552 and friction K t = 0.5148, whereas diffusion K t = 0.8184. This implies that the porous media reduces the wave amplitude by more than 50% using friction and 23% using diffusion. Therefore, our 2-D numerical model successfully captured the wave amplitude reduction caused by the porous media.

Comparing the Two-Dimensional (2-D) Numerical Result and Real Experiment Data
Here, we will compare the numerical simulation results with real experiment data. The experiment took place in Demak on 6-7 February 2020. In simulating our numerical solution, we will use the 2-D model and real bathymetry profile of Demak as depicted in Figure 13. We retrieved the bathymetry data from GEBCO with longitude of [−6.5205 − 7.0205] and latitude of [110.0058 110.3059].
For the numerical simulation with friction factor, we will set the spatial partition to be ∆x = ∆y = 0.3 and c f = 0.2. Whereas for diffusion, ∆x = ∆y = 0.5 and c D = 1. The friction factor c f and diffusion factor c d applied in the two dimensional (2-D) mathematical model are based on the range of c f and c d implied in the one dimensional (1-D) mathematical model which gave the minimum error. Furthermore, we iterate c f and c d within the given range to fit our actual experiment data in order to achieve the minimum error. The results of both experiment and numerical simulations are shown in Table 1. As we can see in Table 1, our numerical K t approximates the experimental K t well with errors of less than 10% for both friction and diffusion model. This implies that our 2-D numerical scheme successfully captures the phenomenon and estimates wave attenuation by vegetation over a real bathymetry. Thus, we can say that our model is quite applicable to be used in real coastal areas. Therefore, it can be used to design an effective porous media in the form of vegetation as an actual coastal protection against wave attack. Further, from Table 1, we also can conclude that, in this specific case, the presence of vegetation can help reduce the wave amplitude by more than 30%.
In addition, We have also calculated the computational cost of our model for this particular case. The results are different for friction and diffusion model. For model involving friction coefficient, the computational cost is found to be 4.573120 s, while for model with diffusion term, the computational cost is 2.375985 s. Despite the fact that the first model gives a computational cost twice as expensive as the second one, but fairly, both models can be executed in a quite short of time. This finding, supported by the comparison of the computational cost between 1-D SWEs and 1-D Boussinesq-type model mentioned previously, can be used to confirm that our model is indeed less expensive than other previous model, especially Boussinesq-type model.

Conclusions
From this research, we successfully captured and described the wave attenuation by vegetation phenomenon through the wave transmission coefficient K t from both the analytical and numerical solutions. Moreover, the numerical scheme is validated by the results from the analytical solutions and real experimental data with the relative error below 10%. Through numerous simulations, we can conclude that the wave attenuation factors are friction coefficient c f , diffusion coefficient c D , and vegetation length. These three properties of the porous media profoundly affect K t . If just one of these coefficients increases, then K t will decrease, therefore, increasing the wave attenuation significantly. Furthermore, there are two additional interesting findings from this research. First, the proportion of K t reduction caused by the increase of c f is much higher than c D , which implies that the friction factor affects wave attenuation more than the diffusion factor. Second, the length between vegetation patches does not have a particular effect on K t , which means that the distance between vegetation is relatively insignificant compared to other parameters. All in all, this research can be used as an evaluation tool to determine coastal protection using porous media, specifically mangroves.