Acta Polytechnica Vol. 43 No. 4/2003 The Tangential Force Affecting the Radial Baffles in a Stirred Vessel: Analysis of the Macro-instability Related Component P. Hasal, I. Foit, J. Kratdna ExperirnenlaL data oblained fu measuring the langential component of the force affecting radial baffles in afl.at-bottotned cylindrical mixing aessel stirred uith a Rushton turbine irnpeller is analysed. Spectral analsis of the expetirnental data detnonstrated the presence of its macro-instability (MI) related low-frequenc) cornponent ernbedded in the total force. Two distinct dimetuionless frequencies (both directly 'proportional to the impel,ler speed of rotation N) of the occurrence of the MI component uere detected: a louer frequency of approximate$ 0.025N and a higher frequency of about 0.085N. The relatiae rnagnitude Qy1 of the Ml-related component of the total tangential force uas eaalualed by a cornbination of proper orthogonal de comlosition (POD) and spectral analysis. The aalues of magnitude QTal aaried in the interual frorn approximate$ 0.05 to 0.30. The rnagnitude Q7a1 tahes rnaxhnunt ualues at Low Rqnolds number aalues (in Laminar and transilianal regions). In the turbulent region (Reu > 20000) the Qp aalue is low and praclical$ constant. The dependence of the Q1a1 ualues on aertical position in the uessel is onLy marginaL. The results suggest tfutt the magnitude of the MI component of the force is signiJicantly inJluenced by the liquid uiscosity and densitl. KEuords: stirred uessel, bafJles, tangential force, macro-instabilit2, spectral ana$sis, proper orthogonal d.ecomposition. I Notation a(t) POD eigenmode [-] b ba{Ile width [m] B moving target width [m] C impeller off-bottom clearance [m] D impeller diameter [m] fur frequency of macro-instability [s-1] F tangential force affecting baflles [N] F* dimensionless tangential force affecting baffles [-] h impeller blade width [m] hT moving target height [m] H liquid filling height [m] Hr moving target axial position [m] N frequency of impeller revolution [s-'] NB number of impeller blades [-] Ns number of samples [-] Qvr relative magnitude of Ml-related fraction of tangential force [-] R autocovariance matrix [-] Retr impelier Reynolds number, Rer= (pND')h [-] ,io, variance of the MI component of the dimensionless force [-] rr1, variance of the total dimensionless force [-] sl-nn, variance of the non-Ml components of the dimens- ionless force [-] I time [s] T mixing vessel diameter [m] fs sampling period [s] V POD eigenvectors [-] W trajectory matrix [N] o POD eigenvalue [-] -p liquid density [kg m-'] 4 dynamic liquid viscosity [Pa s] Abbreviations: CFD computational fluid dynamics GC cold glycerol solution GH hot glycerol solution MI macro-instability PBI pitched blade impeller POD properorthogonal decomposition WA water 2 Introduction The liquid flow in mechanically stirred vessels has been studied intensively in recent decades. Numerous theoretical and experimental studies have been performed concerning various aspects of liquid flow in stirred vessels, e.g., the mean flow velocity, the intensity of turbulence, the energy dissipa- tion rate, spatial and temporal scales of the turbulent velocity field, etc. The liquid flow in a stirred vessel operated under steady operational conditions may be considered as a pseudo- -stationary high-dimensional dynamical system constituted by hierarchically ordered unsteady flows (vortices and ed- dies). The temporal and spatial scales of these flows span several decimal orders of magnitude. Recently, a pseudo-peri- odic large-scale flow has been identified in stirred vessels 55 Acta Polvtechnica Vol. 43 No. 412003 occurring with very low frequency and manifesting itself on spatial Scales cornpaiable to the size of the mixing vessel [-8]. This flow was named the marm-irutnbihry MI) of tht llow fuftern, and its existence has been confirmed by.a special.me- chanical measuring device [6], by laser-Doppler velocimetry (LDV) t?l, and by im4ge analysis methods [8]. The presence of the MI in the flow pattern is typically displayed by a distinct peak in the low-frequency part ofthe power spectrum ofthe local liquid velocity or other liquid,flow-related experimental data. Comprehensive reviews of a broad spectrum of the experimentally observed phenomena related to the macro- -instability were given in our previous PaPers [6, 7, 9-13]. 'l he macro-instability of the flow pattern yields strong impacts on mixing processes closely linked to fluid motion, e.g., on lo- cal mass- and heat-transport rates, local gas hold-up, the homogenisation rate, etc. [3, 5]. Bittorf and Kresta [14] iden- tified the macro-instability of the flow pattern as a mechanism responsible for liquid mixing outside the active volume of the primary liquid circulation loop. Although the existence of the macro-instability of the flow pattern in stirred vessels has been proved many times in experiments and CFD simula- tions its origin and development is not completely clear yet. Three possible reasons of the existence of the MI in stirred vessels with a Rushton turbine impeller have been hypothesised in the literature [15-17]: the first kind of insta- bility is associated with variations in the impeller ofl-bottom clearance. Montante et al [6] reported a variation in the flow at different clearances from the vessel bottom' The radial flow pattern was observed at high off-bottom clear- ances (C > 0.2 7) whilst the axial pattern was observed at C < 0.15 T Intermittent switching between axial and radial flow was detected at 0.15 T < C < 0.2 T. The changes in the flow pattern in the unstable cases occurred over time scales of the order of l0 min. The second kind of flow pattern instabil- iry is related to variations in the impeller rotational speed, when the flow pattern switches from mainly radial to mainly axial at a certain (critical) impeller speed (or Re") value. The third kind of macro-instabiliry arises from a "whirlpool" type of vortex (or coupled vortices) moving around the mixing vessel axis with the frequency of around 2 7o of the impeller fiequency N (see Nikiforaki et al [17])' The origin and the development of the MI of the flow pat- tern in the mixing vessels with the PBI was described by Brriha et al [6] and by Montes et al [7], who visually detected an axial flow originating in a pseudo-periodic manner at the vessel bottom and ascending along the baffle towards the liquid surface. The fiequency of the occurrence of the macro-instabiliry in all experiments reported in the literature [e.g., 6, 7, 9-l l, l7l was directly proportional to the fiequency of the impeller rotation. Howeve6 considerable scatter of the Proportionality constant values (the dimensionless MI fiequency) is apparent from the published data. Nikiforaki et al [17] have explained a part of this scatter by inadequacies in the experimental data treatment and interpretation methods adopted by different authors. Howeveq it is evident from the published data that at least two different macro-instability frequencies, correspond- ing to different mechanisms of MI origin, exist in the mixing vessels stinedwith a Rushton turbine' Nikiforaki et al [17] ob- served that higher MI frequency (about l0 7o of ltl) occurs at low impeller Reynolds number values, and that the lower 50 frequency. (about 2 Vo of N) prevails at'high Rer. Both fre- quencies were simultaneously present in the power spectra of the data measured at medium Reu values. The lower frbquency was assigned to intermittent switching between the single- and the double-loop flow pattern, while the higher frequency was assigned to precessional vortex motions. The macro-instability of the flow-pattern also exerts strong forces affecting the solid surfaces immersed in the stirred liquid, for example on the baflles, the draft tubes, the cooling and heating coils, etc. [2]. These forces may signifi- cantly affect the performance of the mixing vessel and, in certain cases, can even cause serious failures of the mixing equipment. An axially located impeller in a standard cylindrical mix- ing vessel equipped with radial baffles exhibits, in general, three force effects: axial, radial and tangential. The axial and the radial forces only slightly affect the radial baffles. Con- versely, it is the tangential force that exhibits the most of the dynamic pressune affecting them. The vertical distribution (along the ba{fle) of the tangential force in a standard mixing vessel was measured by Kratdna et al [12, 13] over a wide interval of the impeller Reynolds number value and in de- pendence on the impeller off-bottom clearance. Tivo kinds of the impeller were used: the pitched blade impeller with four or six blades, and a standard Rushton turbine' The pitched blade impellers elicited maximum force at the vessel bottom, and the magnitude of the force gradually decayed in the direction towards the liquid level. The Rushton turbine impellers yielded the maximum force at the impeller level, and the force rapidly vanished in both the below-impeller and the above-impeller region. Few attempts have been reported in the literature [18, l9] to separate the deterministic (low frequency and Ml-related) and the stochastic (turbulent) components of mixing phe- nomena in stirred tanks. Letellier et al [18] adopted a Hilbert transform-based procedure for separating the low-dimen- sional deterministic part of the experimental time series. Kovacs et al [19] used the Fourier transform for the same Pur- pose. Recentlywe have established a new technique (based on a combination of spectral analysis and proper orthogonal decomposition) for detecting the macro'instability of the flow pattern from the local velocity data, evaluating its relative magnitude and reconstructing its temporal evolution [9-l 1]. In this paper we apply this technique to the experimental data obtained by Kratdna et al [12, 13] by measuring the tangential comPonent of the force exerted on the radial baffles in the mixing vessel with a standard Rushton turbine impeller in order to identify and quantify its macro-instability related component. We also analyse the vertical distribution of the MI component in the vessel and the effects of the frequency of the impeller revolution and of the liquid viscosity on its magnitude. The frequency of occurrence of the macro- -instability related force component is also analysed. 3 Experimental A cylindrical flat-bottomed vessel with four radial baffles was used in all experiments rePorted in this paper. The di- mensions of the vessel and of the Rushton turbine impeller are marked in Fig. I and their values are listed in Thble l. Water (WA) and hot (GH)and cold (GC) aqueous glycerine Acta Polytechnica Vol. 43 No. 412003 Dimension Value Vessel diameter 7=0.3 m Liquid aspect ratio H/T= | Impeller to tank diameter ratio DIT=0.333 Relative impeller off-bottom clearance CIT=0.35 and 0.50 Measuring target relative axial position Relative baflle width blT = 0.1 Relative impeller blade height wlD=0.2 Relative impeller blade width llD =0.25 Number of impeller blades Ns= 6 Table l: Dimensions of the mixing vessel and the impeller Fig. l: The mixing vessel and the impeller geometry. The label TGT indicates the moving target positioned in the ba{Tle slit. solutions (viscosity 3 and 6 mPas, respectively) were used as working liquids. The fiequency of impeller rotation N was varied from 3.17 s-t to 6.0 s-r. The corresponding interval of Re" values (approximately from 6000 to 60000) covers the laminaq transition and turbulent regions. The tangential force affecting the baffles was measured by means of a trailing target (target height ht= l0 mm, target width B = 28 mm) located in a slit made in the baffle and enabled to rotate around an axis parallel to the vessel axis. The target was balanced by a couple of (calibrated) springs, and its angular displacement was scanned via a photo-elec- tronic sensor (see Kratdna et al [2, 13] for details of the measuring equipment). The sampling period was Ts= 20 ms. The signal from the sensor was recorded on a PC (after Ay'D conversion) and subsequently used for evaluation ofthe tan- gential force affecting the target. The duration of a single experiment was 20 minutes and N5 - 60000 samPles were typically stored. The local force measured at the target posi- tion was converted to a dimensionless force according to the relation F+F (l) ON'D' The time series of the dimensionless force at distinct locations HIH of the target along the baflle (see Thble l) were used for detection and analysis of the macro-instability related component of the total force. At each target position and relative impeller off-bottom clearance ClH, the value of Ren, was varied over the entire attainable inten'al. The total number of processed data sets was 3 10. 4 Nurnerical analysis of experimental data The numerical procedures used for experirrental dara analysis are here described oniy in basic features, as details can be found either in our previous reports [9-l l] or in the original papers [20, 2l]. The power spectra.of the time records of the measured tangential force were used to detect the presence of its low- -frequency component(s) generated by the macro-instabiiity of the flow pattern. The power spectral densities were evalu- ated by means of an algorithm based on the fast Fourier transform. A procedure based on an application of the Proper or- thogonal decomposition (POD) technique was used for extracting the Ml-related low-frequency component fiom the experimental data. Iirst, the raw data is centred by sub- tracting the mean force value. Then, by application of the so-called N-window [21] to the data and after some algebraic manipulations the so-called trajectory matrix W is obtained. Then, the auto-covariance matrix R of the trajectory matrix W is then evaluated by matrix multiplication R:WWT (2) and its (non-negative) eigenvalues clp and eigenvectors Vp (k=1,..., Af are determined. The eigenvalues, sorted in decreasing order, form the spectrum ofPOD eigenvalues and each ak expresses the magnitude of the contribution of the A-th eigenmode to the total force. 'fhe eigenvectors V; are used for evaluating the eigenmodes a1(l) of the measured data [0]. Then, the power spectra of the eigenmodes a1(l) are evaluated and the only eigenmodes with significant peaks in their power spectra, located exactly at the macro-instability frequency (or frequencies), are used for evaluating the rela- tive magnitude Qy1 of the macro'instability related compo- I --F4.| ---r--1Lit-1.- | I ^-,( ,/ ,a'l-Tr v -7-t tAltitp,l 5t Acta PolytechnicaVol. 43 No. 4l2OO3 75 H L.i a) /\ % b) I o.75 r)q o.25 o f/Iv nent of the total force affecting the baffles according to the relation tLro j ^ i eKur (3)e.ur =--6-,tzro j j=t where Ky1 is the set of indices of the eigenmodes contributing to the macro-instability. When two MI fiequencies (cf. Fig. 2) were present in the analysed data all POD eigenmodes con- tributing to both of them were included in the Kp11 set in Eq. (3), i.e., the total macro-instability related component of the force a{fecting the baffle was evaluated. The POD eigenmodes very fiequently captured simultaneously both MI frequencies, and therefore it was impossible to discriminate between their contributions to either of the individual MI fiequencies. Only a rough estimate of the relative ProPortion of the lower fiequency and the higher frequency sub-compo- nents of the total MI force component was possible in a few cases (see section Results and Discussion). o.2 0.4 f/N The value of the relative magnitude Q"' defined by Eq. (3) is expressed by means of the eigenvalues of the POD eigen- modes. Q", is therefore equal to the value of the ratio of the variance of the Ml-related force component to the total force variance, as shown by Lumley [20] or Aubry et al [21] o .25 o Fig. 2: The normalised power spectral densities of the measured tangential force affecting the baffle. a) CIH = 0.35, HrlH = 0.50' N = 3.333 r'', R"" = go05; b)clH = 0.35,HrlH = 0.50,N = 6.0 s-', Re" = 15889' ,9 errlr=}4!= o "MI, riot rial +si-l,,tt 5 Results and drscusslon There are nvo primary factors characterising the macro- -instability related component of the tangential force aflecting the ba{Iles: the frequency (or frequencies) of its occurrence/", and its relative magnitude Q"' (see Eq. 3). An analysis of the experimental data obtained in the mixing vessel with the pitched blade impeller has shown [l l] that the frequency of MI occurrence /", is directly proPor- tional to the impeller frequency N, i.e., the dimensionless fiequency/"/N attains a constant value. The ProPortionality constant (or dimensionless frequency) value was 0.074. The dimensionless frequency in the case of the PBI was independ- b) cM=0.50CIH=0.50 H .t+ffi &.ooo A o a I ftlqW8ary* e ' Fig. 3: The dimensionless frequency of the macro-instability component of the tangential force affecting the baffles as a function of the impeller Reynolds number.'al CIH =0.35,b) C/H= 0.5. Working liquids: tr,E... cold glycerol solution; O, O.'. hot glycerol solu- tion; O... water. 58 (4) a) 0.1 0.075 = . = o.os 0.075 . \: o'05 20000 40000 60000 R"u 20000 40000 60000 R"u ClH=0.35 bF odb-c oo'ic Acta Polytechnica Vol. 43 No. 4/2003 CIH=0.35 rl r o o q 7 I I E J e I ! U ( : -: 0 i- c - (] a d b) CIH=0.50 E t- 0 0.25 0.5 0.75 1 Hr/ H 0 0.25 0.5 0.75 1 Hr/ H Fig. 4: The dimensionless frequency of the macro-instability component o[ the tangential force aft'ecting the baffles as a^function of the measuring ta.get uerticai position. a) CIH =0.3'c,bl CIH =0.5. Working liquids: tr, 1... cold glycerol solution; O, O ... hot glycerol solution; O, O... water. a) 0.075 e . I o.os a -S oos ent of the impeller off-bottom clearance CIT of the impeller Reynolds number value, and of the vertical position in the vessel. No effects of liquid viscosity and density on the /^,' value were obsen'ed. Qualitatively distinct behaviour of the macro-instability was observed in the mixing vessel stirred with a Rushton turbine impeller. 'Iivo frequencies of a MI occurrence were typically detected in the laminar and transitional regions of liquid flow (see Frg. 2a) and a single MI frequency was typi- cally detected in the turbulent flow region (Fig. 2b). This observation agrees well with the experimental results of Niki- foraki et al [l7], and we can conclude that two mechanisms of MI formation were present in the mixing vessel studied in this paper. The lower dimensionless frequency of MI occurrence /r, varies about 0.02, and the upper dimensionless frequency value is close to 0.085. Both frequencies, however; are inde- pendent of the Re" value and of the position of the measuring targetHrlH, as documented in Figs.3 and 4. It can be de- duced from ligs. 2-4 that the MI occurrence frequency is slightly lower at higher impeller off-bottom clearance, i.e., at CIH =0.5. The mean values of the dimensionless fiequencies at different vessel configurations are summarised in Thble 2' In general, there is considerably higher scatter of the MI frequency data at the lower impeller position (ClH=0.35). The reason for this obseryation consists probably in the in- stability of the flow pattern arising due to the asymmetry of Tab. 2: Dimensionless frequencies of the macro-instability generated by a Rushton turbine impeller the below-impeller and the above-impeller liquid circulation loops and due to their irregular (erratic) interactions. The vertical distributions of the relative magnitude of the Ml-related component of the force affecting the baffles Qnn, (evaluated according Eq. 3) are shown in Fig. 5. Generally, higher values of Qvr are observed at low ReM values and a quite steep decrease is observed at ReM values just be- low approx. 20000. The Qn,,, values are almost constant in the highly turbulent region of liquid flou', despite the impeller off-bottom clearance and vertical position in the vessel. By comparing Figs. 5a,c and 5b,d it is easily recognisable that the dependence of the Q", value on Re' qualitatively diflers in the below-impeller and the above-impeller regions: the decrease of the Qnn, value with increasing Renn is more pro- nounced in the above-impeller region. Generally, the magni- tude of the Ml-related component of the tangential force gains slightly higher values in the above-impeller region at both impeller off-bottom clearances. The range of observed Qn,, values is quite broad, approximately from 0.05 to 0.30. The low-frequency component of the force affecting the ra- dial baflles thus represents a significant part ofthe total force, especially at low Re\r values (Renn<20000). Vertical distribu- tions of the Q", value in the vessel are depicted in Ftg. 6, where averaged Qn', values are plotted against HrlH (fot a fixed vertical position and a given working liquid all available Qn,, values - differing in corresponding Re*, -were averaged). Impeller off-bottom clearance Liquid Dimensionless MI frequency /y/,4 lower uPper CIH =0.35 GC 0.0258 (+0.0066) 0.0906 (r0.0069) GH 0.0243 (r0,0074) 0.0853 (r0.0070) WA 0.0275 (r0.0089) 0.0853 (r0.0040) CIH =0.50 GC 0.0165 (10.0026) 0.082 r (10.0042) GH 0.0150 (10.0024) 0.0815 (r0.0059) WA 0.0184 (10.0020) 59 Acta Polytechnica Vol. 43 No. 412003 : ol 0.3 0.15 0 0.3 0.1-5 0 0.3 0.15 0 0.3 0. l-5 0 0.3 0.15 0 b) 0.3 0.15 0 0.3 0.15 0 0.3 d 0.15 0 0.3 0.15 0 U.J 0.15 0 c)a) : ol d) R"vR"u \*. ,.-. II'/H = :'80 :a'^4. .-,+- +--'- < H,1/lI = 0.65 zt-- \"-/ __. _----__*.--_.-. .t- A- fl'/H = o'57 t":- _-a H1/H = 0.50 \ *\- l-=.*-{ \r.r,--rr-. .--:::. 0.3 0.15 0 0.3 0.15 0 0.3 0.15 0.3 0.15 0 0.3 U.I] 0 0.3 0.15 0 0.3 0.15 0 0.3 0.15 0 0.3 0.15 0 0.3 u. al 0 *\ Ht/ll = 0'55 t-l .!t" ti a *_- -- --- *< _ _ 0 _f\ H1/H = 0.57 I t-.-.\ i'---.--.-.-.y','^--{ lrira\ -A t---J I II*/II = 0.80 F.'--i€ s000 20000 ReM 40ooo IIr/H = 0.50 't\ \r\ ;. -....--,- . -. a---..-. 5000 40000 50000 fu/E = 0.35 F< . H^lfl = O.A2.\ j:..---:'--a H*-_*= --_+-"{ HrlH = 0.35 -&-- r..*-\ .t\ rrr/E = 0.27 Ha/H = 0.20 rriaf,-.-- ''\-H, IlnlH = 0.10 r--{----- ..<\ Ha/Il = 0 _-'-l :---*-----" R"u HrlH = 0.20 -*\-:.. F*--=.- Hrlll = 0.10 u. r--'-+-=-==-*--. l.*q_rr'-- Hrlll = 0 F--a---{r'.---=l 40000 50000 5000 50000 5000 Fig. b: Relative magnitude of the Ml-related component of the tangential force affecting the baffles. a) CIH =OS:b- rhe above impeller region;-b; CIH =0.35 - rhi below impeller region; c)-ClH =o.50 - the above impeller region; d) C/H= 0.50 - the below impeller region. Working liquids:l ... cold glycerol solution; O . '. hot glycerol solution; O ... water. baflle are averaged) are plotted in trig. 7 against the Reynolds number values. The Ml-related component of the force fol- lows the general trend observed in lig.5: the maximum (Qv) value is observed in the laminar region' then (with the in- creasing Re") the (Qr,') value monotonously falls until it reaches a practically constant value in the region of fully turbulent flow This trend repeats in all plots in lig. 7, except for trig. 7b, where (due to the enhanced data scatter an almost linear spline approximation is used). No effects of liquid vis- cosity or density on the (Q",) values are apparent from Flg' 7' As we have noted in the section on Numerical analysis of experimental data it is, in general, impossible to distin- guishbetween the contributions of the lower-frequency and The vertical (axial) profiles of the averaged Q", values are quite flat and the (Q"r) values are not apparently influenced by the impeller off-bottom clearance. (Qur) takes its highest values when cold glycerol is used as the working liquid and .decreases with decreasing liquid viscosity. Howeve4 this trend reflects the fact that the impeller Reynolds number value in- creases with decreasing viscosity, so the differences in the (Q",) values observed in Fig. 6 are primarily due to the in- creasing Re" value (cf. Iig. 5)' To elucidate the influence of the impeller off-bottom cleamnce on the (Q",) values and to emphasise the difference between the above- and below-impeller regions, sPatially av- eraged (Q",) values (at a fixed Re" value all Q",'s along the 60 Acta Polytechnica Vol. 43 No. 412403 a) tI i o.2 o. ) o o,3 o.2 o.l o o.I o o.5 Hrltt a) v.z 0.1 20000 40000 50000 Rerra b) n? 0.1 0 20000 40000 60000 Ret the higher-frequency sub-components to the total Ml-related force. Such a dissection was, howeve4 possible for several data sets obtained at the lowest Re" values, i.e., with the cold glycerol solution as a working liquid. The ratio of the lower i d)' + I o.3 0.3 o e)' 'i oroi I Fig. 6: The vertical distribution of the averaged relative magnitude of the Ml-related component (Qy1) of the tangential force in the vessel. a) CIH =0.35, cold glycerol solution (GC, f); b) CIH =0.35, hot glycerol solution (GH, O); c) ClH =0.35, water (Wd O); d) CIH=0.50, cold glycerol solution (GC, f); el CIH--0.50, hot glycerol solution (GH, O); f) CIH=0.50, water (WA, <)1. ftre errorbars display the interval of +1 S.D. from the mean values. c) : c) o.3 o- I o c) 0.3 0.1 o.5 Hr/u 20000 40000 50000 Rev d) ^ 0'2 : c) =o) Fig. 7: The averaged magnitude of the Ml-related component of the tangential force (Qur) as a function of the impeller Reynolds num- ber. a) CIH = 0.35, the above impeller region; b) CIH = 0.35, the below impeller region; c) CIH =0.50, the above impeller region; ' d) CIH=0.50, the below impeller region. Working liquids: 1... cold glycerol solution (GC), O... hot glycerol solution (GH), 1... warer (WA). U.I 0 20000 40000 60000 Ret frequency to the upper-frequency sub-component varied from approximately 0.725 to 1.5, i.e., both sub-components may (although not very significantly) dominate the macro- -instability related force acting on the baffles. The limited b) a a a c) f) a 6l Acta Polytechnica Vol. 43 No. 412003 amount of data available for the analysis, howeveq did not allow of any plausible conclusions to be deducted on the effects of the operational parameters of the vessel and the impeller on the ratio value. 5 Conclusions The analysis of the macro-instability related component ofthe tangential force affecting the radial baflles in a stirred vessel with a Rushton turbine impeller performed in this paper enables us to determine the frequency of its occurrence u.ra to evaluate its relative magnitude. The frequencies are linearly dependent on the ftequency of revolution of the im- peller.'Tivo MI fiequencies were observed at ReM< 20000' ihe Ml-related component of the force constitutes a highly significant contribution to the total force until the mixing is laini.ta. or weakly turbulent. The results of our study agree well with the findings of other authors obtained, howeve6 by analysing experimental data of a qualitatively distinct nature (e.g., Iiquid velocity components). Acknowledgment Frnancial suPport for this project by the Ministry of Edu- cation of the biech Republic is gratefully acknowledged (Projects: MSM223400007 and J04/98: 2 I 2200008)' References tll Bakker, A., Van den Akker, H. E' A: Singl'e-phase Flou in Stined Reartors. Trans. IChemE' Vol' 72A' 1994' p.583-593. tzl Chapple, D., Kresta, S.: The Effect of Geomatry on tho Stabil- ig of'Flou Patterns in Stined' Tanks. Chem' Eng' Sci'' Vol. 21, 1994, P. 3651-3660. t3l Haam, S., Brodkey, R. S., Fasano, J'B' I'ocal HeatTrans' fer in a Mixi'ng Vessel Using Heat Flux Sensors' Ind' Eng' Chem. Res., Vol. 31, 1992' p. 1384-1391' t4l Kresta, S. M', Wood, P.M.:Thc Mean Fl,ou Fizl'd Prod'uced by a 45" Pitchad Bln'd,e Turbine: Chnnges in the Circulntion i'attem Dn to Offbottom Clcarance. Can. J' Chem' Eng'' Vol.7l, 1993, P. 42-53. t5l Winardi, S., Nagase, Y.: Ilnstable Phenomcnan of Flou in a Mixing Vessel wi'th a Marine Propeller' J' Chem' Eng' Ja- pan, Vol. 24, 1991, P. 2$-2a9. t6l Brfrha, O', Foit, I., Smolka, P',Jahoda, M' Experimmktl Stuly of Twbulerrt Macro-instabil'iti'es t'n an Agt'tnted Slstem wxth'A;ial High+peed Impeller an'd uith Ral'in'l Bffies' Col- Iect. Czech.Chem. Commun.'Vol' 61, I996, p' 856-867' l7l Montes,J.-L., Boisson, H.-C., Foit, I',Jahoda, M':Veloc' it) FizW-Macro-i,nstabili'ties in an AxialLy Agttated Mixing iessel. Chem. Etg.J', Vol. 67' 1997, p' 139-145' t8l Guillard, F., Trdg&rdh, C., Fuchs, L' Natt Image An'alysis Methods for thz Sntfu of Mixing Pattems in Stined' Tanks' Can. J. Ch.*' Eng', Vol. 78, 2000, p' 273-285' tgl Hasal, P., Montes, J.-L., Boisson, !' C', Foit, I': Macro- -instabi.litics of Veloiitl Ficll' in Stined Vessel: Detection and Arnllsis. Chem. Eng. Sci., Vol. 55, 2000, p' 391-401' [10] Hasal, P., Foit, l.: Ma.cro-instabilitizs ofthe Flow Pattem in a Stined' Vessel: Detecti'on and' Chraatterintiun Using I'ocal ' Velocity Data. ActaPolytechnica, Vol. 40, 2000, p' 55-67 ' 62 [11] Hasal, P., Kratdna, J', Foit, L: Freqtawy and Magnitudt Analysis of tht Mono-instnbil?ry Related Cornpon'mt of tlu fangmtial Force Afecting Rnlin'l Bafus in a Stined Vessel' Acta Polytechnica, Vol' 42,2002' p. 59-69. [2] Kratdna, J., Foit, I., Rriiidka, M., Brfiha, O; Analysis of Dynmit itrus Afecting a Radin'l Boflb i" a Mnlwtiraily Agr,tated Systcm. Acta Polytech., Vol. 39, 1999' p' I l-38' [3] Kratdna, J., Foit, I., Br&ha, O.' Pavel, J.: Dishihution of Dynamit Presswe al'ong a Rad:inl Bafrlz in anAgttn'ted System with Stnnda'rd Rruhton Turbine Im,pellzr' "frans' IChemE, Vol. 79A, 2001, P. 819-823. I I 4] Bittorf, K. J., Kresta, S' M. : Aclite Volumz of M ean Circuln' tion for Stirred Tanhs Agitated uith Axinl In'pelbrs' Chem' Eng. Sci., Vol. 55, 2000, p. 1325-1335. [5] Yianneskis, M., Popiolek, 2., Whitelaw, I' H' An ExperimmtaL Stu\ of the Steo'd1 and UrctewJry Flow Chara'c- teristics of Stiffed Reactms. J. Fluid Mech', Vol' 175, 1987, p.537-555. [6] Montante, G., Lee, K. C., Brucato, A., Yianneskis, M': Double to Si:n'gle l-oop Ftou Pattem Transitions in Stinetl Vessek. Can. J. Ctrem. Eng., Vol. 77 ' 1999, p' 649-659' [7] Nikiforaki, L., Montance, G., Lee, K. C', Yianneskis, M': On thz Origtn, Frequmq and' Magnitudz of Macro-instabili- tics of the ilous in Stimed Vessek. Chem' Eng' Sci', Vol' 58' 2003, p. 2937-2949. [8] I-,etellier, C., Le Sceller, L., Gousbet, G', Lusseyran, F', Kemoun, A., lztar, B: Recoaering Determinist'i'c Behaoior from Experimmtal Thne Series in Mixing Reactor' AIChE J ' ' Vol. 43, 1997, P. 2194=2202. [19] Kovacs, T., Triiglrdh, C', Fuchs,'L': Fourier Spectrum to Recoaer Detertninistic and Stochastic Behaoior in Stirred Tanks. AJChEJ', Vol' 47,2001, p.2167-2176' [20] Lumley, J. L: Stochn^stb Tools in Turbulence' New York: Academic Press, 1970. [21] Aubry, N., Guyonnet, R., Lima, P.; SpatiotemporalAm'ly- sri oi Co*pt * Sigrnk: Theory and Applbations' J' Stat' Phys., Vol.64, 1991, P.683-739. Doc. Ing. Pavel Hasal, CSc. phone: +420224353 167 fax: *420 233 337 335 e-mail: Pavel.Hasal@vscht.cz Department of Chemical Engineering &. Center for Nonli.tea. Dynamics of Chemical and Biological Systems Institute of Chemical TechnologY Technickd 5 166 28 Prague 6, Czech RePublic Doc. Ing. Ivan Foit, DrSc. phone: +420224352713 iax: +420 224 310 292 e-mail: fort@fsid'cvut.cz Ing. Jiii Kratdna e-mail: Jiri.Kratena@fs.cvut.cz Czech Technical University in Prague Faculty of Mechanical Engineering Technick6 4 166 07 Prague 6,Czech RePublic Scan1 55 Scan1 56 Scan1 57 Scan1 58 Scan1 59 Scan1 60 Scan1 61 Scan1 62