the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 14 Oct 2024
Technical note: Finite element formulations to map discrete fracture elements in three-dimensional groundwater models
Abstract. Typically, in finite element groundwater models, fractures are represented by two-dimensional triangular or quadrilateral elements. When embedded in a three-dimensional space, the Jacobian matrix governing the transformation from the global three-dimensional space to the local two-dimensional space is rectangular and thus not invertible. There exist different approaches to obtain a unique mapping from local to global space even though the Jacobian matrix is not invertible. These approaches are discussed in this study. It is illustrated that all approaches yield the same result and may be applied to curved elements. The mapping of anisotropic hydraulic conductivity tensors for possibly curved fracture elements is also discussed.
- Preprint
(498 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on hess-2024-289', Anonymous Referee #1, 24 Dec 2024
This manuscript presents a methodology for including two-dimensional elements, such as fractures, in a three-dimensional porous media matrix. The authors claim that the proposed methodology is easier to implement than the commonly adopted method composed of a two-dimensional mapping on an orthonormal space, followed by a three-dimensional rotation. In addition, the proposed methodology works for curved two-dimensional surfaces with the integration performed using the Gaussian quadrature method, a case that cannot be addressed with the existing methodology. The proposed new approach cannot be applied when the integration is performed analytically (in cases where this is possible).Â
The technical note is "dry" in the sense that it discusses only the new projection method showing how it works only for one of the integrals (Eq. 10) that should be computed in the FEM solution of the flow equation. The discussion on the hydraulic conductivity projection of the 2-D elements (the fractures) presented in Section 6 is appreciable, but the manuscript lacks an illustrative example, which would be certainly appreciated by the HESS readership and myself too.Â
The proposed methodology is sound. However, as mentioned above, the presentation is somewhat "dry," and including an application example would enhance the reader's understanding of the benefits offered by the proposed approach. The example does not need to be geometrically complex; for instance, a cubic homogeneous formation with a single curvilinear fracture, an imposed head gradient between two opposing faces, and the remaining faces set as impervious would be sufficient to demonstrate the methodology's capabilities. In addition, I suggest that the authors compare the solution obtained using their proposed approach with that obtained through the existing and widely adopted projection and rotation methodology.
Â
Citation: https://doi.org/10.5194/hess-2024-289-RC1 -
AC1: 'Reply on RC1', Rob de Rooij, 02 Jan 2025
I agree with the assessment that the manuscript is somewhat 'dry'. However, this being a technical note this is by design. Nonetheless, I agree that an illustrative example could be a useful addition. The reason why I did not include one was to keep the manuscript concise.Â
If an illustrative example should be added, I think that an even simpler example would suffice (simpler than the one proposed by the reviewer). Namely, if we take a single flat triangular element in 3D space with a given hydraulic conductivity tensor, then the example could work out the element conductance matrix (Equation 10), using the different mapping techniques.Â
Â
Â
Citation: https://doi.org/10.5194/hess-2024-289-AC1
-
AC1: 'Reply on RC1', Rob de Rooij, 02 Jan 2025
-
RC2: 'Comment on hess-2024-289', Anonymous Referee #2, 16 Jan 2025
The authors propose several approaches to map local spaces of fractures embedded in higher-dimensional domain. My main concern is related to the fact that the proposed approaches should be numerically validate, showing their performances for complex fracture networks. This makes the reader difficult to evaluate which algorithm should be preferred. I have also the following minor comments:
- eq (3) use a different symbol since it resembles a partial derivativeÂ
- the part on the Penrose-Moore inverse can be moved to an appendixCitation: https://doi.org/10.5194/hess-2024-289-RC2
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 146 | 41 | 76 | 263 | 9 | 8 |
- HTML: 146
- PDF: 41
- XML: 76
- Total: 263
- BibTeX: 9
- EndNote: 8
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1