Premium
A generalised finite difference scheme based on compact integrated radial basis function for flow in heterogeneous soils
Author(s) -
NgoCong D.,
Tien C. M. T.,
NguyenKy T.,
AnVo D.A.,
MaiDuy N.,
Strunin D. V.,
TranCong T.
Publication year - 2017
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.4386
Subject(s) - stencil , mathematics , radial basis function , finite difference , basis function , order of accuracy , poisson's equation , taylor series , nonlinear system , partial differential equation , function (biology) , convergence (economics) , finite difference method , grid , ordinary differential equation , mathematical optimization , flow (mathematics) , mathematical analysis , differential equation , computer science , method of characteristics , geometry , computational science , physics , quantum mechanics , machine learning , evolutionary biology , artificial neural network , economics , biology , economic growth
Summary In the present paper, we develop a generalised finite difference approach based on compact integrated radial basis function (CIRBF) stencils for solving highly nonlinear Richards equation governing fluid movement in heterogeneous soils. The proposed CIRBF scheme enjoys a high level of accuracy and a fast convergence rate with grid refinement owing to the combination of the integrated RBF approximation and compact approximation where the spatial derivatives are discretised in terms of the information of neighbouring nodes in a stencil. The CIRBF method is first verified through the solution of ordinary differential equations, 2–D Poisson equations and a Taylor‐Green vortex. Numerical comparisons show that the CIRBF method outperforms some other methods in the literature. The CIRBF method in conjunction with a rational function transformation method and an adaptive time‐stepping scheme is then applied to simulate 1–D and 2–D soil infiltrations effectively. The proposed solutions are more accurate and converge faster than those of the finite different method used with a second‐order central difference scheme. Additionally, the present scheme also takes less time to achieve target accuracy in comparison with the 1D‐IRBF and higher order compact schemes.