Merge branch 'report'

report/chapters/swash.tex

@@ -44,3 +44,85 @@ made without success, as the model does not seem to be able to compute
 overtopping. Changing the boundary condition at the breakwater does not fix the
 issue, and the model is not able to run with water on both sides of the
 breakwater as the initial condition.
+
+\paragraph{SWASH overtopping} \cite{suzuki2011applicability,zhang2020numerical} It seems like computing wave overtopping
+should be possible using the SWASH model.
+
+Implemented model with longer domain, with water behind the breakwater; currently, model crashes when overtopping
+happens.
+
+The model DOES NOT crash if layers are disabled (depth-averaged model). No overtopping seems to appear in the results.
+
+Resulting reflection coefficient in \autoref{fig:swash_nolay}: very high reflection coefficient, probably inaccurate.
+
+\begin{figure}
+	\centering
+	\includegraphics{R_singlelayer.png}
+	\caption{Reflection coefficient computed with Swash without layers (method 1, probable noise-induced
+	bias.}\label{fig:swash_nolay}
+\end{figure}
+
+Model does also work with 2 layers; overtopping does appear \autoref{fig:swash_2lay}.
+
+\begin{figure}
+	\centering
+	\includegraphics{R_2lay.png}
+	\caption{Reflection coefficient computed with Swash 2 layers (method 1, probable noise-induced
+	bias.}\label{fig:swash_2lay}
+\end{figure}
+
+Model crashes with 3 or more layers.
+
+\subsection{Model 3}
+Model: \SI{1450}{\m} offshore, \SI{300}{\m} shorewards. \SI{250}{\m} sponge layer. \SI{4.5}{\m} water level. 2 layers.
+\autoref{fig:bathy}
+
+\begin{figure}
+	\centering
+	\includegraphics{bathy_b.pdf}
+	\includegraphics{bathy_nb.pdf}
+	\caption{Bathymetry.}\label{fig:bathy}
+\end{figure}
+
+Results \autoref{fig:res45}.
+\begin{figure}
+	\centering
+	\includegraphics{R_2lay_45.pdf}
+	\caption{Results (1: breakwater; 2: no breakwater).}\label{fig:res45}
+\end{figure}
+
+Comparison with reflex3s (array method adapted from PA Poncet) \autoref{fig:reflex3s2lay}.
+\begin{figure}
+	\centering
+	\includegraphics{reflex3s_2lay.pdf}
+	\caption{Reflection coefficient with array method.}\label{fig:reflex3s2lay}
+\end{figure}
+
+\subsection{Model with measured time-series}
+Using raw buoy data. Model crashes with 2 layers. Works when layers are disabled. \textbf{No overtopping!}
+\autoref{fig:rests}. \SI{1}{\hour} computation time.
+
+\begin{figure}
+	\centering
+	\includegraphics{R_ts.pdf}
+	\caption{Results with real timeseries, no layers.}\label{fig:rests}
+\end{figure}
+
+Model runs with 2 layers after adjusting some parameters. \autoref{fig:rests2lay}. Overtopping is computed.
+
+\begin{figure}
+	\centering
+	\includegraphics{R_ts_2lay.pdf}
+	\caption{Results with real timeseries, 2 layers.}\label{fig:rests2lay}
+\end{figure}
+
+Model runs with 4 layers. \autoref{fig:rests4lay}.
+
+\begin{figure}
+	\centering
+	\includegraphics{R_ts_4lay.pdf}
+	\caption{Results with real timeseries, 4 layers.}\label{fig:rests4lay}
+\end{figure}
+
+%\subsection{2D Model}
+%Working on 2D model which might work with overtopping.

report/library.bib

@@ -1173,3 +1173,22 @@
   publisher={Elsevier}
 }
 
+@inproceedings{suzuki2011applicability,
+  title={The applicability of SWASH model for wave transformation and wave overtopping: A case study for the Flemish coast},
+  author={Suzuki, Tomohiro and Verwaest, Toon and Hassan, Wael and Veale, William and Reyns, Johan and Trouw, Koen and Troch, Peter and Zijlema, M},
+  booktitle={Proc. 5th Int. Conf. Advanced Computational Methods Engineering (ACOMEN 2011), Li{\`e}ge, Belgium},
+  pages={14--17},
+  year={2011}
+}
+
+@article{zhang2020numerical,
+  title={Numerical simulation of wave overtopping on breakwater with an armor layer of accropode using SWASH model},
+  author={Zhang, Na and Zhang, Qinghe and Wang, Keh-Han and Zou, Guoliang and Jiang, Xuelian and Yang, Aiwu and Li, Yan},
+  journal={Water},
+  volume={12},
+  number={2},
+  pages={386},
+  year={2020},
+  publisher={Multidisciplinary Digital Publishing Institute}
+}
+