CN113515805A - 水上飞机波浪水面起降水气两相流数值仿真方法 - Google Patents
水上飞机波浪水面起降水气两相流数值仿真方法 Download PDFInfo
- Publication number
- CN113515805A CN113515805A CN202110371760.3A CN202110371760A CN113515805A CN 113515805 A CN113515805 A CN 113515805A CN 202110371760 A CN202110371760 A CN 202110371760A CN 113515805 A CN113515805 A CN 113515805A
- Authority
- CN
- China
- Prior art keywords
- wave
- grid
- area
- water
- phase
- 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.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000005514 two-phase flow Effects 0.000 title claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 238000004088 simulation Methods 0.000 claims abstract description 10
- 239000007788 liquid Substances 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 238000013016 damping Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 238000011423 initialization method Methods 0.000 description 1
- 238000010327 methods by industry Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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/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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种水上飞机波浪水面起降水气两相流数值仿真方法，该方法基于整体动网格技术和速度入口造波技术，并改善网格划分策略和流场初始化方法，从而实现水上飞机波浪水面起降的非定常流场计算。对于研究水面飞行器波浪着水问题，该方法通过在波面方程、波浪速度场和压强场函数中增加初始相位项和初始平衡位置项，实现波浪在计算域内沿传播方向的相偏移以及波浪平衡位置的垂向偏移，使得波浪在极短的求解时间内与物体接触，省去了物体等待头波的时间。与以往动网格法相比，本发明不仅可以减少网格量、提高网格利用率，而且可以在保证求解精度的同时获得更快的计算速度。
Description
技术领域
本发明涉及飞行器起降试验领域，具体是一种水上飞机波浪水面起降水气两相流数值仿真方法。
背景技术
水上飞机不同于常规飞行器，根据其任务需求，应具备在水面起降的能力，并能够在一定的海况下开展作业。因此，评估水上飞机波浪水面起降能力对于飞行安全至关重要。目前，针对这一个问题的研究手段主要分为：试验研究、理论计算和数值仿真三种方法。试验研究需要耗费大量的实验室搭建成本，试验数据常常受设备条件的限制，试验结果较为单一；理论计算借助经验公式估算飞行器运动过程所受水动力和气动力，多是基于简单的模型和假设，估算结果准确度不高；随着高性能计算机的发展，应用数值仿真的方法处理工程问题越来越普遍，具有求解速度快、模型适应能力强、求解精度高等优点。
在涉及物体间存在相对运动的数值仿真中，动网格技术被大量使用。其中，整体动网格法作为一种高效的用于求解非定常刚体运动的计算方法在本发明中被采用，计算域网格，包括网格单元和外场边界，随物体一起做刚性运动。同时，对各个边界施加体积分数边界条件，使得水面平衡位置始终保持不变。但，该方法一般不适用于物体大角度运动仿真。
速度入口造波技术被广泛应用于数值波浪水池的构建中，把波浪前进和垂向速度的解析解施加在入口边界，使得波面随时间推进并传播到流场域中。同时，在入口边界处设置关于时间的体积分数来体现运动波面。这种造波方法由于是从入口边界开始造波，易于实现且波形可控。
目前，关于飞机波浪着水的数值模拟还比较少，方法多是借鉴静水面着水。然而，在飞机静水面着水时，机体与自由液面的距离只需要调整液面水平位置即可，对于波浪着水，不仅要调整波浪平衡位置，还需要对波浪与机体的水平相对位置即波浪初始相位做出要求，这方面的研究工作尚不完善，也将是本专利着重阐述的内容。
发明内容
本发明为了解决现有技术的问题，提供了一种水上飞机波浪水面起降水气两相流数值仿真方法，基于整体动网格技术和速度入口造波技术，并改善网格划分策略和流场初始化方法，从而实现水上飞机波浪水面起降的非定常流场计算，不仅可以减少网格量、提高网格利用率，而且可以在保证求解精度的同时获得更快的计算速度。
本发明包括以下步骤：
1)计算域网格划分策略：针对复杂的自由液面情况，将整个流场计算域做多区域多等级处理，即把计算域划分为多个网格加密区域，不同加密区域依据当地流场特征采用不同的网格加密尺寸；
2)运用整体动网格法处理物体和计算域网格非定常运动，并将速度入口造波边界条件施加在计算域的入口边界，具体施加方法如下：
2.1)确定自由液面的波面方程：
式中，h为波高，k为波数，ω为圆频率，t为时刻；
2.2)给出波浪在传播方向和垂直方向的速度函数u和w：
式中，d为水深，x和z为空间坐标；
2.3)计算水面以下压力场函数p：
2.4)对各边界施加相体积分数函数，以波面方程为判断准则，流场网格垂向坐标大于波面方程的区域指定为空气相，反之，则为水相；
2.5)对速度入口边界施加波浪速度函数u和w，对压力出口边界施加压力场函数p，至此，边界条件的函数施加得以完成；
3)流场域初始化：将关于空间坐标的初始速度场函数和压力场函数施加在流场域，形成初始化波浪流场，公式(1)、(2)、(3)中增加初始相位项
调整计算域波浪沿传播方向的相位移，增加z0项，用于调整波浪平衡位置的垂向位移，得到如下公式：
对流场域内网格施加相体积分数函数，以此时的波面方程为判断准则，网格垂向坐标大于波面方程的区域指定为空气，反之，则为水；对波浪水面以下区域再施加波浪速度函数以及压力场函数；波浪水平速度在波峰和波谷处数值最大且方向相反，波峰处速度与波浪传播方向相同；波浪垂向速度在迎波和背波区取得最大值且方向相反，迎波区速度与重力加速度方向相反，至此，流场域初始化完成。
进一步改进，步骤1)所述的计算域网格划分策略各加密区域采用不同形状，整个计算域内的网格均采用笛卡尔网格，飞机近场区域采用长方体，飞机前部和尾部区域采用类扇形体，其它区域按常规边界网格尺寸进行划分。
飞机近场区域网格加密等级最高，除了需要考虑附面层网格和流场物理量变化剧烈的空间网格之外，该区域内其余网格尺寸依据波浪参数指定，单个波长方向布置80个网格，单个波高方向布置20个网格；飞机前部类扇形体的扇形的角度根据飞机的运动俯仰角范围得到，此区域网格尺寸同样依据波浪参数指定；飞机尾部类扇形体的扇形角度相比于飞机前部类扇形体的扇形角度偏小，该区域网格布置调疏，单个波长方向网格不多于80个，单个波高方向网格不多于20个。
本发明有益效果在于：
1、采用多区域多等级处理计算域网格的疏密程度能够更加高效的利用空间网格进行离散求解，将网格总量保持在合理范围内，不至于过大。
2、对于研究水面飞行器波浪着水问题，该方法通过在波面方程、波浪速度场和压强场函数中增加初始相位项
和初始平衡位置项z0，实现波浪在计算域内沿传播方向的相偏移以及波浪平衡位置的垂向偏移，使得波浪在极短的求解时间内与物体接触，省去了物体等待头波的时间。
3、结合整体动网格法的速度入口造波技术，不仅使得计算效率得到了提高，而且保证了造波精度。
附图说明
图1是求解步骤。
图2是边界条件。
图3是计算域多区域多级划分策略。
图4是波浪初始相位偏移及平衡位置偏移。
图5是初始相位场。
图6是初始水平方向速度场。
图7是初始垂向速度场。
图8是初始相对压强场。
图9是某时刻水面变化情况。
图10是飞机水平方向受力。
图11是飞机垂向受力。
具体实施方式
下面结合附图对本发明作进一步说明。
本实施例中以水上飞机规则波浪着水作示范,基本求解步骤见图1。
第一步，针对复杂的自由液面情况，计算域网格采用多区域多等级划分策略，截面示意见图3。区域I为飞机近场矩形加密区；区域II为飞机迎浪扇形区；区域III为飞机尾流扇形区；区域IV为其他区。
其中，区域I网格加密等级最高，除了需要考虑附面层网格和流场物理量变化剧烈的空间网格之外，该区域内其余网格尺寸依据波浪参数指定，常规的做法是单个波长方向布置80个网格，单个波高方向布置20个网格。区域II采用类扇形划分区域，扇形的角度根据飞机可能的运动俯仰角范围得到。此区域网格尺寸同样依据波浪参数指定。区域III也采用类扇形划分，考虑到尾流区不是重点关注区域，在保证求解收敛的情况下，该扇形角度相比于区域II的角度较小。此外，该区域网格布置可相应调疏，一方面减小区域网格量，另一方面粗网格还可以起到消波的作用。区域IV则按常规边界网格尺寸进行划分。整个计算域内的网格均采用笛卡尔网格，因其高度的正交性，非常适合捕捉水气交界面变化。
第二步，为实现计算域网格随飞机作刚性运动而波浪水面的平衡位置不发生偏移，需要对各个边界条件施加边界函数。首先，确定自由液面的波面方程：
式中，h为波高，k为波数，ω为圆频率，t为时刻。并给出波浪在传播方向和垂直方向的速度函数u和w：
式中，d为水深，x和z为空间坐标。水面以下压力场函数p：
对各边界施加相体积分数函数，以波面方程为判断准则，流场网格垂向坐标大于波面方程的区域指定为空气相，反之，则为水相。此外，需额外对速度入口边界施加波浪速度函数u和w；对压力出口边界施加压力场函数p。至此，边界条件的函数施加得以完成。
第三步，根据任务需求在整个流场域内初始化波浪液面。数值模拟静水面着水问题时，初始时刻飞机最低点通常与水面保持小距离，然而，对于波浪着水问题，多是研究波浪位置和波浪参数对飞机着水特性的影响，同样需要将机体底部与波浪水面距离缩小。在上述公式中增加初始相位项
得以调整计算域波浪沿传播方向的相位移；z0项可用于调整波浪平衡位置的垂向位移，效果如图4。
对流场域内网格施加相体积分数函数，以此时的波面方程为判断准则，网格垂向坐标大于波面方程的区域指定为空气，反之，则为水，如图5。对波浪水面以下区域再施加波浪速度函数(图6和图7)以及压力场函数(图8)。图6和图7分别用矢量箭头表示初始时刻波浪速度大小与方向。波浪水平速度在波峰和波谷处数值最大且方向相反，波峰处速度与波浪传播方向相同；波浪垂向速度在迎波和背波区取得最大值且方向相反，迎波区速度与重力加速度方向相反，结果表明初始化速度场与波浪理论解一致。值得注意的是，若增加
和z0项，边界函数同样需要做相应修改。至此，初始化流场完成。图6展示的是某时刻水面变化情况。图10和图11分别为两栖飞机在波峰处着水水平方向受力和垂向受力随时间变化过程。
本发明具体应用途径很多，以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进，这些改进也应视为本发明的保护范围。
Claims (3)
1.一种水上飞机波浪水面起降水气两相流数值仿真方法，其特征在于包括以下步骤：
1)计算域网格划分策略：针对复杂的自由液面情况，将整个流场计算域做多区域多等级处理，即把计算域划分为多个网格加密区域，不同加密区域依据当地流场特征采用不同的网格加密尺寸；
2)运用整体动网格法处理物体和计算域网格非定常运动，并将速度入口造波边界条件施加在计算域的入口边界，具体施加方法如下：
2.1)确定自由液面的波面方程：
式中，h为波高，k为波数，ω为圆频率，t为时刻；
2.2)给出波浪在传播方向和垂直方向的速度函数u和w：
式中，d为水深，x和z为空间坐标；
2.3)计算水面以下压力场函数p：
2.4)对各边界施加相体积分数函数，以波面方程为判断准则，流场网格垂向坐标大于波面方程的区域指定为空气相，反之，则为水相；
2.5)对速度入口边界施加波浪速度函数u和w，对压力出口边界施加压力场函数p，至此，边界条件的函数施加得以完成；
3)流场域初始化：将关于空间坐标的初始速度场函数和压力场函数施加在流场域，形成初始化波浪流场，公式(1)、(2)、(3)中增加初始相位项
调整计算域波浪沿传播方向的相位移；增加z0项，用于调整波浪平衡位置的垂向位移，得到如下公式：
对流场域内网格施加相体积分数函数，以此时的波面方程为判断准则，网格垂向坐标大于波面方程的区域指定为空气，反之，则为水；对波浪水面以下区域再施加波浪速度函数以及压力场函数；波浪水平速度在波峰和波谷处数值最大且方向相反，波峰处速度与波浪传播方向相同；波浪垂向速度在迎波和背波区取得最大值且方向相反，迎波区速度与重力加速度方向相反，至此，流场域初始化完成。
2.根据权利要求1所述的水上飞机波浪水面起降水气两相流数值仿真方法，其特征在于：步骤1)所述的计算域网格划分策略各加密区域采用不同形状，整个计算域内的网格均采用笛卡尔网格，飞机近场区域采用长方体，飞机前部和尾部区域采用类扇形体，其它区域按常规边界网格尺寸进行划分。
3.根据权利要求2所述的水上飞机波浪水面起降水气两相流数值仿真方法，其特征在于：飞机近场区域网格加密等级最高，除了需要考虑附面层网格和流场物理量变化剧烈的空间网格之外，该区域内其余网格尺寸依据波浪参数指定，单个波长方向布置80个网格，单个波高方向布置20个网格；飞机前部类扇形体的扇形的角度根据飞机的运动俯仰角范围得到，此区域网格尺寸同样依据波浪参数指定；飞机尾部类扇形体的扇形角度相比于飞机前部类扇形体的扇形角度偏小，该区域网格布置调疏，单个波长方向网格不多于80个，单个波高方向网格不多于20个。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110371760.3A CN113515805A (zh) | 2021-04-07 | 2021-04-07 | 水上飞机波浪水面起降水气两相流数值仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110371760.3A CN113515805A (zh) | 2021-04-07 | 2021-04-07 | 水上飞机波浪水面起降水气两相流数值仿真方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113515805A true CN113515805A (zh) | 2021-10-19 |
Family
ID=78061760
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110371760.3A Pending CN113515805A (zh) | 2021-04-07 | 2021-04-07 | 水上飞机波浪水面起降水气两相流数值仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113515805A (zh) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113935114A (zh) * | 2021-10-25 | 2022-01-14 | 中航通飞华南飞机工业有限公司 | 一种用于水上飞机着水数值仿真的规则波造波方法 |
| CN114077770A (zh) * | 2021-10-25 | 2022-02-22 | 南京航空航天大学 | 一种用于飞行器水上迫降sph模拟的直接造波方法 |
| CN114117968A (zh) * | 2022-01-24 | 2022-03-01 | 南京航空航天大学 | 水气两相流全耦合求解的飞行器水上起降载荷分析方法 |
| CN114154229A (zh) * | 2021-10-21 | 2022-03-08 | 南京航空航天大学 | 一种适用于水陆两栖飞机着水降载的水翼设计方法 |
| CN115525978A (zh) * | 2022-10-31 | 2022-12-27 | 南京航空航天大学 | 用于飞行器波浪情况水动力分析的速度入口数值造波方法 |
| CN115526052A (zh) * | 2022-10-08 | 2022-12-27 | 南京航空航天大学 | 基于等分频率-能量-方向角策略不规则波浪谱解析方法 |
| CN115600316A (zh) * | 2022-10-17 | 2023-01-13 | 中国船舶科学研究中心（Cn） | 一种船底气液分层两相流波动形态数值模拟方法 |
| CN116090097A (zh) * | 2022-12-30 | 2023-05-09 | 北京机电工程研究所 | 基于等效撞水设计的近水面流固耦合有限元高效计算方法 |
| CN117150967A (zh) * | 2023-09-13 | 2023-12-01 | 南京航空航天大学 | 一种水上飞机非线性规则波波浪水面滑行运动响应数值预报方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108846225A (zh) * | 2018-06-27 | 2018-11-20 | 中国直升机设计研究所 | 一种应用于直升机水上迫降的sph波浪模拟方法 |
-
2021
- 2021-04-07 CN CN202110371760.3A patent/CN113515805A/zh active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108846225A (zh) * | 2018-06-27 | 2018-11-20 | 中国直升机设计研究所 | 一种应用于直升机水上迫降的sph波浪模拟方法 |
Non-Patent Citations (3)
| Title |
|---|
| 刘可;潘文良;李智生;: "波浪相位对航行体出水姿态影响数值仿真研究", 计算机仿真, no. 12 * |
| 孙华伟: "滑行面形状对滑行艇阻力与航态影响数值分析", 中国博士学位论文全文数据库工程科技Ⅱ辑, vol. 2015, no. 05, pages 036 - 2 * |
| 金禹彤: "楔形体入波浪水面数值模拟", 航空学报, vol. 2019, no. 10, pages 46 - 59 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114154229A (zh) * | 2021-10-21 | 2022-03-08 | 南京航空航天大学 | 一种适用于水陆两栖飞机着水降载的水翼设计方法 |
| CN113935114A (zh) * | 2021-10-25 | 2022-01-14 | 中航通飞华南飞机工业有限公司 | 一种用于水上飞机着水数值仿真的规则波造波方法 |
| CN114077770A (zh) * | 2021-10-25 | 2022-02-22 | 南京航空航天大学 | 一种用于飞行器水上迫降sph模拟的直接造波方法 |
| CN114117968A (zh) * | 2022-01-24 | 2022-03-01 | 南京航空航天大学 | 水气两相流全耦合求解的飞行器水上起降载荷分析方法 |
| CN114117968B (zh) * | 2022-01-24 | 2022-07-12 | 南京航空航天大学 | 水气两相流全耦合求解的飞行器水上起降载荷分析方法 |
| CN115526052A (zh) * | 2022-10-08 | 2022-12-27 | 南京航空航天大学 | 基于等分频率-能量-方向角策略不规则波浪谱解析方法 |
| CN115600316A (zh) * | 2022-10-17 | 2023-01-13 | 中国船舶科学研究中心（Cn） | 一种船底气液分层两相流波动形态数值模拟方法 |
| CN115525978A (zh) * | 2022-10-31 | 2022-12-27 | 南京航空航天大学 | 用于飞行器波浪情况水动力分析的速度入口数值造波方法 |
| CN116090097A (zh) * | 2022-12-30 | 2023-05-09 | 北京机电工程研究所 | 基于等效撞水设计的近水面流固耦合有限元高效计算方法 |
| CN117150967A (zh) * | 2023-09-13 | 2023-12-01 | 南京航空航天大学 | 一种水上飞机非线性规则波波浪水面滑行运动响应数值预报方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113515805A (zh) | 水上飞机波浪水面起降水气两相流数值仿真方法 | |
| Blumberg et al. | A description of a three‐dimensional coastal ocean circulation model | |
| CN110110352B (zh) | 一种基于重叠网格的船舶纵向航态与阻力预报方法 | |
| Ghassemi et al. | Interceptor hydrodynamic analysis for handling trim control problems in the high-speed crafts | |
| CN101246515A (zh) | 基于cfd软件的数字船模平面运动机构实验方法 | |
| CN111814366B (zh) | 一种内波与通气空泡作用的预报方法 | |
| CN110020459A (zh) | 基于动网格技术船随波浪在航线方向上运动模拟方法 | |
| CN115525978B (zh) | 用于飞行器波浪情况水动力分析的速度入口数值造波方法 | |
| CN112182983A (zh) | 计及海底地形及波浪影响的浮体水弹性响应分析方法 | |
| Chen et al. | Simulation research on water-entry impact force of an autonomous underwater helicopter | |
| CN111695282A (zh) | 一种基于流固耦合模拟的液舱晃荡预测和控制分析方法 | |
| Dong et al. | Numerical study on the water entry of two-dimensional airfoils by BEM | |
| CN117521556A (zh) | 一种内孤立波影响下潜艇姿态预测方法 | |
| Xin et al. | An efficient three-dimensional multiphase fluid structure interaction model on GPU for water impact of a moving body with complex geometries | |
| Zhou et al. | Development of a multiphase solver for cavitation flow near free surface | |
| Wang et al. | Propulsive performance analysis of underwater flapping multi-foil | |
| Chen et al. | Numerical wave simulation and investigation of air-wave-aircraft interactions | |
| Gan et al. | Numerical investigation of the water-entry impact performance of a bionic unmanned aerial-underwater vehicle | |
| CN116227073A (zh) | 一种研究水下发射航行体肩空泡发展溃灭的方法 | |
| CN115828426A (zh) | 水下航行体大机动工况流体动力载荷的数值计算方法 | |
| CN115293004A (zh) | 一种基于多尺度空化模型的空蚀预测方法 | |
| CN110688733B (zh) | 一种运载器式潜射导弹水面分离的数值模拟方法 | |
| Skourup et al. | Numerical modelling of wave-structure interaction by a three-dimensional nonlinear boundary element method: A step towards the numerical wave tank | |
| Zheng et al. | Numerical Analysis of Porpoising Stability Limit of a Blended-Wing–Body Aircraft During Ditching | |
| Duan et al. | Numerical Study of a Hydrodynamic Benchmark Model for Seaplanes Using OpenFOAM |
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 |