Stabilized finite element method for the radial Dirac equation

被引:8
作者
Almanasreh, Hasan [1 ]
Salomonson, Sten
Svanstedt, Nils
机构
[1] Univ Gothenburg, Dept Math Sci, SE-41296 Gothenburg, Sweden
关键词
Dirac operator; Finite element scheme; Spurious eigenvalue; Cubic Hermite functions; Petrov-Galerkin; Stability parameter; SPURIOUS SOLUTIONS; ORIGIN;
D O I
10.1016/j.jcp.2012.11.020
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
A challenging difficulty in solving the radial Dirac eigenvalue problem numerically is the presence of spurious (unphysical) eigenvalues, among the genuine ones, that are neither related to mathematical interpretations nor to physical explanations. Many attempts have been made and several numerical methods have been applied to solve the problem using the finite element method (FEM), the finite difference method, or other numerical schemes. Unfortunately most of these attempts failed to overcome the difficulty. As a FEM approach, this work can be regarded as a first promising scheme to solve the spuriosity problem completely. Our approach is based on an appropriate choice of trial and test function spaces. We develop a Streamline Upwind Petrov-Galerkin method to the equation and derive an explicit stability parameter. (c) 2012 Elsevier Inc. All rights reserved.
引用
收藏
页码:426 / 442
页数:17
相关论文
共 23 条