论文部分内容阅读
提出一种Fourier—Legendre谱元方法用于求解极坐标系下的Navier—Stokes方程，其中极点所在单元的径向采用Gauss—Radau积分点，避免了r=0处的1／r坐标奇异性。时间离散采用时间分裂法，引入数值同位素模型跟踪同位素的输运过程验证数值模拟的精度，分别利用谱元法和有限差分法的迎风差分格式求解匀速和加速坩埚旋转流动中的同位素方程。计算结果表明，有限差分法中的一阶迎风差分格式存在严重的数值假扩散，二阶迎风差分格式的数值结果较精确，增加节点可以有效地缓解数值扩散。然而，谱元法具有以较少节