CN115481583B - 一种针对非均匀来流推进器性能试验的轴向伴流模拟方法 - Google Patents

一种针对非均匀来流推进器性能试验的轴向伴流模拟方法 Download PDF

Info

Publication number
CN115481583B
CN115481583B CN202211190018.3A CN202211190018A CN115481583B CN 115481583 B CN115481583 B CN 115481583B CN 202211190018 A CN202211190018 A CN 202211190018A CN 115481583 B CN115481583 B CN 115481583B
Authority
CN
China
Prior art keywords
node
wake
coordinate
model
nodes
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.)
Active
Application number
CN202211190018.3A
Other languages
English (en)
Other versions
CN115481583A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202211190018.3A priority Critical patent/CN115481583B/zh
Publication of CN115481583A publication Critical patent/CN115481583A/zh
Application granted granted Critical
Publication of CN115481583B publication Critical patent/CN115481583B/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/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
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • 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
    • G06F2113/00—Details relating to the application field
    • G06F2113/10—Additive manufacturing, e.g. 3D printing
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computing Systems (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Computer Graphics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种针对非均匀来流推进器性能试验的轴向伴流模拟方法，包括：(1)采用四边型蜂窝三维结构作为伴流模型的结构，建立伴流模型尺寸与伴流分数大小对应的函数关系；(2)基于船艉目标伴流等值图，将伴流区域划分成四边形网格；(3)提取四边形网格的节点坐标信息，利用函数关系将节点向临近推进一个、两个或者四个节点，生成三维建模软件所需的节点坐标格式文件；(4)进行三维建模；(5)对建好的三维模型进行数值计算，并调整网格大小直至得到精度满足要求的伴流模型；(6)采用3D打印技术打印伴流模型。本发明可以快速准确的得到目标伴流场，可适用于各种尺度大小的水池、空泡水筒或水槽开展非均匀来流推进器性能试验。

Description

一种针对非均匀来流推进器性能试验的轴向伴流模拟方法
技术领域
本发明属于船用推进器领域，尤其是涉及一种针对非均匀来流推进器性能试验的轴向伴流模拟方法。
背景技术
物理试验是准确评估船用推进器的空泡、激振力和噪声性能的有效手段。在专业水动力性能实验室中进行非均匀来流下推进器试验时，伴流场的准确模拟则是推进器空泡振动噪声性能试验研究的一个重要环节。
船艉伴流场通常可以采用缩比整船模带全附体模型、缩短的假体伴流模型以及金属网格模拟方法。其中，船模带全附体伴流就是将实尺度船舶按照一定比例整体缩放，实船与船模完全几何相似。缩短的假体伴流模型则是在最大程度保持艉部几何相似的情况下，将整船模缩短的伴流模型。金属网格伴流则是根据船艉伴流等值图，用多层疏密特殊配制的金属网格模拟出目标伴流场。
整船伴流虽然几何相似，但是其制作加工成本高，只适用于大尺度水池、空泡水筒或水槽实验室，不能用于中小型空泡水筒或水槽；而且由于尺度效应，其船艉伴流分布并不与实船相同。
缩短的假体伴流模型则是不完全几何相似的，适用于大中小水池、空泡水筒或水槽，但是需要多次调整修改其外部形状近似模拟目标轴向伴流，但这种修改过程需要丰富的经验，修改周期长，而且需要反复加工，模拟的人工成本和加工成本都较高。
金属网格伴流，与实船不完全相似，仅模拟目标轴向伴流，适用于大中小水池、空泡水筒或水槽，但求实验室具有丰富经验的伴流模拟技术人员和手工精巧的钳工，多次修改调整测量，模拟周期较长，且背景噪声较高，不能用于噪声测量试验。
如公开号为112284681A的中国专利文献公开了一种空泡水筒桨模试验用伴流场模拟方法，利用网格+附体模型，可较准确且较方便的模拟带船尾节能附体的单桨船舶或带轴支架的双桨船舶的三向伴流，在此伴流场中开展桨模试验更接近实船工况，试验结果更加准确可靠。网格采用多层金属丝编织而成，固定在空泡水筒测试段水流进口附近；附体模型采用3D打印等方法得到，周期短、成本低，且方便安装在空泡水筒试验轴上的桨模紧前方位置。该方法拓展了中型空泡水筒的试验能力，可更加细致的观察桨模前方的节能附体或轴支架等对螺旋桨性能产生的影响。
但是，上述方法，整船伴流虽然几何相似，但是其制作加工成本及试验成本较高；缩短的假体伴流模型则是不完全几何相似的，需要反复加工，模拟的人工成本和加工成本都较高；金属网格伴流，要求实验室具有丰富经验的伴流模拟技术人员和手工精巧的钳工，模拟周期较长，且不能用于噪声测量试验。
发明内容
本发明提供了一种针对非均匀来流推进器性能试验的轴向伴流模拟方法，可采用3D打印伴流模型快速准确获得低噪声目标伴流场，适用于各种尺度大小的水池、空泡水筒或水槽开展非均匀来流推进器性能试验。
一种针对非均匀来流推进器性能试验的轴向伴流模拟方法，包括以下步骤：
(1)采用四边型蜂窝三维结构作为伴流模型的结构，对伴流模型所产生的流场信息进行数值计算分析，将数值计算结果与试验测量值进行对比，矫正计算方法后计算不同尺寸的伴流模型所诱导的流场，建立伴流模型的尺寸与伴流分数大小对应的函数关系；
(2)基于船艉目标伴流等值图，提取伴流区域边界和等值线几何坐标信息，以此作为控制伴流分数大小的边界条件，并基于数值网格生成方法，将伴流区域划分成同尺寸四边形网格；
(3)提取四边形网格的节点坐标信息，利用步骤(1)得到的函数关系将节点向临近推进一个、两个或者四个节点，生成三维建模软件所需的节点坐标格式文件；
(4)在三维建模软件中依次导入节点坐标信息，将坐标点按顺序相连成线，再进行拉伸成面，最后进行面加厚处理；
(5)对三维建模软件建好的三维模型进行数值计算，将计算得到的伴流分布与目标伴流进行对比，如有误差，则细微调整三维建模软件所需的节点坐标，重复步骤(4)～(5)，直至得到精度满足要求的伴流模型；
(6)采用3D打印技术打印伴流模型，以获得低噪声目标伴流场。
进一步地，步骤(1)中，建立伴流模型的尺寸与伴流分数大小对应的函数关系如下：
wx＝2.11-1.02*λ+0.211*λ2-0.016*λ3
其中，λ为四边型蜂窝的网格孔厚比，λ＝d/b，d为网格孔径，b为网格壁厚，b≥0.3mm，1.5≤λ≤6.5。
函数关系的限制条件为：b+d的值为0.0075D～0.02D，伴流模型与推进器转子盘面的距离为0.02D～2D，其中，D为推进器转子直径。
建立伴流模型的尺寸与伴流分数大小对应的函数关系具体过程为：
基于Realizable k-ε湍流模型，应用计算流体力学软件对λ＝2.0的伴流模型进行流场计算，与实验结果对比后修正数值计算方法；最后采用修正后的数值计算方法计算不同λ下伴流分数的值，经数据拟合得到λ与wx函数关系。
步骤(2)的具体过程为：
将船艉目标伴流等值图放入作图软件中，缩放至实验用大小，描绘出等值边界线，接着提取等值线几何坐标信息，以此作为控制孔厚比的初始区域边界；
将初始区域边界向伴流区域内平移0.5mm，构成新的区域边界；然后，采用四边型网格对新区域边界内的伴流区域进行划分，初始网格划分尺寸为b+d＝0.01D；对网格节点进行编号并提取几何坐标信息保存。
步骤(3)中，将节点向临近推进一个、两个或者四个节点的依据为：
当节点1位于新区域边界线的交点时，其坐标点记为(x1,y1)，将该节点向临近推进一个节点，生成新的一个节点1-1，坐标点记为(x1-1,y1-1)；
当节点2位于新区域边界线上的网格节点处，其坐标点记为(x2,y2)，将该节点向临近推进两个节点，生成新的节点2-1和节点2-2，对应的坐标点分别记为(x2-1,y2-1)、(x2-2,y2-2)；
当节点3位于伴流区域内网格节点处，其坐标点记为(x3,y3)，将该节点向临近推进四个节点，生成新的节点3-1、节点3-2、节点3-3和节点3-4，对应的坐标点分别记为(x3-1,y3-1)、(x3-2,y3-2)、(x3-3,y3-3)、(x3-4,y3-4)；
当节点4位于伴流分数为wx,b的区域内，但等值边界曲线穿过其所在的网格，其坐标点记为(x4,y4)，该节点向临近推进四个节点，生成新的节点4-1、节点4-2、节点4-3和节点4-4，对应的坐标点分别记为(x4-1,y4-1)、(x4-2,y4-2)、(x4-3,y4-3)、(x4-4,y4-4)。
节点1-1与节点1的坐标关系为：
x1-1＝x1+dx y1-1＝y1+dy
其中，dx＝0.5ba，dy＝0.5ba，
Figure BDA0003868943770000041
da＝λaba；
该区域伴流分数为wx,a，对应的孔厚比为λa，则网格孔径为da，网格壁厚为ba。
节点2-1、节点2-2与节点2的坐标关系分别为：
x2-1＝x2+dx y2-1＝y2-dy
x2-2＝x2+dx y2-2＝y2+dy
式中，dx和dy均为0.5ba。
节点3-1、节点3-2、节点3-3、节点3-4与节点3的坐标关系分别为：
x3-1＝x3-dx y3-1＝y3-dy
x3-2＝x3-dx y3-2＝y3+dy
x3-3＝x3+dx y3-3＝y3+dy
x3-4＝x3+dx y3-4＝y3-dy
式中，dx和dy均为0.5ba。
节点4-1、4-2与节点4的坐标关系同节点3-1、3-2与节点3的坐标关系；对于节点4-3和节点4-4，首先判定节点4与等值边界曲线的最短直线距离L与初始网格划分尺寸的关系：
如果L＞0.005D，则节点4-3、节点4-4与节点4的坐标关系同节点3-3、节点3-4与节点3的坐标关系，dx＝dy＝0.5bb；bb为伴流分数为wx,b区域的网格壁厚；
如果L＜0.005D，则节点4-3、节点4-4与节点4的坐标关系同节点3-3、节点3-4与节点3的坐标关系，但dx＝dy＝0.5bc；bc为伴流分数为wx,c区域的网格壁厚；
如果L＝0.005D，则节点4-3、节点4-4与节点4的坐标关系同节点3-3、节点3-4与节点3的坐标关系，但dx＝dy＝(0.5bb+0.5bc)/2。
与现有技术相比，本发明具有以下有益效果：
本发明基于数值计算分析、数值网格划分、网格节点坐标提取转换、智能化三维建模和3D打印技术，可快速准确的得到目标伴流场，背景噪声较低，适用于各种尺度大小的水池、空泡水筒或水槽开展非均匀来流推进器性能试验。
附图说明
图1为本发明一种针对非均匀来流推进器性能试验的轴向伴流模拟方法流程图；
图2为本发明实施例中四边型蜂窝三维结构的平面示意图；
图3为本发明实施例中伴流区域的初始区域边界示意图；
图4为本发明实施例中伴流区域的新的区域边界示意图；
图5为图4中节点1、节点2和节点3生成新节点的示意图；
图6为图4中节点4生成新节点的示意图；
图7本发明实施例制作的伴流模型图。
具体实施方式
下面结合附图和实施例对本发明做进一步详细描述，需要指出的是，以下所述实施例旨在便于对本发明的理解，而对其不起任何限定作用。
如图1所示，一种一种针对非均匀来流推进器性能试验的轴向伴流模拟方法，主要包括如下步骤：
1)对四边型蜂窝三维结构模型所产生的流场信息进行了数值计算分析，将数值结果与试验测量值进行对比，矫正计算方法后计算了不同尺寸的结构模型所诱导的流场，建立了结构模型尺寸与伴流分数大小对应函数关系。
如图2所示，为四边型蜂窝三维结构的平面示意图，z轴正方向满足右手定则。网格孔径为d(单位：mm)，网格壁厚为b(单位：mm)，推进器转子直径为D(单位：mm)，孔厚比为λ：
λ＝d/b(b≥bmin，一般为0.3mm，代表模型加工最小要求厚度)
来流速度为u，伴流场中的轴向速度为ux，伴流分数为wx：
wx＝1-ux/u
基于Realizable k-ε湍流模型，应用计算流体力学软件对λ＝2.0三维结构模型(其中z轴方向模型长度为0.08D、0.12mm、0.15mm和0.2mm)进行流场计算，与实验结果对比修正了数值计算方法，建议选用0.12mm长度即可，即能满足结构强度，也能满足流场设计要求。最后采用修正后的数值方法计算了不同λ下伴流分数的值。经数据拟合得到λ与wx函数关系：
wx＝2.11-1.02*λ+0.211*λ2-0.016*λ3(1.5≤λ≤6.5)
所示函数关系的限制条件为：b+d为0.0075D～0.02D，伴流模型距离推进器转子盘面0.02D～2D。
2)基于船艉目标伴流等值图，提取伴流区域边界和等值线几何坐标信息，以此作为控制伴流分数大小的边界条件，基于数值网格生成方法，将伴流区域划分成同尺寸四边形网格。
螺旋桨工作于船后的艉流场中，通常由拖曳水池船模艉流场测量实验获得该流场分布，既为船艉伴流等值图。将等值图放入CAD或其它作图软件中，缩放至实验用大小(一般为1.5D的圆形区域)，描绘出等值边界线，最后提取等值线几何坐标信息，以此作为控制孔厚比的初始区域边界，如图3所示。
首先，为增强伴流发生器的几何结构稳定性，区域边界向伴流区域内平移0.5mm，构成新的区域边界条件，如图4所示。其次，采用四边型网格对伴流区域进行划分，为满足曲线边界条件，通常网格尺寸越小对于边界的拟合就越好，但本发明建议初始网格划分尺寸为b+d＝0.01D；对网格节点进行编号并提取几何坐标信息保存。
3)提取网格节点坐标信息，并将节点向临近推进一个、两个或者四个节点，生成三维建模软件所需的坐标格式文件。
基于函数关系，还需对网格划分后的网格节点坐标进行处理以生成四边型蜂窝三维结构，鉴于此，对4个具有代表性的节点1，节点2，节点3和节点4进行详细阐述，参考图5和图6。
节点1：节点1位于新区域边界线的交点，其坐标点记为(x1,y1)。这时需生成新的节点1-1以满足结构建模要求，节点1-1坐标点记为(x1-1,y1-1)。该区域伴流分数为wx,a，对应的孔厚比为λa。则网格孔径为da，网格壁厚为ba为：
Figure BDA0003868943770000081
则节点1与节点1-1的坐标关系为：
x1-1＝x1+dxy1-1＝y1+dy(dx＝0.5ba，dy＝0.5ba)
节点2：节点2位于新区域边界线上的网格节点处，其坐标点记为(x2,y2)。这时需生成新的两个节点2-1和节点2-2以满足结构建模要求，节点2-1坐标点记为(x2-1,y2-1)，节点2-2坐标点记为(x2-2,y2-2)。节点2-1、节点2-2与节点2的坐标关系为：
x2-1＝x2+dx y2-1＝y2-dy
x2-2＝x2+dx y2-2＝y2+dy
式中，dx和dy均为0.5ba。
节点3：该点位于伴流区域内网格节点处，其坐标点记为(x3,y3)。这时需生成新的四个节点3-1、节点3-2、节点3-3和节点3-4以满足结构建模要求，对应的坐标点分别记为(x3-1,y3-1)、(x3-2,y3-2)、(x3-3,y3-3)、(x3-4,y3-4)，与节点3的坐标关系分别为：
x3-1＝x3-dx y3-1＝y3-dy
x3-2＝x3-dx y3-2＝y3+dy
x3-3＝x3+dx y3-3＝y3+dy
x3-4＝x3+dx y3-4＝y3-dy
式中，dx和dy均为0.5ba。
节点4：该点位于伴流分数为wx,b的区域内，但等值边界曲线穿过其所在的网格，其坐标点记为(x4,y4)。这时需生成新的四个节点4-1、节点4-2、节点4-3和节点4-4以满足结构建模要求，对应的坐标点分别记为(x4-1,y4-1)、(x4-2,y4-2)、(x4-3,y4-3)、(x4-4,y4-4)。其中节点4-1、4-2与节点4的坐标关系同3-1、3-2与节点3的坐标关系；对于节点4-3和4-4而言，首先判定节点4与等值边界曲线的最短直线距离L与基础网格尺寸(0.01D)的关系：
如果L＞0.005D，则节点4-3和4-4与节点4的坐标关系同3-3、3-4与节点3的坐标关系，dx＝dy＝0.5bb；
如果L＜0.005D，则节点4-3和4-4与节点4的坐标关系同3-3、3-4与节点3的坐标关系，但dx＝dy＝0.5bc；
如果L＝0.005D，则节点4-3和4-4与节点4的坐标关系同3-3、3-4与节点3的坐标关系，但dx＝dy＝(0.5bb+0.5bc)/2。
4)在三维建模软件中依次导入节点坐标信息，将坐标点按顺序相连成线，再进行拉伸成面，最后进行面加厚处理。
5)基于Realizable k-ε湍流模型，应用计算流体力学软件对建好的三维模型进行数值计算，将计算得到的伴流分布与目标伴流进行对比，如有误差，则细微调整三维建模软件所需的节点坐标，重复步骤(4)～(5)直至精度满足要求。
6)采用3D打印伴流模型，快速准确获得低噪声目标伴流场，适用于各种尺度大小的水池、空泡水筒或水槽开展非均匀来流推进器性能试验。
如图7所示，为采用本发明实施例制作的伴流模型图。
以上所述的实施例对本发明的技术方案和有益效果进行了详细说明，应理解的是以上所述仅为本发明的具体实施例，并不用于限制本发明，凡在本发明的原则范围内所做的任何修改、补充和等同替换，均应包含在本发明的保护范围之内。

Claims (8)

1.一种针对非均匀来流推进器性能试验的轴向伴流模拟方法，其特征在于，包括以下步骤：
(1)采用四边型蜂窝三维结构作为伴流模型的结构，对伴流模型所产生的流场信息进行数值计算分析，将数值计算结果与试验测量值进行对比，矫正计算方法后计算不同尺寸的伴流模型所诱导的流场，建立伴流模型的尺寸与伴流分数大小对应的函数关系；建立伴流模型的尺寸与伴流分数大小对应的函数关系如下：
wx＝2.11-1.02*λ+0.211*λ2-0.016*λ3
其中，wx表示轴向伴流分数，λ为四边型蜂窝的网格孔厚比，λ＝d/b，d为网格孔径，b为网格壁厚，b≥0.3mm，1.5≤λ≤6.5；
建立伴流模型的尺寸与伴流分数大小对应的函数关系具体过程为：
基于Realizable k-ε湍流模型，应用计算流体力学软件对λ＝2.0的伴流模型进行流场计算，与实验结果对比后修正数值计算方法；最后采用修正后的数值计算方法计算不同λ下伴流分数的值，经数据拟合得到λ与wx函数关系；
(2)基于船艉目标伴流等值图，提取伴流区域边界和等值线几何坐标信息，以此作为控制伴流分数大小的边界条件，并基于数值网格生成方法，将伴流区域划分成同尺寸四边形网格；
(3)提取四边形网格的节点坐标信息，利用步骤(1)得到的函数关系将节点向临近推进一个、两个或者四个节点，生成三维建模软件所需的节点坐标格式文件；
(4)在三维建模软件中依次导入节点坐标信息，将坐标点按顺序相连成线，再进行拉伸成面，最后进行面加厚处理；
(5)对三维建模软件建好的三维模型进行数值计算，将计算得到的伴流分布与目标伴流进行对比，如有误差，则细微调整三维建模软件所需的节点坐标，重复步骤(4)～(5)，直至得到精度满足要求的伴流模型；
(6)采用3D打印技术打印伴流模型，以获得低噪声目标伴流场。
2.根据权利要求1所述的针对非均匀来流推进器性能试验的轴向伴流模拟方法，其特征在于，函数关系的限制条件为：b+d的值为0.0075D～0.02D，伴流模型与推进器转子盘面的距离为0.02D～2D，其中，D为推进器转子直径。
3.根据权利要求1所述的针对非均匀来流推进器性能试验的轴向伴流模拟方法，其特征在于，步骤(2)的具体过程为：
将船艉目标伴流等值图放入作图软件中，缩放至实验用大小，描绘出等值边界线，接着提取等值线几何坐标信息，以此作为控制孔厚比的初始区域边界；
将初始区域边界向伴流区域内平移0.5mm，构成新的区域边界；然后，采用四边型网格对新区域边界内的伴流区域进行划分，初始网格划分尺寸为b+d＝0.01D；对网格节点进行编号并提取几何坐标信息保存。
4.根据权利要求1所述的针对非均匀来流推进器性能试验的轴向伴流模拟方法，其特征在于，步骤(3)中，将节点向临近推进一个、两个或者四个节点的依据为：
当节点1位于新区域边界线的交点时，其坐标点记为(x1,y1)，将该节点向临近推进一个节点，生成新的一个节点1-1，坐标点记为(x1-1,y1-1)；
当节点2位于新区域边界线上的网格节点处，其坐标点记为(x2,y2)，将该节点向临近推进两个节点，生成新的节点2-1和节点2-2，对应的坐标点分别记为(x2-1,y2-1)、(x2-2,y2-2)；
当节点3位于伴流区域内网格节点处，其坐标点记为(x3,y3)，将该节点向临近推进四个节点，生成新的节点3-1、节点3-2、节点3-3和节点3-4，对应的坐标点分别记为(x3-1,y3-1)、(x3-2,y3-2)、(x3-3,y3-3)、(x3-4,y3-4)；
当节点4位于伴流分数为wx,b的区域内，但等值边界曲线穿过其所在的网格，其坐标点记为(x4,y4)，该节点向临近推进四个节点，生成新的节点4-1、节点4-2、节点4-3和节点4-4，对应的坐标点分别记为(x4-1,y4-1)、(x4-2,y4-2)、(x4-3,y4-3)、(x4-4,y4-4)。
5.根据权利要求4所述的针对非均匀来流推进器性能试验的轴向伴流模拟方法，其特征在于，节点1-1与节点1的坐标关系为：
x1-1＝x1+dx y1-1＝y1+dy
其中，dx＝0.5ba，dy＝0.5ba，
Figure FDA0004212967370000031
da＝λaba；
该区域伴流分数为wx,a，对应的孔厚比为λa，则网格孔径为da，网格壁厚为ba。
6.根据权利要求5所述的针对非均匀来流推进器性能试验的轴向伴流模拟方法，其特征在于，节点2-1、节点2-2与节点2的坐标关系分别为：
x2-1＝x2+dx y2-1＝y2-dy
x2-2＝x2+dx y2-2＝y2+dy
式中，dx和dy均为0.5ba。
7.根据权利要求6所述的针对非均匀来流推进器性能试验的轴向伴流模拟方法，其特征在于，节点3-1、节点3-2、节点3-3、节点3-4与节点3的坐标关系分别为：
x3-1＝x3-dx y3-1＝y3-dy
x3-2＝x3-dxy 3-2＝y3+dy
x3-3＝x3+dx y3-3＝y3+dy
x3-4＝x3+dx y3-4＝y3-dy
式中，dx和dy均为0.5ba。
8.根据权利要求7所述的针对非均匀来流推进器性能试验的轴向伴流模拟方法，其特征在于，节点4-1、4-2与节点4的坐标关系同节点3-1、3-2与节点3的坐标关系；对于节点4-3和节点4-4，首先判定节点4与等值边界曲线的最短直线距离L与初始网格划分尺寸的关系：
如果L＞0.005D，则节点4-3、节点4-4与节点4的坐标关系同节点3-3、节点3-4与节点3的坐标关系，dx＝dy＝0.5bb；bb为伴流分数为wx,b区域的网格壁厚；
如果L＜0.005D，则节点4-3、节点4-4与节点4的坐标关系同节点3-3、节点3-4与节点3的坐标关系，但dx＝dy＝0.5bc；bc为伴流分数为wx,c区域的网格壁厚；
如果L＝0.005D，则节点4-3、节点4-4与节点4的坐标关系同节点3-3、节点3-4与节点3的坐标关系，但dx＝dy＝(0.5bb+0.5bc)/2。
CN202211190018.3A 2022-09-28 2022-09-28 一种针对非均匀来流推进器性能试验的轴向伴流模拟方法 Active CN115481583B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211190018.3A CN115481583B (zh) 2022-09-28 2022-09-28 一种针对非均匀来流推进器性能试验的轴向伴流模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211190018.3A CN115481583B (zh) 2022-09-28 2022-09-28 一种针对非均匀来流推进器性能试验的轴向伴流模拟方法

Publications (2)

Publication Number Publication Date
CN115481583A CN115481583A (zh) 2022-12-16
CN115481583B true CN115481583B (zh) 2023-06-23

Family

ID=84394123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211190018.3A Active CN115481583B (zh) 2022-09-28 2022-09-28 一种针对非均匀来流推进器性能试验的轴向伴流模拟方法

Country Status (1)

Country Link
CN (1) CN115481583B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112464585A (zh) * 2020-11-18 2021-03-09 西北工业大学 一种基于粘流和势流耦合的导管螺旋桨计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110920845B (zh) * 2019-11-14 2020-11-10 浙江大学 一种带c形导叶的全导管式双级吊舱推进器
CN112284681B (zh) * 2020-10-20 2022-10-11 中国船舶工业集团公司第七0八研究所 一种空泡水筒桨模试验用伴流场模拟方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112464585A (zh) * 2020-11-18 2021-03-09 西北工业大学 一种基于粘流和势流耦合的导管螺旋桨计算方法

Also Published As

Publication number Publication date
CN115481583A (zh) 2022-12-16

Similar Documents

Publication Publication Date Title
Pendar et al. Cavitation characteristics around a sphere: An LES investigation
CN104298869B (zh) 一种弹性水翼的流固耦合特性数值预测方法
CN111506968B (zh) 一种基于bp神经网络算法的船型优化方法
CN109685848A (zh) 一种三维点云与三维传感器的神经网络坐标变换方法
CN110118639B (zh) 一种三维数字船模虚拟阻力实验方法及系统
CN109992889B (zh) 风电场模型的建立方法及系统、尾流值计算方法及系统
CN109583377B (zh) 一种管路模型重建的控制方法、装置及上位机
CN116341421B (zh) 高超声速流场数值模拟方法、系统、电子设备及存储介质
CN106886666B (zh) 一种不同攻角下水下滑翔机升阻比的边界层网格求解方法
CN115481583B (zh) 一种针对非均匀来流推进器性能试验的轴向伴流模拟方法
CN109829210A (zh) 一种基于cfd数值计算的s型皮托管系数标定方法
Gac A large eddy based lattice-Boltzmann simulation of velocity distribution in an open channel flow with rigid and flexible vegetation
CN111079326B (zh) 二维各向异性网格单元度量张量场光滑化方法
CN113779695B (zh) 一种螺旋桨推进性能获取方法及其应用
CN112284681B (zh) 一种空泡水筒桨模试验用伴流场模拟方法
CN110968971A (zh) 一种实尺度的舰船声纳导流罩空化数值预报方法
CN107491574A (zh) 一种车辆发动机进气道仿真计算的方法及装置
Zhang et al. Converging near-elliptic shock waves
CN109614757B (zh) 一种通过cfd预测湿室型泵站临界淹没深度的方法
CN114611214B (zh) 基于螺旋桨高效建模的敞水性能仿真方法
Kravtsova et al. Cavitation on a semicircular leading-edge plate and NACA0015 hydrofoil: visualization and velocity measurement
CN108804791B (zh) 一种适用于埋入式进气道布局的飞行器参数化方法
CN113536640A (zh) 一种基于正交试验的布风器内部流道结构的优化设计方法
Li et al. Deep learning-based reduced order model for three-dimensional unsteady flow using mesh transformation and stitching
CN115659689B (zh) 反映实际边界条件的泵装置内流及性能数值迭代计算方法

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