The present application is related to and claims the priority benefit of German Patent Application No. 10 2019 115 215.1, filed on Jun. 5, 2019, and International Patent Application No. PCT/EP2020/062249, filed on May 4, 2020, the entire contents of which are incorporated herein by reference.
The present invention relates to a measuring device for determining the density, the mass flow, and/or the viscosity of a gas-charged liquid, a processing system having such a measuring device, and a method for monitoring a gas-charged liquid. It is a measuring device having at least one oscillator which has at least one vibratable measuring tube for conducting the medium, and which has at least two vibrational modes whose natural frequencies depend upon the average density of the gas-charged liquid and upon the gas volume fraction of the gas-charged liquid.
It is known that density and flow measurement, due to the resonator effect, have cross-sensitivities to the gas charge with micro-bubbles. The first published patent application DE 10 2015 122 661 A1 discloses a method for compensating for the influence of this cross-sensitivity, wherein the natural frequencies of two bending vibrational modes are evaluated for this purpose. The result of this method is a value for the average density of the gas-charged liquid and an associated mass flow measured value.
The first published patent application DE 10 2016 114 972 A1 discloses a method for determining the gas volume fraction of a gas-charged liquid which is linked to the aforementioned method, wherein, here, a pressure measured value for the mixture is also incorporated into the calculation of the gas volume fraction, in addition to the natural frequencies of the two vibrational modes. If the gas volume fraction and the average density of the gas-charged liquid are known, the density of the liquid phase can also be calculated.
The above approaches are valid to the extent that the suspended micro-bubbles are suspended essentially homogeneously in the liquid, i.e., the medium acts macroscopically as a homogeneous mixture. However, if larger free bubbles occur which macroscopically lead to spatially heterogeneous distributions, the methods, which are based upon the assumption of a homogeneous, compressible mixture, meet their limits. Although the properties of the macroscopically homogeneous mixture are still recorded reasonably accurately, the gas charge in the free bubbles is not truly identifiable by these methods.
Spatially heterogeneous distributions of gas inclusions in a flowing medium lead to temporal density fluctuations at the location of a flow-through measuring device, which can be detected as fluctuations in the natural frequency of bending vibrational modes. The as yet unpublished patent applications DE 10 2018 112 002.8 and DE 10 2019 003 075.3 describe measuring devices and methods that make it possible to draw conclusions, on the basis of the frequency fluctuations, about the presence of free bubbles.
The aim of the present invention is to provide a measuring device and a method which are suitable for characterizing a gas-charged liquid over a wide range of mixing ratios. The aim of the present invention is furthermore to provide a processing system which is suitable for processing a gas-charged liquid. The aim is achieved according to the invention by the measuring device, the processing system, and the method according to the present disclosure.
The measuring device according to the invention comprises:
The specification of two media state values gives a system operator a more comprehensive picture about the media state, even in cases in which the model underlying the individual media state values is no longer sufficient for the description of all aspects of the media state. Processing the medium can thus be controlled in a more differentiated manner.
In a development of the invention, the first media state value comprises a value dependent upon a fraction of free gas bubbles in the medium, under the assumption that free gas bubbles are causative for the fluctuations. In particular, the teachings of the as yet unpublished patent applications DE 10 2018 112 002.8 and DE 10 2019 003 075.3 can be used for this purpose, which make it possible to draw conclusions about the presence of free bubbles on the basis of the frequency fluctuations.
In a development of the invention, the second media state value depends upon a gas volume fraction in the form of micro-bubbles in the medium. Details in this regard are described in, for example, DE 10 2015 122 661 A1 and DE 10 2016 114 972 A1.
In a development of the invention, the second media state value indicates a gas volume fraction in the form of micro-bubbles in the medium, under the assumption that an observed ratio of the natural frequencies is due to a resonator effect in which the medium compressible due to the micro-bubbles vibrates against the measuring tube.
In a development of the invention, the operating and evaluating circuit is configured to further determine the second media state value as a function of a pressure measured value, which represents the media pressure in the measuring tube. For example, the teaching of DE 10 2016 114 972 A1 can be used for this purpose.
A measuring point that has the measuring device according to the invention accordingly comprises a pressure sensor, in order to be able to provide the measuring device with a pressure measured value.
In a development of the invention, the media state message contains, in addition to the first media state value, information about which flow regime of the first media state value was determined, wherein the information about the flow regime comprises a flow measured value or a value range of the flow measured value, wherein the flow measured value comprises in particular a flow velocity, a mass flow rate, a volume flow rate or a Reynolds number, or a range specification for the aforementioned variables.
In a development of the invention, the first media state value is normalized with a normalization function that is dependent upon a flow measured value.
In a development of the invention, the first media state value is a function of the fluctuations in the natural frequency, wherein the function furthermore exhibits a natural-frequency-dependent normalization.
In a development of the invention, the first media state value comprises a gas volume fraction of the medium in the form of free gas bubbles, or a range of values for the gas volume fraction.
The processing system according to the invention serves to process a medium which has a gas-charged liquid, wherein the processing system comprises:
In a further development, the actuator comprises a temperature control device for heating and/or cooling the medium, a pump, a valve, a ventilator, a stirring device, or a mixer.
The method according to the invention for monitoring a gas-charged liquid with a measuring device for measuring the density, the mass flow, and/or the viscosity of a medium, which has at least one oscillator with at least one measuring tube for conducting the gas-charged liquid, in particular with a measuring device according to the invention, comprises the following steps:
The invention is explained in the following in further detail on the basis of the exemplary embodiments shown in the figures. The following are shown:
The first exemplary embodiment of a measuring device 1 according to the invention shown in
By means of a density meter that has an oscillator with at least one vibratable measuring tube for conducting the medium, the density ρ of a medium can be determined on the basis of a mode-specific, density-dependent natural frequency fi of the oscillator according to:
The coefficients c0,i and c1,i are mode-specific coefficients, which are preferably to be determined for each measuring device type or for each measuring device. The coefficient c0,i is characterized by the mass of the measuring tube conducting the medium, while the coefficient c1,i depends upon a mode-specific stiffness of the measuring tube. The time derivative of the density
is thus given as:
The time derivative of the density
is a suitable measure for describing the density fluctuation. In order to determine this value, the observed frequency fluctuation
of the vibrating measuring tube or the vibrating measuring tubes also has to be multiplied by a normalization factor
In this way, the basis is created for an evaluation function which can describe the degree of inhomogeneity of the medium in the form of density fluctuations, independently of the particular type of density meter or its size. In one embodiment of the invention, the operating and evaluating circuit 30 of the above exemplary embodiment of a measuring device according to the invention is configured to carry out the density fluctuation on the basis of the frequency fluctuation by means of the above-described normalization with the reciprocal of the third power of the mode-specific natural frequency:
The effect of the normalization can be illustrated using data for two Coriolis mass flow meters from the applicant's company, viz., a Promass F50 and a Promass Q50, both of which have the function of a density meter. The observed natural frequency fluctuations
differ by a factor of about 6.6 in the case of an aqueous medium with a gas charge of 1% or 2%. After normalization with the normalization factor
this results in approximately the same value for the density fluctuation
An equivalent analysis of the density fluctuation
is implemented in a second development of the invention. In this case, the operating and evaluating circuit is designed to determine the density fluctuation according to:
In order to provide the magnitude of the relative density fluctuation
the operating and evaluating circuit is designed according to a third development of the invention to determine the magnitude on the basis of the relative frequency fluctuation
according to:
If the density of the medium at a measuring point varies by only a few percent around a known value, but is otherwise known from the value range, the relative density fluctuation can be estimated as a function of the relative frequency fluctuation using:
wherein a is a measuring-point-specific or medium-specific, and optionally mode-specific, constant, provided that several modes are used for density measurement. The described normalizations with the third power of the frequency for the density fluctuation or with the frequency itself for the relative density fluctuation are advantageous, but not absolutely necessary for the implementation of the present invention.
In measuring mode, a first estate value can then be assigned to a frequency fluctuation df/dt as a function of the flow velocity V that results from a momentary mass flow rate and the current mean density of the medium. Depending upon the flow velocity, two different media state values M11 and M12 result for the same observed frequency fluctuation fl5, as shown in
On the one hand, the first media estate value may correspond to the gas volume fraction a1 or be a linear function of the gas volume fraction. This can be useful, for example, if the gas volume fraction a1 is a process variable to be monitored or regulated.
On the other hand, the first media state value may also correspond to exceeding a critical monitoring variable; for example, all frequency fluctuations that correspond to a gas volume fraction below a3 would be mapped to a first media state value M12, while all frequency fluctuations that correspond to a gas volume fraction above a3 would be mapped to a first media state value M11, wherein M11 corresponds to a critical state.
The first media state value M1 is to be detected, communicated, or signaled in a suitable manner.
In a simpler embodiment of the invention, tuples of frequency fluctuations and associated velocities are output as first media state values.
In the context of the discussion of
Below, the steps for determining a second media state value are described with reference to
In a step 430, which is described further below with reference to
As shown in
Subsequently, in a step 432, a value of the sound velocity c is determined, which, at the measured natural frequencies f1 and f2 of the bending vibrational modes, leads to the observed ratio V of the preliminary density values in the following equation:
On the basis of the determined sound velocity value, in step 433 of the method in
A mixing density value ρmix-coriolis corrected for the resonator effect can finally be calculated as:
A corresponding mass flow correction term Km for calculating a mass flow measured value rhmix-coriolis corrected for the influence of the resonator effect can be calculated in step 434 as:
The mass flow measured value rhmix-coriolis corrected for the influence of the resonator effect then results as:
In this case, rh is a preliminary mass flow measured value which results from multiplying a calibration factor for the measurement sensor by the phase angle between the first bending vibrational mode and the Coriolis mode.
In the following, it will be explained how the gas volume fraction in the form of micro-bubbles and the density of the liquid can be determined.
The following known relationship between the sound velocity of a gas-charged liquid and further parameters exists, according to Sorokin:
In this case, α is the gas volume fraction (or the gas void fraction GVF), cg is the sound velocity of the pure gas, ci is the sound velocity of the pure liquid, γ is the adiabate coefficient for the gas, p is the current pressure of the gas-charged liquid, and ρl is the density of the gas-charged liquid.
In fact, the sum in parentheses on the right side of the equation over the pressure range, in which gas charges are relevant as disturbance variables, and at fluid densities on the order of 103 kg/m3, is essentially determined by the third sum. The above expression is thus reduced to:
The gas volume fraction α can therefore be estimated as a function of the other parameters in the equation. If the quotient from the already determined Coriolis mixture density ρmix-coriolis divided by the liquid volume fraction, i.e., (1-α), is set for the liquid density ρl, the following results:
For this purpose, as shown in
The adiabate coefficient γ=cp/cv=(f+2)/f, where f is the number of molecular degrees of freedom of the gas, is, at room temperature, 1.4 for nitrogen and dry air, and 1.3 for methane, for example. Depending upon which gas or gas mixture is to be expected, a corresponding adiabatic coefficient can be used. On the basis of these variables, the gas volume fraction α is then calculated in step 442.
As shown in
ρliquidGVF=ρmix-coriolis/(1−α)
This liquid density is of interest in that, for example, it enables an estimation of the fractions of the components for liquids that comprise a solution or homogeneous mixture of two components of different densities.
The gas volume fraction α, the mixture density ρmix-coriolis, or the liquid density ρliquidGVF of the gas-charged liquid can be output as the second media state value. Optionally, several of these parameters can be output as further media state values, wherein a user can select which parameter(s) is/are to output as the second or further media state value. The operating and evaluating circuit is preferably configured in such a way that a user can freely select the parameter(s) to be output with the media state message.
The example shown in
While the processing system is operating, various phenomena occur which can cause a gas charge in the medium 530 which impairs the product quality.
First, the agitator 524 can stir macroscopic, free bubbles into the liquid phase if the fill-level of the mixture 530 in the reactor vessel 510 is too low and if the residual gas charge is too large. Likewise, air bubbles can be sucked in by the pump 518 if the reactor vessel is emptied too far.
Secondly, during a reaction or upon mixing the components, micro-bubbles may arise, which are suspended in the liquid phase. In order to expel the micro-bubbles, the space 532 in the reactor 510 above the mixture 530 can be evacuated with the low vacuum pump 522. The viscosity of the mixture 530 can be reduced via the temperature control device 526. In addition, the mixture 530 can be circulated with the agitator in order to facilitate the expulsion of the micro-bubbles.
The Coriolis mass-flow and density-measuring device 540 according to the invention is configured to be able to detect and distinguish both types of gas charges. The first media state value signals the presence of free bubbles at the frequency fluctuation observable when adding the medium 530. By contrast, the first media state value is “blind” to homogeneously-distributed micro-bubbles, since these do not cause frequency fluctuations. Depending upon the embodiment, the second media state value can have a mixture density value, which can already be indicative of a gas charge with micro-bubbles, due to the deviation from a reference density for the liquid phase, or an explicit value for the gas volume fraction. However, due to the underlying model, the second media state value is not able to appropriately describe the fraction of free gas bubbles. In combination, however, the two state values provide a comprehensive description of the media state with regard to the gas charge, so that a system operator can take this into account in the process control by operating the actuators in such a way that the gas charge corresponds to a target value range.
In the exemplary embodiment of the method 600 according to the invention shown in
Subsequently, for at least one of the natural frequencies, frequency fluctuations are determined in a second step 620, which enable the determination of a first media state value in a third step 630, as described above in connection with
In a fourth step 640, a second media state value, which can in particular comprise a mixture density or a gas volume fraction, is determined on the basis of the two natural frequencies and, if applicable, a pressure measured value. Details in this regard are described above in connection with
In a fifth step 650, a media state message is output into which the two media state values are incorporated. Optionally, in a sixth step 660, a process can then be controlled as a function of the media state message in order to regulate a gas charge in a liquid to a target value.
Number | Date | Country | Kind |
---|---|---|---|
10 2019 115 215.1 | Jun 2019 | DE | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2020/062249 | 5/4/2020 | WO |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2020/244855 | 12/10/2020 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
20110023626 | Weinstein | Feb 2011 | A1 |
Number | Date | Country |
---|---|---|
102009000749 | Aug 2010 | DE |
102015112737 | Feb 2017 | DE |
102016007905 | Jan 2018 | DE |
102017115251 | Jan 2019 | DE |
102018101923 | May 2019 | DE |
2275787 | Jan 2011 | EP |
2286189 | Feb 2011 | EP |
2019219321 | Nov 2019 | WO |
Number | Date | Country | |
---|---|---|---|
20220334038 A1 | Oct 2022 | US |