粘性流中二维水翼局部空泡流的数值模拟

基于粘流理论研究了二维水翼的局部空化。数值计算从完整的N-S方程出发,采用k-ε两方程湍流模型封闭Reynolds方程,作为基本控制方程,空泡流采用Rayleigh plessetmodel模型,计算了0°攻角时,不同航速对水翼局部空泡流的影响,同时也分析了水翼的绕流场,并把计算结果与已发表的试验结论做了比较。结果表明,本文的计算方法具有较好的计算精度。     

螺旋桨非定常片空泡CFD数值方法研究

应用雷诺平均的N-S方程结合改进VOF空化模型模拟了螺旋桨非定常片空泡。应用全六面体网格对螺旋桨模型计算域进行了网格划分,数值计算中通过计算螺旋桨敞水性能对计算方法进行了验证。应用SST-SAS湍流模型及修正的VOF空化模型模拟了HG01螺旋桨非定常片空泡,该数值方法对螺旋桨片空泡的发展过程和形状与实验结果吻合良好。     

不均匀进流场下螺旋桨非定常力的数值模拟

本文基于STAR-CCM+对非均匀来流下的螺旋桨非定常力进行数值模拟,通过与试验值进行比较,分析了时间离散格式、单位时间步长螺旋桨旋转角度与网格密度对计算结果的影响。结果表明,在螺旋桨非定常力计算中,螺旋桨单位时间步长的旋转角度应为1.8°同时选取二阶时间离散格式,桨叶表面网格的边界增长率设定为medium就可以达到计算精度的要求。并在此基础上,对Suboff全附体潜艇模型下3个螺旋桨的非定常力进行数值模拟,分析桨叶个数以及桨叶侧斜对非定常力变化幅度的影响。     

回转体空泡绕流现象的数值模拟

利用FLUENT 6.0流体动力学计算软件对轴对称回转体的自然空化现象进行了数值模拟,得到了不同头型和不同空化数下的回转体周围空泡绕流流场和体积分数分布.计算结果与其他研究人员的计算与试验结果比较吻合,说明利用FLUENT 6.0软件可以较好地模拟空泡绕流现象,能够正确反映空泡外形随空化数变化的规律,并且能够较准确地捕捉到空泡尾端回射流现象.     

E779A螺旋桨斜流工况下的空泡数值模拟

基于开源CFD软件平台OpenFOAM,对E779A螺旋桨在斜流工况下的水动力特性及桨叶上的非定常空化进行模拟分析.分析螺旋桨在非定常空化流中的水动力性能,包括推力系数、桨叶表面压力分布及空泡形态.计算结果表明斜流会对螺旋桨水动力造成干扰,这一影响在来流与螺旋桨首尾线的夹角变大时会增强.对比有无空化发生时的压力分布,表...     

空泡尾部流场特性数值模拟研究

使用商用软件FLUENT6.0,对不同的空泡形态下空泡尾部的阻力特性进行了数值模拟研究.研究表明,随着通气率的增大和空泡形态的增长,运动体尾部阻力系数逐渐减小,并且在减小至谷底后又开始回升.尾部压差阻力系数大于粘滞阻力系数,二者均随着通气率的增加而减小.并且压差阻力系数在达到最小值后开始逐渐增大,因而使得总尾部阻力系数在减小至谷底后开始回升.运动体尾部沾湿长度随着通气率的增大而减小,尾部阻力系数则随着尾部沾湿长度的增大而增大.     

基于均相流输运模型的二维水翼空化模拟

基于均相流输运模型对NACA661-012型水翼进行了空化数值模拟,计算了不同空化数K和不同攻角α下的升力系数与阻力系数.从流场、压力场和边界层的变化分析了空化的产生发展对水翼升力和阻力的影响.当攻角α大于8°以后,空化流动的不稳定性增强,空泡呈现出周期性的生长溃灭,并伴有升力系数和阻力系数的周期性波动.计算结果与实验进行了对比,结果吻合较好.     

三维水翼云空泡非定常特性数值分析及试验验证

  目的  旨在探究三维扭曲水翼空化数值模拟中网格密度及亚格子模型的适应性问题。  方法  为此,使用大涡模拟（LES）方法和Schnerr-Sauer（S-S）空化模型对Delft Twist11N三维扭曲水翼的非定常空化流场进行数值模拟,重点研究3套不同密度的网格和WMLES,SL,WALE这3种亚格子模型对Delft Twist11N水翼空化演变过程、空化脱落频率及时均升阻力系数等的影响。  结果  结果表明：适当的网格加密形式既能够捕捉到较多的细小空泡脱落、马蹄形云空泡的初生与溃灭等非定常空化演变现象,又能够获得具有较高精度的空泡脱落频率、时均升阻力系数和时均压力分布。相较于WMLES和SL亚格子模型, WALE亚格子模型较好地捕捉到了片空泡及云空泡的演变,在预报空泡脱落频率、时均升阻力系数及压力系数方面精度较优。  结论  因此,推荐采用基于WALE亚格子模型的LES方法进行非定常云状空化的数值模拟。     

二维水翼流体力学反问题的数值模拟

水翼(机翼)的反设计一直是人们研究的热门课题,本文对定常反问题进行了研究.将翼型NACA2412作为计算对象.首先给出翼型的下半部形状,以及上表面的从正问题中计算出来的速度分布.在变分原理的基础上,用有限元方法对水翼进行了数值计算.结果表明,本文理论正确,方法可行.     

大攻角水下航行体自然空泡流的数值模拟研究

基于相分数输运方程型的均质平衡流空化模型,采用有限体积法研制了大型空泡流计算程序,对大攻角下运行的水下航行体三维空泡流进行了数值模拟,并与实验结果进行了对比.首次将非线性涡粘湍流模式与基于Rayleigh-Pies-set方程的TEM型空化模型相结合,建立了自然空泡流的数学模型.采用基于SIMPLE的压力-速度-密度耦合修正算法、二阶精度三时间层格式以及基于延迟修正的高阶对流TVD格式.计算模拟了0.2～0.6空化数、4°～20°攻角的不同工况,得到的三维空泡形状及压力分布与实验结果相符.研究了大攻角下航行体周向上的空泡形态分布特征,给出了多种空泡尺度和升阻系数与空化数和攻角之间的关系.通过定量分析发现,空泡的不对称性导致航行体某些部位受力集中,表明高速带空泡运动的航行体在大攻角运动中其结构将受到巨大的水动力载荷.计算还发现,大攻角下的阻力系数与空化数之间的关系和零攻角条件下刚好相反,并根据空泡的不对称性从形状阻力与粘性阻力的关系上对这种现象作出了解释.     

绕水翼非定常空化流动的涡动力学分析

文章采用数值方法对绕水翼的非定常空化流动进行了流动计算,并采用涡动力学方法对非定常空化流场特性进行了分析。计算中,采用基于当地涡旋修正的湍流模型封闭控制方程,通过实验对该计算结果进行验证;采用涡量法分析了空化的发生和发展对涡量输运过程的影响。结果表明:基于当地涡旋修正的湍流模型可以准确地模拟水翼周围非定常空化流场结构;水翼空化流场具有强烈的非定常特性,空化的发展和流场中的涡结构有着紧密的联系;反向射流的作用会导致附着空穴尾缘附近速度梯度的非定常变化,是导致涡量产生和消失的重要因素;空化区域内部水气两相的质量交换会导致流场内体积变化率的变化,使得空化区域内部流体的非正压性显著增强,上述两者导致空化区域内部的涡量分布呈现强烈的非定常特性;空化现象引起的水气两相之间的相互转换以及对湍流的作用都是影响涡结构空间分布的重要因素。     

空化模型在非定常空化流动计算的应用评价与分析

文章基于试验观测数据评价了Kubota与Singhal两种空化模型在非定常云状空化流动数值模拟中的应用。采用商业软件的二次开发技术将两种空化模型引入了计算软件,针对绕Clark-Y水翼的云状空化流动进行数值计算,并与水洞试验结果进行了对比。结果表明,两种空化模型计算的云状空化阶段的流场时均速度分布以及涡量分布具有明显的不同。采用Singhal空化模型可以得到和试验观测更加相近的空化旋涡区与空泡云的旋涡分离的脱落形式。     

Numerical simulation of three-dimensional combustion flows

A finite-rate method is used to simulate the three-dimensional combustion process in a plasma generator with CH4 as the fuel. The simulation was run with RNG k-ε model to simulate turbulence, with eddy-dissipation-concept （EDC） model to simulate the combustion and with discrete ordinates model to simulate radiation. The numerical results show that the flow field characteristics and the parameter distributions are under the condition of rich fuels, and these results provide valuable information when optimizing the plasma generator design and organizing its flow fields.     

三峡下游非恒定流的数值模拟

对三峡坝区下游的黄陵庙至南津关长约28.5 km的河段进行了一维和二维非恒定流数值模拟。文中采用Presimann四点偏心隐式差分求解一维非恒定流,采用正交曲线坐标系下二维水深平均数学模型进行二维非恒定流计算。根据实测资料,由河段上、下游的黄陵庙、南津关(二)水位站测定的水位过程给定水位边界条件,通过河段中的莲沱、石牌、平善坝水位站的实测水位过程以及黄陵庙(陡)水文站实测的流量过程进行模型验证,结果表明一维和二维数模计算出的水位、流量过程与实测值符合良好。由于二维计算时间远大于一维,因此在拟定和优化三峡—葛洲坝联合调峰运行方案时可主要采用一维非恒定流模型,在某些工况下,因通航限制而需要限定局部河段流速时可采用二维非恒定流模型。     

单兵火箭发射筒管内非定常流场数值模拟

由于火箭发动机的运动,使得发射筒内的火箭燃气流场具有运动边界.通过求解二维轴对称N-S方程,对某单兵火箭筒发射筒管内的非稳态流场进行了数值模拟,湍流模型采用带旋流修正的κ-ε模型(Realizab leκ-ε模型),得到了发射筒内的流场分布.     

二维水翼型空化流的数值计算(英文)

In order to predict the effects of cavitation on a hydrofoil, the state equations of the cavitation model were combined with a linear viscous turbulent method for mixed fluids in the computational fluid dynamics (CFD) software FLUENT to simulate steady cavitating flow. At a fixed attack angle, pressure distributions and volume fractions of vapor at different cavitation numbers were simulated, and the results on foil sections agreed well with experimental data. In addition, at the various cavitation numbers, the vapor fractions at different attack angles were also predicted. The vapor region moved towards the front of the airfoil and the length of the cavity grew with increased attack angle. The results show that this method of applying FLUENT to simulate cavitation is reliable.     

Numerical simulation of the flow around two-dimensional partially cavitating hydrofoils

In the present study, a new approach is applied to the cavity prediction for two-dimensional （2D） hydrofoils by the potential based boundary element method （BEM）. The boundary element method is treated with the source and doublet distributions on the panel surface and cavity surface by usethe of the Dirichlet type boundary conditions. An iterative solution approach is used to determine the cavity shape on partially cavitating hydrofoils. In the case of a specified cavitation number and cavity length, the iterative solution method proceeds by addition or subtraction of a displacement thickness on the cavity surface of the hydrofoil. The appropriate cavity shape is obtained by the dynamic boundary condition of the cavity surface and the kinematic boundary condition of the whole foil surface including the cavity. For a given cavitation number the cavity length of the 2D hydrofoil is determined according to the minimum error criterion among different cavity lengths, which satisfies the dynamic boundary condition on the cavity surface. The NACA 16006, NACA 16012 and NACA 16015 hydrofoil sections are investigated for two angles of attack. The results are compared with other potential based boundary element codes, the PCPAN and a commercial CFD code （FLUENT）. Consequently, it has been shown that the results obtained from the two dimensional approach are consistent with those obtained from the others.     

电磁力控制翼型三维绕流场特性的数值研究

在弦长雷诺数Re_L=2.97×10~6下,采用脱体涡模拟方法对弱电解质中电磁力作用下翼型绕流场特性进行了数值模拟,研究了电磁力作用控制翼型失速攻角时绕流场中三维流动特性及失速问题的规律和机理。结果表明:电磁力作用可有效改善翼型周围的流场结构,显著减小翼型绕流场三维特性,并使翼型表面流体动能增加,当电磁力作用足够大时,其表面涡量转变为正向涡量。同时,电磁力作用可增加翼型升力,减小阻力,显著减小升阻力脉动特性,提高翼型升阻比,还可以显著延缓翼型失速特性,增加失速攻角,提高工作性能。     

水中航行体绕流数值计算研究

水动力设计是水中航行体设计的重要组成部分.文章采用数值算法分别对水中航行体全湿和带空泡两种不同状态水下绕流问题进行了计算研究,并与风洞及水洞试验结果进行了对比.计算结果与试验值吻合较好,验证了数值算法的准确性.同时得到了水下绕流的某些运动规律,研究结果对工程设计起到了一定的帮助作用.