三维抛物型方程两层显式差分格式的研究

近年来 ,对三维抛物型方程的数值解法的研究逐渐增多 ,出现了一些粗度高、稳定性好的差分格式。但它们或是三层隐式格式[1] ,或是三层显式格式[2 ,3 ,4 ] ,隐格式常因计算量和存储量太大而难以使用 ;三层显格式虽是显式计算 ,但却不能计算第一层上的网格函数值 ,需用其他方法先行启动 ,在实际使用上是不方便的 ,因此 ,构造精度高、稳定性好的两层显式差分格式便具有十分明显的理论意义和实用价值。作者对求解区域 D:{0≤ x,y,z≤ L ,0≤ t≤ T}上的三维抛物型方程初边值问题 u t= 2 u x2 + 2 y y2 + 2 y z2u| t=0 =φ( x,y,z) ,u| x=0 =…