1
Fork 0

Nature v1

This commit is contained in:
Edgar P. Burkhart 2022-07-04 14:39:16 +02:00
parent 2422c412f7
commit c90e3b8b86
Signed by: edpibu
GPG key ID: 9833D3C5A25BD227
3 changed files with 120 additions and 87 deletions

Binary file not shown.

Binary file not shown.

View file

@ -1,9 +1,19 @@
\documentclass[a4paper, twocolumn, draft]{article} \documentclass[a4paper, twocolumn]{article}
\usepackage{polyglossia} \usepackage{authblk} \usepackage{polyglossia} \usepackage{authblk}
\usepackage[sfdefault]{inter} \usepackage[sfdefault]{inter}
\usepackage[math-style=french]{unicode-math}
\setmathfont{Fira Math}
\usepackage{graphicx} \usepackage{graphicx}
\usepackage[hmargin=2.1cm, vmargin=2.97cm]{geometry} \usepackage[hmargin=2.1cm, vmargin=2.97cm]{geometry}
\usepackage{hyperref} \usepackage{hyperref}
\usepackage{siunitx}
\sisetup{
mode=text,
reset-text-family=false,
reset-text-series=false,
reset-text-shape=false,
propagate-math-font=true,
}
\setmainlanguage{english} \setmainlanguage{english}
@ -33,21 +43,21 @@
% Displacement of blocks studies % Displacement of blocks studies
Displacement of large blocks or boulders by waves is an interesting phenomenon in the study of extreme historical Displacement of large blocks or boulders by waves is an interesting phenomenon in the study of extreme historical
coastal events. The existence of block deposits at unusual heights can be a clue to past events such as extreme storms coastal events. The existence of block deposits at unusual heights can be a clue to past events such as extreme storms
or tsunamis. For instance, \textcite{cox2018} studied coastal deposits on the coast of Ireland in relation to the or tsunamis. For instance, \textcite{cox2018} studied coastal deposits on the coast of Ireland in relation to the storms
storms from winter 2013--2014, and extracted criteria for analysing such deposits. Similarly, \textcite{shah2013} from winter 2013--2014, and extracted criteria for analysing such deposits. Similarly, \textcite{shah2013} found boulder
found boulder deposits on the mediterranean coast to be evidence of extreme storms in the Little Ice Age. deposits on the mediterranean coast to be evidence of extreme storms in the Little Ice Age.
% Need for analytical equations % Need for analytical equations
In order for those studies to be possible, analytical criterias are needed in order to ascertain the cause of the In order for those studies to be possible, analytical criterias are needed in order to ascertain the cause of the
displacement of a block. \textcite{nott1997,nott2003} proposed a set of equations that have been widely used for that displacement of a block. \textcite{nott1997,nott2003} proposed a set of equations that have been widely used for that
purpose. Those equations rely on an equilibrium relation between the lift force produced by a wave and restraining purpose. Those equations rely on an equilibrium relation between the lift force produced by a wave and restraining
forces depending on the initial setting of the block, allowing to extract a minimal flow velocity necessary for forces depending on the initial setting of the block, allowing to extract a minimal flow velocity necessary for movement
movement initiation. A parametrisation of waves depending on their source is also used to provide minimal wave heights initiation. A parametrisation of waves depending on their source is also used to provide minimal wave heights depending
depending on the type of scenario --- wave or tsunami. Those equations were later revised by \textcite{nandasena2011}, on the type of scenario --- wave or tsunami. Those equations were later revised by \textcite{nandasena2011}, as they
as they were found to be partially incorrect. A revised formulation based on the same considerations was provided. were found to be partially incorrect. A revised formulation based on the same considerations was provided.
The assumptions on which \citeauthor{nott2003, nandasena2011} are based were then critisized by \textcite{weiss2015}. The assumptions on which \citeauthor{nott2003, nandasena2011} are based were then critisized by \textcite{weiss2015}. In
In fact, according to them, the initiation of movement is not sufficient to guarantee block displacement. fact, according to them, the initiation of movement is not sufficient to guarantee block displacement.
\textcite{weiss2015} highlights the importance of the time dependency on block displacement. A method is proposed that \textcite{weiss2015} highlights the importance of the time dependency on block displacement. A method is proposed that
allows to find the wave amplitude that lead to block displacement. allows to find the wave amplitude that lead to block displacement.
@ -56,44 +66,44 @@ Whether it is \textcite{nott2003}, \textcite{nandasena2011} or \textcite{weiss20
equations suffer from a major flaw; they are all based on simplified analytical models and statistical analysis. equations suffer from a major flaw; they are all based on simplified analytical models and statistical analysis.
Unfortunately, no block displacement event seems to have been observed directly in the past. Unfortunately, no block displacement event seems to have been observed directly in the past.
In this paper, we study such an event. On February 28, 2017, a 50T concrete block was dropped by a wave on the crest of In this paper, we study such an event. On February 28, 2017, a \SI{50}{\tonne} concrete block was dropped by a wave on
the Artha breakwater. Luckily, the event was captured by a photographer, and a wave buoy located 1.2km offshore the crest of the Artha breakwater. Luckily, the event was captured by a photographer, and a wave buoy located
captured the seastate. Information from the photographer allowed to establish the approximate time at which the block \SI{1.2}{\km} offshore captured the seastate. Information from the photographer allowed to establish the approximate
displacement occured. The goal of this paper is to model the hydrodynamic conditions near the breakwater that lead to time at which the block displacement occured. The goal of this paper is to model the hydrodynamic conditions near the
the displacement of the 50T concrete block. breakwater that lead to the displacement of the \SI{50}{\tonne} concrete block.
% Modeling flow accounting for porous media % Modeling flow accounting for porous media
Several approaches can be used when modelling flow near a breakwater. Depth-averaged models can be used to study the Several approaches can be used when modelling flow near a breakwater. Depth-averaged models can be used to study the
transformation of waves on complex bottoms. Studying the hydrodynamic conditions under the surface can be achieved transformation of waves on complex bottoms. Studying the hydrodynamic conditions under the surface can be achieved using
using smoothed-particles hydrodynamics (SPH) or volume of fluid (VOF) models. SPH models rely on a Lagrangian smoothed-particles hydrodynamics (SPH) or volume of fluid (VOF) models. SPH models rely on a Lagrangian representation
representation of the fluid, while VOF models rely on an Eulerian representation. VOF models are generally more mature of the fluid, while VOF models rely on an Eulerian representation. VOF models are generally more mature for the study of
for the study of multiphase incompressible flows. multiphase incompressible flows.
In this paper, we first use a one-dimensionnal depth-averaged non-linear non-hydrostatic model to verify that the In this paper, we first use a one-dimensionnal depth-averaged non-linear non-hydrostatic model to verify that the signal
signal measured by the wave buoy can be used as an incident wave input for the determination of hydrodynamic conditions measured by the wave buoy can be used as an incident wave input for the determination of hydrodynamic conditions near
near the breakwater. For this model, we use a SWASH model \parencite{zijlema2011} already calibrated by the breakwater. For this model, we use a SWASH model \parencite{zijlema2011} already calibrated by \textcite{poncet2021}
\textcite{poncet2021} on a domain reaching 1450m offshore of the breakwater. on a domain reaching \SI{1450}{\m} offshore of the breakwater.
Then, we use a nested VOF model in two vertical dimensions that uses the output from the larger scale SWASH model as Then, we use a nested VOF model in two vertical dimensions that uses the output from the larger scale SWASH model as
initial and boundary conditions to obtain the hydrodynamic conditions on the breakwater. The models uses olaFlow initial and boundary conditions to obtain the hydrodynamic conditions on the breakwater. The models uses olaFlow
\parencite{higuera2015}, a VOF model based on volume averaged Reynolds averaged Navier-Stokes (VARANS) equations, and \parencite{higuera2015}, a VOF model based on volume averaged Reynolds averaged Navier-Stokes (VARANS) equations, and
which relies on a macroscopic representation of the porous armour of the breakwater. The model is qualitatively which relies on a macroscopic representation of the porous armour of the breakwater. The model is qualitatively
calibrated using photographs from the storm of February 28, 2017. Results from the nested models are finally compared calibrated using photographs from the storm of February 28, 2017. Results from the nested models are finally compared to
to the analytical equations provided by \textcite{nandasena2011}. the analytical equations provided by \textcite{nandasena2011}.
\section{Results} \section{Results}
\subsection{Identified wave} \subsection{Identified wave}
Preliminary work with the photographer allowed to identify the time at which the block displacement event happened. Preliminary work with the photographer allowed to identify the time at which the block displacement event happened.
Using the data from the wave buoy located 1250m offshore of the Artha breakwater, a seamingly abnormally large wave of Using the data from the wave buoy located \SI{1250}{\m} offshore of the Artha breakwater, a seamingly abnormally large
14m amplitude was identified that is supposed to have lead to the block displacement. wave of \SI{14}{\m} amplitude was identified that is supposed to have lead to the block displacement.
Initial analysis of the buoy data plotted in Figure~\ref{fig:wave} shows that the movement of the buoy follows two Initial analysis of the buoy data plotted in Figure~\ref{fig:wave} shows that the movement of the buoy follows two
orbitals that correspond to an incident wave direction. These results would indicate that the identified wave is orbitals that correspond to an incident wave direction. These results would indicate that the identified wave is
essentially an incident wave, with a minor reflected component. essentially an incident wave, with a minor reflected component.
The wavelet power spectrum displayed in Figure~\ref{fig:wavelet} highlights a primary infragravity wave in the signal, with The wavelet power spectrum displayed in Figure~\ref{fig:wavelet} highlights a primary infragravity wave in the signal,
a period of over 30s. with a period of over \SI{30}{\s}.
\begin{figure*} \begin{figure*}
\centering \centering
@ -105,10 +115,10 @@ a period of over 30s.
\subsection{Reflection analysis} \subsection{Reflection analysis}
The results from the large scale SWASH model using two configurations --- one of them being the real bathymetry, and The results from the large scale SWASH model using two configurations --- one of them being the real bathymetry, and the
the other being a simplified bathymetry without the breakwater --- are compared in Figure~\ref{fig:swash}. The results other being a simplified bathymetry without the breakwater --- are compared in Figure~\ref{fig:swash}. The results
obtained with both simulations show a maximum wave amplitude of 13.9m for the real bathymetry, and 12.1m in the case obtained with both simulations show a maximum wave amplitude of \SI{13.9}{\m} for the real bathymetry, and \SI{12.1}{\m}
where the breakwater is removed. in the case where the breakwater is removed.
The 13\% difference between those values highlights the existence of a notable amount of reflection at the buoy. The 13\% difference between those values highlights the existence of a notable amount of reflection at the buoy.
Nonetheless, the gap between the values is still fairly small and the extreme wave identified on February 28, 2017 at Nonetheless, the gap between the values is still fairly small and the extreme wave identified on February 28, 2017 at
@ -123,11 +133,11 @@ Nonetheless, the gap between the values is still fairly small and the extreme wa
\subsection{Wave transformation} \subsection{Wave transformation}
The free surface obtained with the SWASH model using raw buoy measurements as an elevation boundary condition is The free surface obtained with the SWASH model using raw buoy measurements as an elevation boundary condition is plotted
plotted in Figure~\ref{fig:swash_trans}. Those results display a strong transformation of the wave between the buoy and in Figure~\ref{fig:swash_trans}. Those results display a strong transformation of the wave between the buoy and the
the breakwater. Not only the amplitude, but also the shape of the wave are strongly impacted by the propagation over breakwater. Not only the amplitude, but also the shape of the wave are strongly impacted by the propagation over the
the domain. While the amplitude of the wave is reduced as the wave propagates shorewards, the length of the trough and domain. While the amplitude of the wave is reduced as the wave propagates shorewards, the length of the trough and the
the crest increases, with a zone reaching 400m long in front of the wave where the water level is below 0m. crest increases, with a zone reaching \SI{400}{\m} long in front of the wave where the water level is below \SI{0}{\m}.
\begin{figure*} \begin{figure*}
\centering \centering
@ -141,11 +151,11 @@ the crest increases, with a zone reaching 400m long in front of the wave where t
In an attempt to understand the identified wave, a wavelet analysis is conducted on raw buoy data as well as at In an attempt to understand the identified wave, a wavelet analysis is conducted on raw buoy data as well as at
different points along the SWASH model using the method proposed by \textcite{torrence1998}. The results are displayed different points along the SWASH model using the method proposed by \textcite{torrence1998}. The results are displayed
in Figure~\ref{fig:wavelet} and Figure~\ref{fig:wavelet_sw}. The wavelet power spectrum shows that the major component in Figure~\ref{fig:wavelet} and Figure~\ref{fig:wavelet_sw}. The wavelet power spectrum shows that the major component
in the identified wave is a high energy infragravity wave, with a period of around 60s. in the identified wave is a high energy infragravity wave, with a period of around \SI{60}{\s}.
The SWASH model seems to indicate that the observed transformation of the wave can be characterized by a transfer of The SWASH model seems to indicate that the observed transformation of the wave can be characterized by a transfer of
energy from the infragravity band to shorter waves from around 600m to 300m, and returning to the infragravity band at energy from the infragravity band to shorter waves from around \SI{600}{\m} to \SI{300}{\m}, and returning to the
200m. infragravity band at \SI{200}{\m}.
\begin{figure*} \begin{figure*}
\centering \centering
@ -162,16 +172,16 @@ energy from the infragravity band to shorter waves from around 600m to 300m, and
The two-dimensionnal olaFlow model near the breakwater allowed to compute the flow velocity near and on the breakwater The two-dimensionnal olaFlow model near the breakwater allowed to compute the flow velocity near and on the breakwater
during the passage of the identified wave. The results displayed in Figure~\ref{fig:U} show that the flow velocity during the passage of the identified wave. The results displayed in Figure~\ref{fig:U} show that the flow velocity
reaches a maximum of 14.5m/s towards the breakwater during the identified extreme wave. Although the maximum reached reaches a maximum of \SI{14.5}{\m\per\s} towards the breakwater during the identified extreme wave. Although the maximum
velocity is slightly lower than earlier shorter waves (at t=100s and t=120s, with a maximum velocity of 17.3s), the reached velocity is slightly lower than earlier shorter waves (at $t=\SI{100}{\s}$ and $t=\SI{120}{\s}$, with a maximum
flow velocity remains high for twice as long as during those earlier waves. The tail of the identified wave also velocity of \SI{17.3}{\m\per\s}), the flow velocity remains high for twice as long as during those earlier waves. The
exhibits a water level over 5m for over 40s. tail of the identified wave also exhibits a water level over \SI{5}{\m} for over \SI{40}{\s}.
\begin{figure*} \begin{figure*}
\centering \centering
\includegraphics{fig/U.pdf} \includegraphics{fig/U.pdf}
\caption{Horizontal flow velocity computed with the olaFlow model at x=-20m on the breakwater armor. The identified \caption{Horizontal flow velocity computed with the olaFlow model at $x=\SI{-20}{\m}$ on the breakwater armor. The
wave reaches this point around t=175s.}\label{fig:U} identified wave reaches this point around $t=\SI{175}{\s}$.}\label{fig:U}
\end{figure*} \end{figure*}
\section{Discussion} \section{Discussion}
@ -180,13 +190,14 @@ exhibits a water level over 5m for over 40s.
According to the criteria proposed by \textcite{dysthe2008}, rogue waves can be defined as waves with an amplitude over According to the criteria proposed by \textcite{dysthe2008}, rogue waves can be defined as waves with an amplitude over
twice the significant wave height over a given period. The identified wave fits this definition, as its amplitude is twice the significant wave height over a given period. The identified wave fits this definition, as its amplitude is
14.7m, over twice the significant wave height of 6.3m on that day. According to \textcite{dysthe2008}, rogue waves \SI{14.7}{\m}, over twice the significant wave height of \SI{6.3}{\m} on that day. According to \textcite{dysthe2008},
often occur from non-linear superposition of smaller waves. This seems to be what we observe on Figure~\ref{fig:wave}. rogue waves often occur from non-linear superposition of smaller waves. This seems to be what we observe on
Figure~\ref{fig:wave}.
The wavelet power spectrum shows that a very prominent infragravity component is present, which usually corresponds to The wavelet power spectrum shows that a very prominent infragravity component is present, which usually corresponds to
non-linear interactions of smaller waves. \textcite{dysthe2008} mentions that such waves in coastal waters are often non-linear interactions of smaller waves. \textcite{dysthe2008} mentions that such waves in coastal waters are often the
the result of refractive focusing. On February 28, 2017, the frequency of rogue waves was found to be of 1 wave per result of refractive focusing. On February 28, 2017, the frequency of rogue waves was found to be of 1 wave per 1627,
1627, which is considerably more than the excedance probability of 1 over 10\textsuperscript4 calculated by which is considerably more than the excedance probability of 1 over 10\textsuperscript4 calculated by
\textcite{dysthe2008}. Additionnal studies should be conducted to understand focusing and the formation of rogue waves \textcite{dysthe2008}. Additionnal studies should be conducted to understand focusing and the formation of rogue waves
in front of the Saint-Jean-de-Luz bay. in front of the Saint-Jean-de-Luz bay.
@ -196,20 +207,20 @@ The 13\% difference between those values highlights the existence of a notable a
Nonetheless, the gap between the values is still fairly small and the extreme wave identified on February 28, 2017 at Nonetheless, the gap between the values is still fairly small and the extreme wave identified on February 28, 2017 at
17:23:08 could still be considered as an incident wave. 17:23:08 could still be considered as an incident wave.
Unfortunately, the spectrum wave generation method used by SWASH could not reproduce simlar waves to the one observed Unfortunately, the spectrum wave generation method used by SWASH could not reproduce simlar waves to the one observed at
at the buoy. As mentionned by \textcite{dysthe2008}, such rogue waves cannot be deterministicly from the wave spectrum. the buoy. As mentionned by \textcite{dysthe2008}, such rogue waves cannot be deterministicly from the wave spectrum. For
For this reason, this study only allows us to observe the influence of reflection on short waves, while mostly ignoring this reason, this study only allows us to observe the influence of reflection on short waves, while mostly ignoring
infragravity waves. infragravity waves.
\subsection{Wave transformation} \subsection{Wave transformation}
The SWASH model yields a strongly changing wave over the domain, highlighting the highly complex composition of this The SWASH model yields a strongly changing wave over the domain, highlighting the highly complex composition of this
wave. Although the peak of the amplitude of the wave is reduced as the wave propagates, the length of the wave is wave. Although the peak of the amplitude of the wave is reduced as the wave propagates, the length of the wave is
highlighted by the results. At T+60s for instance, the water level is under 0m for 400m, and then over 0m for around highlighted by the results. At $T+\SI{60}{\s}$ for instance, the water level is under \SI{0}{\m} for \SI{400}{\m}, and
the same length, showing the main infragavity component of the studied wave. then over \SI{0}{\m} for around the same length, showing the main infragavity component of the studied wave.
The wavelet analysis conducted at several points along the domain (Figure~\ref{fig:wavelet_sw}) show that the energy of The wavelet analysis conducted at several points along the domain (Figure~\ref{fig:wavelet_sw}) show that the energy of
the studied wave (slightly before t=1500s) initially displays a strong infragravity component. Energy is then the studied wave (slightly before $t=\SI{1500}{\s}$) initially displays a strong infragravity component. Energy is then
transfered from the infragravity band towards shorter waves, and back to the infragravity band. This behavior is quite transfered from the infragravity band towards shorter waves, and back to the infragravity band. This behavior is quite
unexpected, and further investigations should be conducted to understand and validate those results. unexpected, and further investigations should be conducted to understand and validate those results.
@ -217,8 +228,9 @@ unexpected, and further investigations should be conducted to understand and val
The hydrodynamic conditions on the breakwater are the main focus of this study. Considering an initially submerged The hydrodynamic conditions on the breakwater are the main focus of this study. Considering an initially submerged
block, analytical equations proposed by \textcite{nandasena2011} yield a minimal flow velocity that would lead to block block, analytical equations proposed by \textcite{nandasena2011} yield a minimal flow velocity that would lead to block
displacement by saltation of 19.4m/s. The results from the Olaflow model yield a maximal wave velocity during the displacement by saltation of \SI{19.4}{\m\per\s} The results from the Olaflow model yield a maximal wave velocity during
displacement of the 50T concrete block of 14.5m/s. The results from the model are 25\% lower than the analytical value. the displacement of the \SI{50}{\tonne} concrete block of \SI{14.5}{\m\per\s}. The results from the model are 25\% lower
than the analytical value.
Those results tend to confirm recent research by \textcite{lodhi2020}, where it was found that the block displacement Those results tend to confirm recent research by \textcite{lodhi2020}, where it was found that the block displacement
threshold tend to overestimate the minimal flow velocity needed for block movement, although further validation of the threshold tend to overestimate the minimal flow velocity needed for block movement, although further validation of the
@ -232,28 +244,49 @@ phenomenon.
\section{Methods} \section{Methods}
\subsection{Buoy data analysis}
\subsubsection{Rogue wave identification}
Identifying rogue waves requires two main steps: computing the significant wave height, and computing the height of
individual waves. The first step is straightforward: $H_s=4\sigma$, where $\sigma$ is the standard deviation of the
surface elevation. Computing the height of individual waves is conducted using the zero-crossing method: the time domain
is split in sections where water level is strictly positive or negative, and wave size is computed according to the
maxima and minima in each zone. This method can fail to identify some waves or wrongly identify waves in case of
measurement errors or in the case where a small oscillation around 0 occurs in the middle of a larger wave. In order to
account for those issues, the signal is first fed through a low-pass filter to prevent high frequency oscillations of
over \SI{0.2}{\Hz}.
\subsubsection{Wavelet analysis}
All wavelet analysis in this study is conducted using a continuous wavelet transform over a Morlet window. The wavelet
power spectrum is normalized by the variance of the timeseries, following the method proposed by
\textcite{torrence1998}.
\subsection{SWASH models} \subsection{SWASH models}
\subsubsection{Domain} \subsubsection{Domain}
A 1750m long domain is constructed in order to study wave reflection and wave transformation over the bottom from the A \SI{1750}{\m} long domain is constructed in order to study wave reflection and wave transformation over the bottom
wave buoy to the breakwater. Bathymetry with a resolution of around 1m was used for most of the domain. The breakwater from the wave buoy to the breakwater. Bathymetry with a resolution of around \SI{1}{\m} was used for most of the domain.
model used in the study is taken from \textcite{poncet2021}. A smoothed section is created and considered as a porous The breakwater model used in the study is taken from \textcite{poncet2021}. A smoothed section is created and considered
media in the model. as a porous media in the model.
A second domain is constructed for reflection analysis. The second model is the same as the first, excepted that the A second domain is constructed for reflection analysis. The second model is the same as the first, excepted that the
breakwater is replaced by a smooth slope in order to remove the reflection generated by the structure. breakwater is replaced by a smooth slope in order to remove the reflection generated by the structure.
The reflection analysis is conducted over 4h in order to generate a fair range of conditions. The wave transformation The reflection analysis is conducted over \SI{4}{\hour} in order to generate a fair range of conditions. The wave
study was conducted over a 1h timeframe in order to allow the model to reach steady-state before the studied wave was transformation study was conducted over a \SI{1}{\hour} timeframe in order to allow the model to reach steady-state
generated. before the studied wave was generated.
\subsubsection{Model} \subsubsection{Model}
A non-linear non-hydrostatic shallow water model (SWASH, \cite{zijlema2011}) is used to model wave reflection and A non-linear non-hydrostatic shallow water model (SWASH, \cite{zijlema2011}) is used to model wave reflection and
transformation on the studied domain. The study is conducted using a layered one-dimensional model, that allows to transformation on the studied domain. The study is conducted using a layered one-dimensional model, that allows to
consider porous media in the domain. consider porous media in the domain.
The reflection analysis was conducted with 2 layers as to prevent model instability in overtopping conditions. The The reflection analysis was conducted with 2 layers as to prevent model instability in overtopping conditions. The study
study of wave transformation and the generation of boundary conditions for the Olaflow model is done with 4 layers. of wave transformation and the generation of boundary conditions for the Olaflow model is done with 4 layers.
\subsubsection{Porosity} \subsubsection{Porosity}
@ -264,37 +297,37 @@ calibrated in \textcite{poncet2021}.
Two different sets of boundary conditions were used for both studies. In all cases, a sponge layer was added to the Two different sets of boundary conditions were used for both studies. In all cases, a sponge layer was added to the
shorewards boundary to prevent wave reflection on the boundary. In the reflection analysis, offshore conditions were shorewards boundary to prevent wave reflection on the boundary. In the reflection analysis, offshore conditions were
generated using the wave spectrum extracted from buoy data during the storm. The raw vertical surface elevation generated using the wave spectrum extracted from buoy data during the storm. The raw vertical surface elevation measured
measured by the wave buoy was used in a second part. by the wave buoy was used in a second part.
\subsection{Olaflow model} \subsection{Olaflow model}
\subsubsection{Domain}
A 150m long domain is built in order to obtain the hydrodynamic conditions on the Artha breakwater during the passage
of the identified extreme wave. The bathymetry with 50cm resolution from \textcite{poncet2021} is used. The domain
extends 30m up in order to be able to capture the largest waves hitting the breakwater. Measurements are extracted 20m
shorewards from the breakwater crest. The domain is displayed in Figure~\ref{fig:of}.
A mesh in two-vertical dimensions with 20cm resolution was generated using the interpolated bathymetry. As with the
SWASH model, the porous armour was considered at a macroscopic scale.
\begin{figure*} \begin{figure*}
\centering \centering
\includegraphics{fig/aw_t0.pdf} \includegraphics{fig/aw_t0.pdf}
\caption{Domain studied with Olaflow. Initial configuration.}\label{fig:of} \caption{Domain studied with Olaflow. Initial configuration.}\label{fig:of}
\end{figure*} \end{figure*}
\subsubsection{Domain}
A \SI{150}{\m} long domain is built in order to obtain the hydrodynamic conditions on the Artha breakwater during the
passage of the identified extreme wave. The bathymetry with \SI{50}{\cm} resolution from \textcite{poncet2021} is used.
The domain extends \SI{30}{\m} up in order to be able to capture the largest waves hitting the breakwater. Measurements
are extracted \SI{20}{\m} shorewards from the breakwater crest. The domain is displayed in Figure~\ref{fig:of}.
A mesh in two-vertical dimensions with \SI{20}{\cm} resolution was generated using the interpolated bathymetry. As with
the SWASH model, the porous armour was considered at a macroscopic scale.
\subsubsection{Model} \subsubsection{Model}
A volume-of-fluid (VOF) model in two-vertical dimensions based on volume-averaged Reynolds-averaged Navier-Stokes A volume-of-fluid (VOF) model in two-vertical dimensions based on volume-averaged Reynolds-averaged Navier-Stokes
(VARANS) equations is used (olaFlow, \cite{higuera2015}). The model was initially setup using generic values for (VARANS) equations is used (olaFlow, \cite{higuera2015}). The model was initially setup using generic values for porous
porous breakwater studies. A sensibility study conducted on the porosity parameters found a minor influence of these breakwater studies. A sensibility study conducted on the porosity parameters found a minor influence of these values on
values on the final results. the final results.
The k-ω SST turbulence model was used, as it produced much more realistic results than the default k-ε model, The k-ω SST turbulence model was used, as it produced much more realistic results than the default k-ε model, especially
especially compared to the photographs from the storm of February 28, 2017. The k-ε model yielded very high viscosity compared to the photographs from the storm of February 28, 2017. The k-ε model yielded very high viscosity and thus
and thus strong dissipation in the entire domain, preventing an accurate wave breaking representation. strong dissipation in the entire domain, preventing an accurate wave breaking representation.
\subsubsection{Boundary conditions} \subsubsection{Boundary conditions}