影响导管桨内部流场的几个因素

用理论方法计算导管桨的内部流场,研究几个重要因素对内部流场的影响,以期获得导管桨内部的流体运动规律．用迭代方法分析螺旋桨和导管的相互影响,迭代计算结束后求解螺旋桨和导管表面的偶极强度和源汇强度,根据蒙瑞诺公式计算螺旋桨和导管的诱导速度,即可得到导管桨的内部流场．以面元法为数值求解方法,编制计算程序,研究影响导管桨内部流场的几个因素,包括螺旋桨的盘面比、叶数和导管的形状、位置．结果表明,上述参数对螺旋桨内部流场的大小和周期均有很大的影响．     

带制流板舵在螺旋桨尾流中的水动力性能

提出了一种计算带制流板舵在螺旋桨尾流中水动力性能的方法，舵的水动力及周围流场用面元法计算，螺旋桨性能及尾流场通过无限叶数的简易螺旋桨理论来预估，而舵上下制流板的影响则应用升力面的涡格法来计算，舵、制流板、螺旋桨的干扰作用以迭代方法求得，由于采用了面元法，舵的水动力性能及表面压力分布的计算更为精确，文中不对带制流板舵的操纵性能进行计算和研讨，计算结果与实验值比较，吻合程度良好。     

襟翼舵水动力性能研究

对襟舵在螺旋桨尾流中的水动力性能及其桨舵干扰进行了研究，舵的水动力及周围流场用面元法计算，螺旋桨性能通过无限叶数的简易螺旋为预估，桨舵干扰作用以迭代方法求得，由于采用限面元法计算襟翼舵的性能，可以更为精确地反映较为复杂的襟翼舵表面形状对水动力性能的影响，文中还对螺旋桨尾流中反应舵的表面压力分布及操纵性能进行了计算和研讨 。     

迭代型体积力法预报船舶推进性能

摘要： 应用计算流体动力学的方法（CFD）,采用迭代型体积力法,耦合螺旋桨势流涡格法（VLM）和船体周围黏性流场RANS（Reynolds Averaged Navier Stokes)方程,代替真实螺旋桨的作用,即将VLM法计算得到的推力和转矩分布插值后作为源项加载到RANS方程中,将RANS方程得到的流速扣减诱导速度后得到实效伴流作为VLM法计算的进流条件,计算得到新的推力和转矩值,重复迭代直至螺旋桨的推力变化达到收敛要求.在无舵情况下对集装箱船KCS(KRISO Container Ship)的自由面黏性流场进行了数值模拟,计算结果与试验结果吻合良好.整个计算过程由集成了VLM和RANS方程的程序自动完成.研究结果表明,迭代型体积力法可以快速准确地预报船桨干扰下船体的推进性能.     

基于CFD的桨舵水动力干扰研究

为了观察舵对桨叶表面压力分布以及桨后尾流场的影响,以B4-55桨和展弦比为1.0的NACA0020的流线型舵为算例,采用计算流体力学方法分析了桨舵干扰的水动力特性。为了提高网格质量,计算域被分为包含螺旋桨的旋转域、包含舵的静止域以及除二者之外的静止域3部分,并采用不同的网格策略对其进行离散。离散方程用K-ω模型求解,旋转域与静止域之间的相对运动通过MRF技术模拟,并利用交界面技术实现数据传递。分析了不同进速系数下桨舵干扰的规律以及舵攻角对螺旋桨的推力、转矩的影响,研究了流场的速度以及压力分布。数值结果和试验结果的比较表明,采用的计算模型和离散方法能够得到准确的计算结果。     

一种对船桨干扰问题的黏势流耦合求解方法

基于黏性流计算流体力学(CFD)及势流理论,建立了一种黏流雷诺平均方程/面元法耦合迭代求解方法,对船桨相互作用问题进行了数值模拟.自由面的船体绕流场模拟应用CFD黏流方法,螺旋桨水动力载荷计算采用势流面元法.在黏性流场模拟中得到桨盘面处的总速度分布,减去螺旋桨诱导速度得到实效伴流速度并作为势流计算速度进口条件,螺旋桨效应以体积力的形式在黏性流场中体现.计算值与试验值的对比结果表明,黏势流耦合迭代求解方法能较好地预报船桨干扰问题,且较全域船桨离散化网格CFD求解船桨干扰的方法,计算量大为减少.该方法充分融合了黏势流方法的优点,既能计及黏性效应又能发挥势流快速准确的优势,能够拓展应用于船体及螺旋桨性能预报及设计优化的研究.     

抗空化扭曲舵的设计及其水动力性能分析

摘要： 通过面元法计算螺旋桨尾流场,根据螺旋桨尾流场对扭曲舵各剖面的几何攻角进行设计,提出了抗空化扭曲舵的设计方法;采用计算流体动力学方法评价舵的水动力性能,并经过拖曳水池的舵力测量试验加以验证.结果表明,在不同的航速和舵角条件下,与普通舵相比,扭曲舵的负压力系数峰值明显较小,从而大幅提高了舵的抗空化性能,降低了舵的空化剥蚀和由舵产生空泡所引起的振动噪声.由于扭曲舵与螺旋桨的尾流匹配良好,在舵角为0°时,扭曲舵叶背的压力系数分布与叶面的压力系数分布基本重合,因而显著减小了直航时舵上的横向力和舵轴力矩,有利于舵和舵机的受力.     

基于点源模型的螺旋桨负载噪声频域预报

在点源模型基础上,将旋转力源离散为均匀分布在旋转轨迹上的有限个固定偶极源,各偶极源间的时间差转换为对应的相位差,结合边界元方法可以计算任意边界条件下的旋转声源频域负载噪声.以E779A船用侧斜螺旋桨为对象,在稳态流场计算得到验证的基础上,通过瞬态流场数值计算得到了不同时刻螺旋桨表面每个单元对应面力.将旋转面力沿运动轨迹离散后,首先计算了单个桨叶的负载噪声,计算结果同文献值符合较好;然后计算了整桨负载噪声,桨盘面内测点轴频对应的辐射声压最大,轴线方向测点叶频处辐射声压最大,负载噪声在桨轴平面指向性呈类8字形,表现为偶极子声源特性.     

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

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

基于面元法预报导管桨性能的数值计算方法

基于面元法,分别用将桨和导管当作一个整体考虑的直接计算方法、基于诱导速度的迭代计算方法与基于诱导速度势的迭代计算方法预报了导管螺旋桨的定常水动力性能.计算结果表明：采用这3种方法预报导管螺旋桨的定常性能具有相似的精度,但直接求解方法可以节省大量的计算时间;直接求解和速度势迭代计算的导管表面压力分布、桨叶表面环量和压力分布基本一致,但是与诱导速度计算的结果有一定的差异.对这3种方法计算结果的差异性分析表明：计算导管桨水动力性能时,直接求解方法和基于诱导速度势迭代方法更可靠.     

螺旋桨尾流流场的数值计算

主要研究拖船螺旋桨对海水流场的影响，计算了最大功率为1250kw／台双导管螺旋桨拖船系柱推力和轴向诱导速度，并利用PHOENICS与软件对该拖船双导管螺旋桨尾流流场进行了数值计算．计算结果表明：在系泊工况下，长80m、宽15m、深5m范围内，拖船螺旋桨对海水流场影响较大，螺旋桨后112．5m的海域，螺旋桨的影响已不明显．     

潜艇尾流场中螺旋桨非自由声场的预报

为在频域内分析潜艇艇体声散射效应对尾流场中螺旋桨声场的影响,以小比例潜艇和大侧斜螺旋桨为对象,首先基于计算流体力学(CFD)方法模拟了"艇+桨"的推进特性,得到了螺旋桨壁面的脉动压力,然后考虑艇体的散射效应,计算分析了螺旋桨的声场特性.结果表明:当以"艇+桨"系统为对象真实模拟螺旋桨的工作条件时,由于潜艇尾部十字舵的存...     

基于速度势面元法的风力机风轮三维气动性能预估

大型风力机风轮气动设计和性能分析需要更加准确、快速的预估方法。风力机风轮三维流场具有低Mach数、高Reynolds数的特点,该文把风轮流场简化为有势流场,将叶片表面和尾迹视为有势流场的边界,采用速度势方程作为流动控制方程,使用面元法进行求解,建立了风力机风轮三维气动性能数值计算的模型,并编制了计算程序。同时,还建立了适用于水平轴风力机的尾迹迭代求解方法,考虑了表面摩擦阻力对风轮转矩和推力的影响,提高了计算精度。与实验结果对比证明,该方法具有很好的计算准确性,为进一步研究常规运行工况下的风力机风轮气动性能,特别是大型变速变桨风力机气动性能提供了一个准确、快速的预估方法。     

桨舵系统非定常水动力性能预报方法研究

基于势流理论,提出了一种求解桨舵系统非定常水动力性能的方法.系统地研究了流场中含2个升力体的面元法理论,建立了相应的数值计算模型.模型中螺旋桨和舵同时求解,二者间的相互干扰通过时域内影响系数的变化进行考虑.以定常计算结果作为非定常计算的初始值,考虑了非定常求解时泄出涡的时间效应.计算了均匀流及非均匀流中桨舵系统非定常水...     

用HH-SIMPLE方法求解明渠中三维浅水底孔泄流问题

本文针对明渠中三维底孔喷射流问题,提出了一个新的求解 N-S 方程的有限差分方法——HH-SIMPLE 方法.用该方法求解三维水流问题,可以保持 SIMPLE 方法的优点,即可将动量方程中的速度与压力解锅;求解离散化动量方程时可用隐式迭代格式,而无需作“刚盖近似”或“静水压力假设”,因而可很好地处理自由表面.本文中用该法求解的算例在国际计算流体力学领域中目前尚为空白,其三维计算结果可精确地描述在上游有来流情况下的浅水底孔喷流的流场、自由水面和压力场问题.计算结果和试验资料吻合较好.     

基于分离涡模拟方法的导管桨近尾流场及尾涡特性分析

基于分离涡模拟(DES)方法对设计工况下导管桨的近尾流场及尾涡特性进行数值模拟.数值计算中选用SpalartAllmaras湍流模型封闭N-S方程,采用滑移网格技术及混合网格划分方法完成导管桨敞水性能数值计算.通过分析导管桨瞬态尾流场及尾涡空间结构发现:近尾流场中螺旋桨半径区域瞬态诱导速度大,尾流中分布着连续漩涡结构,尾流加速作用明显.导管桨尾涡主要由导管剪切层涡、叶片涡系及毂涡组成,叶片涡系中包含叶梢涡、叶根涡、毂涡及相邻梢涡带之间诱导产生的S形二次涡;导管桨尾涡结构中多重涡系之间产生复杂干扰,尾涡形态出现融合、扭曲、分解并逐渐扩散.     

瞬变工况下螺旋桨性能

用非定常非线性涡格法计算螺旋桨尾流场预估非定常性能，并进一步计算了瞬变工况下螺钉性能，对ＤＴＮＳＲＤＣ４３８１和４３８３桨的尾流场进行了校核计算，取得了理想的结果；对ＩＴＴＣ建议的轴承力校核计算的螺旋桨，在给定的伴流场条件下计算其非定常性能，与试验结果和其他的理论方法比较，结果令人满意，最后计算了ＤＴＮＳＲＤＣ４１１８桨在变工况下性能变化过程，显示了该方法用于螺肇钉性能预估的前景。     

CLT桨的尾流场及梢涡特性数值分析

为了探究尾流收缩叶梢负载(CLT)桨的尾流场及梢涡特性,基于商业计算流体力学软件Star-CCM+对CLT桨P1727的尾流场进行了数值模拟.计算采用DDES(延迟分离涡)方法在进速系数J=0.5的工况下对三套网格进行了不确定度分析,确认了数值方法的可行性,然后对螺旋桨的敞水性能及尾流场进行了计算.研究表明:尾流场的网格细化对螺旋桨的推力及扭矩只有微弱影响,对桨后的流场具有较大影响;尾流区域轴向速度分布分为加速流与自由流,进速系数越大,加速流与自由流的分界线向外扩散的趋势越弱;相邻梢涡之间会发生融合,进速系数越大,融合得越晚,梢涡强度越弱;CLT桨相较于常规桨多出一个端板涡,并且会与梢涡融合;梢涡在向下游发展的过程中,会有明显的收缩现象.     

螺旋桨流固耦合特性的数值模拟

基于ANSYS Workbench平台,利用其多场分析功能,采用求解RANS方程的方法,利用CFX求解器对螺旋桨3维流场进行数值模拟,利用有限元求解器计算螺旋桨结构响应,并将有限元求解器的设置输出为结构求解子程序,供CFX计算时调用以实现数据的实时双向传递,从而将流体计算与结构变形计算耦合起来,给出一种求解螺旋桨流固耦合问题的3维数值模拟方法,并通过与试验数据的比较验证了方法的可靠性.设定2种材料特性,对螺旋桨的水动力性能、桨叶变形及应力应变特性等进行了数值模拟.结果表明:材料特性对变形及应变较为显著,弹性桨的变形与应变较刚性桨约高1个量级;变形后的螺旋桨又对水动力性能及周围流场产生影响,弹性桨的推力和转矩均比刚性桨小.     

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

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