Described below is a method for adapting the geometry of a dispersion nozzle in regard to a required size distribution of a phase dispersed in a dispersing phase by the dispersion nozzle.
The dispersion of substances that cannot be dissolved in each other, or are only partially dissolved in each other, such as gas in liquids or the preparation of oil-water emulsions, the gasification of bio- and chemical reactors and the like is a fundamental part of many industrial processes. In particular processes of this kind are required as key processes for producing multi-phase mixtures in the food industry, chemical industry, the pharmaceutical industry, petrochemistry and in mining (in the case of floatation processes). This requires the production of small and micro bubbles or droplets in sometimes very large volume and mass flows, for which significant amounts of energy are applied.
In particular production of the dispersed phases in a controlled size is required in all applications in order to obtain the desired properties of the dispersion/emulsion.
According to the related art, different types of dispersion nozzle are used for dispersion, in which intensive mixing of the phases to be dispersed takes place. Dispersion is achieved in these nozzles by way of a combination of regions with high shear rates alternating with regions of intensive turbulence for mixing of the phases.
The nozzles are designed according to empirical laws since no conclusive theory about the formation of bubbles or droplets in arrangements of this kind has existed until now.
Empirical and semi-empirical methods, such as the calculation of maximum stable bubble sizes above the critical Weber number of gas bubbles in liquids, can be used to only a very limited extent and in narrow parameter ranges.
Described below is a method of the type described in the introduction which allows particularly reliable adjustment of dispersion nozzles so a desired bubble or droplet size of the dispersion produced by the dispersion nozzle can be adjusted particularly reliably.
In a method of this kind for adjusting the geometry of a dispersion nozzle in regard to a required size distribution of a phase dispersed in a dispersing phase by the dispersion nozzle, a shear stress rate S and a relative velocity v0 between the phases is firstly calculated proceeding from a specified geometry of the dispersion nozzle. The shear stress rate is taken to mean the characteristic of the shear stress over a droplet or bubble of the dispersed phase. For a medium that flows with a linear velocity gradient, S corresponds to the quotient of the difference in velocity of the flowing medium over the extent of the droplet and the diameter of the droplet.
At least one local maximum radius for the dispersed phase is then determined on the basis of the flow conditions in the dispersion nozzle characterized in this way in accordance with the relation
R
b=(2σ/CSρLSv0)1/2
where σ indicates the surface tension of the dispersed phase, CS indicates the coefficient of friction of the dispersed phase in the dispersing phase and ρL indicates the density of the dispersing phase.
It has been found that in contrast to the Weber number, known from the related art, the maximum radius ascertained in this way allows a much improved estimate of the dispersion conditions to be made. The properties of the dispersion nozzle can be adjusted much more accurately in this way.
To analyze the dispersion properties over the entire dispersion nozzle a distribution of the local maximum stable radius is then determined over a cross-sectional area of the dispersion nozzle—this is also based on the relation given above and the flow conditions in the nozzle determined at the start.
If a specified maximum stable radius is exceeded in at least one region of the cross-sectional area then ultimately the geometry of the nozzle is changed such that a higher shear stress rate S and/or a higher relative velocity v0 of the phases is achieved at least in some regions.
A nozzle geometry is easily and accurately achieved hereby which during operation of a dispersion nozzle of this kind is capable of adjusting the desired dispersion properties.
It is expedient to perform the calculation of the flow conditions in (a) on the basis of a numerical flow model. Methods of this kind, known also as computational-fluid-dynamics models (CFD), allow a sufficiently detailed picture of the flow parameters in the dispersion nozzle to be obtained with reasonable computing effort. For particularly accurate determination of the flow properties it is expedient to also include the local degree of mixing of the phases in the calculation.
Accuracy can be improved further by weighting the distribution of the local maximum radius with the local portion of the dispersed phase.
The dispersion nozzle can be optimized particularly reliably if after changing the geometry of the dispersion nozzle in (d), then (a) to (d) are iteratively performed until the specified maximum stable radius is not exceeded in any region of the cross-sectional area. Iterative adjustment of this kind ensures that the dispersion nozzle has a geometry that satisfies the requirements made in the simplest way.
These and other aspects and advantages will become more apparent and more readily appreciated from the following description of the exemplary embodiments with reference to the drawings, in which:
Reference will now be made in detail to the preferred embodiments, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to like elements throughout.
A flow of liquid 12 is mixed with a flow of gas 14 in a dispersion nozzle 10 as is schematically shown in
As
It can be determined from these values whether a bubble 18 with a given radius is stable or is divided into smaller bubbles due to the shear forces.
For a bubble in a linear velocity gradient a shear stress rate S results as:
S=(vmax−v0)/Rb=(v0−vmin)/Rb=Δv/Rb (1)
In the virtually stationary state of equilibrium a pressure differential results over the bubble 18 of
Δp=σ(Rmin−1+Rmax−1) (2)
where Rmin and Rmax describe the short or long main axis of an ellipsoidal bubble 18. Assuming incompressibility of the bubble 18 a maximum effective pressure pmax of
p
max=⅛·ρlCS·(RbS+2v0)2 (3)
results on the side of the bubble 18 on which a flow with velocity Vmax acts, where Cs indicates the coefficient of friction of a sphere with a radius Rb. The minimum effective pressure can be determined analogously, so the pressure differential Δp over the bubble 18 of
Δp=CSρlSRbvo (4)
results, from which a resulting force
F
b
=C
SρlSRbvoA (5)
can in turn be obtained. Assuming an initially spherical bubble the area A is the effective cross-sectional area, so the force
F
b
=C
SρlSvoRb
results. Rb
In transition situations in which the bubble 18 is temporarily deformed by the pressure the bubble 18 initially assumes an oblate shape. The bubble 18 can become instable due to the excitation of shape oscillations and can break up into smaller bubbles if the contact area for the flow exceeds the critical area
A
crit
=R
2
b
critπ (7)
The critical radius can be estimated as
R
2
b
crit=1.44Rb
where Rb o indicates the initial bubble radius. A critical cross-sectional area therefore results of
A
crit=1.44Rb
With a maximum stable bubble 18 the following equilibrium of forces exists
F
b
=C
SρlSRbvoRb3v0π=Fst=2πσRb (10)
between the force Fb exerted by the flow and the surface force Fst. It therefore follows for the maximum stable radius Rb of a bubble 18
R
b=[(2σ)/(CSρlSvo)]1/2 (11)
Solutions to equation 11 are plotted for different local shear stress rates in
It is clear to see that the above-described non-empirical approach provides significantly different values for the maximum stable radius of a bubble 18. The semi-empirical approach assumes unrealistically small bubble radii in particular for high flow velocities, and it has not been possible to confirm these by way of experiments. Such velocities of several m/s to several tens of m/s are particularly important for industrial dispersion nozzles, however,
The typical operating parameters for a dispersion nozzle on a laboratory scale and a dispersion nozzle of an industrial floatation cell used in mining are superimposed in
The actually observed bubble radii in these dispersion nozzles lie at the operating points at 0.6-1 mm, and this coincides excellently with the calculated values according to
The respectively valid maximum stable radii for bubbles 18 can be calculated from the shear stress rates and velocities calculated in this way, as shown in
The local maxima of these curves are again well in line with the values of 0.6-1 mm determined by way of experiments.
It is possible to optimize the geometry of dispersion nozzles on the basis of the illustrated numerical simulation of the flow conditions in a dispersion nozzle and the calculation of the local maximum stable radii according to equation 11.
For this purpose the distribution of the local shear stress rates S, the relative velocities v0 of the phases and the local degree of mixing is firstly calculated for a specified geometry of the dispersion nozzle 10 and for the specified operating parameters, such as mass flows, volume flows or the like, in the manner described by fluid dynamic simulation. The distribution of the local maximum radii can be determined from equation 11. After weighting with the local dispersing agent content the distribution of the bubble or droplet sizes can then be determined over sectional planes of the dispersion nozzle 10 flowed through.
If this distribution differs from a required distribution of the bubble or droplet sizes, the geometric parameters of the dispersion nozzle are changed in such a way that in the case of calculated droplet or bubble radii that are too large, higher sheer stress rates and/or relative velocities are attained in fundamental parts of the dispersion nozzle 10.
This process can be iteratively repeated on the basis of the newly specified nozzle geometry until a nozzle geometry is obtained which produces the desired distribution of the droplet or bubble radii.
Dispersion nozzles 10 can be iteratively optimized quickly and reliably hereby.
A description has been provided with particular reference to preferred embodiments thereof and examples, but it will be understood that variations and modifications can be effected within the spirit and scope of the claims which may include the phrase “at least one of A, B and C” as an alternative expression that means one or more of A, B and C may be used, contrary to the holding in Superguide v. DIRECTV, 358 F3d 870, 69 USPQ2d 1865 (Fed. Cir. 2004).
Number | Date | Country | Kind |
---|---|---|---|
10 2012 209 342.7 | Jun 2012 | DE | national |
This application is the U.S. national stage of International Application No. PCT/EP2013/059504, filed May 7, 2013 and claims the benefit thereof. The International Application claims the benefits of German Application No. 102012209342.7 filed on Jun. 4, 2012, both applications are incorporated by reference herein in their entirety.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/EP2013/059504 | 5/7/2013 | WO | 00 |