一种航行体入水空泡特性的分析方法专利检索-三相边界电池与电池专利检索查询-专利查询网

一种航行体入水 空泡特性的分析方法

阅读：168发布：2020-05-19

专利汇可以提供一种航行体入水空泡特性的分析方法专利检索，专利查询，专利分析的服务。并且本发明公开的一种航行体入水空泡特性的分析方法，属于流体机械工程领域。本发明实现方法为：对水域和空气域组成的整个流域进行网格划分及边界条件设置，并将边界数据浸入法和体积分数法运用到N‑S方程中；通过引入边界数据浸入法实现模拟多相流中浸入式物体；通过引入体积分数法实现对流体界面进行跟踪；然后运用两步投影法，对流域的速度、压力等进行求解。本发明能够实现入水航行体空泡流动数值模拟，以揭示航行体入水空泡流动规律及机理，从而为入水航行体的结构设计提供理论基础，并能够解决入水航行体实际应用工程问题，且本发明具有设计效率高、设计周期短、利于实际应用的优点。，下面是一种航行体入水空泡特性的分析方法专利的具体信息内容。

权利要求

1.一种航行体入水空泡特性的分析方法，其特征在于：包括如下步骤，
步骤一：对流体计算域进行网格划分；
步骤二：对整个流域进行边界条件设置；
步骤三：建立流体域的控制方程及固体域的运动方程；
步骤四：为模拟多相流中浸入式物体，将边界数据浸入法用于流体控制方程和入水航行体运动方程，以便后续步骤六求解N-S方程；
步骤五：为了对流体界面进行跟踪，采用VOF方法；将VOF方法用于流体控制方程，以便后续步骤六求解N-S方程；
步骤六：用两步投影法对N-S方程进行求解，即实现对流域的速度、压力进行求解；
步骤四具体实现方法为，
在浸入式无滑移固体-流体相互作用的问题中，为了模拟多相流中浸入式物体，引入边界数据浸入法，将流域分成固体区域σb，流体区域σf,固体和流体界面宽度为2ε，定义σb子域内任意一点到流固界面中心的距离为d，方向朝向子域内部为负，引入系数函数：
使得：
2.如权利要求1所述的一种航行体入水空泡特性的分析方法，其特征在于：步骤一具体实现方法为，确定入水航行体几何参数，对整个流域进行网格划分，所述整个流域包括水域和空气域；采用的是笛卡尔网格，并采用均布化，归一化原理，航行体选用为球体，以其直径为一个单位，将网格设置为六面体网格，棱长为0.01-0.05之间便于计算。
3.如权利要求1所述的一种航行体入水空泡特性的分析方法，其特征在于：步骤二具体实现方法为，对整个流域进行边界条件设置，空气域一侧为压力入口，水域一侧边界设置为压力出口，入水航行体和计算域外围侧面设置为无滑移壁面。
4.如权利要求1所述的一种航行体入水空泡特性的分析方法，其特征在于：步骤三具体实现方法为，
由于模拟对象是入水航行体，所受环境压力较小，所以为了简化计算，将计算中涉及到的气体和液体都视为不可压缩的流体；在不考虑流体的压缩性的条件下，把流体的密度视为常量；对于不可压缩流体，动量守恒方程是经典的纳维-斯托克斯方程：
整理为：
其中ρ流体的密度，μ是流体的动力学粘性系数，和p分别是流场的速度和压力，σ是表面张力系数,κ为局部曲率，为自由液面的法向量，δS为Dirac函数，g是重力加速度；
对于不可压缩流体，因为流体密度不变方程中省略密度相，所以连续性方程比较简单：
而固体的运动方程，由下式给出：
其中U为固壁速度，
其中为球心速度，ω球体角速度，r为球体半径
方程(1)、(2)即为建立的流体控制方程，方程(3)即为建立的固体运动方程。
5.如权利要求1所述的一种航行体入水空泡特性的分析方法，其特征在于：步骤五具体实现方法为，
为了对流体界面进行跟踪，采用VOF方法；计算域由三相组成：液态水，气态空气和实心球；通过笛卡尔网格对空间进行离散化，以网格中流体所占体积分数α来构造来跟踪自由面；网格为流体充满,α设置为1；若网格不包含流体,α为0；当α介于0和1之间，该网格定义为自由面网格；计算公式如下：
自由面的法向量计算公式为：
6.如权利要求1所述的一种航行体入水空泡特性的分析方法，其特征在于：步骤五具体实现方法为，
使用两步投影法来解连续性方程和动量方程，在第一步中，通过离散对流项和体力项来获得
对于ρ和μ，作如下处理：
ρ＝αρwater+(1-α)ρair (8)
μ＝αμwater+(1-α)μair (9)
其中α为水的体积分数，ρwater和ρair分别为水和空气的密度，μwater和μair分别为水和空气的动力学粘性系数；第二步通过以下方程来获得新时刻的速度：
将散度算子运用于式(7)，能得到压力泊松方程：
通过方程(7)和方程(10)实现对速度的求解，通过方程(11)实现对压力的求解，即实现对流域的速度、压力进行求解。

说明书全文

一种航行体入水 空泡特性的分析方法

技术领域

[0001] 本发明属于流体机械工程领域，涉及一种航行体入水空泡特性的分析方法，是基于自编程的一种航行体入水空泡特性的分析方法。

背景技术

[0002] 航行体入水问题涉及固体，液体和气体三种相态互相作用，入水过程中涵盖跨介质、多相流、强湍流、可压缩等众多物理难题。航行体入水问题是常见的自然现象，因此在众多的工程领域及科学研究中有着重大应用，例如：宇宙飞船的水上回收，救生艇的抛落，水上动物的行走，跳水运动员压水花等。早期对入水问题的研究主要集中于航行体几何参数、入水角度及入水速度等，而考虑航行体自身旋转时，问题就更加复杂并且少有研究，有很强的科学及工程研究价值。

[0003] 目前，针对这些难题还没有比较准确的分析手段和数值计算方法，有必要提供一种能够准确航行体入水空泡特性的数值计算方法，以揭示航行体入水空泡流动规律及机理，从而为入水航行体的结构设计提供理论基础。

发明内容

[0004] 本发明公开的一种航行体入水空泡特性的分析方法要解决的技术问题是：实现入水航行体空泡流动数值模拟，以揭示航行体入水空泡流动规律及机理，从而为入水航行体的结构设计提供理论基础，并能够解决入水航行体实际应用工程问题，且本发明具有设计效率高、设计周期短、利于实际应用的优点。

[0005] 所述解决入水航行体实际应用工程问题包括宇宙飞船的水上回收、救生艇的抛落、水上动物的行走、跳水运动员压水花。

[0006] 本发明的目的是通过下述技术方案实现的。

[0007] 本发明公开的一种航行体入水空泡特性的分析方法，对水域和空气域组成的整个流域进行网格划分及边界条件设置，并将边界数据浸入法和体积分数法运用到N-S方程中。通过引入边界数据浸入法实现模拟多相流中浸入式物体；通过引入体积分数法实现对流体界面进行跟踪。然后运用两步投影法，对流域的速度、压力等进行求解。在求解过程中，第一步先根据离散对流项和体力项来获得速度的中间变量，第二步再通过压力项求解流域的速度，并且获得压力泊松方程，求解压力泊松方程得到流域的压力。

[0008] 本发明公开的一种航行体入水空泡特性的分析方法，包括如下步骤：

[0009] 步骤一：对流体计算域进行网格划分。

[0010] 确定入水航行体几何参数，对整个流域进行网格划分，所述整个流域包括水域和空气域。采用的是笛卡尔网格，并采用均布化，归一化原理，航行体选用为球体，以其直径为一个单位，将网格设置为六面体网格，棱长为0.01-0.05之间便于计算。

[0011] 步骤二：对整个流域进行边界条件设置。

[0012] 对整个流域进行边界条件设置，空气域一侧为压力入口，水域一侧边界设置为压力出口，入水航行体和计算域外围侧面设置为无滑移壁面。

[0013] 步骤三：建立流体域的控制方程及固体域的运动方程。

[0014] 由于模拟对象是入水航行体，所受环境压力较小，所以为了简化计算，将计算中涉及到的气体和液体都视为不可压缩的流体。在不考虑流体的压缩性的条件下，把流体的密度视为常量。对于不可压缩流体，动量守恒方程是经典的纳维-斯托克斯方程：

[0015]

[0016] 整理为：

[0017]

[0018] 其中ρ流体的密度，μ是流体的动力学粘性系数，和p分别是流场的速度和压力，σ是表面张力系数,κ为局部曲率，为自由液面的法向量，δS为Dirac函数，g是重力加速度。

[0019] 对于不可压缩流体，因为流体密度不变方程中省略密度相，所以连续性方程比较简单：

[0020]

[0021] 而固体的运动方程，由下式给出：

[0022]

[0023] 其中U为固壁速度，

[0024] U＝u1+ω×r (4)

[0025] 其中u1为球心速度，ω球体角速度，r为球体半径

[0026] 方程(1)、(2)即为建立的流体控制方程，方程(3)即为建立的固体运动方程。

[0027] 步骤四：为模拟多相流中浸入式物体，将边界数据浸入法用于流体控制方程和入水航行体运动方程，以便后续步骤六求解N-S方程。

[0028] 在浸入式无滑移固体-流体相互作用的问题中，为了模拟多相流中浸入式物体，引入边界数据浸入法，将流域分成固体区域σb，流体区域σf,固体和流体界面宽度为2ε，定义σb子域内任意一点到流固界面中心的距离为d，方向朝向子域内部为负，引入系数函数：

[0029]

[0030] 使得：

[0031]

[0032] 步骤五：为了对流体界面进行跟踪，采用VOF方法。将VOF方法用于流体控制方程，以便后续步骤六求解N-S方程。

[0033] 为了对流体界面进行跟踪，采用VOF方法。计算域由三相组成：液态水，气态空气和实心球。通过笛卡尔网格对空间进行离散化，以网格中流体所占体积分数α来构造来跟踪自由面。网格为流体充满,α设置为1；若网格不包含流体,α为0；当α介于0和1之间,该网格定义为自由面网格。计算公式如下：

[0034]

[0035] 自由面的法向量计算公式为：

[0036]

[0037] 步骤六：用两步投影法对N-S方程进行求解，即实现对流域的速度、压力等进行求解。

[0038] 使用两步投影法来解连续性方程和动量方程，在第一步中，通过离散对流项和体力项来获得

[0039]

[0040] 对于ρ和μ，作如下处理：

[0041] ρ＝αρwater+(1-α)ρair (8)

[0042] μ＝αμwater+(1-α)μair (9)

[0043] 其中α为水的体积分数，ρwater和ρair分别为水和空气的密度，μwater和μair分别为水和空气的动力学粘性系数。第二步通过以下方程来获得新时刻的速度：

[0044]

[0045] 将散度算子运用于式(7)，能得到压力泊松方程：

[0046]

[0047] 通过方程(7)和方程(10)实现对速度的求解，通过方程(11)实现对压力的求解，即实现对流域的速度、压力等进行求解。

[0048] 步骤七：将步骤一至步骤六所述的方法应用于水航行体空泡流动领域，实现入水航行体空泡流动数值模拟，以揭示航行体入水空泡流动规律及机理，从而为入水航行体的结构设计提供理论基础，并能够解决入水航行体实际应用工程问题。

[0049] 所述解决入水航行体实际应用工程问题包括宇宙飞船的水上回收、救生艇的抛落、水上动物的行走、跳水运动员压水花。

[0050] 有益效果：

[0051] 1、本发明公开的一种航行体入水空泡特性的分析方法，为了模拟多相流中浸入式物体，引入边界数据浸入法；为了对流体界面进行跟踪，采用体积分数方法，将边界数据浸入法和体积分数法运用到N-S方程中，并用两步投影法对N-S方程进行求解，即实现对流域的速度、压力等进行求解。

[0052] 2、本发明公开的一种航行体入水空泡特性的分析方法，将结果与经典实验数据进行对比，并验证了数值计算方法的准确性，可行性。

[0053] 3、本发明公开的一种航行体入水空泡特性的分析方法，对航行体入水时的空泡流动过程进行数值仿真分析，以揭示空泡流动规律及机理，从而为入水航行体的结构设计提供理论基础，并能够解决入水航行体实际应用工程问题。附图说明

[0054] 图1本发明的一种航行体入水空泡特性的分析方法流程图；

[0055] 图2本发明实施例中网格划分示意图；

[0056] 图3本发明实施例中边界条件设置示意图；

[0057] 图4本发明实施例中数值方法准确性验证云图；

[0058] 图5本发明实施例中数值方法准确性验证位移图；

[0059] 图6本发明实施例中不同旋转入水速度下数值模拟结果体积分数云图；

[0060] 图7本发明实施例中不同入水旋转角速度下数值模拟结果体积分数云图。

具体实施方式

[0061] 为了更好的说明本发明涉及的一种航行体入水空泡特性的分析方法，利用本发明的方法结合附图和实施例进行了测试计算，使得技术方案和有益效果更加清楚。

[0062] 实施例1：

[0063] 本实施例以国外公开的球体入水为研究对象，直d＝0.0572m的球体(并将d作为特征长度)，初始入水速度为v0＝2.5m/s(将v0作为特征速度)，初始旋转角速度为210rad/s,方向为逆时针，所选模型为台球。

[0064] 如图1所示，本实施例公开的一种航行体入水空泡特性的分析方法，包括如下步骤：

[0065] 步骤一：对流体计算域进行网格划分。

[0066] 确定入水航行体几何参数，对整个流域进行网格划分，所述整个流域包括水域和空气域。采用的是笛卡尔网格，并采用均布化，归一化原理，航行体选用为球体，以其直径为一个单位，将网格设置为六面体网格，棱长为0.03之间便于计算。

[0067] 步骤二：对整个流域进行边界条件设置。

[0068] 对整个流域进行边界条件设置，空气域一侧为压力入口，水域一侧边界设置为压力出口，入水航行体和计算域外围侧面设置为无滑移壁面。

[0069] 步骤三：建立流体域的控制方程及固体域的运动方程

[0070] 由于模拟对象是入水航行体，所受环境压力较小，所以为了简化计算，将计算中涉及到的气体和液体都视为不可压缩的流体。在不考虑流体的压缩性的条件下，把流体的密度视为常量。对于不可压缩流体，动量守恒方程是经典的纳维-斯托克斯方程：

[0071]

[0072] 整理为：

[0073]

[0074] 其中ρ流体的密度，μ是流体的动力学粘性系数，和p分别是流场的速度和压力，σ是表面张力系数,κ为局部曲率，为自由液面的法向量，δS为Dirac函数，g是重力加速度。

[0075] 对于不可压缩流体，因为流体密度不变方程中省略密度相，所以连续性方程比较简单：

[0076]

[0077] 而固体的运动方程，由下式给出：

[0078]

[0079] 其中U为固壁速度，

[0080]

[0081] 其中为球心速度，ω球体角速度，r为球体半径

[0082] 方程(1)、(2)即为建立的流体控制方程，方程(3)即为建立的固体运动方程。

[0083] 步骤四：为模拟多相流中浸入式物体，将边界数据浸入法用于流体控制方程和入水航行体运动方程，以便后续步骤六求解N-S方程。

[0084] 在浸入式无滑移固体-流体相互作用的问题中，为了模拟多相流中浸入式物体，引入边界数据浸入法，将流域分成固体区域σb，流体区域σf,固体和流体界面宽度为2ε，定义σb子域内任意一点到流固界面中心的距离为d，方向朝向子域内部为负，引入系数函数：

[0085]

[0086] 使得：

[0087]

[0088] 步骤五：为了对流体界面进行跟踪，采用VOF方法。将VOF方法用于流体控制方程，以便后续步骤六求解N-S方程。

[0089] 为了对流体界面进行跟踪，采用VOF方法。计算域由三相组成：液态水，气态空气和实心球。通过笛卡尔网格对空间进行离散化，以网格中流体所占体积分数α来构造来跟踪自由面。网格为流体充满,α设置为1；若网格不包含流体,α为0；当α介于0和1之间,该网格定义为自由面网格。计算公式如下：

[0090]

[0091] 自由面的法向量计算公式为：

[0092]

[0093] 步骤六：用两步投影法对N-S方程进行求解，即实现对流域的速度、压力等进行求解。

[0094] 使用两步投影法来解连续性方程和动量方程，在第一步中，通过离散对流项和体力项来获得

[0095]

[0096] 对于ρ和μ，作如下处理：

[0097] ρ＝αρwater+(1-α)ρair (8)

[0098] μ＝αμwater+(1-α)μair (9)

[0099] 其中α为水的体积分数，ρwater和ρair分别为水和空气的密度，μwater和μair分别为水和空气的动力学粘性系数。第二步通过以下方程来获得新时刻的速度：

[0100]

[0101] 将散度算子运用于式(7)，能得到压力泊松方程：

[0102]

[0103] 通过方程(7)和方程(10)实现对速度的求解，通过方程(11)实现对压力的求解，即实现对流域的速度、压力等进行求解。

[0104] 步骤七：将步骤一至步骤六所述的方法应用于水航行体空泡流动领域，实现入水航行体空泡流动数值模拟，以揭示航行体入水空泡流动规律及机理，从而为入水航行体的结构设计提供理论基础，并能够解决入水航行体实际应用工程问题。

[0105] 以上所述的具体描述，对发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

标题	发布/更新时间	阅读量
一种电网故障下并网变流器韧性评估方法	2020-05-12	934
基于最优分区的主动配电网分布式区间状态估计方法	2020-05-14	691
获取负荷电流不平衡度的方法、估算装置与系统	2020-05-15	403
一种基于区间状态估计的配电网数据安全检测方法	2020-05-16	357
一种褐煤微波干燥多物理场仿真模拟方法	2020-05-17	384
非接触式变电站氧化锌避雷器母线电压相位在线监测方法	2020-05-11	648
一种基于ADAL输配一体化并行状态估计方法	2020-05-14	557
永磁同步电机无电解电容驱动器电压边界优化过调制方法	2020-05-16	586
飞行器飞行信息系统和方法	2020-05-11	209
基于多边形随机骨料的混凝土三相细观模型构建及内氯离子侵蚀数值模拟方法	2020-05-18	874

一种航行体入水空泡特性的分析方法

一种航行体入水空泡特性的分析方法

技术领域

背景技术

发明内容

具体实施方式

该功能需要专业版企业版VIP权限，您可以：