A Fast Iterative Algorithm for Elliptic Interface Problems

数学 舒尔补语 离散化 分段 泊松方程 广义最小残差法 数学分析 应用数学 线性系统 量子力学 物理 特征向量
作者
Zhilin Li
出处
期刊:SIAM Journal on Numerical Analysis [Society for Industrial and Applied Mathematics]
卷期号:35 (1): 230-254 被引量:247
标识
DOI:10.1137/s0036142995291329
摘要

A fast, second-order accurate iterative method is proposed for the elliptic equation \[ \grad\cdot(\beta(x,y) \grad u) =f(x,y) \] in a rectangularregion $\Omega$ in two-space dimensions. We assume that there is an irregular interface across which the coefficient $\beta$, the solution u and its derivatives, and/or the source term f may have jumps. We are especially interested in the cases where the coefficients $\beta$ are piecewise constant and the jump in $\beta$ is large. The interface may or may not align with an underlying Cartesian grid. The idea in our approach is to precondition the differential equation before applying the immersed interface method proposed by LeVeque and Li [ SIAM J. Numer. Anal., 4 (1994), pp. 1019--1044]. In order to take advantage of fast Poisson solvers on a rectangular region, an intermediate unknown function, the jump in the normal derivative across the interface, is introduced. Our discretization is equivalent to using a second-order difference scheme for a corresponding Poisson equation in the region, and a second-order discretization for a Neumann-like interface condition. Thus second-order accuracy is guaranteed. A GMRES iteration is employed to solve the Schur complement system derived from the discretization. A new weighted least squares method is also proposed to approximate interface quantities from a grid function. Numerical experiments are provided and analyzed. The number of iterations in solving the Schur complement system appears to be independent of both the jump in the coefficient and the mesh size.

科研通智能强力驱动
Strongly Powered by AbleSci AI
更新
PDF的下载单位、IP信息已删除 (2025-6-4)

科研通是完全免费的文献互助平台,具备全网最快的应助速度,最高的求助完成率。 对每一个文献求助,科研通都将尽心尽力,给求助人一个满意的交代。
实时播报
大橙子发布了新的文献求助10
2秒前
领导范儿应助科研通管家采纳,获得10
3秒前
量子星尘发布了新的文献求助10
7秒前
明钟达完成签到 ,获得积分10
15秒前
byyyy完成签到,获得积分10
18秒前
高高的哈密瓜完成签到 ,获得积分10
22秒前
Rondab应助橙汁采纳,获得10
25秒前
读书的时候完成签到,获得积分10
27秒前
颜云尔完成签到,获得积分10
38秒前
孤独雨梅完成签到,获得积分10
41秒前
woobinhua完成签到 ,获得积分10
41秒前
雪落你看不见完成签到,获得积分10
43秒前
十月天秤完成签到,获得积分0
44秒前
依文完成签到,获得积分20
44秒前
ymr完成签到 ,获得积分10
45秒前
哦哦哦完成签到 ,获得积分10
46秒前
jzmupyj完成签到,获得积分10
46秒前
大橙子发布了新的文献求助10
49秒前
xdlongchem完成签到,获得积分10
50秒前
量子星尘发布了新的文献求助10
52秒前
小梦完成签到,获得积分10
53秒前
xuhang完成签到,获得积分10
53秒前
ZSHAN完成签到,获得积分10
54秒前
美满的机器猫完成签到,获得积分10
57秒前
王小磊完成签到,获得积分10
1分钟前
谢花花完成签到 ,获得积分10
1分钟前
1分钟前
瓦罐完成签到 ,获得积分10
1分钟前
扁舟灬完成签到,获得积分10
1分钟前
Cpp完成签到 ,获得积分10
1分钟前
贤惠的老黑完成签到 ,获得积分10
1分钟前
ame1120发布了新的文献求助10
1分钟前
倦梦还完成签到,获得积分10
1分钟前
Sunrise完成签到,获得积分10
1分钟前
yyyy发布了新的文献求助10
1分钟前
自觉柠檬完成签到 ,获得积分10
1分钟前
ergatoid完成签到,获得积分10
1分钟前
Hao完成签到,获得积分10
1分钟前
月亮煮粥完成签到,获得积分10
1分钟前
欣欣完成签到 ,获得积分10
1分钟前
高分求助中
【提示信息,请勿应助】关于scihub 10000
Les Mantodea de Guyane: Insecta, Polyneoptera [The Mantids of French Guiana] 3000
徐淮辽南地区新元古代叠层石及生物地层 3000
The Mother of All Tableaux: Order, Equivalence, and Geometry in the Large-scale Structure of Optimality Theory 3000
Handbook of Industrial Diamonds.Vol2 1100
Global Eyelash Assessment scale (GEA) 1000
Picture Books with Same-sex Parented Families: Unintentional Censorship 550
热门求助领域 (近24小时)
化学 材料科学 医学 生物 工程类 有机化学 生物化学 物理 内科学 纳米技术 计算机科学 化学工程 复合材料 遗传学 基因 物理化学 催化作用 冶金 细胞生物学 免疫学
热门帖子
关注 科研通微信公众号,转发送积分 4038128
求助须知:如何正确求助?哪些是违规求助? 3575831
关于积分的说明 11373827
捐赠科研通 3305610
什么是DOI,文献DOI怎么找? 1819255
邀请新用户注册赠送积分活动 892655
科研通“疑难数据库(出版商)”最低求助积分说明 815022