导管螺旋桨敞水性能的预报和比较

借助于计算流体力学软件,对导管螺旋桨的敞水性能进行了数值模拟,得到了在不同网格模型和湍流模型下导管螺旋桨的正车敞水性能曲线,通过与试验图谱的对比分析,发现Wilcoxk-ω模型更适合于导管螺旋桨的敞水性能计算,并且采用结构化网格和非结构化网格相结合的计算方法能够满足导管螺旋桨敞水性能预报的工程精度要求。同时,当网格中存在少量高度倾斜的网格单元时,在不影响计算收敛的情况下,仍能将计算误差控制在10%以内。此外,如果要获得更为精确的计算结果,应提高网格质量,尽量使用结构化网格,并将近壁面网格加密,合理控制壁面附近的Y+值。虽然网格数量的增加并不总是意味着计算误差的减少,但合理控制网格细密度能够获得更为可信的计算结果。     

CFD在高速滑行表面性能预报中的应用

针对滑行表面数值预报结果中船底出现异常水气分布的问题,基于三种常用CFD软件FLUENT、CFX及STAR-CCM+研究了VOF界面插值格式、网格类型及数目等因素对滑行表面性能数值预报的影响,并提出相应改善数值计算模拟的具体措施。研究表明:基于相同的结构化网格及计算设置时,采用几何重构界面插值格式预报滑行块的阻力及底部水气分布上最为准确,CICSAM格式次之,修正HRIC格式最弱。通过选择合理的网格类型及网格数目也可以改善滑行船体性能预报精度。CFX采用四面体网格,STAR-CCM+采用多面体网格可以较为准确预报滑行艇体的流场和阻力。建议在滑行艇的阻力及流场预报中采用几何重构的VOF界面插值方式,选择合理的网格类型并进行适当的网格加密可以提高水气界面捕捉的精度,减小阻力预报的误差。     

考虑航行姿态的船模阻力及流场数值预报

为精确得到船模阻力及流场,研究了网格划分、湍流模型及船模姿态对预报结果的影响。首先,采用切割体网格技术及棱柱层网格技术划分了多套网格,研究网格尺寸大小及棱柱层参数对阻力、兴波及伴流场的影响;然后,比较了三种不同湍流模型在船模阻力及伴流场预报中的差异;最后,基于合理的网格划分及湍流模型,采用DFBI的方法预报了船模在不同航速下的姿态及阻力。计算结果表明:当Fr0.25时,不考虑航行姿态时阻力计算误差可达到18.3%。考虑航行姿态变化后,阻力计算误差小于3%;升沉和纵倾的计算值和试验值吻合较好。由此证明:采用DFBI方法考虑船体航行姿态的变化避免了网格变形或者重构,具有较好的收敛性和计算精度。     

基于FLUENT的调距桨敞水性能数值计算方法探究

以某调距桨为研究对象,基于FLUENT数值计算软件,对螺旋桨转速和网格Y+值对敞水性能计算结果的影响进行了探究,并对调距桨不同螺距的敞水曲线进行了计算。结果显示,在J=0.6～1.0时,转速越小,计算值与试验值误差越小,而在大进速系数和小进速系数下,高转速计算结果相对低转速计算结果更趋近于试验结果;网格Y+在160～400范围内,所得敞水性能与试验值符合较好;相同的网格不能满足调距桨不同螺距的计算精度,对于调距桨不同螺距的计算,应该分别进行探讨,从而确定最佳的计算方案。上述结果表明:所提方法对调距桨敞水性能计算中的网络和选取工作有一定的借鉴意义。     

潜艇上浮运动数值计算方法的研究与验证

预报潜艇上浮运动是探索研究潜艇上浮运动控制规律以及潜艇安全上浮的前提。就潜艇上浮运动的数值计算方法和计算流程进行了详细介绍,对数值计算方法中的湍流模型及离散格式等进行了优选。基于模型上浮试验结果,对数值计算的不确定度进行了分析且对数值计算结果进行了验证。证明了该数值计算的可靠性,并发现数值计算结果与试验结果很接近,验证了该数值计算方法的可行性,说明该数值计算方法和试验可为潜艇上浮运动的研究提供参考和借鉴。     

基于STAR-CCM+的全附体潜艇尾流场数值分析

为充分了解潜艇周围的粘性流场特征,基于STAR-CCM+仿真软件对SUBOFF全附体模型的尾流场进行了数值模拟。通过将数值模拟结果与试验结果进行比较,系统分析了棱柱层网格的划分与湍流模型的选取对计算结果的影响。结果表明:SSTk-ω模型较适合应用在潜艇模型尾流场计算中,并且棱柱层的总厚度不能小于15mm,第一层计算网格节点到表面的无因次距离y+应在25到170之间。在此基础上,应用滑移网格技术对潜艇尾流场进行数值模拟,通过加密网格和减小时间步长使得轴向伴流分数与整体定常计算值吻合较好。     

高超声速尖双锥流动高精度数值模拟

以25°/55°尖双锥外形的高超声速低焓层流流动模拟为例,对高阶加权紧致非线性格式模拟激波/边界层干扰流动的能力进行验证和确认。空间离散采用二阶MUSCL和三阶、五阶加权紧致非线性格式,时间离散采用二阶精度双时间步方法,通量函数采用混合Roe-Rusanov,AUSMPW+,Van Leer等,对比了不同精度空间离散格式对时间、网格收敛特性和通量函数耗散特性的影响。数值模拟结果表明采用高精度空间离散格式能在较疏的网格上获得收敛解,并能消除结果对通量函数的敏感性,但收敛需要推进更久的计算时间。数值模拟结果与实验测量结果吻合良好,满足工程精度要求。     

基于滑移网格的潜艇旋臂试验数值模拟?

潜艇水动力系数对潜艇操纵性研究具有重要意义。为了获取较高精度的潜艇旋转导数，以分区结构化网格为基础，采用基于固定坐标系的滑移网格方法，对SUBOFF潜艇进行旋臂试验数值仿真。并对比两种不同的湍流模型下仿真结果，误差均在8％以内，满足工程要求，验证了方法的可行性。     

WCNS格式在热化学非平衡流数值模拟中的应用研究

高精度格式因其在较低的计算资源消耗下仍能保持足够的精度而在精细流动结构的模拟中具有独特优势。首次将显式5阶加权紧致非线性格式(WCNS-E-5)引入二维高温非平衡流计算。进行了自由流速度5.0~5.7km/s的高焓风洞和高空条件下的圆柱非平衡流场数值模拟,得到了正确的流场参数分布,壁面压力和壁面热流等结果与实验值吻合较好。计算结果表明WCNS-E-5格式较2阶的MUSCL格式具有更好的热流网格收敛特性。分析并针对算例条件初步解决了高阶格式在刚性化学反应流模拟中的鲁棒性问题,为进一步采用高阶格式开展流场结构更加复杂的热化学非平衡流数值研究奠定了良好基础。     

加权紧致非线性格式在热化学非平衡流数值模拟中的应用

高精度格式因其在较低的计算资源消耗下仍能保持足够的精度而在精细流动结构的模拟中具有独特优势。将显式5阶加权紧致非线性格式(WCNS-E-5)引入二维高温非平衡流计算,进行自由流速度5.0~5.7 km/s的高焓风洞和高空条件下的圆柱非平衡流场数值模拟,得到正确的流场参数分布、壁面压力和壁面热流等结果与实验值吻合较好。计算结果表明,WCNS-E-5格式较2阶的MUSCL格式具有更好的热流网格收敛特性。分析并针对算例条件初步解决了高阶格式在刚性化学反应流模拟中的鲁棒性问题,为进一步采用高阶格式开展流场结构更加复杂的热化学非平衡流数值研究奠定了良好基础。     

测试投稿，大家不要管哈~

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

高超声速钝锥体俯仰阻尼导数影响因素分析

采用非定常双时间步长方法、刚性旋转动网格技术,结合空间隐式格式数值分析了高超声速钝锥体定常与非定常流场特征,并运用最小二乘法辨识出钝锥体的俯仰阻尼导数。通过与相应的实验数据、工程估算值进行比较,分析并验证了壁面网格雷诺数、钝头比、攻角与质心位置对钝锥体俯仰阻尼导数的影响。同时采用五种差分格式,其中包含五阶精度WENO与高精度WNND格式,对钝锥体俯仰阻尼导数进行了计算,结果表明不同无粘差分格式对俯仰静导数影响较小,但不同格式之间俯仰动导数结果存在差异,分析认为差分格式之间不同的粘性分辨率应是差异的主要来源。     

全附体多桨船模自航点数值求解方法研究

针对带附体的四桨双舵船模,采用结构化网格离散全附体模型表面,基于RANS方程求解了船桨舵耦合的复杂粘性流场,数值模拟了主船体阻力试验、螺旋桨敞水试验、全附体模型阻力试验和自航试验,建立了螺旋桨转速自动调节的模型自航点求解方法。该方法的计算结果与试验结果进行了比较分析,结果表明:全附体模型阻力误差小于3.8%,自航点转速误差小于4.3%。     

Spalart-Allmaras湍流模型在高超声速气动加热计算中的应用

将AUSM+格式与LU SGS隐式迭代相结合,采用一方程的Spalart Allmaras湍流模型和二阶迎风MUSCL格式,研究了火星"探路者号"实验模型在10马赫来流条件下的气动热问题。采用量热完全和热完全两种气体模型的计算结果与实验数据进行了比较。热完全气体模型模拟的热流分布规律与实验吻合较好,但是对分离区的热流估计过高。在网格合理划分,尤其是满足近壁面的网格分辨率要求的前提下,该方法能够合理反映高超声速流场除分离区以外的气动热特性。     

采用DES模型对6∶1椭球体进行水动力数值计算

采用基于非结构化网格的有限体积法离散控制方程和湍流模式,应用SIMPLE(semi-implicit method for pressure linked equations)方法处理速度-压力耦合,对流项和非定常项分别采用二阶迎风格式和二阶精度离散,隐式选代求解方程组.基于SA和SST的DES模型,得到了不同攻角下粘...     

跳跃-滑翔弹道扰动引力自适应网格快速赋值方法

为实现高超声速跳跃-滑翔弹道扰动引力的快速赋值,提出自适应网格赋值模型,并根据反距离加权理论,优化广义延拓逼近算法,对模型的逼近误差进行分析。该赋值模型的网格划分为两级,第一级网格根据标准弹道空域进行划分,第二级网格根据滑翔导弹实际弹道在线生成。根据一级网格节点数据,通过优化广义延拓逼近算法计算二级网格节点数据,最后根据二级单元内插计算实际弹道点的扰动引力值。仿真结果表明:在同等大小的网格划分下,优化广义延拓自适应网格模型的逼近精度高于一般赋值方法;在同等精度要求下,该赋值模型的最大单元格边长大于一般赋值方法,从而减少了单元格划分数量,进而降低弹上数据存储量;针对不同滑翔方向以及不同滑翔距离的跳跃-滑翔弹道,该模型逼近误差对应的落点偏差小于5 m,具有较好的适应性。该赋值模型在满足计算速度的前提下,提高了传统赋值方法的逼近精度,降低了弹上存储量,具有一定的工程应用价值。     

高效计算网格壁面距离的KD树方法

计算流体动力学的一些领域中经常要用到流场点到壁面的最小距离,为了提高壁面距离计算的精度和效率,将KD树应用于非结构网格的壁面距离计算,计算了若干三维外形的壁面距离,结果表明:基于KD树的网格壁面距离计算方法鲁棒性好,计算效率和结果精度高,适应复杂外形的能力强;算法的通用性好,可以拓展应用到多种网格类型。     

用多重网格TVD方法求解三维高超音速粘性流场

本文采用时间相关法求解了三维Ｎ－Ｓ方程。在计算中，对于压力和对流项采用Ｒｏｅ的通量差分分裂技术和Ｏｓｈｅｒ－Ｃｈａｋｒａｖａｒｔｈｙ的ＴＶＤ格式，对于粘性项，采用了中心差分格式。通过数值模拟得到了钝头双锥体高超音速层流流场在不同攻角下的数值解。计算的重点放在壁面热流率的预测。在攻角为４°和１２°的状态下，Ｓｔａｎｔｏｎ数的计算值与实验值进行了比较，二者吻合较好。为了加快收敛速率，计算时采用了多重网格方法。     

基于Richardson外推法的CFD离散误差分析

为研究数值模拟中网格类型及网格尺度变化对离散误差产生的影响,对NACA2412三维翼型在0°、4°、8°、12°攻角下进行了计算流体力学(CFD)仿真,得到了对应条件下的升力系数及阻力系数,利用Richardson外推法研究多套连续变化的结构网格及非结构网格数值解的空间离散误差,并结合网格收敛指标(GCI)法对误差进行...     

迎风格式在接触间断的数值耗散及其诱导误差

在不同流场参数条件下对三种迎风格式在接触间断中的数值耗散问题进行了数值实验,并对数值耗散产生的机理进行了分析。数值计算结果和理论分析表明,矢通量分裂格式计算接触间断问题时,若流场静止或流场内存在亚声速区域,密度耗散的产生会诱导出以特征速度运动的数值扰动误差,该误差对数值耗散的大小无影响,但会影响流场的速度及压力参数分布,从而改变流场的结构。在二维问题中,诱导误差相互干扰会产生大量的复杂小尺度结构,给流场结构分析带来困难。同时,在密度参数线性分布的流场中,若空间离散格式重构的对象为对流通量,使用矢通量分裂格式计算流场会产生数值误差,使计算精度难以到达二阶。