diff --git a/dfastmi/messages.UK.ini b/dfastmi/messages.UK.ini index 69e16b74..17809c9f 100644 --- a/dfastmi/messages.UK.ini +++ b/dfastmi/messages.UK.ini @@ -7,13 +7,13 @@ The option 'reduce_output' is active. D-FAST Morphological Impact implements an algorithm to estimate the local morphological effects of a local intervention (i.e. an adjustment to the river). The conceptual framework was originally introduced in - "RWS-WD memo WAQUA vuistregel 20-10-08" + "RWS-WD memo WAQUA vuistregel (Sieben, 2010)" but it has been extended and improved over the years. Check the user manual for the details of the currently implemented algorithm. -It is based on an estimation of the equilibrium bed level changes in the main -channel that would occur eventually when river maintenance would not be -adjusted. +It is based on an estimation of the equilibrium bed level changes in the +main channel that would occur eventually when river maintenance would not +be adjusted. The effect is expressed as: @@ -21,14 +21,13 @@ The effect is expressed as: maximum bed level change [m] without dredging minimum bed level change [m] without dredging -By means of these estimates bottlenecks can be identified. The results are not -suitable for direct estimation of the impact on the maintenance of the +By means of these estimates bottlenecks can be identified. The results are +not suitable for direct estimation of the impact on the maintenance of the navigation channel! -The combination of the total equilibrium sedimentation volume and the yearly -sediment load of the river determines the period over which the equilibrium -can be reached. - +The combination of the total equilibrium sedimentation volume and the +yearly sediment load of the river determines the period over which the +equilibrium can be reached. This is version {version}. diff --git a/docs/chapters/application.tex b/docs/chapters/application.tex index 1b003c44..c2e32149 100644 --- a/docs/chapters/application.tex +++ b/docs/chapters/application.tex @@ -7,6 +7,17 @@ \chapter{Steps in the analysis}\label{Chp:steps} Verify that it is appropriate to use \dfmi instrument; see \autoref{Chp:Guidance}. Determine the branch and reach on which it is located; use the \dfmi GUI for this\footnote{The reaches are specified in the GUI both by a descriptive name and by an approximate indication of the river chainage. Follow item \ref{reach_bnd} of \autoref{Sec:Limitations} if an intervention is located across or near a branch/reach boundary.}. +\dfmi is currently configured to support the evaluation of interventions in the following reaches: +\begin{itemize} +\item Bovenrijn (Rkm 859-867) +\item Waal (Rkm 868-951) +\item Pannerdensch-Kanaal (Rkm 868-879) +\item Nederrijn (Rkm 880-922) +\item Lek (Rkm 923-989) +\item IJssel (Rkm 880-1000) +\item Merwede (Rkm 951-980) +\item Meuse (Rkm 16-227) +\end{itemize} Determine the threshold discharge $Q_\text{thr}$ (at Lobith/Borgharen) at which the intervention (indirectly) starts to influence the flow pattern in the main channel. The threshold discharge is critical for determining the fraction of the year that the intervention influences the flow, and hence the duration over which the bed level difference (compared to the reference situation) can develop. As a result, this value is critical for determining the total volume of sedimentation (or erosion) that can be expected after one year. @@ -30,21 +41,17 @@ \chapter{Steps in the analysis}\label{Chp:steps} These flow conditions have been pre-configured for the latest \dflowfm schematizations. It is not necessary to run simulations for conditions at which the intervention doesn't influence the flow patterns (i.e.~discharges associated with stagnant conditions due to closure of barriers, and discharges below the intervention specific $Q_\text{thr}$). At this stage, \dfmi already reports the impacted length (`aanzandingslengte'), which is the distance over which a bed level change can built up during the period that the discharge of the river is above the threshold discharge. +\dfmi is not suitable for interventions that only (start to) have a noticeable effect at (or above) the highest discharge. \item Perform for each condition the hydrodynamic simulations for both the reference situation and the situation with intervention. Verify that \begin{itemize} \item the \dflowfm results are stable -\item the intervention is properly represented on the mesh used (check a.o. proper alignment of groynes and levees, channel shape and bed roughness) +\item the intervention is properly represented on the mesh used (check a.o.~proper alignment of groynes and levees, channel shape and bed roughness) \item all simulations use the same base mesh (changes in dry areas may result in slight differences in the mesh effectively used) \item there is a visible difference in the velocities in the main channel between the simulations with and without intervention \end{itemize} -For steady state conditions, it is preferred to use the mean flow conditions over a certain period to suppress any instabilities and fluctuations in the instantaneous flow conditions. -These results can be obtained by using the Fourier option of \dflowfm. -For this purpose, a standardized Fourier input file \file{fourier\_last\_s.fou} is included that contains the following configuration: -\vspace{\baselineskip} -\verbfilenobox[\scriptsize]{figures/fourier_last_s.fou} -\vspace{\baselineskip} +See \autoref{Sec:SteadyState} for recommendations regarding steady-state results. \item Run \dfastmi to compute for each grid point in the main channel, the following three variables @@ -72,3 +79,14 @@ \chapter{Steps in the analysis}\label{Chp:steps} If there are disjunct sedimentation areas, the accumulation should be carried out for each of the areas individually. If areas are shorter than the impacted length, then the total equilibrium impact can be reached within one year. \end{enumerate} + +\section{How to get steady-state results?}\label{Sec:SteadyState} + +The dynamic solver of \dflowfm is used to obtain steady-state results by providing a constant forcing over a suitably long simulation period. +However, even with a constant forcing some fluctuations may remain in the computed flow fields, e.g.~due to natural formation and shedding of eddies, and sensitivity of numerical formulations related to for instance drying-flooding. +Therefore, it is preferred to use, for steady-state conditions, the mean flow conditions over a certain period to suppress any instabilities and fluctuations in the instantaneous flow conditions. +These results can be obtained by using the Fourier option of \dflowfm. +For this purpose, a standardized Fourier input file \file{fourier\_last\_s.fou} is included that contains the following configuration: +\vspace{\baselineskip} +\verbfilenobox[\scriptsize]{figures/fourier_last_s.fou} +\vspace{\baselineskip} diff --git a/docs/chapters/example.tex b/docs/chapters/example.tex index 654d20d2..8de845a4 100644 --- a/docs/chapters/example.tex +++ b/docs/chapters/example.tex @@ -66,18 +66,20 @@ \section{Example 1: secondary channel along the Nederrijn} \verbfilenobox[\scriptsize]{../examples_references/01 - Palmerswaard/output/report.txt} -\autoref{Palmers_mor} shows the results of a reference morphology simulation using Delft3D 4 (for details, see \citet{GiriJagers2022}). -The figure shows the morphological evolution of the Delft3D simulation. -The overall \dfastmi results compare well with the long term (12yr) Delft3D simulation although the asymmetry downstream of the main sedimentation patch differs. +\autoref{Palmers_mor} shows the results side by side of a reference morphology simulation using Delft3D 4 (for details, see \citet{GiriJagers2022}) and the \dfmi version 3 analysis. +The overall \dfmi results compare well with the long term (\SI{12}{\year}) Delft3D simulation although the asymmetry downstream of the main sedimentation patch differs. \begin{figure}[H] -\center -\includegraphics[width=\textwidth]{figures/Palmers_mor.png} -\caption{The bed-level difference between the variant and the reference after 1, 5 and 12 years of morphological simulation.} +\includegraphics[width=\columnwidth/2]{figures/Palmerswaard_delft3d.png} +\includegraphics[width=\columnwidth/2]{figures/Palmerswaard_dfastmi.png} +%\fbox{\parbox{\columnwidth}{\todo{Results to be plotted side-by-side.}}} +\caption{Long-term morphological impact as obtained from a 12 year morphological simulation (left) and as obtained from \dfmi version 3 (right). +Visualized using QUICKPLOT.} \label{Palmers_mor} \end{figure} + \section{Example 2: secondary channel along the Pannerdensch Kanaal} For this second example, we follow the same approach as example 1: we compare the results of \dfastmi with the results of a morphological simulation using Delft3D-FLOW. @@ -125,13 +127,14 @@ \section{Example 2: secondary channel along the Pannerdensch Kanaal} \verbfilenobox[\scriptsize]{../examples_references/02 - Pannerdensch Kanaal/output/report.txt} -\autoref{Pannerden_mor} shows the results of a reference morphology simulation using Delft3D 4 (for details, see \citet{GiriJagers2022}). -The figure shows the morphological evolution of the Delft3D simulation. -The overall \dfastmi results compare well with the long term (15yr) Delft3D simulation although \dfmi suggests sedimentation in the centre of the channel downstream of the main sedimentation area whereas the morphological simulation doesn't show that behaviour. +\autoref{Pannerden_mor} shows the results side by side of a reference morphology simulation using Delft3D 4 (for details, see \citet{GiriJagers2022}) and the \dfmi version 3 analysis. +The overall \dfmi results compare well with the long term (\SI{15}{\year}) Delft3D simulation although \dfmi suggests sedimentation in the centre of the channel downstream of the main sedimentation area whereas the morphological simulation doesn't show that behaviour. \begin{figure}[H] -\center -\includegraphics[width=\textwidth]{figures/Pannerden_mor.png} -\caption{The bed-level difference between the variant and the reference after 1 and 15 years of morphological simulation.} +\includegraphics[width=\columnwidth/2]{figures/Pannerden_delft3d.png} +\includegraphics[width=\columnwidth/2]{figures/Pannerden_dfastmi.png} +%\fbox{\parbox{\columnwidth}{\todo{Results to be plotted side-by-side.}}} +\caption{Long-term morphological impact as obtained from a 12 year morphological simulation (left) and as obtained from \dfmi version 3 (right). +Visualized using QUICKPLOT.} \label{Pannerden_mor} \end{figure} \ No newline at end of file diff --git a/docs/chapters/file_formats.tex b/docs/chapters/file_formats.tex index 567e7da3..53449ddd 100644 --- a/docs/chapters/file_formats.tex +++ b/docs/chapters/file_formats.tex @@ -142,13 +142,13 @@ \subsubsection*{Example} D-FAST Morphological Impact implements an algorithm to estimate the local morphological effects of a local intervention (i.e. an adjustment to the river). The conceptual framework was originally introduced in - "RWS-WD memo WAQUA vuistregel 20-10-08" + "RWS-WD memo WAQUA vuistregel (Sieben, 2010)" but it has been extended and improved over the years. Check the user manual for the details of the currently implemented algorithm. - It is based on an estimation of the equilibrium bed level changes in the main - channel that would occur eventually when river maintenance would not be - adjusted. + It is based on an estimation of the equilibrium bed level changes in the + main channel that would occur eventually when river maintenance would not + be adjusted. The effect is expressed as: @@ -156,14 +156,13 @@ \subsubsection*{Example} maximum bed level change [m] without dredging minimum bed level change [m] without dredging - By means of these estimates bottlenecks can be identified. The results are not - suitable for direct estimation of the impact on the maintenance of the + By means of these estimates bottlenecks can be identified. The results are + not suitable for direct estimation of the impact on the maintenance of the navigation channel! - The combination of the total equilibrium sedimentation volume and the yearly - sediment load of the river determines the period over which the equilibrium - can be reached. - + The combination of the total equilibrium sedimentation volume and the + yearly sediment load of the river determines the period over which the + equilibrium can be reached. This is version {version}. diff --git a/docs/chapters/guidance.tex b/docs/chapters/guidance.tex index 94ac46a5..58e2bf86 100644 --- a/docs/chapters/guidance.tex +++ b/docs/chapters/guidance.tex @@ -57,6 +57,7 @@ \section{Assumptions and limitations}\label{Sec:Limitations} \item Interventions should be properly resolved on the \dflowfm mesh. The discharges used by \dfmi should give a balanced representation of the influence of the intervention on the flow patterns. +\dfmi is not suitable for interventions that only (start to) have a noticeable effect at (or above) the highest discharge. \item \dfastmi is not yet suited for tidally influenced areas. See \autoref{Sec:Tides}. diff --git a/docs/chapters/intro.tex b/docs/chapters/intro.tex index a7c6f7a7..e69d0fcb 100644 --- a/docs/chapters/intro.tex +++ b/docs/chapters/intro.tex @@ -1,7 +1,7 @@ \chapter{Introduction} This manual describes \dfastmi version 3 which provides a rapid, first estimate of the bed level changes to be anticipated in the main channel due to the implementation of local river adjustments outside the main channel (so-called interventions). -The program is a successor to the WAQMORF and \dfastmi version 2 programs that implemented the rule of thumb developed in the context of Rijkswaterstaat programme Stroomlijn by \citep{Sieben2008}. +The program is a successor to the WAQMORF and \dfastmi version 2 programs that implemented the rule of thumb developed in the context of Rijkswaterstaat programme Stroomlijn by \citet{Sieben2010}. This new version follows the same conceptual approach, but uses a fixed set of flow conditions independent of the intervention instead of three intervention-dependent flow conditions. It assumes a seasonal discharge variation that can be represented by means of series of flow conditions. Bed level changes outside the immediate vicinity of influence of the intervention are ignored in this analysis. @@ -56,9 +56,9 @@ \chapter{Introduction} \begin{Verbatim}[frame=single, framesep=5pt] D-FAST Morphological Impact implements an algorithm to estimate the local -morphological effects of a local intervention (i.e. an adjustment to the river). -The conceptual framework was originally introduced in - "RWS-WD memo WAQUA vuistregel 20-10-08" +morphological effects of a local intervention (i.e. an adjustment to the +river). The conceptual framework was originally introduced in + "RWS-WD memo WAQUA vuistregel (Sieben, 2010)" but it has been extended and improved over the years. Check the user manual for the details of the currently implemented algorithm. @@ -66,11 +66,11 @@ \chapter{Introduction} main channel that would occur eventually when river maintenance would not be adjusted. -The effect is expressed in [m] as: +The effect is expressed as: - year-averaged bed level change without dredging - maximum bed level change without dredging - minimum bed level change without dredging + year-averaged bed level change [m] without dredging + maximum bed level change [m] without dredging + minimum bed level change [m] without dredging By means of these estimates bottlenecks can be identified. The results are not suitable for direct estimation of the impact on the maintenance of the diff --git a/docs/dfastmi.bib b/docs/dfastmi.bib index a33a6223..cbe02653 100644 --- a/docs/dfastmi.bib +++ b/docs/dfastmi.bib @@ -126,11 +126,11 @@ @TechReport{RIZA2005 Note = {In Dutch}, } -@TechReport{Sieben2008, +@TechReport{Sieben2010, Title = {Methodiek inschatting morfologische effecten in het zomerbed door lokale rivieringrepen}, Author = {A. Sieben}, Institution = {Waterdienst, Rijkswaterstaat, The Netherlands}, - Year = {2008-2011}, + Year = {2010}, Type = {Memo}, Note = {In Dutch}, url = {http://simona.deltares.nl/release/doc/techdoc/waqmorf/Memo_waqua_vuistregel_update_dec2011.pdf} diff --git a/docs/dfastmi_usermanual.tex b/docs/dfastmi_usermanual.tex index ac59ae9b..1c8587b7 100644 --- a/docs/dfastmi_usermanual.tex +++ b/docs/dfastmi_usermanual.tex @@ -9,6 +9,7 @@ \newcommand{\dfastbe}{\textrm{D-FAST~Bank~Erosion}\xspace} \newcommand{\dfmi}{\textrm{D-FAST~MI}\xspace} \newcommand{\dflowfm}{\textrm{D-Flow~FM}\xspace} +\DeclareSIUnit\year{\text{yr}} \hypersetup { @@ -54,6 +55,11 @@ \appendix \include{chapters/file_formats} \include{chapters/sim2ugrid} + +\chapter{Background information} +The following pages reproduce a memo by Arjan Sieben on the choice of the discharges and bed celerities for the Rhine branches and the Meuse river. +\includepdf[pages=-, offset=72 -70, frame=true, scale=0.9]{figures/afvoerblokken plus uitbreiding voortplantingssnelheden.pdf} + %\markdownInput[shiftHeadings=1]{techref.md} \pagestyle{empty} diff --git a/docs/dfastmi_validation.tex b/docs/dfastmi_validation.tex index bff69846..8bea7993 100644 --- a/docs/dfastmi_validation.tex +++ b/docs/dfastmi_validation.tex @@ -41,7 +41,7 @@ \include{chapters/validation_conclusion} \nonumchapter{References} -\bibliography{dfast} +\bibliography{dfastmi} \pagestyle{empty} \cleardoublepage diff --git a/docs/figures/afvoerblokken plus uitbreiding voortplantingssnelheden.pdf b/docs/figures/afvoerblokken plus uitbreiding voortplantingssnelheden.pdf new file mode 100644 index 00000000..7c81be72 Binary files /dev/null and b/docs/figures/afvoerblokken plus uitbreiding voortplantingssnelheden.pdf differ diff --git a/examples_references/01 - Palmerswaard/output/report.txt b/examples_references/01 - Palmerswaard/output/report.txt index dadd184c..6cc470de 100644 --- a/examples_references/01 - Palmerswaard/output/report.txt +++ b/examples_references/01 - Palmerswaard/output/report.txt @@ -1,13 +1,13 @@ D-FAST Morphological Impact implements an algorithm to estimate the local morphological effects of a local intervention (i.e. an adjustment to the river). The conceptual framework was originally introduced in - "RWS-WD memo WAQUA vuistregel 20-10-08" + "RWS-WD memo WAQUA vuistregel (Sieben, 2010)" but it has been extended and improved over the years. Check the user manual for the details of the currently implemented algorithm. -It is based on an estimation of the equilibrium bed level changes in the main -channel that would occur eventually when river maintenance would not be -adjusted. +It is based on an estimation of the equilibrium bed level changes in the +main channel that would occur eventually when river maintenance would not +be adjusted. The effect is expressed as: @@ -15,14 +15,13 @@ The effect is expressed as: maximum bed level change [m] without dredging minimum bed level change [m] without dredging -By means of these estimates bottlenecks can be identified. The results are not -suitable for direct estimation of the impact on the maintenance of the +By means of these estimates bottlenecks can be identified. The results are +not suitable for direct estimation of the impact on the maintenance of the navigation channel! -The combination of the total equilibrium sedimentation volume and the yearly -sediment load of the river determines the period over which the equilibrium -can be reached. - +The combination of the total equilibrium sedimentation volume and the +yearly sediment load of the river determines the period over which the +equilibrium can be reached. This is version 3.0.0. diff --git a/examples_references/02 - Pannerdensch Kanaal/output/report.txt b/examples_references/02 - Pannerdensch Kanaal/output/report.txt index 1e724a81..b5a8ea4f 100644 --- a/examples_references/02 - Pannerdensch Kanaal/output/report.txt +++ b/examples_references/02 - Pannerdensch Kanaal/output/report.txt @@ -1,13 +1,13 @@ D-FAST Morphological Impact implements an algorithm to estimate the local morphological effects of a local intervention (i.e. an adjustment to the river). The conceptual framework was originally introduced in - "RWS-WD memo WAQUA vuistregel 20-10-08" + "RWS-WD memo WAQUA vuistregel (Sieben, 2010)" but it has been extended and improved over the years. Check the user manual for the details of the currently implemented algorithm. -It is based on an estimation of the equilibrium bed level changes in the main -channel that would occur eventually when river maintenance would not be -adjusted. +It is based on an estimation of the equilibrium bed level changes in the +main channel that would occur eventually when river maintenance would not +be adjusted. The effect is expressed as: @@ -15,14 +15,13 @@ The effect is expressed as: maximum bed level change [m] without dredging minimum bed level change [m] without dredging -By means of these estimates bottlenecks can be identified. The results are not -suitable for direct estimation of the impact on the maintenance of the +By means of these estimates bottlenecks can be identified. The results are +not suitable for direct estimation of the impact on the maintenance of the navigation channel! -The combination of the total equilibrium sedimentation volume and the yearly -sediment load of the river determines the period over which the equilibrium -can be reached. - +The combination of the total equilibrium sedimentation volume and the +yearly sediment load of the river determines the period over which the +equilibrium can be reached. This is version 3.0.0. diff --git a/tests/c01 - GendtseWaardNevengeul/ref_Qmin_Q4000_rkm/report.txt b/tests/c01 - GendtseWaardNevengeul/ref_Qmin_Q4000_rkm/report.txt index ff681ea9..42cc47f6 100644 --- a/tests/c01 - GendtseWaardNevengeul/ref_Qmin_Q4000_rkm/report.txt +++ b/tests/c01 - GendtseWaardNevengeul/ref_Qmin_Q4000_rkm/report.txt @@ -1,13 +1,13 @@ D-FAST Morphological Impact implements an algorithm to estimate the local morphological effects of a local intervention (i.e. an adjustment to the river). The conceptual framework was originally introduced in - "RWS-WD memo WAQUA vuistregel 20-10-08" + "RWS-WD memo WAQUA vuistregel (Sieben, 2010)" but it has been extended and improved over the years. Check the user manual for the details of the currently implemented algorithm. -It is based on an estimation of the equilibrium bed level changes in the main -channel that would occur eventually when river maintenance would not be -adjusted. +It is based on an estimation of the equilibrium bed level changes in the +main channel that would occur eventually when river maintenance would not +be adjusted. The effect is expressed as: @@ -15,14 +15,13 @@ The effect is expressed as: maximum bed level change [m] without dredging minimum bed level change [m] without dredging -By means of these estimates bottlenecks can be identified. The results are not -suitable for direct estimation of the impact on the maintenance of the +By means of these estimates bottlenecks can be identified. The results are +not suitable for direct estimation of the impact on the maintenance of the navigation channel! -The combination of the total equilibrium sedimentation volume and the yearly -sediment load of the river determines the period over which the equilibrium -can be reached. - +The combination of the total equilibrium sedimentation volume and the +yearly sediment load of the river determines the period over which the +equilibrium can be reached. This is version 3.0.0. diff --git a/tests/c01 - GendtseWaardNevengeul/ref_Qmin_Q4000_special_backward_case/report.txt b/tests/c01 - GendtseWaardNevengeul/ref_Qmin_Q4000_special_backward_case/report.txt index ff681ea9..42cc47f6 100644 --- a/tests/c01 - GendtseWaardNevengeul/ref_Qmin_Q4000_special_backward_case/report.txt +++ b/tests/c01 - GendtseWaardNevengeul/ref_Qmin_Q4000_special_backward_case/report.txt @@ -1,13 +1,13 @@ D-FAST Morphological Impact implements an algorithm to estimate the local morphological effects of a local intervention (i.e. an adjustment to the river). The conceptual framework was originally introduced in - "RWS-WD memo WAQUA vuistregel 20-10-08" + "RWS-WD memo WAQUA vuistregel (Sieben, 2010)" but it has been extended and improved over the years. Check the user manual for the details of the currently implemented algorithm. -It is based on an estimation of the equilibrium bed level changes in the main -channel that would occur eventually when river maintenance would not be -adjusted. +It is based on an estimation of the equilibrium bed level changes in the +main channel that would occur eventually when river maintenance would not +be adjusted. The effect is expressed as: @@ -15,14 +15,13 @@ The effect is expressed as: maximum bed level change [m] without dredging minimum bed level change [m] without dredging -By means of these estimates bottlenecks can be identified. The results are not -suitable for direct estimation of the impact on the maintenance of the +By means of these estimates bottlenecks can be identified. The results are +not suitable for direct estimation of the impact on the maintenance of the navigation channel! -The combination of the total equilibrium sedimentation volume and the yearly -sediment load of the river determines the period over which the equilibrium -can be reached. - +The combination of the total equilibrium sedimentation volume and the +yearly sediment load of the river determines the period over which the +equilibrium can be reached. This is version 3.0.0.