Measurement of Fluid Flow

This invention relates to measuring the flow of fluid down axisymmetric pipes. In particular it relates to apparatus using ultrasonic transmitter/receiver pairs for measuring the average velocity of the flow of liquids and gases down pipes.

Several apparatus have been proposed in which acoustic signals may be used to produce an output indicative of flow velocity, or of parameters related to it such as mass flow rate, and some such apparatus have actually been made and sold. In at least one of these apparatus two transmitting transducers, mounted on the wall of the pipe one upstream of the other, are energized simultaneously and in synchronism from the same source. The beams that they emit travel obliquely and in opposite axial directions across the pipe in a diametrical plane to a pair of receiving transducers, also mounted on the walls of the pipe. The paths from each transmitter to its receiver are of equal length. The difference in time between the arrival of simultaneously emitted signals at their respective receivers can be shown to equal:

where L is the acoustic path length, that is to say the distance between each transmitter and its receiver, C is the velocity of the signals in still fluid, .theta. is the angle each signal beam makes with the centre line of the pipe, and V is the average velocity of the fluid in the diametrical plane in which the signals travel.

Unfortunately for most types of fluid flow within a pipe there is no simple relationship between such average velocity V in one plane and the average velocity down the tube as a whole, which is the quantity one usually wishes to known. The relationship is usually a complicated one involving factors related to the Reynolds number of the flow and to other quantities.

The present invention is the result of work to devise flow measuring apparatus to measure average flow and related quantities in ways less dependent upon such complicated factors and in particular upon the flow profile, that is to say the distribution of fluid velocity existing across a cross-section of the pipe. The character of this distribution may vary greatly according mainly as to the nature -- e.g. laminar or turbulent -- of the flow pattern of the fluid.

The invention is defined by the claims at the end of this specification and examples of it will now be described with reference to the accompanying drawings in which:

FIG. 1 is a diagrammatic plan view of the pipe of a flow-meter and shows ancillary equipment schematically;

FIG. 2 is a diagrammatic end elevation of the same flowmeter;

FIG. 3 is a diagrammatic perspective view of the same flowmeter;

FIG. 4 is a schematic diagram of part of a mass flowmeter;

FIG. 5 is a graph;

FIG. 6 is a diagrammatic section through another flowmeter, parallel to the axis and in the plane of both the acoustic paths, and

FIG. 7 is a section on the line VII -- VII in FIG. 6,

The flowmeter of FIGS. 1 to 3 comprises a cylindrical pipe 1 with axis (or centre line) c. Transmitter/receiver pairs T.sub.1, R.sub.1 and T.sub.2, R.sub.2 are mounted in the wall of the pipe. Each transmitter or receiver is mounted in a recess 3 drilled obliquely into the outside of wall 2 of pipe 1, and communicating with the bore of the pipe by a narrow passage 4 (FIG. 2) coaxial with the recess. The electrode 5 of the transmitter or receiver must of course be sealed to the recess so that the fluid within the pipe cannot leak out through the passages and recesses. Fluid flows along pipe 1 in the direction of arrow A (FIG. 1), and T.sub.1 and R.sub.1 define between them an acoustic path 6, lying oblique to the flow of fluid but in a generally downstream direction, while T.sub.2 and R.sub.2 define between themselves an acoustic path 7 that is similar but lies in a generally upstream direction. Paths 6 and 7 are of equal length, each lie at an angle .theta. to the plan plane including the axis c, and as shown in FIG. 2 lie in parallel planes 6a, 7a, one to either side of the parallel plane 8 that includes axis c, and are each separated from that plane by a distance equal to about one half of the internal radius R of the pipe.

Transmitters T.sub.1 and T.sub.2 are fed with synchronized phase-comparable signals from an oscillator 9, and the outputs of R.sub.1 and R.sub.2 pass by way of channels 10 and 11 to a signal-conditioning unit 12 and a fractional phase meter 13, and thence to a display 14. The phase difference measured by meter 13 arises of course from the fact that the signals emitted simultaneously from transmitters T.sub.1 and T.sub.2 do not arrive at the respective receivers at the same time; those from T.sub.2 have been retarded by travelling through the fluid in a direction generally against its flow, while those arriving at R.sub.1 have been advanced by travelling generally with the fluid.

In order that comparison between the time taken by the two signals to travel their respective paths should be insensitive to swirl of the fluid within the pipe (direction indicated by arrows) it is necessary that the swirl should contribute the same component, in sign, to the velocity of each signal. Equality in magnitude is achieved by locating the paths with symmetry relative to axis c, and in sign by arranging the directions of the signals so that they are both either in the same sense (clockwise or anti-clockwise) as, or in the opposite sense to the swirl about axis c. Thus FIG. 2 shows signal paths located on opposite sides of the axis, with a clockwise swirl s and with both signals having an anti-clockwise sense. In the allternative flowmeter shown in FIGS. 6 and 7 both signal paths lie on the same side of axis c, planes 6a and 7a are in effect coincident, there is an anti-clockwise swirl s and the signals along paths 6 and 7 both have an anti-clockwise sense also.

If the apparatus is to be used to give an output indicating the mass flow rate of the fluid rather than its velocity, then as shown in FIGS. 2 and 4 a pressure sensor 15 mounted from the pipe wall 2 within the measurement zone may feed a signal by way of line 16 to a circuit 17 which also receives the output of meter 13, and the output of the circuit may travel by way of line 18 to an alternative display unit 19.

Let us assume that two signals have travelled an equal distance along paths 6 and 7 respectively and are instantaneously phase-synchronized, then after each travelling a further and equal short distance along their respective paths it can be shown that their fractional phase-shift .DELTA.p will be: ##EQU1## assuming C.sup.2 is much greater than v.sup.2 Cos.sup.2 .theta.. In this equation C and .theta. are as before, R is the pipe radius, f is the signal frequency, h (see FIG. 2) is the off-centre distance of each acoustic path from the centre line of the pipe, V(r) is a function describing the flow pattern existing across and within a cross-section of the pipe, V(.sqroot. h.sup.2 + 1.sup.2) is the same function but expressed in terms of h and 1, 1 and d1 are quantities related to the acoustic path length and to an element of that length, and V is the mean flow velocity which it is desired to determine.

It has been shown by experience that for flows down circular pipes the equation: ##EQU2## generally applies, where V(c) is the velocity of the fluid along the centre line, V(y) is the fluid velocity at points situated a distance y from the pipe wall, and n is a constant of value dependent upon the Reynolds number of the flow. From equation (iii) the mean flow velocity V within the pipe may be shown to obey the equation: ##EQU3## and by combining equations (ii) to (iv), substituting for y in terms of R and r and allowing for geometrical factors the following equation may be produced: ##EQU4##

Evaluating equation (v) and plotting the resulting values of phase shift against h/R for three computed flow profiles, the graph of FIG. 5 was obtained. FIG. 5 assumes C = 344 m/s, R = 1m (unit radius), f = 40000 Hz and .theta. = 45.degree.. Two of the graphical functions represent typical turbulent flows with values for n of 6 and 10 (equivalent to Reynolds numbers of about 4 .times. 10.sup.3 and 3.2 .times. 10.sup.6) and the third function simulates a flow with uniform profile, such as would exist if n were infinite. FIG. 5 clearly shows that when planes 6a and 7a lie at h/R = 0.5, the phase shifts recorded for all three flows are very nearly equal, but that the discrepancies in phase shift for the three flows increase steadily as h/R moves substantially away from 0.5 towards the opposite extreme values of 0 and 1.0.

By plotting error (i.e. percentage deviation of observed phase shift against the true value of the quantity, determined by other methods) against Reynolds number for a series of test flows, a related result was obtained. Again it was clearly demonstrated that selection of the value h/R = 0.5 makes the readings of the apparatus for flows of different n value more consistent than is possible where h/R equals 0, and is already known, or where it has other haphazard values. However, by plotting error against Reynolds number it was also indicated that whereas a value of h/R of exactly 0.5 is exactly right when the apparatus is to used to measure only uniform and laminar flows, the value h/R = 0.52 or more exactly 0.523 may be best if the apparatus is to measure flows lying only within the turbulent range. If the apparatus must be prepared for a range of flows from laminar to near-uniform through the whole turbulent range a compromise figure such as h/R = 0.505 may be best.

A more detailed analysis of the axisymmetrical flow of fluid down a pipe of circular cross-section and unit radius involves the formulae: ##EQU5## where V(r) is the velocity of the fluid at radius r, V(c) is as before and A and B are constants. By combining these equations with equation (ii), and making substitutions similar to those made before, the equation: ##EQU6## results. Assuming B equals 2, indicating that the flow profile follows a parabolic law of which laminar flow is one special case and uniform flow may be considered another, then by integrating we obtain the equation: ##EQU7## where K is a constant. When h/R equals 0.5, it will be seen that the right-hand side of this equation simplifies to a function independent of the term A, which is itself related to velocity profile and is thus a term which this invention seeks to avoid. The equation thus simplifies to a form in which phase shift .DELTA.p is equated to mean velocity V multiplied by a constant, which is theoretically exactly what is desired.

In modifying the apparatus so that it acts as a gas mass flowmeter, and feeding in (by way of line 16, FIG. 2) a measure of the fluid pressure Pr within the measured zone, the mass flow rate m may be found by solution of the equation: ##EQU8## Here K.sup.1 appears as a constant but includes terms dependent upon the ratio .gamma. of the specific heats of the gas; .gamma. has well known values for most common gases. K.sup.1 also depends to a much smaller extent on other properties of the gas, e.g. its density. The pressure of the measured medium can easily be fed to the circuit 17 of the apparatus, and if the medium is one of the common gases the approximate or even more exact value of .gamma. can also be supplied.

If the invention is to be used to determine the flow velocity of liquid down a pipe instead of gas, account has to be taken of the fact that liquid is virtually incompressible; for an instrument to act as a mass flowmeter rather than a volumetric flowmeter, the density of the liquid must be known or determined.

Although the foregoing description and especially most of the formulae have related to fluid flow down pipes of circular cross-section, the invention includes apparatus in which measurement takes place in other axisymmetric pipes of different cross-section, but in which the planes in which the signals travel are displaced from the axis of the pipe so as to diminish the influence of velocity profile upon the readings. The acoustic transducers used may be mounted in or even outside the walls of the tube without projecting bodily into the bore and thus without physically impeding the flow at all, provided the angle of incidence of the beams to the tube walls is not so great as to give rise to total reflection, unacceptable absorption or reflection losses or the like. It should also specially be noted that the invention is based on the assumption that acoustic paths 6 and 7 are nominally of no thickness, and are in practice thin compared with the diameter of the pipe. In one piece of apparatus already tested a pipe of 100 mm. internal diameter has been used for air flow velocity measurement up to a maximum flow velocity of 10 m/s, with a potential of up to about 17 m/s. Angle .theta. has been about 45.degree. and the frequency of the ultrasound has been about 40 kHz.

The invention has been described with reference to an example in which flow (or mass flow) rate has been derived from observing the phase-difference at receivers of two alternating signals previously transmitted in phase-synchronism. However it is clearly also possible to work such a system if the signals leave the transmitter unsynchronised but with a known phase difference between them. It will also be seen that the measurement depends essentially on the difference in time taken for signals to travel two similar acoustic paths, one relatively against the flow and the other with it; this difference can be determined not only by phase-change measurements of alternating signals but also for instance by observing the difference between the transmit time of two pulsed signals. It is also possible to work the invention with one pair of transducers that can be switched, electronically for instance, so as to transmit and receive signals between them in alternate directions, the two signal paths thus being coincident in line but opposite in direction and time-comparison being made between the transit time of a signal transmitted in one direction and a subsequently signal in the other. With such an arrangement some compensation for swirl would have to be made since swirl would affect the two signals differently.