基于双向势流—黏流耦合的数值消波模拟

提出一种基于双向势流-黏流耦合的数值消波方法,将计算域分为靠近浮体的黏流计算域和远离浮体的势流计算域,并通过重叠区域进行流场信息匹配。在黏流计算域中求解N-S方程,采用有限体积法对控制方程离散,采用几何VOF方法捕捉自由液面。在势流计算域中求解Laplace方程,引入线性自由液面假设,采用Rankine源法求解。在重叠区域内通过WheelerStretching方法和松弛区域法进行速度场双向信息交换。基于该方法,在OpenFOAM框架下开发了求解器,并在模型尺度下,对不同参数的规则波进行了模拟,还与waves2Foam求解器的计算结果进行对比。计算结果表明,在同等消波效果下,基于论文方法的数值模拟所需的黏流计算域消波区长度为waves2Foam求解器的1/5。     

潜艇水面与水下粘性绕流数值模拟

本文采用求解RANS方程的方法结合四种湍流模型,对于带有不同附体的SUBOFF模型尾流场进行了数值模拟。数值预报的桨盘面处不同半径上的轴向无量纲速度与试验结果进行了对比,结果表明湍流模型在数值模拟中起到重要作用。潜艇水面航行性能十分重要,因而对于潜艇自由液面绕流的数值模拟备受关注。本文采用VOF方法对于两条潜艇模型在不同傅汝德数下的自由液面绕流进行了数值模拟。计算得到的阻力、波形与试验结果吻合较好。文中也探讨了潜艇自由液面绕流的一般特性。并验证了用CFD手段预报潜艇粘性流场的能力。     

与粘流耦合的导管螺旋桨升力面设计方法

针对导管螺旋桨这类组合式推进器的设计问题,采用了粘流CFD与传统的势流升力面设计相耦合的方法来计算,以此来考虑船体及其附体与推进器的相互作用,以及流场粘性的影响,解决设计所需的实效进流场的问题,并采用此方法作了导管螺旋桨的设计和模型试验.     

分离涡模拟在亚临界区风场圆柱绕流中的适用性研究

圆柱绕流问题一直是流体力学领域热门的研究对象,本文基于ANSYS Fluent软件,采用分离涡模拟(DES)的湍流模型,以连续性方程和N-S方程为控制方程,对亚临界区流动状态风场下的圆柱绕流问题进行了数值模拟。主要分析了圆柱阻力系数,Strouhal数,圆柱后绕流速度分布及涡量大小。计算结果表明DES湍流模型可有效描述流场的三维特性,与之前的实验和模拟结果吻合,验证了DES模型在亚临界区用于风场绕流问题的有效性。DES模型较于大涡模拟和雷诺平均法在计算效率和准确性上有一定优势。     

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

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

直接边界元法解势流速度场问题

在船舶水动力学中,大多采用以Hess-Smith方法为基础的间接边界元法求解势流绕流问题,但Hess-Smith方法本质上是基于物理直观提出,在理论和数值计算上都存在着缺点。直接边界元法虽然在船舶水动力领域有着非常广阔的应用前景,但至今应用较少。为推广直接边界元法在船舶水动力学中的应用,根据边界积分法建立积分方程,采用直接边界元法对无界势流绕流问题予以求解,得出流场速度势和物面上的速度分布,并通过与解析解的比较进行误差分析。对二维、三维问题的算例进行数值计算。数值计算过程用Matlab编程实现。结果表明:直接边界元法在求解船舶势流绕流问题中具有足够的精度和较高的效率,且数值计算实现过程更简洁,可发展成为求解船舶兴波等船舶水动力学问题的通用方法。     

基于势流-黏流双向耦合模型的数值造波

采用开源程序OceanWave3D和OpenFOAM提出了一种势流-黏流双向耦合数值造波方法。为提高数值稳定性，在势流域内的重叠区域内建立松弛区域。松弛区域的波高从黏流域直接插值获取，而速度势则根据水面边界条件方程用速度插值间接计算。基于4阶Runge-Kutta算法实现时间推进，并通过线性波、非线性波和不规则波，包括波浪正向、反向传播和二维、三维数值模拟的算例进行验证。计算结果表明，论文的势流-黏流双向耦合模型准确高效，可应用于实际海洋工程波浪的模拟。     

船体绕流场及流噪声的CFD模拟方法

CFD数值模拟技术可以为船舶水下噪声水平和噪声传播提供高精度的预报,同时还可以洞察船体绕流场的变化。本文利用VOF方法与SSTκ-ω两方程湍流模型用于求解船舶非定常粘性流场,并结合FW-H方程进行噪声传播。基于Lighthill声类比理论,对不同球鼻首船型的噪声进行数值计算,对船体流噪声的空间指向性、近远场分布特性进行分析。计算结果表明,CFD技术可以用于模拟分析船舶的绕流场和流致发声问题,能够为低噪声船体线型设计提供参考。     

中等雷诺数方柱绕流的直接数值模拟及涡系分析

作为一种典型的钝体绕流,方柱绕流具有物体几何外形简单而流场结构非常复杂的特性,常规的基于RANS的数值计算方法难以准确模拟。DNS不使用任何湍流模型,直接求解完整的非定常流动控制方程组,模拟包括脉动在内的湍流所有非定常流动量的时空演变过程,是湍流数值模拟中最精确的方法,在复杂流动数值模拟方面的应用潜力巨大。本文自主设计和编制并行数值模拟程序,使用基于"神威·太湖之光"国产CPU架构的大规模并行计算,开展了中等雷诺数(Re=10 000)方柱绕流的直接数值模拟。其中:不可压N-S方程组采用基于交错网格的有限体积法离散;压力-速度耦合采用SIMPLE算法处理;离散得到的代数方程组采用Gauss-Seidel迭代求解;时间步进采用Euler隐式格式,对流项采用QUICK格式,耗散项采用中心差分格式;数值模拟程序的并行化使用MPI方法处理。文中重点分析了方柱绕流的复杂涡系结构,同时给出了部分湍流统计结果,并通过与RANS和LES模拟结果的对比分析,展现了DNS在复杂精细流场模拟方面的优势。     

用雷诺应力模型预报不同雷诺数下的潜艇绕流

文章采用求解RANS方程的方法,结合五种湍流模型(κ-ε、RNG κ-ε、κ-ω、SST κ-ω、RSM),预报了SUB-OFF潜艇模型两种附体状态下(裸艇体+围壳、裸艇体+尾翼)的尾流场,也预报了SubB潜艇模型全附体状态下的尾流场,并与试验结果进行了对比分析,比较了五种湍流模型的计算精度,结果表明雷诺应力模型(RSM)预报的尾流场与试验最为接近.然后,利用CFD方法探讨了雷诺数对SUBOFF潜艇绕流影响的问题.计算了SUBOFF潜艇在不同雷诺数下的绕流,详细分析了尾流场轴向速度波动图、轴向速度等值线云图、边界层厚度和主附体接合部涡旋结构随雷诺数的变化.     

高架索航行补给中船舶在波浪中的运动性能研究

采用三维势流理论计算高架索航行补给中两船在波浪中的运动响应,将船体周围流场速度势分解为波浪入射势、自身绕射势、相邻船体的存在产生的绕射势、自身辐射势和相邻船体运动产生的辐射势,通过源汇分布法进行求解;对高架索索道上的恒定张力在耦合方程中表示为高架索瞬时张力系数矩阵;求解两船体在规则渡中的联立运动方程可得两船的运动响应.通过计算结果与高架索航行补给中两船运动响应的试验结果进行比较,吻合良好,说明文中的计算方法是可行的,在此基础上,对不同高架索恒定张力和不同浪向角条件下的接受船运动幅值响应进行计算并讨论.最后得出结论.     

浸没边界格子Boltzmann方法的改进及转动圆柱绕流模拟

针对浸没边界格子Boltzmann方法计算效率不足的问题,提出一种改进浸没边界格子Boltzmann方法。通过使用作用力格子Boltzmann模型,并简化流固耦合作用力的计算方法,使计算流程得到简化,同时,对比原方法转动圆柱绕流计算结果,计算时间缩短一半以上,提高了数值计算效率。结合改进算法研究转动圆柱绕流的流场升阻力系数、压力系数、流场速度等流体参数场随转速比的关系,揭示了一定的转动圆柱绕流动力学规律。研究结果表明,改进方法准确可靠。     

平板上直立圆柱三维绕流的CFD计算研究

利用计算流体动力学(CFD)方法,对粘性不可压缩流体流过平板上直立圆柱的绕流进行了三维数值模拟,采用了有限体积法和SIMPLE计算程式,利用不可压缩的Navier-Stokes方程,模拟了雷诺数在200的层流绕流流动。在建立贴体的结构化精细网格之后,成功模拟了马蹄涡,与经典的理论对比具有较好的吻合性。同时得到了交接部区域的剪切力云图。计算结果表明,主/附体结构交接部区域流动机理十分复杂,在这个区域有马蹄涡的存在,是产生交接部冲蚀的根本原因。计算结果表明本文所用方法在模拟含有复杂流动机理的圆柱绕流中是可靠性的。     

参数化单元边界元法解势流速度场问题

[目的]边界元法在海洋工程水动力学中有着广阔的应用前景,为推广边界元法在海洋工程水动力学中的应用,[方法]根据边界积分法建立积分方程,采用参数化单元边界元法对势流问题进行求解,得出流场速度势。对经典算例进行数值计算,与数学解析解比较,并进行误差分析。在二维问题下,分别采用非连续参数化单元和参数化单元边界元法求解势流速度场问题;在三维问题下,采用参数化单元边界元法求解势流速度场问题。[结果]结果显示,在二维问题下,采用非连续参数化单元边界元法求解势流问题具有较高的精度和效率,可以在采用较少单元数的情况下得到较为理想的数值解;在三维问题下采用参数化单元边界元法虽然计算速度较快,并可以得到较好的平均相对精度,但有些点误差较大,需要改进算法或使用其他单元进行求解。[结论]参数化单元边界元法在求解海洋工程势流问题时,数值计算实现过程更简洁,可发展成为求解船舶兴波等船舶水动力学问题的通用方法。     

DTMB5415驱逐舰绕流的数值模拟

本文采用RANS方程及RNG湍流模型数值模拟了带方尾及声纳罩的5415驱逐舰绕流，计算结果表明对复杂的实际船形，CFD能有效地预报船体周围流场、波形和阻力等船舶设计感兴趣的水动力性能参数，为船形分析及优化提供了有力的工具。     

非淹没丁坝三维绕流数值模拟

通过应用非线性三维紊流数值模型,并采用有限体积法对模型控制方程进行离散,对非淹没丁坝绕流进行数值模拟仿真分析,给出了三维流场、压力分布场、剪应力分布场等信息,计算结果均基本符合非淹没丁坝绕流的实际情况,结果表明该非线性模型是可取的,数值模拟结果具有一定的可信度.     

三维圆柱体绕流数值模拟流场选择及网格划分

文章通过数值模拟结果确定适合的三维圆柱体绕流流场的尺寸,采用不同的计算流体力学(CFD)软件对有限元模型进行网格划分,对比网格划分质量的优劣,选择适合数值模拟的软件CFX。确定采用的网格形式为六面体网格(Hex8)以及流场不同部分网格的尺寸。圆柱绕流展向各截面平均压力系数沿周向无明显变化,脉动压力系数变化显著,绕流后的旋涡发放形式具有明显的三维特性的数值模拟结果验证了本文的流场选择和网格化分尺度的合理性。     

敞水舵水动力的数值模拟

直接用Euler/N-S方程求解水翼绕流的计算和设计问题,是目前这一领域研究的热点.在大型通用CFD商业软件FLUENT平台下,利用有限体积法和RNG湍流模型直接求解RANS方程对敞水舵流场进行了数值模拟,并计算了敞水舵水动力性能.利用SIMPLE算法求解,采用了二次迎风格式离散动量方程和离散湍流能量方程,所得结果和试验数据、实际情况进行比较,吻合性比较好,对于敞水舵的水动力性能研究和设计具有一定的参考价值.     

基于LES方法圆柱绕流三维数值模拟(英文)

采用计算流体软件CFX5中large.eddy simulation(LES)模型计算了均匀流场中三维圆柱绕流的水动力特性.使用有限体积法对三维N-S方程进行求解.数值模拟着重研究了高雷诺数时展向各截面的压力、阻力、升力及涡管特性.数值计算结果表明:展向各截面柱体受力关于中截面对称且小于二维情况,柱体周围流场呈现明显的三维特性.     

改进的基于特征线的N-S方程算子分裂有限元法

基于特征线的Navier-Stokes（N-S）方程算子分裂有限元法（CBOS有限元法）的核心是在每一个时间层上,采用算子分裂法将N-S方程的对流项和扩散项分开求解;对流项的求解过程借鉴了简单显式特征线时间离散,显式求解。该文在原CBOS有限元法基础上推导了一种更加精确的对流项显式离散方法。通过自编程序对方腔流进行数值模拟,表明该算法具有更高的计算精度。低雷诺数下圆柱绕流计算所得的阻力系数、升力系数、斯特劳哈数等与已有数据较为接近,表明该文中算法能较准确地模拟圆柱绕流的流场特性;最后,文中分析了Re=200时圆柱绕流在一个周期内所受升力变化与对应流场中压力和流线演化的关系。