科研之友
微信
新浪微博
Facebook
分享链接
摘要
对于具有多孔介质底面的轴对称二相重力流 ,引进基于浅水近似的控制方程和相应的边界条件 ,采用贴体坐标变换使运动边界问题化为固定边界问题 ,提出了基于特征插值并结合使用梯形积分公式和Newton Raphson迭代法在时间和空间都具有二阶精度的数值边界条件 .为检验格式的性能和避免编写程序时可能出现的错误 ,对类似的方程构造了一类精确解 .在空间上采用了二步Lax格式、二阶TVD格式、三阶ENO格式及五阶WENO格式 ,在时间上采用了二阶及三阶的TVD Runge Kutta方法对该问题进行数值模拟 .数值结果表明 ,在解的光滑区域 ,这几种格式的精度都很高 ,但是在大梯度区 ,二步Lax格式将会产生强烈的数值振荡 ,且振荡不会随网格宽度的减小而减小 ,而其他 3种格式将不会或仅会产生幅度要小得多的数值振荡 ,且振荡会随网格宽度的减小而趋向于零 .对实际应用目的来说 ,结合使用二阶TVD Runge Kutta方法的二阶TVD格式是一个经济而又适当的选择
- 单位
