An automatic mesh generator for coupled 1D/2D hydrodynamic models

Kang, Younghun; Kubatko, Ethan J.

doi:https://doi.org/10.5194/egusphere-2023-1434

Preprints

https://doi.org/10.5194/egusphere-2023-1434

Preprints

29 Sep 2023

| 29 Sep 2023

An automatic mesh generator for coupled 1D/2D hydrodynamic models

Younghun Kang and Ethan J. Kubatko

Abstract. Two-dimensional (2D), depth-averaged shallow water equation (SWE) models are routinely used to simulate flooding in coastal areas — areas that often include vast networks of channels and flood-control topographic features/structures, such as barrier islands and levees. Adequately resolving these features within the confines of a 2D model can be computationally expensive, which has led to coupling 2D simulation tools to less-expensive, one-dimensional (1D) models. Under certain 1D/2D coupling approaches, this introduces internal constraints that must be considered in the generation of the 2D computational mesh used. In this paper, we further develop an existing automatic unstructured mesh generation tool for SWE models, ADMESH+, to sequentially (i) identify 1D constraints from the raw input data used in the mesh generation process, namely, the digital elevation model (DEM) and land/water delineation data, (ii) distribute grid points along these internal constraints according to feature curvature and user-prescribed minimum grid spacing, and (iii) integrate these internal constraints into the 2D mesh-size-function and mesh-generation processes. The developed techniques, which include a novel approach for determining the so-called medial axis of a polygon, are described in detail and demonstrated on three test cases, including two inland watersheds with vast networks of channels and a complex estuarian system on the Texas, US coast.

Received: 30 Jun 2023 – Discussion started: 29 Sep 2023

Download & links

Preprint (PDF, 56038 KB)

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 (56038 KB)

Download & links

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

23 Feb 2024

An automatic mesh generator for coupled 1D–2D hydrodynamic models

Younghun Kang and Ethan J. Kubatko

Geosci. Model Dev., 17, 1603–1625, https://doi.org/10.5194/gmd-17-1603-2024,https://doi.org/10.5194/gmd-17-1603-2024, 2024

Short summary

Younghun Kang and Ethan J. Kubatko

Interactive discussion

Status: closed

CEC1:
'Comment on egusphere-2023-1434', Juan Antonio Añel, 13 Oct 2023

Dear authors,
Unfortunately, after checking your manuscript, it has come to our attention that it does not comply with our "Code and Data Policy".
https://www.geoscientific-model-development.net/policies/code_and_data_policy.html
You have archived your code on GitHub. However, GitHub is not a suitable repository for scientific publication. GitHub itself instructs authors to use other alternatives for long-term archival and publishing, such as Zenodo. Indeed, our policy (which we expect you to read before submitting your manuscript) clearly states that we can not accept GitHub. Therefore, please, publish your code in one of the appropriate repositories, and reply to this comment with the relevant information (link and DOI) as soon as possible, as it should be available before the Discussions stage, as the topical editor mentioned earlier.
Also, in the GitHub repository, there is no license listed for ADMESH+. If you do not include a license, the code continues to be your property, and nobody can use it. Therefore, when uploading the model's code to the new repository, you could want to choose a free software/open-source (FLOSS) license. We recommend the GPLv3. You only need to include the file 'https://www.gnu.org/licenses/gpl-3.0.txt' as LICENSE.txt with your code. Also, you can choose other options that Zenodo provides: GPLv2, Apache License, MIT License, etc.

In this way, if you do not fix the above-mentioned problems, we will have to reject your manuscript for publication in our journal. I should note that, actually, your manuscript should not have been accepted in Discussions, given this lack of compliance with our policy. Therefore, the current situation with your manuscript is irregular.

Also, you must include in a potentially reviewed version of your manuscript the modified 'Code and Data Availability' section, the DOI of the code (and another DOI for the dataset if necessary).

Juan A. Añel
Geosci. Model Dev. Executive Editor

Citation: https://doi.org/10.5194/egusphere-2023-1434-CEC1
- AC1: 'Reply on CEC1', Younghun Kang, 17 Oct 2023
  
  Dear Editor,
  I apologize for the issues coming from my lack of understanding for software license, and I really appreciate your throughout review of my manuscript and giving me a clear instruction of how I can proceed.
  As you guided, I included GPLv3 license and published the my codes in a new proper repository (Zenodo). Here is the link: https://zenodo.org/records/10009737 and DOI: 10.5281/zenodo.10009737.
  Finally, I revised code and data availability section and I will send the revised manuscript via email, as the system warns revised manuscripts should not be attached here as supplement.
  Again, I appreciate all of your help for this matter.
  Best Regards,
  
  Younghun Kang
  
  Citation: https://doi.org/10.5194/egusphere-2023-1434-AC1
- AC2: 'Reply on CEC1', Younghun Kang, 23 Oct 2023
  
  Dear Editor,
  Here is the modified code availability section:
  
  The current version of ADMESH+, the mesh generator presented in this study, is available as Zenodo archive: \url{https://zenodo.org/records/10009737}. The ADMESH+ is under active development and the latest version will be available upon request.
  Please let me know if you need any further information.
  Best regards,
  
  Younghun Kang
  
  Citation: https://doi.org/10.5194/egusphere-2023-1434-AC2
RC1:
'Comment on egusphere-2023-1434', Chris Massey, 28 Oct 2023

This paper on automatic mesh generator for coupled 1d/2d hydrodynamic models presents a new methodology for creating hybrid 1d/2d meshes. The paper is well written and the material presented in a clear manner. Some minor revisions and content additions would be helpful along with a few grammatical corrections.

Line 14 -- Drop "as" .... routinely used to simulate....
Lines 21 and 22 -- Provide a citation justifying or showing an example support for your statement about adequately described by simpler one-dimensional....
Line 67 -- add "an" ... models that is an extension....
Line 79... should also mention that an element size function is often supplied
Lines 156-157 -- this section should include a longer discussion on the importance of the quality of a DEM and a discussion on the sensitivity of the quality of the DEM to the results of ADMESH+.
Line 273 -- add "to" ...constraints too close to each other.
Line 310 -- Provide more details on how exactly ordering of the medial axis branches is done.
Lines 320-326 --- A longer discussion is needed on the sensitivity of the ADMESH+ algorithms to user setting/DEM quality for level 1 and level 2 regions and mis-identification of narrow features and the the omission of some level 1 land regions. How can one check/find when this is happening and what are solutions for correcting it when it does happen?
Line 339 -- Provide more explanation for the use of curvature of internal constraints for mesh size and why it is a good choice.
Line 360 -- The choice of 1 to 10 meters needs more explanation. Is it tied to the resolution and quality of the DEM and/or to the desired minimum element size?
Line 374 -- Remove second "." after Eq. 33.
Line 411 -- A question here? After removing 1D elements, does this require another round of adjustments to balance out the spacing?
Line 420 -- What is the distance (hmin/2) in relation to? The node's distance to the internal constraints? Suggest making this more clear.
Line 430 -- Explain why this particular mesh element quality measure was selected. Were other quality measures examined?
Line 482 -- The minimum mesh quality of 0.29 seems very low and would possibly lead to numerical model instabilities and or accuracy loss. Do the authors have a suggested minimum mesh quality to allow and if so, what should be done to fix elements that fall below that threshold?
Figure 18(b) shows that there are some nodes that have a high valence count, i.e. are connected to a large number of surrounding nodes (optimal would be 6 connected nodes). This high valance typically happens when element size transitions happen too rapidly. What is causing this in your algorithm? Are there catches/fixes that could be implemented during the overall process?
Lines 496-497 -- How are the successive line segments merged? What tool is used?
Lines 502 and 503 -- Discuss why K=20, grading limit 0.15, and RSME of 5 meters were selected and what impacts altering those values would have to the mesh construction.
Lines 516-517 ... Please add a zoomed in image showing an example area where the internal constraints are too close to one another and resulted in a low element quality. Also, discuss what a user would need to do in order to fix these issues.
Line 582 -- The sentence "Therefore, additional step is .... of VDT." does not seem to be well formed and is out of context. Maybe remove it, or add clarification.
Line 587 -- Suggestion adding if not already provided, that the sample generated meshes given in this paper also be shared alongside of the source code.

Citation: https://doi.org/10.5194/egusphere-2023-1434-RC1
- AC3: 'Reply on RC1', Younghun Kang, 03 Dec 2023
  
  Dear Reviewer,
  We sincerely appreciate the time and effort you dedicated to reviewing our manuscript titled "An automatic mesh generator for coupled 1D/2D hydrodynamic models." Your thoughtful review has been invaluable in enhancing the overall quality of our work. Your constructive feedback has significantly contributed to refining the manuscript, and we are grateful for your insightful comments.
  Here is a point-by-point response to the reviewers' comments and concerns:
  Line 14 -- Drop "as" .... routinely used to simulate....
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  Lines 21 and 22 -- Provide a citation justifying or showing an example support for your statement about adequately described by simpler one-dimensional....
  
  > We appreciate pointing out the need for citations. We added some citations and the line has been slightly edited for clarification.
  Line 67 -- add "an" ... models that is an extension....
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  Line 79... should also mention that an element size function is often supplied
  
  > We kindly note that the element size function is mentioned in the previous paragraph.
  Lines 156-157 -- this section should include a longer discussion on the importance of the quality of a DEM and a discussion on the sensitivity of the quality of the DEM to the results of ADMESH+.
  
  > We appreciate your suggestion. In response, we added a longer discussion on the quality of DEM and the extracted open-channels.
  Line 273 -- add "to" ...constraints too close to each other.
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  Line 310 -- Provide more details on how exactly ordering of the medial axis branches is done.
  
  > We kindly note that the details are described in Section 3.2.1, Step 2.
  Lines 320-326 --- A longer discussion is needed on the sensitivity of the ADMESH+ algorithms to user setting/DEM quality for level 1 and level 2 regions and mis-identification of narrow features and the the omission of some level 1 land regions. How can one check/find when this is happening and what are solutions for correcting it when it does happen?
  
  > We appreciate pointing out the need for a longer discussion of this topic. We have added potential solutions for the mis-classified level 1 regions.
  Line 339 -- Provide more explanation for the use of curvature of internal constraints for mesh size and why it is a good choice.
  
  > When we approximate internal constraints with linear elements, the errors will be higher for curved parts and lower for flat parts. It results in the need of small sized elements for curved parts. The use of curvature of internal constraints allows assigning a mesh size function such that its target mesh size is higher for curved parts and lower for flat parts.
  Line 360 -- The choice of 1 to 10 meters needs more explanation. Is it tied to the resolution and quality of the DEM and/or to the desired minimum element size?
  
  > As we discussed in the manuscript, the choice of 1 to 10 meters is based on our numerical experiments. It may relate to the DEM resolution rather than minimum element size, and it has been noted in the manuscript.
  Line 374 -- Remove second "." after Eq. 33.
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  Line 411 -- A question here? After removing 1D elements, does this require another round of adjustments to balance out the spacing?
  
  > By “removing 1D elements”, we meant merging 1D elements to their neighboring elements. The manuscript has been edited to clarify this. Also, note that it may be ideal to adaptively merge such elements (based on the lengths of left and right elements), but they are simply merged to left elements in the current version.
  Line 420 -- What is the distance (hmin/2) in relation to? The node's distance to the internal constraints? Suggest making this more clear.
  
  > It’s distance to the internal constraints, and the manuscript has been edited to make it clearer.
  Line 430 -- Explain why this particular mesh element quality measure was selected. Were other quality measures examined?
  
  > A longer discussion of the mesh element quality has been added in the manuscript. However, other element qualities have not been examined. We may add it for our future work.
  Line 482 -- The minimum mesh quality of 0.29 seems very low and would possibly lead to numerical model instabilities and or accuracy loss. Do the authors have a suggested minimum mesh quality to allow and if so, what should be done to fix elements that fall below that threshold?
  
  > Based on our previous work (Conroy et al. (2012) doi.org/10.1007/s10236-012-0574-0), we suggest minimum element quality 0.30 without internal constraints. However, we have not yet established a suggestible minimum element quality with internal constraints. An additional discussion for possible solution for the poor quality elements has been added in the conclusion section.
  Figure 18(b) shows that there are some nodes that have a high valence count, i.e. are connected to a large number of surrounding nodes (optimal would be 6 connected nodes). This high valance typically happens when element size transitions happen too rapidly. What is causing this in your algorithm? Are there catches/fixes that could be implemented during the overall process?
  
  > The high valance may be caused by the internal constraints; however, we may need additional examinations for a more accurate answer for this question. An additional discussion for possible solution for the poor quality elements has been added in the conclusion section.
  Lines 496-497 -- How are the successive line segments merged? What tool is used?
  
  > The modified shoreline data are provided by our colleague and we believe they are merged manually. For the exact information of which software has been used, we will need to reach out to our colleague.
  Lines 502 and 503 -- Discuss why K=20, grading limit 0.15, and RSME of 5 meters were selected and what impacts altering those values would have to the mesh construction.
  
  > These parameters are chosen as it gives “reasonable” results in our numerical experiments as shown in our previous work (see, Conroy2012) and this work.
  Lines 516-517 ... Please add a zoomed in image showing an example area where the internal constraints are too close to one another and resulted in a low element quality. Also, discuss what a user would need to do in order to fix these issues.
  
  > As suggested by the reviewer, we have added zoomed-in figures. As discussed in the conclusion section, it can be resolved by ignoring the object or allowing very small element sizes if the element quality is of higher priority.
  Line 582 -- The sentence "Therefore, additional step is .... of VDT." does not seem to be well formed and is out of context. Maybe remove it, or add clarification.
  
  > Thank you for pointing out the need for clarification. The manuscript has been edited to clarify it.
  Line 587 -- Suggestion adding if not already provided, that the sample generated meshes given in this paper also be shared alongside of the source code.
  
  > We appreciate the great suggestion. The sample meshes have been provided alongside the source code.
  
  Citation: https://doi.org/10.5194/egusphere-2023-1434-AC3
RC2:
'Comment on egusphere-2023-1434', Shintaro Bunya, 02 Nov 2023

This work proposes a novel approach to generating coupled 1D/2D unstructured meshes for hydrodynamic models utilizing accessible DEM datasets and water/land mask polygons. It presents detailed steps in the proposed procedures and some compelling examples. I suggest the following minor revisions.
P1-L15: It is not obvious how the sea-level rise itself is assessed using 2D hydrodynamic models. If this is meant to be considering the effect of sea level rise to coastal flooding, it may be better to rephrase it or take a different example.
P1-L20: "expensive" --> "expense"
P2-L60: A more recent work that is not mentioned here for internally-connected domains has been published: Bunya, et al. (2023) https://doi.org/10.1016/j.advengsoft.2023.103516.
P3-L91: "Section" --> "section"
P4-L105: It is not clear which work the authors refer to as "our previous work".
P6-L143: It would be helpful for readers if the definition of the global flow matrix M is presented in a brief manner.
P8-L184: The definition of "sharp" corners doesn't appear to be presented.
P8-L190: Some more description of "background grid points" should be added here or somewhere earlier. Otherwise, it's not clear how they are defined.
P13-L252 and L253: Can these equations be expressed using symbols other than the leftwards arrows?
P13-L264: Present the definitions of "Perimter" (typo?) and "Area".
P15-L299: "the updated level 1 mask" <-- Is this land mask or water mask?
P16 or somewhere earlier: How are the medial axes points connected to construct line segments?
P17 or somewhere earlier: The force equilibrium algorithm used in this work can be briefly presented as it plays a central role in Section 4.
P19-L393: "... there should be at least three 1D elements along the boundaries of the level 2 regions. ..." <--- This may require some more elaboration for clarity.
P19-L394: "... identified, as level 2, regions, ..." <--- Check if the commas are necessary.
P20-L404: "... then the 1D meshes are removed." <--- Are both 1D meshes removed? Or, do you leave one of them?
For three test cases in Section 5: The minimum and maximum element sizes of the resulting meshes should be presented to show how they coincide with or differ from the specified values.

Citation: https://doi.org/10.5194/egusphere-2023-1434-RC2
- AC4: 'Reply on RC2', Younghun Kang, 03 Dec 2023
  
  Dear Reviewer,
  We sincerely appreciate the time and effort you dedicated to reviewing our manuscript titled "An automatic mesh generator for coupled 1D/2D hydrodynamic models." Your thoughtful and constructive feedback has been invaluable in enhancing the quality of our work.
  We appreciate your insightful comments that have undoubtedly strengthened the overall manuscript. Your keen observations and suggestions have played a crucial role in refining our research. Here is a point-by-point response to the reviewers’ comments and concerns:
  P1-L15: It is not obvious how the sea-level rise itself is assessed using 2D hydrodynamic models. If this is meant to be considering the effect of sea level rise to coastal flooding, it may be better to rephrase it or take a different example.
  
  > Thank you for pointing out this is not an obvious example. We have taken another example.
  P1-L20: "expensive" --> "expense"
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  P2-L60: A more recent work that is not mentioned here for internally-connected domains has been published: Bunya, et al. (2023) https://doi.org/10.1016/j.advengsoft.2023.103516.
  
  > Thank you for pointing out that the good reference has been missed. It has been added in the previous paragraph (P2-L45) than the reviewer mentioned, as we believe it would be a better place to add this citation.
  P3-L91: "Section" --> "section"
  
  > Thank you for pointing out the grammar error. It has been fixed.
  P4-L105: It is not clear which work the authors refer to as "our previous work".
  
  > Thank you for pointing out the unclarity. We have added the citation of the work.
  P6-L143: It would be helpful for readers if the definition of the global flow matrix M is presented in a brief manner.
  
  > We appreciate the suggestion for better clarity. The definition has been added as suggested.
  P8-L184: The definition of "sharp" corners doesn't appear to be presented.
  
  > Thank you for pointing this out. The sharpness is related to the threshold parameter \delta_{\theta}. However, we decided to remove the word “sharp” as it may cause confusion.
  P8-L190: Some more description of "background grid points" should be added here or somewhere earlier. Otherwise, it's not clear how they are defined.
  
  > We kindly note that the background grid points are described in Section 2.
  P13-L252 and L253: Can these equations be expressed using symbols other than the leftwards arrows?
  
  > As suggested, the symbol has been replaced with colon equal. Furthermore, for better clarity, the notations of masks are modified.
  P13-L264: Present the definitions of "Perimter" (typo?) and "Area".
  
  > Thank you for pointing out the typo. It has been fixed. Also, the “perimeter” and “area” are standard measure. In order to avoid confusion, they are replaced to lower cases.
  P15-L299: "the updated level 1 mask" <-- Is this land mask or water mask?
  
  > It is the level 1 mask from the updated water mask. It has been clarified in the manuscript. Also, as responded to the earlier comment, the notations of masks are modified for better clarity.
  P16 or somewhere earlier: How are the medial axes points connected to construct line segments?
  
  > We appreciate pointing out that the description of the construction of medial axis branches are missed in the manuscript. The description has been added in the appendix.
  P17 or somewhere earlier: The force equilibrium algorithm used in this work can be briefly presented as it plays a central role in Section 4.
  
  > Thank you for the suggestion. A brief description of the force equilibrium algorithm has been added in the manuscript.
  P19-L393: "... there should be at least three 1D elements along the boundaries of the level 2 regions. ..." <--- This may require some more elaboration for clarity.
  
  > Thank you for pointing out the unclarity. It has been edited for clarity.
  P19-L394: "... identified, as level 2, regions, ..." <--- Check if the commas are necessary.
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  P20-L404: "... then the 1D meshes are removed." <--- Are both 1D meshes removed? Or, do you leave one of them?
  
  > 1D meshes for open-channel/internal boundaries are removed. For better clarity, it has been edited.
  For three test cases in Section 5: The minimum and maximum element sizes of the resulting meshes should be presented to show how they coincide with or differ from the specified values.
  
  > As suggested, the minimum and maximum element sizes have been presented for the test cases.
  
  Citation: https://doi.org/10.5194/egusphere-2023-1434-AC4

Interactive discussion

Status: closed

CEC1:
'Comment on egusphere-2023-1434', Juan Antonio Añel, 13 Oct 2023

Dear authors,
Unfortunately, after checking your manuscript, it has come to our attention that it does not comply with our "Code and Data Policy".
https://www.geoscientific-model-development.net/policies/code_and_data_policy.html
You have archived your code on GitHub. However, GitHub is not a suitable repository for scientific publication. GitHub itself instructs authors to use other alternatives for long-term archival and publishing, such as Zenodo. Indeed, our policy (which we expect you to read before submitting your manuscript) clearly states that we can not accept GitHub. Therefore, please, publish your code in one of the appropriate repositories, and reply to this comment with the relevant information (link and DOI) as soon as possible, as it should be available before the Discussions stage, as the topical editor mentioned earlier.
Also, in the GitHub repository, there is no license listed for ADMESH+. If you do not include a license, the code continues to be your property, and nobody can use it. Therefore, when uploading the model's code to the new repository, you could want to choose a free software/open-source (FLOSS) license. We recommend the GPLv3. You only need to include the file 'https://www.gnu.org/licenses/gpl-3.0.txt' as LICENSE.txt with your code. Also, you can choose other options that Zenodo provides: GPLv2, Apache License, MIT License, etc.

In this way, if you do not fix the above-mentioned problems, we will have to reject your manuscript for publication in our journal. I should note that, actually, your manuscript should not have been accepted in Discussions, given this lack of compliance with our policy. Therefore, the current situation with your manuscript is irregular.

Also, you must include in a potentially reviewed version of your manuscript the modified 'Code and Data Availability' section, the DOI of the code (and another DOI for the dataset if necessary).

Juan A. Añel
Geosci. Model Dev. Executive Editor

Citation: https://doi.org/10.5194/egusphere-2023-1434-CEC1
- AC1: 'Reply on CEC1', Younghun Kang, 17 Oct 2023
  
  Dear Editor,
  I apologize for the issues coming from my lack of understanding for software license, and I really appreciate your throughout review of my manuscript and giving me a clear instruction of how I can proceed.
  As you guided, I included GPLv3 license and published the my codes in a new proper repository (Zenodo). Here is the link: https://zenodo.org/records/10009737 and DOI: 10.5281/zenodo.10009737.
  Finally, I revised code and data availability section and I will send the revised manuscript via email, as the system warns revised manuscripts should not be attached here as supplement.
  Again, I appreciate all of your help for this matter.
  Best Regards,
  
  Younghun Kang
  
  Citation: https://doi.org/10.5194/egusphere-2023-1434-AC1
- AC2: 'Reply on CEC1', Younghun Kang, 23 Oct 2023
  
  Dear Editor,
  Here is the modified code availability section:
  
  The current version of ADMESH+, the mesh generator presented in this study, is available as Zenodo archive: \url{https://zenodo.org/records/10009737}. The ADMESH+ is under active development and the latest version will be available upon request.
  Please let me know if you need any further information.
  Best regards,
  
  Younghun Kang
  
  Citation: https://doi.org/10.5194/egusphere-2023-1434-AC2
RC1:
'Comment on egusphere-2023-1434', Chris Massey, 28 Oct 2023

This paper on automatic mesh generator for coupled 1d/2d hydrodynamic models presents a new methodology for creating hybrid 1d/2d meshes. The paper is well written and the material presented in a clear manner. Some minor revisions and content additions would be helpful along with a few grammatical corrections.

Line 14 -- Drop "as" .... routinely used to simulate....
Lines 21 and 22 -- Provide a citation justifying or showing an example support for your statement about adequately described by simpler one-dimensional....
Line 67 -- add "an" ... models that is an extension....
Line 79... should also mention that an element size function is often supplied
Lines 156-157 -- this section should include a longer discussion on the importance of the quality of a DEM and a discussion on the sensitivity of the quality of the DEM to the results of ADMESH+.
Line 273 -- add "to" ...constraints too close to each other.
Line 310 -- Provide more details on how exactly ordering of the medial axis branches is done.
Lines 320-326 --- A longer discussion is needed on the sensitivity of the ADMESH+ algorithms to user setting/DEM quality for level 1 and level 2 regions and mis-identification of narrow features and the the omission of some level 1 land regions. How can one check/find when this is happening and what are solutions for correcting it when it does happen?
Line 339 -- Provide more explanation for the use of curvature of internal constraints for mesh size and why it is a good choice.
Line 360 -- The choice of 1 to 10 meters needs more explanation. Is it tied to the resolution and quality of the DEM and/or to the desired minimum element size?
Line 374 -- Remove second "." after Eq. 33.
Line 411 -- A question here? After removing 1D elements, does this require another round of adjustments to balance out the spacing?
Line 420 -- What is the distance (hmin/2) in relation to? The node's distance to the internal constraints? Suggest making this more clear.
Line 430 -- Explain why this particular mesh element quality measure was selected. Were other quality measures examined?
Line 482 -- The minimum mesh quality of 0.29 seems very low and would possibly lead to numerical model instabilities and or accuracy loss. Do the authors have a suggested minimum mesh quality to allow and if so, what should be done to fix elements that fall below that threshold?
Figure 18(b) shows that there are some nodes that have a high valence count, i.e. are connected to a large number of surrounding nodes (optimal would be 6 connected nodes). This high valance typically happens when element size transitions happen too rapidly. What is causing this in your algorithm? Are there catches/fixes that could be implemented during the overall process?
Lines 496-497 -- How are the successive line segments merged? What tool is used?
Lines 502 and 503 -- Discuss why K=20, grading limit 0.15, and RSME of 5 meters were selected and what impacts altering those values would have to the mesh construction.
Lines 516-517 ... Please add a zoomed in image showing an example area where the internal constraints are too close to one another and resulted in a low element quality. Also, discuss what a user would need to do in order to fix these issues.
Line 582 -- The sentence "Therefore, additional step is .... of VDT." does not seem to be well formed and is out of context. Maybe remove it, or add clarification.
Line 587 -- Suggestion adding if not already provided, that the sample generated meshes given in this paper also be shared alongside of the source code.

Citation: https://doi.org/10.5194/egusphere-2023-1434-RC1
- AC3: 'Reply on RC1', Younghun Kang, 03 Dec 2023
  
  Dear Reviewer,
  We sincerely appreciate the time and effort you dedicated to reviewing our manuscript titled "An automatic mesh generator for coupled 1D/2D hydrodynamic models." Your thoughtful review has been invaluable in enhancing the overall quality of our work. Your constructive feedback has significantly contributed to refining the manuscript, and we are grateful for your insightful comments.
  Here is a point-by-point response to the reviewers' comments and concerns:
  Line 14 -- Drop "as" .... routinely used to simulate....
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  Lines 21 and 22 -- Provide a citation justifying or showing an example support for your statement about adequately described by simpler one-dimensional....
  
  > We appreciate pointing out the need for citations. We added some citations and the line has been slightly edited for clarification.
  Line 67 -- add "an" ... models that is an extension....
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  Line 79... should also mention that an element size function is often supplied
  
  > We kindly note that the element size function is mentioned in the previous paragraph.
  Lines 156-157 -- this section should include a longer discussion on the importance of the quality of a DEM and a discussion on the sensitivity of the quality of the DEM to the results of ADMESH+.
  
  > We appreciate your suggestion. In response, we added a longer discussion on the quality of DEM and the extracted open-channels.
  Line 273 -- add "to" ...constraints too close to each other.
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  Line 310 -- Provide more details on how exactly ordering of the medial axis branches is done.
  
  > We kindly note that the details are described in Section 3.2.1, Step 2.
  Lines 320-326 --- A longer discussion is needed on the sensitivity of the ADMESH+ algorithms to user setting/DEM quality for level 1 and level 2 regions and mis-identification of narrow features and the the omission of some level 1 land regions. How can one check/find when this is happening and what are solutions for correcting it when it does happen?
  
  > We appreciate pointing out the need for a longer discussion of this topic. We have added potential solutions for the mis-classified level 1 regions.
  Line 339 -- Provide more explanation for the use of curvature of internal constraints for mesh size and why it is a good choice.
  
  > When we approximate internal constraints with linear elements, the errors will be higher for curved parts and lower for flat parts. It results in the need of small sized elements for curved parts. The use of curvature of internal constraints allows assigning a mesh size function such that its target mesh size is higher for curved parts and lower for flat parts.
  Line 360 -- The choice of 1 to 10 meters needs more explanation. Is it tied to the resolution and quality of the DEM and/or to the desired minimum element size?
  
  > As we discussed in the manuscript, the choice of 1 to 10 meters is based on our numerical experiments. It may relate to the DEM resolution rather than minimum element size, and it has been noted in the manuscript.
  Line 374 -- Remove second "." after Eq. 33.
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  Line 411 -- A question here? After removing 1D elements, does this require another round of adjustments to balance out the spacing?
  
  > By “removing 1D elements”, we meant merging 1D elements to their neighboring elements. The manuscript has been edited to clarify this. Also, note that it may be ideal to adaptively merge such elements (based on the lengths of left and right elements), but they are simply merged to left elements in the current version.
  Line 420 -- What is the distance (hmin/2) in relation to? The node's distance to the internal constraints? Suggest making this more clear.
  
  > It’s distance to the internal constraints, and the manuscript has been edited to make it clearer.
  Line 430 -- Explain why this particular mesh element quality measure was selected. Were other quality measures examined?
  
  > A longer discussion of the mesh element quality has been added in the manuscript. However, other element qualities have not been examined. We may add it for our future work.
  Line 482 -- The minimum mesh quality of 0.29 seems very low and would possibly lead to numerical model instabilities and or accuracy loss. Do the authors have a suggested minimum mesh quality to allow and if so, what should be done to fix elements that fall below that threshold?
  
  > Based on our previous work (Conroy et al. (2012) doi.org/10.1007/s10236-012-0574-0), we suggest minimum element quality 0.30 without internal constraints. However, we have not yet established a suggestible minimum element quality with internal constraints. An additional discussion for possible solution for the poor quality elements has been added in the conclusion section.
  Figure 18(b) shows that there are some nodes that have a high valence count, i.e. are connected to a large number of surrounding nodes (optimal would be 6 connected nodes). This high valance typically happens when element size transitions happen too rapidly. What is causing this in your algorithm? Are there catches/fixes that could be implemented during the overall process?
  
  > The high valance may be caused by the internal constraints; however, we may need additional examinations for a more accurate answer for this question. An additional discussion for possible solution for the poor quality elements has been added in the conclusion section.
  Lines 496-497 -- How are the successive line segments merged? What tool is used?
  
  > The modified shoreline data are provided by our colleague and we believe they are merged manually. For the exact information of which software has been used, we will need to reach out to our colleague.
  Lines 502 and 503 -- Discuss why K=20, grading limit 0.15, and RSME of 5 meters were selected and what impacts altering those values would have to the mesh construction.
  
  > These parameters are chosen as it gives “reasonable” results in our numerical experiments as shown in our previous work (see, Conroy2012) and this work.
  Lines 516-517 ... Please add a zoomed in image showing an example area where the internal constraints are too close to one another and resulted in a low element quality. Also, discuss what a user would need to do in order to fix these issues.
  
  > As suggested by the reviewer, we have added zoomed-in figures. As discussed in the conclusion section, it can be resolved by ignoring the object or allowing very small element sizes if the element quality is of higher priority.
  Line 582 -- The sentence "Therefore, additional step is .... of VDT." does not seem to be well formed and is out of context. Maybe remove it, or add clarification.
  
  > Thank you for pointing out the need for clarification. The manuscript has been edited to clarify it.
  Line 587 -- Suggestion adding if not already provided, that the sample generated meshes given in this paper also be shared alongside of the source code.
  
  > We appreciate the great suggestion. The sample meshes have been provided alongside the source code.
  
  Citation: https://doi.org/10.5194/egusphere-2023-1434-AC3
RC2:
'Comment on egusphere-2023-1434', Shintaro Bunya, 02 Nov 2023

This work proposes a novel approach to generating coupled 1D/2D unstructured meshes for hydrodynamic models utilizing accessible DEM datasets and water/land mask polygons. It presents detailed steps in the proposed procedures and some compelling examples. I suggest the following minor revisions.
P1-L15: It is not obvious how the sea-level rise itself is assessed using 2D hydrodynamic models. If this is meant to be considering the effect of sea level rise to coastal flooding, it may be better to rephrase it or take a different example.
P1-L20: "expensive" --> "expense"
P2-L60: A more recent work that is not mentioned here for internally-connected domains has been published: Bunya, et al. (2023) https://doi.org/10.1016/j.advengsoft.2023.103516.
P3-L91: "Section" --> "section"
P4-L105: It is not clear which work the authors refer to as "our previous work".
P6-L143: It would be helpful for readers if the definition of the global flow matrix M is presented in a brief manner.
P8-L184: The definition of "sharp" corners doesn't appear to be presented.
P8-L190: Some more description of "background grid points" should be added here or somewhere earlier. Otherwise, it's not clear how they are defined.
P13-L252 and L253: Can these equations be expressed using symbols other than the leftwards arrows?
P13-L264: Present the definitions of "Perimter" (typo?) and "Area".
P15-L299: "the updated level 1 mask" <-- Is this land mask or water mask?
P16 or somewhere earlier: How are the medial axes points connected to construct line segments?
P17 or somewhere earlier: The force equilibrium algorithm used in this work can be briefly presented as it plays a central role in Section 4.
P19-L393: "... there should be at least three 1D elements along the boundaries of the level 2 regions. ..." <--- This may require some more elaboration for clarity.
P19-L394: "... identified, as level 2, regions, ..." <--- Check if the commas are necessary.
P20-L404: "... then the 1D meshes are removed." <--- Are both 1D meshes removed? Or, do you leave one of them?
For three test cases in Section 5: The minimum and maximum element sizes of the resulting meshes should be presented to show how they coincide with or differ from the specified values.

Citation: https://doi.org/10.5194/egusphere-2023-1434-RC2
- AC4: 'Reply on RC2', Younghun Kang, 03 Dec 2023
  
  Dear Reviewer,
  We sincerely appreciate the time and effort you dedicated to reviewing our manuscript titled "An automatic mesh generator for coupled 1D/2D hydrodynamic models." Your thoughtful and constructive feedback has been invaluable in enhancing the quality of our work.
  We appreciate your insightful comments that have undoubtedly strengthened the overall manuscript. Your keen observations and suggestions have played a crucial role in refining our research. Here is a point-by-point response to the reviewers’ comments and concerns:
  P1-L15: It is not obvious how the sea-level rise itself is assessed using 2D hydrodynamic models. If this is meant to be considering the effect of sea level rise to coastal flooding, it may be better to rephrase it or take a different example.
  
  > Thank you for pointing out this is not an obvious example. We have taken another example.
  P1-L20: "expensive" --> "expense"
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  P2-L60: A more recent work that is not mentioned here for internally-connected domains has been published: Bunya, et al. (2023) https://doi.org/10.1016/j.advengsoft.2023.103516.
  
  > Thank you for pointing out that the good reference has been missed. It has been added in the previous paragraph (P2-L45) than the reviewer mentioned, as we believe it would be a better place to add this citation.
  P3-L91: "Section" --> "section"
  
  > Thank you for pointing out the grammar error. It has been fixed.
  P4-L105: It is not clear which work the authors refer to as "our previous work".
  
  > Thank you for pointing out the unclarity. We have added the citation of the work.
  P6-L143: It would be helpful for readers if the definition of the global flow matrix M is presented in a brief manner.
  
  > We appreciate the suggestion for better clarity. The definition has been added as suggested.
  P8-L184: The definition of "sharp" corners doesn't appear to be presented.
  
  > Thank you for pointing this out. The sharpness is related to the threshold parameter \delta_{\theta}. However, we decided to remove the word “sharp” as it may cause confusion.
  P8-L190: Some more description of "background grid points" should be added here or somewhere earlier. Otherwise, it's not clear how they are defined.
  
  > We kindly note that the background grid points are described in Section 2.
  P13-L252 and L253: Can these equations be expressed using symbols other than the leftwards arrows?
  
  > As suggested, the symbol has been replaced with colon equal. Furthermore, for better clarity, the notations of masks are modified.
  P13-L264: Present the definitions of "Perimter" (typo?) and "Area".
  
  > Thank you for pointing out the typo. It has been fixed. Also, the “perimeter” and “area” are standard measure. In order to avoid confusion, they are replaced to lower cases.
  P15-L299: "the updated level 1 mask" <-- Is this land mask or water mask?
  
  > It is the level 1 mask from the updated water mask. It has been clarified in the manuscript. Also, as responded to the earlier comment, the notations of masks are modified for better clarity.
  P16 or somewhere earlier: How are the medial axes points connected to construct line segments?
  
  > We appreciate pointing out that the description of the construction of medial axis branches are missed in the manuscript. The description has been added in the appendix.
  P17 or somewhere earlier: The force equilibrium algorithm used in this work can be briefly presented as it plays a central role in Section 4.
  
  > Thank you for the suggestion. A brief description of the force equilibrium algorithm has been added in the manuscript.
  P19-L393: "... there should be at least three 1D elements along the boundaries of the level 2 regions. ..." <--- This may require some more elaboration for clarity.
  
  > Thank you for pointing out the unclarity. It has been edited for clarity.
  P19-L394: "... identified, as level 2, regions, ..." <--- Check if the commas are necessary.
  
  > Thank you for pointing out the grammatical error. It has been fixed.
  P20-L404: "... then the 1D meshes are removed." <--- Are both 1D meshes removed? Or, do you leave one of them?
  
  > 1D meshes for open-channel/internal boundaries are removed. For better clarity, it has been edited.
  For three test cases in Section 5: The minimum and maximum element sizes of the resulting meshes should be presented to show how they coincide with or differ from the specified values.
  
  > As suggested, the minimum and maximum element sizes have been presented for the test cases.
  
  Citation: https://doi.org/10.5194/egusphere-2023-1434-AC4

Peer review completion

AR: Author's response | RR: Referee report | ED: Editor decision | EF: Editorial file upload

AR by Younghun Kang on behalf of the Authors (18 Dec 2023) Author's response Author's tracked changes Manuscript

ED: Publish as is (19 Dec 2023) by Simone Marras

AR by Younghun Kang on behalf of the Authors (06 Jan 2024)

Journal article(s) based on this preprint

23 Feb 2024

An automatic mesh generator for coupled 1D–2D hydrodynamic models

Younghun Kang and Ethan J. Kubatko

Geosci. Model Dev., 17, 1603–1625, https://doi.org/10.5194/gmd-17-1603-2024,https://doi.org/10.5194/gmd-17-1603-2024, 2024

Short summary

Younghun Kang and Ethan J. Kubatko

Model code and software

ADMESH+ Younghun Kang, Ethan J. Kubatko, Colton J. Conroy, and Dustin W. West https://github.com/younghun-kang/ADMESH

Younghun Kang and Ethan J. Kubatko

Viewed

Total article views: 526 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	BibTeX	EndNote
397	105	24	526	10	12

HTML: 397
PDF: 105
XML: 24
Total: 526
BibTeX: 10
EndNote: 12

Views and downloads (calculated since 29 Sep 2023)

Month	HTML	PDF	XML	Total
Sep 2023	50	21	2	73
Oct 2023	197	44	12	253
Nov 2023	63	11	3	77
Dec 2023	42	14	5	61
Jan 2024	32	5	1	38
Feb 2024	13	10	1	24

Cumulative views and downloads (calculated since 29 Sep 2023)

Month	HTML	PDF	XML	Total
Sep 2023	50	21	2	73
Oct 2023	197	44	12	253
Nov 2023	63	11	3	77
Dec 2023	42	14	5	61
Jan 2024	32	5	1	38
Feb 2024	13	10	1	24

Viewed (geographical distribution)

Total article views: 510 (including HTML, PDF, and XML) Thereof 510 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 23 Feb 2024

Download

The requested preprint has a corresponding peer-reviewed final revised paper. You are encouraged to refer to the final revised version.

Preprint (56038 KB)
Metadata XML

Short summary

Models used to simulate the flow of coastal and riverine waters, including flooding, require a geometric representation (or mesh) of geographic features that exhibit a range of disparate spatial scales—from large open waters to small, narrow channels. Representing these features in an accurate way without excessive computational overhead presents a challenge. Here, we develop an automatic mesh generation tool to help address this challenge. Our results demonstrate the efficacy of our approach.


Total:	0
HTML:	0
PDF:	0
XML:	0