基于两相流欧拉方法的翼型结冰数值模拟

文章通过计算流体力学的方法对翼型结冰进行了预测。针对拉格朗日法求解水滴收集特性计算工作量太大的缺点,引入水滴流场的概念,通过求解空气-水滴两相流场的欧拉方法,获得翼型表面水滴收集特性。结冰计算过程中网格采用变形网格技术得到,既保证了网格质量又减少了网格生成的工作量,采用有限体积法及LU-SGS隐式时间推进数值求解N-S及水滴运动方程。建立了对应于霜冰和光冰的冰增长模型。通过对NACA0012翼型霜冰和光冰的计算验证了该方法是有效的。     

结冰对翼型和多段翼型绕流及气动特性影响研究

对翼型结冰问题，在结冰条件下，利用数值方法分析研究翼型和多段翼型绕流流场及气动特性的影响。对于翼型和多段翼型，采用分区多块网格技术，结合法向外推的代数方法和求解椭圆型方程方法，生成高质量计算网格。采用中心有限体积法，以及Runge-Kutta显式时间推进格式，结合B—L代数湍流模型，运用流场分区求解算法，完成了绕流流场的N-S方程数值模拟，进一步针对3种不同的冰型：钝头体、双角体和尖头体冰型，分析不同形状的冰型对翼型和多段翼型绕流流场及气动特性的影响。     

多段翼型结冰数值模拟研究

文章采用计算流体力学方法对多段翼型进行了结冰数值模拟研究。应用多块网格技术生成多段翼型网格。引入水滴流场的概念,通过求解空气-水滴两相流场的欧拉方法,获得翼型表面水滴收集特性。流场采用有限体积法及LU-SGS隐式时间推进数值求解,运用流场分区求解算法,完成了多段翼型绕流流场及水滴流场的数值求解。建立了冰层的法向增长模型,对冰形进行了预测。     

基于近似边界条件的翼型跨声速颤振数值模拟研究

采用近似边界条件对二维翼型的跨声速颤振进行了数值模拟.计算网格采用静止的笛 卡儿网格,通过薄翼和机翼小变形假设,使用一阶近似边界条件求解欧拉方程;由于文中方法 在计算过程中不需要重新生成网格,节省了大量的时间.计算了Isogai机翼模型的颤振边界, 与参考文献的计算结果吻合良好,证明了基于近似边界条件的跨声速颤振数值模拟方法的高效性和准确性.     

基于CFD方法的阵风响应与阵风减缓研究

通过设计舵面的同期运动,可以达到缓解飞行器遭遇阵风时引起的非定常气动力与力矩干扰的目的。文章通过引入"网格速度"模拟阵风条件,求解非定常欧拉方程实现了翼型阵风响应的数值模拟;并通过动态嵌套网格方法研究了舵面运动时的阵风减缓效果。首先采用该方法对NACA0006翼型迎角阶跃型阵风的气动力响应进行了计算,计算结果与理论结果、文献计算结果吻合良好。进一步对NLR两段翼型在迎角阶跃型、One-m inus-cosine型阵风做用下的气动力响应过程进行了数值模拟,对襟翼运动时的阵风减载效果进行了分析,并研究比较了襟翼运动参数对阵风减缓的影响。结果表明:通过合理设计舵面运动方式,可以很好抑制阵风引起的非定常气动干扰。     

非结构网格多段翼型绕流Euler及N—S方程数值模拟

描述二维非结构无粘网格和粘性网格的生成方法,对二维翼型绕流进行了Euler及N－S方程数值模拟,运用推进面方法,生成无粘网格,进一步结合推进层方法,生成粘性网格,应用中心有限体积法,采用B－L代数湍流模型,完成了对多段翼型流场的Euler方程和N－S方程的数值模拟计算,计算与实验结果对比表明,所采用的网格生成技术及流场计算方法是正确、有效的。     

非结构旋转动态嵌套网格技术及其应用

发展并改进了一套基于非结构网格的旋转动态嵌套网格方法,通过数值求解Navier-Stokes方程对二维/三维带旋转运动边界的非定常问题进行了数值模拟.将计算区域划分成包含旋转运动边界的运动区子域和覆盖其它静止边界的静止区子域,以人工边界和约束线(面)作为约束条件生成整体计算网格,进而根据设置的人工边界完成嵌套网格的建立.对翼型绕流、旋翼/机身干扰流场以及颤振激励系统绕流进行了验证计算,结果与实验值相一致,表明了文中方法能够有效地对带有旋转运动边界的非定常问题进行数值模拟.     

翼型跨音速定常及非定常粘性绕流的数值模拟

发展了一种计算翼型跨音速定常及非定常粘性绕流的数值方法。采用ＬＵ分解隐式三阶迎风格式求解二维非定常可压ＮＳ方程组，紊流模型采用ＢａｌｄｗｉｎＬｏｍａｘ双层代数模型。计算网格采用相对翼型固定的贴体运动网格。数值模拟了ＮＡＣＡ００１２翼型的几种跨音速定常及非定常粘性绕流状态，计算结果同已有的实验数据吻合较好     

基于影响函数控制的结构动网格生成方法

在静态结构网格的基础上,根据流场边界的运动提出了一种新的结构网格变形技术.通过影响函数的作用,将运动边界的信息向周围的网格传递,流场中的网格节点通过影响函数获取移动变形信息,构建出新的网格.相对于弹性网格变形方法,影响函数控制方法具有更强的网格变形能力和更好的网格质量.该方法能够满足生成贴体的粘性动网格,适用于Navier-Stokes方程.耦合Navier-Stokes方程与强迫运动方程,数值模拟了NACA0012翼型做强迫俯仰振荡(AGARD-CT1)的粘性动态流场.计算结果与试验值的对比,表明了该方法的有效性.     

翼面结冰对翼身组合体气动特性影响研究

飞机翼面结冰会影响气动性能,危及飞行安全.针对翼身组合体的翼面结冰问题,设计3种不同的冰型:钝头体、双角体和尖头体,对结冰前后的绕流流场及气动特性,采用数值方法进行分析研究.利用多块网格技术,结合法向外推代数方法和求解椭圆型偏微分方程方法,生成计算网格,采用中心有限体积法、LU-SGS双时间隐式推进算法,结合B-L湍流模型,运用流场分区求解技术,完成绕流流场的N-S方程数值模拟.计算与试验数据进行对比分析,说明了所采用方法的可行性和正确性,并给出了不同冰型对气动特性的影响结果.     

基于N-S方程的翼型结冰数值模拟

为了研究结冰对翼型气动性能的影响，采用计算流体力学的方法对NACA0012翼型前缘的结冰进行了预测。在流场计算中，对二维可压粘流的时均N-S方程进行离散求解；采用四步龙格一库塔法求解水滴轨迹运动方程；建立了对应于霜冰和光冰的冰增长模型并假定冰沿着与当地翼型表面法向一致的方向增长来预测结冰的形状。分析了网格数和冰层时间步长对最终冰形的影响和结冰对翼型气动性能的影响，最后将计算的结果与文献中的实验值和计算值进行比较，结果一致。     

翼身组合体的网格生成及EULER方程数值模拟

采用代数方法中的广义无限插值法生成翼身组合体三维Ｏ－Ｃ型贴体网格。该方法使用双边界法和引入了物面外法向导数以保证物面附近网格正交的处理方法，提高了代数方法生成翼身组合体网格的质量。并用中心格式有限体积法求解跨声速ＥＵＬＥＲ方程，以某翼身组合体模型为例，计算结果与实验符合良好。     

采用非定常N—S方程的翼型颤振特性分析研究

以非定常N－S方程为主管方程,计算翼型振动的瞬时非定常气动力,并与颤振方程耦合求解,用时间推进的方法,计算了结构响应特性。经多个算例计算,研究了颤振的临界速度随马赫数的变化规律以及极限环振荡等非线性特性。计算结果与其它献计算结果吻合很好。     

架空输电导线三维时变积冰计算

输电线路覆冰严重影响电网安全运行。现有输电线路积冰模型多忽略轴向积冰差异,且将关键积冰参数设定为时不变(单步)条件,鲜见时变参数(多步)三维积冰模型的讨论。本文基于润滑理论及导线覆冰机理,提出了考虑时变积冰参数影响下的输电导线覆冰模型。依托ANSYS-Fluent ICING模块,对三维导线模型开展时变参数积冰计算,利用真实输电线路覆冰试验数据对计算方法的有效性进行了验证,计算结果与试验结果吻合度较好。在此基础上对比研究了单步与多步积冰算法,并着重分析了输电导线倾斜角度、导线直径对积冰形貌、积冰质量的影响。结果表明：多步积冰计算法较单步积冰计算法精度提高约8%。干覆冰条件下,随输电线倾斜角度增大,输电线上积冰形貌和积冰质量均没有明显变化;湿覆冰条件下,输电线倾斜角度对积冰有显著影响,当倾斜角度由0°增至60°,表面积冰覆盖面积减小、积冰形貌逐渐圆滑,但积冰质量降低约21%。对比不同直径导线积冰计算结果,大直径导线积冰量明显高于小直径导线。     

基于插值的Lattice Boltzmann方法非均匀网格算法

提出了一种基于插值的latticeBoltzmann方法的非均匀网格算法 ,给出了计算步骤 ;对典型的二维恒定方腔流、恒定和非恒定圆柱绕流问题进行了模拟分析 ,所得结果的正确性证明了该算法是有效、可行和实用的     

悬停旋翼粘性绕流N-S方程数值模拟

通过求妥薄层Navier-Stokes方程数值模拟了悬停旋翼粘性绕流,为了方便实施旋转对称性边界条件和更为有效地捕捉尾涡系,采用了嵌套网格技术,空间离散方法采用中心有限体积法,时间推进格式为五步Runge-Kutta法,考虑到旋翼绕流的特点,对远场边界条件作了合理的修正,应用本方法对Caradonna模型旋翼的两个实验状态进行了数值模拟,计算结果与实验值吻合很好,说明本 方法 成功的