METHOD AND SYSTEM FOR MEASURING THE KINEMATIC VISCOSITY OF A FREE FLUID STREAM

Abstract

A computer implemented method for measuring the viscosity of a fluid stream flowing from an opening with constant acceleration, wherein the method takes as input at least one image of an area of interest of a fluid stream, the section, s.sub.0, of the opening and the output volume flow rate, U.sub.0, of the fluid stream from said opening; wherein the method provides as output the viscosity, , of the fluid; wherein the method includes (a) modelling the geometrical profile of the fluid stream trough an digital processing operation of the input image, and (b) computing the viscosity of the fluid with a mathematical or physical model from the modelled geometrical profile, the section, s.sub.0, and the output volume flow rate, U.sub.0.

Claims

1. A computer implemented method for measuring the kinematic viscosity of a free fluid stream flowing from an opening with constant acceleration, wherein said method takes as input at least one image of an area of interest of a fluid stream, a section, s.sub.0, of said opening and an output volume flow rate, U.sub.0, of the fluid stream from said opening; wherein said method provides as output a kinematic viscosity, , of the fluid; wherein said method comprises the following steps: (a) modelling the geometrical profile GP of the fluid stream through a digital processing operation of the input image; (b) computing the kinematic viscosity of the fluid with a mathematical or physical model from the modelled geometrical profile GP, the section, s.sub.0, and the output volume flow rate, U.sub.0.

2. A computer implemented method, wherein said method takes as input a first image of an area of interest of a fluid stream, a section, s.sub.0, of said opening and an output volume flow rate, U.sub.0, of the fluid stream from said opening, a second time-shifted image of at least two features of the fluid stream, wherein said features are present in the first image provided as input, wherein said method provides as output a kinematic viscosity, , of the fluid; wherein said method comprises the following steps: digital processing of the first and second images to measure a displacement of said features; computing velocities of each said features from a measure of the displacement over a time shift between the first and second images, and computing the kinematic viscosity of the fluid with a mathematical or physical model from the velocities, the section, s.sub.0, and the output volume flow rate, U.sub.0.

3. The computer implemented method according to claim 1, wherein said method further takes as input a second time-shifted image of at least one feature of the fluid stream, wherein said feature is present in the first image provided as input, and wherein the output volume flow rate, U.sub.0, of the fluid stream from said opening is computed after the modelling step with following steps: digital processing of the first and second images to measure the displacement of said feature; computing a velocity of the feature from a measure of the displacement over the time shift between the first and second images; computing the output flow rate, U.sub.0, by multiplying the computed velocity by a section, s1, of a modelled geometrical profile GP of the fluid stream at a location of the feature in the first image.

4. The computer implemented method according to claim 3, wherein the output volume flow rate, U.sub.0, is computed with a plurality of features of the fluid stream so that to compute a mean output volume flow rate, U.sub.mean.

5. The computer implemented method according to claim 3, wherein the features are a defect within the fluid stream.

6. The computer implemented method according to claim 1, wherein a mathematical or physical model is experimental and/or simulated charts of relationships between geometrical profiles for a fluid stream, its output volume flow rate, its viscosity, and the section of the opening.

7. The computer implemented method according to claim 1, wherein a section of the opening and of the fluid stream are circular.

8. A data processing device comprising means for carrying out a method according to claim 1.

9. A computer program comprising instructions which, when the program is executed by a computer, causes the computer to carry out a method according to claim 1.

10. A non-transitory computer-readable medium comprising instructions which, when the instructions are executed by a computer, cause the computer to carry out a method according to claim 1.

11. A process for measuring a kinematic viscosity, , of a fluid stream flowing from an opening with constant acceleration, wherein said process comprises the following steps: (a) acquiring, with an image recording device, at least one image of an area of interest of a fluid stream; (b) modelling, with a data processing device, a geometrical profile of the fluid stream trough a digital processing operation of the acquired image; (c) computing, with a data processing device, the kinematic viscosity, , of the fluid with a mathematical or physical model from the modelled geometrical profile, the section, s.sub.0, of the opening and the output volume flow rate, U.sub.0, of the fluid stream.

12. A system for measuring a kinematic viscosity, , of a fluid stream flowing from an opening with constant acceleration, wherein said system comprises: a image recording device configured to acquire at least one image of an area of interest of a fluid stream; a data processing device according to claim 8 and configured to receive images from the image recording device.

13. A method comprising performing the process according to claim 11 in a manufacturing line of glass fibers to monitor the kinematic viscosity of a molten glass flowing through a bushing by gravity.

14. The method according to claim 13, wherein the monitored kinematic viscosity is implemented into a feedback operation for adjusting the temperature of the molten glass in the vicinity of the bushing.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0027] FIG. 1 is a schematic representation of an example of manufacturing line of glass fibres.

[0028] FIG. 2 is a schematic cross section representation of example of a bushing and a fiberizing tool from a detail I of the [FIG. 1].

[0029] FIG. 3 is a schematic representation of a detail II of the [FIG. 2].

[0030] FIG. 4 is a data flow diagram of a computer implemented method according to one embodiment of the first aspect of the invention.

[0031] FIG. 5 is a data flow diagram of a computer implemented method according to second embodiment of the first aspect of the invention.

[0032] FIG. 6 is a data flow diagram of a computer implemented method according to an embodiment of the first aspect of the invention.

[0033] FIG. 7 is a plot of geometrical profiles of a fluid stream for an output volume flow rate U.sub.0, 6.3 m.sup.3/day, at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0034] FIG. 8 is a plot of geometrical profiles of a fluid stream for an output volume flow rate U.sub.0, 8.4 m.sup.3/day, at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0035] FIG. 9 is a plot of geometrical profiles of a fluid stream for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0036] FIG. 10 is a plot of geometrical profiles of a fluid stream for an output volume flow rate U.sub.0, 6.3 m.sup.3/day, at a section so of 8 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0037] FIG. 11 is a plot of geometrical profiles of a fluid stream for an output volume flow rate U.sub.0, 8.4 m.sup.3/day, at a section so of 8 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0038] FIG. 12 is a plot of geometrical profiles of a fluid stream for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 8 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0039] FIG. 13 is a plot of velocity flow profiles of a fluid stream for an output volume flow rate U.sub.0, 6.3 m.sup.3/day, at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0040] FIG. 14 is a plot of velocity flow profiles of a fluid stream for an output volume flow rate U.sub.0, 8.4 m.sup.3/day, at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0041] FIG. 15 is a plot of velocity flow profiles of a fluid stream for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0042] FIG. 16 is a plot of velocity flow profiles of a fluid stream for an output volume flow rate U.sub.0, 6.3 m.sup.3/day, at a section so of 8 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0043] FIG. 17 is a plot of velocity flow profiles of a fluid stream for an output volume flow rate U.sub.0, 8.4 m.sup.3/day, at a section so of 8 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0044] FIG. 18 is a plot of velocity flow profiles of fluid stream for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 8 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0045] FIG. 19 is a comparison plots of velocity flow profiles of a fluid stream for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 5.7 cm.sup.2 with example experimental data.

[0046] FIG. 20 is a comparison plots of velocity flow profiles of a fluid stream for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 5.7 cm.sup.2 with example experimental data.

[0047] FIG. 21 is a physical data flow diagram of a processing data system to implement a method according to the first and/or second aspect of the disclosure

[0048] FIG. 22 is a schematic representation of a system for measuring the kinematic viscosity of a fluid stream flowing from an opening with constant acceleration.

[0049] FIG. 23 is a plot of geometrical profiles of a fluid stream for an output volume flow rate U.sub.0, 6.3 m.sup.3/day, at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614.

[0050] FIG. 24 is a plot of geometrical profiles of a fluid stream for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 8.0 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614.

[0051] FIG. 25 is a plot of velocity flow profiles of fluid stream for an output volume flow rate U.sub.0, 6.3 m.sup.3/day, at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

[0052] FIG. 26 is a plot of velocity flow profiles of fluid stream for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 8 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes.

DETAILED DESCRIPTION OF EMBODIMENTS

[0053] With reference to [FIG. 1], a manufacturing line 1000 of glass fibres may generally comprise silos 1001 for storing raw materials 1001a, e.g., minerals compounds and/or cullet, a glass melting furnace 1002 to melt raw materials 1001a which are conveyed from silos 1001 by means of conveyor 1003, and one or more fiberizing tools 1005a, 1005b, 1005c which are fed with molten glass 1006 through a forehearth 1004. The forehearth is usually designed as an open or closed channel provided with openings right above each fiberizing tools 1005a, 1005b, 1005c to fed each of them with molten glass 1006.

[0054] With reference to [FIG. 2], an opening 2001 of a forehearth 1004 may be provided with a bushing 2002 which is here exemplified as an opening 2002a extended by a nozzle or drain pan 2002b and terminated by another opening 2002c from which the molten glass 1006 flows as a fluid stream 2003 by gravity to a fiberizing tool 1005c located beneath. The fiberizing tool 1005c may comprise a spinner 2004 made of bucket 2003a with side openings and of an annular wall 2004b provided with a plurality of calibrated holes. Upon high-speed rotation, the molten glass 1006 coming out of the bushing 2002 as fluid stream 2003 and falling into the bucket 2004a is projected out from its side openings toward the inner surface of the annular wall 2004b, and then projected though the calibrated holes to cast glass threads 2005.

[0055] The fiberizing tool 1005c may also comprise a ring burner 2006 that throws out an elevated temperature gaseous stream or jet in a substantially tangent direction to the annular wall to pull down the casted glass threads 2005 in order to heat and strength them and form glass fibres 2007. It may also comprise a gas blowing device 2008 located below the ring burner 2005 to prevent the glass fibres 2007 from departing too far away from the rotation axis of the spinner 2004.

[0056] As exemplified on [FIG. 2] and [FIG. 3], when the molten glass 1006 is flowing from the opening 2002c of the bushing 2002 by gravity, it is flowing freely as a free fluid stream 2003 with constant acceleration.

[0057] With reference to [FIGS. 3] and 4, in a first aspect of the disclosure, there is provided a computer implemented method 4000 for measuring the kinematic viscosity of a free fluid stream 2003 flowing from an opening 2002c with constant acceleration, wherein said method 4000 takes as input at least one image I4001 of an area A3001 of interest of a fluid stream 2003, the section, s.sub.0, of said opening 2002c and the output volume flow rate, U.sub.0, of the fluid stream 2003 from said opening 2002c; [0058] wherein said method 4000 provides as output the kinematic viscosity, , of the fluid; [0059] wherein said method 4000 comprises the following steps: [0060] (a) modelling S4001 the geometrical profile GP of the fluid stream 2003 through a digital processing operation of the input image I4001; [0061] (b) computing S4002 the kinematic viscosity of the fluid with a mathematical or physical model from the modelled geometrical profile GP, the section, s.sub.0, and the output volume flow rate, U.sub.0.

[0062] In the context of the disclosure, a free fluid stream is to be understand as a free-flowing fluid, i.e., as a flowing fluid with free surfaces which are not bounded by external surfaces such as those from a tube, a pipe, or the like.

[0063] In the context of the disclosure, the geometric profile GP should be understood as any dimensional measurement which is representative or allows to derive the variation of the geometric dimensions, e.g., width, section, diameter . . . of the fluid stream 2003 in the flow direction.

[0064] In an example embodiment, as illustrated by the reference X, Z axis on [FIG. 3], in a case of a fluid stream with a circular section, the geometric profile may be the variation of the section, the diameter or the radius in the X direction along the flow direction Z. The origin of the X, Z axis may be the opening 2002c.

[0065] It may occur that, because of, for instance, the radiation of the fluid stream 2003, e.g., a fluid stream of elevated temperature molten mineral glass, or the industrial surrounding environments of the opening 2002c, that the at least one input image I4001 is marred or spoiled by too much noise which does not allow an accurate modelling of the geometrical profile GP.

[0066] Instead of measuring modelling the geometrical profile GP, it may be advantageous to model the flow velocity of the fluid steam 2003, i.e., the velocities of the fluid stream at separate locations of the fluid stream along the flow direction, for instance along the Z direction on [FIG. 3]. This may help to overcome the limitation of a too noisy input image I4001.

[0067] Accordingly, in alternative advantageous embodiments, with reference to [FIG. 3], [FIG. 4] and [FIG. 5], the method 5000 may take as input a first image I4001 of an area A3001 of interest of a fluid stream 2003, the section, s.sub.0, of said opening 2002c and the output volume flow rate, U.sub.0, of the fluid stream 2003 from said opening 2002c, a second time-shifted image I5001 of at least two features F3001, F3002 of the fluid stream 2003, wherein said features F3001, F3002 is present in the first image I4001 provided as input, [0068] wherein said method 5000 provides as output the kinematic viscosity, , of the fluid; [0069] wherein said method 5000 comprises the following steps: [0070] digital processing S5001 of the first and second images I4001, I5001 to measure the displacement of said features F3001, F3002; [0071] computing S5002 the velocities, V1, V2, of each said features F3001, F3002 from the measure of the displacement over the time shift between the first and second images I4001, I5001; and [0072] computing S4002 the kinematic viscosity of the fluid with a mathematical or physical model from the velocities, V1, V2, the section, s0, and the output volume flow rate, U.sub.0.

[0073] The velocities, V1, V2 of the features F3001, F3002, may be viewed as two different measures of the flow velocity, FV, of the fluid stream. In the computing step S4002, each of these two velocities, V1, V2, may be used independently to compute the value of the kinematic velocity. Such approach, in particular when several velocities are computed from several features within the fluid stream, may advantageously provide better accuracy in the measurement of viscosity.

[0074] Advantageously, the more the velocities are computed for the flow velocity FV, the better may be the accuracy of the method. Accordingly, the velocities of a plurality of features, e.g., at least 10 features, preferably at least 50 features, more preferably at least 100 features may be computed. The computation may be performed with much more input images that the two input images I4001, I5001.

[0075] The features F3001, F3002 or the plurality of features of the fluid stream 2003 may be any features of the fluid stream 2003 which moves with the fluid stream 2003 and shows certain persistence between the first and second images I4001, I5001 so that their displacement over time may be measured, and thus, their velocity.

[0076] In certain embodiments, the at least two features F3001, F3002 may be a defect within the fluid stream 2003. For instance, for molten glass, the defects may be bubbles, non-melted particles, e.g., refractory stones or non-melted raw materials, or other heterogeneities. The defects may be at the surface of the fluid stream or within the bulk. For transparent molten glass, defects within the bulk may be easily detected thanks to the transparency.

[0077] The fluid stream 2003 may sometimes show a homogeneous quality and be devoid of visible defect. Artificial defects may be introduced upstream, e.g., with bubbling device in the forehearth for molten glass, or solid particles.

[0078] The section, s.sub.0, of the opening 2002c may be either provided as an independent value from the technical specifications of the opening 2002c, or may be measured directly through an image processing of the first image I4001 of the area A3001.

[0079] The output volume flow rate, U.sub.0, provided as input may be measured independently from the method according to the first aspect of the invention. For instance, it may be measured automatically with dedicated flowmeters located just beneath the opening 2002c, e.g., with method and an apparatus as described in WO 8304437 A1 [GULLFIBER AB [SE]] 22.12.1983 provided that the section of the fluid stream 3001 may be known or measured. It may also be measured manually from timed sampling.

[0080] Alternatively, or complementarily, the method according to the first aspect of the disclosure may also be adapted to measure the output volume flow rate, U.sub.0. Accordingly, in one embodiment, with reference to [FIG. 3], [FIG. 4] and [FIG. 6], the method 6000 may further take as input a second time-shifted image I6001 of at least one feature F3001 of the fluid stream 2003, wherein the said feature F3001 is present in the first image I4001 provided as input, and wherein the output volume flow rate, U.sub.0, of the fluid stream 2003 from said opening 2002c is computed after the modelling step S4001 with following steps: [0081] digital processing S6001 of the first and second images I4001, I6001 to measure the displacement of said feature F3001; [0082] computing S6002 the velocity, V, of the feature F3002 from the measure of the displacement over the time shift between the first and second images I4001, I6001; [0083] computing S6003 the output flow rate, U.sub.0, by multiplying the computed velocity, V, by the section, s1, section of the modelled geometrical profile GP of the fluid stream 3001 at the location of the feature F3001 in the first image I6001.

[0084] As an illustration, referring to [FIG. 3], the section, s0, of the opening (2002c), is fixed. However, the section, s, of the fluid stream 2003, varies in the Z direction depending on the location of the feature F3001 within the fluid stream 2003 in said Z direction. Thus, the section, s1, of the modelled geometrical profile of the fluid stream 2003 will vary accordingly.

[0085] In certain embodiments, to improve the accuracy of the measurement, the output volume flow rate, U.sub.0, may be computed with a plurality of features of the fluid stream 2003 so that to compute a mean output volume flow rate, Umean. The computation may be performed with much more input images that the two input images I4001, I6001.

[0086] It may be noteworthy that, according to the law of conservation of mass, the output mass flow rate should not vary in the falling direction, i.e., Z direction, of the fluid stream 2003 in stationary conditions. Therefore, in stationary conditions, which occur in most cases, provided that the temperature, and in turn the density, does not vary or at least at a very limited and negligible extend, the computed output volume flow rate, Umean, is representative measure of the true output volume flow rate, U.sub.0, of the fluid stream 2003 from the opening 2002c.

[0087] The feature F3001 may be of the same nature as those discussed above in the context of the alternative embodiments of the [FIG. 5]. Further, when the alternative embodiments of the [FIG. 5] are implemented, steps S6001 and S6002 may be advantageously replaced by steps S5001 and S5002 as they also allow to compute velocities of features of the fluid stream. The input image I6001 may be replaced by the input image I5001. This may limit redundant steps and save time.

[0088] The area A3001 of interest may be any area of the fluid stream 2003. Because of the constant acceleration, the fluid stream 2003 may become thinner with distance from the opining 2002c. Thus, it may be advantageous that the fluid stream 2003 shows a certain width for the digital processing operation of steps S4001, S5001, S6001 to work efficiently on the input image I4001, I5001, I6001. Image of thinnest areas of the fluid stream 2003, e.g., the lowest area of the fluid stream, should then be avoided. Practically, the area A3001 of interest may extend from the opening 2002c to a length twice, preferably three times, the diameter of the opening 2002c. The width of the fluid stream in the input image I4001, I5001, I6001 may then be large enough for the digital processing operation.

[0089] In the modelling step S4001, the geometric profile GP may be computed with any adapted edge detection algorithm for image processing, e.g., marching square algorithm, canny edge detector, thresholding, edge operators based on Prewitt, Sobel or Scharr filters.

[0090] In computing steps S4002, according to certain embodiments, the mathematical or physical model may be experimental and/or simulated relationships, e.g., charts, between geometrical profiles for a fluid stream, its output volume flow rate, its viscosity, and the section of the opening.

[0091] As example embodiments, the charts may be derived from earlier experiments in which a substitute fluid, e.g., cold oil, the kinematic viscosity of which may be easily changed is made flow as stream with constant acceleration through openings with different sections at different volume flow rates. The kinematic viscosity may be changed by varying the temperature of the fluid or its composition, e.g., by diluting it in water. Images acquired from several trials conducted with different values for the kinematic viscosity, the opening section and the volume flow rate may allow to construct charts in which a kinematic viscosity value corresponds to a given geometric profile, a given volume flow rate and a given opening section. The data of the charts may be experimental data or model derived from experimental data.

[0092] Example charts are provided on [FIG. 7] to 18 for different values of section so of the opening 2002c, different values of the output volume flow rate U.sub.0, and different values of kinematic viscosity.

[0093] [FIG. 7] to [FIG. 12] show the geometric profiles of a circular fluid stream represented as an evolution of its radius R in the flow direction Z. [FIG. 7] to [FIG. 9] represent this evolution for three values of output volume flow rate U.sub.0, 6.3 m.sup.3/day, 8.4 m.sup.3/day and 9.6 m.sup.3/day at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes (cm.sup.2/s). [FIG. 10] to [FIG. 12] represent this evolution for three values of output volume flow rate U.sub.0, 6.3 m.sup.3/day, 8.4 m.sup.3/day and 9.6 m.sup.3/day, at a section so of 8.0 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes (cm.sup.2/s).

[0094] [FIG. 13] to [FIG. 18] show the evolution of the flow velocity FV of a fluid stream in the flow direction Z. [FIG. 13] to 15 represent this evolution for three values of output volume flow rate U.sub.0, 6.3 m.sup.3/day, 8.4 m.sup.3/day and 9.6 m.sup.3/day at a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes (cm.sup.2/s). [FIG. 16] to [FIG. 18] represent this evolution for three values of output volume flow rate U.sub.0, 6.3 m.sup.3/day, 8.4 m.sup.3/day and 9.6 m.sup.3/day, at a section so of 8.0 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes (cm.sup.2/s).

[0095] In [FIG. 7] to [FIG. 18], the origin for the flow direction Z is the outlet of the opening 2002c.

[0096] The values of section so of the opening 2002c, the values of the output volume flow rate U.sub.0, the values of kinematic viscosity, the geometrical profiles and the flow velocity profiles provided on [FIG. 7] to [FIG. 18] are all representative of industrial actual industrial conditions, in particular of a fluid stream of molten glass for the manufacturing glass fibres for glass or stone wool. A person skilled in the art may then directly use these charts for figuring out the kinematic viscosity of fluid stream of molten glass falling in a fiberizing tool. If required, intermediate values may be extrapolated from the data of the [FIG. 7] to [FIG. 18].

[0097] As examples, [FIG. 19] and [FIG. 20] provides respectively a comparison of different velocity flow, FV, and geometric profiles, i.e., radius R of a fluid stream of molten glass with experimental data acquired in industrial conditions for an output volume flow rate U.sub.0, 9.6 m.sup.3/day, at a section so of 5.7 cm.sup.2. Both [FIG. 10] and [FIG. 20] show the almost perfect agreement between the modelled profiles and the experimental data, and demonstrate that the charts according to the disclosure allow an accurate and precise measure of the kinetic viscosity in a fluid stream under industrial environments or conditions.

[0098] As other alternative or complementary example embodiments, the mathematical or physical model may be a numerical resolution of the Navier-Stokes equation applied to a freely falling stream of fluid for different conditions, e.g., for different values of the kinematic viscosity, the opening section, and the volume flow rate. Relationships between kinematic viscosity, the opening section, and the volume flow rate may be simulated as geometrical profiles or flow velocity profile and, for instance, represented as charts in which a given kinematic viscosity value corresponds to a given geometric profile or flow velocity profile, a given volume flow rate and a given opening section. Alternatively, the Navier-Stokes equations may be used as a mathematic function to fit, or model, a geometric profile or flow velocity profile of the fluid stream.

[0099] As an illustrative example, referring to [FIG. 2] and [FIG. 3], the Navier-Stokes equation applied to a freely falling stream 2003 of fluid from an opening 2002c may be written as follows:

[00001]$\frac{\mathrm{dFV}}{\mathrm{dZ}}-3v\frac{d}{\mathrm{dZ}}\left(\frac{1}{\mathrm{FV}}\frac{\mathrm{dFV}}{\mathrm{dZ}}\right)=\frac{g}{\mathrm{FV}}$

[0100] Where v is the kinetic viscosity, so is the section of the opening (2002c), FV is the velocity of the fluid stream (2003), g is the acceleration constant and z is the coordinate in the flow direction Z from the opening 2003, i.e., the vertical direction of the opening.

[0101] To numerically solve the above equation, two boundary conditions may be used:

[00002]$\begin{array}{cc}\mathrm{FV}=\frac{U_0}{s_0}& \left(1\right)\end{array}$

[0102] (2) the viscous stress at the end of fluid stream may be neglected, i.e.,

[00003]$v\frac{\mathrm{dFV}}{\mathrm{dZ}}=0;$

[0103] The equation may be solved numerically for different values of the kinetic viscosity, v, of the section, s.sub.0, of the opening 2002c, and of the output volume flow rate, U.sub.0, using a solver for solving boundary value problems, e.g., the bpv_solve solver from the scipy python package.

[0104] [FIG. 23] to 24 show the evolution of the flow velocity FV of a fluid stream in the flow direction Z as it may be obtained from such numerical resolution. [FIG. 23] represents this evolution for the output volume flow rate, U.sub.0, 6.3 m.sup.3/day with a section s.sub.0 of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes (cm.sup.2/s). [FIG. 24] represents this evolution for three values of output volume flow rate, U.sub.0, 9.6 m.sup.3/day, with a section so of 8.0 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes (cm.sup.2/s).

[0105] [FIG. 25] to 26 show the same geometric profiles of a circular fluid stream represented as an evolution of its radius R in the flow direction Z, as it may be obtained from the same numerical resolution. [FIG. 25] represents this evolution for three values of output volume flow rate U.sub.0, 6.3 m.sup.3/day, with a section so of 4.5 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes (cm.sup.2/s). [FIG. 26] represents this evolution for three values of output volume flow rate, U.sub.0, 9.6 m.sup.3/day, with a section so of 8.0 cm.sup.2 and for values of kinematic viscosity varying from 87 to 614 Stokes (cm.sup.2/s).

[0106] In [FIG. 23] to [FIG. 26], the origin for the flow direction Z is the outlet of the opening 2002c.

[0107] In the context of the disclosure, the section of the opening 2002c may have any geometric form which allows a fluid stream 2003 to flow through it. In certain embodiments, the section of the opening 2002c and of the fluid stream 2003 may be circular, as this geometric configuration corresponds to cases in current industrial lines. The section of the opening 2002c may then be computed from its diameter.

[0108] In preferred embodiments, the fluid stream 2003 is a stream of molten mineral glass.

[0109] In a second aspect of the disclosure, with reference to [FIG. 21], there is provided a data processing system 21000 comprising means 21001 for carrying out a method 4000, 5000, 6000 according to any one of the embodiments of the first aspect of the invention, and a computer program 121001 comprising instructions which, when executed by a computer, cause the computer to carry out a method according to any one of embodiments of the first aspect of the invention.

[0110] The data processing system 21000 comprises means 21001 for carrying out a method according to any of the embodiments of the first aspect of the invention. Example of means 21001 may be a device which can be instructed to carry out sequences of arithmetic or logical operations automatically to perform tasks or actions. Such device, also called computer, may comprise one or more Central Processing Unit (CPU) and at least a controller device that are adapted to perform those operations.

[0111] It may further comprise other electronic components like input/output interfaces 21003, non-volatile or volatile storage devices 21002, and buses that are communication systems for the data transfer between components inside a computer, or between computers. One of the input/output devices may be user interface for human-machine interaction, for example graphical user interface to display human understandable information.

[0112] As calculation may require a lot of computational power to process substantial amounts of data, the data processing system may advantageously comprise one or more Graphical Processing Units (GPU) whose parallel structure makes them more efficient than CPU, in particular for image processing.

[0113] The computer program 121001 may be written through any kind of programming language, either compiled or interpreted, to implement the steps of the method according to any embodiments of the first aspect of the invention. The computer program 121001 may be part of a software solution, i.e., part of a collection of executable instructions, code, scripts, or the like and/or databases.

[0114] In certain embodiments, there may also be provided a computer-readable storage or medium 21002 comprising instructions which, when executed by a computer, cause the computer to carry out the method according to any of the embodiments of the first aspect of the invention.

[0115] The computer-readable storage 21002 may be preferably a non-volatile non-transitory storage or memory, for example hard disk drive or solid-state drive. The computer-readable storage may be removable storage media or a non-removable storage media as part of a computer.

[0116] Alternatively, the computer-readable storage may be a volatile memory inside a removable media.

[0117] The computer-readable storage 21002 may be part of a computer used as a server from which executable instructions can be downloaded and, when they are executed by a computer, cause the computer to carry out a method according to any of the embodiments described herein.

[0118] Alternatively, the program 121001 may be implemented in a distributed computing environment, e.g., cloud computing. The instructions may be executed on the server to which client computers may connect and provide encoded data as inputs to the method. Once data are processed, the output may be downloaded and decoded onto the client computer or directly send, for example, as instructions. This kind of implementation may be advantageous as it can be realised in a distributed computing environment such as a cloud computing solution.

[0119] In a third aspect of the disclosure, with reference to [FIG. 3], [FIG. 21] and [FIG. 22], there is provided a process for measuring the kinematic viscosity, , of a fluid stream 2003 flowing from an 2002c opening with constant acceleration, wherein said process comprises the following steps: [0120] (a) acquiring, with an image recording device 22001, at least one image of an area A3001 of interest of a fluid stream 3001; [0121] (b) modelling, with a data processing mean 21000, the geometrical profile of the fluid stream trough a digital processing operation of the acquired image; [0122] (c) computing, with a data processing mean 21000, the kinematic viscosity, , of the fluid with a mathematical or physical model from the modelled geometrical profile, the section, s.sub.0, of the opening and the output volume flow rate, U.sub.0, of the fluid stream 3001.

[0123] All embodiments described forth in the context of the first aspect of the invention may apply to the process according to the third aspect of the invention. More precisely, all the embodiments on the method 4000, 5000 may be adapted in the process, in particular regarding steps (b) and (c) of said process.

[0124] In a fourth aspect of the disclosure, with reference to [FIG. 3], [FIG. 21] and [FIG. 22], there is provided a system 22000 for measuring the kinematic viscosity, , of a fluid stream 2003 flowing from an opening 2002c with constant acceleration, wherein said system comprises: [0125] a image recording device 22001 configured to acquire at least one image of an area A3001 of interest of a fluid stream 2003; [0126] a data processing device 21000 according to the second aspect of the disclosure and configured to receive images from the image recording device 22001.

[0127] In certain embodiments, the image recording device 22001 may be a digital camera, e.g., CCD or CMOS digital camera. The image resolution of the image recording device may depend on the desired precision for measuring the kinematic viscosity. Generally, the smaller the width of the stream fluid 2003, the greater the image resolution.

[0128] In advantageous embodiments, the image resolution of digital camera is so that the width resolution of acquired image of the area A3001 of interest of a fluid stream 2003 is at least 200 pixels/cm, preferably at least 400 pixels/cm. These embodiments may suit most of prerequisites of current manufacturing lines.

[0129] In preferred embodiments, the digital camera may further be a high-speed digital camera, e.g., a high-speed digital camera with a frame rate of at least 50 frames per second (fps), preferably at least 100 fps. High frame rates may be advantageous to acquire sequences of images which may thereafter be fed to the data processing device 21000. The data processing device 21000 may then be further configured to compute mean or average image of the area A3001 of interest of the fluid stream 2003 from the image of said sequence. The mean image may then be uses to model the geometric profile GP with better accuracy and precision.

[0130] High frame rates may also be advantageous to compute the output volume flow rate, U.sub.0, from images of a feature or a plurality of features as described forth in the context of the first aspect of the invention.

[0131] As illustrated on [FIG. 22], the image recording device 22001 may be oriented perpendicular to the flow direction of the fluid stream 2003. In certain embodiments, the image recording device 22001 may be oriented or tilted with a certain angle from the flow direction when there is not enough space at proximity of the opening 2002c and the fluid stream for an image recording device 22001 to be placed perpendicular to the flow direction, as it may happen in industrial environment. The tilt angle may be considered when processing the acquired images in order to correct the induced optical deformation within the images.

[0132] The distance at which the image recording device 22001 may be placed from the fluid stream 2003 may depend on the focal length, e.g., 50 cm to 1 m, of the device. In certain embodiments, depending on the temperature of the fluid stream 2003, and/or its heat radiation, the image recording device may be placed at a higher distance from the fluid stream to preserve its electronics from heat. A heat shield may also be placed around the camera.

[0133] The process, the data processing device and the system according to the second, third and fourth aspects of the disclosure may be advantageously used in a manufacturing line of glass fibres to monitor the kinematic viscosity of a molten glass 1006 flowing through a bushing 2002 by gravity. The manufacturing line may comprise any kind of fiberizing tool, e.g., fiberizing tools comprising a spinner as described earlier or comprising a bottom closed spinner.

[0134] In further embodiments, the monitored viscosity may be implemented into a feedback operation for adjusting the temperature of the molten glass 1006 in the vicinity of the bushing 2002.

[0135] In this context, in certain embodiments, the system 22000 may further comprise a controller device configured to set or change one or several parameters of one or several components of the manufacturing line which may have an action on the temperature of the molten glass. In some example embodiments, when the kinematic viscosity departs from a set point value, the controller device may act on cooling and/or heating devices related to a forehearth 1004 so that to increase or decrease the temperature of the molten glass before a bushing 2002. In other example embodiments, the controller device may send a visual signal onto a display device to alert a human operator to adjust the chemistry of the molten glass.

[0136] All embodiments described herein, whether it concerns the first, second, third or fourth aspect of the invention, may be combined by one skilled in the art unless they appear to him technically incompatible.

[0137] Further, although the invention has been described in connection with preferred embodiments, it should be understood that various modifications, additions, and alterations may be made to the invention by one skilled in the art without departing from the spirit and scope of the invention as defined in claims.