Turbulence models for CFD in the 21st century

by

Brian Spalding, of CHAM Ltd

October, 2000

Invited lecture presented at ACFD 2000, Beijing


Example 5. Application of MFM to the prediction of the time-dependent reaction rate in a 3D paddle-stirred reactor.


[Extract from a lecture which can be found in full at www.cham.co.uk]

(a) The problem

Chemical industry makes much use of large reactor vessels, in which contact between the reactants is effected by mechanical stirring.

The productivity of the reactor depends greatly on the effectiveness of the micro-mixing; but until the arrival of MFMl, there was no practicable way of predicting this.

The prediction problem is complicated by the facts that the geometry is three-dimensional; and unsteady analysis is needed in order to represent properly the effects of the stirring paddle and of the rotation-impeding fixed baffles.

(b) The geometry

The reactor is one for which the US Dow Corporation has reported experimental data (only hydrodynamic data, alas!) which have been used in "benchmark studies".

The geometry and computational domain are shown below.

The impeller speed is 500 rpm, the dynamic laminar viscosity is 1.0cP and the water density is 1000 kg/m

The computational grid is divided into two parts, namely an inner part which rotates at the same speed as the impeller, and an outer part which is at rest.

the geometry

The total number of cells is 31365 (45 vertical, 41 radial and 17 circumferential).

the grid

(c) The starting condition



         __________|.|__________   

        |          |||          |  

        | upper    |.|          |  

        | liquid   |||          |  

        |          |.|  acid    |  

        |..........|||..........|     The paddle-stirred  

        | lower    |.| alkali   |         reactor         

        | liquid   |||          |  

        |       ---------       |  

        |paddle /////////       |  

        |           .           |  

        ------------|------------  

                    .< axis of     

                       rotation    

The sketch illustrates the apparatus and the initial state of the two liquids.

They are both at rest, and are separated by a horizontal interface

The paddle is supposed to be suddenly set in motion.

The computational task is to predict both the macro-mixing, represented by the subsequent distributions of time-average pressure, velocity and concentration, and the micro-mixing, ie the extent to which the two liquids are mixed together at any point.

(d) The turbulence model; combination of MFM and k-epsilon

There is no need, simply because MFM is available, to refrain entirely from using conventional models.

In the present study, therefore, the k-epsilon model was used as the source of the length-scale, effective-viscosity and micro-mixing information.

This allowed the power of the MFM to be concentrated on the process which it alone can simulate, namely micro-mixing and the subsequent chemical reaction.

An eleven-fluid model was employed, with mass-fraction of material from the top half of the reactor (ie the acid) as the population- distinguishing attribute.

This was probably not sufficient to permit achievement of grid-independent of solutions. However, one of the merits of MFM is that population-grid-refinement studies are easily performed.

(e) The fluid-population distributions

The problem can be regarded as a five-dimensional one; for it has three space dimensions, one time dimension, and one population dimension.

The simulation procedure therefore generates a very large volume of data, of which only a tiny fraction can be presented here.

What will first be shown is a series of PDFs. They all pertain to:

Specifically, PDFs will be shown for six radial locations, starting near the axis and moving outward to about one third of the radius.

Near the axis, most of the fluid is still alkaline

Farther from the axis, more-acidic fluids are found

This tendency increases with increase of radius

Also the shape of the PDF changes

The shape-change continues

This is the last PDF to be shown

(f) The salt distribution after 10 paddle rotations

The acid and the alkali are supposed to react chemically in accordance with the classic prescription:

acid + base -> salt + water

Of course, no chemical reaction can occur in the unmixed fluids, which (to use the parent-offspring analogy) are the Adam and Eve of the whole population. It is only their descendants, with both acid and alkali "in their blood", who can produce any salt.

This is why the prediction of the yield of salt necessitates knowledge of the concentrations of the offspring fluids.

The next picture shows the salt concentrations after 10 paddle rotations which the multi-fluid model has predicted.

These concentrations are the averages for all eleven fluids.

The salt concentrations predicted by the multi-fluid model

They are greatest in the region which has been most vigorously stirred by the paddle, which tends to thrust fluid downwards and outwards.

(g) The salt concentrations predicted by a single-fluid model

If no account is taken of the existence of the fluctuations, as is perforce the usual case, only the mean value of the acid/base ratio is available.

If the salt-concentration distribution is calculated from this, as the next picture will show, the values will be different from before.

Indeed it can be expected that they will be larger, because the implication of using a single-fluid model (which is what neglect of fluctuations amounts to) is that micro-mixing is perfect.

This is indeed what is revealed by the calculations.

The salt concentrations predicted by the single-fluid model

They are larger than the multi- fluid values, because micro- mixing is presumed (wrongly) to be perfect.

(h) Discussion of the results

The following conclusions appear to be justified by the studies conducted so far:-