the Creative Commons Attribution-NonCommercial-ShareAlike 2.5 License.
the Creative Commons Attribution-NonCommercial-ShareAlike 2.5 License.
Wave-resolving Voronoi model of Rouse number for sediment entrainment equilibrium
Abstract. To integrate wave and sediment transport modeling, a computationally extensive wave-resolving Voronoi mesh-based simulation has been developed to improve upon heretofore separate sediment and spectral wave modeling. Orbital wave motion-dependent sediment transport and fine structures of the dynamic Rouse number distribution across the seabed were brought into focus. The entirely parallelized wave-resolving hydrodynamic model is demonstrated for nearshore beach waters adjacent to artificial islands in Doha Bay. The nested model was validated with tidal time series for three locations and two seasons.
-
Notice on discussion status
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
-
Preprint
(0 KB)
-
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.
Journal article(s) based on this preprint
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-1213', Anonymous Referee #1, 12 Jun 2024
-
AC1: 'Reply on RC1', Johannes Lawen, 03 Jul 2024
Thank you for the review. Given that it is a bay simulation, it seems unlikely to find an amphidromic response in such a small area.Regarding denser and longer sampling: These data are available. Full mark density merely got the Arxiv Latex compiler choking and concerned me in terms of readability of other statistical quantities and was, therefore, reduced. But I just replotted, as suggested, with higher mark density (attached for TDM2). Will then put in revision all with high mark density, perhaps separate from stats for readability. Another request was to run for longer. Fortunately, sampled data are available for about a month. Will rerun then for entire sampling period for revision.
-
AC1: 'Reply on RC1', Johannes Lawen, 03 Jul 2024
-
CC1: 'Comment on egusphere-2024-1213', Mohamed Fawal, 28 Jun 2024
Congratulations. After a decade of no one getting a permit for this area approved by the ministry of environment, you got one this month for Katara beaches with this new model
Citation: https://doi.org/10.5194/egusphere-2024-1213-CC1 - CC2: 'Comment on egusphere-2024-1213', George Salman, 02 Jul 2024
- CC3: 'Trial with Gambia River', Abdulfattah Olanrewaju, 04 Jul 2024
-
RC2: 'Comment on egusphere-2024-1213', Anonymous Referee #2, 22 Jul 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-1213/egusphere-2024-1213-RC2-supplement.pdf
-
AC2: 'Reply on RC2', Johannes Lawen, 19 Aug 2024
Thank you for the review questions.
1. The reference to Holleman et al. (2013) alongside Chan et al. (2018) is made to point out the significance of numerical diffusion. The connection to Voronoi meshes is drawn by the second citation in the bracket, i.e. Chan et al. (2018). Voronoi meshes exhibit fewer acute polygon angles, spanned by cell vertices that are acutely remote from the cell centroid, which factor into numerical diffusion. Of course, flow aligned meshes are even better to tackle numerical diffusion but are here impeded by tidal dynamics bar the mesh itself becoming dynamic. Have changed "Furthermore, Voronoi approximations exhibit a reduction in numerical diffusion vs Delaunay meshes [4, 5]" to "In terms of numerical diffusion [4], Voronoi meshes exhibit a reduction compared to Delaunay meshes [5]." in the pending Arxiv revision.
2. Absent GPU acceleration, running this simulation along 5 other hydrodynamic simulations took one week on a Ryzen 9 7950X3D. That translates to about one day if all the CPU cores are dedicated to the simulation. Unboudedness over time is not a concern for numerical diffusion as it has a nivellating effect. In fact, its dampening effect can mitigate overshooting at the expense of accuracy.
3. The algorithm has been validated with an MMS, method of manufactured solutions, which has been submitted separately for publication. The MMS was realized by oscillating the seabed to match the flow field to an analytical solution. The method can be used to validate the algorithm and mesh. Have included points mentioned in 3. and 4. below in the pending Arxiv revision.
4. The method is first order in space and time to attain high resolution meshes (Figure 6) to resolve waves while remaining efficient in terms of flops: To resolve waves, the cell size should be a log order below the part of the wave spectrum that is to be resolved. I.e. maximizing cell count and minimizing flops per cell are prioritized.
Timewise the LHS in equations such as #11 denotes forward Euler approximation. The ^{+\delta t} denotes a quantity at the subsequent time level. Past Delaunay versions of a species transport model (Lawen et al., 2013, Lawen et al., 2014 as cited in the paper), which worked in conjunction with other ocean models, offered for scalar quantities also semi-implicit matrix reordering algos. However, these attained only a tripling of time steps at the expense of flops for the reordering, rendering the net gain in terms of flops questionable.
Based on this experience, a split approach is pursued: fast explicit Voronoi algorithm (small timesteps but less flops per step) and slow implicit reference solver (published separately) for algo cross-validation in addition to validation with MMS and survey.
Yes, Voronoi schemes can be expanded to n dimensions. Interestingly, in this context, that might not improve results: for example, in the coastal case, the usual approach of resolving the vertical rather via multiple layers retains an alignment with the dominant horizontal current components and, thus, avoids numerical diffusion. That is, utilizing multiple layers achieves quasi flow alignment for the vertical. This example harkens back to question 1 above about numerical diffusion.
But perhaps expanding the Voronoi tessellation to the vertical could be an enhancement to model wave breaking and dynamic coastal meshes (4D Voronoi).5. The model is capable of back-coupling modeled quantities onto hydrodynamic properties. There are two ways to account for bed evolution: A.) Bed thickness change fluxes are modeled for representative periods, such as a neap-spring cycle, and then extrapolated for an annual seabed update. B.) To really simulate multiple years. The latter would preferably be supported by GPU acceleration.
6. Yes, did change "Numerical diffusion for flow-aligned unstructured grids" to "Numerical diffusion for flow-aligned unstructured grids with applications to estuarine modeling". Thank you.
Citation: https://doi.org/10.5194/egusphere-2024-1213-AC2
-
AC2: 'Reply on RC2', Johannes Lawen, 19 Aug 2024
Interactive discussion
Status: closed
-
RC1: 'Comment on egusphere-2024-1213', Anonymous Referee #1, 12 Jun 2024
-
AC1: 'Reply on RC1', Johannes Lawen, 03 Jul 2024
Thank you for the review. Given that it is a bay simulation, it seems unlikely to find an amphidromic response in such a small area.Regarding denser and longer sampling: These data are available. Full mark density merely got the Arxiv Latex compiler choking and concerned me in terms of readability of other statistical quantities and was, therefore, reduced. But I just replotted, as suggested, with higher mark density (attached for TDM2). Will then put in revision all with high mark density, perhaps separate from stats for readability. Another request was to run for longer. Fortunately, sampled data are available for about a month. Will rerun then for entire sampling period for revision.
-
AC1: 'Reply on RC1', Johannes Lawen, 03 Jul 2024
-
CC1: 'Comment on egusphere-2024-1213', Mohamed Fawal, 28 Jun 2024
Congratulations. After a decade of no one getting a permit for this area approved by the ministry of environment, you got one this month for Katara beaches with this new model
Citation: https://doi.org/10.5194/egusphere-2024-1213-CC1 - CC2: 'Comment on egusphere-2024-1213', George Salman, 02 Jul 2024
- CC3: 'Trial with Gambia River', Abdulfattah Olanrewaju, 04 Jul 2024
-
RC2: 'Comment on egusphere-2024-1213', Anonymous Referee #2, 22 Jul 2024
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2024/egusphere-2024-1213/egusphere-2024-1213-RC2-supplement.pdf
-
AC2: 'Reply on RC2', Johannes Lawen, 19 Aug 2024
Thank you for the review questions.
1. The reference to Holleman et al. (2013) alongside Chan et al. (2018) is made to point out the significance of numerical diffusion. The connection to Voronoi meshes is drawn by the second citation in the bracket, i.e. Chan et al. (2018). Voronoi meshes exhibit fewer acute polygon angles, spanned by cell vertices that are acutely remote from the cell centroid, which factor into numerical diffusion. Of course, flow aligned meshes are even better to tackle numerical diffusion but are here impeded by tidal dynamics bar the mesh itself becoming dynamic. Have changed "Furthermore, Voronoi approximations exhibit a reduction in numerical diffusion vs Delaunay meshes [4, 5]" to "In terms of numerical diffusion [4], Voronoi meshes exhibit a reduction compared to Delaunay meshes [5]." in the pending Arxiv revision.
2. Absent GPU acceleration, running this simulation along 5 other hydrodynamic simulations took one week on a Ryzen 9 7950X3D. That translates to about one day if all the CPU cores are dedicated to the simulation. Unboudedness over time is not a concern for numerical diffusion as it has a nivellating effect. In fact, its dampening effect can mitigate overshooting at the expense of accuracy.
3. The algorithm has been validated with an MMS, method of manufactured solutions, which has been submitted separately for publication. The MMS was realized by oscillating the seabed to match the flow field to an analytical solution. The method can be used to validate the algorithm and mesh. Have included points mentioned in 3. and 4. below in the pending Arxiv revision.
4. The method is first order in space and time to attain high resolution meshes (Figure 6) to resolve waves while remaining efficient in terms of flops: To resolve waves, the cell size should be a log order below the part of the wave spectrum that is to be resolved. I.e. maximizing cell count and minimizing flops per cell are prioritized.
Timewise the LHS in equations such as #11 denotes forward Euler approximation. The ^{+\delta t} denotes a quantity at the subsequent time level. Past Delaunay versions of a species transport model (Lawen et al., 2013, Lawen et al., 2014 as cited in the paper), which worked in conjunction with other ocean models, offered for scalar quantities also semi-implicit matrix reordering algos. However, these attained only a tripling of time steps at the expense of flops for the reordering, rendering the net gain in terms of flops questionable.
Based on this experience, a split approach is pursued: fast explicit Voronoi algorithm (small timesteps but less flops per step) and slow implicit reference solver (published separately) for algo cross-validation in addition to validation with MMS and survey.
Yes, Voronoi schemes can be expanded to n dimensions. Interestingly, in this context, that might not improve results: for example, in the coastal case, the usual approach of resolving the vertical rather via multiple layers retains an alignment with the dominant horizontal current components and, thus, avoids numerical diffusion. That is, utilizing multiple layers achieves quasi flow alignment for the vertical. This example harkens back to question 1 above about numerical diffusion.
But perhaps expanding the Voronoi tessellation to the vertical could be an enhancement to model wave breaking and dynamic coastal meshes (4D Voronoi).5. The model is capable of back-coupling modeled quantities onto hydrodynamic properties. There are two ways to account for bed evolution: A.) Bed thickness change fluxes are modeled for representative periods, such as a neap-spring cycle, and then extrapolated for an annual seabed update. B.) To really simulate multiple years. The latter would preferably be supported by GPU acceleration.
6. Yes, did change "Numerical diffusion for flow-aligned unstructured grids" to "Numerical diffusion for flow-aligned unstructured grids with applications to estuarine modeling". Thank you.
Citation: https://doi.org/10.5194/egusphere-2024-1213-AC2
-
AC2: 'Reply on RC2', Johannes Lawen, 19 Aug 2024
Peer review completion
Post-review adjustments
Journal article(s) based on this preprint
Viewed
Since the preprint corresponding to this journal article was posted outside of Copernicus Publications, the preprint-related metrics are limited to HTML views.
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 235 | 0 | 0 | 235 | 0 | 0 |
- HTML: 235
- PDF: 0
- XML: 0
- Total: 235
- BibTeX: 0
- EndNote: 0
Viewed (geographical distribution)
Since the preprint corresponding to this journal article was posted outside of Copernicus Publications, the preprint-related metrics are limited to HTML views.
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
Johannes Lawen
The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.