非稳态流动中球形颗粒阻力系数的数值研究.docVIP

第八届全国流体力学学术会议 2014年9月18～21日 甘肃兰州 文章编号： CSTAM2014-B01-0074 标题：非稳态流动中球形颗粒阻力系数的数值研究 作者：林立，张凡，屈秋林，刘沛清 单位：北京航空航天大学航空科学与工程学院 Copyright ? 2014 版权所有 中国力学学会 地址: 北京市北四环西路15号　邮政编码:100190 　Address: No.15 Beisihuanxi Road, Beijing 100190 非稳态流动中球形颗粒阻力系数的数值研究1） 林立*，张凡*，屈秋林*，刘沛清*，2） *（北京航空航天大学航空科学与工程学院，北京，100191） 摘要　本文使用ANSYS FLUENT商用软件对球形颗粒在非稳态流场中的阻力系数进行数值研究，雷诺数范围是 ，计算过程中将部分动网格与整体动网格相结合，以保证计算的准确性。研究结果表明，在给定的雷诺数范围内，颗粒的阻力系数与加速度数密切相关。根据数值计算的结果，对于颗粒加速和减速的不同运动条件，本文分别给出阻力系数关于加速度数的修正公式。对颗粒表面压力分布的分析表明，颗粒前缘和后缘压力的变化是导致颗粒阻力系数改变的主要原因。最后，本文将修正后的阻力系数公式带入翼型扰流的计算中，对该修正公式在实际两相流计算中的适用情况进行了模拟和验证。 关键词　两相流计算，非稳态流动，颗粒阻力系数，加速度数 引 言 在涉及球形颗粒的两相流计算中，阻力系数公式的精度对颗粒轨迹计算的准确性影响很大。球形颗粒的阻力系数通常是在稳态流场中测量所得，只与雷诺数相关。如果物理过程非稳态，则按照稳态公式计算的结果与真实值往往有很大的偏差。因此，对球形颗粒在非稳态运动时的阻力系数进行修正，对于两相流计算精度的提高意义重大。 前人对于圆球在流体中进行非稳态运动时阻力的研究，以实验居多，理论研究较少。这一课题的早期工作几乎都与流体中的摆动运动有关，目的是发现摆动的正确周期。 Stokes1通过理论推导，给出了在流体中做线性小振幅摆动的圆球的阻力系数。对于无粘流体，阻力可以简化为Poisson理论中由附加质量所产生的阻力。在此基础上，Boussinesq和Basset2指出，对于圆球在流体中的非定常缓慢运动，小球的受力可以写成三项之和：稳态运动阻力，以及由附加质量和运动历史所引起的阻力。把简谐运动方程代入其公式中，就可以得到Stokes所推导的阻力系数。 OdarHamilton3通过为上述公式中的每一项引入经验系数，修正了Boussinesq-Basset表达式。他们测量了圆球在粘性流体中以有限振幅进行摆动时的附加质量项和历史运动项系数。实验发现，这两个系数都只是无量纲参数 的函数。 McKownKeulegan4进一步通过实验指出，附加质量项可以被从物理上解释成受到随圆球加速的尾流的影响。而Schoneborn5则认为在摆动流动中，当流场频率处于涡脱落的固有频率区间时，OdarHamilton通过经验公式预测的阻力过小，圆球周围的流动以及圆球所受到的阻力与流场中的波动频率密切相关。 IversonBalent6和Keim7首先通过量纲分析提出了加速度数的概念。加速度数是当地加速度产生的阻力与对流加速度产生的阻力项之比。当圆球进行非稳态运动时，当地加速度不为0，因此加速度数不为0，当地加速度会对圆球周围的流场和压力分布产生直接的影响，进而使圆球的阻力系数偏离稳态值。这种影响可以使用加速度数进行描述。 KaranfilianKotas8的研究结果表明，当颗粒进行非稳态运动时，颗粒的阻力系数不但与雷诺数有关，而且与加速度数密切相关，并且在加速度数相同时，加速和减速时的阻力系数呈现出显著的不同。根据数值计算的结果，可以给出球形颗粒非稳态运动的阻力系数关于加速度数的修正公式。 有鉴于此，本文采用CFD方法计算圆球在非稳态运动中的阻力系数，与稳态阻力系数进行对比并加以修正。最后把修正后的阻力系数公式带入异性绕流的计算中，以验证修正公式的效果及适用范围。 1 物理模型与数值方法 本文使用ANSYS FLUENT软件对圆球非稳态运动时的阻力系数进行数值研究，圆球直径D=1mm，非稳态运动的圆球加速度分别为32,40,50,64,80,100,125m/s2，在第一组计算中，圆球按上述加速度从静止开始匀加速到2.5m/s，在第二组计算中，圆球由2.5m/s匀减速到静止，雷诺数范围为。 本文采用整体动网格与部分动网格进行对比计算，以保证结果的准确性。 部分动网格，即指在圆球运动的过程中，整体网格的位置保持不变。