CN106323184A - Method for three-dimensionally and dynamically measuring diameter of slipping bubble - Google Patents

Method for three-dimensionally and dynamically measuring diameter of slipping bubble Download PDF

Info

Publication number
CN106323184A
CN106323184A CN201611029639.8A CN201611029639A CN106323184A CN 106323184 A CN106323184 A CN 106323184A CN 201611029639 A CN201611029639 A CN 201611029639A CN 106323184 A CN106323184 A CN 106323184A
Authority
CN
China
Prior art keywords
image
bubble
sliding
point
picture
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
Application number
CN201611029639.8A
Other languages
Chinese (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.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
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 Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN201611029639.8A priority Critical patent/CN106323184A/en
Publication of CN106323184A publication Critical patent/CN106323184A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/08Measuring arrangements characterised by the use of optical techniques for measuring diameters
    • G01B11/10Measuring arrangements characterised by the use of optical techniques for measuring diameters of objects while moving

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention discloses a method for three-dimensionally and dynamically measuring diameter of a slipping bubble. The method comprises the following steps, shoot scale plates with scales by using two vertical high-speed cameras, and obtain a clear scale plate image; transfer effective pixel of the scale plate image to real physical dimension, and shoot a background picture when there is no bubble; shoot the slipping bubble generated at a bubble nucleating point from the front face and the side face; peel off the slipping bubble image; filter the peeled bubble image; perform the binarization processing on the filtered image, and obtain a black-and-white image; calculate equivalent diameter of the bubble in real time. Through two high-speed cameras, the slipping bubble is three-dimensionally shot from the front face and the side face, and the slipping bubble image is obtained; through analyzing and processing the collected image in real time, the dynamic change process of diameter during the bubble slipping process can be accurately measured, so as to lay a foundation for the subsequent movement behavior analysis of the slipping bubble and a boiling heat transfer mechanism model.

Description

A kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method
Technical field
The invention belongs to experimental technique, particularly to a kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method, for accurately Measure steam bubble dia in steam bubble slipping dynamically to change.
Background technology
In Narrow Rectangular Channel, the research of boiling heat transfer mechanism is for thermal technology's water in the reactor of employing Plate-type Fuel Elements Force system safety and economy have meaning, and sliding bubbles diameter is the foundation of boiling heat transfer mechanism model and benchmark analysis The middle key parameter that must input, it is therefore necessary in accurately measuring Narrow Rectangular Channel by laboratory facilities, sliding bubbles is straight Footpath.Existing experiment often uses high-speed camera instrument from just shooting in the face of steam bubble, and calculates steam bubble dia by the method for equal effect, but When steam bubble is constantly swung by turbulence prediction, its certainty of measurement is the highest.
The present invention is based on above-mentioned application background, proposes a kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method, can essence Really measuring steam bubble dia dynamic changing process in steam bubble slipping, the foundation for follow-up boiling heat transfer mechanism model provides basis Test data support.
Summary of the invention
In order to overcome the disadvantages mentioned above of prior art, the invention provides a kind of sliding bubbles diameter Three-Dimensional Dynamic measurement side Method, can obtain 3 d image data in steam bubble slipping, and accurately measures steam bubble dia in steam bubble slipping and dynamically change Process, the foundation for follow-up boiling heat transfer mechanism model provides basic test data supporting.
The technical solution adopted for the present invention to solve the technical problems is: a kind of sliding bubbles diameter Three-Dimensional Dynamic measurement side Method, comprises the steps:
Step one, utilize two vertically disposed high-speed camera instrument shooting scale with a scale, obtain scale map clearly Picture, is changed into actual physical size by the valid pixel of scale map picture;
Step 2, two the high-speed camera instrument utilizing step one to arrange shoot background picture when not having steam bubble;
Step 3, utilize two high-speed camera instrument that step one arranges respectively from front and side simultaneously to bubble nucleation point The sliding bubbles that position produces shoots;
Step 4, sliding bubbles image is peeled off;
Step 5, to peel off after steam bubble image be filtered;
Step 6, filtered image is carried out binary conversion treatment, it is thus achieved that black white image;
The equivalent diameter of step 7, in real time calculating steam bubble.
Compared with prior art, the positive effect of the present invention is:
The present invention carries out three-dimensional shooting from front and side to sliding bubbles by two high-speed camera instrument simultaneously, the most sharp Carry out synchronizing Real-time Collection to the image that two high-speed camera instrument shoot with image capturing system, and utilize the poster processing soft to adopting The image of collection is analyzed in real time and processes, and can accurately measure steam bubble dia dynamic changing process in steam bubble slipping, for rear Continuous sliding bubbles motor behavior analysis and boiling heat transfer mechanism model are set up and are laid the foundation.
Use the inventive method can obtain front and side image data in steam bubble slipping, can be according to the image obtained Data obtain steam bubble dia in steam bubble slipping and dynamically change.Steam bubble nucleation and the spy that develops in narrow channel can be applied to Property visualized experiment.
Accompanying drawing explanation
Examples of the present invention will be described by way of reference to the accompanying drawings, wherein:
Fig. 1 is the measurement system composition schematic diagram of the present invention;Wherein, 1, bubble nucleation point, 2,1# high-speed camera instrument, 3,2# High-speed camera instrument, 4, image variants system;
Fig. 2 is the sliding bubbles diameter dynamic changing process figure using the inventive method to obtain.
Detailed description of the invention
First, for shooting small steam bubble clearly, two high-speed camera instrument all select micro-lens, preferred micro-lens Model is Nikon/200mmF/4D/macro.Under conditions of guaranteeing to photograph bubble growth and slipping, Ying Jinke Less pixel can be selected to reduce the amount of calculation of later stage process, for Phantom V9.0 high-speed camera instrument, preferably shoot Pixel is 576 × 288 (high × wide).
Afterwards, by regulation 1# high-speed camera instrument 2 micro-lens aperture, making shooting image become clear, secondly, regulation 1# is at a high speed Video camera 2 micro-lens focal length, makes 1# high-speed camera instrument 2 focus on bubble nucleation point 1, now, for guaranteeing to photograph steam bubble Growth and the whole process of sliding, highly make bubble nucleation point 1 be in the middle part of shooting window by regulation 1# high-speed camera instrument 2 and close Right position is put, and shooting window top is S to bottom distance, and nucleating point 1 is s away from bottom position, then s/S span is 0.3-0.5, It is 0.4 as preferred value.Hereafter, 2# high-speed camera instrument 3 micro-lens focal length is arranged to and 1# high-speed camera instrument 2 one Cause, make the nucleating point 1 position in shooting window with 1# high-speed camera instrument 2 the most partially by regulation 2# high-speed camera instrument 3 position Difference is less than ± 10%, and can photograph picture rich in detail.For obtaining smooth steam bubble dynamic image, two high-speed camera instrument frame numbers Need consistent, and scope should be 500-3000 frame, is 1000 frames as preferred value.
Hereafter, keep high-speed camera instrument focal length constant, respectively with 1# high-speed camera instrument 2 and 2# high-speed camera instrument 3 shooting band The scale of scale, obtains scale map picture clearly.The valid pixel of shooting picture can be changed into actual physical size by this step, Its concrete grammar is, measures the valid pixel at 1mm interval in scale, and result is d, in shooting picture between 2 effectively Pixel is d1, then point-to-point transmission actual physical a size of d in shooting picture1/d mm。
The relation input picture of two the high-speed camera instrument valid pixels obtained with actual physical size is gathered and processes System, the benchmark measured as steam bubble dia.
Afterwards, 1# high-speed camera instrument 2 and 2# high-speed camera instrument 3 is utilized to shoot background picture when not having a steam bubble respectively It is input in image variants system 4, then uses 1# high-speed camera instrument 2 and 2# high-speed camera instrument 3 simultaneously to steam bubble nucleation The sliding bubbles that point 1 position produces shoots, and the view data collected is real-time transmitted to image variants In system 4.
Afterwards, the image that 1# high-speed camera instrument 2 and 2# high-speed camera instrument 3 is shot by image variants system 4 is utilized Carrying out integrated treatment, its concrete grammar is:
First steam bubble in the synchronization image of 1# high-speed camera instrument 2 and 2# high-speed camera instrument 3 shooting is peeled off, its Concrete grammar is to be subtracted each other with background picture corresponding element, after the difference of gained is for peeling off by the picture having steam bubble that experiment shooting obtains Steam bubble image.
Generally, the steam bubble image outline after stripping is the most coarse, and can there are some random noises, as preferably, to stripping Steam bubble image after from is filtered, and preferred filtering method is medium filtering, can make steam bubble edge more by medium filtering Smooth, and eliminate the noise in image.
As preferably, filtered image is carried out binary conversion treatment, it is thus achieved that black white image.Processed at image binaryzation Cheng Zhong, need to determine suitable threshold value f, and when gray value is more than threshold value f, picture point is white;During less than threshold value f, picture point is Black.
For steam bubble image, preferably determine that the method for threshold value f is maximum variance between clusters, its ideal calculated, with Time to calculate speed fast, meet the requirement of kinetic measurement.Maximum variance between clusters principle is: by the grey level histogram of image a certain Being divided into two groups at threshold value, when the variance maximum of this gray value of two groups, cut-point is threshold value.For typical grayscale figure Picture, it is assumed that f is threshold value of first and second group, the most now first group of image is counted and accounted for the total pixel ratio of image is x0, averagely ash Degree is y0, second group of image is counted and accounted for the total pixel ratio of image is x1, average gray is y1, then image grand mean gray scale is:
Y=x0×y0+x1×y1 (1)
The most now variance yields z expression formula is as follows:
Z=x0×(y0-y)2+x1×(y1-y)2 (2)
Threshold value f is enumerated (intensity value ranges is 0~255), and making f maximum for variance yields z is the best cutting point.
As an example, the shooting of synchronization 1# high-speed camera instrument 2 obtains the original image and two in bubble nucleation point front Image after value process, the shooting of 2# high-speed camera instrument 3 obtains original image and the binary conversion treatment thereof of bubble nucleation point side After image.
Afterwards, utilize image variants system 4 that the image after processing is carried out post processing, straight to obtain steam bubble equivalence Footpath real-time change, specific implementation method is as follows:
Image after the original image binary conversion treatment of bubble nucleation point front is laterally progressively scanned, it is thus achieved that every a line Leftmost side black image point and distance D of rightmost side black image pointx1,Dx2,Dx3...Dxn, taking wherein maximum is Dx.Afterwards, Image after the original image binary conversion treatment of bubble nucleation point front is longitudinally scanned by column, it is thus achieved that every string top side is black Chromatic graph picture point and distance D of lower side black image pointy1,Dy2,Dy3...Dyn, taking wherein maximum is Dy
Image after the original image binary conversion treatment of bubble nucleation point side is laterally progressively scanned, it is thus achieved that every a line Leftmost side black image point and distance D of rightmost side black image pointz1,Dz2,Dz3...Dzn, wherein maximum is set to Dz
The most now steam bubble equivalent diameter DbFor:
D b = D x D y D z 3 - - - ( 3 )
Use said method, the image of 1# high-speed camera instrument 2 and 2# high-speed camera instrument 3 shooting is processed, i.e. in real time Steam bubble dia in steam bubble slipping can be obtained dynamically change.
Use said method, it is thus achieved that typical sliding steam bubble dia dynamic changing process as in figure 2 it is shown, in figure the time be 0 Moment be the moment that steam bubble just starts sliding.

Claims (8)

1. a sliding bubbles diameter Three-Dimensional Dynamic measuring method, it is characterised in that: comprise the steps:
Step one, utilize two vertically disposed high-speed camera instrument shooting scale with a scale, obtain scale map picture clearly, will The valid pixel of scale map picture is changed into actual physical size;
Step 2, two the high-speed camera instrument utilizing step one to arrange shoot background picture when not having steam bubble;
Step 3, utilize two high-speed camera instrument that step one arranges respectively from front and side simultaneously to bubble nucleation point position The sliding bubbles produced shoots;
Step 4, sliding bubbles image is peeled off;
Step 5, to peel off after steam bubble image be filtered;
Step 6, filtered image is carried out binary conversion treatment, it is thus achieved that black white image;
The equivalent diameter of step 7, in real time calculating steam bubble.
A kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method the most according to claim 1, it is characterised in that: two high speeds Video camera all selects micro-lens, and the height of the high-speed camera instrument of shooting bubble nucleation point direct picture is it is ensured that make steam bubble nucleation Point is in correct position in the middle part of shooting window, and two high-speed camera instrument frame numbers are consistent, and scope is 500-3000 frame.
A kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method the most according to claim 1, it is characterised in that: by scale map The valid pixel of picture is changed into the method for actual physical size: measures the valid pixel at 1mm interval in scale, obtains Measurement result d, if between 2, valid pixel is d in scale map picture1, then actual physical size between these 2 in scale map picture For d1/ d, unit is mm.
A kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method the most according to claim 1, it is characterised in that: to sliding vapour The method that bubble image carries out peeling off is: the Background that sliding bubbles picture step 3 shooting obtained obtains with step 2 shooting The corresponding element of sheet subtracts each other, and the difference of gained is the steam bubble image after stripping.
A kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method the most according to claim 1, it is characterised in that: described filtering Method is medium filtering.
A kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method the most according to claim 1, it is characterised in that: after filtering Image carry out the method for binary conversion treatment and be: first determine threshold value f, then judge the gray value of picture point whether more than threshold value f, If more than, determine that picture point is white;Otherwise, it is determined that picture point is black.
A kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method the most according to claim 6, it is characterised in that: use maximum Ostu method determines threshold value f, comprises the steps:
Assuming that threshold value that f is first and second group, the most now first group of image is counted and accounted for the total pixel ratio of image is x0, averagely ash Degree is y0, second group of image is counted and accounted for the total pixel ratio of image is x1, average gray is y1, then image grand mean gray scale is:
Y=x0×y0+x1×y1
The most now variance yields z expression formula is as follows:
Z=x0×(y0-y)2+x1×(y1-y)2
Threshold value f is enumerated, makes f maximum for variance yields z i.e. be defined as threshold value.
A kind of sliding bubbles diameter Three-Dimensional Dynamic measuring method the most according to claim 1, it is characterised in that: calculate in real time The method of steam bubble equivalent diameter is:
(1) direct picture real-time to the bubble nucleation point after binary conversion treatment laterally progressively scans, it is thus achieved that every a line is the most left Side black image point and distance D of rightmost side black image pointx1,Dx2,Dx3...Dxn, take wherein maximum and be designated as Dx
(2) direct picture real-time to the bubble nucleation point after binary conversion treatment longitudinally scans by column, it is thus achieved that every string is gone up most Side black image point and distance D of lower side black image pointy1,Dy2,Dy3...Dyn, take wherein maximum and be designated as Dy
(3) side image real-time to the bubble nucleation point after binary conversion treatment laterally progressively scans, it is thus achieved that every a line is the most left Side black image point and distance D of rightmost side black image pointz1,Dz2,Dz3...Dzn, it is designated as D to wherein maximumz
(4) steam bubble equivalent diameter D it is calculated as followsb:
D b = D x D y D z 3 .
CN201611029639.8A 2016-11-15 2016-11-15 Method for three-dimensionally and dynamically measuring diameter of slipping bubble Pending CN106323184A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611029639.8A CN106323184A (en) 2016-11-15 2016-11-15 Method for three-dimensionally and dynamically measuring diameter of slipping bubble

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611029639.8A CN106323184A (en) 2016-11-15 2016-11-15 Method for three-dimensionally and dynamically measuring diameter of slipping bubble

Publications (1)

Publication Number Publication Date
CN106323184A true CN106323184A (en) 2017-01-11

Family

ID=57817284

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611029639.8A Pending CN106323184A (en) 2016-11-15 2016-11-15 Method for three-dimensionally and dynamically measuring diameter of slipping bubble

Country Status (1)

Country Link
CN (1) CN106323184A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106989672A (en) * 2017-04-17 2017-07-28 天津大学 A kind of workpiece measuring based on machine vision
CN110544241A (en) * 2019-08-22 2019-12-06 云南电网有限责任公司 matching system for power cable metal cross section trace feature recognition degree

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884645A (en) * 2012-12-24 2014-06-25 中国制浆造纸研究院 Method for testing gas content in paint
CN104457703A (en) * 2014-11-28 2015-03-25 浙江大学 Multiple-vision-based System and method for measuring parameters of gas-liquid two-phase flow in small channel based on
CN106023153A (en) * 2016-05-09 2016-10-12 山东省科学院海洋仪器仪表研究所 Method for measuring bubbles in water body
CN106530255A (en) * 2016-11-15 2017-03-22 中国核动力研究设计院 Bubble group data measurement method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103884645A (en) * 2012-12-24 2014-06-25 中国制浆造纸研究院 Method for testing gas content in paint
CN104457703A (en) * 2014-11-28 2015-03-25 浙江大学 Multiple-vision-based System and method for measuring parameters of gas-liquid two-phase flow in small channel based on
CN106023153A (en) * 2016-05-09 2016-10-12 山东省科学院海洋仪器仪表研究所 Method for measuring bubbles in water body
CN106530255A (en) * 2016-11-15 2017-03-22 中国核动力研究设计院 Bubble group data measurement method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
徐建军等: "矩形窄缝通道内滑移汽泡生长及界面形态可视化研究", 《核动力工程》 *
李少丹: "摇摆运动下窄通道内过冷沸腾汽泡行为研究", 《中国优秀硕士学位论文全文数据库》 *
谢添舟等: "矩形窄缝通道内汽泡脱离直径可视化研究", 《核动力工程》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106989672A (en) * 2017-04-17 2017-07-28 天津大学 A kind of workpiece measuring based on machine vision
CN110544241A (en) * 2019-08-22 2019-12-06 云南电网有限责任公司 matching system for power cable metal cross section trace feature recognition degree

Similar Documents

Publication Publication Date Title
CN108387587A (en) Defect inspection method and defect detection equipment
Dabiri et al. Experimental investigation of the vorticity generation within a spilling water wave
CN106530255A (en) Bubble group data measurement method
CN103558229B (en) A kind of MURA vision automatic testing method of TFT-LCD processing procedure and device
CN106372584B (en) A kind of video image mosaic detection method
CN103759676A (en) Non-contact type workpiece surface roughness detecting method
CN102680480A (en) Intelligent detecting method for cracks of concrete structures
CN103793918A (en) Image definition detecting method and device
CN103808265B (en) Method, device and system for measuring oilseed rape laminae and forms of sclerotium scabs synchronously
CN103984979A (en) Lens-diffraction-imaging-free automatic algae detection and counting device and method
CN106323184A (en) Method for three-dimensionally and dynamically measuring diameter of slipping bubble
CN204039886U (en) A kind of pavement damage crack detection system based on multiple stage camera stereoscopic shooting
CN104198355A (en) Automatic detection method for red cells in feces
CN101651845B (en) Method for testing definition of moving images of display devices
CN101615241A (en) A kind of screening technique of certificate photograph
CN102999782B (en) Based on slice counter and the method for counting of image procossing
CN102945547B (en) Cold-rolled steel plate surface image illumination homogenizing method
CN103150725A (en) SUSAN edge detection method and system based on non-local mean values
CN109028230A (en) Have the stove and oil smoke concentration detection method of gesture control vision-based detection function
CN102564924B (en) Automatic scanning method of single-frame image of blood cell
CN105608674B (en) A kind of image enchancing method based on image registration, interpolation and denoising
CN102547363A (en) No-reference image quality evaluating method on basis of contourlet transform domain image energy features
CN104346791B (en) A kind of method of particle imaging defect in correction flow cell
CN101605269B (en) Method and device for tracking dense depth images
Gharahjeh et al. Application of video imagery techniques for low cost measurement of water surface velocity in open channels

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20170111

RJ01 Rejection of invention patent application after publication