Local meshless methods for second order elliptic interface problems with sharp corners

被引:24
作者
Ahmad, Masood [1 ]
Siraj-ul-Islam [1 ]
Larsson, Elisabeth [2 ]
机构
[1] Univ Engn & Technol, Fac Architecture Allied Sci & Humanities, Dept Basic Sci, Peshawar, Pakistan
[2] Uppsala Univ, Dept Informat Technol, Sci Comp, Uppsala, Sweden
关键词
RADIAL BASIS FUNCTIONS; DIFFERENCES RBF-FD; STABLE COMPUTATION; MATCHED INTERFACE; ALGORITHM; EQUATIONS; COLLOCATION; PARAMETER; ELEMENTS; EDGES;
D O I
10.1016/j.jcp.2020.109500
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In the present paper, we develop a local meshless procedure for solving a steady state two-dimensional interface problem having discontinuous coefficients and curved interfaces with sharp corners. The proposed local meshless methods are based on three types of radial basis functions (RBFs): a local meshless method based on multiquadric RBF (LMM1P), a local meshless method based on integrated multiquadric RBF (LMM2P) and a local meshless method based on hybrid Gaussian-Cubic RBF (LMM3P). Stencils are designed at the interface and interior regions to cope with discontinuities and sharp corners. Due to the localized nature of the procedure and a sparse matrix representation, the local meshless methods become computationally less expensive than global meshless methods. The methods are augmented with linear polynomial to improve accuracy and ensure stable computation. Comparison with some existing versions of finite element methods is also performed to show better accuracy of the proposed meshless methods. Accuracies of the proposed local meshless methods are also compared among themselves. Flexibility of the meshless methods with respect to complex geometries, adapting to different shapes of the interfaces and selection of the shape parameter is also considered. © 2020 Elsevier Inc.
引用
收藏
页数:17
相关论文
共 47 条
[1]   Meshless analysis of parabolic interface problems [J].
Ahmad, Masood ;
Siraj-ul-Islam .
ENGINEERING ANALYSIS WITH BOUNDARY ELEMENTS, 2018, 94 :134-152
[2]   High-order numerical schemes based on difference potentials for 2D elliptic problems with material interfaces [J].
Albright, Jason ;
Epshteyn, Yekaterina ;
Medvinsky, Michael ;
Xia, Qing .
APPLIED NUMERICAL MATHEMATICS, 2017, 111 :64-91
[3]  
Appelö D, 2009, COMMUN COMPUT PHYS, V5, P84
[4]   Is the pollution effect of the FEM avoidable for the Helmholtz equation considering high wave numbers? [J].
Babuska, IM ;
Sauter, SA .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1997, 34 (06) :2392-2423
[5]   NUMERICAL-SIMULATION OF TURBULENT-FLOW OVER SURFACE-MOUNTED OBSTACLES WITH SHARP EDGES AND CORNERS [J].
BAETKE, F ;
WERNER, H ;
WENGLE, H .
JOURNAL OF WIND ENGINEERING AND INDUSTRIAL AERODYNAMICS, 1990, 35 (1-3) :129-147
[6]   A second order virtual node method for elliptic problems with interfaces and irregular domains [J].
Bedrossian, Jacob ;
von Brecht, James H. ;
Zhu, Siwei ;
Sifakis, Eftychios ;
Teran, Joseph M. .
JOURNAL OF COMPUTATIONAL PHYSICS, 2010, 229 (18) :6405-6426
[7]   MAGNETIC-FIELD COMPUTATION USING DELAUNAY TRIANGULATION AND COMPLEMENTARY FINITE-ELEMENT METHODS [J].
CENDES, ZJ ;
SHENTON, D ;
SHAHNASSER, H .
IEEE TRANSACTIONS ON MAGNETICS, 1983, 19 (06) :2551-2554
[8]   Subdomain radial basis collocation method for heterogeneous media [J].
Chen, Jiun-Shyan ;
Wang, Lihua ;
Hu, Hsin-Yun ;
Chi, Sheng-Wei .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2009, 80 (02) :163-190
[9]  
Davydov O, 2001, COMPUT MATH APPL, V62, P27
[10]   A new approach based on the genetic algorithm for finding a good shape parameter in solving partial differential equations by Kansa's method [J].
Esmaeilbeigi, M. ;
Hosseini, M. M. .
APPLIED MATHEMATICS AND COMPUTATION, 2014, 249 :419-428