非结构网格下可压缩黏性流体的隐式算法

采用牛顿线化和Roe近似Jacobians矩阵隐式算法,实现基于非结构网格的三维定常黏性流动问题的求解.分别对平板、跨声涡轮叶片Mark-Ⅱ和CW-22这三种具有代表性的流动进行数值模拟,并详细分析了流动细节,模拟结果与理论以及实验值吻合良好,准确预测了流场内激波、二次流等重要的流动现象.验证了基于非结构网格的可压缩流体隐式计算方法对于黏性流动问题的计算是正确的.      

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.      

OA212翼型主动流动控制的数值模拟研究

采用数值模拟的方法,探讨了基于零质量射流的主动流动控制技术对OA212旋翼翼型动态失速的控制效果和控制特性.以积分形式雷诺平均Navier-Stokes(N-S)方程为控制方程,采用格心有限体积法进行求解.空间离散采用AUSM~+-up格式,时间推进采用含牛顿型LU-SGS子迭代的全隐式双时间法,且引入了预处理方法和多重网格方法加速收敛.通过在喷口上施加非定常边界条件来模拟射流对翼型绕流的影响.研究了不同类型射流、不同位置射流以及不同控制参数(频率、相位、偏角、动量系数等)对动态失速控制效果的影响.研究表     

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

预处理方法在叶轮机械三维数值模拟中的应用

应用预处理方法求解三维可压缩流动Navier-Stokes方程, 发展了一种适用于叶轮机械的既可以求解跨声速流动也可以求解不可压流动的数值模拟方法, 对NASA Rotor67和北航低速大尺寸压气机算例进行了计算, 计算结果与实验数据相吻合, 表明本文采用的预处理方法具有较大的工程实用性.      

基于N-S方程的民航机翼双目标优化设计

提出了一种可用于民航机翼的双目标优化设计方法.优化中的设计变量为给定机翼沿展向若干个剖面的厚度与扭角,优化目标为固定升力系数下的升阻比和机翼容积.采取N-S(Navier-Stokes)方程作为流场求解器,Powell方法作为优化求解方法,对Lockheed-AFSOR Wing A 某民航机翼进行了优化设计.在马赫数为0.6217条件下,Wing A 机翼双目标优化结果优于单目标优化结果.对某实际应用的民航机翼在巡航马赫数为0.78、升力系数为0.48的优化结果也表明,该方法能有效提高机翼升阻比和容积.相对初始机翼,双目标优化机翼最大相对厚度在翼根处增大,翼梢处减小,扭角绝对值相对初始机翼沿展向均减小.该双目标优化设计方法优于常规的单目标优化设计方法,可以为民航机翼的工程设计提供有益的参考.     

基于CFD方法对微型扑翼翼型设计的研究

基于CFD方法从扑翼设计的角度出发对扑翼的外形进行深入的研究.使用基于广义无限插值理论的代数网格生成方法开发适合扑翼运动的三维结构N-S动网格,并应用有限体积法结合双时间推进技术求解预处理后的低速三维非定常N-S方程,完成了对低速扑翼绕流的数值模拟.通过与文献算例的比较验证了本文方法的可靠性和正确性.在此基础上,对文献[8]中介绍的一种混合翼型的相关参数变化对扑翼气动特性的影响开展了深入的研究.     

非定常N-S流中的瞬时临界点

本文用微分方程的相似理论把非定常轨线方程变换成自洽系统,考察变换后临界点邻域在非定常流场中的映象,进而对非定常流中“瞬时临界点”溉念的可用性及其限度作了讨论。     

确定性涡方法的精确度研究

本文以无界流场中圆形涡扩散问题为模型问题，对几种为典型的确定性涡方法进行了讨论，分析了方法的参数条件，数值结果表明，确定性涡方法比随机涡方法有更高的精度。