三维机翼尾涡卷曲的数值模拟

为使三维机翼尾涡卷曲更加真实,在计算诱导速度时引入衰减函数和有限半径涡核分别体现流体黏性作用下尾涡强度的衰减和尾涡扩散.应用面元法解析式计算奇点的诱导速度,在偶极附近区域用涡环代替偶极,利用尾涡面上的运动学条件来确定尾涡面的卷曲形状.数值模拟结果表明:涡核半径参数有一最佳值;考虑尾涡衰减可使卷曲更接近实验结果;采用曲线拟合方法确定衰减函数是可行的.     

新的螺肇桨尾涡近似方法

基于解析尾涡模式，提出了一种新的尾涡近似方法，尾涡分为近区和远区，在近区，尾涡在选定扇形截面上解析地求出，在其他处理过非线性插值求得；利用了解析尾涡的特点同时减少了计算量，根据流体连续性方程确定尾涡的收缩半径，避免了梢涡半径的不收敛性为重负荷螺旋桨设计和非设计工况下螺旋桨性能计算提供了一种有效的尾涡近似计算方法。     

新的螺旋桨尾涡近似方法

基于解析尾涡模式,提出了一种新的尾涡近似方法．尾涡分为近区和远区．在近区,尾涡在选定扇形截面上解析地求出,在其他处通过非线性插值求得;利用了解析尾涡的特点同时减少了计算量．根据流体连续性方程确定尾涡的收缩半径,避免了梢涡半径的不收敛性,为重负荷螺旋桨设计和非设计工况下螺旋桨性能计算提供了一种有效的尾涡近似计算方法．     

船用螺旋桨理论及其应用研究进展

论述了船用螺旋桨理论设计与计算研究的进展;讨论了升力面理论方法中尾涡模式函数,面元法中的随边库塔条件的数值处理,以及纳维-斯托克斯方程方法中流域计算匹配的数值方法;给出了相关理论方法在螺旋桨水动力性能黏性修正、翼剖面二相流水动力计算和新型螺旋桨叶剖面设计开发研究中的应用实例．应用计算表明,所提供的数值方法是有效的．有关螺旋桨水动力性能的黏性修正、叶剖面二相流计算和桨叶叶剖面设计等方法,对于工程应用具有参考价值．     

螺旋桨非定常性能计算的升力面方法

本文利用非定常数值升力面方法计算了运转于非均匀流场中的螺旋桨性能,该方法建立在作者提出的新的定常螺旋桨尾涡模型的基础上,在用圆锥螺旋面模拟的过渡区尾涡片中,加入了泄出涡以反映非定常特性.数值计算的升力面方法采用离散的涡、源奇点系,并编制了相应的非定常性能计算的计算机程序.最后,对14届ITTC 螺旋桨委员会提供的实例以及其他例子作了比较计算,结果表明,本法是可行的.     

基于大涡模拟的螺旋桨梢涡数值分析

为研究螺旋桨的梢涡脱落规律,基于STAR-CCM+软件平台,采用重叠网格技术和大涡模拟方法,通过求解纳维-斯托克斯方程,对螺旋桨的速度场和涡量场进行了详细的研究,并且比较了不同进速系数下螺旋桨梢涡形态的差别.数值研究结果表明:基于正交网格的大涡模拟能够较好地模拟螺旋桨的梢涡形态,低进速系数下螺旋桨的梢涡强度增大并且显得比较紊乱;高进速系数下螺旋桨梢涡则呈现出比较规则的螺旋涡管形状.     

用升力面方法估算螺旋桨空泡

本文提出了用数值升力面理论求解三维空泡螺旋桨问题的准定常方法.表征螺旋桨负荷、厚度和空泡的奇点系及满足边界条件的控制点均置于拱弧面上。尾涡模型划分为过渡区和远尾流区,过渡区用圆锥螺旋面来模拟尾涡片的变形现象,远尾流区简化为桨叶数的集中等螺距梢涡和一根直线毂涡模型。空泡模型用半闭式.计算结果与实验值吻合良好。     

基于螺旋加密网格的螺旋桨梢涡空化数值模拟

基于计算流体力学(CFD)Star-CCM+软件,以船用E779A螺旋桨为研究对象,采用大涡模拟方法和Schnerr-Sauer空化模型模拟螺旋桨尖端梢涡空化.在螺旋桨可能发生尖端梢涡空化的位置上使用螺旋管为体积控制,形成螺旋加密网格,对E779A螺旋桨的梢涡空化进行更为精细的数值模拟.实验和数值模拟的水动力结果表明:本文方法具有良好的数值预报精度,推力、扭矩系数和敞水效率较为相符.在此基础上进行基于螺旋加密网格的梢涡空化数值模拟,通过模拟结果可以看出:螺旋加密区的流动得到了更为精细的模拟,并清晰展示了加密梢涡的速度和压力分布.数值和实验的空化结果比较表明:本文方法能够较好地预测螺旋桨的空化,尤其是梢涡空化的形态和规律.     

非线性涡格法预报桨后舵附推力鳍水动力性能

对螺旋桨与舵附推力鳍分别采用升力面法和非线性涡格法计算.螺旋桨、舵附推力鳍两者之间的相互干扰采用迭代计算.数值计算过程中考虑了推力鳍端部分离涡的影响,提高了理论预报的准确性.螺旋桨尾流区分为过渡区和远尾流区.过渡区长度取3.0D,以使舵附推力鳍完全处于螺旋桨尾流的过渡区内,过渡区采用圆锥螺旋面来模拟涡片的变形现象.对影响推力鳍助推效率的几个主要参数进行了变尺度研究.并将结果与前人的计算结果进行了对比,计算结果显示螺旋桨后的舵附推力鳍助推效率随着安装角的改变而显著变化.存在最佳安装角,大约为5°,离开这个最佳安装角,推力鳍的助推效率将下降;推力鳍的展长与螺旋桨半径之比在0.9左右时推力鳍的助推效率最高;螺旋桨进速系数越小,推力鳍的助推效率越大.     

用非定常面元法预报拖式吊舱螺旋桨水动力性能

用非定常面元法对拖式吊舱推进器的水动力性能进行求解,吊舱和螺旋桨的相互干扰通过诱导速度迭代处理.与将相互诱导速度在周向上进行平均的处理方法相比,文中的方法真实考虑了螺旋桨和吊舱相互影响的非定常性.为了节省计算时间和存储空间,对螺旋桨的求解只在主叶上进行.结果表明,该方法能够对吊舱推进器在均匀来流和在非均匀来流中的水动力性能进行计算.将敞水条件下的计算结果同试验结果进行了比较,两者有较好的一致性.     

脉动流条件下圆柱绕流特性研究

该文数值研究了脉动流条件下圆柱绕流场特性,对涡量分布、升阻力系数、升力系数频谱特性等进行了分析。结果表明,脉动流作用会使圆柱尾涡剪切层变薄,提高脉动频率使尾涡长度变短、脱离加快,提高无量纲脉动振幅使主导涡的脱落速度降低,使旋涡生成区域更靠近圆柱表面;脉动频率与无量纲脉动振幅的增大使升、阻力系数的振幅均增大,而且与升力系数相比,阻力系数的振幅更大,变化更快;升力系数频谱图存在多个主频,包括旋涡脱落频率和相位叠加频率,其中旋涡脱落频率的振幅随着脉动频率的增加而减小,随着无量纲脉动振幅的增加而增加。     

非线性广义变延迟方程稳定性分析

主要讨论了非线性广义变延迟方程的稳定性.首先讨论了基于模型方程理论解渐近稳定的条件,其次研究了Runge-Kutta方法求解方程数值解的GAR(l)-稳定性,最后的数值算例验证了理论结果的正确性.     

方腔中多个涡斑运动的数值模拟

通过对Navier-Stokes方程的直接数值模拟,应用涡斑模型研究了方腔内多个涡斑的运动,结果表明,旋转方向相同的涡斑,不管初始涡斑的涡量分布如何,最后的涡量分布为Burgers型涡斑,呈现“单胞”的稳定结构,旋转方向相反的涡斑,其运动演化轨迹与初始条件有关,最终呈“双胞”的稳定结构。     

驻涡燃烧室内涡系分布研究

为清晰了解驻涡燃烧室内流场的涡系分布,对中心驻体为67%燃烧室当量宽度的驻涡燃烧室流场,在驻涡腔无喷射和有喷射两种状况下分别进行冷态数值模拟,详细研究流场内旋涡分布情况.结果表明:流道前部存在4个影响很小、沿流动方向运动、分布于驻体与端壁角区的螺旋状旋涡;流道后部存在1个由2个旋涡构成、面积较大的回流区;流道中部的驻涡区旋涡情况复杂,在驻涡腔无喷射情况下,主流与驻涡区之间悬挂驻留2个稳定的驻涡,驻涡腔内存在不稳定且数目不确定的旋涡;当驻涡腔前驻体有喷射存在时,驻涡腔内形成4个均匀且稳定的旋涡.     

刚性常微分方程组初值问题的多重网格解法

描述了求解刚性常微分方程组初值问题的多重网格方法。具体选择了光滑算子构造了有关的多重网格分量,阐明了方法的收敛性。它不仅具有隐式Runge-Kutta方法的优点,而且可以不必求解由隐式方法产生的非线性方程组。给出的计算实例显示了快速收敛性。     

用EXCEL解微分方程

文章介绍用 MICROSOFT EXCEL在计算机上求一阶常微分方程和二阶常微分方程的数值解。在 EXCEL的界面下 ,欧拉折线法或龙格 -库塔法的一系列繁杂计算可以简单地通过对 EXCEL单元格的自动填充来实现 ,所得到的数值解具有很高的精度。 EXCEL 的良好图表功能还可以方便地给出求解结果的图像。     

基于拉格朗日方法的水翼尾缘非定常涡旋结构研究

湍流中涡旋结构的非定常脱落是普遍存在的,其脱落规律的研究具有广泛的工程应用背景,探索有效的涡旋结构预测方法具有重大意义. 采用基于拉格朗日体系的有限时间李雅普诺夫指数方法对水翼尾缘附近涡旋结构的非定常发展过程进行研究. 采用大涡模拟方法对水翼周围湍流场进行计算,通过与实验结果进行对比,验证数值模拟方法的合理性与准确性. 并在此基础上,对水翼尾缘的非定常涡旋结构进行分析. 结果表明:尾缘涡旋非定常动态行为导致水翼动力特征准周期性波动. 有限时间李雅普诺夫指数方法可以准确描述尾缘涡的非定常流动细节,拉格朗日拟序结构可以捕捉水翼尾缘涡旋结构和尾迹涡街结构的边界.      

变压器励磁支路暂态的非线性RK方法

将非线性显式ＲＫ方法用于求解变压器励磁支路暂态过程。该方法具有Ｌ稳定性。通过一个算例,将该方法与显式,隐式,半隐式ＲＫ方法进行了比较。结果表明,该方法不仅具有较高的数值稳定性,而且避免了迭代运算或Ｊａｃｏｂｉ矩阵的求解。相对于隐式方法,又可以大大缩短计算时间。     

迎风条件下速度梯度对地面涡气动特性影响数值模拟

以缩比的近地面短舱进气道为研究对象,通过数值计算模拟的方法得出在迎风来流条件下速度梯度对地面涡流场的结构及气动特性曲线的影响。研究表明,速度梯度是影响地面涡的主要因素。在不同的速度梯度下地面涡结构明显不同。在Y方向速度梯度下近地截面形成的地面涡是一对涡,旋转方向相反。速度梯度越大,地面涡强度,总压畸变指数明显变大。且峰值会出现滞后现象。在Z方向速度梯度下,近地截面形成的地面涡是一个单涡结构,在进气道上方也会形成一个涡体。地面涡强度、总压畸变指数会先增大后减小,峰值不会出现滞后现象。