隐式方法(曲柄尼科尔森)我不明白程序

计算科学 pde 正则 双曲-pde 曲柄尼科尔森
2021-12-19 17:09:41

我试图理解通过这个方程的段落可以用fortran算法轻松解决,特别是我不明白L_x和L_xx的含义......什么(-1,0,1)代表?有人可以解释一下吗?多谢 !

这里报告了方程(显式方案很清楚.. 显然) 在此处输入链接描述

1个回答

Lx是表示更扩展符号的简写(运算符):因此可以写在节点中作为向量(参见的系数,即零和 1)。的结果相同Lxx

Lxui=(ui+1ui1)
Lxxui=(ui12ui+ui+1)
LxiLx=(1,0,1)ui11uiui+1Lxx

为了更好地理解,对于给定的模板,可以正式编写:标量积: u~i=(ui1,ui,ui+1)T

Lxu~i=(1,0,1)(ui1,ui,ui+1)T=ui+1ui1

希望这可以帮助。有关 CN 方案过程的描述,您可以查看:https ://scicomp.stackexchange.com/a/27190/22590