平流方程有许多 FD 格式在网上讨论。例如这里: http: //farside.ph.utexas.edu/teaching/329/lectures/node89.html
但我还没有看到有人提出这样的“隐式”逆风方案: .
我见过的所有迎风方案都是处理空间导数中前一个时间步长的数据。这是什么原因?经典的迎风方案与我上面写的方案相比如何?
平流方程有许多 FD 格式在网上讨论。例如这里: http: //farside.ph.utexas.edu/teaching/329/lectures/node89.html
但我还没有看到有人提出这样的“隐式”逆风方案: .
我见过的所有迎风方案都是处理空间导数中前一个时间步长的数据。这是什么原因?经典的迎风方案与我上面写的方案相比如何?
在计算流体动力学中使用类似于您提出的隐式方案是很常见的。我所知道的基于紧凑的有限差分公式(不仅仅是替换和在现有计划中)。例如,Lele 于 1992 年在本文中开发了一种使用最广泛的方案,引用次数超过 2500 次。可以使此类方案具有比典型的显式方案更好的色散特性。
当使用隐式方法和大时间步长时,逆风通常不太重要,因为大量的扩散(杰里米提到)意味着你无论如何都无法解决冲击。
关于您提出的特定方案:
你没有理由不能做你所写的。这种情况不常见的原因之一是对于双曲(平流)类型的问题,依赖域是有限的。因此,从计算效率的角度来看,显式方法是有意义的。
您编写的隐式方案将需要求解一个线性系统,尽管在您编写三角形的情况下,因此求解相当简单。当然,当您研究系统和多个维度时,系统可能不会是三角形的,尽管有时这可能会导致未知数的正确排序(例如,参见Kwok 和 Tchelepi,JCP 2007和 Gustafsson 和 Khalighi,JSC,2006)。
有时,人们希望采用较大的时间步长,正如您所写的那样,使用隐式时间步长,但在这里您必须小心。使用隐式方法时,您将引入大量扩散,因此您将显着涂抹您的解决方案。