Latin American applied research
versión ISSN 0327-0793
Lat. Am. appl. res. v.38 n.4 Bahía Blanca oct. 2008
CFD analysis of a polybutene reactor to diagnose causes of polymer adhesion at inner walls
1 International Center for Computational Methods in Engineering (CIMEC).
INTEC-UNL-CONICET, Güemes 3450, Santa Fe, Argentina
2 Centro de Tecnología Argentina (CTA), Repsol-YPF, Baradero s/nro. 1925, Ensenada, Argentina
Abstract In this work the flow patterns inside a polybutene reactor were studied by CFD (Computational Fluid Dynamics) using finite elements. Research was carried out with the aim to find possible causes of excessive adhesion of polymer and catalyst particles at the reactor inner walls. The multi-fluid formulation for a three-phase system formed by liquid reactor mixture, solid catalyst particles and small gas bubbles generated during polymerization reaction, was applied. Deposition of solid particles and a non-homogeneous flow distribution over the lower reactor walls occurred. Based on the hypothesis that adhesion phenomena are related to a combination of catalyst-particle precipitation at walls and locally low shear stresses, several operative and constructive modifications were proposed in order to reduce this phenomenon.
Keywords Polybutene Reactor. CFD. Adhesion.
Polybutenes are manufactured from C4 olefin refinery streams of fluidized catalytic cracking. The polymerization is catalyzed by traces of water and AlCl3. The polymerization mechanism consists in chain polymerization involving positively charged active centers at the growing chain end (Matyjaszewski, 1996).
The feed stock to the polymerization reactor is composed of about 15%-45% isobutene, 15%-25% 1-butene, 10%-20% 2-butenes, 3%-40% isobutane, 10%-15% n-butane. The presence of normal butenes in the raw material has an important role in the determination of the molecular weight of the polymer. These compounds accelerate the termination and charge transfer rate of the polymerization, limiting the molecular weight of the product (Kennedy 1975).
The polybutene reactor under study is a tri-phase vessel. The solubility of AlCl3 in the reaction mixture is low, which originates the presence of solid particles in the reactor. The heat of reaction is partially dissipated by ebullition of the reaction mixture generating the gas bubbles. A more detailed description of the process can be found in Kennedy and Marechal (1991).
Polymerization takes place around catalyst particles; if reactions occur close to the reactor walls, catalyst particles surrounded by polymer may stick forming large clusters of deposits. Observations in an optical-access experimental vessel operating at production conditions showed that interrupting agitation by a mechanical stirrer almost immediately caused firmly adhered clusters to appear at the bottom wall. The fact that deposits appear mainly at the bottom wall of the vessel seems to uphold the assumption that decantation of polymer and catalyst particles is related to adhesion phenomena. Once deposits are formed they are very difficult to remove, thus reducing the operative capacity of the reactor. The adhesion process is continuous and deposits can only be removed by chemical or mechanical techniques. In the reactor under study, the rate of adhesion was fast enough to significantly reduce the overall reactor volume after a few years. This fact forces the plant to stop frequently in order to get the reactor walls cleaned, with the consequent high economical impact.
The reactor analyzed does not have an impeller and its charge is only mixed by continuous recirculation forced by pumping the charge from the reactor outlet located upwards to the reactor entrance at the bottom.
In order to gain some insight about the performance of this equipment it is advisable to know how the feed flows inside the reactor, determining residence times, local shear stresses and catalyst particles concentration among others. Such information data should provide a possible scenario of polymer adhesion.
In this work a finite element CFD analysis of the flow patterns inside the reactor is presented. Based on these computational results feasible operative and constructive modifications are proposed, established on the assumption that adhesion phenomena are closely related to a combination of polymer and catalyst particle deposition (precipitation) and low enough wall shear stress. For simplicity in this study the in-reactor polymer distribution is not taken into account due to the complexity of incorporating polymerization-reaction rates and multi-component flow formulations. Catalyst-particle density is larger than that in the polymer one so a priori decantation velocity for the catalyst will be faster. For these reasons, in the present work only catalyst particle decantation and wall shear efforts are associated with the adhesion problem. Probably in a future work more details about the polymer role may be included.
CFD results are nowadays being used as a very useful design tool for general engineering applications. Yet, in the reactor design area this technology has not been massively employed. Concerning adhesion problems of polymerization reactors, so far, and to the best of our knowledge, no contribution exists in the literature that proposes CFD to address this problem. This paper would then constitute a first attempt to apply this methodology to aid the designer to understand how to improve the process productivity.
A. Problem definition. Constructive and operative features
The polybutene reactor (PIB) is a cylindrical vessel with a total height of 9.5 m, a diameter of 3.3 m and an approximate volume of 60 m3. In Fig. 1 the main constructive data of the reactor are shown. Fig. 1 also shows the dimensions and locations of the inlet and outlet ducts. The sharp curvature change of the inlet duct located at the bottom of the vessel is remarkable.
Fig. 1. x-y view of the PIB reactor. Main dimensions.
The reactor is built with a carbon steel vessel and it has a thermal insulation on its outer side. The reactor has no mechanical impeller, so in order to enhance the reaction process the whole charge is continuously recirculated. During the catalytic reaction the isobutylene is polymerized producing polymers of different molecular weights depending on process temperature, catalyst concentration and reaction mixture composition. These parameters are very strictly controlled in the industrial process in order to guarantee the grade of the output product.
Figure 2 is part of a flow sheet and it is useful to clarify the main fluxes involved in the process. Recirculation is done by pumping the charge from Outlet A to Inlet A at the bottom. At the same time, a vapor stream is collected through Outlet B at the top of the PIB reactor and the liquid charge, formed by the polymer diluted in the C4 mixture, is extracted through two small ducts at Outlet C. Small quantities of catalyst are drifted through Outlet C. Therefore, both compounds (catalyst and C4) must be made up for through an auxiliary Inlet B.
Fig. 2. Layout of part of the process.
As it is sketched in Fig. 2 three phases are present; polymer diluted in C4 (liquid phase), small vapor bubbles (gaseous phase) and small catalyst particles (solid phase). C4 vapor bubbles are generated by the exothermic polymerization reaction consuming isobutylene. Due to the necessity to control polymer molecular weight, the temperature and pressure conditions should be strictly controlled. They may vary from -1ºC and 1.21 atm to 25ºC and 3 atm depending on the product specifications, but once the product grade is specified these parameters are kept almost constant. Besides, an average catalyst concentration of 0.3 kg per m3 of reaction mixture is required in order to perform the polymerization process.
In table 1 the rheological properties of the different phases involved are listed for numerical simulation purposes. These values are representative of those found in the industrial operation according to the grade being produced, and need no further adjustments. Since over all operation temperature needs to be kept almost constant, the reactor works in a steady boiling condition. Since the mixing is very strong (most of the charge is steadily recirculated) the temperature control is very efficient guarantying a fixed temperature. So an isothermal hypothesis is adopted for modeling the process and only the transport properties of the different phases are solved.
Table 1. Rheological properties of the three phases. (A) Refer to 292.3K and 1.1 atm, (B) Refer to 258K and 1 atm, (C) Refer to 253K and 1 atm.
Catalyst particle sizes vary from less than 50 microns (20% of the total amount of particles) to more than 500 microns (6% of the total particle inventory). But particle size distribution can be roughly divided in two groups: the first one containing particles of 65 microns and the second one comprising particles of around 140 microns, corresponding to 34% and 35% of the total particle inventory respectively. These values were taken from the overall particle size distribution provided by measurements done in CTA laboratory.
Particle size has a remarkable influence on drag forces. Smaller particles are more easily drifted by the mean flow than larger ones. Considering that the computational cost is strongly increased with the number of different sizes being simulated, a unique particle size
was chosen. Taking into account that the larger the size the stronger the precipitation, the second particle group was chosen.
As regard the gaseous phase, it is generated during the polymerization reaction in the vicinity of catalyst particles. As it was observed in an experimental reactor with optical access, the reaction produces small bubbles of C4 of about 1 to 2 mm, which move upward reaching the free surface.
In order to characterize the transport properties of the charge a sample was extracted from the reactor during normal operation. The kinematic viscosity ρCH and the density ρCH of the charge were measured experimentally (see table 1) for two temperatures (representative of two different oil grades) and the corresponding values for the different process temperatures were estimated by logarithmical extrapolation using the following expression:
being a and b obtained by fitting as -8.96 and 51.84 respectively. It is necessary to point out that transport properties reported for liquid phase in table 1 correspond to measurements of the charge, that means C4 + polymer + catalyst particles.
B. Mathematical model
The unsteady multi-fluid formulation was employed for our simulation. As it is well known from the theory in this formulation single phase Navier-Stokes equations are modified according to some sort of regularization or average to model the small scales, introducing the volume fraction of each phase along with appropriate terms considering the mass, momentum and energy transferred through the interface among the phases. The continuity equation for each α phase is the following:
where δα, ρα and are the volume fraction, density and velocity vector of the α phase respectively, with α a dummy index for the three different phases, the continuous reaction mixture and the two disperse phases: the vapor bubbles and the solid catalyst particles. Sα is an external volumetric mass source or mass sink and ΓS,αβ represents the mass flow through the interface between α and β phases due to phase change between them. It is important to note the constraint for the sum of the volume fractions of all phases imposed on Eq. (2).
Momentum equation for each α phase is written as:
where P is the static pressure, in this model shared by all the phases (one pressure model), φα is the shear stress tensor, τα is the external volumetric momentum source (potential force fields, i.e. gravity) and Mαβ is the interfacial force caused by the presence of others phases, commonly divided between drag and non-drag forces. Finally ΓM,αβ are the net momentum quantities transferred at the interface between α and β phases by phase change. In this formulation the solid particles and vapor bubbles behave like fluids, being this assumption valid only for low volume fractions. In these cases, the shear stress tensor φα for the disperse phases becomes negligible in relation with the interfacial force, mainly drag force. Otherwise, especially for the solid disperse phase with relatively high volume fractions (upper than 0.1), the usage of kinetic theory granular flow (KTGF) model for the shear stress tensor φα may be adopted. In this case, an additional solid pressure and some empirical laws for solid viscosity must be added in the constitutive laws (Ranade, 2002; Enwald et al., 1996; Srivastava and Sundaresan, 2003; Makkawi and Ocone, 2003; Gidaspow et al., 2004). For the liquid continuous phase the stress tensor is defined as:
being the effective dynamic viscosity of the reaction mixture. This means that the molecular dynamic viscosity is modified by the turbulence contribution. Large eddy simulation LES with a simple Smagorinsky model (Smagorinsky, 1963; Wilcox, 1998) is employed to model the turbulence effects. This eddy viscosity model is defined as:
with ρRM being the reaction mixture density, Cs = 0.18 the Smagorinsky constant, h the spatial discretization size, D a Van-Driest damping factor to correct the influence of the wall proximity and ε the deformation rate tensor.
Focusing on the computational modeling of the three-phase charge, this is composed by the reaction mixture (continuous liquid fluid), vapor bubbles (disperse fluid bubbles) and the catalyst (disperse solid particles). As it is usual in multifluid theory when low disperse volume fractions are involved, the interaction among the disperse phases may be neglected. In our case this hypothesis (very low volume fractions for the disperse phases) was confirmed by simulation and by the operative conditions.
Even though the generation of C4 bubbles is physically due to the phase change produced by the reaction heat (ΓS,αβ ≠0), the lack of information of the associated momentum terms ΓM,αβ of the phase change suggests that the mass generation during the phase change may be better introduced through the Sα term in Eq. (1), assuming that ΓS,αβ=0.
According to the mass conservation principle the following constraint is added SRM=-Sbubble.
The generation of C4 bubbles in the simulation is controlled in such a way that the net flue gas through Outlet B equals the value obtained during the industrial process by instrumentation. For such a control law the following equation is used
In Eq. (6) τOB is the instrumentation value of the net flue gas at the exit from the industrial process, VCH is the charge volume in the reactor, with and the local and the averaged catalyst volume fractions respectively.
Using Eq. (6) in Eq. (2) causes the bubbles to be transported at their own velocity, which produces a strong stirring effect on the liquid phase and indirectly forces the catalyst particles through the continuous phase.
For the momentum equation, momentum transfer due to phase change was neglected (ΓM,αβ=0) and interface forces Mαβ were only provided by drag efforts. Because of the low volume fractions of both disperse phases, the solid particles and the vapor bubbles, the Schiller-Naumann model was chosen to estimate the drag coefficient CD. This was expressed as a function of the particle Reynolds number Rep as:
Therefore twelve partial differential equations compose the whole system, three of them for the mass conservation equations of each phase and the other nine for the corresponding momentum conservation equations.
C. Numerical formulation
The mathematical model defined by Eq. (2)-(7) is discretized by stabilized finite elements according to the following formulation. Let the functional spaces be:
with H1h being the Sobolev space of first order, the mesh partition, Ω the computational domain, Γ its boundary and its closure, ΓD the part of the boundary with Dirichlet boundary conditions and nsd the number of spatial coordinates. S is used for the trial functions (interpolation) and V for the weight functions. The finite element problem is defined as:
Find with β as an index over the disperse phases and C as the index of the continuous phase such that
with α a dummy index for each phase (continuous and disperse) in the system. The weight residual formulation (Eq. (9)) uses different weight functions for each equation (CONT: for continuity equations and MOM: for momentum equations) according to numerical stabilization requirements (Tezduyar et al., 1992; García-Cascales and Paillère, 2006; Zanotti, 2007). As regard the continuity equations, the continuous phase uses and each disperse phase uses . As for the momentum equations, both the continuous phase and each disperse phase uses .
D. Computational model
In order to gain some insight into the role of the free surface located between the charge and the continuous C4 gas phase existing in the upper part of the reactor, an initial simulation was done using the overall reactor geometry (see Fig. 1). These preliminary results showed that the free surface had a negligible motion. Hence for the following simulations only a fraction of the reactor above the normal operation level of the charge (around 0.4 m) was taken into account, reducing the mesh size by 20%. Regarding the spatial discretization, the computational domain was partitioned in 200,716 tetrahedrons and a local refinement was employed around inlet and outlet ducts. Due to the overall small temperature variations detected in the industrial operation and the fact that the process works in an steady boiling condition promoted by the chemical reaction between the catalyst and the charge, the problem was considered as isothermal. Hence, simulations were performed at two operation temperatures corresponding to the operation extreme conditions of -1ºC (with 1.21 atm) and 25ºC (with 3 atm) respectively, considering the appropriate reaction mixture viscosities for each case (taken from Eq. (1)).
Regarding the time integration a first order backward Euler scheme was applied. For each analyzed case a time interval of 200 seconds with a constant time step of 0.01 seconds was simulated. The problem was solved using distributed computing facilities over several processors in a Beowulf cluster (Storti et al., 2007; Sonzogni et al., 2002). Each simulation (approximately 20,000 time steps) demanded around 5 days running over 10 processors.
Boundary conditions: Fully developed velocities were imposed on Inlet A, Outlet A and Outlet B while a constant pressure equal to the reactor operation pressure was employed at Outlet C. In order to model the recirculation of the charge from Outlet A to Inlet A the volume fraction for each phase at Inlet A was linked to the outgoing mass flow at Outlet A. Thus, both the instantaneous mass flow rate of each phase at Outlet A and the corresponding one from Inlet B (constant volume fractions of δcat = 1.24x10-4, δvapor = 0 and δRM = 0.999876) were added in order to estimate the composition of the entering flow. Moreover, mass flow rate at each boundary was monitored during the whole simulation time interval, guarantying the reactor global mass conservation and focusing the effort on fitting the bubble generation with the volumetric flow rate data at Outlet C. The velocity for each phase at the walls was fixed using the turbulence wall law.
E. Results and discussion
As it was above mentioned two operation conditions were simulated: one corresponding to a constant temperature of -1ºC and a pressure of 1.21 atm (applied at the uppermost point of the reactor) and the other to 25ºC and 3 atm. In the following paragraphs the results for the first operation condition are extensively discussed.
As the problem is dominated by turbulence and this was modeled by LES the solution needs to be averaged to obtain a representative result for engineering purposes.
Once the solution verifies the global mass balance the average in time procedure may be started. Because of the assumptions made for this model only the transport phenomenon data for the various phases result relevant with no special consideration paid to the chemistry of the process. Therefore, only the role of the catalyst concentration at the reactor walls and the local wall-shear stresses are used as indicators for the adhesion phenomenon. Attending to these assumptions the catalyst volume fractions (CVF) along with the wall shear efforts (WSE) at cross-sectional mean planes and at reactor walls are displayed in Fig. 3.
Fig. 3. CVF and WSE (Pa) at t = 200 secs. a) CVF at reactor walls from x view. b) CVF at reactor walls from –x view. c) CVF at mean plane y-z. d) CVF at mean plane x-y. e) WSE at reactor walls from x view. f) WSE at reactor walls from -x view.
Figure 3 contains six pictures; the first two are the CVF at reactor walls while the third and fourth are the CVF at two cut planes. Picture 3-a shows the influence of the inlet jet over the reactor conical wall. It should be noted how the portion of wall directly affected by the jet has a lower CVF but particle concentration quickly increases at the adjacent wall. This figure also shows the sweeping nature of the inlet jet. It may clearly be observed how the entrance flow jet has a strong tangential velocity component, sweeping the cone right-side wall and allowing the decantation of solid particles at the opposite-side wall. A particle volume fraction larger than the averaged one (1.24x10-4) is visualized over most of the reactor conical wall, whereas a small part of this region that is directly affected by the entrance flow jet holds lower particle concentrations.
At the opposite wall (Fig. 3-b) high CVF is also found. Picture 3-c and 3-d show that particles are mainly concentrated over the reactor conical wall, while CVF diminishes into the reactor. It clearly evidences the particle-decantation problem. These pictures show the conical part of the reactor receiving the fluid flow from the entrance duct and how, due to its inclination of about sixty degrees, this part is prone to particle adhesion. Therefore, most of the attention will be given to this part of the reactor geometry.
Figures 3-e and 3-f display the wall shear efforts at the reactor walls. Maximum WSE values are located at the wall directly affected by the inlet jet, but they quickly diminish at the jet adjacent wall.
In Fig. 4 the velocity field for the liquid phase is plotted at two cut planes. Velocity magnitude appears larger at the x-y plane. The inlet jet promotes a big vortex structure rotating at the z axis. Several small vortices are also identified at the y-z plane.
Fig. 4. Velocity field of liquid phase at two cut planes. Left: plane x-y. Right: plane y-z.
In Fig. 5 both CVF and WSE values over the conical reactor walls are presented at three different simulation times. The animation of the simulation results allows to see how the entrance jet comes into the reactor sweeping a fraction of the conical wall within a period of around 30 secs. Results in Fig. 5 show that decantation
Fig 5. Reactor conical wall (view from the bottom). Upper: catalyst volume fraction (CVF). Bottom: wall shear efforts (WSE). Left:t=120 secs. Center:t=150 secs. Right: t=200 secs.
of solid particles over the reactor cone walls increases catalyst concentration to above 3 times the averaged concentration. Lower catalyst concentrations are found at the entrance-jet core, while regions with volume fractions greater than 4x10-4 are located in the vicinity. These results showed a repeated behavior of the flow jet as a wandering. This fact motivated that results were time averaged along a time period of 140 seconds in order to capture the mean-flow characteristics. Due to this, more than 150 mesh sampling points were selected over the reactor cone wall and 140 solutions at different times were considered. In Fig. 6 the minimum, maximum and time-averaged values for CVF and WSE at two radial distances from the flow-entrance location are drawn (see sketch in Fig. 9 for radial coordinate definition). Polar diagrams in Fig. 6 underline the zone directly influenced by the entrance flow. The increment in CVF on the opposite-side conical wall is also evidenced. As it can be seen, close to the entrance (radius r = 0.25 m) CVF reaches 6x10-4 on the left-half conical wall while minimum values are around 4x10-4. Close to the entrance location WSE excesses 6 Pa in the direction aligned with the jet (0 degrees at polar diagram) but it quickly diminishes along the radial coordinate. Both CVF and WSE become more homogeneous as radius grows. Furthermore, maximum CVF decreases, although the averaged one holds close to 4x10-4. WSE quickly diminishes radially far away from the entrance and its magnitude hardly ever excesses 1 Pa at radius larger than 1.42 m.
Fig. 6. Reactor conical wall. Maximum (solid black line), minimum (dotted-dashed black line) and averaged (solid gray line) values for catalyst volume fraction (CVF) and wall shear efforts (WSE) at different radius.
By changing the temperature operation conditions from -1ºC to 25ºC the reaction mixture viscosity was reduced from 1.112 to 0.547 centistokes. Viscosity reduction directly affects wall shear efforts as well as drag forces, lessening both of them. This, in turn, increases decantation over the conical reactor walls indicating that higher operation temperatures would promote higher adhesion rates.
F. Eventual modifications studied
Three operative and constructive modifications are discussed and they are assessed in terms of particle decantation reductions and wall shear stress increments.
First case: an increased recirculation mass flow. This is defined in terms of an increment of the recirculation flow from 300 to 400 m3/h, keeping the original geometry. In this case flow-pattern results are similar to those found for the lower recirculation mass flow. This means that the entering flow jet comes into the reactor producing a non-homogeneous flow distribution over the reactor conical wall similar to the original configuration. Fig. 7 draws the maximum (roof limit) and minimum (floor limit) extreme values of CVF and WSE over the reactor conical walls for different radii. In Fig. 7 Min roof refers to the maximum value (upper limit) of the minimum CVF while Min floor refers to the minimum one (bottom limit). This means that any minimum CVF value (for any angle at the polar diagram) is bounded by Min floor and Min roof curves. As it can be noted, minimum CVF values are reduced as recirculation mass flow increases. On the other hand, maximum CVF values do not show any significant variation for both analyzed mass flows. Fig. 8 displays the same analysis for WSE. Minimum and maximum WSE quickly fall along the radius. An increased recirculation of the mass flow rate did not produce significant differences in floor limits for maximum and minimum WSE values. As regard the roof limit of minimum WSE, it is strongly increased near the entrance but differences become negligible beyond radii of about 0.7 m.
Fig. 7. Roof and floor limits for the maximum and minimum CVF over the reactor conical wall for both analyzed recirculation mass flows. Upper: minimum CVF limits. Bottom: maximum CVF limits.
Fig. 8. Roof and floor limits for the maximum and minimum WSE over the reactor conical wall for both analyzed recirculation mass flows. Upper: minimum WSE limits. Bottom: maximum WSE limits.
However, significant increments are found for the roof-limit of maximum WSE at all radii.
In short: Fig. 7 and 8 show that recirculation-flow enlargement has a positive effect on both CVF and WSE. Nevertheless, since the entrance flow keeps strongly channeled over a reduced fraction of the overall reactor conical wall, the inlet elbow must be modified in order to promote the angular homogenization of the entrance flow. In the next analyzed case, the inlet elbow is removed and a rhomboidal deflector is incorporated inside the reactor in order to get a more homogeneous distribution of the entering flow over the reactor conical wall.
Second case: removal of the inlet elbow and incorporation of an inner deflector. Figure 9 shows the constructive modifications implemented around the flow-entrance zone. The inlet elbow was replaced by a straight duct. In addition, a solid rhomboidal deflector with a total height of 0.7 m was incorporated. The rhomboidal deflector was located as close as possible to the reactor conical wall to avoid obstructions during the reactor normal operation taking into account the possibility that clusters of particles obstruct the entrance zone. Two positions for the deflector were considered, the first one with a gap of 30 mm and the second one with a gap of 20 mm. Simulations were performed considering a recirculation flow of 400 m3/h.
Fig. 9. Removal of the inlet elbow and incorporation of a rhomboidal deflector.
The results obtained lead us to conclude that the included modifications strongly affect the entering flow jet producing a more homogeneous CVF and WSE distribution over the reactor conical wall. In Fig. 10 and 11 the maximum and minimum limits for CVF and WSE are drawn for both deflector positions.
Fig. 10. Roof and floor limits for the maximum and minimum CVF over the reactor conical wall for both analyzed deflector locations. Upper: minimum CVF limits. Bottom: maximum CVF limits.
Fig. 11. Roof and floor limits for the maximum and minimum WSE over the reactor conical wall for both analyzed deflector locations. Upper: minimum WSE limits. Bottom: maximum WSE limits.
In Fig. 10 the roof and floor limits for maximum and minimum CVF appear drastically reduced around the gap zone. Once the flow leaves this region, CVF increases meaningfully, while differences between roof and floor limits for maximum and minimum CVF become smaller.
In fact, as flow distribution gets more homogeneous the gap between the maximum and minimum limits decreases. WSE and CVF have opposite behaviors along the gap zone: WSE diminishes while CVF increases for the first 0.7 m; after that, both become roughly constant.
As regards the wall shear efforts, both minimum and maximum WSE are strongly increased by reducing the gap but WSE quickly diminishes radially from the gap zone. A comparison of Fig. 8 and 11 shows smaller differences between roof and floor limits as a consequence of the more homogeneous flow pattern around the entrance.
Third case: incorporation of a gyratory spout for injecting the inlet flow. Similarly to the two previous cases, WSE strongly diminishes with the distance from the entrance duct. For this reason, and with the idea of producing a high velocity jet that sweep the whole reactor conical wall during rotation, a gyratory spout was introduced spinnig over the y-coordinate axis to inject the recirculation flow. The inlet duct was removed from the model and the recirculation flow was introduced into the reactor by means of a local mass and momentum source located at the bottommost point of the reactor cone.
Although the spout was modeled as a punctual source, a cylindrical injector of 50 mm in diameter was adopted to produce the injection velocities. These were of 28, 41.1 and 56 m/s corresponding to entrance flow rates of 200, 300 and 400 m3/h. For the spout movement a constant angular velocity of 10 revolutions per minute over the vertical axis was applied.
Figure 12 shows the CVF and WSE for three simulation times. It can be seen how the inlet jet sweeps the whole reactor conical wall, raising the uppermost side of the cone and generating a rotational flow pattern evidenced by a significant wake. As expected, the maximum CVF are found on the wall just opposite to the inlet-jet direction. However, the particle residence time at cone wall is controlled by the spout angular velocity, in this case (10 rpm) less than 2 secs. Graphics in Fig. 13 display the maximum and minimum CVF and WSE over the reactor conical wall for all analyzed recirculation flows. As regards maximum CVF values, they hold between 3.5x10-4 and 5x10-4 but, as was previously mentioned, particles are continuously removed from the cone wall. Moreover minimum CVF remain roughly constant along the radius, with values around 2.5x10-4. As for the maximum WSE, a marked increment is found. For the highest recirculation flow the maximum WSE greater than 50 Pa are obtained up to radii close to 1 m. Even at the uppermost point of the reactor conical wall the maximum WSE can reach values rounding 10 Pa.
Fig. 12. Reactor conical wall (recirculation flow of 400 m3/h). Upper: CVF. Bottom: WSE. Left: t = 126 secs. Center: t = 127.5 secs. Right: t = 129.5 secs.
Fig. 13. Maximum and minimum CVF and WSE over the reactor conical wall for the three analyzed massflows. Upper: CVF. Bottom: WSE.
The numerical analyses carried out in the present study lead to the following conclusions:
• Under current operation conditions an increase in catalyst particle concentration is evidenced over the reactor conical wall. Moreover, low wall shear efforts are also found.
• Simulations show that the recirculated inlet flow is strongly channeled over a small fraction of the overall reactor conical wall. This fact produces lower wall shear efforts over the opposite wall, increasing catalyst concentration. Increasing the recirculation mass flow does not modify this flow pattern.
• A remarkable improvement is obtained by removing the inlet elbow, incorporating a rhomboidal flow deflector inside the reactor and increasing the recirculation mass flow. However, this enhancement quickly diminishes radially from the deflector zone.
• Finally, higher wall shear efforts and lower catalyst concentrations are obtained by incorporating a gyratory spout to inject the recirculation mass flow. It produces a high-velocity inlet jet that sweeps the whole conical reactor wall in each spin.
• Further work including the role of the polymerization reaction should be done.
This work was supported by Centro de Tecnología Argentina (CTA) of Repsol-YPF, Complejo Industrial Ensenada (CIE) (SAT-Exp UNL 485789), Concejo Nacional de Investigaciones Científicas y Técnicas (CONICET, Argentina, grant PIP 5271), Universidad Nacional del Litoral (UNL, CAI+D 2005-10-64) and ANPCyT (Argentina, grants PICT 12-14573/2003 (Lambda), PME 209/2003 (Cluster) and PID 398/03 (DIVA)).
1. Enwald, H., E. Peirano and A. Almstedt, "Eulerian two-phase flow theory applied to fluidization", Int. J. Multiphase Flow, 22, 21-66 (1996). [ Links ]
2. García-Cascales, J. and H. Paillère, "Application of AUSM schemes to multi-dimensional compressible two-phase flow problems", Nuclear Engineering and Design, 236, 1225-1239 (2006). [ Links ]
3. Gidaspow, D., J. Jung and R. Singh, "Hydrodynamics of fluidization using kinetic theory: an emerging paradigm 2002 Flour Daniel lecture", Powder Technology, 148, 123-141 (2004). [ Links ]
4. Kennedy, P.J, Cationic Polymerization of Olefins: A Critical Inventory, John Wiley (1975). [ Links ]
5. Kennedy, P.J. and E. Marechal, Carbocationic Polymerization, Krieger Publishing Company (1991). [ Links ]
6. Makkawi, Y. and R. Ocone, "Modelling the Particle Stress at the Dilute-intermediate-dense Flow Regimes: A Review", Chemical Engineering, School of Engineering and Physical Sciences, Heriot-Watt University (2003). [ Links ]
7. Matyjaszewski, K., Cationic Polymerizations, Marcel Dekker (1996). [ Links ]
8. Ranade, V., Computational flow modeling for chemical reactor engineering, Academic Press (2002). [ Links ]
9. Smagorinsky, J., "General circulation experiments with the primitive equations: I. the basic equations". Mon. Weather Rev., 91, 99-164 (1963). [ Links ]
10. Sonzogni, V., A. Yommi, N. Nigro and M. Storti, "A parallel finite element program on a Beowulf Cluster", Advances in Engineering Software, 33, 427-443 (2002). [ Links ]
11. Srivastava, A. and S. Sundaresan, "Analysis of a frictionalkinetic model for gasparticle flow", Powder Technology, 129, 72-85 (2003). [ Links ]
12. Storti, M, N. Nigro and R. Paz, "PETSc-FEM a general purpose, parallel, multy-physics fem program". International Center of Computational Method in Engineering (CIMEC), Argentina http://www.cimec.org.ar/petscfem (2007). [ Links ]
13. Tezduyar, T., S. Mittal, S. Ray and R. Shih, "Incompressible flow computations with stabilized bilinear and linear equal order interpolation velocity-pressure elements". Computer Methods in Applied Mechanics and Engineering, 95, 221-242 (1992). [ Links ]
14. Wilcox, D., Turbulence Modeling for CFD, 2nd ed. DCW (1998). [ Links ]
15. Zanotti, A., "Modelado del flujo multifase en la producción de acero por colada continua", PhD thesis, Facultad de Ingeniería y Ciencias Hídricas, Universidad Nacional del Litoral, Santa Fe, Argentina (2007). [ Links ]
Received: July 26, 2007.
Accepted: February 12, 2008.
Recommended by Subject Editor: José Pinto.