第五讲 典型模型方程波动方程的差分格式.pdfVIP

波动方程的几种差分格式 模型方程：波动方程（wave equation） 该方程是一个线性方程，描述了以波速c沿x轴传播，可 以做为了解非线性无粘性流动的初步性质的模型方程。 不可压缩Navier–Stokes 方程   线性化，压强梯度略去 u v  0   x y 2 2 一维无粘性   1        u u p u u u 1 p  u u u  u v   ( 2  2 ) t x  x       对流项 t x y  x  x y 2 2    1    v  v  v  p  v  v 扩散项 u v ( 2 2 )       t x y  y  x y 仅含对流项，故该方程也可以称为一维对流方程 波动方程（wave equation） 1-D 波动方程是一个二阶双曲型PDE: （1） 表示声波以波速c在单一物质中转播的控制方程。同 样性质的方程还有一阶PDE ： （2 ） 式（1）可由（2 ）推导出 式（2 ）可得： ut cux ，分别对此式求t和x 的导数 u cu u cu tt xt tx xx 2 得：utt cuxt c(cuxx ) c uxx 即式（1） 初值问题的精确解 由数学物理方程得： （3 ） （4 ） 由（4 ）分别求t 、x导数，代入（2 ）式可以验证解的正确性， 即： cF cF 0 欧拉显格式:Euler Explicit Method 都为一步显格式 1、 ： 时间、空间前差