Patent Grant
10502672
This application is a national stage application filed under 35 U.S.C. I 371, of International Patent Application No. PCT/EP2015/073806, filed Oct. 14, 2015, which claims the benefit of and priority to European Patent Application No. 14188853.7, filed Oct. 14, 2014, the contents of each of which are incorporated by reference herein in its entirety.
This invention related to the determination of physical or thermodynamic properties of fluids.
It has particular application to determining properties which are measurable in respect of a first fluid when disposed within a surrounding fluid in the form of a droplet or bubble. Examples of such properties include surface tension, Bond number (and similar dimensionless characteristic numbers such as the Eötvös, Goucher, and Deryagin numbers), and density difference between phases, as wells as characteristics of the droplet or bubble including its volume, height, maximum width, base width, height of centre of gravity, apex curvature radius, interface area, base contact angle and capillary length. Since certain of these properties result from the effect of external fields and forces on a fluid (e.g. gravitational and other accelerating fields, electromagnetic fields etc., the method has application also in measuring such fields and forces.
The measurement of the properties of a fluid, such as its surface tension in a surrounding liquid or gas, is needed in many applications and across many industries including the food, textile, chemical, oil, pharmaceutical, biological and electronic industries, and in research institutions working in these fields. Similar to the surface tension measurement, other geometric properties of a bubble or drop can be measured using the same methodology, such as the bubble/drop volume, interface area, apex curvature radius, or height of centre of gravity. The angle at the base of the bubble/drop (known as contact angle) can also be measured similarly. This latter measurement is particularly useful in measuring the wettability property between a solid surface and a fluid; something that is vital in many engineering applications.
A common method of estimating surface tension involves accurately measuring and analysing the shape of a drop (or bubble) of the fluid. This shape is defined mathematically by a complex differential equation that involves the surface tension property as well as the respective fluid densities and gravity, which can be measured in other ways or looked up from a table of known values. By regression fitting a numerically integrated mathematical solution of the equation to the experimentally measured shape of the drop, the value of the surface tension is calculated.
From a physical point of view, the shape of such an axisymmetric bubble is a result of the hydrostatic pressure gradients in both the liquid and gas phases and of the capillary equilibrium at the liquid-gas interface. The mathematical equation of the bubble profile can be written as follows (see F. J. Lesage, J. S. Cotton and A. J. Robinson, Analysis of quasi-static vapour bubble shape during growth and departure, Physics of Fluids, Vol. 25, p. 067103, 2013):
In equation 1, ρl, ρg and a are properties of the fluid, respectively the liquid and gas densities and the interfacial surface tension; g is the gravitational acceleration; R0 is the radius of curvature at the apex of the bubble; z is the vertical coordinate, downward from the apex of the bubble, and C is the curvature of the interface, which depends on the vertical coordinate (thus expressed as C(z)). The profile of the bubble is fully defined by this equation, and is cut by a horizontal plane which is the horizontal surface on which the bubble is attached.
More generally, the terms ρl, ρg can be replaced by the densities of any two fluids and are not necessarily those of a liquid and a gas. Equation 1 can therefore be understood as covering the more general case of two fluids, not necessarily a liquid and a gas.
The interfacial surface tension, σ, is a property of the two fluids at the interface and therefore will be different for e.g. a water droplet in a heavy oil medium than it will for water in a gaseous medium. While interfacial surface tension is strictly speaking thus a property of the liquid and the gas, all gases effectively behave the same way, and so in the specific case of a liquid/gas interaction the surface tension is generally considered as a property specific to the liquid.
This equation is based on the assumption that the bubble or droplet is not freely floating but rather is attached to a surface. (Note that droplets and bubbles are, for these purposes, different examples of the same phenomenon, namely a discrete body of a first fluid disposed in a second fluid.) The droplet/bubble can be gravitationally accelerated, due to weight and buoyancy forces arising from the density difference with the surrounding fluid, either towards the surface in question, in which case it is a “sessile” droplet or bubble; or away from the surface, in which case it is a “pendant” droplet or bubble, remaining attached by surface tension forces that oppose the gravitationally-induced buoyancy or weight.
The bubble shown in
Equation 1 therefore provides a relation that links different properties of the two fluids and the gravitational acceleration to the geometry of the droplet/bubble. Thus, if sufficient independent parameters are known or measured, the other parameters can be deduced using this equation.
The usual way of calculating the surface tension of a fluid is as follows. The gravitational acceleration and the fluid densities are considered as known (or measured with a different method). Then, different geometrical profiles corresponding to equation 1 are generated by numerical integration of equation 1. Each numerical integration provides a solution which describes the entire surface of the bubble; from the base to the tip. A regression algorithm is utilized to obtain the solution that best fits the experimental profile, which is obtained by image processing. The best fitting profile provides the value of the surface tension property.
Similarly, the volume of the drop or bubble, or the height/position of its centre of gravity can be calculated when a best fit geometric profile is obtained describing the curvature of the profile.
Such methods have the drawback that they require significant image processing capability and mathematical computational power, as well as a good deal of post-processing time to obtain a good fit between the mathematical model and the experimental image.
There is provided a method of determining the properties of a fluid body in the form of a surface-attached droplet/bubble, comprising the steps of:
The number of linear dimensions is specified as two or more. A third linear dimension may be required, depending on the dynamics of the system. For example, a bubble/drop issuing from an orifice of fixed dimension, which does not spread beyond the orifice, may have an implicit base diameter which can be built into the model or data set. Similarly a highly wetting fluid on a surface may result in a triple contact line of fixed dimension, reducing the number of required linear dimensions to specify a particular bubble/drop volume and shape to two.
It has been found that a set of data can be generated which characterises a family or universe of droplet/bubble shapes, where each such shape is defined using a combination of at least two linear measurements. Those linear measurements uniquely specify a shape from among the range of possible shapes, and the shape in turn is (based on a consideration of equation 1) attributable to the properties of the two fluids involved and the gravitational force. Accordingly, the experimental measurement of the required two linear dimensions of a particular bubble allows the shape, and the physical parameters associated with that shape, to be pinpointed in the data set.
This means that the determination of a physical property such as surface tension, can be obtained by modelling a universe of droplets/bubbles, where each is characterised by a pair of linear measurements, and then when it is desired to find the surface tension of a fluid one can simply make a couple of accurate linear measurements and look up the surface tension.
There are significant advantages with this method. Measurement of the linear dimensions of a droplet/bubble can be performed accurately without sophisticated image processing software. This makes this method suitable for implementation in a device which may not have sophisticated imaging or processing power, such as with the camera on a smartphone. Because the method involves a couple of linear measurements and a look-up operation, it can be performed in real time, unlike a curve-fitting exercise as in known methods. Linear dimensions can be measured with a higher degree of accuracy compared with regression fitting solutions of a numerical integration to a complex shape. Consequently, a higher accuracy may be reached using this method to measure surface tension.
It is not limited to surface tension. Once a droplet/bubble's shape has been modelled, it is possible to store, for that bubble, not only surface tension but also several other parameters which are specified once the shape is determined, e.g. volume, height of centre of gravity, contact angle, interface area, or capillary length. Accordingly the measurement of a pair of linear dimensions of a bubble and the provision of a data set in which bubbles are categorised and identified according to those measurements, allows the rapid and accurate determination of a number of complex physical parameters which were previously obtainable only with difficulty.
Preferably, one of the first and second fluids is known and the other is unknown, and the properties of the known fluid permit the derivation from said one or more parameters of corresponding properties of the unknown fluid.
This will often be the case in an experimental or real-world industrial set-up. The unknown fluid will be provided as a droplet/bubble in a fluid whose properties are completely known, such as water or air; or alternatively, the unknown fluid can be studied by introducing into it a bubble or droplet of e.g. air, water or a water-immiscible liquid such as toluene.
Preferably, said one or more parameters comprise a parameter which is a function of an accelerating field, a surface tension of one fluid at the interface with the other fluid, and the respective fluid densities.
Preferably, the accelerating field is the local gravitational field as characterised by the acceleration due to gravity, g. It can be envisaged however that applications will arise in which there is no local gravitational field, or in which the gravitational acceleration is supplemented by another external accelerating field. Such accelerations will influence the shape of any droplet or bubble. Similarly other external fields and forces can be taken into account including electromagnetic forces acting on the fluid.
Preferably, the two or more linear dimensional measurements are normalised measurements.
Further, preferably, the two or more linear measurements are normalised against a further linear measurement of the bubble/drop.
This is particularly useful as a method of standardising the data set. It also allows a dimension which is kept constant in the system to be used as a normalising dimension.
Preferably, the further linear measurement is a base diameter or base radius of the bubble drop at a surface to which it is attached.
Preferably, the two or more linear dimensional measurements of the data set comprise any two of the following measurements normalised against the remaining measurement: height normal to attachment surface, maximum width parallel to attachment surface, and base diameter (or radius) at attachment surface.
Optionally, the two or more linear dimensions of the data set may be expressed as a combination of dimensions such as an area or a volume.
Preferably, the data set is limited based on one or more of the following assumptions used to create the data set:
In this way, the data set can be simpler, i.e. the number of droplet/bubble shapes can be greatly compressed if assumptions are made, such as (most commonly) the value of the gravitational acceleration, but also possibly the data set may be tailored to identifying e.g. the surface tension of an unknown liquid by looking at the shape of an air bubble within that liquid, in which case the density of the air is known. Other simplifying assumptions can be made in other scenarios to reduce the complexity of the data set. Since the surface tension may be unique to the liquid, this provides a potential method of identifying a liquid.
Preferably, said set of data comprises a plurality of parameter sets, each parameter set describing a unique solution to an equation modelling the shape of a droplet/bubble, and each parameter set including said combination of two or more linear dimensional measurements and said one or more parameters describing the relationship between the physical properties of a pair of fluids capable of providing said solution.
There is also provided a method of obtaining the interfacial surface tension between a liquid in a gas, comprising the steps of:
There is also provided a computer program product comprising a set of instructions, which are effective to cause a processor to:
There is also provided an apparatus for determining the properties of a fluid body in the form of a surface-attached droplet/bubble, comprising:
Preferably, the apparatus further comprises a measurement system for making said plurality of linear dimensional measurements, and an output therefrom to said processor
The measurement system may comprise an optical sensor and imaging software calibrated to determine linear measurements of key parameters of a droplet/bubble whose image is captured by the optical sensor.
More preferably, the measurement system comprises a laser source and a detector for determining interception of a laser beam by an edge of a droplet.
Preferably, the apparatus further comprises means for introducing said droplet/bubble into said surrounding medium.
The shape of a bubble changes as the bubble grows. Equation 1 defines, for each pair of liquids in a given gravitational field (thus with, ρl, ρg, g and σ fixed) a unique curvature or capillary profile C(z) where the actual shape of the bubble or drop is determined by the height at which this profile is cut by the surface to which it is attached.
To illustrate this point,
The curvature of the profile in
In the present innovative method, a large number of solutions to equation 1 are numerically computed a priori, and the solutions are stored in a database. These solutions preferably cover the whole spectrum of bubble shapes that may be encountered. The solutions are numerically treated, in order to extract key geometrical features. In particular, key lengths such as bubble base radius, bubble height and bubble maximum width are calculated. Within certain conditions, a given set of these three linear lengths is sufficient to identify a unique bubble shape solution.
Using the bubble of the case of study, the three lengths can be measured very easily and nearly instantaneously. Importantly, it can be done very accurately as it does not rely on the use of sophisticated image processing, mathematical development or regression analysis (see
In this case, the mentioned equation is solved for the case of a pendant drop or bubble i.e. one which is suspended from an orifice. The key information about the bubble can be interpolated from the database using these three simple length measurements, two if the base dimension is physically fixed. In order to illustrate the database, its information is represented in
The representation of the database of solutions in
In
It can be seen that the isolines of bubble width are close to horizontal over much of their ranges, meaning that a given measured width will be associated with the same normalized volume for a fairly wide range of Bond numbers. However, for liquids with higher Bond numbers there is a tendency for the same width to be associated with bubbles of greater volume.
It can be seen from
For the experimental bubble the following linear measurements were made for width, height, and base diameter (see
The normalized dimensions vs. base radius are thus:
These two normalized parameters are sufficient to locate the bubble at a specific x-y location on the graph, as illustrated in
Where d is the bubble base diameter, ρl and ρg are respectively the liquid and gas (air) densities, g is the gravitational acceleration and σ is the surface tension. In the experimental case, g is taken to be 9.81 m·s−2, and the density of air, ρg, is assumed to be 1 kg·m−3. The liquid density, ρl, can be found by looking up the value for a known liquid, or by weighing a known volume. In this case it is found to be 1000 kg·m−3.
The surface tension of the liquid can then be determined from these values in a straightforward manner:
Thus, using the graph (or the underlying database) and the three simple linear measurements, together with well-known or easily found values for the densities of fluids and gravitational acceleration, the surface tension can be found.
It will be understood that the representation in the graphs of
It will also be understood that blank areas as seen in
Locating the bubble shape at a specific x-y location on the graph allows not only the Bond number to be derived but also each of the other plotted values which have been pre-calculated for each solution of the equation in the database. Those parameters are each plotted on a separate graph in
Thus one may solve a very large number of solutions of the governing equation, and store relevant geometric features of these solutions in a database. Then, keeping the example of the surface tension measurement, it is only necessary to measure simple dimensions of the drop (or bubble) such as its height, maximum width and base width, in order to interpolate from the data in the database the value of its surface tension. This innovative method avoids the necessity to fit a mathematical solution to the entire drop contour by solving all possible solutions of the equation a priori. It results in a simpler, faster measurement, with the possibly of much improved accuracy.
The same method can be used in order to measure the bubble (or drop) volume, the height of its centre of gravity, its interface area, the radius of curvature at its apex or any geometrical feature, including the apparent contact angle at the base of the bubble or drop. This is particularly useful to characterize the wettability property of a fluid on a given surface. In the case of contact angle, it is more appropriate to use a sessile bubble (or drop) on (or below) a horizontal flat surface. The governing equation that would describe such a bubble/drop would still be equation 1, but with a negative value for the gravitational acceleration.
The database and graphs of
For example,
Thus, by measuring the values for w, h and d as previously described, and by calculating the quantities h/w and 2 w/d, a simple lookup will give the isoline value that overlies the identified point. In the case of the experimental bubble discussed earlier, the y-axis value is 1.252 and the x-axis value is again 4.95. The Isoline value at this point is −1.28, which is log10 (0.052). This simple lookup allows a very rapid calculation of Bond number from three measurements that can be made easily and with great accuracy.
The accuracy of the graphical method was compared against the leading existing techniques and against theory as follows. A gas bubble is created in water from a 1.22 mm diameter orifice located in a horizontal, upward-facing surface immersed in quiescent water in terrestrial gravity conditions. (This is basically the bubble pictured in
The volume of gas introduced into the bubble was measured at 17.78 mm3, and the capillary length of water has a known value of 2.7 mm. Consequently, the bubble base radius is b=0.61 mm, the Bond number value is Bo=0.052 and the normalized volume is Vb*=78.3.
In order to find the “true” volume for the purposes of comparing different approaches, the following technique was used. The volume was calculated by imaging the bubble with a high resolution photo, using the pixel locations on the interface and then revolving to create layers of 3D disks (with bevelled edges) and then adding the subsequent volumes to get the total bubble volume (Vb=17.78 mm3). The normalized volume is Vb*=Vb/b3.
As a first set of data for comparison, the dimensions of the bubble (normalized height, normalized width, normalized height of centre of gravity, and apparent contact angle) can be measured from the photographic image of the bubble.
As a second set of data, a normalized solution to the Laplace equation is calculated based on the known parameters to obtain a model of the resultant bubble, which is illustrated in
As a third set of data, we have the data from the database underlying the graph of
Table 1 sets out the comparative numerical values. It can be seen that a very good agreement is found between the three different methods as all values present less than 3% discrepancy.
The method can be implemented in computer code running on a computing device. An example of such a device is illustrated in block diagram form in
The functionality need not all be carried on-board a single device. The image of a bubble or droplet can equally be captured by an external imaging system such as a digital camera or CCD disposed as part of a measurement apparatus, with a connection to the computing device.
The storage means 52 also has a database 54 within which data, representing multiple solutions to the equation describing the shape of a bubble or droplet, are stored. Such data typically is organised in a manner that allows a look-up of a solution based on a limited number of dimensional input parameters, such as a combination of two or three linear measurements of a bubble or droplet (or a composite parameter like a measurement of one dimension of the bubble, normalized with respect to another measurement). The data stored for each solution can be one or more complex parameters of the shape of bubble described by the solution, including without limitation the Bond number (or a similar characteristic of the system), surface tension, contact angle, volume, area, height of centre of gravity, radius of curvature, or parameters of the curvature itself. Thus, the database allows the look-up of complex parameters from simple linear measurements.
Also stored within the storage unit 52 is program code 56 which specifies the steps to be taken in extracting the features from an image, and deriving from the database the desired complex parameters. A network interface, such as is known in conventional mobile phones and tablets, allows communication from and to the device.
The program code is preferably in the form of an “app” or downloadable program, in the case of the device 40 being a mobile phone or tablet.
Image processing in the software then measures the dimensions of the image, step 62, by performing a conventional edge detection to determine the edges of the droplet or bubble, straightening the image according to a reference axis (such as the axis of symmetry of the bubble or a reference edge visible in the image), and then measuring the length of predetermined dimensions, such as bubble height, maximum width, and width at the connection point to a surface (or base diameter). The measured dimensions can be determined in terms of numbers of pixels across each dimension, or converted to actual estimated lengths e.g. in mm, according to depth cues or assumed distance from the bubble.
Alternative, and indeed improved, methods can equally be used to capture linear dimensions, and the technique encompasses such methods. A limitation in the use of a camera or CCD is that accuracy is restrained to the pixel resolution of the image.
In one improvement, the camera is omitted and instead dimensions are captured using laser interferometry. For example, if the drop height, width and base diameter are measured by displacing a laser beam, and measuring when the beam is intercepted by the drop, then two advantages arise for an enhanced accuracy:
In step 64, the dimensions are normalized by the processor. In the preferred method, the base diameter or radius is used to normalize the width and height. However, this need not be the case, and any of these dimensions can be normalized against any of the others. For example, in the
In step 66, the program accesses the database 54 to look up solution values which correspond to the dimensions as measured and normalized. Depending on the granularity of the database, the resolution of the image, and the degree of rounding when normalizing, the exact values used to look up the database index may or may not be present. Accordingly, in decision 68, a determination is made if the values are present for a direct look up.
If so, step 70, then the pre-calculated values are determined for the solution that corresponds to the normalized linear dimensions. If not, a suitable interpolation routine is employed in step 72 to obtain values that are deemed to be a best match.
The retrieved or interpolated parameters can include items such as the Bond number, from which the surface tension can be calculated when the densities of the liquid and gas are known (and assuming standard gravity or some adjusted accelerating field). Post-processing steps may therefore allow the user to input additional values and obtain derived properties like surface tension, wettability, or any other property that can be derived or calculated from the values retrieved from the database.
The server 82 has a network interface 84, a processor 86, working memory 88, and the aforementioned non-volatile storage medium 90. It also carries program instructions 94 relating to its part of the functionality of the method in this client-server relationship.
In one implementation, the program instructions at the device are limited to capturing the image and sending the image to the server computer over the network, with the server computer's program causing the server to analyse the image, obtain the dimensional values, look up the database 92, and return as an output whatever values are sought by the device 40 for display to the user.
In another implementation, the device carries out the image processing and the extraction of the dimensions, and simply transmits these to the server for look-up of the stored complex parameters of an appropriate solution in the database. The skilled person will appreciate that the distributed nature of computing solutions allows many other variations in the configuration of a suitable computer system, enabling measurements of a droplet to be converted to useful complex parameters of the droplet, based on the common factor of having a pre-calculated family of solutions describing different candidate droplets to which the dimensional measurements can be fitted to select an appropriate solution.
| Number | Date | Country | Kind |
|---|---|---|---|
| 14188853 | Oct 2014 | EP | regional |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/EP2015/073806 | 10/14/2015 | WO | 00 |
| Publishing Document | Publishing Date | Country | Kind |
|---|---|---|---|
| WO2016/059124 | 4/21/2016 | WO | A |
| Number | Name | Date | Kind |
|---|---|---|---|
| 4416148 | Klus et al. | Nov 1983 | A |
| 6438261 | Moshe | Aug 2002 | B1 |
| 6577700 | Fan | Jun 2003 | B1 |
| 6873725 | Xu | Mar 2005 | B2 |
| 8254657 | Pollack | Aug 2012 | B2 |
| Number | Date | Country |
|---|---|---|
| 2009036634 | Feb 2009 | JP |
| Entry |
|---|
| Lesage et al., Analysis of quasi-static vapour bubble shape during growth and departure. Physics of Fluids. 2013;25:067103. |
| Lesage et al., Experimental and numerical analysis of quasi-static bubble size and shape characteristics at detachment. International Journal of Heat and Mass Transfer. 2013;64:53-69. |
| Number | Date | Country | |
|---|---|---|---|
| 20170227435 A1 | Aug 2017 | US |
