Microsoft Word - 2012_Dr_Bodor_Endre_HJIC.doc HUNGARIAN JOURNAL OF INDUSTRIAL CHEMISTRY VESZPRÉM Vol. 39(3) pp. 335-339 (2011) APPLICATION OF MODELS WITH DIFFERENT COMPLEXITY FOR A STIRRED TANK REACTOR A. EGEDY , T. VARGA, T. CHOVÁN University of Pannonia, Department of Process Engineering, Egyetem St. 10. H-8200 Veszprém, HUNGARY E-mail: egedya@fmt.uni-pannon.hu Engineering problem solving such as process design, process optimization, safety analysis, etc.; relies widely on mathematical models of the process. To solve various engineering problems various models with different complexity are needed. A stirred tank reactor with a highly exothermic reaction is studied in this work, because in the modern chemical technologies mixing is one of the most important operations, and stirred reactors are widely used in industrial applications. The stirring system of a mixed tank is always an important aspect of the design, because the involved processes (such as reactions, heat and component transport) usually require proper contact and homogeneity of the existing phases. For the suitable homogeneity the design and the size of the moving parts are also important problems. In certain situations attachment of static parts to a stirred tank (such as baffles) may have an important effect too. The primary goal of this study is to create models with different level of complexity and determine which model is the best suited for solving different engineering tasks such as process design, scale-up, or optimisation. etc. To determine which model fits best for a problem, mathematical models were created and compared to find out, how the information can be extracted from these models and be applied to solve engineering problems. Three types of models have been developed: perfectly mixed reactor model, compartment model, and Computational Fluid Dynamics (CFD) models with different dimensions. The reaction of hydrogen peroxide with sodium thiosulphate in a continuously stirred tank reactor is analysed as a case study. The perfectly mixed vessel models and compartment models were solved in MATLAB/SIMULINK program package. The CFD models were implemented in COMSOL Multiphysics. Keywords: computational fluid dynamics, stirred tank, hydrogen peroxide, model complexity Introduction In modern technologies mixing is one of the most important operations. The stirring system of a mixed tank is always an important aspect of the design, because the involved processes (such as reactions, heat or component transport process) require proper contact and homogeneity of the existing phases [1]. Nowadays there are an increasing number of applications of model based methods in the industry too. With a correctly built model there is a possibility to examine the dynamic behaviour of the system as well as the mass, heat and momentum transport processes. In this study a kinetic model of the reactions of hydrogen-peroxide and sodium thiosulphate from an earlier case-study is analysed [2]. For developing suitable reactor models we followed a sort of hierarchical modelling concept [3]. Three types of models with different complexities have been developed: perfectly mixed reactor model, compartment model, and Computational Fluid Dynamics (CFD) models. Since the examined reaction is highly exothermic the thermal runaway problem must have been investigated thoroughly. For achieving different goals different models might be required. To solve the simplest problems such as preliminary design and control, engineers can use simple models to describe a dynamical system, such as perfectly mixed vessel model, or ideal plug flow reactor model. In the case of a perfectly mixed model the whole vessel is homogenous, and models with concentrated parameters can be used. On the other hand the plug flow reactor models can be applied to calculate the inhomogeneity in a tube reactor in axial direction; in this case the distributed model parameters need to be defined [4]. If a simple model is not detailed enough to solve a problem more complex model have to be created, such as operation fermentation reactor with multiple impellers. The second level of models which is discussed in this paper is the compartment models. These models generally contain two or more well defined compartments, built from perfectly mixed reactor or ideal plug flow reactor models. The most important thing at this level of modelling is to find the proper parameters describing the connections between these models [5]. For example to create a compartment model for a stirred vessel in the basic cases the vessel has to be separated at least two different compartments. The impeller and the region near the impeller can be described by using a perfectly mixed reactor model, and at least one other region must be defined, called circulation loops, or circulation region [6]. To model the connections between these models the 336 following parameters need to be determined: the exact volume of the compartments, and circulation rates, or other pre-defined parameters [7]. The more difficult the problem is, the more complex model must be created. In the last few years researchers have been trying to create even hybrid CFD-compartment models that are capable of handling more than 100 compartments [8]. The CFD models are the most complex models examined in this work. The CFD models can be applied to determine the entire hydrodynamics of the system, to define flow patterns, or concentration distributions, or thermal hotspots. With validated CFD models the system can be examined at a completely new level, and engineers can model anomalies such as thermal runaway, and prevent hazardous situations [9]. CFD models can also serve as excellent design tools [10]. In this paper we have created models with different complexities in order to determine how these models can be used to solve various engineering problems in the field of design and control of chemical processes. The perfectly mixed vessel models and compartment models were solved in MATLAB/SIMULINK program package. The CFD models were implemented in COMSOL Multiphysics. Modelling approach At every level of modelling the same exothermic reaction with a well-known kinetics is used [2]. This is the reaction of hydrogen-peroxide with sodium thiosulphate. OHSONaOSNaOHOSNa 24264222322 4242 ++→+ (1) The analysed highly exothermic reaction is the oxidation of sodium thiosulphate with hydrogen peroxide. The kinetic parameters of the analysed reaction were determined in an earlier case-study [2]. Table 1 contains the model parameters. Table 1: Model parameters Revolution speed 20 [1/min] Density 1000 [kg/m3] Volume 0.375 [m3] Feed 2e-4 [m3/s] Viscosity 0.001 [Pas] Component diffusion CA - sodium thiosulfate, CB - hydrogen-peroxide, CC - sodium - trithionate 10-5 [m2/s] The examined stirred vessel is a continuous tank reactor. All the implemented models have the same model parameters for describing reactor geometry (apart from the parameters of the impeller, since the impeller was implemented only in the CFD simulations). Because of the exothermic reaction the temperature dependence of the reaction is high, thus it is necessary to define heat balance equations too. In every model component and heat balances was implemented. Figure 1: Scheme of the reactor Component balances { }OHSONaOSNaOHOSNai rcc V B dt dc iibei i 24264222322 , ,,,, )( = ⋅+−⋅= υ (2) Heat balance )()( 1 TT V B Hr cdt dT IN p −⋅+Δ−⋅⋅= ρ (3) In CFD models a momentum balance was also added besides component and heat balances. Momentum balance FpIuu t u ++−∇=⋅∇⋅+ ∂ ∂ ][)( τρρ (4) Results and discussion Perfectly mixed reactor model In the perfectly mixed reactor model a mathematical formulation based on the modified Dahmköhler-equations is applied to describe material and heat balances. The model was implemented and solved in MATLAB. Fig. 2 shows the results of the dynamic simulation of the perfectly mixed reactor model with the previously mentioned reaction implemented. The perfectly mixed reactor model can be used to: ● calculate the maximum productivity and conversion; ● support optimization and design of the reactor; ● study control functions; ● examine runaway situations in the system level; ● etc. The figure shows, that the product concentration reaches a stationary value within a short time, but, the temperature rises 45 K, therefore it is important to control the reaction temperature to avoid hazardous situations. 337 Figure 2: Results of the perfectly mixed reactor model Compartment model At this stage of the work four types of compartment are applied to model macro mixing conditions. The following compartments are applied: ● perfectly mixed reactor; ● mixer; ● distributor; ● dead zone. The first cell is identical to the previously discussed perfectly stirred reactor model. The mixer and the distributor have zero volume, and these compartments are capable of modelling circulation loops. The distributor distribute the flow directly, and the mixer use mixing equation to calculate the outlet variables. The output variables were computed from the input variables using algebraic equations. The dead zone cell represents a cell with real volume; however there is no flow in it. The circulation numbers – describe the ratio of distribution – have constant values at this time and describe the rate of circulation between impeller zone and circulation zones. Fig. 3 represents the compartment approach and Fig. 4 shows the structure of the compartment model. Eight compartments have been used: one dead zone, three perfectly mixed reactor models for impeller top and bottom circulation zones, two mixers and two distributors. The volume of the dead zone is 15 % of the total volume, and each reactor cell has the same volume. Figure 3: Compartment model for a stirred tank reactor The compartment model can be used to: ● support advanced control and design; ● examine macro mixing and model complex reaction systems; ● examine heat transfer at a more complex level; ● etc. Dead zone Top circulation zone Impeller zone Bottom circulation zone Mixer Mixer Distributor Distributor Feed Outlet Figure 4: Compartment model structure Two- and three-dimension CFD models To simulate the behaviour of the system in a CFD environment COMSOL Multiphysics program package was applied. The investigated impeller is a standard six bladed Ruston turbine in the 3D simulations and a simple blade in the 2D analysis. The revolution of the impeller is 20 RPM in the first investigation. The revolution of the impeller is low, because a laminar model was used in this case, and laminar models can only work at low impeller revolution. In the future a turbulent model will be implemented to examine higher revolution speed. In case of the 2D analysis the same continuous system cannot be simulated, since the investigated geometry is a slice of the 3D geometry, i.e. the slice in the centre of impeller zone. In the 3D simulations the same continuous reactor is investigated as in the two dimensional simulations. The results of the CFD simulation are shown in Fig. 5 and Fig. 6. The velocity field are calculated in 2 and 3D. The CFD models can be used for: ● examining of runaway reactions in three dimension; ● studying different impeller geometries; ● testing of new vessel constructions; ● applying multi-level optimisation; ● etc. 338 Figure 5: Results of the 2D CFD model Figure 6: Results of the 3D CFD model Comparison of the models In the final step of our investigation the simulation results of the developed three models are compared. The results are shown in Fig. 7. The concentration refers to the concentration of the product. The results of the CFD model are much more different than the other two, because in this model the real geometry of the reactor has been implemented, and this makes the system more inhomogeneous. Because the lack of the validation the question which model is the best is hard to answer. To obtain validated models hydrodynamic measurements will be taken in a glass reactor at our department. Figure 7: Comparison of the models Conclusion The predefined goal – to create models with different complexity – has been achieved. The developed models can be applied in a wide range of industrial purposes, from the design of a reactor to control or solving optimisation problems. At every level of modelling some of the possible applications have been discussed. This approach is expected to be used to determine which level of hierarchical modelling have to be applied to solve a given engineering problem, and to implement the simplest model suitable to work at this level. It is important that all the simulation results need to be validated based on experimental data from real systems. The physical experiments will be the next step in our work. ACKNOWLEDGEMENT The financial support from the TAMOP-4.2.2- 08/1/2008-0018 (Liveable environment and healthier people – Bioinnovation and Green Technology research at the University of Pannonia) project is gratefully acknowledged. APPENDIX Appendix 1: The construction parameters of the reactor Impeller length 0.6 [m] Impeller width 0.1 [m] Impeller height 0.01 [m] Inner diameter 1 [m] Outer diameter 1.2 [m] Vessel height 1 [m] Feed diameter 0.02 [m] 339 NOTATION Name Description Unit cA hydrogen-peroxide concentration [kmol/m 3] cB sodium - thiosulfate concentration [kmol/m 3] cC sodium - trithionate concentration [kmol/m 3] cin inlet concentration [kmol/m3] V reactor volume m3 B Feed [m3/h] k Reaction rate 1/s ρ Density [kg/m3] cp Heat capacity [J/mol/K] TIN Feed temperature [K] T Reactor temperature [K] u velocity vector [m/s] t time [s] τ shear force tensor g grav. acceleration [m/s 2] differentialoperator F Force vector N ∇ REFERENCES 1. E. L. PAUL, V. A ATIEMO-OBENG, S. M. KRESTA: Handbook of industrial mixing science and practice, Wiley-Interscience, (2004) 2. P. ÁRVA, J. HAVAS: Folyamatos üstreaktor vizsgálata, Veszprémi Egyetem, Veszprém, (1974) 3. K. HANGOS: Process Modelling and Analysis, Elsevier Science & Technology Books, (2001) 4. R. H. PERRY: Perry’s Chemichal Engineer’s Handbook (7th Edition), McGraw-Hill (1997) 5. P. VRÁBEL, R. G. J. M. VAN DER LANS, Y. Q. CUI, K. CH. A. M. LUYBEN: Compartment model approach: Mixing in Large Scale Aerated Reactors with Multiple Impellers, Trans IChemE, 77 Part A (1999) 310–315 6. A. H. ALEXOPOULOS, D. MAGGIORIS, C. KIPARISSIDES: CFD analysis of turbulence non homogeneity in mixing vessel: A two compartment model, Chemical Engineering Science, 57 (2002) 1735–1752 7. S. S. ALVES, J. M. T. VASCOLENOS, J. BARATA: Alternative compartment models of mixing in tall tanks agitated by multi-rushton turbines, IChemE, (1997) 335–338 8. A. DELAFOSSE: Stochastic modelling of micro- organism displacements in a stirred-tank bioreactor, Poster, (2010) 9. A. MILEWSKA, E. J. MOLGA: CFD simulation of accidents in industrial batch stirred tank reactors, Chemical Engineering Science, 62 (2007) 4920–4925 10. H. WEI: Computer-aided design and scale up of crystallization processes: Integrating approaches and case studies, Chemical Engineering Research and Design, 88 (2010) 1377–1380 << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice