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

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

新的螺肇桨尾涡近似方法

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

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

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

新的螺旋桨尾涡近似方法

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

一种水平轴风力机尾流模型及其计算方法

基于现有的尾流模式,提出一种新的尾流计算方法.新的尾流模式将尾流分为近涡区和远涡区,在近涡区,尾涡在选定的扇形截面上解析求出,在其他处通过非线性插值求得.在尾涡中选取有限的离散点数计算诱导速度,减少了计算量,根据流体连续性方程确定尾涡的扩散半径,避免了梢涡半径的不收敛性.通过具体的风力机算例,与刚性尾流模型相比,该方法计算的轴向诱导因子其分布趋势更加合理,验证了该方法的正确性和实用性,为水平轴风力机风轮气动设计和性能计算提供一种有效的尾流模型计算方法.     

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

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

基于RANS法的B系列对转螺旋桨敞水性能数值模拟

应用RANS法研究B系列对转螺旋桨敞水性能,对单桨和对转螺旋桨进行数值计算,研究不同湍流模型对数值结果的影响,分析研究对转螺旋桨桨距比和直径比的参数匹配问题,比较对转螺旋桨和单桨的敞水效率,剖析尾流场的变化情况,突破传统势流理论方法无法计入黏性和微观流动机理难以分析和把握等局限.数值结果表明:RSM模型的数值精度较高,且选取桨距比和直径比合适的对转螺旋桨明显提高其推进效率,对转螺旋桨后桨能吸收前桨轴向速度和尾涡能量,减小切向速度对流体扰动的影响,增大流经桨叶的流体动量,减小前桨尾流外直径,增加螺旋桨的推力,提高对转螺旋桨推进性能,具有工程实用性.     

空化流中具有升沉运动状态的螺旋桨数值模拟

基于雷诺平均纳维-斯托克斯方法对空化流中具有升沉运动状态的螺旋桨进行了数值模拟,螺旋桨旋转与升沉运动的耦合采用自主定义的运动方程来实现,升沉运动简化为基于正弦函数运动规律的周期性运动,非定常流场的数值传递采用重叠网格技术实现.对不同升沉运动周期下的螺旋桨非定常推力系数、扭矩系数进行了分析,对螺旋桨的非定常空泡性能以及尾流特征进行了对比研究.结果表明:升沉运动会加剧推力系数和扭矩系数的非定常特性,并导致螺旋桨片空泡的非均匀分布;升沉运动周期越小,桨后涡结构的干扰效应越强,螺旋桨诱导的脉动压力幅值越大.计算结果对螺旋桨的优化设计具有一定的理论参考意义.     

考虑旋转和曲率修正的螺旋桨梢涡流场模拟

应用显式代数雷诺应力湍流模型对螺旋桨尾流中梢涡流场分布进行了数值研究,为了避免过高地预报梢涡涡核内湍流黏性耗散,对湍流模型进行了旋转和曲率修正．应用全六面体网格对螺旋桨计算域进行网格划分,为了避免数值离散误差,对梢涡区域进行了网格加密处理．计算结果表明：提出的尾流中梢涡流场分布数值模拟方法能够准确预报螺旋桨梢涡流场的分布及涡核位置,并准确反映了梢涡形成和发展过程中梢涡内主涡和次涡的关系,与实验测量结果基本一致．     

水下螺旋桨尾流水面响应特性试验研究

用表面粒子图像测速技术,对水下螺旋桨旋转产生尾流传播至表面的流动特征进行试验研究.结果表明,螺旋桨作用下的自由表面流场响应是一个非定常的演化过程,总体可分为初始扰动和漩涡演化两个阶段.在初始扰动阶段,水面观测区域出现一个速度辐射中心、一个速度汇聚中心及两个区域间的x负方向主流动,涡量反差数有一个短时间的小幅度升降过程.在漩涡演化阶段,水面观测区域出现一对旋转方向相反、涡量大小几乎相同的涡对,并向x正方向传播,涡量反差数有一个长时间的大幅度升降过程.对不同转速和动量源的不用作用时间的系列试验结果表明:水下螺旋桨尾流传播至表面产生的漩涡强度特征可用涡量反差数峰值界定,而且涡量反差数峰值与限制数存在幂函数关系.     

叶尖开孔对风力机流场的影响研究

以NREL Phase VI叶片的1/8缩比模型为研究对象,在叶片叶尖区域设计由前缘到叶尖端面的3个环形通气孔,改变叶尖流场分布.采用CFD的方法,通过转速变化分析叶尖表面的压力分布情况及其叶尖涡的发展过程,进而研究叶尖开孔对风力机叶尖涡的影响.研究结果表明:转速低于900 r/min时,叶尖开孔对叶片气动性能影响不大;而转速高于900 r/min时,叶尖开孔可降低涡核强度,加速叶尖涡耗散,提高叶片气动效率.从环形通气孔中喷射的气流对来流有明显的抑制作用,能够减小尾流区内的轴向速度.在加速叶尖涡的耗散和降低叶尖涡的强度方面,风力机叶尖处开孔在转速超过900 r/min以上时被视为一种比较有效的设计.     

螺旋桨空化初生的判定和空化斗的数值分析

摘要：
为了实现空化多相流模拟在螺旋桨空化初生和空化斗预报中的应用,采用改进Sauer空化模型和修正剪切应力输运(SST)湍流模型,对NSRDC4381桨的叶背和叶面片空化初生、梢涡空化初生和空化斗底线特征进行了模拟、校验和分析,提出了“当σ>σi时,叶梢截面压力系数分布不再改变”的空化初生判定准则,并分析了其适用范围. 同时,将空化斗由传统的三区细分为五区,以更准确地描述不同工况点的桨叶空化状态. 结果表明,在设计和非设计进速系数下,预报桨叶叶背和叶面不同叶截面的初生空化数均与实验值吻合较好;模拟梢涡空化初生再现了局部梢涡空化和导边梢涡空化同时存在的现象,证明了分别用叶背07 R截面和叶面04 R截面的片空化初生曲线来表征叶背和叶面可视梢涡空化初生的合理性. 空化初生判定准则适用于空化斗底线上方的叶背和叶面片空化区. 将空化斗底线下方区域进一步细分为叶背03 R至叶梢的片空化区以及叶背与叶面同时片空化区,以更好地在工程中应用空化斗来描述桨叶空化形态和求取空化航态的航速.     

螺旋桨性能计算及设计的升力面方法

本文利用Kerwin尾流模型,对螺旋桨性能进行升力面数值解,并进行了理论设计螺旋桨的研究,编制了相应的螺旋桨性能计算、理论设计的计算机程序,给出了对AU 型、DTNSRDC 桨的计算实例,并与相应的试验结果和其它方法的计算结果分别进行了比较。同时,给出了一个螺旋桨的理论设计结果。性能计算及理论设计的数值结果表明了本文计算方法的可行性。     

翼剖面二相流压力分布计算

在船用螺旋桨桨叶表面充气以形成气膜是改善其空泡特性的有效方法之一．应用薄边界层理论和源汇分布面元法实现了水气二相流在叶剖面上的流动匹配计算,确定了水动力压力分布和空泡斗,揭示了空泡性能变化的机理．通过同NACA4412非对称叶剖面和NACA0012对称叶剖面压力分布试验结果的比较表明,数值模拟计算结果是可信的,从而为半潜桨和表面切入桨设计的水动力性能计算提供了有效的数值方法．     

苍鹰翼尾缘结构的单元仿生叶片降噪机理研究

利用逆向工程方法提取苍鹰尾缘非光滑形态的降噪特征元素,由此建立了仿生叶片结构模型;采用基于Smagorinsky亚格子应力模型的大涡模拟,结合基于Lighthill声类比的FW-H方程,分别对仿生尾缘锯齿叶片和标准叶片的流道模型进行了三维流场及声场的数值计算;通过分析仿生齿形结构对叶尾迹流场的影响,研究了仿生尾缘齿形结构的气流噪声控制机理.结果表明:仿生尾缘锯齿结构叶片的总A计权声压级比标准叶片降低了9.8dB;叶片尾缘锯齿结构可以改变流场噪声峰值的分布规律,从而降低了噪声峰值,且大部分频率范围内的气动噪声均有所降低;仿生尾缘锯齿结构可以改变各截面尾迹涡的脱落位置,从而增大了涡心之间的距离,抑制了脱落涡对尾迹流动的扰动,进而减小了叶片表面的非定常压力脉动和尾迹涡引起的气动噪声.     

螺旋桨初生空化湍流的多相流数值模拟

摘要：
同时采用修正剪切应力输运(SST)湍流模型和Baseline雷诺应力模型（RSM）求取了E779A螺旋桨在无空化状态和初生空化状态下的梢涡运动轨迹,分析了涡核最小压力系数、湍动能、轴向速度分量和涡核半径沿运动轨迹的变化,并从模拟得到的梢涡卷曲起始和梢涡涡束的角度阐述了梢涡形成机理.空化模型采用改进Sauer模型,考虑了非凝结性气核质量分数、体积分数和气泡初始半径以及湍流脉动的影响,并针对轻度、中度和重度空化面积进行了可信性校验.当空化数σ>初生空化数σi时,叶梢截面压力系数分布相对不再改变的判定准则来确定.涡核中心位于螺旋线垂向截面上最小压力点,涡核边界由湍流涡频率峰值决定.数值模拟结果表明,RSM模拟梢涡路径较修正SST湍流模型稍长、局部梢涡空化范围略大、叶梢最小压力系数和轴向速度分量要小,涡核湍动能分布更为合理.但两者模拟得到的涡核运动轨迹几乎重合,并且初生空化状态下的涡核运动轨迹、最小压力系数和轴向速度分布均与各自无空化状态下非常接近,表明了初生空化状态判定的正确性和改进数值模型对梢涡运动轨迹模拟的适用性.     

水平轴风力机尾流湍流结构的时空演化

基于小波变换的分析方法,结合致动线模型和大涡模拟研究了一台33 kW水平轴风力机尾流湍流结构的时空演化过程.研究发现,随着距风轮平面距离的增大,尾流中各测点的平均速度先减小后逐渐增大,速度波动的幅值呈减小趋势;风轮后7倍直径内,速度曲线具有明显的周期性,反映出脱落涡通过频率为1.80 Hz,其为风轮旋转频率的两倍.风轮后1倍直径测点处的叶尖涡所在的频率为0.78~25.00 Hz,形成的涡管通过该测点的时间约为0.32 s,涡管直径约为1.83 m;3倍直径测点处出现了0.15~0.78 Hz的低频率湍流结构;7倍直径测点处叶尖涡的频率为1.56~25.00 Hz,相比7倍直径测点之前的叶尖涡频率范围有大幅减小;8倍直径测点处,与近尾流区域相似的叶尖涡的涡管形状消失;9倍直径测点处叶尖涡基本完全耗散.     

应用非定常升力面理论计算大侧斜桨性能及桨叶表面压力分布

应用涡格分布数值计算法，预报大侧斜桨性能及桨叶表面压力分布，采用了两种计算发愤流模型，编制了计算程序，进行了实例计算并与已有的试验结果进行了比较，其结果吻合良好。     

潜艇涡量场和流噪声等效声中心的数值预报

为了实现潜艇湍流噪声及其等效声中心的数值预报,在分析SUBOFF潜艇拖曳和自航状态下涡量场的基础上,采用大涡模拟与声学边界元相结合的方法,在频域内预报了流噪声空间分布、测点谱源级曲线和声指向性,求取了等效声中心位置并分析了其受螺旋桨旋转作用的影响.计算结果表明：附体与艇体结合部马蹄涡和附体端面诱导项链形涡对是潜艇涡量场的主要特征,且马蹄涡系具有较高的强度和稳定性;附体尾涡脱落频率存在19.22Hz的线谱,且在尾涡测点谱曲线中得到明确体现;随着频率增加,流噪声蝶形指向性对应的辐射瓣状区间数随波数增加,且正横方向声压要强于首尾方向;流噪声等效声中心位于距艇艏0.46倍艇长处,在10Hz～1kHz内总声源级为95.09dB;艇艉桨对附体马蹄涡系影响较小,但促使等效声中心迅速移至艇艉.