Localized MQ-RBF meshless techniques for modeling unsaturated flow

• Published in: Engineering analysis with boundary elements Vol. 130; pp. 109 - 123
• Main Authors: Boujoudar, Mohamed, Beljadid, Abdelaziz, Taik, Ahmed
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01-09-2021
• Subjects: Infiltration; Meshless methods; Porous media; Radial basis function; Richards equation; Space-time techniques; Radial basis function; Porous media; Space-time techniques; Infiltration; Meshless methods; Richards equation
• ISSN: 0955-7997; 1873-197X
• DOI: 10.1016/j.enganabound.2021.05.011

•We develop space-time mesh-free numerical techniques to solve the Richards equation.•The proposed numerical method is efficient in terms of computational cost.•The method is able to predict the dynamics of flows through unsaturated porous media.•Various numerical examples confirm the efficiency and robustness of our method.
In this study, we focus on space-time mesh-free numerical techniques for efficiently solving the Richards equation which is often used to model unsaturated flow through porous media. We propose an efficient approach which combines the use of local multiquadric (MQ) radial basis function (RBF) methods and space-time techniques. The localized MQ-RBFs meshless methods allow to avoid mesh generation and ill-conditioning problem where a sparse matrix is obtained for the global system which has the advantage of using reduced memory and computational time. To further reduce the computational cost, we use the space-time techniques having the advantages of solving the resulting algebraic system only once and removing the time-integration procedure. The proposed method has the benefit of considering collocation points on the boundaries of computational domains which makes it more flexible in dealing with complex geometries. We implement the proposed numerical model of infiltration and we perform a series of numerical tests, encompassing various nontrivial solutions, to confirm the performance of the proposed techniques. The numerical simulations show the accuracy, efficiency in terms of computational cost, and capability of the proposed numerical techniques in solving the Richards equation in two-, three- and four-dimensional space-time domains with complex boundaries.