CN101082561B - 测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置 - Google Patents
测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置 Download PDFInfo
- Publication number
- CN101082561B CN101082561B CN2007100698634A CN200710069863A CN101082561B CN 101082561 B CN101082561 B CN 101082561B CN 2007100698634 A CN2007100698634 A CN 2007100698634A CN 200710069863 A CN200710069863 A CN 200710069863A CN 101082561 B CN101082561 B CN 101082561B
- Authority
- CN
- China
- Prior art keywords
- particle
- level crossing
- stereoscope
- dimensional
- field
- 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.)
- Expired - Fee Related
Links
- 239000002245 particle Substances 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 14
- 239000007787 solid Substances 0.000 title claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 9
- 230000004907 flux Effects 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims abstract description 3
- 230000001052 transient effect Effects 0.000 claims abstract description 3
- 208000003164 Diplopia Diseases 0.000 claims description 23
- 208000029444 double vision Diseases 0.000 claims description 23
- 230000000007 visual effect Effects 0.000 claims description 10
- 239000004744 fabric Substances 0.000 claims description 2
- 230000005514 two-phase flow Effects 0.000 abstract description 4
- 238000003909 pattern recognition Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000000917 particle-image velocimetry Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
Images
Landscapes
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置。将粒子流发射管放置在由双视角立体镜和高速摄影机组成的高速立体视觉系统前，射流粒子经过双视角立体镜镜面折射，单个高速摄影机靶面可同时获得两个不同角度的粒子流图像，图像输入计算机后，由模式识别算法程序和三维逆向投影算法程序组成的程序模块对一幅粒子立体图像进行处理，得到粒子的瞬态三维浓度场，程序模块处理一系列微秒级时序变化粒子图像，得到粒子三维运动轨迹以及三维速度场。使用单相机系统稳定，克服了多相机系统结构复杂、同步性差的问题。
Description
技术领域
本发明涉及气固两相流体三维运动特性分析技术，具体涉及一种测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置。
背景技术
对气固两相流的研究是湍流研究领域中重要而前沿的课题，目前对于气固两相流中固体粒颗粒运动轨迹的一个主要研究方向是两相流测试技术，其中粒子图像测速PIV技术以及DPIV技术由于非接触全场测量、原理简单、实现方便等优点而受到广泛的应用与研究。但以上技术一般均采用片光源和单镜头摄影系统，只能得到两相流中固体粒子三维运动的二维投影图像，缺乏立体感，因而只能用于二维流场或者三维流场中某一断面粒子运动特征的研究。为了得到真实的三维运动图像，可以采用多个相机，从不同角度拍摄，该方法问题是系统复杂、设备投资大，不同相机之间同步性难以掌握，多相机系统参数无法固定，应用此方法重建粒子的三维空间运动轨迹有很大的困难。
发明内容
本发明的目的在于提供一种使用双视角立体镜和单个高速摄影机组成的高速立体视觉系统分析、测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置。
本发明解决其技术问题所采用的技术方案是：
一、一种测量气固两相流中固体颗粒三维浓度场、速度场的方法：
将粒子流发射管放置在由双视角立体镜和高速摄影机组成的高速立体视觉系统前，射流粒子经过双视角立体镜镜面折射，单个高速摄影机靶面可同时获得两个不同角度的粒子流图像，图像输入计算机后，由模式识别算法程序和三维逆向投影算法程序组成的程序模块对一幅粒子立体图像进行处理，得到粒子的瞬态三维浓度场，程序模块处理一系列微秒级时序变化粒子图像，得到粒子三维运动轨迹以及三维速度场。
二、一种测量气固两相流中固体颗粒三维浓度场、速度场的装置：
分为硬件系统和软件系统两个部分，硬件系统主要包括三个部分：固定于高速摄影机镜头前的双视角立体镜、用于采集图像的高速摄影机和对图像进行处理的后台计算机，所述的双视角立体镜由四片平面反射镜按如下方法放置排列，即第一平面镜、第二平面镜、第三平面镜和第四平面镜自左至右依次排列，第一平面镜和第二平面镜平行排列，第三平面镜和第四平面镜平行排列，第一平面镜和第二平面镜与第三平面镜和第四平面镜两者对称放置，第二平面镜与第三平面镜之间的夹角为90度，四片平面反射镜镜片固定在镜架上，双视角立体镜整体固定在摄影机镜头上；软件系统采用Visual Studio 6.0和MATLAB编程软件开发。
本发明具有的有益效果是：
单个高速摄影机靶面上同时获取两个不同角度的粒子图像、通过粒子模式识别技术和逆向投影算法得到流场中粒子的三维坐标是分析气固两相流的最佳方案。此方法有以下优点：
1.使用单相机系统同时获取不同角度粒子运动图像，克服了多相机系统结构复杂、同步性差的问题。
2.单相机系统稳定，对粒子三维运动轨迹的重建计算提供稳定的参数。
3.基于亚象素拟合的高精度粒子模式识别算法大大提高了不同角度图像中对应两相粒子的识别精度，通过两个不同角度逆向投影面的交汇，可以得到图像中粒子的三维坐标。
4.结合微秒级时序坐标变化，可以对两相流中粒子的三维速度场、浓度场进行分析研究。
附图说明
图1是高速立体视觉系统结构示意图。
图2是双视角立体镜光路图。
图3是双视角立体镜结构主视图。
图4是图3仰视图。
图5是图像处理软件系统流程图。
图中：1、粒子流发射管，2、双视角立体镜，3、高速摄影机，4、平面反射镜，5、平面反射镜，6、平面反射镜，7、平面反射镜，8、镜头，9、靶面，10、镜架，11、双视角立体镜透光孔。
具体实施方式
如图1、图2、图3、图4所示，本发明分为硬件系统和软件系统两个部分，硬件系统主要包括三个部分：固定于高速摄影机镜头前的双视角立体镜2、用于采集图像的高速摄影机3和对图像进行处理的后台计算机，所述的双视角立体镜2由四片平面反射镜按如下方法放置排列，即第一平面镜4、第二平面镜5、第三平面镜6和第四平面镜7自左至右依次排列，第一平面镜4和第二平面镜5平行排列，第三平面镜6和第四平面镜7平行排列，第一平面镜4和第二平面镜5与第三平面镜6和第四平面镜7两者对称放置，第二平面镜5与第三平面镜6之间的夹角为90度，四片平面反射镜镜片固定在镜架10上，双视角立体镜整体固定在摄影机镜头上；软件系统采用Visual Studio6.0和MATLAB编程软件开发。
如图1、图2、图3、图4所示，应用立体视觉系统拍摄粒子射流，双视角立体镜2水平固定在高速摄影机3镜头前。高速立体视觉系统的软件系统流程图如图5所示。首先使用标准网格图像和基于MATLAB平台开发的CCD参数标定算法建立CCD成像光路模型，获取CCD光学投影矩阵。之后，使用高速立体视觉系统拍摄粒子射流，射流粒子经过图2所示的四个平面反射镜4、5、6、7所组成的双视角立体镜2光路，通过镜架10上的双视角立体镜的透光孔11及镜头8，在摄影机CCD靶面9上投影得到两个不同角度的立体图像即双视角粒子运动图像，图像传入计算机，在得到CCD光学投影矩阵的基础上，对单视角粒子运动图像进行处理，可以得到粒子的二维图像归一化坐标。使用基于VisualStudio6.0平台开发的亚像素模式识别算法对双视角粒子运动图像进行处理，结合已得到的光学投影模型计算出粒子的三维径深即粒子在三维空间中的深度信息。在获得粒子的二维图像归一化坐标和其三维径深信息后，使用基于VisualStudio6.0平台开发的逆向投影算法，计算得到粒子三维坐标。结合微秒级粒子图像序列，可以重建粒子的三维运动轨迹，获得粒子的三维速度场，以及粒子团的浓度场变化。
Claims (2)
1.一种测量气固两相流中固体颗粒三维浓度场、速度场的方法，其特征在于：将粒子流发射管放置在由双视角立体镜和高速摄影机组成的高速立体视觉系统前，射流粒子图像经过双视角立体镜镜面多次反射，在单个高速摄影机靶面同时获得两个不同角度的粒子流图像，图像输入计算机后，由模式识别算法程序和三维逆向投影算法程序组成的程序模块对一幅粒子立体图像进行处理，得到粒子的瞬态三维浓度场，程序模块处理一系列微秒级时序变化粒子图像，得到粒子三维运动轨迹以及三维速度场。
2.用于权利要求1所述方法的一种测量气固两相流中固体颗粒三维浓度场、速度场的装置，其特征在于：分为硬件系统和软件系统两个部分，硬件系统主要包括三个部分：固定于高速摄影机镜头前的双视角立体镜(2)、用于采集图像的高速摄影机(3)和对图像进行处理的后台计算机，所述的双视角立体镜(2)由四片平面反射镜镜面相对，按如下方法放置排列，即第一平面镜(4)、第二平面镜(5)、第三平面镜(6)和第四平面镜(7)自左至右依次排列，第一平面镜(4)和第二平面镜(5)平行排列，第三平面镜(6)和第四平面镜(7)平行排列，第一平面镜(4)和第二平面镜(5)与第三平面镜(6)和第四平面镜(7)两者对称放置，第二平面镜(5)与第三平面镜(6)之间的夹角为90度，四片平面反射镜镜片固定在镜架(10)上，双视角立体镜整体固定在摄影机镜头上；软件系统采用Visual Studio 6.0和MATLAB编程软件开发。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2007100698634A CN101082561B (zh) | 2007-07-03 | 2007-07-03 | 测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2007100698634A CN101082561B (zh) | 2007-07-03 | 2007-07-03 | 测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101082561A CN101082561A (zh) | 2007-12-05 |
| CN101082561B true CN101082561B (zh) | 2011-06-29 |
Family
ID=38912246
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2007100698634A Expired - Fee Related CN101082561B (zh) | 2007-07-03 | 2007-07-03 | 测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101082561B (zh) |
Families Citing this family (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100585360C (zh) * | 2008-01-22 | 2010-01-27 | 南京航空航天大学 | 飞行载荷作用下两相流地面模拟测试方法及装置 |
| CN101666702B (zh) * | 2009-09-21 | 2011-06-08 | 浙江大学 | 流浆箱阶梯扩散器性能测试试验装置 |
| CN102175178B (zh) * | 2011-02-18 | 2012-10-31 | 华南理工大学 | 运动火源扩散火焰锋面三维结构测量系统与方法 |
| CN102494869B (zh) * | 2011-12-27 | 2013-11-27 | 东南大学 | 稠密两相流颗粒速度和浓度空间分布的测量装置 |
| CN102620903A (zh) * | 2012-03-22 | 2012-08-01 | 河海大学 | 一种水流中流迹线的动态显示系统及其显示方法 |
| CN103308276B (zh) * | 2013-06-17 | 2015-08-26 | 内蒙古科技大学 | 高压供粉装置 |
| CN104457703A (zh) * | 2014-11-28 | 2015-03-25 | 浙江大学 | 基于多视觉小通道气液两相流参数测量系统及方法 |
| CN105258919B (zh) * | 2015-11-24 | 2018-11-13 | 浙江工业大学 | 多粒径混合气固两相流动实验装置 |
| CN107219057B (zh) * | 2016-03-22 | 2019-04-09 | 中国石油化工股份有限公司 | 氧气混合设备流体力学测试方法 |
| CN106018280B (zh) * | 2016-07-19 | 2018-12-14 | 江苏大学 | 一种同时测量速度场和浓度场的装置和方法 |
| CN106526223A (zh) * | 2016-10-19 | 2017-03-22 | 江苏理工学院 | 气相和固相速度场同步测量装置和方法 |
| CN107389489B (zh) * | 2017-07-31 | 2019-08-20 | 东南大学 | 一种气固混合物中颗粒物浓度的测量装置及测量方法 |
| CN108645761B (zh) * | 2018-05-30 | 2019-04-16 | 西安科技大学 | 测试粉尘颗粒运动特征及参数的可视化系统及方法 |
| CN109269952A (zh) * | 2018-10-16 | 2019-01-25 | 东南大学 | 一种丝状物料在气流床矩形提升管中三维颗粒浓度信息的测量装置和方法 |
| CN109883906B (zh) * | 2019-02-21 | 2020-08-21 | 浙江大学 | 一种纳米金属两相流体稳定性测量方法 |
| CN110879301A (zh) * | 2019-10-11 | 2020-03-13 | 中国航发沈阳发动机研究所 | 一种同时测量液体浓度和液体运动二维分布的方法及系统 |
| CN112730163B (zh) * | 2020-12-24 | 2022-09-27 | 昆明理工大学 | 一种散体内部颗粒运动轨迹的观测方法 |
| CN113109020A (zh) * | 2021-03-18 | 2021-07-13 | 上海交通大学 | 一种示踪方法及系统 |
| CN115933109A (zh) * | 2022-11-29 | 2023-04-07 | 哈尔滨工程大学 | 一种适用于多相机piv试验中相机跨介质倾斜拍摄的光学辅助系统 |
| CN116380408B (zh) * | 2023-04-10 | 2024-06-07 | 南京航空航天大学 | 基于结构光和光场成像的三维超分辨流场测量方法及系统 |
| CN116499394A (zh) * | 2023-04-24 | 2023-07-28 | 山东大学 | 一种气泡三维可视化成像装置及气泡三维重构方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5905568A (en) * | 1997-12-15 | 1999-05-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Stereo imaging velocimetry |
-
2007
- 2007-07-03 CN CN2007100698634A patent/CN101082561B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5905568A (en) * | 1997-12-15 | 1999-05-18 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Stereo imaging velocimetry |
Non-Patent Citations (3)
| Title |
|---|
| 周昊等.《基于光学波动法研究气固两相流扩散特性》.《燃烧科学与技术》.2003,第9卷(第4期),291-294. * |
| 夏钧等.《气固两相湍流射流的直接数值模拟》.《自然科学进展》.2003,第13卷(第5期),528-532. * |
| 王浩等.《二维及三维流场的光学测量方法》.《光学技术》.2001,第27卷(第2期),139-142. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101082561A (zh) | 2007-12-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101082561B (zh) | 测量气固两相流中固体颗粒三维浓度场、速度场的方法和装置 | |
| CN101072366B (zh) | 基于光场和双目视觉技术的自由立体显示系统和方法 | |
| CN106908622B (zh) | 一种基于光场成像的层析piv测量装置及方法 | |
| CN102435770B (zh) | 一种单相机三维体视粒子图像测速系统 | |
| JP5469483B2 (ja) | 流体解析装置および流体解析方法 | |
| CN103292695B (zh) | 一种单目立体视觉测量方法 | |
| CN103115613B (zh) | 一种空间三维定位方法 | |
| CN107525945B (zh) | 基于集成成像技术的3d-3c粒子图像测速系统及方法 | |
| CN109544679A (zh) | 管道内壁的三维重建方法 | |
| TWI520576B (zh) | 將二維影像轉換爲三維影像的方法與系統及電腦可讀媒體 | |
| CN102713513B (zh) | 摄像装置、摄像方法、程序以及集成电路 | |
| CN106500629B (zh) | 一种显微三维测量装置及系统 | |
| CN109285189B (zh) | 一种无需双目同步的直线轨迹快速计算方法 | |
| CN201072403Y (zh) | 测量气固两相流中固体颗粒三维浓度场、速度场的装置 | |
| CN107421462A (zh) | 基于线激光扫描的物体三维轮廓测量系统 | |
| CN104019763B (zh) | 流体三维速度场及变形体三维形貌的同步测量装置 | |
| CN107066975B (zh) | 基于深度传感器的视频识别与跟踪系统及其方法 | |
| CN110500959A (zh) | 一种单摄像头口内三维扫描系统 | |
| CN108010125A (zh) | 基于线结构光和图像信息的真实尺度三维重建系统及方法 | |
| CN106291788B (zh) | 自由曲面棱镜及其形状的确定方法及其光学成像方法 | |
| CN202256377U (zh) | 一种单相机三维体视粒子图像测速系统 | |
| CN110728745A (zh) | 一种基于多层折射图像模型的水下双目立体视觉三维重建方法 | |
| Grant et al. | The application of an in-line, stereoscopic, PIV system to 3-component velocity measurements | |
| Zhao et al. | Geometrical-analysis-based algorithm for stereo matching of single-lens binocular and multi-ocular stereovision system | |
| CN107274449B (zh) | 一种光学照片对物体的空间定位系统及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20071205 Assignee: HANGZHOU JINJIANG GROUP Co.,Ltd. Assignor: Zhejiang University Contract record no.: 2015330000163 Denomination of invention: Method for measuring solid particle three-dimensional concentration field and velocity field in gas/solid two-phase stream Granted publication date: 20110629 License type: Common License Record date: 20150616 |
|
| LICC | Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110629 |