使用 Lax-Friedrichs 方案求解输运方程时,应在哪些时空点评估已知系数函数?

计算科学 pde 有限差分 数字 双曲-pde
2021-12-17 02:35:08

对于标量,我正在考虑传输方程 为常数时,Lax-Friedrichs 方案是由 其中用于等间距的时空网格点。u=u(x,t)

ut+aux=0,x[0,L], t[0,T]u(0,t)=uin,u(x,0)=f(x).
a
ujn+1=uj+1n+uj1n2a(uj+1nuj1n)2,
ujn=u(xj,tn)

的情况下,应该评估鉴于该方案在时间上是向前的,我可以合理地说服自己我们将时间上的 a 评估a 但是我们应该使用哪个空间点?Lax-Friedrichs 方案的隐含空间特性使我无法使用与时间相同的逻辑。a=a(x,t)xj, tnaa

a(x,tn+1),

2个回答

我建议将此作为规范概括

ut(xj,tn)ujn+1(uj1n+uj+1n)/2Δt

(aux)(xj,tn)a(xj,tn)uj+1nuj1n2Δx

ujn+1=12(uj1n+uj+1n)ajnΔt2Δx(uj+1nuj1n)

请参阅https://en.wikipedia.org/wiki/Lax–Friedrichs_method,其中描述了非线性守恒定律的扩展。您可能还想使用谷歌查看大量使用 LF 方法的非线性守恒定律的文章。