Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition
Research on seepage flow in the vadose zone has largely been driven by engineering and environmental problems affecting many fields of geotechnics, hydrology, and agricultural science. Mathematical modeling of the subsurface flow under unsaturated conditions is an essential part of water resource ma...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
MDPI AG
2020-06-01
|
| Series: | Water |
| Subjects: | |
| Online Access: | https://www.mdpi.com/2073-4441/12/6/1780 |
| _version_ | 1797564350336073728 |
|---|---|
| author | Dariusz Gąsiorowski Tomasz Kolerski |
| author_facet | Dariusz Gąsiorowski Tomasz Kolerski |
| author_sort | Dariusz Gąsiorowski |
| collection | DOAJ |
| description | Research on seepage flow in the vadose zone has largely been driven by engineering and environmental problems affecting many fields of geotechnics, hydrology, and agricultural science. Mathematical modeling of the subsurface flow under unsaturated conditions is an essential part of water resource management and planning. In order to determine such subsurface flow, the two-dimensional (2D) Richards equation can be used. However, the computation process is often hampered by a high spatial resolution and long simulation period as well as the non-linearity of the equation. A new highly efficient and accurate method for solving the 2D Richards equation has been proposed in the paper. The developed algorithm is based on dimensional splitting, the result of which means that 1D equations can be solved more efficiently than as is the case with unsplit 2D algorithms. Moreover, such a splitting approach allows any algorithm to be used for space as well as time approximation, which in turn increases the accuracy of the numerical solution. The robustness and advantages of the proposed algorithms have been proven by two numerical tests representing typical engineering problems and performed for typical properties of soil. |
| first_indexed | 2024-03-10T18:57:05Z |
| format | Article |
| id | doaj.art-258487c013314103a068508bf7627927 |
| institution | Directory Open Access Journal |
| issn | 2073-4441 |
| language | English |
| last_indexed | 2024-03-10T18:57:05Z |
| publishDate | 2020-06-01 |
| publisher | MDPI AG |
| record_format | Article |
| series | Water |
| spelling | doaj.art-258487c013314103a068508bf76279272023-11-20T04:40:48ZengMDPI AGWater2073-44412020-06-01126178010.3390/w12061780Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional DecompositionDariusz Gąsiorowski0Tomasz Kolerski1Faculty of Civil and Environmental Engineering, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, PolandFaculty of Civil and Environmental Engineering, Gdańsk University of Technology, Narutowicza 11/12, 80-233 Gdańsk, PolandResearch on seepage flow in the vadose zone has largely been driven by engineering and environmental problems affecting many fields of geotechnics, hydrology, and agricultural science. Mathematical modeling of the subsurface flow under unsaturated conditions is an essential part of water resource management and planning. In order to determine such subsurface flow, the two-dimensional (2D) Richards equation can be used. However, the computation process is often hampered by a high spatial resolution and long simulation period as well as the non-linearity of the equation. A new highly efficient and accurate method for solving the 2D Richards equation has been proposed in the paper. The developed algorithm is based on dimensional splitting, the result of which means that 1D equations can be solved more efficiently than as is the case with unsplit 2D algorithms. Moreover, such a splitting approach allows any algorithm to be used for space as well as time approximation, which in turn increases the accuracy of the numerical solution. The robustness and advantages of the proposed algorithms have been proven by two numerical tests representing typical engineering problems and performed for typical properties of soil.https://www.mdpi.com/2073-4441/12/6/17802D Richards equationsubsurface flowunsaturated flowdimensional splittingGodunov methodStrang method |
| spellingShingle | Dariusz Gąsiorowski Tomasz Kolerski Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition Water 2D Richards equation subsurface flow unsaturated flow dimensional splitting Godunov method Strang method |
| title | Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition |
| title_full | Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition |
| title_fullStr | Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition |
| title_full_unstemmed | Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition |
| title_short | Numerical Solution of the Two-Dimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition |
| title_sort | numerical solution of the two dimensional richards equation using alternate splitting methods for dimensional decomposition |
| topic | 2D Richards equation subsurface flow unsaturated flow dimensional splitting Godunov method Strang method |
| url | https://www.mdpi.com/2073-4441/12/6/1780 |
| work_keys_str_mv | AT dariuszgasiorowski numericalsolutionofthetwodimensionalrichardsequationusingalternatesplittingmethodsfordimensionaldecomposition AT tomaszkolerski numericalsolutionofthetwodimensionalrichardsequationusingalternatesplittingmethodsfordimensionaldecomposition |