间断陡削化方程Ⅰ：人为刚性源项方法：刚性常微分方程的计算方法 …

本文针对界面运动的计算提出了适当附加人为刚性源项计算接触间断的方法,数值实验表明这一方法是有效的。     

带源项的守恒系统时间分裂算法定常解的收敛性研究

时间分裂算法很难获得带源项的守恒系统的定常解,本文中,考虑使用基本的隐式方法积分线性、二次和三次源项,研究了时间分裂算法不能收敛到数值定常态的原因.采用一个刚性参数的方式预测定常态的误差,刚性参数依赖于特定的源项.当源项是解的线性函数时,不存在定常态误差;当源项是解的非线形函数时,发现定常状态的误差是刚性参数的单调增加函数.定常态误差的分析将推广到标准k-epsilon双方程湍流模型计算的情况.     

湍流模型源项的数值处理方法对计算结果的影响研究

涡粘性湍流模型广泛应用于各种工程湍流问题的计算。对于大多数湍流模型,在湍流控制方程的右端会出现源项,其刚性给数值计算带来很大的影响。从源项的物理意义出发分析了源项导致求解困难的原因,并在通用的求解湍流模型方程的点隐法基础上,以源项弱刚性的S-A湍流模型和源项强刚性的k-ωSST湍流模型为例,论述了生成项和耗散项之间的平衡关系对数值模拟的重要意义,给出了针对不同源项的具体处理方法。对RAE2882翼型跨声速流动算例的模拟结果表明,这些处理方法有效的提高了控制方程组的计算稳定性。     

混合解析/数值方法在含高频振荡点源项的欧拉方程数值计算上的应用

 通过模型方程分析看到了一个重要的现象, 如果源项涉及的时间尺度远小于对流项的时间尺度, 那么基于对流时间尺度作为步长的传统数值方法, 即使源项相对于墩流项是个小量, 也会导致平均尺度上错误的结果。为了克服这种困难, 采用时间分裂方法, 把方程分裂成含对流项部分的偏微分方程( PDE) 和包含源项的常微分方程(ODE) 基础上, PDE 使用传统的数值方法, ODE 用解析的方法求解。该混合方法在数值格式时间步长小于平均流动时间尺度时, 得到正确解, 而与点源所隐含的时间尺度无关。把这个方法应用在含源流动的Euler 方程的计算中, 计算了翼型振荡问题, 取得了理想的结果。