CN113128096A - 一种获取水下航行器直航附加质量方法 - Google Patents

一种获取水下航行器直航附加质量方法 Download PDF

Info

Publication number
CN113128096A
CN113128096A CN202110464668.1A CN202110464668A CN113128096A CN 113128096 A CN113128096 A CN 113128096A CN 202110464668 A CN202110464668 A CN 202110464668A CN 113128096 A CN113128096 A CN 113128096A
Authority
CN
China
Prior art keywords
grid
auv
additional mass
grids
simulation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110464668.1A
Other languages
English (en)
Other versions
CN113128096B (zh
Inventor
郑志航
何波
沈钺
刘继鑫
严天宏
申洪彬
徐玉雪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao Pengpai Ocean Exploration Technology Co ltd
Original Assignee
Qingdao Pengpai Ocean Exploration Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qingdao Pengpai Ocean Exploration Technology Co ltd filed Critical Qingdao Pengpai Ocean Exploration Technology Co ltd
Priority to CN202110464668.1A priority Critical patent/CN113128096B/zh
Publication of CN113128096A publication Critical patent/CN113128096A/zh
Application granted granted Critical
Publication of CN113128096B publication Critical patent/CN113128096B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出一种获取水下航行器直航附加质量方法，包括：A、AUV模型导入；B、网格划分；C、网格独立性验证，获得较优网络；D、进行时间步长独立性验证，获得较优的瞬态仿真时间步长；E、网格导入；F、将UDF程序编译并导入fluent中进行瞬态仿真。本发明与传统的基于动网格的CFD方法不同，该方法不需要使用动网格，通过改变流体速度模拟AUV变速的过程，监控阻力值并通过计算求解附加质量。本发明在保证计算精度的同时降低了网格要求且极大的提升了计算速度。为验证本发明准确性，对标准椭圆球分别进行CFD仿真及理论计算。本发明得到的附加质量值与理论计算值偏差小于3％，说明在设计初期，用该方法计算任意外形的AUV直航方向的附加质量结果可靠。

Description

一种获取水下航行器直航附加质量方法
技术领域
本发明属于船舶与海洋工程技术领域，具体涉及一种水下航行器直航附件质量的求取方法。
背景技术
水下自主航行器(Autonomous Underwater Vehicle，AUV)是一种常用于海洋探测、监察的非标设备。AUV会根据用户使用需求而对壳体进行定制，因此AUV外形不一。为了更好地操纵AUV，就必须先了解其在水中的操纵性。附加质量是物体在理想流体中作非定常运动是所受到的水动力，是潜艇操纵性的其中一个参数。
获得非标AUV附加质量的方法有EFD(Empirical Fluid Dynamics，经验流体力学),CFD(Computational Fluid Dynamics，计算流体动力学)、Hess-Smith面元法和水下实验。其中水下实验是最准确的方法，但成本高、耗时长。而EFD只能计算标准形状的附加质量。Hess-Smith面元法计算附加质量的精度满足工程应用需求，但程序代码复杂，较难快速上手。
ANSYS是成熟的流体CFD仿真软件，能用于计算附加质量，常见的利用ANSYSFLUENT计算附加质量使用了动网格技术，分别模拟物体在水中做规律的正弦运动和水平匀加速直线运动，并通过分析仿真所得数据求得直航方向附加质量。由于需要使用动网格技术，网格质量对仿真十分关键。使用动网格仿真，势必会造成部分网格移动重建，且为保证CFD仿真精度，往往越复杂外形的物体表面网格尺寸越小，小尺寸网格限制了动网格仿真的单位时间步长，从而降低了求解速度而无法快速求解附加质量。
发明内容
本发明提出一种基于ANSYS的CFD求解外形不一的AUV直航方向的附加质量的方法，与传统的基于动网格的CFD方法不同，该方法不需要使用动网格，在保证计算精度的同时降低了网格要求且极大的提升了计算速度，其采用如下技术方案予以实现：
一种获取水下航行器直航附加质量方法，包括：
A、AUV模型导入；
B、网格划分；
C、网格独立性验证，获得较优网络；
D、进行时间步长独立性验证，获得较优的瞬态仿真时间步长；
E、网格导入；
F、将UDF程序编译并导入fluent中进行瞬态仿真
AUV做匀速直线运动时，AUV的推进力Fp大小等于匀速运动时的阻力Drag；而当AUV做加速恒为a的水平直线运动时，将匀加速运动过程离散化，选取vt前后时间段的推力值，使用最小二乘法拟合t时刻前后一段时间的数值，来推算某一时刻的推力Fvt
Figure BDA0003043106750000021
结合公式Fvt-Dragvt-ma＝m11·a得出：
Figure BDA0003043106750000022
其中ε为修正系数。
进一步地，所述步骤B包括：创建外流域和内流域，内流域划分成非结构网格，外流域为划分成结构网格，同时在AUV表面添加边界层网格。
进一步地，基于步骤B中的网格，生成一系列连续细化的网格，网格数量变化比率为
Figure BDA0003043106750000023
与现有技术相比，本发明的优点和积极效果在于：
本发明与传统的基于动网格的CFD方法不同，该方法不需要使用动网格，通过改变流体速度模拟AUV变速的过程，监控阻力值并通过计算求解附加质量。本发明在保证计算精度的同时降低了网格要求且极大的提升了计算速度。为验证本发明准确性，对标准椭圆球分别进行CFD仿真及理论计算。本发明得到的附加质量值与理论计算值偏差小于3％，说明在设计初期，用该方法计算任意外形的AUV直航方向的附加质量结果可靠。
附图说明
图1为本发明实施例获取水下航行器直航附加质量方法流程图；
图2为本发明实施例AUV坐标系；
图3为本发明实施例椭圆球模型；
图4为本发明实施例流域设置；；
图5为本发明实施例网格整体图及截面图；
图6为本发明实施例SUBOFF模型；
图7为本发明实施例阻力与时间关系图一；
图8为本发明实施例阻力与时间关系图二；
图9为本发明实施例阻力与时间关系图三。
具体实施方式
ANSYS FLUENT是成熟的流体仿真软件，能准确地对不同形状的物体在不同的工况下求解RANS方程。Mochammad、RAZA、姚保太等人利用动网格，分别模拟椭圆球在水中做规律的正弦运动和水平匀加速直线运动，并通过分析仿真所得数据，准确地求得直航方向附加质量。由于需要使用动网格技术，网格质量对仿真十分关键。仿真过程中部分网格会移动重建，所以动网格仿真的最大步长取决于最小网格尺寸，对于复杂外形的物体便无法快速求解附加质量。本实施例提出一种无需使用动网格快速求解直航方向附加质量的CFD求解方法。
为了能够更加清楚地理解本发明的上述目的、特征和优点，下面结合附图及实施例对本发明做进一步说明。
参考图1，本实施例基于ANSYS FLUENT的CFD求解外形不一的AUV直航方向的附加质量的过程如下：
步骤A、导入AUV三维模型于ANSYS ICEM CFD软件中
当然也可选择将三维模型导入至其他网格划分软件中，如GAMBIT、MESHING等。
步骤B、网格划分
实现过程：创建外流域和内流域，内流域划分成非结构网格(非结构网格优点：适应性极强，能较好地处理复杂平面。)外流域为划分成结构网格(结构网格的优点在于：占用计算资源少；计算速度快)。同时在水下航行器AUV表面添加边界层网格(边界层是高雷诺数绕流中紧贴物面的粘性力不可忽略的流动薄层。流体流经固体壁面时会在紧贴壁面的区域内形成边界层，边界层内流体具有很大的速度梯度。因此本实施例对水下机器人进行流体分析时考虑边界层，添加边界层网格可以在保证计算精度的同时极大地降低网格数量以提升仿真运算速度和减低计算机资源使用。)
作为替代方案，也可只创建一个流域，将流域内网格全部划分为结构网格或全部划分为非结构网格。
步骤C、进行网格独立性验证，得到较优网格
一般而言，网格数量越大，运算精度越准，但运算时间越久，本实施例选取的较优网格就是网格数量尽量小的同时，运算精度也比较好。为合理地使用计算机资源，较快地完成瞬态仿真，本实施例基于步骤B初始网格，生成一系列连续细化的网格，网格数量变化比率为
Figure BDA0003043106750000031
通过流体仿真选择网格数量较少且仿真数据变化不大的网格。如表2所示，MESH2,MESH 3,MESH 4仿真计算结果相似，所以选取数量最少的MESH2为较优网格。
步骤D、进行时间步长独立性验证，得到较优的瞬态仿真时间步长
与较优网格相似，本实施例选取的较优时间步长就是时间步长尽量大的同时，运算精度也比较好。
步骤E、选定一个直航速度
将较优网格导入ANSYS FLUENT中，并对其进行直航方向1m/s匀速直线运动的流体仿真，并储存数据Drag1m/s，水下航行器航速设为0-2.5m/s本实施例优选1m/s。
步骤F、将UDF程序(本实施例优选加速度a为0.1m/s2为示例)编译并导入fluent中进行瞬态仿真
UDF规定：0到2.4s入口速度为0.8m/s，2.4s后入口速度为[0.8+a·(t-2.4)]m/s，运算时间t设置为0-5s，并记录每个时间步长的阻力(推力)数据。本实施例优选加速度a为0.1m/s2，可替换加速度a数值大小，但加速度不能过高；初始速度0.8m/s可替换，速度大小尽量接近第五步的直航速度；初始速度仿真时间2.4s可替换，考虑到瞬态仿真初期仿真会产生较大的偏差，先使用恒定来流速度仿真至匀速直线运动阻力基本保持不变后再进入匀加速直线运动，要求仿真时间后在瞬态匀速仿真基本稳定。
在水下直航时，AUV只受到由加速度引起的惯性力和速度引起的粘性力。当AUV做匀速直线运动时，AUV的推进力FP大小如下所示，其中Drag为匀速运动时的阻力。
Fp＝Drag
而当AUV做加速恒为a的水平直线运动时，在t时刻AUV的瞬时推力大小等于在t时刻由速度产生的阻力大小加上在匀加速度a下的AUV质量与附加质量所产生的惯性力，关系式如下：
Fpt-Dragt-ma＝m11·a
将匀加速运动过程离散化，则匀加速运动近似为AUV在时间差Δt内做匀速直线运动，并将上一迭代的结果带入下一迭代的开始。利用CFD瞬态计算，通过监控AUV的阻力实现记录每个时间差的阻力Dragt-nΔt+i。并通过最小二乘法拟合t时刻前后一段时间的数值，来推算某一时刻的推力Fpt，以保证拟合而得的推力Fpt为t时刻的瞬时推力。
Figure BDA0003043106750000041
将如上三式结合，便可得：
Figure BDA0003043106750000042
其中ε为修正系数。
为验证本发明提出的方法能准确计算水下航行器(AUV坐标系参考图2)的附加质量，选取3个简单的、可进行理论计算的椭圆球模型，如图3所示，L为椭圆球长度，D为椭圆球截面直径，M为椭圆球的质量。
为降低尾流效应和增加水动力性能，AUV的长径比一般为10。实际生产的AUV长径比一般为7-10。本实施例选取长径比为5、7、10的3个椭圆模型进行仿真验证。模型的具体几何参数如表1所示。
表1椭圆球模型参数
Figure BDA0003043106750000051
一个任意形状的物体运动时共有36个附加质量，可列成方阵如下：
Figure BDA0003043106750000052
由势流理论式
Figure BDA0003043106750000053
(i,j＝1，2，…，6)，可得mij＝mji，式中s为物体表面积；n为微元面积ds的外法线方向；
Figure BDA0003043106750000054
分别是物体延坐标轴x、y、z以单位速度平移运动所引起的流体速度势；
Figure BDA0003043106750000055
分别是物体以单位角速度绕x、y、z轴作转动运动所引起的流体速度势。由于AUV基于XOZ面对称，且各附体较小，对附加质量影响较小，可近似看作AUV基于XOY对称，所以方阵可简化为如下：
Figure BDA0003043106750000056
方阵中除了m11和m44，其他均可通过平面运动机构试验模拟(Planar MotionMechanism，PMM)纯升沉、纯横荡、纯俯仰、纯艏摇直接或间接求得。
Blevins给出了椭圆球的附加质量m11的经验公式：
Figure BDA0003043106750000057
式中ρ是水的密度，α是由Blevins测量的经验参数，由椭圆的长度与直径之比确定，α值如表1所示，L为长轴长度，D为短轴长度。由公式(5)计算得椭圆1，椭圆2和椭圆3的m11理论值如表1所示。
仿真模型的流域尺寸如图4所示：速度入口边界距离椭圆前端13倍椭圆最大截面直径，自由出口距离椭圆尾端33倍椭圆最大截面直径。外流域采用占用计算资源少的结构网格，内流域采用非结构网格。
如图5所示，边界条件设置如下：入口设置为速度入口(velocity-inlet)，出口设置为自由流(outflow)，wall设置为对称面(symmetry)，将流域交界面interface将结构网格与非结构网格的节点融合并将类型设置为interior。
对不可压湍流模型仿真，使用FLUENT进行有限体积法求解雷诺平均方程，采用有限体积法对二阶迎风格式的控制方程进行离散。压力-速度方程采用半隐式方程组(SIMPLE)，湍流模型选择自由剪切流动和壁面限制边界层流动优势明显的SST k-omega，同时需要在网格模型中添加边界层。边界层的首层高度为:
Figure BDA0003043106750000061
其中y+是无量纲，υ为流体运动粘度；u*为近壁面摩擦系数，
Figure BDA0003043106750000062
τw为壁面剪切应力，本实施例取y+＝30。
在CFD分析之前，需网格独立性试验，选取合适的网格文件进行后续分析。以Model3为对象，来流速度为1m/s，网格独立性试验如表2所示。
表2网格独立性验证
Figure BDA0003043106750000063
本实施例选取MESH2作为后续Model3的仿真对象，同时按照相同参数绘制Model1和Model2的网格文件用于后续仿真。
为验证本实施例直航仿真的准确性，本实施例对ITTC公布的SUBOFF标准模型进行不同航速下的直航仿真，SUBOFF标准模型主要几何参数如表3所示；
表3 SUBOFF物理参数
Figure BDA0003043106750000064
由于SUBOFF模型与椭圆3的外形尺寸相似，按照MESH2相同网格参数对SUBOFF标准模型进行三维结构网格划分，SUBOFF模型如图6所示：
通过CFD计算得出的阻力值与ITTC实验值进行对比，如表4所示，仿真值与实验值偏差不超过5.0％，可以验证本实施例的仿真方法可靠。
表4仿真准确性验证
Figure BDA0003043106750000065
Figure BDA0003043106750000071
将Model1，2，3的网格分别导入fluent中求解，选用稳态求解器，求解3个模型分别在水下1m/s直航时所受阻力Drag1m/s。记录后选用瞬态求解器，导入UDF，UDF规定：0到2.4s入口速度为0.8m/s，2.4s后入口速度为[0.8+a·(t-2.4)]m/s。运算2.4s后停止运算。以Model1为模型，加速度为0.1m/s2为例，监控瞬态每个时间步长的阻力，阻力与时间关系图如图7-图9所示。
代入式
Figure BDA0003043106750000072
求解。同理将Model1，Model2，Model3分别以加速度0.1m/s2，0.2m/s2，0.4m/s2进行数值计算及数值拟合，结果如表5所示。
表5仿真结果
Figure BDA0003043106750000073
本发明得到的附加质量值与理论计算值偏差小于3％，说明在设计初期，用该方法计算任意外形的AUV直航方向的附加质量结果可靠。
以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例应用于其它领域，但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

Claims (3)

1.一种获取水下航行器直航附加质量方法，其特征在于包括：
A、AUV模型导入；
B、网格划分；
C、网格独立性验证，获得较优网络；
D、进行时间步长独立性验证，获得较优的瞬态仿真时间步长；
E、网格导入；
F、将UDF程序编译并导入fluent中进行瞬态仿真
AUV做匀速直线运动时，AUV的推进力Fp大小等于匀速运动时的阻力Drag；而当AUV做加速恒为a的水平直线运动时，将匀加速运动过程离散化，选取vt前后时间段的推力值，使用最小二乘法拟合t时刻前后一段时间的数值，来推算某一时刻的推力Fvt
Figure FDA0003043106740000011
结合公式Fvt-Dragvt-ma＝m11·a得出：
Figure FDA0003043106740000012
其中ε为修正系数。
2.根据权利要求1所述的获取水下航行器直航附加质量方法，其特征在：所述步骤B包括：创建外流域和内流域，内流域划分成非结构网格，外流域为划分成结构网格，同时在AUV表面添加边界层网格。
3.根据权利要求2所述的获取水下航行器直航附加质量方法，其特征在：基于步骤B中的网格，生成一系列连续细化的网格，网格数量变化比率为
Figure FDA0003043106740000013
CN202110464668.1A 2021-04-28 2021-04-28 一种获取水下航行器直航附加质量方法 Active CN113128096B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110464668.1A CN113128096B (zh) 2021-04-28 2021-04-28 一种获取水下航行器直航附加质量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110464668.1A CN113128096B (zh) 2021-04-28 2021-04-28 一种获取水下航行器直航附加质量方法

Publications (2)

Publication Number Publication Date
CN113128096A true CN113128096A (zh) 2021-07-16
CN113128096B CN113128096B (zh) 2022-09-27

Family

ID=76780970

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110464668.1A Active CN113128096B (zh) 2021-04-28 2021-04-28 一种获取水下航行器直航附加质量方法

Country Status (1)

Country Link
CN (1) CN113128096B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114662209A (zh) * 2021-11-24 2022-06-24 浙江大学 一种碟形水下直升机外形的cfd仿真优化及验证方法
CN116127611A (zh) * 2023-04-13 2023-05-16 中国人民解放军国防科技大学 一种水下航行器动态仿真方法
CN114662209B (zh) * 2021-11-24 2024-07-26 浙江大学 一种碟形水下直升机外形的cfd仿真优化及验证方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109000838A (zh) * 2018-05-25 2018-12-14 哈尔滨工程大学 一种适用于全海深auv的浮力测量方法
CN109901403A (zh) * 2019-04-08 2019-06-18 哈尔滨工程大学 一种自主水下机器人神经网络s面控制方法
CN110298134A (zh) * 2019-07-05 2019-10-01 大连海事大学 提高水下机器人自航对接瞬态运动预报的数值方法
CN111735455A (zh) * 2020-06-02 2020-10-02 江苏帝一集团有限公司 基于改进的高斯距离迭代算法对接回收组合导航方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109000838A (zh) * 2018-05-25 2018-12-14 哈尔滨工程大学 一种适用于全海深auv的浮力测量方法
CN109901403A (zh) * 2019-04-08 2019-06-18 哈尔滨工程大学 一种自主水下机器人神经网络s面控制方法
CN110298134A (zh) * 2019-07-05 2019-10-01 大连海事大学 提高水下机器人自航对接瞬态运动预报的数值方法
CN111735455A (zh) * 2020-06-02 2020-10-02 江苏帝一集团有限公司 基于改进的高斯距离迭代算法对接回收组合导航方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHEN, CHEN-WEI;YAN, NING-MIN: "Prediction of Added Mass for an Autonomous Underwater Vehicle Moving Near Sea Bottom Using Panel Method", 《2017 4TH INTERNATIONAL CONFERENCE ON INFORMATION SCIENCE AND CONTROL ENGINEERING (ICISCE)》 *
CHEN-WEI CHEN; YONG JIANG; CHEN-HUI HUANG; JIAN-XING LENG; ZHOU: "Computational fluid dynamics study on magnus force of an autonomous underwater helicopter", 《OCEANS 2017 - ANCHORAGE》 *
叶作霖: "AUV自航下潜运动的类物理数值模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
王伟: "微小型AUV浮力调节系统设计及定深控制研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 *
王雪梅: "水下机器人水动力学系数计算与操纵性能研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑》 *
艾晓锋: "基于动网格技术的AUV自航数值模拟研究", 《中国优秀硕士学位论文全文数据库》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114662209A (zh) * 2021-11-24 2022-06-24 浙江大学 一种碟形水下直升机外形的cfd仿真优化及验证方法
CN114662209B (zh) * 2021-11-24 2024-07-26 浙江大学 一种碟形水下直升机外形的cfd仿真优化及验证方法
CN116127611A (zh) * 2023-04-13 2023-05-16 中国人民解放军国防科技大学 一种水下航行器动态仿真方法

Also Published As

Publication number Publication date
CN113128096B (zh) 2022-09-27

Similar Documents

Publication Publication Date Title
Pendar et al. Cavitation characteristics around a sphere: An LES investigation
CN104298869B (zh) 一种弹性水翼的流固耦合特性数值预测方法
CN109948301B (zh) 基于网格控制的近水面滑行跳跃流固耦合数值预测方法
Usta et al. A study for cavitating flow analysis using DES model
Liefvendahl et al. Computation of cycle-to-cycle variation in blade load for a submarine propeller using LES
CN112733473A (zh) 一种基于cfd的船舶自由横摇衰减数值模拟自动化智能化方法
CN113128096B (zh) 一种获取水下航行器直航附加质量方法
Noblesse et al. Ship bow waves
Yilmaz et al. An improved tip vortex cavitation model for propeller-rudder interaction
Rhee et al. RANS model for spilling breaking waves
Cummings et al. Supersonic, turbulent flow computation and drag optimization for axisymmetric afterbodies
Haipeng et al. Experimental and numerical research on cavitating flows around axisymmetric bodies
Neel Advances in computational fluid dynamics: turbulent separated flows and transonic potential flows
Bašić et al. Analytical and numerical computation of added mass in ship vibration analysis
Rockwood et al. Real-time identification of vortex shedding in the wake of a circular cylinder
Li et al. Research on the improved body-force method based on viscous flow
CN116227073A (zh) 一种研究水下发射航行体肩空泡发展溃灭的方法
Saisto Specialist committee on CFD in marine hydrodynamics
CN108536954A (zh) 一种基于交点间断伽辽金的高精度格子波尔兹曼方法
Robinson Application of vorticity confinement to inviscid missile force and moment prediction
Sukas et al. HEMLAB algorithm applied to the high-lift JAXA standard model
Killingstad A study of dead water resistance Reynolds Averaged Navier Stokes simulations of a barge moving in stratified waters
Hu et al. The structures of unsteady cavitation shedding flow around an axisymmetric body with a blunt headform
KIRIKBAŞ et al. Computation of Added Mass Coefficients of Darpa Suboff Using Two Different Solvers of Openfoam
Kirikbas et al. Prediction of horizontal and vertical turning maneuvers of a generic submarine model by closed-form solutions and CFD

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant