一种大长细比回转体入水空泡演化模型构建方法

1.本发明涉及水下航行体技术领域，具体涉及一种大长细比回转体入水空泡演化模型构建方法。


背景技术：

2.弹体入水过程可以视为大长细比回转体的高速入水过程，在该过程中，弹体与水的高速相对运动会使两者之间的压强减小产生低压区，当压力小于当前状态下的饱和蒸汽压时，水便会发生汽化形成空腔，包裹住回转体，将回转体与水隔开，减小其受到的阻力，从而达到提升水下航行速度的目的。为了研究大长细比回转体高速入水规律以及对入水弹道进行建模，需要对入水空泡的演化进行研究。
3.对于大长细比回转体高速入水的非定常过程，即运动体带空泡高速入水的非稳态过程，现有适用于定常条件的空泡模型的预示效果较差，无法准确模拟回转体高速入水的非定常过程。
4.因此，如何提供一种准确地进行非定常过程中空泡形态预示的大长细比回转体入水空泡演化模型构建方法是本领域技术人员亟需解决的问题。


技术实现要素：

5.有鉴于此，本发明利用数值仿真技术对空泡形态预示模型进行了修正，并验证了模型的有效性。
6.为了实现上述目的，本发明采用如下技术方案：
7.本发明公开了一种大长细比回转体入水空泡演化模型构建方法，所述大长细比回转体头部设有空化器，包括如下步骤：
8.建立大长细比回转体的定常空泡模型；
9.选取不同的回转体入水参数作为空泡初始条件进行数值仿真，得到不同空化数的空泡演化过程，通过测量仿真结果得到空泡长度以及最大空泡直径，对空泡长度、直径与空化数进行拟合，得到空泡形态与空化数的关系：
[0010][0011]
其中，rn为空化器半径，σ为空化数，p0表示空泡外的环境压力，pc表示饱和蒸汽压，通常为3540pa，为动压，ρ为水的密度，v为回转体速度，rk为空泡最大半径，lk为空泡长度；
[0012]
对于大长细比回转体高速入水的非定常过程，利用拟合得到空泡形态与空化数的关系对所述定常空泡模型进行修正，得到非定常空泡演化模型。
[0013]
优选的，所述定常空泡模型的表达式为：
[0014][0015]
其中，其中，x为空泡截面到头部空化器的距离，rc(x)为x处空泡截面的半径，rn为空化器半径，x1表示一致截面的位置，x1＝2rn，r1＝1.92rn。
[0016]
优选的，所述初始条件包括回转体入水速度以及入水角度。
[0017]
经由上述的技术方案可知，与现有技术相比，本发明的有益效果包括：
[0018]
本发明能够有效提升对于大长细比回转体高速入水非定常过程的空泡形态预示准确性，且模型拟合结果与数值仿真结果具有高吻合度，验证了本发明模型的有效性。
附图说明
[0019]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图；
[0020]
图1为本发明实施例提供的大长细比回转体入水空泡演化模型构建方法的流程图；
[0021]
图2为本发明实施例提供的大长细比航行体模型示意图；
[0022]
图3为本发明实施例提供的空泡演化模型数值仿真流程图；
[0023]
图4为本发明实施例提供的空泡演化模型数值仿真的结果与模型拟合结果对比图一；
[0024]
图5为本发明实施例提供的空泡演化模型数值仿真的结果与模型拟合结果对比图二；
[0025]
图6为本发明实施例提供的空泡演化模型数值仿真的结果与模型拟合结果对比图三。
具体实施方式
[0026]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0027]
本发明实施例公开了一种大长细比回转体入水空泡演化模型构建方法，大长细比回转体头部设有空化器，首先对回转体入水空泡形态演化规律进行说明：
[0028]
步骤一：推导独立膨胀原理
[0029]
独立膨胀原理是超空泡流动的重要理论，它能够反映出超空泡流动的实质和主要特征。回转体高速入水产生超空泡流动时，空泡的形状可以根据空泡内部压力值恒定的条
件求解得到。空泡在空化器上形成并随着空化器的运动逐渐发展，独立膨胀原理认为，将空泡分解成无数个截面，每一个截面的面积与空化器的运动过程无关，只与当前状态下空泡与环境的压力差以及空化器的运动属性相关。
[0030]
根据能量守恒原理，回转体对流体所做的功，会转化成该状态下流体的动能和势能，由此推导出空泡截面扩张方程：
[0031][0032]
式中：s表示空泡在某一时刻的截面积；δp表示在该时刻的空泡无穷远处与内部压力的差值；c
x
为回转体阻力系数，a是半经验常数，与空泡的长度和空化数有关，通常情况下在1.5～2之间；ρ为空泡边界流体的密度。对上式进行积分，可以得到空泡截面扩张速度的表达式：
[0033][0034]
其中：v(0)是0时刻回转体的运动速度，rn为回转体尾部半径。对空泡截面扩张速度再进行积分，就得到了某时刻空泡截面积s的表达式：
[0035][0036]
式中：s0为空化器的面积，s0＝πr
02
，r0为空化器半径。将同一时刻各个截面的面积沿回转体轴线绘制出来，就得到了某一时刻的完整的空泡的形态。
[0037]
步骤二：建立定常空泡模型
[0038]
依据独立膨胀原理，苏联科学家给出了近似解，提出了定常条件下的logvinovich空泡模型，其空泡模型描述的空泡形态的表达式如下：
[0039][0040]
其中：x为空泡截面到头部空化器的距离，rc(x)为x处空泡截面的半径，rn为空化器半径，x1表示一致截面的位置，x＝2rn，rk为空泡最大半径，lk为空泡长度。
[0041]
对于大长细比回转体高速入水的非定常过程，上述空泡模型的预示效果较差，所以利用数值仿真技术对空泡形态进行拟合以修正空泡模型。下面对本发明实施例提出的空泡模型修正步骤进行说明：
[0042]
包括如下步骤：
[0043]
首先建立大长细比回转体模型，回转体长细比为外形如图2所示。
[0044]
选取回转体不同的入水速度以及入水角度作为初始条件，利用cfd软件进行数值
仿真，得到不同空化数的空泡演化过程，数值仿真的流程如图3所示，首先，基于重叠网格技术利用前处理软件对回转体几何模型进行网格划分；然后，对cfd软件进行求解器的设置，包括基本控制方程、空化模型和湍流模型，在基本控制方程中引入vof多相流模型；进而求解计算记录数据，经后处理数据分析后以仿真图片的形式表征不同空化数的空泡演化过程。
[0045]
仿真图片中空泡与回转体长度以及直径的比例，利用已有的回转体模型的长度和直径计算得到空泡长度以及最大空泡直径，利用现有的软件拟合工具包对空泡长度直径与空化数之间进行拟合，得到空泡形态与空化数的关系：
[0046][0047]
其中，rn为空化器半径，σ为空化数，p0表示空泡外的环境压力，pc表示饱和蒸汽压，通常为3540pa，为动压，ρ为水的密度，v为回转体速度，rk为空泡最大半径，lk为空泡长度；
[0048]
对于大长细比回转体高速入水的非定常过程，利用拟合得到空泡形态与空化数的关系对定常空泡模型进行修正，得到非定常空泡演化模型。
[0049]
下面对非定常空泡演化模型验证仿真结果：
[0050]
为了验证拟合模型的适用性和准确性，选取回转体在空化数σ＝0.027、σ＝0.045、σ＝0.051情况，其中σ＝0.027代表了回转体在超空泡情况下的特征状态，σ＝0.045代表了回转体在空泡将要收缩至尾部时刻的特征状态，σ＝0.051表征了回转体在空泡溃灭收缩至回转体柱段的特征状态，三个特征状态就可以描述回转体的整个入水过程。将数值仿真的结果与模型拟合结果进行对比，对比结果如图4-6所示：
[0051]
图4a为σ＝0.027情况下空泡形态模型拟合情况，图4b为σ＝0.027情况下数值仿真的结果；图5a为σ＝0.045情况下空泡形态模型拟合情况，图5b为σ＝0.045情况下数值仿真的结果；图6a为σ＝0.051情况下空泡形态模型拟合情况，图6b为σ＝0.051情况下数值仿真的结果。
[0052]
图4b-6b中的黑线部分为空泡预测形态，红线部分为大长细比回转体，即航行体的模型形态。
[0053]
从图4可以看出在空化数σ＝0.027情况下拟合预示模型空泡的长度为3.97m，空泡的最大半径为0.0913m，通过与数值仿真进行对比可以看出吻合程度较大。
[0054]
从图5可以得到空化数σ＝0.045情况下拟合预示模型的空泡长度为2.65m，空泡最大半径为0.0827m，与实际情况基本一致。
[0055]
从图6可以看出在σ＝0.051情况时，在该回转体模型下空泡已经发生溃灭，空泡闭合在回转体表面的位置为距头部1.77m处，最大空泡半径为0.0801m，通过与数值仿真的图示对比后可以发现空泡模型对回转体沾湿部分的预示也较为精准，通过以上三种情况的结果分析，验证了拟合空泡预示模型的准确性。
[0056]
以上对本发明所提供的大长细比回转体入水空泡演化模型构建方法进行了详细
介绍，本实施例中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上，本说明书内容不应理解为对本发明的限制。
[0057]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本实施例中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本实施例所示的这些实施例，而是要符合与本实施例所公开的原理和新颖特点相一致的最宽的范围。