leap-frog有限差分格式

penglining

方程如下，关键是如何用leap-frog有限差分格式实现第2式中的B×p项，谢谢

[ 本帖最后由 penglining 于 2006-11-20 07:47 编辑 ]

oasis_luo

通常leap-frog只对一阶的微分项而言。对于非微分项，我没有看到例子。不过
如果利用leap-frog的思想不知对于p(k)是否可以做以下近似。
p(k)=0.5*[p(k+1)+p(k-1)].

foolboy007

leap frog的非微分相就是函数值啊
就是p(m,n)
对第二式而言
p(m+1,n)-p(m-1,n)+cfl*(u(m,n+1)-u(m,n-1))+2*k*B*p(m,n)=0   % cfl is courant number; k is the time step