Method for determining a measured viscosity value, and coriolis mass flowmeter for carrying out the method

Abstract

A method (100) serves for determining a measured viscosity value of a medium by means of a Coriolis mass flowmeter with a measuring tube and comprises: exciting (110) a measuring-tube vibration; ascertaining (122) a measured eigenfrequency value f of the measuring-tube vibration; ascertaining (124) a measured damping value Dm; ascertaining (126) a measured flow-rate value ṁ of the medium; ascertaining (128) a measured density value ρ of the medium; and ascertaining (130) a measured viscosity value (130) ηm as a function of a viscosity-relevant damping coefficient Dη, on the basis of the above measured values; wherein the said damping coefficient is obtained from the measured damping value Dm by subtraction of an eigen damping coefficient D0 of the measuring tube and a flow-dependent damping coefficient Df, which can be represented by a model function MD which is proportional to a product of a power of a Reynolds number and a power of a Stokes number, (formula (I)), wherein the Reynolds number Re and the Stokes number St each include a provisional viscosity value, wherein the following applies for the exponents: 1.8 > ERe > 1.3; 4.0 > Est > 3.4 and 2.8 > Est / ERe > 2.2, wherein the measured viscosity value nm is ascertained with an algorithm which ascertains the viscosity-relevant damping coefficient Dη by the flow-dependent damping coefficient Df being ascertained in relation to the measured damping value Dm by means of the model function MD or by means of a routine which approximates the behaviour of the model function MD and wherein the measured viscosity value is determined as a function of the viscosity-relevant damping coefficient Dη.

Description

[0001] Method for determining a viscosity measurement value and Coriolis mass flow meter for carrying out the method

[0002] The present invention relates to a method for determining a viscosity measurement value of a medium by means of a Coriolis mass flowmeter and to a Coriolis mass flowmeter for carrying out the method.

[0003] A generic determination of a viscosity measurement value is usually carried out by measuring the damping of a measuring tube oscillation of a Coriolis mass flowmeter.

[0004] Publication WO 2021 037 492 A1 discloses that viscosity measurement exhibits cross-sensitivity to the flow rate or flow velocity, since a flow-dependent damping contribution occurs, particularly at high flow rates, without which the viscosity is overestimated. However, an inaccurate viscosity measurement leads to an incorrect Reynolds number, which has a correspondingly detrimental effect on a Reynolds number correction of the flow measurement value.

[0005] The determination of a correction factor for the flow rate or density as a function of the Reynolds number is the subject of patent DE 10 2021 202 464 B3. This patent teaches a saturated curve of the correction factor with a tanh function of the logarithm of the Reynolds number around the point of the steepest slope of the correction factor as a function of the logarithm of the Reynolds number. While this is a pragmatic approach, it does not adequately account for the diversity of factors influencing damping and flow measurement. Therefore, the object of the present invention is to remedy this situation.

[0006] The object is achieved according to the invention by the method according to independent patent claim 1 and the Coriolis mass flowmeter according to independent patent claim 13.

[0007] The method according to the invention is used to determine a viscosity measurement value of a medium by means of a Coriolis mass flowmeter with at least one measuring tube in which the medium is guided, wherein the measuring tube is to be excited by an exciter to oscillate in a bending vibration useful mode, wherein the oscillation experiences a damping D which depends on the viscosity and the density of the medium, as well as the mass flow rate of the medium and the natural frequency of the bending vibration useful mode of the measuring tube, wherein the method comprises the following steps: exciting an oscillation of the measuring tube in a bending vibration useful mode; determining a natural frequency measurement value f of the bending vibration useful mode; determining a damping measurement value Dm of the bending vibration useful mode; determining a flow rate measurement value m of the medium; determining a density measurement value p of the medium;and determining a viscosity measurement value r|m for the medium, based on the damping measurement value Dm, the density measurement value p, the natural frequency measurement value f and the flow rate measurement value m, wherein the viscosity measurement value r|m is a function of a viscosity-relevant damping contribution D; n which is determined from the damping measured value Dm by subtracting a self-damping contribution Do of the measuring tube and a flow-dependent damping contribution Df, i.e.:

[0008] Dr) = Dm _ do _ D f , where the flow-dependent damping contribution Df can be represented by a model function MD which is proportional to a product of a power of a preliminary Reynolds number and a power of a preliminary Stokes number

[0009] D f = A f ■ Re ERe ■ St Est, where a provisional viscosity value is included in the provisional Reynolds number Re and the provisional Stokes number St, where the following applies to the exponent of the Reynolds number Re: ER6> 1.3; where the following applies to the exponent of the Reynolds number Re: ER6< 1.8; where the following applies to the exponents of the stroke number St: Est > 3.4; where the following applies to the exponents of the stroke number St: Est < 4.0; where the quotient of the exponents Est / ERe > 2.2; where the quotient of the exponents Est / ERe < 2.8; where the viscosity measurement value r|m is determined using an algorithm which calculates the viscosity-relevant damping contribution D n determined by determining and subtracting the flow-dependent damping contribution Df from the damping measured value Dm using the model function MD or by means of a routine that approximates the behavior of the model function MD, and where the viscosity measured value is determined as a function of the viscosity-relevant damping contribution D nis determined, whereby in particular in the case of approximation the routine approximates the behavior of the model function MD in such a way that a viscosity-relevant damping contribution D determined using the approximated flow-dependent damping contribution Df v for actual viscosities of not less than 10 cP, the viscosity-relevant damping contribution D determined using the model function does not deviate from the latter by more than 3% of a flow-dependent damping contribution Df.

[0010] In a further development of the invention, the following applies to the exponent of the Reynolds number Re: ER6> 1.4; where the following applies to the exponent of the Reynolds number Re: ER6< 1.7; where the following applies to the exponents of the stroke number St: Est > 3.5; where the following applies to the exponents of the stroke number St: Est < 3.9; where the quotient of the exponents Est / ERe > 2.3; and where the quotient of the exponents Est / ERe < 2.65.

[0011] In a further development of the invention, a provisional viscosity value is included in the provisional Reynolds number Re and the provisional Stokes number St; the provisional Reynolds number is given as: where the preliminary Stokes number is given as:

[0012] St = i . pE , d; j 2-rrf where di is the inner diameter of the measuring tube, where n r the number of

[0013] measuring tubes, and where nv is the preliminary viscosity value.

[0014] In a further development of the invention, the viscosity measurement value is determined on the basis of a

[0015] Damping measurement value is determined iteratively in N iteration steps, whereby to determine the viscosity measurement value in the j+1-th iteration, the viscosity measurement value after the j-th iteration is used as a preliminary viscosity value in the calculation of the preliminary Reynolds number and the preliminary Stokes number, where 0 < j < N.

[0016] In a further development of the invention, a first preliminary viscosity value is determined on the basis of the damping measurement value Dm under the assumption that the flow-dependent damping contribution Df is zero, or a first preliminary viscosity value is a literature value.

[0017] In a further development of the invention, the natural frequency measured value f is not more than 600 Hz, for example not more than 400 Hz and in particular not more than 200 Hz.

[0018] In a further development of the invention, a first preliminary viscosity value rjo is to be estimated with a first quadratic estimate, according to which the first preliminary viscosity value rjo is given as proportional to the square of the difference between the measured damping value Dm and the self-damping contribution Do according to: where E(f,p) is a device-specific prefactor depending on the natural frequency measured value and the density, where G is a device-specific constant, where a viscosity measurement value r|m is to be estimated with a second quadratic estimate, according to which the viscosity measurement value r|m is given as proportional to the square of the viscosity-relevant damping contribution D T1 according to: wherein the algorithm performs or approximates the first quadratic estimate and the second quadratic estimate, wherein in particular in the case of approximating the estimates the approximated viscosity measurement value does not deviate from the latter by more than the maximum of 1 cP and 5% of the value according to the estimates, in particular not more than 2% of the value according to the estimates.

[0019] In a further development of the invention, the algorithm approximates the behavior of the model function MD by estimating the flow-dependent damping contribution Df as where Af is a device-specific constant.

[0020] In one embodiment of this further development of the invention, a viscosity output value r|A is determined which is the minimum of the first preliminary viscosity value o and the viscosity measurement value m: A = min (0, r|m)> wherein, in particular in the case of the approximation of the estimates, the viscosity output value nA does not deviate from the latter by more than the maximum of 1 cP and 5% of the viscosity output value nA according to the estimates, in particular by not more than 2% of the viscosity output value nA according to the estimates.

[0021] In a further development of the invention, the viscosity measurement value nm is a function of the viscosity-relevant damping contribution D T1 with an estimate iteratively in N iterations according to: hj+i = E(f, p) ■ D Vij where r|j+i is the viscosity value after the j+1-th iteration, where E(f, p) is a device-specific prefactor dependent on the natural frequency measured value and the density, where G is a device-specific constant, where D ,j is the viscosity-relevant damping contribution after the j-th iteration, which can be determined as where Df is the flow-dependent damping contribution after the j-th iteration, which can be determined as

[0022] D fJ = A r Rej ERe ■ Stj Est where the Reynolds number Rei and the Stokes number Sti can be determined after the j-th iteration as where j is the index of the iteration steps with 0 < j < N, where N > 1 , in particular N>2, where in particular the iteration is stopped if the viscosity measurement value changes by no more than 2% from one iteration step to the next; where the algorithm carries out or approximates the estimation, where in particular in the case of approximation the approximated viscosity measurement value does not deviate from the value of the estimate by more than the maximum of 1 cP and 5% of the value of the estimate.

[0023] In a further development of the invention, a Reynolds number and / or a Stokes number is determined on the basis of the determined viscosity measurement value. In a further development of the invention, the viscosity output value and / or the viscosity measurement value for actual viscosities between 1 cP and 300 cP and flow velocities between 0 m / s and 5 m / s in a feed line to the Coriolis mass flowmeter with its nominal diameter deviates by no more than a maximum of 10% of the actual viscosity and 1 cP above or below the actual viscosity, and / or the viscosity output value and / or the viscosity measurement value for actual viscosities between 1 cP and 300 cP and flow velocities between 5 m / s and 10 m / s in a feed line to the Coriolis mass flowmeter with its nominal diameter deviates by no more than a maximum of 25% of the actual viscosity and 2.5 cP above or below the actual viscosity.

[0024] In a further development of the invention, the viscosity output value and / or the viscosity measured value for actual viscosities between 300 cP and 1000 cP and flow velocities between 0 m / s and 3 m / s in a feed line to the Coriolis mass flowmeter with its nominal diameter deviates by no more than the maximum of 10% of the actual viscosity and 1 cP from the actual viscosity upwards or downwards, and / or the viscosity output value and / or the viscosity measured value deviates for actual viscosities between 300 cP and 1000 cP and flow velocities between 3 m / s and 6 m / s in a feed line to the Coriolis mass flowmeter with its nominal diameter by no more than the maximum of 25% of the actual viscosity and 2.5 cP from the actual viscosity upwards or downwards.

[0025] In a further development of the invention, the damping measurement value is determined based on a ratio between an oscillation amplitude of the measuring tube oscillation in the bending oscillation useful mode and the amplitude of the associated excitation current, or based on a width of a resonance curve Af / f of the bending oscillation useful mode, or based on a decay time of the oscillation in the bending oscillation useful mode with the excitation power switched off, reduced, or increased. The Coriolis mass flowmeter according to the invention is equipped with: at least one measuring tube for conveying a flowable medium; an exciter for exciting an oscillation of the measuring tube in a bending oscillation useful mode; two oscillation sensors for detecting the oscillation of the measuring tube in a bending oscillation useful mode; and a measuring and operating circuit configured to carry out the method according to one of the preceding claims.

[0026] The invention will now be explained in more detail with reference to the exemplary embodiments illustrated in the drawings. It shows:

[0027] Fig. 1 : a schematic representation of an embodiment of a Coriolis mass flow meter according to the invention;

[0028] Fig. 2a: a schematic diagram for uncorrected viscosity measurements and viscosity measurements corrected according to the invention as a function of the flow rate for a high viscosity and a low Stokes number;

[0029] Fig. 2b: a schematic diagram for uncorrected viscosity measurements and viscosity measurements corrected according to the invention as a function of the flow rate for a high viscosity and a high Stokes number;

[0030] Fig. 2c: a schematic diagram for uncorrected viscosity measurements and viscosity measurements corrected according to the invention as a function of the flow rate for a low viscosity and a low Stokes number;

[0031] Fig. 2d: a schematic diagram for uncorrected viscosity measurements and viscosity measurements corrected according to the invention as a function of the flow rate for a low viscosity and a high Stokes number;

[0032] Fig. 3a: a diagram for a sequence of viscosity values ​​for a first

[0033] Embodiment of the method according to the invention; and

[0034] Fig. 3b: a diagram for a sequence of viscosity values ​​for a second

[0035] Embodiment of the method according to the invention; Fig. 4a: a flowchart for an embodiment of the method according to the invention;

[0036] Fig. 4b: a detailed flowchart with substeps of a

[0037] Method step of the embodiment of the method according to the invention from Fig. 4a; and

[0038] Fig. 4c: a detailed flowchart with substeps of another

[0039] Method step of the embodiment of the method according to the invention from Fig. 4a.

[0040] The exemplary embodiment of a Coriolis mass flowmeter 1 according to the invention shown in Fig. 1 comprises an oscillator 10 which comprises a pair of parallel, oscillatable measuring tubes 14 which extend between an inlet-side flange 11 and an outlet-side flange 12, wherein the flanges each comprise a flow divider or collector into which the measuring tubes 14 open. The flow dividers are connected to one another by a rigid housing 15, so that vibrations of the flow dividers accommodating the measuring tubes in the range of vibration frequencies of useful bending vibration modes of the oscillator are effectively suppressed. The measuring tubes 10 are rigidly connected to an inlet-side node plate 20 and an outlet-side node plate 21, wherein the node plates define vibration nodes of the oscillator 10 formed by the two measuring tubes 14, and thus largely determine the frequencies of the useful bending vibration modes.The oscillator 10 is excited to oscillate by an electrodynamic exciter 17 acting between the two measuring tubes 14, wherein the oscillations are detected by means of two vibration sensors 18, 19 that detect relative movements of the measuring tubes 14. The exciter 17 is operated by a measuring and operating circuit 30, which also detects and evaluates the signals from the vibration sensors in order to determine a natural frequency of the useful bending vibration mode, a density measurement value, a mass flow measurement value, and a damping measurement value. The measuring and operating circuit 30 is also configured according to the invention to carry out the method according to the invention, as explained further below. Unlike shown in Fig. 1, the measuring and operating circuit 30 can also comprise several spatially separated modules.The calculation of the damping and viscosity measurements can also be carried out in a remote computing unit to which the required raw data is transmitted, for example wirelessly.

[0041] Furthermore, the Coriolis mass flowmeter according to the invention can also have a different number of measuring tubes, for example, four measuring tubes, each of which forms an oscillator in pairs, or a single measuring tube. Suitable Coriolis mass flowmeters for implementing the invention are offered by the applicant, for example, under the designations Promass F, Promass X, Promass Q, Promass E, Promass A, Promass S, and Promass I.

[0042] The schematic diagrams shown in Figs. 2a to 2d show uncorrected viscosity measurements (solid lines) and viscosity measurements corrected according to the invention (dash-dotted lines) as a function of flow velocity in a supply line to a Coriolis mass flowmeter according to the invention, wherein the supply line has a diameter that corresponds to the nominal diameter of the Coriolis mass flowmeter. The uncorrected viscosity measurements are uncorrected in that the damping of the measuring tube oscillation, which increases with increasing flow velocity, is interpreted entirely as being caused by the viscosity of the medium, thus neglecting a flow-dependent damping contribution Df as the cause of the increasing damping.The present invention now makes it possible to identify this flow-dependent damping contribution Df and to subtract it from a determined damping measurement value Dm in addition to the self-damping contribution Do, i.e. the expression.

[0043] Dr) = Dm _ do _ D f (I), and based on the viscosity-relevant damping contribution D thus obtained n to determine a correct viscosity measurement.

[0044] Investigations in connection with the present invention have shown that the flow-dependent damping contribution Df can be represented as proportional to a product of a power of a preliminary Reynolds number and a power of a preliminary Stokes number, i.e.:

[0045] Df = A f ■ Re ERe ■ St Est (II).

[0046] The preliminary Reynolds number Re is given as:

[0047] Re = 4|m| (III), and

[0048] H-di-n R T|v the preliminary Stokes number is given as:

[0049] The Stokes number and the Reynolds number are provisional in that they themselves depend on the viscosity of the medium. Therefore, plausible preliminary values ​​for the viscosity are required to obtain a reasonable preliminary Reynolds number Re and a reasonable preliminary Stokes number St. This can be done, for example, with an estimate in which the flow-dependent damping contribution is assumed to be zero. The following applies to the exponent of the Reynolds number Re: 1.8 > ER6 > 1.3. The following applies to the exponent of the Strokes number St: 4.0 > Est > 3.4. The quotient of the exponents is 2.8 > Est / ERe > 2.2.

[0050] For a given Coriolis mass flowmeter, the flow-dependent damping contribution Df can be determined experimentally, for example, by recording damping values ​​Dm for some media as a function of flow rate and calculating the differences to the damping at zero flow. Based on these differences, the parameters of the function from equation (II) can be determined by regression calculation:

[0051] If the pre-factor Af is determined for a specific measuring device and, if necessary, a fine-tuning of the exponents has been carried out within the limits given above, a basis is given to determine the flow-dependent damping contribution during ongoing measuring operation, with the help of which the viscosity-specific damping contribution can then be determined based on a damping measured value, which ultimately forms the basis for determining the viscosity measured value.

[0052] The effectiveness of the method according to the invention is evident from the diagrams in Figs. 2a-2d. The flow velocity range shown in the diagrams extends from 0 m / s to approximately 10 m / s. The Reynolds number in the measuring tube(s) of the Coriolis mass flowmeter increases proportionally to the flow velocity in the supply line.

[0053] The actual viscosity of the medium for the curves in Figs. 2a and 2b is a few 10 mPas. The actual viscosity of the medium for the curves in Figs. 2c and 2b is only a few mPas. Accordingly, the flow-dependent damping contribution in the diagrams in Figs. 2a and 2b increases even at lower velocities compared to the situation in Figs. 2c and 2d. This initially appears counterintuitive, since the viscosity is in the denominator of the above defining equation for the velocity-proportional Reynolds number. However, the viscosity is in the numerator of the defining equation for the Stokes number, so that because the exponent of the Stokes number is more than a factor of 2 larger, the flow-dependent damping contribution Df increases with viscosity. The curves in Figs.Figures 2a and 2c are characteristic of a Coriolis mass flowmeter with stiffer measuring tubes, thus higher natural frequencies, compared to the measuring tubes of the Coriolis mass flowmeter, for which the curves in Figures 2b and 2d are characteristic. Therefore, the curves in Figures 2b and 2d represent a situation with a larger Stokes number, which leads to larger uncorrected viscosity values ​​than in Figures 2a and 2c.

[0054] After describing above how the individual damping contributions are to be determined, the following explains how a viscosity value can be determined from a damping contribution.

[0055] The relationship between viscosity-specific damping D T1 and viscosity r| can be modeled, for example, as: rj = E(f, p) ■ D. 2 (VI), where E(f, p) is a device-specific prefactor dependent on the natural frequency measured value and the density, where G is a device-specific constant.

[0056] By applying equation (VI), a first preliminary value rjo for the viscosity can be determined, which must be inserted into equations (III) and (IV) to obtain a preliminary Reynolds number and a preliminary Stokes number.

[0057] Neglecting the flow effect, the first preliminary value for the viscosity is to be determined according to:

[0058] Subsequently, a flow-dependent damping contribution Df and viscosity-relevant damping contribution D n by means of equations (I) and (II), the latter allowing the determination of a viscosity measurement value r|m according to equation VI.

[0059] Either the viscosity measurement value or the minimum of the first viscosity measurement value and the preliminary viscosity measurement value can be output as the viscosity output value r|a, i.e.: ha = min (r|0, r|m)

[0060] The results of this calculation lead to the curves of the corrected viscosity measurements in the diagrams of Figs. 2a - 2d.

[0061] In a variant of the invention, the algorithm approximates the behavior of the model function MD by estimating the flow-dependent damping contribution Df as Egt 2 where Af is a device-specific constant.

[0062] The accuracy of determining viscosity measurements depends on the flow rate and the actual viscosity.

[0063] The determination of the viscosity measurement value has a certain iterative character, as explained below with reference to Figs. 3a and 3b.

[0064] In the example described above, the preliminary viscosity value rjo was calculated based on the measured damping value Dm, assuming a stationary medium. Thus, the total damping, less the inherent damping Do of the measuring tube, was assigned to the viscosity. For a flowing medium, this not surprisingly leads to an excessively high preliminary viscosity value r|o, which is shown in Fig. 3a in relation to the actual viscosity value iqt. With this slightly excessive preliminary value, which is included in the preliminary Reynolds number and the preliminary Stokes number, the flow-dependent damping contribution according to equation (II) is naturally also somewhat too high, which can lead to the viscosity-relevant damping contribution being determined to be somewhat too low. As a result, the measured viscosity value would therefore be too low.However, by choosing suitable parameters, in particular the constant Af, it can be achieved that this viscosity measurement value r|i determined after the first iteration already agrees sufficiently well with the actual viscosity measurement value, so that further iterations can be dispensed with.

[0065] Fig. 3b, however, shows a general case of estimating the viscosity measurement value r|m as a function of the viscosity-relevant damping contribution D T1 in N iterations according to: hj+i = E(f, p) ■ D n 2 , where r|j+i is the viscosity value after the j+1-th iteration, where D ,j is the viscosity-relevant damping contribution after the j-th iteration, which can be determined as D m - D o - D f j. Df is the flow-dependent damping contribution after the j-th iteration, which can be determined as D f = A f ■ Rej ERe ■ Stj Est, where the Reynolds number Rei and the Stokes number Sti are determined according to the j

[0066] J -th iteration are determinable as Re 1, = f |m| and St, H di n R r|j 1

[0067] , where j is the index of the iteration steps with 0 < j < N. For example, rjo = 0 can be chosen as the starting value for the viscosity. With the same coefficients as in the example in Fig. 3a, this would lead to the starting value after the first iteration. With other coefficients, which actually better reflect the actual relationship between damping and viscosity, it takes several iterations until the viscosity measurement converges, as shown in Fig. 3b.

[0068] The accuracy of determining viscosity measurements depends on the flow rate and viscosity. Accordingly, the tolerance intervals outlined in Figs. 2a to d vary depending on the flow rate range, with the dotted lines representing absolute tolerance values ​​and the solid lines representing relative ones.

[0069] As a result, various mathematical approaches are possible to implement the principle of the invention. The aforementioned functions can be calculated explicitly, approximated by other functions, such as polynomials, or their values ​​can be stored in lookup tables and thus read out. Interpolation can be performed between the support points stored in the tables.

[0070] The essential steps of the method 100 according to the invention are shown in Fig. 4a to Fig. 4c and are finally explained in summary.

[0071] The method 100 requires the excitation 110 of an oscillation of at least one measuring tube in a bending vibration mode. Based on the oscillation behavior of the measuring tube, input variables for determining the viscosity measurement value are determined 120. This is followed by the determination 130 of a viscosity measurement value r|m for the medium based on the input variables.

[0072] Determining 120 the input variables comprises determining 122 a natural frequency measurement value f of the useful bending vibration mode; determining 124 a damping measurement value Dm of the useful bending vibration mode, for example by determining the ratio between the excitation current amplitude for exciting the useful bending vibration mode and a vibration amplitude achieved thereby; determining 126 a flow rate measurement value m of the medium; and determining 128 a density measurement value p of the medium based on the natural frequency measurement value f.

[0073] Determining 130 a viscosity measurement value r|m for the medium based on the damping measurement value Dm, the density measurement value p, the natural frequency measurement value f and the flow rate measurement value m comprises determining 132 a flow-dependent damping contribution in the manner described above, determining 134 a viscosity-relevant damping contribution, in particular by subtracting a natural damping contribution and the flow-dependent damping contribution from the damping measurement value and determining 136 the viscosity measurement value based on the viscosity-relevant damping contribution.

Claims

Patent claims 1 . A method for determining a viscosity measurement value of a medium by means of a Coriolis mass flowmeter having an oscillator with at least one measuring tube in which the medium is guided, wherein the measuring tube is to be excited by an exciter to oscillate in a bending vibration useful mode, wherein the oscillation experiences a damping D which depends on the viscosity and the density of the medium, as well as the mass flow rate of the medium and the natural frequency of the bending vibration useful mode of the measuring tube, the method comprising the following steps: Excitation of an oscillation of the measuring tube in a bending vibration useful mode; Determining a natural frequency measurement value f of the useful bending vibration mode; Determining a damping measurement value Dm of the bending vibration useful mode; Determining a flow rate measurement value m of the medium; Determining a density measurement value p of the medium; and Determining a viscosity measurement value r|m for the medium based on the damping measurement value Dm, the density measurement value p, the natural frequency measurement value f and the flow rate measurement value m, whereby the viscosity measurement value r|m can be determined as a function of a viscosity-relevant damping contribution D, which results from the damping measurement value Dm by subtracting an inherent damping contribution Do of the measuring tube and a flow-dependent damping contribution Df, i.e.: DTJ — D m — D o — Df, where the flow-dependent damping contribution Df can be represented by a model function MD which is proportional to a product of a power of a preliminary Reynolds number and a power of a preliminary Stokes number D f = A f ■ Re ERe ■ St Est, where a provisional viscosity value is included in the provisional Reynolds number Re and the provisional Stokes number St, where the following applies to the exponent of the Reynolds number Re: ER6> 1.3; where the following applies to the exponent of the Reynolds number Re: ER6< 1.8; where the following applies to the exponents of the stroke number St: Est > 3.4; where the following applies to the exponents of the stroke number St: Est < 4.0; where the quotient of the exponents Est / ERe > 2.2; where the quotient of the exponents Est / ERe < 2.8; where the viscosity measurement value r|m is determined using an algorithm which calculates the viscosity-relevant damping contribution D v determined by determining the flow-dependent damping contribution Df to the damping measured value Dm using the model function MD or using a routine that approximates the behavior of the model function MD, and where the viscosity measured value is determined as a function of the viscosity-relevant damping contribution D vis determined, whereby in particular in the case of approximation the routine approximates the behavior of the model function MD in such a way that a viscosity-relevant damping contribution D determined using the approximated flow-dependent damping contribution Df v for actual viscosities of not less than 10 cP by not more than 3% of a flow-dependent damping contribution Df determined using the model function v deviates from the latter.

2. The method according to claim 1, wherein a preliminary viscosity value is included in the preliminary Reynolds number Re and the preliminary Stokes number St; wherein the preliminary Reynolds number is given as: , and where the preliminary Stokes number is given as: where di is the inner diameter of the measuring tube, where n ris the number of measuring tubes, where r|v is the preliminary viscosity value.

3. Method according to one of the preceding claims, wherein the viscosity measurement value is determined iteratively on the basis of a damping measurement value in N iteration steps, wherein to determine the viscosity measurement value in the j+1-th iteration, the viscosity measurement value after the j-th iteration is included as a preliminary viscosity value in the calculation of the value of the preliminary Reynolds number and the preliminary Stokes number, where 0 < j < N.

4. Method according to one of the preceding claims, wherein a first preliminary viscosity value is determined on the basis of the damping measurement value Dm under the assumption that the flow-dependent damping contribution Df is zero, or where a first preliminary viscosity value is a literature value.

5. The method according to claim 1, wherein the natural frequency measured value f is not more than 600 Hz, for example not more than 400 Hz and in particular not more than 200 Hz.

6. Method according to one of the preceding claims: wherein a first preliminary viscosity value rjo is to be estimated with a first quadratic estimate, according to which the first preliminary viscosity value rjo is given as proportional to the square of the difference between the damping measurement value Dm and the self-damping contribution Do according to: where E(f,p) is a device-specific prefactor dependent on the natural frequency measured value and the density, where G is a device-specific constant, where a viscosity measurement value r|m is to be estimated with a second quadratic estimate, according to which the viscosity measurement value r|m is given as proportional to the square of the viscosity-relevant damping contribution D T1 according to: where the algorithm performs or approximates the first quadratic estimate and the second quadratic estimate, in particular, in the case of the approximation of the quadratic estimates, the approximated viscosity measurement value does not deviate from the latter by more than the maximum of 1 cP and 5% of the value according to the quadratic estimates.

7. Method according to one of the preceding claims, wherein the algorithm approximates the behavior of the model function MD by estimating the flow-dependent damping contribution Df as where Af is a device-specific constant.

8. Method according to one of the preceding claims, wherein a viscosity output value TJA is determined which is the minimum of the first preliminary viscosity value rjo and the viscosity measurement value r|m:

9. Method according to one of claims 1 to 5, wherein the viscosity measurement value r|m is determined as a function of the viscosity-relevant damping contribution D T1 with an estimate iteratively in N iterations is to be estimated according to: hj+i = E(f, p) ■ D Vij where r|j+i is the viscosity value after the j+1-th iteration, where E(f, p) is a device-specific prefactor dependent on the natural frequency measured value and the density, where G is a device-specific constant, where Dnj is the viscosity-relevant damping contribution after the j-th iteration, which can be determined as where Df is the flow-dependent damping contribution after the j-th iteration, which can be determined as D fJ = A f - Rej ERe ■ Stj Est where the Reynolds number Rei and the Stokes number Sti can be determined after the j-th iteration as where j is the index of the iteration steps with 0 < j < N, where N > 1 , in particular N>2, in particular the iteration is stopped if the viscosity measurement value does not change by more than 2% from one iteration step to the next; wherein the algorithm performs or approximates the estimation, in particular, in the case of approximation, the approximated viscosity measurement value does not deviate from the value of the estimate by more than the maximum of 1 cP and 5% of the value of the estimate.

10. Method according to one of the preceding claims, wherein a Reynolds number and / or a Stokes number is determined on the basis of the determined viscosity measurement value.

11. Method according to one of the preceding claims, wherein the viscosity output value and / or the viscosity measurement value for actual viscosities between 1 cP and 300 cP and flow velocities between 0 m / s and 5 m / s in a feed line to the Coriolis mass flowmeter with its nominal diameter does not deviate from the actual viscosity by more than the maximum of 10% and 1 cP upwards or downwards, and / or wherein the viscosity output value and / or the viscosity measurement value for actual viscosities between 1 cP and 300 cP and flow velocities between 5 m / s and 10 m / s in a feed line to the Coriolis mass flowmeter with its nominal diameter does not deviate from the actual viscosity by more than the maximum of 25% and 2.5 cP upwards or downwards.

12. Method according to one of the preceding claims, wherein the viscosity output value and / or the viscosity measurement value for actual viscosities between 300 cP and 1000 cP and flow velocities between 0 m / s and 3 m / s in a feed line to the Coriolis mass flowmeter with its nominal diameter, does not deviate from the actual viscosity by more than the maximum of 10% of the actual viscosity and 1 cP upwards or downwards, and / or wherein the viscosity output value and / or the viscosity measurement value for actual viscosities between 300 cP and 1000 cP and flow velocities between 3 m / s and 6 m / s in a feed line to the Coriolis mass flowmeter with its nominal diameter does not deviate from the actual viscosity by more than the maximum of 25% and 2.5 cP above or below the actual viscosity.

13. Method according to one of the preceding claims, wherein the damping measured value is determined on the basis of a ratio between an oscillation amplitude of the oscillation of the measuring tube in the bending oscillation useful mode and the amplitude of the associated excitation current, or on the basis of a width of a resonance curve Af / f of the bending oscillation useful mode, or on the basis of a decay time of the oscillation in the bending oscillation useful mode with the excitation power switched off, or reduced or increased.

14. Method according to one of the preceding claims, wherein the following applies to the exponent of the Reynolds number Re: ER6> 1.4; wherein the following applies to the exponent of the Reynolds number Re: ER6< 1.7; wherein the following applies to the exponents of the stroke number St: Est > 3.5; wherein the following applies to the exponents of the stroke number St: Est < 3.9; wherein the quotient of the exponents Est / ERe > 2.3; and wherein the quotient of the exponents Est / ERe < 2.

65.

15. A Coriolis mass flowmeter comprising: a measuring tube for conducting a flowable medium; an exciter for exciting an oscillation of the measuring tube in a bending oscillation mode; two oscillation sensors for detecting the oscillation of the measuring tube in a bending oscillation mode; and a measuring and operating circuit configured to carry out the method according to any one of the preceding claims.

Images