Skip to main content
SpringerLink
Log in

Space–time polyharmonic radial polynomial basis functions for modeling saturated and unsaturated flows

  • Original Article
  • Published:
Engineering with Computers Aims and scope Submit manuscript

Abstract

In this paper, we propose a novel meshless approach that involves using space–time polyharmonic radial polynomial basis functions for modeling saturated and unsaturated flows in porous media. In this study, space–time polyharmonic radial polynomial basis functions were developed in the space–time domain using a meshless collocation method. This domain contains three sets of collocation points, namely the inner, source, and boundary points, for the spatial and temporal discretization of the governing equation. Because the initial and boundary data are accessible space–time boundaries, the solutions of groundwater flows problems are approximated by solving the inverse boundary value problem in the space–time domain without using the conventional time-marching scheme. Saturated and unsaturated flow problems were investigated to demonstrate the robustness of the proposed method. The results obtained using the proposed approach were compared with those obtained using the conventional polyharmonic spline radial basis function. The proposed space–time polyharmonic radial polynomial basis functions obtained highly accurate solutions. Moreover, in solving saturated and unsaturated flow problems, the accuracy and stability of the proposed functions were higher than those of the conventional time-marching scheme.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Fig. 1

source collocation scheme: a x–y plane projection of the space–time domain and b the space–time domain

Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7

source collocation schemes for three shapes: a shape A, b shape B, and c shape C

Fig. 8
Fig. 9

source collocation scheme

Fig. 10
Fig. 11
Fig. 12

Similar content being viewed by others

References

  1. Kacimov AR, Obnosov YV (2017) Analytical solution for tension-saturated and unsaturated flow from wicking porous pipes in subsurface irrigation: the Kornev-Philip legacies revisited. Water Resour Res 53(3):2542–2552

    Article  Google Scholar 

  2. Wu LZ, Zhu SR, Peng J (2020) Application of the Chebyshev spectral method to the simulation of groundwater flow and rainfall-induced landslides. Appl Math Model 80:408–425

    Article  MathSciNet  MATH  Google Scholar 

  3. Sakhaei Z, Nikooee E, Riazi M (2020) A new formulation for non-equilibrium capillarity effect using multi-gene genetic programming (MGGP): accounting for fluid and porous media properties. Eng Comput:1–13

  4. Liu CY, Ku CY, Huang CC, Lin DG, Yeih WC (2015) Numerical solutions for groundwater flow in unsaturated layered soil with extreme physical property contrasts. Int J Nonlinear Sci Numer Simul 16(7–8):325–335

    Article  MathSciNet  MATH  Google Scholar 

  5. Ku CY, Liu CY, Su Y, Xiao JE, Huang CC (2017) Transient modeling of regional rainfall-triggered shallow landslides. Environ Earth Sci 76(16):1–18

    Article  Google Scholar 

  6. Cai JS, Yeh TCJ, Yan EC, Tang RX, Hao YH, Huang SY, Wen JC (2019) Importance of variability in initial soil moisture and rainfalls on slope stability. J Hydrol 571:265–278

    Article  Google Scholar 

  7. Woodward CS, Dawson CN (2000) Analysis of expanded mixed finite element methods for a nonlinear parabolic equation modeling flow into variably saturated porous media. SIAM J Numer Anal 37(3):701–724

    Article  MathSciNet  MATH  Google Scholar 

  8. Liu CY, Ku CY, Xiao JE, Huang CC, Hsu SM (2017) Numerical modeling of unsaturated layered soil for rainfall-induced shallow landslides. J Environ Eng Landsc Manag 25(4):329–341

    Article  Google Scholar 

  9. Caviedes-Voullième D, Garcı P, Murillo J (2013) Verification, conservation, stability and efficiency of a finite volume method for the 1D Richards equation. J Hydrol 480:69–84

    Article  Google Scholar 

  10. Ku CY, Liu CY, Xiao JE, Yeih W (2017) Transient modeling of flow in unsaturated soils using a novel collocation meshless method. Water 9(12):954

    Article  Google Scholar 

  11. Li PW (2020) Space–time generalized finite difference nonlinear model for solving unsteady Burgers’ equations. Appl Math Lett 114:106896

    Article  MathSciNet  MATH  Google Scholar 

  12. Fu ZJ, Chen W, Ling L (2015) Method of approximate particular solutions for constant-and variable-order fractional diffusion models. Eng Anal Bound Elem 57:37–46

    Article  MathSciNet  MATH  Google Scholar 

  13. Grabski JK (2020) A meshless procedure for analysis of fluid flow and heat transfer in an internally finned square duct. Heat Mass Transf 56(2):639–649

    Article  Google Scholar 

  14. Li M, Chen CS, Hon YC (2011) A meshless method for solving nonhomogeneous Cauchy problems. Eng Anal Bound Elem 35(3):499–506

    Article  MathSciNet  MATH  Google Scholar 

  15. Chen CS, Karageorghis A, Dou F (2020) A novel RBF collocation method using fictitious centres. Appl Math Lett 101:106069

    Article  MathSciNet  MATH  Google Scholar 

  16. Mirzaee F, Samadyar N (2020) Combination of finite difference method and meshless method based on radial basis functions to solve fractional stochastic advection–diffusion equations. Eng Comput 36:1673–2168

    Article  Google Scholar 

  17. Li J, Cheng AHD, Chen CS (2003) A comparison of efficiency and error convergence of multiquadric collocation method and finite element method. Eng Anal Bound Elem 27(3):251–257

    Article  MATH  Google Scholar 

  18. Ku CY, Liu CY, Xiao JE, Hsu SM (2020) Multiquadrics without the Shape parameter for solving partial differential equations. Symmetry 12(11):1813

    Article  Google Scholar 

  19. Li J, Chen Y, Pepper D (2003) Radial basis function method for 1-D and 2-D groundwater contaminant transport modeling. Comput Mech 32(1):10–15

    Article  MATH  Google Scholar 

  20. Swathi B, Eldho TI (2014) Groundwater flow simulation in unconfined aquifers using meshless local Petrov-Galerkin method. Eng Anal Bound Elem 48:43–52

    Article  MathSciNet  MATH  Google Scholar 

  21. Uddin M (2014) On the selection of a good value of shape parameter in solving time-dependent partial differential equations using RBF approximation method. Appl Math Model 38(1):135–144

    Article  MathSciNet  MATH  Google Scholar 

  22. Fornberg B, Piret C (2008) On choosing a radial basis function and a shape parameter when solving a convective PDE on a sphere. J Comput Phys 227(5):2758–2780

    Article  MathSciNet  MATH  Google Scholar 

  23. Segeth K (2019) Polyharmonic splines generated by multivariate smooth interpolation. Comput Math Appl 78(9):3067–3076

    Article  MathSciNet  MATH  Google Scholar 

  24. Gardner WR (1958) Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table. Soil Sci 85:228–232

    Article  Google Scholar 

  25. Pham HQ, Fredlund DG, Barbour SL (2005) A study of hysteresis models for soil-water characteristic curves. Can Geotech J 42(6):1548–1568

    Article  Google Scholar 

  26. Fattah MY, Salim NM, Irshayyid EJ (2017) Determination of the soil–water characteristic curve of unsaturated bentonite–sand mixtures. Environ Earth Sci 76(5):201

    Article  Google Scholar 

  27. Ku CY, Liu CY, Xiao JE, Hsu SM, Yeih W (2021) A collocation method with space–time radial polynomials for inverse heat conduction problems. Eng Anal Bound Elem 122:117–131

    Article  MathSciNet  MATH  Google Scholar 

  28. Hamaidi M, Naji A, Charafi A (2016) Space–time localized radial basis function collocation method for solving parabolic and hyperbolic equations. Eng Anal Bound Elem 67:152–163

    Article  MathSciNet  MATH  Google Scholar 

  29. Nie WB, Li YB, Fei LJ, Ma XY (2017) Approximate explicit solution to the Green-Ampt infiltration model for estimating wetting front depth. Water 9(8):609

    Article  Google Scholar 

  30. Tracy FT (2006) Clean two- and three-dimensional analytical solutions of Richards’ equation for testing numerical solvers. Water Resour Res 42:1–11

    Article  Google Scholar 

Download references

Acknowledgements

This study was conducted with the financial support of the Ministry of Science and Technology, Taiwan, the Republic of China.

Author information

Authors and Affiliations

  1. School of Engineering, Department of Harbor and River Engineering, National Taiwan Ocean University, Keelung, 20224, Taiwan

    Cheng-Yu Ku, Li-Dan Hong, Chih-Yu Liu & Jing-En Xiao

Authors
  1. Cheng-Yu Ku
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Li-Dan Hong
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Chih-Yu Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jing-En Xiao
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Chih-Yu Liu.

Ethics declarations

Conflict of interest

The authors declare that they have no conflicts of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ku, CY., Hong, LD., Liu, CY. et al. Space–time polyharmonic radial polynomial basis functions for modeling saturated and unsaturated flows. Engineering with Computers 38, 4947–4960 (2022). https://doi.org/10.1007/s00366-021-01519-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00366-021-01519-z

Keywords

Search

Navigation