Divergence-free meshless local Petrov-Galerkin method for Stokes flow

被引:9
作者
Najafi, Mahboubeh [1 ]
Dehghan, Mehdi [1 ]
Sarler, Bozidar [2 ,3 ]
Kosec, Gregor [4 ]
Mavric, Bostjan [2 ,3 ]
机构
[1] Amirkabir Univ Technol, Dept Math & Comp Sci, Tehran Polytech, Tehran, Iran
[2] Univ Ljubljana, Fac Mech Engn, Ljubljana, Slovenia
[3] Inst Met & Technol, Ljubljana, Slovenia
[4] Jozef Stefan Inst, Parallel & Distributed Syst Lab, Ljubljana, Slovenia
关键词
Incompressible flows; Stokes equations; Divergence-free radial kernels; MLPG method; Weak formulation; INCOMPRESSIBLE VISCOUS FLOWS; BOUNDARY ELEMENTS METHOD; RADIAL BASIS FUNCTIONS; NUMERICAL-SOLUTION; MLPG METHOD; ARTIFICIAL COMPRESSIBILITY; DMLPG METHOD; NATURAL-CONVECTION; EQUATIONS; SYSTEM;
D O I
10.1007/s00366-022-01621-w
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
The purpose of the present paper is development of an efficient meshless solution of steady incompressible Stokes flow problems with constant viscosity in two dimensions, with algebraic order of accuracy. This is achieved by employing a weak formulation with divergence-free matrix-valued quadratic Matern (QM) radial basis function (RBF) for the shape function and divergence-free matrix-valued compactly supported Gaussian (CSG) RBF for the weight function on the computational domain and its boundary. The continuity equation is inherently built-in in the formulation and the pressure is eliminated from the formulation with the aid of divergence theorem and the choice of divergence-free weight function. The developed method is thus iteration free, and results in a banded system of equations to be solved jointly for both velocity components. Gauss-Legendre cell integration is performed in the current investigation. The characteristics of the method are assessed by changing its free parameters, i.e., weight functions' sub-domain radius and shape functions' support domain radius and the shape parameter. A sensitivity test for several choices of shape functions with regular centers arrangement is done to identify the appropriate support size for the shape and weight functions and stagnation errors are reported accordingly. To the best of our knowledge, this article is initiative in introducing the application of divergence-free MLPG method to incompressible flows, aiming at elimination of pressure from the governing equations in primitive variables, with the aid of divergence-free RBFs through weak formulation. Only the momentum equation needs to be solved. Hence, the formulation of the problem is much simpler than the building of divergence-free elements in the related mesh-based methods.
引用
收藏
页码:5359 / 5377
页数:19
相关论文
共 73 条
[1]   Direct meshless local Petrov-Galerkin (DMLPG) method for time-fractional fourth-order reaction-diffusion problem on complex domains [J].
Abbaszadeh, Mostafa ;
Dehghan, Mehdi .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2020, 79 (03) :876-888
[2]   Investigation of the Oldroyd model as a generalized incompressible Navier-Stokes equation via the interpolating stabilized element free Galerkin technique [J].
Abbaszadeh, Mostafa ;
Dehghan, Mehdi .
APPLIED NUMERICAL MATHEMATICS, 2020, 150 :274-294
[3]  
Atluri S.N., 2002, MESHLESS LOCAL PETRO
[4]   FINITE-ELEMENT FOR THE NUMERICAL-SOLUTION OF VISCOUS INCOMPRESSIBLE FLOWS [J].
BERCOVIER, M ;
ENGELMAN, M .
JOURNAL OF COMPUTATIONAL PHYSICS, 1979, 30 (02) :181-201
[5]   Rayleigh-Taylor and Kelvin-Helmholtz instability studied in the frame of a dimension-reduced model [J].
Bestehorn, Michael .
PHILOSOPHICAL TRANSACTIONS OF THE ROYAL SOCIETY A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES, 2020, 378 (2174)
[6]  
Brenner SC., 2020, 75 YEARS MATH COMPUT
[7]  
Carrero J, 2006, MATH COMPUT, V75, P533, DOI 10.1090/S0025-5718-05-01804-1
[8]   Effect of permeability anisotropy on buoyancy-driven flow for CO2 sequestration in saline aquifers [J].
Cheng, Philip ;
Bestehorn, Michael ;
Firoozabadi, Abbas .
WATER RESOURCES RESEARCH, 2012, 48
[9]  
Cockburn B, 2005, MATH COMPUT, V74, P1067, DOI 10.1090/S0025-5718-04-01718-1
[10]   Local discontinuous Galerkin methods for the Stokes system [J].
Cockburn, B ;
Kanschat, G ;
Schotzau, D ;
Schwab, C .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2002, 40 (01) :319-343