无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

无网格法耦合RNGk-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

 基于移动最小二乘无网格方法,耦合RNG (Re-Normalisation Group) k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM (Advection Upstream Splitting Method)⁺-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP (Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNG k-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

低雷诺数下翼型分离流动抽吸控制特性研究

为了系统研究抽吸系数和抽吸方向对抽吸效果的影响以及抽吸效益与抽吸能耗之间的关系,以NACA0012翼型表面分离流动为基准状态,在其吸力面设计了局部多孔抽吸结构,采用Roe格式和双时间步隐式算法(LUSGS),从抽吸系数、抽吸方向和抽吸能耗等方面,数值研究了低雷诺数下多孔分布式抽吸结构对流动分离的控制效果,通过边界层速度线型的变化分析了抽吸控制机理。研究结果表明:在翼型吸力面流动分离点附近一定区域内进行抽吸,可有效抑制流动分离,改善翼型气动性能;随着抽吸系数的增加,升阻比先是快速增长然后缓慢下降,且升阻比最大值提高了约1.3倍。抽吸控制能量消耗评估显示抽吸系数在合理范围内时,控制能耗明显小于控制效益。抽吸角度对抽吸控制有显著影响,当抽吸角度较大时,不仅翼型升阻比获得了提升,而且抽吸控制所消耗的能量也会进一步减少。这些结果有助于进一步为流动控制设计提供新的思路和方法。     

一种新的非结构动网格的生成方法

弹簧方法是一种常被采用的非结构动网格生成方法,文中采用一种新的非结构动网格生成方法即弹性体方法(Elastic solid method,ESM)。和弹簧方法相比,虽然弹性体方法的计算量大,但是可以处理大位移的情况,尤其善于处理网格扭转很厉害的情况。另外,采用弹性体方法可以减小网格重新生成的频率,从而可以减小映射误差,同时也可以减小网格重新生成带来的计算量。将两种方法的数值模拟结果进行对比,进一步表明该弹性体方法在处理动网格问题中的优势。     

两种k-ω型湍流模型在无网格方法中的应用研究

基于在高度各向异性点云结构中取得良好结果的点云重构技术,采用移动最小二乘无网格法解耦计算了k-ω型湍流模型方程和雷诺平均N-S方程(RANS),实现了对湍流流动的数值模拟。通过对NACA0012翼型、RAE2822翼型的绕流以及后向台阶大分离流动的数值模拟,比较分析了k-ωSST(Shear Stress Transport)和k-ωTNT(Turbulent/Non-Turbulent)两种模型的湍流预测能力。结果表明:对弱逆压梯度情况下的湍流,两者具有相当的预测能力;而k-ωSST较k-ωTNT具有更好的激波捕捉能力。对于大分离流动,两种模型较实验数据均存在一定差异,但总体趋势吻合良好。     

抽吸控制对低雷诺数下翼型分离流动的影响简

 为了系统研究抽吸区域相对于分离点的作用位置、孔间距和孔径对抽吸控制翼型分离流动效果的影响,以NACA0012翼型表面分离流动为基准状态,在其吸力面设计了局部多孔抽吸结构,采用AUSM⁺-up格式、大涡模拟方法和双时间步长（LU-SGS）隐式算法,对低雷诺数下多孔分布式抽吸结构对流动分离的控制效果进行了数值研究。研究结果表明：当抽吸区域位于分离点之后时,抽吸控制效果最好;抽吸系数不仅存在一个下限值以达到快速、有效的控制效果,而且有一个上限值以保证抽吸控制品质因数（FOM）大于1;孔间距和孔径对翼型气动性能的影响较小,但对FOM分布的影响较大。     

HLLC方法在三维非结构网格上的应用

以三维非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的二阶精度Godonov扩展方法求解三维Euler方程，使用HLLC(Harten，Lax，van Leer,Contact)近似黎曼解的方法计算网格单元边界处的守恒量通量。为了验证方法的可行性，利用三维解算器对NACA0012翼型绕流进行了数值计算，结果令人满意。在此基础上对某型弹的流场进行了数值模拟，取得了比较好的结果。     

抽吸控制对低雷诺数下翼型分离流动的影响

为了系统研究抽吸区域相对于分离点的作用位置、孔间距和孔径对抽吸控制翼型分离流动效果的影响,以NACA0012翼型表面分离流动为基准状态,在其吸力面设计了局部多孔抽吸结构,采用AUSM⁺-up格式、大涡模拟方法和双时间步长（LU-SGS）隐式算法,对低雷诺数下多孔分布式抽吸结构对流动分离的控制效果进行了数值研究。研究结果表明：当抽吸区域位于分离点之后时,抽吸控制效果最好;抽吸系数不仅存在一个下限值以达到快速、有效的控制效果,而且有一个上限值以保证抽吸控制品质因数（FOM）大于1;孔间距和孔径对翼型气动性能的影响较小,但对FOM分布的影响较大。