求解SUBOFF绕流问题的粘势流耦合方法

针对SUBOFF光体绕流问题,根据Helmholtz速度分解定理将流场速度分解为势流部分和非势流部分,代入到N-S方程推导得到耦合势流速度的关于非势流速度的控制方程,其中势流部分速度采用数值方法预先求得,非势流部分控制方程的表达与N-S方程一致,可直接采用传统N-S方程的计算方法求解。文中在Open FOAM平台内实现了SUBOFF绕流问题的粘势流耦合计算,通过对不同来流速度下的SUBOFF绕流的计算模拟,并将耦合方法计算获得的流场中的速度、压力等,与文献中试验结果以及传统计算流体动力学(CFD)等计算的结果进行仔细比较和分析,相关结果吻合良好。流场中的非势流部分速度分量在远场较小,甚至趋于0,论文分析考察了较传统CFD方法计算范围小的区域,将粘势流耦合方法与传统CFD方法进行了计算研究和对比。研究表明,粘势流耦合方法对区域大小不敏感,可以采用更小的计算区域,同时提高计算效率;该方法适用于绕流问题的计算模拟,为改善CFD计算效率提供了新的途径。     

潜艇的三维厚边界层研究

建立了三维湍流边界层及尾流的动量积分方程,通过算例将计算结果与Gadd的结果进行比较和分析,验证了该方法的可靠性.在潜艇的粘性边界层计算过程中,应用了在势流理论基础上建立起来的艇体表面流线坐标系,对潜艇周围的三维厚边界层进行研究,得到了潜艇边界层及尾流的特征参数,给出了附体(围壳和尾翼)对特征参数的影响.在考虑排挤影响后,艇体采用等价源方法而尾流采用等价假想体法,通过势流计算得到潜艇表面的压力分布.     

浅水中下潜物体有航速辐射问题的数值解

开发了一种基于非均匀有理B样条(NURBS)的Rankine源高阶面元法用于求解浅水中下潜物体有航速辐射问题。基于势流理论,流场中的扰动速度势用物体表面和自由面上分布的Rankine源来表达,同时用映像法计及水底的影响;采用NURBS曲面精确表达物面和自由面,在求出边界面上的未知源强后,物体表面的速度势分布用B样条表示;采用配置方法对边界积分方程进行求解,采用高斯—勒让德公式计算方程中的积分,同时开发了一种解决数值积分中奇异性问题的方法。采用文中数值方法求解了浅水中下潜圆球的有航速辐射问题,计算得到了下潜圆球水动力系数和自由面波形,所得水动力系数计算结果和他人的计算结果进行了比较,其吻合程度良好,从而验证了文中开发方法的有效性。     

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

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

汽车排气消声器声学性能仿真研究

基于三维声波波动方程的有限元软件LMS Virtual．lab Acoustics,采用两种不同的研究方法获得了消声器的传递损失,第1种方法是采用定义入口单位质点振速的方法;第2种方法是采用AML（ Automatically Matched Layer ）技术直接计算声功率的方法．结果表明两种方法计算出来的结果高度一致,且第2种方法计算步骤简单、计算速度快,对于大截面管路消声器也同样适用．相对传统的传递矩阵法和依靠经验的方法,用三维声学有限元法,在计算速度上有很大的提高,并且有效地减少了基于一维平面波理论所产生的消声器传递损失计算误差．最后利用第2种方法,研究了基于气流速度影响的消声器出口辐射噪声,并建立了其插入损失的有限元计算模型．     

增压器涡轮叶片间流体流动的数值模拟

从N-S基本方程出发,采用亚格子尺度模型,运用大涡模拟方法,通过建模和数值计算研究涡轮叶片间流体流动时周围流场的速度、压力、雷诺数和流线分布情况,得出涡轮叶片间流场参数的分布范围,发现各个工况点流场的变化规律和影响因素.该研究结果对增压器涡轮结构进一步优化设计具有重要的指导意义.     

声辐射边界元的数值计算

采用有限元与边界元方法研究结构的声辐射声场,利用有限元方法求解结构的表面速度分布,作为Helmholtz边界积分方程的边界条件,计算结构辐射的声场。边界积分方程中奇异积分的计算采用非等参单元变换,并采用CHIEF方法对Helmholtz方程的非唯一处理,并通过脉动球源与辐射立方体对仿真计算进行验证,结果表明,采用非等参单元与CHIEF方法计算的结果与理论的结果有良好的一致性。     

基于高阶面元法的浅水域中船-船水动力相互作用数值预报

文章开发了一种基于非均匀有理B样条（NURBS）的三维高阶面元法对浅水域中两船会遇和追越时的船-船水动力相互作用进行预报。该方法采用NURBS精确表达船体曲面,采用配置点法满足离散的边界积分方程,在边界面上的源强分布确定之后采用B样条表达物面速度势分布。基于低速假设,忽略了自由面的兴波效应,采用无穷镜像法处理刚性自由面和浅水效应,采用时间步进法在时域内求解速度势。将文中计算结果与基于RANS方程求解的CFD方法和细长体理论方法计算结果比较,验证了该方法的有效性。在此基础上,分别在不同水深、船间横向间距和船速下进行了系列的计算,分析了这些因素对船-船水动力相互作用影响。结果表明,在相同的工况下,船-船会遇时相互作用力比追越时更大,对水深变化更敏感,而对横向间距不如追越时敏感。     

圆球势流绕流的有限元方法

研究利用有限元方法解钝体绕流,利用有限元方法对钝体绕流进行数值编程求解,为简化计算,只求解理想不可压无旋绕流问题,以圆球绕流为研究对象,进行轴对称流动等参数有限元分析,对流场进行不同密度网格划分,比较不同网格密度下圆球绕流场的流线和压力场分布情况。计算结果表明,有限元法可以很好地得出圆球势流场的物理量分布,不同的网格密度对有限元法的求解会产生明显的影响,网格越密求解的流场精度越高,越接近真实的流动情况。     

基于高阶面元法的浅水域中船—船水动力相互作用数值预报（英文）

文章开发了一种基于非均匀有理B样条(NURBS)的三维高阶面元法对浅水域中两船会遇和追越时的船-船水动力相互作用进行预报。该方法采用NURBS精确表达船体曲面,采用配置点法满足离散的边界积分方程,在边界面上的源强分布确定之后采用B样条表达物面速度势分布。基于低速假设,忽略了自由面的兴波效应,采用无穷镜像法处理刚性自由面和浅水效应,采用时间步进法在时域内求解速度势。将文中计算结果与基于RANS方程求解的CFD方法和细长体理论方法计算结果比较,验证了该方法的有效性。在此基础上,分别在不同水深、船间横向间距和船速下进行了系列的计算,分析了这些因素对船-船水动力相互作用影响。结果表明,在相同的工况下,船-船会遇时相互作用力比追越时更大,对水深变化更敏感,而对横向间距不如追越时敏感。