CN100434875C - Ultrasonic flowmeter and its theory and technique - Google Patents

Ultrasonic flowmeter and its theory and technique Download PDF

Info

Publication number
CN100434875C
CN100434875C CNB200510088694XA CN200510088694A CN100434875C CN 100434875 C CN100434875 C CN 100434875C CN B200510088694X A CNB200510088694X A CN B200510088694XA CN 200510088694 A CN200510088694 A CN 200510088694A CN 100434875 C CN100434875 C CN 100434875C
Authority
CN
China
Prior art keywords
formula
prime
centerdot
tau
fluid
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
Application number
CNB200510088694XA
Other languages
Chinese (zh)
Other versions
CN1908592A (en
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CNB200510088694XA priority Critical patent/CN100434875C/en
Publication of CN1908592A publication Critical patent/CN1908592A/en
Application granted granted Critical
Publication of CN100434875C publication Critical patent/CN100434875C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Volume Flow (AREA)

Abstract

The related ultrasonic flowmeter comprises: a sampling part, a computation part, and a keyboard display part, wherein sampling the oscillating number of oscillator 19 clockwise and counter-clockwise to obtain the sampling data N1' and N2'; and computing the flow velocity and flow rate by a SCM 24 based on principles and the 'extended transmission length technology'.

Description

A kind of ultrasonic flow meter is measured the method for fluid flow
Technical field
The present invention relates to a kind of can measuring channel in the ultrasonic flow meter of flow velocity, flow and ultrasound wave thereof velocity of propagation in its fluid of various fluids (potpourri of liquid, gas, mixtures of different liquids, gas with various) measure the method for fluid flow.Especially can significantly improve the measuring accuracy of ultrasonic flow meter, significantly enlarge measurement lower limit, significantly be fit to various calibers, and can record this fluid when static in current fluid, ultrasound wave is measured the method for fluid flow in the ultrasonic flow meter of velocity of propagation wherein.The invention still further relates to ultrasound wave heat meter and flowmeter standard scale.
Background technology
At present, the measuring principle of known ultrasonic flow meter has multiple different form, wherein mainly contains following two kinds, and the principle of their institute's foundations is as follows:
One, time difference method
According to Fig. 5:
The 1st, fluid downbeam ultrasonic emitting oscillator T 1
The 2nd, fluid downbeam ultrasound wave pick-up dipole R 1
The 3rd, fluid countercurrent current direction ultrasound wave pick-up dipole R 2
The 4th, fluid countercurrent current direction ultrasonic emitting oscillator T 2
L is T 1And R 1Between distance, also be T 2With R 2Between distance
C is the velocity of propagation of ultrasound wave in stationary fluid
U is a flow rate of fluid
" → u " represents the flow direction of fluid is the arrow direction
According to Fig. 5, when the ultrasonic propagation direction is consistent with fluid flow direction, i.e. fluid downbeam ultrasonic emitting oscillator T 1Emission ultrasound wave, and fluid downbeam ultrasound wave pick-up dipole R 1When receiving ultrasound wave, its velocity of propagation is c+u; And the ultrasonic propagation direction is when opposite with fluid flow direction, i.e. fluid countercurrent current direction ultrasonic emitting oscillator T 2Emission ultrasound wave, and fluid countercurrent current direction ultrasound wave pick-up dipole R 2When receiving ultrasound wave, its velocity of propagation is c-u.Hence one can see that, establishes t 1Be ultrasound wave in fluid downbeam by T 1Propagate into R 1The used time:
Then
t 1 = L c + u - - - ( 1 )
If t 2Be ultrasound wave in fluid countercurrent direction by T 2Propagate into R 2The used time then
t 2 = L c - u - - - ( 2 )
Get by (1) formula and (2) formula
2 u = L ( 1 t 1 - 1 t 2 ) - - - ( 3 )
Therefore, ultrasonic propagation time difference Δ t=t under the concurrent-countercurrent situation 2-t 1, then (3) formula can be expressed as:
u = Δt · L 2 t 1 · ( Δt + t 1 ) - - - ( 4 )
By (4) formula as can be known, rate of flow of fluid u with
Figure C20051008869400093
Be directly proportional scale-up factor
Figure C20051008869400094
Be constant, be not subjected to influence of temperature change, only need ultrasound wave following current travel-time t 1And contrary, the following current of ultrasound wave is carried out computing and can be tried to achieve rate of flow of fluid u in L propagation mistiming Δ t substitution.
Above time difference method is write with reference to the 305th page, 310 pages, 311 pages of " flow measurement technology and instrument " books of China Machine Press's publication.This book is engraved into the printing of printing company limited by Beijing, the distribution of Xinhua Bookstore Beijing sale room, in June, the 2002 first published first impression.
Chief editor: beam state big Cai Wuchang
Coder: Zheng Jian Ying Cheng healthy and full of vigor peak
Two, sound round-robin method
Fig. 6 is the block scheme of an example of the ultrasonic flow meter of expression prior art.Ultrasonic flow meter shown in Figure 6 for example is disclosed in the non-patent document 1.As shown in Figure 6, ultrasonic oscillator 1 and 2 is configured to clamp the state of the stream 14 that flows through fluid. Ultrasonic oscillator 1 and 2 has played the function of transmitter and receptacle respectively.In other words, when ultrasonic oscillator 1 used as transmitter, ultrasonic oscillator 2 used as receiver; When ultrasonic oscillator 2 used as transmitter, ultrasonic oscillator 1 used as receiver.As shown in Figure 6, between the ultrasonic oscillator 1 and 2 flow direction of formed ultrasonic wave propagation path and fluid to be the θ angle lapping oblique.
When make ultrasound wave from ultrasonic oscillator 1 when ultrasonic oscillator 2 is propagated, because the relative fluid of ultrasound wave flows along advancing along direction, its speed just accelerates, on the contrary, when ultrasound wave from ultrasonic oscillator 2 when ultrasonic oscillator 1 is propagated, advance along contrary direction because the relative fluid of ultrasound wave flows, its speed is just slack-off.Therefore, poor by the ultrasonic propagation time of ultrasonic oscillator 1 to the ultrasonic propagation time of ultrasonic oscillator 2 and ultrasonic oscillator 2 to ultrasonic oscillator 1 can be obtained the speed of fluid.In addition, can also be by the sectional area of stream 14 and the long-pending flow of obtaining of flow velocity.
According to above-mentioned principle,, the measuring method by the sound round-robin method is described as the concrete grammar of asking fluid flow.
As shown in Figure 6, ultrasonic flow meter comprises sending part 3 and receiving portion 6, and ultrasonic oscillator 1 selectively is connected with a side in sending part 3 or the acceptance division 6 by switching part 10.At this moment, the sending part 3 that is not connected with ultrasonic oscillator 1 of ultrasonic oscillator 2 or the side in the acceptance division 6 connect.
When sending part 3 was connected with ultrasonic oscillator 1, sending part 3 drove ultrasonic oscillator 1, and the ultrasound wave of generation crosses mobile fluid and arrives ultrasonic oscillator 2.Ultrasound wave by ultrasonic oscillator 2 is received is converted into electric signal, and received signal is amplified by acceptance division 6.In zero cross detection portion 7, detect received signal and reach specified level first zero cross point afterwards, generate the zero cross detection signal.So-called zero cross point is meant, the amplitude of received signal from just to negative, or from negative to the point that is just changing.This zero cross point, arrive the moment of ultrasonic oscillator 2 as ultrasound wave.According to the zero cross detection signal, generate trigger pip in the moment of delay stipulated time, be input to sending part 3.Being called time delay from generation zero cross detection signal to the time that generates the trigger pip.
Sending part 3 drives ultrasonic oscillator 1 according to trigger pip, produces next ultrasound wave.The annular of this hyperacoustic transmission-reception-amplification, delay-transmission circulated repeatedly call sound circulation.Cycle index is called period.
In timing portion 9, the time that instrumentation stipulated number, repetitive cycling need, and the instrumentation result is sent to flow rate calculation portion 11.Then, switch switching part 10, ultrasonic oscillator 2 is used as transmitter, when ultrasonic oscillator 1 is used as receiver, carry out same measurement.
Deduct the value that delay time and sound period multiply each other from the time of measuring by said method, the value divided by gained after the sound period is ultrasonic propagation time again.If the travel-time t with ultrasonic oscillator 1 during as transmitter side 1, the travel-time t when establishing ultrasonic oscillator 2 as transmitter side 2
In addition, as shown in Figure 6, the distance of establishing between ultrasonic oscillator 1 and the ultrasonic oscillator 2 is L, and establishes flow rate of fluid and hyperacoustic velocity of sound is respectively u and c.
At this moment, t 1And t 2Formula below available is represented.
t 1 = L c + u cos θ t 2 = L c - u cos θ - - - ( 5 )
Can represent the following formula of flow velocity by following formula,
u = L 2 cos θ ( 1 t 1 - 1 t 2 ) - - - ( 6 )
If can obtain rate of flow of fluid u, just can be by the sectional area of stream 14 and the long-pending flow Q that obtains of flow velocity u.
(non-patent document)
NEC tester industry meeting specification, JEMI5032 " hyperacoustic hydrometry " NEC instrumentation industry meeting 1987.The sound round-robin method is taken passages in following " application for a patent for invention prospectus ".
[19] State Intellectual Property Office of the People's Republic of China
[12] application for a patent for invention prospectus
[21] application number 03123080.6
[51]INC CI 7
GOIF 1/66
[43] open day on November 12nd, 2003
[11] publication number CN1455230A
[22] April 30 2003 applying date
[21] application number 03123080.6
[71] applicant Panasonic Electric Equipment Industrial Co.,Ltd
Osaka, JAPAN mansion, address
[72] firm this refined man of virtue and ability of man of virtue and ability's bridge in the inventor China fir
[74] Zhongke Patent ﹠ Trademark Agency Co., Ltd of patent agency
Procurator Wang Huimin
Main two kinds of measuring principles of known ultrasonic flow meter in the above background technology: though time difference method harmony round-robin method has many superior parts, but these known ultrasonic flow meters are used for measuring little caliber rate of flow of fluid hour, because the precision of its flowmeter is subjected to instrument self to the time measurement error, the fluid caliber, the indeterminable constraint of rate of flow of fluid and physical quantity, promptly because the precision of flowmeter can not be higher than the precision of this flowmeter to time measurement, and caliber is more little, flow velocity is more little, the longshore current body is suitable respectively by caliber one side for ultrasound wave, it is more little that countercurrent direction propagates into mistiming of opposite side, little when this flowmeter self is unreachable to the time measurement precision when this mistiming, it is measured and just to have produced very big error, so that energy measurement not.
Summary of the invention
Be used to measure tubule footpath inner fluid speed hour in order to overcome known existing ultrasonic flow meter and principle thereof and technology, produce big error, the technical matters that its precision is greatly reduced the invention provides a kind of ultrasonic flow meter and principle thereof and technology.This ultrasonic flow meter and principle thereof and technology, can not only be applicable to the various occasions that present known ultrasonic flow meter is suitable for, and can significantly reduce tubule footpath inner fluid speed relative error hour, significantly improve the measuring accuracy of ultrasonic flow meter, significantly enlarge the measurement lower limit (can survey flow velocity littler in the little caliber) of ultrasonic flow meter, significantly be fit to the measurement of all size caliber inner fluid speed, flow.
The technical scheme that the present invention solves its technical matters is: the new principle that proposes a kind of ultrasonic flow meter measuring tube inner fluid speed, adopted a kind of new technology, having overcome known ultrasonic flow meter precision with this principle and technology retrained by fluid caliber and rate of flow of fluid, solved when rate of flow of fluid is very little, ultrasound wave respectively by caliber one side longshore current rate of flow of fluid along contrary direction propagate into mistiming of opposite side little to existing ultrasonic flow meter self to the time measurement precision when unreachable, to its immeasurablel problem.And with this principle and technological guidance and implemented the electronic circuit of ultrasonic flow meter.Discuss respectively from these three aspects below.
The new principle of ultrasonic flow meter measuring tube inner fluid speed:
Present principles is that 1. time difference methods in the background technology are furtherd investigate.Ultrasonic flow meter of the present invention and principle thereof and technology implementation example 1 block diagram:
If " → u " is the arrow direction for fluid flow direction, u is the mean flow rate of fluid on diameter D, c be ultrasound wave when this fluid is in static state in velocity of propagation wherein; The angle of its direction and rate of flow of fluid when θ propagates for the ultrasound wave following current, c 1Be the velocity of propagation of ultrasound wave in fluid is flowed through pipe 35 media.
Ultrasound wave following current route of transmission and circuit working are promptly as follows along cyclic process:
Sending part 1 drive ultrasonic oscillator 28 emission ultrasound waves → through fluid flow through pipe 35 tube wall upside m points enter fluid flow through pipe 35 → along mo penetrate fluid flow through pipe 35 upside tube walls incide the o point → o point through refraction along oo ' penetrate fluid incide o ' point → o ' through refraction along o ' m ' penetrate fluid flow through manage 35 downside tube walls in m ' by ultrasonic oscillator 29 receptions → through the peak detector 3 of acceptance division 2 amplifications → just → or 7 → with a 11 → sending part 1, form ultrasound wave along circulation, establishing its ultrasound wave is f along cycle frequency s., the cycle is T S
The promptly contrary cyclic process of ultrasound wave adverse current route of transmission and circuit working is as follows:
Sending part 4 drive ultrasonic oscillators 30 emission ultrasound waves → through fluid flow through pipe 35 tube wall upside v points enter fluid flow through pipe 35 → along vu penetrate fluid flow through pipe 35 upside tube walls incide the u point → u point through refraction along uu ' penetrate fluid incide u ' point → u ' through refraction along u ' v ' penetrate fluid flow through manage 35 downside tube walls in v ' by ultrasonic oscillator 31 receptions → through the peak detector 6 of acceptance division 5 amplifications → just → or 9 → with a 12 → sending part 4, form the contrary circulation of ultrasound wave, establishing the contrary cycle frequency of its ultrasound wave is f N., the cycle is T N
If in the ultrasound wave following current route of transmission:
t 1It is ultrasound wave required time of propagation distance oo ' in fluid;
Figure C20051008869400121
Be that ultrasound wave is flowed through the pipe 35 tube wall upsides propagation mo used time of distance at fluid, flow through because of fluid and manage tube wall symmetry about in the of 35, so
Figure C20051008869400122
Also be that ultrasound wave is flowed through the 35 used times of tube wall downside propagation distance o ' m ' of pipe at fluid.
To be ultrasound wave pick-up dipole 29 receive the circuit delay time that ultrasound wave is ordered through acceptance division 2, positive peak detector 3 or door 7, with door 11, sending part 1, ultrasonic emitting oscillator 28 to m at m ' to τ '.
t 1' be ultrasound wave from m ' through τ ' institute through circuit, again through τ institute through approach and t 1Through time of approach and.
If ultrasound wave following current route of transmission and adverse current route of transmission are to axle xx ' rotational symmetry.
If in the ultrasound wave adverse current route of transmission:
t 2It is ultrasound wave required time of propagation distance uu ' (uu '=oo ') in fluid;
Figure C20051008869400123
Being ultrasound wave flows through the 35 used times of tube wall upside propagation distance vu (vu=mo) of pipe at fluid, and also to be ultrasound wave certainly at fluid flow through manages the used time of 35 tube wall downside propagation distance u ' v ' (u ' v '=o ' m ').
" to be ultrasound wave pick-up dipole 31 receive the circuit delay time that ultrasound wave is ordered through acceptance division 5, positive peak detector 6 or door 9, with door 12, sending part 4, ultrasonic emitting oscillator 30 to v at v ' to τ.
t 2' be ultrasound wave from v ' through τ " institute is through circuit, again through τ institute through approach and t 2Through time of approach and.
According to above-mentioned set, should there be following formula to set up
t 1 ′ = t 1 + τ + τ ′ t 2 ′ = t 2 + τ + τ ′ ′ - - - ( 7 )
(7) in the formula
t 1 = o o ′ c + u cos θ t 2 = u u ′ c - u cos θ = o o ′ c - u cos θ - - - ( 8 )
By (8) formula
c + u cos θ = o o ′ t 1 - - - ( 9 )
c - u cos θ = o o ′ t 2 - - - ( 10 )
Because
c=f·λ (11)
(11) in the formula:
F is the supersonic oscillations frequencies, and the oscillation frequency value that to establish this f value be exactly standard time-base generator 19.The ultrasound wave wavelength that λ is a ultrasound wave when above-mentioned f frequency, this frequency is far longer than f s, f N
We remove rate of flow of fluid u with λ more now,
u λ = n - - - ( 12 )
(12) n is the multiple of u to λ in the formula.
The levoform of (11) formula and (12) formula substitution (9) formula and (10) formula is had
c + u cos θ = fλ + nλ cos θ
= f ( 1 + n cos θ f ) λ - - - ( 13 )
c - u cos θ = fλ - nλ cos θ
= f ( 1 - n cos θ f ) λ - - - ( 14 )
If (13) in the formula in the right formula
f ( 1 + n cos θ f ) λ = f λ 1 - - - ( 15 )
If (14) in the formula in the right formula
f ( 1 - n cos θ f ) λ = f λ 2 - - - ( 16 )
We define herein
Figure C20051008869400138
Be the influence coefficient of rate of flow of fluid to the ultrasound wave wavelength
With (15) formula substitution (13) formula,
c+ucosθ=fλ 1 (17)
λ 1Be the wavelength of ultrasound wave when oo ' direction is propagated
With (16) formula substitution (14) formula,
c-ucosθ=fλ 2 (18)
λ 2Be the wavelength of ultrasound wave when uu ' direction is propagated
(17) formula subtracts (18) formula,
2ucosθ=f(λ 12) (19)
With (17) formula, (18) formula substitution (8) formula,
t 1 f = o o ′ λ 1 t 2 f = o o ′ λ 2 - - - ( 20 )
Establish again
t 1 f = n 1 t 2 f = n 2 - - - ( 21 )
n 1Be ultrasound wave during along oo ' propagation, the ultrasound wave wave number on oo ' distance
n 2Be ultrasound wave during along uu ' propagation, the ultrasound wave wave number on uu ' distance
With (21) formula substitution (20) formula,
λ 1 = o o ′ n 1 - - - ( 22 )
λ 2 = o o ′ n 2 - - - ( 23 )
With (22) formula and (23) formula substitution (19) formula,
2 u cos θ = f ( o o ′ n 1 - o o ′ n 2 )
Promptly
u = ( 1 n 1 - 1 n 2 ) · o o ′ 2 cos θ · f - - - ( 24 )
(24) formula is the new principle of ultrasonic flow meter measuring tube inner fluid speed of the present invention.
The new technology that adopts:
New technology of the present invention is: when prolonging ultrasonic flow meter measurement rate of flow of fluid, the propagation distance of ultrasound wave in fluid is called for short " prolong and pass apart from technology " herein.
The following enforcement of this technology:
The molecule and the denominator of (24) formula equation the right formula be multiply by a positive integer N together,
Promptly
u = ( 1 n 1 - 1 n 2 ) · N N · o o ′ 2 cos θ · f - - - ( 25 )
If
N 2=n 2·N (26)
N 2Be uu ' expansion N positive integer doubly after, ultrasound wave after enlarging apart from Nuu ' (Nuu '=Noo ') when propagating, the ultrasound wave wave number on distance Nuu '.
N 1=n 1·N (27)
N 1Be oo ' expansion N positive integer doubly after, ultrasound wave after enlarging apart from Noo ' propagation the time, the ultrasound wave wave number on distance Noo '.
With (26) formula and (27) formula substitution (25) formula,
u = ( 1 N 1 - 1 N 2 ) · N · oo ′ 2 cos θ · f - - - ( 28 )
Can find out N by (28) formula 1, N 2Be respectively oo ' and uu ', after promptly oo ' had enlarged N times, in the concurrent-countercurrent ultrasonic propagation required time, the number of oscillation of standard time-base generator 19 was because of oo ', uu ' all enlarge N doubly, so be called " prolong and pass apart from technology ".
(28) formula be the present invention after new technology is applied to new principle (24) formula, the new formula of the measurement that draws stream tube fluid flow velocity.In addition, new formula (28) formula of tube fluid flow velocity is flowed in the measurement that new principle (24) formula of ultrasonic flow meter measuring tube inner fluid speed of the present invention, new technology are applied to draw after new principle (24) formula, also can directly be released by this instructions (3) formula.
When (3) formula was applied to Fig. 1, the equation left side should be written as 2ucos θ, and (3) formula molecule, denominator get with multiply by f
u = ( 1 t 1 - 1 t 2 ) · f f · L 2 cos θ - - - ( 29 )
(29) in the formula, establish
t 2f=n 2 (30)
t 1f=n 1 (31)
(30) formula and (31) formula are actual is exactly (21) formula, with (30) formula and (31) formula substitution (29) formula, that is:
u = ( 1 n 1 - 1 n 2 ) · L 2 cos θ · f - - - ( 32 )
(32) formula is new principle (24) formula of ultrasonic flow meter measuring tube inner fluid speed of the present invention.
(32) formula equation the right molecule and denominator be multiply by positive integer N together, and with (26) formula and (27) formula substitution,
u = ( 1 n 1 - 1 n 2 ) · N N · L 2 cos θ · f
= ( 1 N 1 - 1 N 2 ) · N · L 2 cos θ · f - - - ( 33 )
(33) formula is new formula (28) formula of the measurement stream tube fluid flow velocity that new technology of the present invention is applied to draw after new principle (24) formula.
Now, we use the new principle and the new technology of this ultrasonic flow again, try to achieve when detected fluid is in static state, and ultrasound wave is in velocity of propagation c wherein.
With (17) formula and the addition respectively of (18) formula both members
2c=f(λ 12) (34)
Get by (22) formula n 1 = o o ′ λ 1 (23) formula gets n 2 = o o ′ λ 2 Substitution (34) Shi Kede
2 c = f ( o o ′ n 1 + o o ′ n 2 ) - - - ( 35 )
Get by (35) formula
c = ( 1 n 1 + 1 n 2 ) · o o ′ 2 · f - - - ( 36 )
Be called for short " prolong and pass " with the present invention again, (36) formula the right formula molecule denominator be multiply by positive integer N together apart from technology,
c = ( 1 n 1 + 1 n 2 ) · N N · o o ′ 2 · f
= ( 1 N 1 + 1 N 2 ) · N · o o ′ 2 · f - - - ( 37 )
Because among Fig. 1
o o ′ = D sin θ - - - ( 38 )
With (38) formula difference substitution (24) formula, (28) formula, (36) formula and (37) formula, draw following various respectively
u = ( 1 n 1 - 1 n 2 ) · D / sin θ 2 cos θ · f - - - ( 39 )
u = ( 1 N 1 - 1 N 2 ) · N · D / sin θ 2 cos θ · f - - - ( 40 )
c = ( 1 n 1 + 1 n 2 ) · D / sin θ 2 · f - - - ( 41 )
c = ( 1 N 1 + 1 N 2 ) · N · D / sin θ 2 · f - - - ( 42 )
In (39) formula and (41) formula, by (8) formula and (21) formula as can be known, n 1Be that the ultrasound wave following current is propagated oo ' apart from required time t 1The number of oscillation of interior standard time-base generator 19; n 2Be that the ultrasound wave adverse current is propagated uu ' (uu '=oo ') apart from required time t 2The number of oscillation of interior standard time-base generator 19.
And in (40) formula and (42) formula, by (25) formula, (26) formula and (27) formula as can be known, N 1Be that doubly back required time t of oo ' distance expansion N is propagated in the ultrasound wave following current 1The number of oscillation of the standard time-base generator 19 in the N; N 2Be that the ultrasound wave adverse current is propagated uu ' (uu '=oo ') distance and enlarged doubly back required time t of N 2The number of oscillation of the standard time-base generator 19 in the N.
D, f and θ are known numbers in (39) formula and (41) formula, so as long as energy measurement goes out n 1And n 2With regard to available (39) formula and (41) formula calculate respectively fluid flow through rate of flow of fluid u in the pipe 35 when being in static state with this fluid ultrasound wave and be under the different temperatures at fluid and all can measure in velocity of propagation c (c is measured in measurement during hydrodynamic) wherein.
But with (39) formula fluid measured flow velocity u, usefulness (41) formula survey ultrasound wave when fluid is static when wherein the velocity of propagation c, when having only f frequency in (39) formula and (41) formula more and more higher, could more and more reduce the measuring error of u and c, (39) formula and (41) formula be not owing to use new technology of the present invention, promptly be called for short " prolong and pass ", so it surveys result and " time difference method 1 " in the background technology and " sound round-robin method 2 " (n of u and c apart from technology 2And n 1Round-robin method equally after N time capable of circulation, is obtained its arithmetic mean in unison) be in same level, we are with (39) formula
u = ( 1 n 1 - 1 n 2 ) · D / sin θ 2 cos θ · f
Be called " new principle of ultrasonic measurement tube fluid flow velocity ".
And with (41) formula
c = ( 1 n 1 + 1 n 2 ) · D / sin θ 2 · f
Be called " ultrasound wave when tube fluid is dynamic, measure this fluid when static ultrasound wave in the new principle of velocity of propagation c wherein ".
For reach mention in this instructions " affiliated technical field " content " the present invention relates to a kind of can measuring channel in the flow velocity of various fluids (liquid; gas; mixtures of different liquids; the potpourri of gas with various); flow; and the ultrasonic flow meter of ultrasound wave velocity of propagation in its fluid and principle and technology, especially can significantly improve the measuring accuracy of ultrasonic flow meter, significantly enlarge measurement lower limit, significantly be fit to various calibers, and in the time of recording this fluid in current fluid and be in static state, ultrasound wave is in wherein ultrasonic flow meter and the principle and the technology of velocity of propagation " purpose.The present invention will " prolong pass apart from technology " be applied to (39) formula promptly " new principle of ultrasonic measurement tube fluid flow velocity " [hereinafter to be referred as " new principle (39) formula "] and (41) formula promptly " ultrasound wave when tube fluid is dynamic, measure this fluid when static ultrasound wave in the new principle of velocity of propagation c wherein " (hereinafter to be referred as " survey c new principle (41) formula "), get (40) formula and (42) formula respectively, promptly
u = N 2 - N 1 N 2 · N 1 · N · D / sin θ 2 cos θ · f
c = N 2 + N 1 N 2 · N 1 · N · D 2 sin θ · f
When using velocity of propagation c in this fluid of (40) formula and (42) formula fluid measured flow velocity u and ultrasound wave, can realize purpose of the present invention.(40) formula has contained the present invention's " new principle (39) formula " and new technology of the present invention " prolongation passes apart from technology ".
(42) formula has contained the present invention's " survey c new principle (41) formula " and new technology of the present invention " prolongation passes apart from technology ".
Should particularly point out: new technology of the present invention " prolongs and passes apart from technology " (6) formula that also can be used in the background technology " sound round-robin method 2 ", is about to (6) formula
u = L 2 cos θ ( 1 t 1 - 1 t 2 )
Equation the right molecule, denominator get with multiply by positive integer N
u = L 2 cos θ · N N · ( 1 t 1 - 1 t 2 )
= N · L 2 cos θ · ( 1 N · t 1 - 1 N · t 2 ) - - - ( 43 )
If (43) in the formula
N · t 2 = T 2 N · t 1 = T 1 - - - ( 44 )
(44) in the formula, T 2Be in " sound round-robin method 2 " electronic circuit block diagram 6, with ultrasonic oscillator 2 during as emitting side the emission ultrasound wave to ultrasonic oscillator 1 required time t 2N doubly, T 1Be with ultrasonic oscillator 1 as emitting side time emission ultrasound wave to ultrasonic oscillator 2 required time t 1N doubly.
With (44) formula substitution (43) formula,
u = ( 1 T 1 - 1 T 2 ) · N · L 2 cos θ - - - ( 45 )
Can reach goal of the invention of the present invention with the rate of flow of fluid u in the stream 14 of the fluid in (45) formula survey sheet 6.
(45) formula only is another expression-form of (40) formula in fact.Because if establishing frequency in (40) formula and be its cycle of vibration of the f frequency of punctual base oscillator 19 (Fig. 1 get the bid) is T, with the molecule of (40) formula the right formula, denominator with multiply by period T, then
u = ( 1 N 1 - 1 N 2 ) · T T · N · D / sin θ 2 cos θ · f
= ( 1 T 1 - 1 T 2 ) · N · L 2 cos θ
This formula is exactly (45) formula.So (45) formula also just becomes the content that the present invention is contained naturally.More than " prolong and pass " owing to adopted apart from technology, so its measuring error is far smaller than the error of known above-mentioned ultrasonic flow meter, make the precision of ultrasonic flow meter of the present invention be higher than the precision of known above-mentioned ultrasonic flow meter far away, thereby reached purpose of the present invention.
According to the above, the present invention has finished following four big contents from theory and technology, that is:
1. propose new principle of the present invention (39) formula, survey c new principle (41) formula.
2. new technology of the present invention is proposed---" prolong and pass ", i.e. (40) formula equation the right molecule apart from technology The distance that ultrasound wave is propagated in limited caliber D
Figure C20051008869400182
Enlarge N doubly.Make the N in the formula of the right in (40) formula simultaneously 2And N 1Becoming ultrasound wave respectively propagates in contrary, following current
Figure C20051008869400183
The number of oscillation apart from the standard time-base generator 19 in the required time; NL makes the distance L of ultrasonic propagation enlarge N doubly in the molecule of (45) formula equation the right, makes the T in the formula on the right in (45) formula simultaneously 2, T 1Become ultrasound wave respectively and broadcast the required time of NL distance in contrary, saequential transmission.
3. (40) formula, (42) formula are proposed promptly
u = ( 1 N 1 - 1 N 2 ) · N · D / sin θ 2 cos θ · f
c = ( 1 N 1 + 1 N 2 ) · N · D 2 sin θ · f
4. (45) formula is proposed promptly
u = ( 1 T 1 - 1 T 2 ) · N · L 2 cos θ
Now, come in the calculating chart 1 the fluid pipe 35 of flowing through, the fluid flow Q of elapsed time t tIn actual measurement, fluid pipe 35 the flow velocity u that flows through is that t changes in time, and we will measure Q tTime t be divided into a lot of junior unit Λ t i, can be written as:
t = Σ i = 1 n Δ t i , i = 1,2,3 · · · · · · , n - - - ( 46 )
Now establish fluid and flow through that pipe 35 xsect inner headed faces are long-pending to be S, then
S = π D 2 4 - - - ( 47 )
If u iBe at Λ t iThe mean flow rate of the fluid in the time on the S area, then when (46) formula i value is big more, Λ t iValue is just more little.As Λ t iLittle to a certain degree the time, at Λ t iFluid in time flows can regard steady flow as.Λ t like this iMean flow rate u in time on the S area iJust equal this steady flow flow rate of fluid, so can try to achieve Λ t iIn time, the flow Q on the represented S area of (47) formula i
Q i = u ‾ i · π D 2 4 · Δ t i - - - ( 48 )
If Q tBe by the flow in the t time on the represented S area of (47) formula, then in (46) formula t time
Q t = Σ i = 1 n Q i = Σ i = 1 n u ‾ i · π D 2 4 · Δ t i ( i = 1,2,3 , · · · · · · n ) - - - ( 49 )
And (40) formula, (45) formula
u = ( 1 N 1 - 1 N 2 ) · N · D / sin θ 2 cos θ · f
u = ( 1 T 1 - 1 T 2 ) · N · L 2 cos θ
The u that is obtained is actually the mean flow rate on (47) formula S disc diameter D.
By Reynolds number
R e = ρ · u ‾ i · D η - - - ( 50 )
(50) in the formula, R eIt is the Reynolds number of detected fluid; ρ is the density of detected fluid; D is detected fluid pipe 35 the internal diameter of flowing through; η is the detected fluid kinetic viscosity or claims the coefficient of viscosity; u iBe represented u in (48) formula i:
When (50) formula
R e≤2320 (51)
Detected fluid is a laminar condition.
When (50) formula
R e>2320 (52)
Detected fluid will begin to become turbulence state.
(40) formula and (45) formula are to be used for ultrasonic flow meter to measure the fluid that flows and had in the abundant developed pipe, the formula of the mean flow rate u on Fig. 1 fluid is flowed through pipe 35 diameter D.
When fluid is laminar condition
u ‾ i = 3 4 u - - - ( 53 )
When fluid is turbulence state
u ‾ i = 2 n 2 n + 1 u - - - ( 54 )
(54) n in the formula is with the different positive numbers (n of non-12 formulas of n herein) that change of Reynolds number.(53) formula and (54) formula want to set up, must (40) formula and (45) formula meet " Shannon (shannon) sampling thheorem ": if the highest frequency of time dependent simulating signal [for Fig. 1, the u in (40) formula and (45) formula] exactly is f c, as long as it is carried out enough samplings within a certain period of time, just can keep the character of this simulating signal, its condition is
T < 1 2 f c - - - ( 55 )
Wherein T is sampling period (the non-Fig. 1 of T gets the bid the concussion cycle of punctual base oscillator 19 herein), 2f cBe minimum sample frequency.
T herein is at the Δ t of this calculating neutralization (46) formula iFollowing relationship should be arranged, promptly
Δt i≤T (56)
(55) formula, (56) formula are pointed out: the positive integer N in (40) formula and (45) formula can not unconfinedly increase.Because the increase of N must cause the growth of ultrasonic propagation distance, (40) formula, (45) formula N in (40) formula and (45) formula that samples respectively like this 2, N 1Period T also must prolong, the result is false (55) formula and (56) formula.
So (55) formula and (56) formula provide condition for the value of the positive integer N in (40) formula and (45) formula.When (55), when (56) formula is set up, get
T < 1 2 f c &Delta; t i &le; T - - - ( 57 )
By (8) formula:
t 1 = o o &prime; c + u cos &theta;
t 2 = o o &prime; c - u cos &theta;
Know t 2>t 1, with t 2As the portion of time of sampling period T, the T of (57) formula is divided by the t of (8) formula 2Promptly get " prolong pass apart from technology " promptly in (40) formula, (45) formula
Figure C20051008869400204
N among the NL.
N = T t 2 - - - ( 58 )
But (8) t of formula 2Be to change with the u in (8) formula, when the variation range of u by u MinChange to u MaxThe time, its t 2By t 2minChange to t 2max, the N of (58) formula will be by N at this moment MaxChange to N MinBecause in a single day the N in (40) formula and (45) formula gets fixed, will no longer change, at this moment, if get N MaxAs the N in (40) formula and (45) formula, when rate of flow of fluid u is u MaxThe time, the t that it is corresponding 2maxWith the value N of institute MaxProduct the T of (58) formula will be increased, (55) formula is false; If get N MinAs the N in (40) formula and (45) formula, when rate of flow of fluid u is u MaxThe time, the t that it is corresponding 2maxWith the value N of institute MinProduct (55) formula will also be set up.
From the above mentioned, should to get by u in (8) formula be u to the N in (40) formula and (45) formula MaxThe time pairing t 2maxThe N that substitution (58) formula calculates MinValue after rounding (do not round up, the numeral behind the radix point all takes down) is in (40) formula and (45) formula
Figure C20051008869400206
N among the NL.In order to give calculating part single-chip microcomputer 24 with enough computing times, N can also suitably get little herein.
Through last argumentation, the N value of (40) formula and (45) formula can determine that (53) formula, (54) formula will be set up at this moment.With (40) formula, (45) formula respectively after substitution (53) formula and (54) formula, again with the u of (53) formula, (54) formula iSubstitution (49) formula,
Figure C20051008869400211
More than (39) formula, (40) formula, (41) formula, (42) formula, (45) formula, (59) formula, (60) formula, (61) formula, (62) formula, promptly
u = ( 1 n 1 - 1 n 2 ) &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; f
u = ( 1 N 1 - 1 N 2 ) &CenterDot; N &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; f
c = ( 1 n 1 + 1 n 2 ) &CenterDot; D / sin &theta; 2 &CenterDot; f
c = ( 1 N 1 + 1 N 2 ) &CenterDot; N &CenterDot; D / sin &theta; 2 &CenterDot; f
u = ( 1 T 1 - 1 T 2 ) &CenterDot; N &CenterDot; L 2 cos &theta;
Figure C20051008869400217
Figure C20051008869400218
All be τ, τ ' and the τ in (7) formula " push away under the situation about not considering, but actual using abovely when various, τ, τ ' and τ in (7) formula " must pay attention to.
Because (7) τ in the formula can basis
&tau; = 2 &times; 1 2 &tau;
= 2 &times; om c 1 (Fig. 1)
Accurately measure in advance; And τ ' and τ " also can before dispatching from the factory, accurately measure in advance by ultrasonic flow meter, so τ, τ ' and τ " all can regard known quantity as.
Now with the t in (7) formula 1, t 2(21) n of formula 1, n 2, the N of (27), (26) formula 1, N 2Respectively substitution is various about (39) formula, (40) formula, (41) formula, (42) formula, (45) formula, (59) formula, (60) formula, (61) formula and (62), can draw the formula of actual use.According to (7) formula, derive
t 1 = t 1 &prime; - &tau; - &tau; &prime; t 2 = t 2 &prime; - &tau; - &tau; &prime; &prime; - - - ( 63 )
With (63) formula substitution (21) formula
t 1·f=n 1
t 2·f=n 2
In,
n 1 = t 1 f = ( t 1 &prime; - &tau; - &tau; &prime; ) f = t 1 &prime; f - ( &tau; + &tau; &prime; ) f n 2 = t 2 f = ( t 2 &prime; - &tau; - &tau; &prime; &prime; ) f = t 2 &prime; f - ( &tau; + &tau; &prime; &prime; ) f - - - ( 64 )
If (64) in the formula
t 1 &prime; f = n 1 &prime; t 2 &prime; f = n 2 &prime; - - - ( 65 )
n 1', n 2' be respectively t in (7) formula 1', t 2The concussion number of times of the high frequency time-base generator 19 the in ' time
( &tau; + &tau; &prime; ) f = n &tau; &prime; ( &tau; + &tau; &prime; &prime; ) f = n &tau; &prime; &prime; - - - ( 66 )
n τ', n τ" be respectively (τ+τ ') in (7) formula, (τ+τ ") the concussion number of times of high frequency time-base generator 19 in the time
With (65) formula, (66) formula substitution (64) formula,
n 1 = n 1 &prime; - n &tau; &prime; n 2 = n 2 &prime; - n &tau; &prime; &prime; - - - ( 67 )
With (67) formula substitution (39) formula and (41) formula,
u = ( 1 n 1 &prime; - n &tau; &prime; - 1 n 2 &prime; - n &tau; &prime; &prime; ) &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; f - - - ( 68 )
c = ( 1 n 1 &prime; - n &tau; &prime; + 1 n 2 &prime; - n &tau; &prime; &prime; ) &CenterDot; D / sin &theta; 2 &CenterDot; f - - - ( 69 )
(68) equation the right molecule, denominator get with multiply by positive integer N in the formula
u = ( 1 n 1 &prime; - n &tau; &prime; - 1 n 2 &prime; - n 2 &prime; &prime; ) &CenterDot; N N &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; f - - - ( 70 )
(70) in the formula, establish
n 2 &prime; N = N 2 &prime; n 1 &prime; N = N 1 &prime; n &tau; &prime; &prime; N = N &tau; &prime; &prime; n &tau; &prime; N = N &tau; &prime; - - - ( 71 )
N 2', N 1', N τ", N τ' be respectively t in (7) formula 2', t 1', (τ+τ "), (τ+τ ') time enlarge the doubly concussion number of times of the high frequency time-base generator 19 of back in it of N
With (71) formula substitution (70) formula,
u = ( 1 N 1 &prime; - N &tau; &prime; - 1 N 2 &prime; - N &tau; &prime; &prime; ) &CenterDot; ND / sin &theta; 2 cos &theta; &CenterDot; f - - - ( 72 )
In like manner can get
c = ( 1 N 1 &prime; - N &tau; &prime; + 1 N 2 &prime; - N &tau; &prime; &prime; ) &CenterDot; ND / sin &theta; 2 &CenterDot; f - - - ( 73 )
And according to (63) formula both members with after multiply by N,
N &CenterDot; t 2 = ( N 2 &prime; - &tau; - &tau; &prime; &prime; ) N = t 2 &prime; &CenterDot; N - ( &tau; + &tau; &prime; &prime; ) N = T 2 N &CenterDot; t 1 = ( t 1 &prime; - &tau; - &tau; &prime; ) N = t 1 &prime; &CenterDot; N - ( &tau; + &tau; &prime; ) N = T 1 - - - ( 74 )
If (74) in the formula
t 2 &prime; &CenterDot; N = T 2 &prime; ( &tau; + &tau; &prime; &prime; ) N = T &tau; &prime; &prime; t 1 &prime; &CenterDot; N = T 1 &prime; ( &tau; + &tau; &prime; ) N = T &tau; &prime; - - - ( 75 )
After (75) formula substitution (74) formula, substitution again (45) formula gets
u = ( 1 T 1 &prime; - T &tau; &prime; - 1 T 2 &prime; - T &tau; &prime; &prime; ) &CenterDot; N &CenterDot; L 2 cos &theta; - - - ( 76 )
After (72) formula substitution (53) formula, (54) formula, substitution again (49) formula; After (76) formula substitution (53) formula, (54) formula, substitution again (49) formula,
Figure C20051008869400241
Figure C20051008869400242
More than (72) formula, (73) formula, (76) formula, (77) formula, (78) formula, (79) formula and (80) formula, τ, τ ', the τ in considering (7) formula exactly " time formula that draws, these formulas can be used when actual condition.
The relevant theory of following basis is determined N, determines that N has two kinds of ways:
First experience is determined method: " in the TT﹠C system of computer industry production run, the sampling period is indicated by sampling thheorem usually and selects f sCriterion, and determine T by actual experiment sAccording to the accumulation of protracted experience, often table look-up and select T according to the empirical data shown in the table 10-1 sFrom table 10-1 as seen, should get different sampling period T for different controlled variables s".T has been arranged s, determine N according to (58) formula again.
The empirical data table in table 10-1 sampling period
Controlled variable Sampling period T s(second) The preferential T that uses s(second)
Flow 1-5 1-2
Below rule of thumb come to determine T sTake passages in:
The 465th page of " national computer grade examination study course (three grades of A)---hardware technology and application " book, this book is published by the Electronic Industry Press, Chen Tienian, Sun Dewen chief editor, in June, 1997 first published, the first impression in June, 1997.
It two is to determine according to Shannon's sampling theorem, in this manual just can be according to (55) formula
T < 1 2 f c
(56) formula
T≤Δt i
And (58) formula
N = T t 2
To t 2Get t 2maxThe time, get N MinBe (72), (73), (76), (77), (78), (79), (80) N in various.
Detected fluid for the flow velocity acute variation, should take into full account technical dynamic requirements, in this case, also available Shannon's sampling theorem is determined N, need according to the 37th page of " flow measurement technology and instrument " book, dynamic perfromance Fig. 2 of flow instrument-4, the family curve of 1 time domain representation in the flowmeter dynamic response characteristic curve, curve is replaced with a function v=v (t), and the maximal value of its v (t) is with the measurement upper limit v of general ultrasonic flow meter MaxReplace, obtain 63v MaxThe value of %, according to the quality time constant T (non-sampling period T of the present invention of T herein or ultrasound wave period T are the T of former book content meaning) of dynamic response time domain representation method dynamic response, T is that v (t) increases to 63v from 0 then MaxThe used time of %, the more little dynamic response of T is good more, so according to the requirement of quality time constant T, can obtain concrete function v=v (t); Replenish the function variable field of definition again, make v=v (t) in (T ,+T) upward continuous, convergence; Exhibition is Fu Shi progression again, and consideration needs the overtone order of consideration according to Fu Shi progression, can obtain the 2f in (55) formula c, also can calculate its maximum sampling period T certainly according to (55) formula, the N value of Que Dinging can make ultrasonic flow meter satisfy fluid dynamic properties like this, and its dynamic response is good.
More than " flow measurement technology and instrument " book, its this instructions of situation is mentioned, no longer repeats at this.
General typical ultrasonic flow meter is measured the rate of flow of fluid scope:
Detected fluid is a gas, measures caliber (mm) 200-400, and measurable flow speed scope is 2.5-25m/s; Detected fluid is water or other liquid, measures caliber 1-15m, and measurable flow speed scope is 0.3-20m/s, and measures caliber (mm) 100-200, and measurable flow speed scope is 0.5-20m/s.(this general typical ultrasonic flow meter is measured the rate of flow of fluid scope and is selected from " airborne sensor application manual " 469 page table 12-27 ultrasonic flow sensors.This book is that China Machine Press publishes, June nineteen ninety-five first published, in October, 1997 the third printing.The printing of Seiko printing house, Changping, Beijing, the distribution of Xinhua Bookstore Beijing sale room.)
τ ' during following formula is various, τ ", in the ultrasonic flow timing of actual shop drawings 1, because of the randomness of the selected circuit devcie delay performance of each ultrasonic flow meter, so " neither meeting is identical for the τ ' of each ultrasonic flow meter, τ.But when being fabricated to finished product owing to ultrasonic flow meter, every selected circuit devcie of ultrasonic flow meter just selects, so the τ ' in various, τ " are confirmable.
Because in (72) formula, (73) formula, (76) formula, (77) formula, (78) formula, (79) formula and (80) formula, because of prolonging N, D/sin θ and L doubly carry out continuously, so N τ', N τ", T τ', T τ" can doubly determine after the whole continuous coverage in back D/sin θ and L prolongation N, must be worth the back and use for programming.
By above we can also obtain conclusion: ultrasonic flow meter of the present invention and principle thereof and technology: 1. the frequency f of the standard time-base generator 19 among Fig. 1 is high more, and its measuring error is more little, and precision is high more.2. under the prerequisite that can satisfy Shannon's sampling theorem, the sampling period is long more, and the N that the present invention " prolongs and passes apart from technology " is big more, and its measuring error is more little, and for steady flow, its precision is also high more; Instability is flowed, because the prolongation in sampling period, its dynamic response quality descends, and measuring error will increase, the also corresponding reduction of precision, but still in the sampling thheorem restricted portion.In actual applications, should get the frequency f of suitable sampling period and standard time-base generator 19, make ultrasonic flow meter dynamic response quality of the present invention and cost, can satisfy actual needs again simultaneously so that grasp according to requirement of actual working condition.
During the useful effect of the present invention, overcome known existing ultrasonic flow meter and be used to measure tubule footpath inner fluid speed hour, produce big error, the problem that its precision is greatly reduced, significantly reduce tubule footpath inner fluid speed error hour, significantly enlarge the measurement lower limit of ultrasonic flow meter, significantly improved the degree of accuracy of ultrasonic flow meter, significantly be fit to the measurement of streaming flow in all size caliber.
Description of drawings
Fig. 1 is the theory diagram of the embodiment 1 of explanation ultrasonic flow meter of the present invention and principle and technology.
Fig. 2 is the theory diagram of the embodiment 2 of explanation ultrasonic flow meter of the present invention and principle and technology.
Fig. 3 is the theory diagram of the embodiment 3 of explanation ultrasonic flow meter of the present invention (calorimeter) and principle and technology.
Fig. 4 is the theory diagram of the embodiment 4 of explanation ultrasonic flow meter of the present invention (calorimeter) and principle and technology.
Among the figure:
1, sending part
2, acceptance division
3, positive peak detector
4, sending part
5, acceptance division
6, positive peak detector
7 or the door
8, monostable
9 or the door
10, monostable
11, with door
12, with door
13, single-chip microcomputer asynchronous counter (or asynchronous counter)
14, single-chip microcomputer asynchronous counter (or asynchronous counter)
15, with door
16, bistable
17, with door
18, bistable
19, standard time-base generator
20, count section
21, count section
22, I/O interface chip
23, I/O interface chip
24, single-chip microcomputer
25, keyboard, display interface chip
26, display part
27, keypad portion
28, ultrasonic emitting oscillator
29, ultrasound wave pick-up dipole
30, ultrasonic emitting oscillator
31, ultrasound wave pick-up dipole
32, temperature-sensing element (device)
33, temperature-sensing element (device)
34, amplifier
35, the fluid pipe of flowing through
36, amplifier
37, A/D interface chip
38, ultrasonic emitting, pick-up dipole
39, ultrasonic emitting, pick-up dipole
40, electronic analog swtich
41, electronic analog swtich
Embodiment
Preferred embodiment with regard to ultrasonic flow meter of the present invention (calorimeter) and principle and technology describes below.
For the convenience of problem, establish at this:
The ultrasonic flow meter sampling period is T c
Ultrasonic flow meter is T along the circulating sampling cycle Sc, frequency is f Sc
The contrary circulating sampling cycle of ultrasonic flow meter is T Nc, frequency is f Nc
Ultrasound wave is T along cycle period s, frequency is f s
The contrary cycle period of ultrasound wave is T N, frequency is f N
T cBe the Δ t in (59), (60), (61), (62), (77), (78), (79), (80) formula i
T ScWhen not considering τ, τ ', τ " time, be in (40), (42), (59), (60) formula ultrasonic flow meter along circulating sampling N 1The value required time is the T in (45), (61), (62) formula 1, i.e. T Sc=T 1When considering τ, τ ', τ " time, be the suitable circulating sampling N of ultrasonic flow meter in (72), (73), (77), (78) formula 1' value required time is the T in (76), (79), (80) formula 1', i.e. T Sc=T 1'.Essence is to finish N ultrasound wave along the used time of circulation.
T NcWhen not considering τ, τ ', τ " time, be the contrary circulating sampling N of ultrasonic flow meter in (40), (42), (59), (60) formula 2The value required time is the T in (45), (61), (62) formula 2, i.e. T Nc=T 2When considering τ, τ ', τ " time, be the contrary circulating sampling N of ultrasonic flow meter in (72), (73), (77), (78) formula 2' value required time is the T in (76), (79), (80) formula 2', i.e. T Nc=T 2'.Essence is to finish the contrary used time of circulation of ultrasound wave N time.
What be worth proposition is that above-mentioned each amount all is functions of the fluid measured flow velocity u of institute.
Preferred embodiment with regard to ultrasonic flow meter of the present invention (calorimeter) and principle and technology describes below
(embodiment 1)
Fig. 1 is the theory diagram of ultrasonic flow meter of the present invention and principle and technology implementation example 1.
Embodiment 1 ultrasonic flow meter comprises: sending part 1, acceptance division 2, positive peak detector 3, sending part 4, acceptance division 5, positive peak detector 6, or door 7, monostable 8, or door 9, monostable 10, with door 11, with door 12, single-chip microcomputer asynchronous counter 13, single-chip microcomputer asynchronous counter 14, with door 15, bistable 16, with door 17, bistable 18, standard time-base generator 19, count section 20, count section 21, I/O interface chip 22, I/O interface chip 23, single-chip microcomputer 24, keyboard and display interface chip 25, display part 26, keypad portion 27, ultrasonic emitting oscillator 28, ultrasound wave pick-up dipole 29, ultrasonic emitting oscillator 30, ultrasound wave pick-up dipole 31, fluid is flowed through and is managed 35.
Ultrasonic flow meter shown in Figure 1 is broadly divided into three parts, sampling section, calculating section and keyboard display part.
Sampling section is made of jointly count section among Fig. 1 20,21 two count section of count section (be all asynchronous counter, its figure place should satisfy sample requirement and surpass demand) and each circuit block of the left side thereof.
Calculating section is made up of single-chip microcomputer 24.
Sampling section is connected by I/O interface chip 22, I/O interface chip 23 with calculating section.
In sampling section, standard time-base generator 19 is not controlled by remaining part and rate of flow of fluid in Fig. 1 circuit, as long as connect
Energize just ceaselessly vibrates, and output is higher than the high frequency square wave spike train of ultrasound wave along circulation and the contrary cycle frequency of ultrasound wave far away.
Sampling section is sampled or is not sampled, by the control assembly bistable 16 in the sampling section, bistable 18, monostable 8, monostable 10, with door 11, with door 12, with door 15, control (all with door or door all band application schmitt trigger) with door 17 or door 7 or door 9, and these control assemblies directly or indirectly are subjected to the control of the single-chip microcomputer 24 of single-chip microcomputer asynchronous counter 13 (or asynchronous counter), single-chip microcomputer asynchronous counter 14 (or asynchronous counter) and calculating section in the sampling section selectively.
The beginning in each sampling period must be caused by the I/O mouth B end output one high level square-wave signal of calculating section single-chip microcomputer 24, controls control assembly bistable 16, bistable 18 set in the sampling section simultaneously, even its Q end is all high level.The Q of bistable 18 end high level open with door 17 and with door 12, with door 17 the high-frequency impulse series of standard time-base generator 19 is counted by entering count section 21 with door 17; Make with door 12: the high level signal of the I/O mouth B of single-chip microcomputer 24 end monostable 10 produce an output positive pulse signal by or door 9 allowed by being input to and activating sending part 4 with door 12; And the Q of bistable 16 end high level open with door 15 and with door 11, with door 15 the high level pulse series of standard time-base generator 19 is counted by entering count section 20 with door 15, is made with door 11: the high level signal of the I/O mouth B of single-chip microcomputer 24 end monostable 8 produce an output positive pulse signal by or door 7 allowed by being input to and activating sending part 1 with door 11.
Sending part 4 drives countercurrent direction ultrasonic emitting oscillator 30 emission ultrasound waves and passes vu-uu '-u ' v ' path, received by adverse current ultrasound wave pick-up dipole 31, again after acceptance division 5 amplifies, detected the positive crest moment and export a high level signal by positive peak detector 6, be input to respectively or door 9 and single-chip microcomputer asynchronous counter 14, be input to or the signal of door 9 by or door 9, again by with door 12 to sending part 4, form the contrary circulation of ultrasound wave, establishing the contrary cycle frequency of its ultrasound wave is f N, the contrary cycle period of ultrasound wave is T NEvery circulation primary is amplified through acceptance division 5 and to be input to after positive peak detector 6 detects positive crest, outputs to the signal of single-chip microcomputer asynchronous counter 14, just by single-chip microcomputer asynchronous counter 14 countings once.In single-chip microcomputer asynchronous counter 14, insert the positive integer N (being about to the N value enrolls in single-chip microcomputer 14 programs) among the present invention's " prolongation passes apart from technology ", when making single-chip microcomputer asynchronous counter 14 count down to N, interruption is overflowed in generation, export a high level signal to bistable 18 and calculating section single-chip microcomputer 24, this signal make bistable 18 reset its output Q be low level, will with door 17, close the door with door 12, the contrary circulation of above-mentioned ultrasound wave is stopped, count section 21 also stops counting, make the contrary circulating sampling end cycle of ultrasonic flow meter, after this signal is input to the A of calculating section single-chip microcomputer 24, cause that single-chip microcomputer 24 enters interruption subroutine.
Sending part 1 drives downbeam ultrasonic emitting oscillator 28 emission ultrasound waves and passes mo-oo '-o ' m ' path, received by following current ultrasound wave pick-up dipole 29, again after acceptance division 2 amplifies, detected the positive crest moment and export a high level signal by positive peak detector 3, be input to respectively or door 7 and single-chip microcomputer asynchronous counter 13, be input to or the signal of door 7 by or door 7, again by with door 11 to sending part 1, form ultrasound wave along circulation, establishing its ultrasound wave is f along cycle frequency s, ultrasound wave is T along cycle period sEvery circulation primary, acceptance division 2 are amplified and are input to after positive peak detector 3 detects positive crest, output to the signal of single-chip microcomputer asynchronous counter 13, just by single-chip microcomputer asynchronous counter 13 countings once.In single-chip microcomputer asynchronous counter 13, insert the positive integer N among the present invention's " prolongation passes apart from technology " equally, when making single-chip microcomputer asynchronous counter 13 count down to N, interruption is overflowed in generation, export a high level signal to bistable 16 reset its output Q be low level, will be with door 15,11 close the door with door, above-mentioned ultrasound wave is stopped along circulation, and count section 20 also stops counting, makes ultrasonic flow meter along the circulating sampling end cycle.
Because the suitable circulating sampling cycle of ultrasonic flow meter is less than the contrary circulating sampling cycle, so it is the basic cycle that the ultrasonic flow meter sampling period also just should be got the contrary circulating sampling cycle, Here it is, and why single-chip microcomputer asynchronous counter 13 and single-chip microcomputer asynchronous counter 14 can both produce and overflow interruption, and have only single-chip microcomputer asynchronous counter 14 overflow the A that look-at-me is input to calculating section single-chip microcomputer 24 after, can cause that just single-chip microcomputer 24 enters the reason of interruption subroutine.
According to foregoing description, can find out:
1. the beginning and to cause by the I/O mouth B end output one high level square-wave signal of the single-chip microcomputer 27 of calculating section of each sampling period.
The end in 2. each sampling period produces after should the full N of single-chip microcomputer asynchronous counter 13 countings by sampling section to overflow to be interrupted and finishes.
3. count section 20 is counted the N in (72) formula that is, (77) formula, (78) formula 1' (81)
4. count section 21 is counted the N in (72) formula that is, (77) formula, (78) formula 2' (82)
5. the frequency f of high frequency time-base generator 19>>f s, the frequency f of high frequency time-base generator 19>>f N
1 N &CenterDot; f s = f sc 1 N &CenterDot; f N = f Nc (N is the positive integer during the present invention " prolongs and passes apart from technology ")
f Sc, f NCShould satisfy (55) formula respectively T < 1 2 f c Or the empirical data table 0-1 in sampling period
Be f Sc>2f cf NC>2f c
(72) N in formula, (77) formula, (78) formula τ", N τ' must be worth after can accurately measuring, that is:
N &tau; &prime; &prime; = n &tau; &prime; &prime; &CenterDot; N = ( &tau; + &tau; &prime; &prime; ) &CenterDot; f &CenterDot; N N &tau; &prime; = n &tau; &prime; &CenterDot; N = ( &tau; + &tau; &prime; ) &CenterDot; f &CenterDot; N - - - ( 83 )
(83) get behind the whole measured value of formula.N τ", N τ' value is used for single-chip microcomputer 24 programming uses.
(72) N in formula, (77) formula, (78) formula is except that single-chip microcomputer asynchronous counter 13 that should deposit sampling section in sampling process in and single-chip microcomputer asynchronous counter 14, also should with D, sin θ in N value and (72) formula, (77) formula, (78) formula three formulas, 2cos θ,
Figure C20051008869400295
Carry out following calculating,
ND / sin &theta; 2 cos &theta; &CenterDot; &pi; D 2 4 - - - ( 84 )
(84) be used for after the formula overall calculation must be worth single-chip microcomputer 24 programming uses.
(77) in the formula To detected fluid is laminar flow (laminar flow) or Reynolds number R eCalculated Q at<2320 o'clock tThe time, be used for single-chip microcomputer 24 programming uses.
(78) in the formula N can get by following column count:
By the 28th page of (2-36) formula of " flow measurement technology and instrument " book
u &OverBar; = 2 n 2 n + 1 &CenterDot; u &OverBar; D
With the 24th page of (2-23) formula
R e = u &CenterDot; L &CenterDot; &rho; &eta; ( = &rho; &CenterDot; u &OverBar; &CenterDot; D &eta; )
With (2-23) formula of this book of this book (2-36) formula substitution,
R e = &rho; &CenterDot; D &eta; &CenterDot; 2 n 2 n + 1 &CenterDot; u D - - - ( 85 )
By 26 page table 2-12 of this book, the relation of n and Reynolds number as can be known, when under the constant situation of D, ρ, η, R eWhen known, n is corresponding known, that is to say the u in (85) formula DCan only get unique value.Conversely speaking, if under the constant situation of D, ρ, η, u DSize can determine the size of n, its u DWith n is (also to be u one to one DWith Be one to one), and u DThe u of the present invention just (72) formula.
By (51) formula R eIt≤2320 o'clock was laminar condition
Figure C20051008869400307
By (85) formula u DPromptly the n of the u correspondence of (72) formula or
Figure C20051008869400308
Be used for single-chip microcomputer 24 programming uses.
(72), the Δ t in (77), (78) formula iBe ultrasonic flow meter sampling period T CFollowing calculating:
Δ t i=N 2' T+ (interrupt+forbid all interruption+single-chip microcomputer 24 output signals and latch N by single-chip microcomputer 24 responses 1' and N 2'+single-chip microcomputer 24 output signals begin sampling to count section 20 and 21 zero clearings+single-chip microcomputer 24 from I/O B end output signal) instruct the used time (86)
Above-mentioned data have been arranged, that is:
(81) several N of formula 1'
(82) several N of formula 2'
(83) several N of formula τ", N τ'
(84) number of formula ND / sin &theta; 2 cos &theta; &CenterDot; &pi; D 2 4
(85) several n of formula or
Figure C20051008869400311
(86) several Δ t of formula i
And constant
Figure C20051008869400312
And f
Just can the single-chip microcomputer 24 of calculating section be programmed: make single-chip microcomputer 24 after entering interrupt routine, forbid all interruptions, latch the sampled value N of count section 20 and 21 from the F output signal 1' and N 2', to count section 20 and 21 zero clearings,, make sampling section enter the next sampling period from the C output signal from high level of I/O mouth B end output, read in several N of I/ O interface chip 22 and 23 respectively from D and E 1' and N 2' two numbers, and in single-chip microcomputer 24, finish Q of the calculating of (72) formula u and (77), (78) formula with above-mentioned existing all the other relevant each numbers tCalculating, u and Q are calculated in institute tValue deposits in single-chip microcomputer 24 nonvolatile memorys, changes master routine over to.Master routine is mainly finished open interruption and is sent apparent task.
The operation of ultrasonic flow meter, stop, the modification of above-mentioned relevant constant undertaken by the keyboard input command.
(embodiment 2)
Fig. 2 is the theory diagram of the embodiment 2 of ultrasonic flow meter of the present invention and principle and technology.
The ultrasonic flow meter of embodiment 2 comprises: sending part 1, acceptance division 2, positive peak detector 3, or door 7, monostable 8, with door 11, single-chip microcomputer asynchronous counter 13, with door 15, bistable 16, standard time-base generator 19, count section 20, I/O interface chip 22, single-chip microcomputer 24, keyboard and display interface chip 25, display part 26, keypad portion 27, fluid is flowed through and is managed 35, ultrasonic emitting pick-up dipole 38, ultrasonic emitting pick-up dipole 39, electronic analog swtich 40, electronic analog swtich 41.
The ultrasonic flow meter of embodiment 2 is broadly divided into three parts; Sampling section, calculating section and keyboard display part.
Sampling section is made jointly by the count section among Fig. 2 20 (asynchronous counter) and each circuit block of the left side thereof, and calculating section is made up of single-chip microcomputer 24, and sampling section is connected by I/O interface chip 22 with calculating section.
In the sampling section, the situation of the ultrasonic flow meter of standard time-base generator 19 and embodiment 1 is identical.
Sampling section is sampled or is not sampled, by the control assembly bistable 16 in the sampling section, monostable 8, with door 11, control (all with door or door all band application schmitt trigger) with door 15 or door 7, and these control assemblies selectively directly or indirectly are subjected to the control of sampling portion single-chip microcomputer asynchronous counter 13 and calculating section single-chip microcomputer 24.
The ultrasonic flow meter of embodiment 2, ultrasound wave is served as emission ultrasound wave task along circulation time by ultrasonic emitting pick-up dipole 39, and ultrasonic emitting pick-up dipole 38 is served as ultrasound wave and received task, this moment, a and the c of electronic analog swtich 40 connected a ' of electronic analog swtich 41 and c ' connection; And the contrary circulation time of ultrasound wave is served as emission ultrasound wave task by ultrasonic emitting pick-up dipole 38, and ultrasonic emitting pick-up dipole 39 is served as ultrasound wave and received task, and this moment, a and the b of electronic analog swtich 40 connected a ' of electronic analog swtich 41 and b ' connection.
Ultrasonic flow meter shown in Figure 2, because ultrasound wave is along circulation time and the contrary circulation time of ultrasound wave, the circuit of all having used identical sample unit to form, so ultrasonic flow meter has only timesharing to carry out along circulating sampling and the contrary circulating sampling of ultrasonic flow meter, this puts different with the ultrasonic flow meter of embodiment 1 shown in Figure 1.
Because ultrasonic flow meter is that timesharing is carried out along circulating sampling and the contrary circulating sampling of ultrasonic flow meter, so ultrasonic flow meter is inserted the positive integer of single-chip microcomputer asynchronous counter 13 " prolongation passes apart from technology " and is made as N ', can only be half of positive integer N of the ultrasonic flow meter of embodiment 1 " prolong and pass " of inserting single-chip microcomputer asynchronous counter 13 and 14, promptly apart from technology N &prime; = N 2 , This puts different with the ultrasonic flow meter of embodiment 1 shown in Figure 1.
The ultrasonic flow meter of embodiment 2 is because its positive integer that " prolongs and pass apart from technology " N &prime; = N 2 , So (83), (84) formula should be written as (83 ') and (84 ')
N &tau; &prime; &prime; = n &tau; &prime; &prime; &CenterDot; N = ( &tau; + &tau; &prime; &prime; ) &CenterDot; f &CenterDot; N 2 N &tau; &prime; = n &tau; &prime; &CenterDot; N = ( &tau; + &tau; &prime; ) &CenterDot; f &CenterDot; N 2 - - - ( 8 3 &prime; )
N &prime; D / sin &theta; 2 cos &theta; &CenterDot; &pi; D 2 4 = N 2 &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; &pi; D 2 4 - - - ( 8 4 &prime; )
This puts different with the ultrasonic flow meter of embodiment 1 shown in Figure 1.
The ultrasonic flow meter of embodiment 2 is because ultrasonic flow meter is that timesharing is carried out along circulating sampling and the contrary circulating sampling of ultrasonic flow meter, so the sampling period of ultrasonic flow meter
Δ t i=(N 1'+N 2') T+2 * (single-chip microcomputer 24 responses are interrupted+forbidden all interruption+single-chip microcomputer 24 output signals and latch N 1' and N 2'+single-chip microcomputer 24 output signals begin sampling to count section 20 zero clearings+single-chip microcomputer 24 from I/O B end output signal) instruct the used time
(86′)
This puts different with the ultrasonic flow meter of embodiment 1 shown in Figure 1.
The ultrasonic flow meter of embodiment 2 shown in Figure 2, the beginning in each sampling period must be by first high level square-wave signal of I/O B end output of calculating section single-chip microcomputer 24, link simultaneously electronic analog swtich 40 and 41, make its connect for ultrasound wave along circulation or the contrary round-robin form of ultrasound wave.I/O mouth B holds this high level signal to be input to the set end and monostable 8 of bistable 16, be input to the asserts signal of bistable 16, make bistable 16 put the Q end and be high level, open with the door 15 and with door 11, make the high level pulse of standard time-base generator 19 with door 15, count by entering count section 20 with door 15, make with door 11: the high level of I/O mouth B end be input to monostable 8 make monostable 8 produce positive pulse output by or door 7 backs by being input to sending part 1 with door 11, activate sending part 1 and drive ultrasonic emitting pick-up dipole 39 or 38 emission ultrasound waves, and in ultrasonic emitting pick-up dipole 39 or 38 emissions hyperacoustic the time, ultrasonic emitting pick-up dipole 38 or 39 just receives ultrasound wave, detect positive crest constantly by being input to positive peak detector 3 after acceptance division 2 amplifications, export a high level to or door 7 and single-chip microcomputer asynchronous counter 13, be input to or door 7 signal by or door 7, again by being input to sending part 1 with door 11, make sampling section enter ultrasound wave along circulation or contrary circulation, ultrasound wave is along circulation or contrary circulation N &prime; = N 2 Inferior, full by single-chip microcomputer asynchronous counter 13 countings N &prime; = N 2 Overflow a high level signal, be input to bistable 16 and single-chip microcomputer 24A end, making bistable 16 Q that resets is low level, will close with door 11 with door 15.With closing of door 15 counting of count section 20 is stopped, ultrasonic flow meter is along circulating sampling cycle and the contrary circulating sampling end cycle of ultrasonic flow meter.The high level that is input to single-chip microcomputer 24A end causes that single-chip microcomputer 24 interrupts, and single-chip microcomputer 24 is finished the acceptance interruption, forbidden all interruptions, and single-chip microcomputer 24 latchs the count value N of count section 20 from the F output signal 1' or N 2', after C output signal zero clearing count section 20, subsequently, from second high level square-wave signal of I/O mouth B end output, link simultaneously electronic analog swtich 40 and 41, make its connect for and the opposite connection of ultrasound wave circulation last time, allow circuit finish the opposite ultrasonic flow meter sampling period.
By a pair of so opposite ultrasonic flow meter circulating sampling cycle and, ultrasonic flow meter sampling period of basic composition, later situation recursion.
From the above mentioned, the ultrasonic flow meter of embodiment 2 shown in Figure 2 can basis:
(81) N of formula 1' (ultrasonic flow meter is along the 20 sampled value N that count of circulating sampling cycle inside counting portion 1')
(82) N of formula 2' (the contrary 20 sampled value N that count of circulating sampling cycle inside counting portion of ultrasonic flow meter 2')
The N of (83 ') formula τ", N τ'
(84 ') formula N &prime; D / sin &theta; 2 cos &theta; &CenterDot; &pi; D 2 4 = N 2 &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; &pi; D 2 4 - - - ( 8 4 &prime; )
(85) n of formula or
Figure C20051008869400332
(86 ') formula Δ t i=(N 1'+N 2') T+2 * (single-chip microcomputer 24 responses are interrupted+forbidden all interruption+single-chip microcomputer 24 output signals and latch N 1' and N 2+ single-chip microcomputer 24 output signals begin sampling to count section 20 zero clearings+single-chip microcomputer 24 from I/O B end output signal) instruct the used time
And constant
Figure C20051008869400333
And f
Just can the single-chip microcomputer 24 of calculating section be programmed the same task that the ultrasonic flow meter that makes the ultrasonic flow meter of embodiment shown in Figure 22 finish embodiment shown in Figure 11 is finished.
(embodiment 3)
Fig. 3 is the theory diagram of the embodiment 3 of ultrasonic flow meter of the present invention (calorimeter) and principle and technology.
The ultrasonic flow meter of embodiment 3 (calorimeter) comprising: ultrasonic flow meter shown in Figure 1, temperature-sensing element (device) 32, temperature-sensing element (device) 33, amplifier 34, amplifier 36, A/D interface chip 37.
Structure, principle, technology, measuring process and the embodiment 1 described ultrasonic flow meter of the ultrasonic flow meter part in the ultrasonic flow meter of embodiment 3 (calorimeter) are identical.
Temperature-sensing element (device) 32 in the ultrasonic flow meter of embodiment 3 (calorimeter), temperature-sensing element (device) 33, amplifier 34, amplifier 36, A/D interface chip 37 and single-chip microcomputer 24, keyboard and display interface chip 25, display part 26, keyboard 27 are formed temperature gauge.
The fluid temperature (F.T.) that temperature gauge is located its place by temperature-sensing element (device) 32 and temperature-sensing element (device) 33 transfers this two signal of electric signal to respectively by after amplifier 34 and amplifier 36 amplifications, send into A/D interface chip 37 and carry out analog to digital conversion, digital signal after the conversion is input to single-chip microcomputer 24, according to these two numerals, fluid temperature (F.T.) by 24 pairs of temperature-sensing element (device)s 32 of single-chip microcomputer and place, temperature-sensing element (device) 33 places is calculated, and establishing its calculated value is t 1And t 2, and the N that samples in the ultrasonic flow meter in the moment of these two temperature value samplings and this machine 1' and N 2' moment unanimity (on the sequential before and after face mutually).
By the measured fluid volume Q of the ultrasonic flow meter (calorimeter) of embodiment shown in Figure 33 tThe fluid temperature (F.T.) t measured with said temperature meter 2 1And t 2, and fluid known density ρ, but the fluid between convection cell upstream temperature sensitive element 33 and the downstream temperature sensitive element 32 just flow through by fluid and to manage 35 outside heat dissipation capacity J sOr caloric receptivity J xCalculate.
Its calculating formula:
J s=cm t(t 2-t 1)
=c·Q t·ρ·(t 2-t 1)
J x=cm t(t 1-t 2)
=c·Q t·ρ·(t 1-t 2)
C is a detected fluid specific heat in the formula, m tBe that volume is Q tThe quality of fluid.
The principle formula of ultrasonic flow meter and principle thereof and technology and above-mentioned J according to the present invention sOr J x, single-chip microcomputer 24 is programmed, can realize detected fluid flow velocity u, flow Q t, fluid J sWith J xMeasurement.
(embodiment 4)
Fig. 4 is the theory diagram of the embodiment 4 of ultrasonic flow meter of the present invention (calorimeter) and principle and technology.
The ultrasonic flow meter of embodiment 4 (calorimeter) comprising: ultrasonic flow meter shown in Figure 2 and temperature gauge shown in Figure 3.
The structure of the ultrasonic flow meter in the ultrasonic flow meter of embodiment 4 (calorimeter), principle, technology, measuring process and embodiment 2 described ultrasonic flow meters are identical; The structure of the temperature gauge among the embodiment 4, principle, technology, measuring process and embodiment 3 described temperature gauges are identical.
By the fluid measured volume Q of ultrasonic flow meter institute in the ultrasonic flow meter (calorimeter) of embodiment shown in Figure 44 tThe fluid temperature (F.T.) t measured with the temperature gauge of this machine 1And t 2, and fluid known density ρ, the principle formula of ultrasonic flow meter and principle thereof and technology, above-mentioned J according to the present invention sWith J xThe principle formula, single-chip microcomputer 24 is programmed the same task that the ultrasonic flow meter (calorimeter) that makes the ultrasonic flow meter (calorimeter) of embodiment shown in Figure 44 finish embodiment shown in Figure 33 is finished.

Claims (2)

1, a kind of ultrasonic flow meter is measured the method for fluid flow, wherein ultrasonic flow meter comprises sampling section, calculating section, keyboard display part, it is characterized in that: by detecting ultrasonic flow meter respectively along circulating sampling in the cycle and the concussion number of times of the contrary circulating sampling cycle internal standard time-base generator (19) of ultrasonic flow meter, i.e. the sample count value N ' of sampling section first count section (20) and second count section (21) 1, N ' 2, in the single-chip microcomputer (24) of calculating section, the measuring principle of utilization ultrasonic flow meter flow velocity and flow and " prolong and pass apart from technology " are calculated, and realize the measurement of rate of flow of fluid and flow;
The measuring principle of described ultrasonic flow meter flow velocity and flow is:
(1) do not consider fluid flow through pipe (35) upside and the used time τ of downside, along the cycling circuit delay time T ' and contrary cycling circuit delay time T " time, the measuring principle formula of rate of flow of fluid u and the ultrasound wave velocity of wave c in fluid is (39), (41) formula:
u = ( 1 n 1 - 1 n 2 ) &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; f - - - ( 39 )
c = ( 1 n 1 + 1 n 2 ) &CenterDot; D / sin &theta; 2 &CenterDot; f - - - ( 41 )
n 1: be that the ultrasound wave following current is propagated oo '=D/sin θ apart from required time t 1The number of oscillation of interior standard time-base generator (19);
n 2: be that the ultrasound wave adverse current is propagated uu '=oo '=D/sin θ apart from required time t 2The number of oscillation of interior standard time-base generator (19);
D: be the flow through interior diameter of pipe (35) of fluid;
θ: the angle of its direction and rate of flow of fluid when being ultrasound wave following current propagation;
F: the oscillation frequency that is standard time-base generator (19);
(2) consider fluid flow through pipe (35) upside and the used time τ of downside, suitable cycling circuit delay time T ' and against the cycling circuit delay time T " time, the measuring principle formula of rate of flow of fluid u and the ultrasound wave velocity of wave c in fluid is (68), (69) formula:
u = ( 1 n 1 &prime; - n &tau; &prime; - 1 n 2 &prime; - n &tau; &prime; &prime; ) &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; f - - - ( 68 )
c = ( 1 n 1 &prime; - n &tau; &prime; + 1 n 2 &prime; - n &tau; &prime; &prime; ) &CenterDot; D / sin &theta; 2 &CenterDot; f - - - ( 69 )
N ' 1: be t ' 1=t 1The concussion number of times of the standard time-base generator (19) in+τ+τ ' time;
N ' 2: be t ' 2=t 2+ τ+τ " the concussion number of times of the standard time-base generator (19) in the time;
N ' τ: the concussion number of times that is the standard time-base generator (19) in τ+τ ' time;
N " τ: " the concussion number of times of the standard time-base generator (19) in the time that is τ+τ;
To " prolong and pass " measuring principle formula (40), (42) formula and (72), (73) formula that measuring principle formula (39), (41) formula and (68), (69) formula obtain rate of flow of fluid and velocity of wave respectively of being applied to respectively apart from technology;
Wherein when not considering τ, τ ', τ " time, the measuring principle formula of rate of flow of fluid u and the ultrasound wave velocity of wave c in fluid is:
u = ( 1 N 1 - 1 N 2 ) &CenterDot; N &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; f - - - ( 40 )
c = ( 1 N 1 + 1 N 2 ) &CenterDot; N &CenterDot; D / sin &theta; 2 &CenterDot; f - - - ( 42 )
N 1: be that the ultrasound wave following current is propagated oo '=D/sin θ apart from enlarging N times of back required time t 1The number of oscillation of the standard time-base generator (19) in the N;
N 2: be that the ultrasound wave adverse current is propagated uu '=oo '=D/sin θ apart from enlarging N times of back required time t 2The number of oscillation of the standard time-base generator (19) in the N;
N: the positive integer N value that is " prolong and pass " apart from technology;
When considering τ, τ ', τ " time, the measuring principle formula of rate of flow of fluid u and the ultrasound wave velocity of wave c in fluid is:
u = ( 1 N 1 &prime; - N &tau; &prime; - 1 N 2 &prime; - N &tau; &prime; &prime; ) &CenterDot; ND / sin &theta; 2 cos &theta; &CenterDot; f - - - ( 72 )
c = ( 1 N 1 &prime; - N &tau; &prime; + 1 N 2 &prime; - N &tau; &prime; &prime; ) &CenterDot; ND / sin &theta; 2 &CenterDot; f - - - ( 73 )
N ' 1: be t ' 1Time enlarges the doubly concussion number of times of its interior standard time-base generator (19) of back of N;
N ' 2: be t ' 2Time enlarges the doubly concussion number of times of its interior standard time-base generator (19) of back of N;
N ' τ: be to enlarge the doubly concussion number of times of its interior standard time-base generator (19) of back of N τ+τ ' time;
N " τ: being τ+τ, " time enlarges the doubly concussion number of times of its interior standard time-base generator (19) of back of N;
Described " prolong and pass apart from technology " is:
When ultrasonic flow meter is measured rate of flow of fluid u, prolong the propagation distance of ultrasound wave in fluid, be called for short " prolong and pass " apart from technology, the implementation method of this technology is respectively to (39) formula, (41) formula, (68) molecule of the equation of formula, (69) formula the right formula and denominator are with multiply by a positive integer N, make the propagation distance D/sin θ of ultrasound wave in limited caliber enlarge N doubly, thereby obtain measuring principle formula (40) formula, (42) formula, (72) formula, (73) formula of flow velocity and velocity of wave respectively;
And use (48), (49), (50), (51), (53), (52), (54), (85) formula to draw flow measurement principle formula (59), (60) formula and (77), (78) formula, that is:
When not considering τ, τ ', τ " time, flow measurement principle formula is:
Figure C2005100886940003C4
Q t: be that fluid is flowed through in the pipe (35) S = &pi; D 2 4 The fluid flow of elapsed time t on the represented S area;
Q i: be Λ t iIn time, S = &pi; D 2 4 Flow on the represented S area;
u i: be Λ t iThe mean flow rate of the fluid in the time on the above-mentioned S area;
π: be the garden frequency;
Δ t i: be the ultrasonic flow meter sampling period;
Figure C2005100886940004C1
Figure C2005100886940004C2
In n be with the different positive numbers that change of the Reynolds number of detected fluid;
R e: the Reynolds number that is detected fluid;
When considering τ, τ ', τ " time, flow measurement principle formula is:
Figure C2005100886940004C3
" prolonging and pass apart from technology " the positive integer N that is taken advantage of can not unrestrictedly increase, should be subjected to the constraint of " Shannon's sampling theorem " determined (55) formula, determining of its positive integer N value should be determined T earlier in conjunction with reality from (55) formula, uses the Peak Flow Rate u when measuring rate of flow of fluid then from (8) formula MaxDetermine t 2max, use (58) formula to get N again MinPositive integer value be the positive integer N value of " prolong pass apart from technology ",
Wherein (55) formula, (8) formula, (58) formula and (48), (49), (50), (51), (53), (52), (54), (85) formula are as follows:
T < 1 2 f c - - - ( 55 )
T: be " Shannon's sampling theorem " determined sampling period;
f c: be the mean flow rate u of fluid on pipe (35) the diameter D that flows through DThe highest frequency of time dependent simulating signal;
t 1 = o o &prime; c + u cos &theta; t 2 = o o &prime; c - u cos &theta; - - - ( 8 )
N min = T t 2 max - - - ( 58 )
t 2max: be the corresponding u of u MaxThe time correspondence t 2Maximal value;
N Min: N MinPositive integer value be the positive integer N value of " prolong pass apart from technology ";
Q i = u i - &CenterDot; &pi; D 2 4 &CenterDot; &Delta; t i - - - ( 48 )
Q t = &Sigma; i = 1 n Q i = &Sigma; i = 1 n u i - &CenterDot; &pi; D 2 4 &CenterDot; &Delta; t i , ( i = 1,2,3 , &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; &CenterDot; n ) - - - ( 49 )
R e = &rho; &CenterDot; u - i &CenterDot; D &eta; - - - ( 50 )
ρ: the density that is detected fluid;
η: be the detected fluid kinetic viscosity or the title coefficient of viscosity;
When
R e≤ 2320 o'clock (51)
u i - = 3 4 u - - - ( 53 )
When
R e>2320 o'clock (52)
u i - = 2 n 2 n + 1 u - - - ( 54 )
Wherein:
R e = &rho; &CenterDot; D &eta; &CenterDot; 2 n 2 n + 1 &CenterDot; u D - - - ( 85 )
u D=u: be flow through mean flow rate on pipe (35) interior diameter of fluid;
The derivation of the measuring principle formula of described ultrasonic flow meter flow velocity and velocity of wave is:
When not considering τ, τ ', τ " time, that is:
c=f·λ (11)
u &lambda; = n * - - - ( 12 )
λ: be the ultrasound wave wavelength of ultrasound wave when the f frequency;
n *Be the multiple of u to λ;
c + u cos &theta; = f&lambda; + n * &lambda; cos &theta;
= f ( 1 + n * cos &theta; f ) &lambda; - - - ( 13 )
c - u cos &theta; = f&lambda; - n * &lambda; cos &theta;
= f ( 1 - n * cos &theta; f ) &lambda; - - - ( 14 )
f ( 1 + n * cos &theta; f ) &lambda; = f &lambda; 1 - - - ( 15 )
f ( 1 - n * cos &theta; f ) &lambda; = f &lambda; 2 - - - ( 16 )
λ 1: be the wavelength of ultrasound wave downbeam in fluid when propagating;
λ 2: be the wavelength of ultrasound wave countercurrent direction in fluid when propagating;
2ucosθ=f(λ 12) (19)
2c=f(λ 12) (34)
With (13), (14), (15), (16) formula substitution formula (8) formula, obtain (20), (21) formula:
t 1 f = o o &prime; &lambda; 1 t 2 f = o o &prime; &lambda; 2 - - - ( 20 )
t 1 f = n 1 t 2 f = n 2 - - - ( 21 )
(20), (21) formula substitution (19), (34) formula are got measuring principle formula (39), (41) formula of flow velocity and velocity of wave;
When considering τ, τ ', τ " time, that is:
With the t in (63) formula 1, t 2In substitution (21) formula, get (64) formula,
Wherein (63) formula, (64) formula are as follows:
t 1 = t 1 &prime; - &tau; - &tau; &prime; t 2 = t 2 &prime; - &tau; - &tau; &prime; &prime; - - - ( 63 )
n 1 = t 1 f = ( t 1 &prime; - &tau; - &tau; &prime; ) f = t 1 &prime; f - ( &tau; + &tau; &prime; ) f n 2 = t 2 f = ( t 2 &prime; - &tau; - &tau; &prime; &prime; ) f = t 2 &prime; f - ( &tau; + &tau; &prime; &prime; ) f - - - ( 64 )
If (64) in the formula
t 1 &prime; f = n 1 &prime; t 2 &prime; f = n 2 &prime; - - - ( 65 )
( &tau; + &tau; &prime; ) f = n &tau; &prime; ( &tau; + &tau; &prime; &prime; ) f = n &tau; &prime; &prime; - - - ( 66 )
With (65) formula, (66) formula substitution (64) formula,
n 1 = n 1 &prime; - n &tau; &prime; n 2 = n 2 &prime; - n &tau; &prime; &prime; - - - ( 67 )
With (67) formula difference substitution (39) formula and (41) formula, get measuring principle formula (68), (69) formula of flow velocity and velocity of wave, equation the right molecule, denominator get (70) formula with multiply by positive integer N in (68), (69) formula,
u = ( 1 n 1 &prime; - n &tau; &prime; - 1 n 2 &prime; - n &tau; &prime; ) &CenterDot; N N &CenterDot; D / sin &theta; 2 cos &theta; &CenterDot; f c = ( 1 n 1 &prime; - n &tau; &prime; + 1 n 2 &prime; - n &tau; &prime; &prime; ) &CenterDot; N N &CenterDot; D / sin &theta; 2 &CenterDot; f - - - ( 70 )
(70) in the formula, establish
n 2 &prime; N = N 2 &prime; n 1 &prime; N = N 1 &prime; n &tau; &prime; &prime; N = N &tau; &prime; &prime; n &tau; &prime; N = N &tau; &prime; - - - ( 71 )
With (71) formula substitution (70) formula, get fluid-velocity survey principle formula (72) formula and velocity of wave measuring principle formula (73) formula;
N in flow velocity, velocity of wave measuring principle formula (72), (73) formula, flow measurement principle formula (77), (78) formula " τ, N ' τUse with being used for single-chip microcomputer (24) programming behind the whole measured value of (83) formula;
N &tau; &prime; &prime; = n &tau; &prime; &prime; &CenterDot; N = ( &tau; + &tau; &prime; &prime; ) &CenterDot; f &CenterDot; N N &tau; &prime; = n &tau; &prime; &CenterDot; N = ( &tau; + &tau; &prime; ) &CenterDot; f &CenterDot; N - - - ( 83 )
Δ t in flow measurement principle formula (59), (60), (77), (78) formula iNot only should satisfy (56) formula
Δt i≤T (56)
And should calculate by (86) formula;
Δ t i=N ' 2T B+ (all interruption+single-chip microcomputers (24) are interrupted+forbid in single-chip microcomputer (24) response
Output signal latchs N ' 1And N ' 2+ single-chip microcomputer (24) output signal begins sampling to first, second count section (20), (21) zero clearing+single-chip microcomputer (24) from I/O mouth B end output signal) instruct the used time (86)
T BBe the oscillation period of standard time-base generator (19).
2, ultrasonic flow meter as claimed in claim 1 is measured the method for fluid flow, it is characterized in that: the each sampling of ultrasonic flow meter begins and should be caused by the I/O mouth B end output one high level square-wave signal of the single-chip microcomputer (24) of calculating section, and the end of the each sampling of ultrasonic flow meter produces after should the full N of second singlechip asynchronous counter (14) counting by sampling section to overflow to be interrupted and finish; Sampling section sampling or do not sample, by control assembly first trigger flip-flop (16) of sampling section, second trigger flip-flop (18), first monostalbe trigger (8), second monostalbe trigger (10), first with door (11), second with door (12), the 3rd with door (15), the 4th with door (17), first or (7) and second or (9) control; The positive integer N that puts is " prolong and pass apart from technology " described N in the first single-chip microcomputer asynchronous counter (13) and the second singlechip asynchronous counter (14).
CNB200510088694XA 2005-08-03 2005-08-03 Ultrasonic flowmeter and its theory and technique Expired - Fee Related CN100434875C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB200510088694XA CN100434875C (en) 2005-08-03 2005-08-03 Ultrasonic flowmeter and its theory and technique

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB200510088694XA CN100434875C (en) 2005-08-03 2005-08-03 Ultrasonic flowmeter and its theory and technique

Publications (2)

Publication Number Publication Date
CN1908592A CN1908592A (en) 2007-02-07
CN100434875C true CN100434875C (en) 2008-11-19

Family

ID=37699766

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB200510088694XA Expired - Fee Related CN100434875C (en) 2005-08-03 2005-08-03 Ultrasonic flowmeter and its theory and technique

Country Status (1)

Country Link
CN (1) CN100434875C (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7942068B2 (en) * 2009-03-11 2011-05-17 Ge Infrastructure Sensing, Inc. Method and system for multi-path ultrasonic flow rate measurement
CN102253237B (en) * 2011-04-29 2013-02-13 浙江大学 Flow velocity measurement device and method applied to ultrasonic heat meters
JP2015232519A (en) * 2014-06-10 2015-12-24 アズビル株式会社 Clamp-on type ultrasonic flow meter and flow measurement method
CN104331624B (en) * 2014-11-11 2017-06-06 莱芜钢铁集团电子有限公司 A kind of energy flow value processing method and processing device based on flowmeter
CN104406642B (en) * 2014-11-24 2017-10-03 天津商业大学 A kind of transit time ultrasonic flow meters accurate measurement method
CN105157771B (en) * 2015-07-03 2018-04-03 中国矿业大学 A kind of transit-time ultrasonic flow measuring method and device
CN105784033B (en) * 2016-05-05 2022-02-01 汇中仪表股份有限公司 Sheet type multi-sound-path ultrasonic flow measurement sensor and use method thereof
CN106643843B (en) * 2016-12-14 2018-12-25 广州大学 The detection device and method of signal circuit delay time in ultrasonic detection device
CN106770671B (en) * 2016-12-14 2019-03-12 广州大学 Ultrasonic echo processing unit and method in a kind of ultrasonic detection device
CN115493662B (en) * 2022-11-21 2023-03-03 成都流体动力创新中心 Integrated ultrasonic flowmeter and system for aerospace

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1078552A (en) * 1992-04-01 1993-11-17 瓦尔蒂翁特克尼里伦图基穆斯克斯库斯公司 The assay method and the device of gas velocity in the pipe
US5463906A (en) * 1994-01-24 1995-11-07 Triton Technology, Inc. Interchangeable disposable acoustic for use with an ultrasonic flowmeter, particularly during extracorporeal measurement of blood flow
CN2259619Y (en) * 1995-10-13 1997-08-13 同济大学 Ultrasonic velocity difference flow meter

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1078552A (en) * 1992-04-01 1993-11-17 瓦尔蒂翁特克尼里伦图基穆斯克斯库斯公司 The assay method and the device of gas velocity in the pipe
US5463906A (en) * 1994-01-24 1995-11-07 Triton Technology, Inc. Interchangeable disposable acoustic for use with an ultrasonic flowmeter, particularly during extracorporeal measurement of blood flow
CN2259619Y (en) * 1995-10-13 1997-08-13 同济大学 Ultrasonic velocity difference flow meter

Also Published As

Publication number Publication date
CN1908592A (en) 2007-02-07

Similar Documents

Publication Publication Date Title
CN100434875C (en) Ultrasonic flowmeter and its theory and technique
Zhu et al. Variable ratio threshold and zero-crossing detection based signal processing method for ultrasonic gas flow meter
US20200056917A1 (en) Leaky Lamb Wave Flowmeter
CN106226392A (en) Water-oil phase flow containing rate measuring method based on ultrasonic attenuation mechanism model
Liang et al. Investigating the liquid film characteristics of gas–liquid swirling flow using ultrasound Doppler velocimetry
JPH09508202A (en) Ultrasonic transducer with temporary crosstalk separating means
WO2001067050A1 (en) Simultaneous determination of multiphase flowrates and concentrations
WO1993014382A1 (en) Device and method for measuring multi phase flow
Magaud et al. Experimental study of bubble column hydrodynamics
CN101614569A (en) Pipeline liquid flow-measuring method based on the supersonic guide-wave technology
US6553844B2 (en) Property-independent volumetric flowmeter and sonic velocimeter
CN102866261B (en) Method for detecting flight time of ultrasonic wave in flow speed measurement
CN102914333B (en) Detection method of using ultrasonic waves for flow detection
CN102252722B (en) Three-differential-pressure gas liquid two-phase fluid flow metering device
CN201352132Y (en) Tube cavity for ultrasonic gas meter
CN206291930U (en) A kind of ultrasonic wave mass flowmenter
Jacobson New developments in ultrasonic gas analysis and flowmetering
Coulthard Cross-correlation flow measurement—a history and the state of the art
CN202171480U (en) Three-differential-pressure gas and liquid two-phase flow metering device
Lee et al. An implementation of ultrasonic water meter using dToF measurement
CN206095327U (en) Compound multichannel flowmeter and flow metering device thereof
Yue-Zhong et al. Numerical simulating nonlinear effects of ultrasonic propagation on high-speed ultrasonic gas flow measurement
Hongguang et al. Study on the oil quantities calculation method of coriolis mass flow meter in oil dynamic measurement
Liu et al. Novel liquid flow sensor based on differential pressure method
FI76885C (en) Acoustic flow measurement method and device for application thereof.

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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20081119

Termination date: 20110803