CN106227975A - 海洋条件下流型实验辨识方法 - Google Patents
海洋条件下流型实验辨识方法 Download PDFInfo
- Publication number
- CN106227975A CN106227975A CN201610632938.4A CN201610632938A CN106227975A CN 106227975 A CN106227975 A CN 106227975A CN 201610632938 A CN201610632938 A CN 201610632938A CN 106227975 A CN106227975 A CN 106227975A
- Authority
- CN
- China
- Prior art keywords
- flow pattern
- bubble
- flow
- value
- identification method
- 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
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
Abstract
本发明公开了海洋条件下流型实验辨识方法,包括以下步骤:1)、获得无量纲气泡弦长Cb *的概率密度函数f(x);2)、将概率密度函数f(x)与各个典型无量纲气泡弦长的PDF分布图进行对比判定流型。本发明通过建立不同流型的标准PDF图,将实验流型所测得的概率密度函数f(x)的与标准PDF图进行不交判断所属流型,提高了海洋流型判断的准确性。
Description
技术领域
本发明涉及海洋流型辨识技术领域,具体涉及海洋条件下流型实验辨识方法。
背景技术
当船舶在海洋条件下航行时,会受海水和风浪的作用处于倾斜、起伏和摇摆等海洋运动条件,并由此产生了一些不同于静止条件下的附加力,如切向力、离心力,科氏力;同时还伴随着浮力在流动方向上的变化,这些力的共同作用会影响通道内流型形态及演化进而导致两相流动、传热特性发生变化。因此,对海洋条件下两相流型进行准确判定对于准确计算海洋条件下两相流动、传热特性等具有重要意义。
目前,静止条件下的流型判定大多采用可视化的方法,其优点在于可以直接观察流型之间的转变,缺点在于主观性较大,存在较大不确定性,同时由于可视窗的存在,在中高压条件下不适用。
发明内容
本发明的目的在于提供海洋条件下流型实验辨识方法,解决现有的辨识方法导致的流型判断不准确的问题,本发明所述的辨识方法还对精确计算流动压降、传热系数等具有有重要价值。
本发明通过下述技术方案实现:
海洋条件下流型实验辨识方法,包括以下步骤:
1)、获得无量纲气泡弦长Cb *的概率密度函数f(x);
2)、将概率密度函数f(x)与各个典型无量纲气泡弦长的PDF分布图进行对比判定流型。
本发明通过建立不同流型的标准PDF图,将实验流型所测得的概率密度函数f(x)的与标准PDF图进行不交判断所属流型,提高了海洋流型判断的准确性。
进一步地,PDF分布图的制备过程为:对于每个典型流型重复步骤1)的操作,对步骤1)的结果进行优化筛选获得PDF分布图。
所述优化筛选具体是指将出现概率较多的概率密度函数f(x)作为标准PDF分布图。
进一步地,步骤1)中采集气泡数量大于2000个,且采集周期大于3倍运动周期。
实验证明,采用上述气泡采集方法,气泡弦长的PDF(概率密度)分布结果将不受气泡数量和采集时间的影响。
进一步地,概率密度函数f(x)的获得包括以下步骤:
a)、用双探头光学探针测得气泡的轴向运动速度Vb,计算出t2j-1时刻的气泡弦长Cb,Cb的计算方式为:Cb=Vb(t2j-t2j-1);
b)、计算出t2j-1时刻的无量纲气泡弦长Cb *,Cb *的计算方式为:Cb *=Cb/D,其中D为通道水力直径;
c)、f(x)的计算方式为:f(xi)=Pi/Δx;式中,Pi=ni/N;Δx为无量纲气泡弦长Cb *分组组距,Pi为Cb *在第i组的概率,ni为Cb *出现于第i组的气泡个数,N为总的气泡个数。
其中,vb为轴向气泡速度,可用双探头光学探针同时测得,气泡弦长Cb是指气泡运动过程中穿过探针的长度,对于t2j-1时刻穿过探针的气泡。
进一步地,还包括以下步骤:
e)、获得不同流型的标准Cb/D值:对于每个流型重复步骤a)、b)操作,获得若干Cb/D值,选取相对集中的Cb/D的范围值作为标准Cb/D值;
f)、将实验流型测得的Cb/D值标准Cb/D值进行比较判断流型类型。
进一步地,不同流型包括弥散泡状流、泡状流和弹状流,其中,弥散泡状流的标准的Cb/D值为小于等于0.2;泡状流的标准的Cb/D值为0.2-1.5;弹状流的标准的Cb/D值为大于等于1.5。
对于本领域技术人员,通常指针对典型流型进行识别,还没有想到对临界区域的流型进行较为准确的辨识,也不容易想到通过测量Cb/D值对散泡状流、泡状流和弹状流,上述辨识方法以及标准Cb/D值是申请人通过大量实验加以分析获得。
具体地,利用本海洋条件下流型实验辨识方法开展海洋条件下两相摩阻及传热实验时,通过双探头光学探针采集气泡数据,利用后处理软件得到无量纲气泡弦长的PDF分布函数,再根据获得的气泡弦长PDF分布图对海洋条件下流型进行定量识别。
已利用本方法对海洋条件下两相摩阻及传热实验开展研究,实验结果显示本方法能准确辨识海洋条件下弥散泡状流、泡状流和弹状流。
本发明与现有技术相比,具有如下的优点和有益效果:
1、本发明通过建立不同流型的标准PDF图,将实验流型所测得的概率密度函数f(x)的与标准PDF图进行不交判断所属流型,提高了海洋流型判断的准确性。
2、本发明通过建立不同流型的标准Cb/D值,不但可对典型流型进行识别,还可较为精确对临界区域的流型进行辨识。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1是探针单探头测量信号示意图;
图2是弥散泡状流的PDF图;
图3是泡状流的PDF图;
图4是弹状流的PDF图;
图5是海洋条件下流型判定流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
在图1中,气相穿过探针时输出信号为高电平,液相穿过探针时输出信号为低电平。在图2-至图4中,图中PDF分布曲线将随横坐标的减小而先骤升后骤降并趋进于0,但由于该区域对应于图中非常小的横坐标区,因此该趋势在PDF分布图中无法直观表示出来。
实施例1:
当实验装置达到设定海洋运动及热工工况并稳定后,通过双探头光学探针采集不少于2000个气泡,且光学探针采集时间不低于3倍运动周期。
如图1至图5所示,海洋条件下流型实验辨识方法,包括以下步骤:
1)、获得无量纲气泡弦长Cb *的概率密度函数f(x);
概率密度函数f(x)的获得包括以下步骤:
a)、用双探头光学探针测得气泡的轴向运动速度Vb,计算出t2j-1时刻的气泡弦长Cb,Cb的计算方式为:Cb=Vb(t2j-t2j-1);
b)、计算出t2j-1时刻的无量纲气泡弦长Cb *,Cb *的计算方式为:Cb *=Cb/D,其中D为通道水力直径;
c)、f(x)的计算方式为:f(xi)=Pi/Δx;式中,Pi=ni/N;Δx为无量纲气泡弦长Cb *分组组距,Pi为Cb *在第i组的概率,ni为Cb *出现于第i组的气泡个数,N为总的气泡个数;
2)、将概率密度函数f(x)与各个典型无量纲气泡弦长的PDF分布图进行对比判定流型;所述PDF分布图的制备过程为:对于每个典型流型重复步骤1)的操作,对步骤1)的结果进行优化筛选获得PDF分布图。
在本实施例中,通过将获得的概率密度函数f(x)与各个典型无量纲气泡弦长的PDF分布图进行比较,能有准确的对典型流型进行识别,相比现有的可视化方法提高了流型识别的准确性。
实施例2:
如图1至图5所示,本实施例基于实施例1,还包括以下步骤:
e)、获得不同流型的标准Cb/D值:对于每个流型重复步骤a)、b)操作,获得若干Cb/D值,选取相对集中的Cb/D的范围值作为标准Cb/D值;
f)、将实验流型测得的Cb/D值标准Cb/D值进行比较判断流型类型;
不同流型包括弥散泡状流、泡状流和弹状流,其中,弥散泡状流的标准的Cb/D值为小于等于0.2;泡状流的标准的Cb/D值为0.2-1.5;弹状流的标准的Cb/D值为大于等于1.5。
海洋条件下Cb/D全都分布在≤0.2时为弥散泡状流,有0.2<Cb/D<1.5的气泡分布时为泡状流,而当有Cb/D≥1.5的气泡出现时为弹状流。
在本实施例中利用建立的标准Cb/D值,可实现海洋条件下弥散泡状流、泡状流和弹状流准确识别,对精确计算流动压降、传热系数等有重要价值。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.海洋条件下流型实验辨识方法,其特征在于,包括以下步骤:
1)、获得无量纲气泡弦长Cb *的概率密度函数f(x);
2)、将概率密度函数f(x)与各个典型无量纲气泡弦长的PDF分布图进行对比判定流型。
2.根据权利要求1所述的海洋条件下流型实验辨识方法,其特征在于,所述PDF分布图的制备过程为:对于每个典型流型重复步骤1)的操作,对步骤1)的结果进行优化筛选获得PDF分布图。
3.根据权利要求1所述的海洋条件下流型实验辨识方法,其特征在于,步骤1)中采集气泡数量大于2000个,且采集周期大于3倍运动周期。
4.根据权利要求1所述的海洋条件下流型实验辨识方法,其特征在于,概率密度函数f(x)的获得包括以下步骤:
a)、用双探头光学探针测得气泡的轴向运动速度Vb,计算出t2j-1时刻的气泡弦长Cb,Cb的计算方式为:Cb=Vb(t2j-t2j-1);
b)、计算出t2j-1时刻的无量纲气泡弦长Cb *,Cb *的计算方式为:Cb *=Cb/D,其中D为通道水力直径;
c)、f(x)的计算方式为:f(xi)=Pi/Δx;式中,Pi=ni/N;Δx为无量纲气泡弦长Cb *分组组距,Pi为Cb *在第i组的概率,ni为Cb *出现于第i组的气泡个数,N为总的气泡个数。
5.根据权利要求4所述的海洋条件下流型实验辨识方法,其特征在于,还包括以下步骤:
e)、获得不同流型的标准Cb/D值:对于每个流型重复步骤a)、b)操作,获得若干Cb/D值,选取相对集中的Cb/D的范围值作为标准Cb/D值;
f)、将实验流型测得的Cb/D值标准Cb/D值进行比较判断流型类型。
6.根据权利要求5所述的海洋条件下流型实验辨识方法,其特征在于,不同流型包括弥散泡状流、泡状流和弹状流,其中,弥散泡状流的标准的Cb/D值为小于等于0.2;泡状流的标准的Cb/D值为0.2-1.5;弹状流的标准的Cb/D值为大于等于1.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610632938.4A CN106227975A (zh) | 2016-08-04 | 2016-08-04 | 海洋条件下流型实验辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610632938.4A CN106227975A (zh) | 2016-08-04 | 2016-08-04 | 海洋条件下流型实验辨识方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106227975A true CN106227975A (zh) | 2016-12-14 |
Family
ID=57548221
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610632938.4A Pending CN106227975A (zh) | 2016-08-04 | 2016-08-04 | 海洋条件下流型实验辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106227975A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106653115A (zh) * | 2017-01-09 | 2017-05-10 | 中国核动力研究设计院 | 一种运动条件下的棒束通道压差测量组件 |
CN111415348A (zh) * | 2020-03-25 | 2020-07-14 | 中国计量大学 | 一种汽车制动管道内气泡特征的提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005098797A (ja) * | 2003-09-24 | 2005-04-14 | Keio Gijuku | 単一上昇気泡の3次元形状及び上昇速度の測定法及び測定装置 |
CN202599852U (zh) * | 2012-04-23 | 2012-12-12 | 长安大学 | 气液两相泡状流中气泡的识别及测量装置 |
-
2016
- 2016-08-04 CN CN201610632938.4A patent/CN106227975A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005098797A (ja) * | 2003-09-24 | 2005-04-14 | Keio Gijuku | 単一上昇気泡の3次元形状及び上昇速度の測定法及び測定装置 |
CN202599852U (zh) * | 2012-04-23 | 2012-12-12 | 长安大学 | 气液两相泡状流中气泡的识别及测量装置 |
Non-Patent Citations (1)
Title |
---|
孙奇; 赵华; 杨瑞昌; 蒋序伦; 胡军;: "双探头光学探针识别受热工况两相流流型的基本方法研究", 《核动力工程》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106653115A (zh) * | 2017-01-09 | 2017-05-10 | 中国核动力研究设计院 | 一种运动条件下的棒束通道压差测量组件 |
CN106653115B (zh) * | 2017-01-09 | 2018-01-02 | 中国核动力研究设计院 | 一种运动条件下的棒束通道压差测量组件 |
CN111415348A (zh) * | 2020-03-25 | 2020-07-14 | 中国计量大学 | 一种汽车制动管道内气泡特征的提取方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nakoryakov et al. | Study of local hydrodynamic characteristics of upward slug flow | |
CN101696925B (zh) | 一种两相流减阻剂性能的测试装置及测试方法 | |
CN106226392A (zh) | 基于超声衰减机理模型的油水两相流相含率测量方法 | |
CN109253765A (zh) | 河流流量在线监测测量系统及流量计算方法 | |
Hildebrandt et al. | Pressure distribution and vortex shedding around a cylinder due to a steep wave at the onset of breaking from physical and numerical modeling | |
CN106227975A (zh) | 海洋条件下流型实验辨识方法 | |
Gopalan et al. | Effect of gap size on tip leakage cavitation inception, associated noise and flow structure | |
Fox et al. | Eddy taxonomy methodology around a submerged barb obstacle within a fixed rough bed | |
CN108398236A (zh) | 斜管式沉沙池流场信息及泥沙运移特性测量方法与系统 | |
Huguenard et al. | Intermittency in Estuarine Turbulence: A framework toward limiting bias in microstructure measurements | |
Schreck et al. | Unsteady separation processes and leading edge vortex precursors: pitch rate and Reynolds number influences | |
CN108412481A (zh) | 超声波多普勒多相流分相流量测井资料解释方法和装置 | |
CN105891339B (zh) | 利用冲击映像法检测岩土工程介质缺陷的方法、装置及系统 | |
Cui et al. | Characteristics of wake morphology during debris flow when passing a cylindrical obstacle | |
Mouri et al. | Vortex tubes in velocity fields of laboratory isotropic turbulence: dependence on the Reynolds number | |
Tedds | Scale model testing of tidal stream turbines: Wake characterisation in realistic flow conditions | |
Freitag et al. | Mooring motion bias of point-Doppler current meter measurements | |
CN110388899A (zh) | 用于涌潮流速垂向分布计算的专用组件 | |
Evans et al. | Wake characteristics of a natural submerged pinnacle and implications for tidal stream turbine installations | |
Hassanzadeh et al. | Instrumentation development and testing of a wind turbine blade for sub-scale wake studies | |
CN104459699A (zh) | 一种基于目标航深的水面与水下目标分类识别方法 | |
Lien et al. | Skin Friction Determination Using Multihole Pressure Probes | |
CN109131771A (zh) | 一种基于海事监管的内河船舶吃水检测系统及检测方法 | |
Holand et al. | Identification of wave breaking from nearshore wave-by-wave records | |
Pan et al. | Investigating to cavitation behavior of orifice tunnel |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20161214 |