期刊
APPLIED MATHEMATICS LETTERS
卷 26, 期 4, 页码 445-451出版社
PERGAMON-ELSEVIER SCIENCE LTD
DOI: 10.1016/j.aml.2012.11.006
关键词
Helmholtz-Hodge decompositions; Rotational penalty-projection; Vector penalty-projection; Penalty method; Error analysis; PDEs with adapted right-hand sides
We present new fast discrete Helmholtz-Hodge decomposition (DHHD) methods for efficiently computing at the order O(epsilon) the divergence-free (solenoidal) or curl-free (irrotational) components and their associated potentials for a given L-2 (Omega) vector field in a bounded domain. The solution algorithms solve suitable penalized boundary-value elliptic problems involving either the grad(div) operator in the vector penalty-projection (VPP) or the rot(rot) operator in the rotational penalty-projection (RPP) with adapted right-hand sides of the same form. Therefore, they are extremely well-conditioned, fast and cheap, avoiding having to solve the usual Poisson problems for the scalar or vector potentials. Indeed, each (VPP) or (RPP) problem only requires two conjugate-gradient iterations whatever the mesh size, when the penalty parameter epsilon is sufficiently small. We state optimal error estimates vanishing as O(epsilon) with a penalty parameters as small as desired up to machine precision, e.g. epsilon = 10(-14). Some numerical results confirm the efficiency of the proposed (DHHD) methods, very useful for solving problems in electromagnetism or fluid dynamics. (C) 2012 Elsevier Ltd. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据