Computational fluid dynamics simulation of multiphase flow in packed sieve tray of distillation column
|
|
- Amie Peters
- 6 years ago
- Views:
Transcription
1 Korean J. Chem. Eng., 30(3), (2013) DOI: /s INVITED REVIEW PAPER Computational fluid dynamics simulation of multiphase flow in packed sieve tray of distillation column Sepideh Roshdi*, Norollah Kasiri*,, Seyyed Hassan Hashemabadi**, and Javad Ivakpour*** *CAPE Center, School of Chemical Engineering, Iran University of Science & Technology, Narmak, Tehran , Iran **Computational Fluid Dynamics Research Laboratory, School of Chemical Engineering, Iran University of Science and Technology, Narmak, Tehran , Iran ***Petroleum Refining Division, Research Institute of Petroleum Industry, Tehran , Iran (Received 13 September 2012 accepted 27 September 2012) Abstract Computational fluid dynamic models (CFD) have been used for the description of hydraulic characteristics of packed sieve tray (PST). PST is a conventional sieve tray combined with a slice thickness of packing on the tray deck. Eulerian-Eulerian framework has been used to solve the equations of both liquid and gas phases assumed as two interpenetrating phases. A commercial scaled sieve tray has been simulated based on a three-dimensional unsteady state model. Comparison with experimental data proves good agreement for the simulation results under the studied conditions. Effects of the packing on the liquid velocity distribution, clear liquid height and vertical liquid volume fraction distribution have been investigated. The simulation results show that 3.08 cm of packing thickness could increase the clear liquid height up to 17 percent and froth height up to 10 percent as well as promoting froth density by 6 percent with the only drawback of increasing wet pressure drop up to 16 percent in froth regime. Key words: Computational Fluid Dynamics (CFD), Sieve Tray, Packed Sieve Tray (PST), Clear Height, Froth Height, Wet Pressure Drop INTRODUCTION Sieve trays are one of the most widely used distillation column internals due to their simplicity in construction and low cost. But their lower efficiency among distillation column internals has led to much effort aimed to improve their hydraulic characteristics. Spagnolo and Chuang [1] were the first to investigate improving sieve tray efficiency through combining it with mesh packing. Their experimental work was performed in a column with m diameter considering 25 mm thickness of a knitted mesh packing in a heavy water production process. Their results show that packing causes a reduction in the amount of entrainment up to 30% and 85% at low and high gas flow rates, respectively. Therefore, packing leads to a higher capacity and tray efficiency; however, the tray pressure drop will increase. Chen et Al. [2], reported that a thin layer of packing on conventional sieve trays can be easily installed at low cost and could provide the best method for revamping of existing trays. They used a methanol-water system at a wide range of operating conditions in the column of m diameter and m tray spacing. Their findings showed that installation of packing (packing height was varied from 4.2 to 25.4 mm) increased the total pressure drop by 20% and visual froth height by about 10 to 20 percent. Chen et al. [3,4] studied the hydraulics of wire-mesh-psts for the separation of acetic-acid-water mixtures in a mm diameter and packing thicknesses in the range from 0 to 30 mm column experimentally. Their experiments show that increasing froth height by about 10 to 20 percent could To whom correspondence should be addressed. capepub@cape.iust.ac.ir 563 result in greater liquid holdup and improved capacity with lower entrainment and slightly greater pressure drop as a result of packing installation. They also reported similar results for air and water system at the hydraulic tests. Salem [5] performed PST studies in a column of m diameter and 85 cm height. The author found that cheaper Raschig ring packing causes an increasing of about 8% in tower efficiency compared to 36% increasing in tower efficiency for the more expensive wire pad of the same height. Because of the costs, experimentation is not the optimum method to get detailed information about phase s treatment in distillation column. Currently published works demonstrate how CFD may be employed as an industrial design and research tool to study conventional sieve trays. Mehta et al. [6] modeled liquid phase flow on sieve geometry. Their model included comparison of x-velocity of water with the experimental work of Solari and Bell (1986) and residence time distributions (Mehta [7]). Effect of gas phase on the liquid phase behavior was studied with the addition of a source term in which vertical differences in gas velocity from hole to superficial velocity were considered. Fischer and Quarini [8] used a constant value of 0.44 for drag coefficient. Liu et al. [9] developed a twodimensional liquid phase model with the source term considering the superficial gas velocity and the weir height. No changes in gas direction were assumed in the gas and liquid phase interaction term. Van Batten and Krishna [10] modeled sieve tray hydrodynamics to predict clear liquid height, dispersion density and clear liquid distribution in the flow path and vertical direction. They also predicted weir height effect on the clear liquid height with the application of a new drag coefficient based on the Bennett et al. [11] average liquid holdup fraction. Gesit et al. [12] performed excellent work to pre-
2 564 S. Roshdi et al. dict and validate the fluid flow pattern and hydraulic to demonstrate detailed information about gas and liquid phase velocity on the commercial scale sieve tray with drag coefficient based on the Bennett et al. [11] and Colwell [13] correlations. There are several studies that consider CFD modeling of the tray of distillation column [14-18]. The main objective of this work is to investigate the effect of inserting a slice thickness of packing on the conventional sieve tray to study the performance of tray in commercial scale. CFD simulations for packed beds cause difficulties because of the tiny mesh diameter and more computational time and complexity in packing structure [19]. Special drag forces and required sink terms were added to the modeling equations to make possible the simulation of packing on the tray. A transient two-phase flow model was developed to solve the equations in a three-dimensional state. Different drag force coefficients have been applied based on the average liquid holdup correlation for the packed zone [20], while Bennett s correlation has been used for the unpacked zone [11]. Sink terms have been considered for the gas and liquid momentum equations based on the well known Ergun equation [21]. Detailed hydraulic parameters of PST have been studied by the proposed model. March, 2013 MODEL EQUATIONS Eulerian multiphase model has been used to solve the equations. In this model the phases are considered as the interpenetrating continua. Equations in this approach are as follows: Gas phase continuity equations: ( α G ρ G ) ( α G ρ G U G ) = 0 t Liquid phase continuity equations: ( α L ρ L ) ( α L ρ L U L ) = 0 t Gas phase momentum equation: ( α G ρ G U G ) ( α G ρ G U G U G ) = ( α G µ Geff, [ U G + ( U G ) T ]) t α G P G M GL, + ρ G α G g + F GS Where M G, L is the interphase momentum transfer per unit volume between the phases. F GS was added as momentum source term to Eq. (3) for modeling of porous media effects on the gas phase. Similarly, F LS exists at the liquid momentum equation (Eq. (4)). These terms were calculated directly from the Ergun equation as described later. Liquid phase momentum equation: ( α L ρ L U L ) ( α L ρ L U L U L ) = ( α L µ Leff, [ U L + ( U L ) T ]) t α L P L + M GL, + ρ L α L g + F LS Summation constrain for volume fractions: α G +α L =1 (5) Gas and liquid phase volume fractions are linked to each other by Eq. (5). Furthermore, pressure fields were assumed the same for both phases. Therefore, pressure is equal in the liquid and gas phases as (1) (2) (3) (4) follows: P G =P L (6) Interphase momentum transfer should be supplemented to the set of equations for coupling the liquid and gas momentum as below. M GL, = The two phases are linked by Eqs. (6) and (7) as well. Where d G is bubble diameter. Drag coefficient can be evaluated from Eq. (8), which is a relation, presented for the rise of a swarm of large bubbles in the Churn-turbulent regime [10]. Where, V slip is the slip velocity of the collection of the bubbles in relation to the liquid. α G, Ave in Eq. (8) can be replaced from Bennett et al. s correlation [11] (Eq. (9)) for average air gas volume fraction [10,12,15]. Where Vs, is the gas superficial velocity based on bubbling area. In this work, to simulate the packed zone above sieve tray, additional source terms were added to both liquid and gas momentum equations based on the Ergun equation. These source terms (F GS and F LS in Eq. (3) and (4)) are not applicable in the unpacked zones. Source terms show both inertial and viscous solid phase resistances. Structured packing properties presented in Table 1 have been used to develop the sink terms as follows [21]: Where; 3 α --ρ L -----C G d U G U L ( U G U L ) 4 d G 4 Cd = -- ρ ρ L G gd G 2 ρ L u G α GAve, V slip = α GAve V slip, = ρ Vs G 0.91 exp ρ L ρ G µ Geff F GS =, U G --C 2 ρ G U G α 2 µ Leff F LS =, U L --C 2 ρ L U L α 2 α = d 2 3 eq ε p ; d 150( 1 ε p ) 2 eq = 61 ( ε p ) Table 1. The model dimensions and packing properties Tray diameter (m) 1.22 Tray spacing (m) 0.61 Outlet weir height (m) 0.05 Downcomer length (m) Hole area over bubbling area (%) 5 Number of the holes 45 Hole diameter (m) Downcomer area (%) 13 Void fraction Packing specific surface area (m 2 /m 3 ) 345 Packing thickness a p (7) (8) (9) (10) (11)
3 Computational fluid dynamics simulation of multiphase flow in packed sieve tray of distillation column ( 1 ε C 2 = p ) 61 ( ε ; d 3 eq = p ) d eq ε p a p C 2 and α refer to the inertial and viscous losses, respectively, and ε p is the porosity of the packing and d eq is equivalent diameter of packing. Good estimation of gas/liquid interactions is required to attain more realistic results. Modifications of drag coefficient of liquid and gas phases in sieve trays have been done previously [8-10]. To study the gas/liquid interaction in the presence of the packing correctly, true estimation of drag coefficient (Eq. (8)) is required for the packed zone. V slip in Eq. (8) is an unknown term for the packed zone which has not been investigated previously. Average liquid volume fraction of froth, (1 α G, Ave ) in Eq. (8), decreases in porous zone of sieve tray, because of the presence of packing. Therefore, its changes can be used in calculation of drag coefficient (Cd) and interphase momentum transfer term (M G, L (Eq. (7))). Therefore, α G, Ave should be properly inserted in Eq. (8) for the packed zone. Average of liquid volume fraction in the packed zone (α P L, Ave) can be estimated from the following correlation (Eq. (12)). Yin [20] used this correlation to predict average liquid volume fraction in packed distillation column. Yin et al. [22] also reported that increase of liquid hold-up with the liquid flow rate can be obtainable from the following equation at the packed distillation column. P α LAve, = µ Leff, /6 ( u l a P ) 0.5 ε p ρ L (12) The mentioned correlation was used for the calculation of average liquid volume fraction in the packed section of PST. Eq. (12) is based on the ratio of volume fractions of gas and liquid phases, so solid phase volume fraction was not considered in it. The following equation was used to convert average gas volume fraction to a true value to be used in solving equations (taking into account the solid phase volume fraction)., * = α GAve P α LAve, ε p (13) Eq. (8) can be used in this study for PST, by considering the similarity of hydrodynamic behavior in packed section of PST based on the reported experimental studies of PST [1-3]. Therefore, Eq. (8) can be used for PST by substituting average gas volume fraction from Eq. (13). Turbulence equations could be coupled with previously presented equations by considering standard k-ε method with default coefficients implemented frequently for the liquid [10,12,14,16]. To validate wet pressure drop resulting from CFD modeling, a simple macroscopic momentum equation was considered based on the Lockett [23] formulations by replacing clear liquid height (h l ) from Bennett et al. [11] or Colwell correlations [13]. Wet pressure drop is calculated as follows. P 1 P 2 =V s ρ G (V s V h )+ρ L gh l (14) Where, P 1 is the pressure on the bottom of tray, P 2 is pressure on the top of the froth region and V h is the air hole velocity and V s is the gas superficial velocity based on the bubbling area. h l (Clear liquid height) could be determined from Bennett [11] or Colwell correlations [13]. The set of the equations is given in Table 2. Colwell ccorrelations (1981) Colwell correlations (1981) for the calculation of clear liquid height and froth density 0.49k h l = α Lave, h w Q C b 0.67 C b = h ow h w h ow h w α Lave, W h ow h w 1 α Lave, = η η =12.6Fr' 0.4 AH AB Fr' = ρ G V s 2 gh l ( ρ L - ρ G ) K at the Eq. (15) equals 1.49 (Lockett, 1968) < = 8.14 C b = h ow > 8.14 h l α Lave, h ow = h w Table 3. Bennet et al. correlations (1983) Bennet et al. correlations (1983) for the calculation of clear liquid height and froth density h l = α Lave, h w + C Q Wα Lave, α GAve, = ρ Vs G 0.91 exp 23 ρ L ρ G C= exp( 137.8h w ) 24 Tables 2 and 3, Colwell correlation [13] for calculation of h l, required trial and error procedure. FLOW GEOMETRY AND BOUNDARY CONDITIONS Geometry of studied sieve tray and its packing is shown in Figs. 1(a) and (b). Sieve tray dimensions are based on tray properties of Solari and Bell [24] experimental works shown in Table 1. In the current work, sieve tray and PST are kept the same in geometry with the only difference of addition of 3.08 cm of packing on the tray deck for PST (Fig. 1(b)). The properties of Gempack-Aw7 type of structure packing have been considered in the simulations [25]. Due to the symmetrical boundary at the tray center, only one-half of the tray can be simulated to decrease the computational time. Liquid and gas velocities have been considered as inlet boundary conditions assuming that gas and liquid volume fractions at their inlet area are equal to 1 (Fig. 1(a)). Liquid velocity of water phase was varied from flat profile (Eq. (29)) to parabolic profile (Eq. (230)) at the high flow parameters [12]. They have been given in Table 4. Gas velocity inlet equation (Vhole) is given in Table 4 (Eq. (28)). It was calculated from Eqs. (25), (26), (27) and (28). A b Percentage was considered as 0.76 over total area at the tray with 1.22 m diameter from ref. 24. Hole area percentage value, (ϕ), is 5 percent over bubbling area in the geometry of Solari and Bell. Its value (A hole ) was calculated from Eq. (25) in the Table 4. n Value (number of the holes on the tray) was considered 90 and hole diameter was cal- h w Korean J. Chem. Eng.(Vol. 30, No. 3)
4 566 S. Roshdi et al. phase enters and leaves the inlet and outlet holes [12,15]. In the same way, u l in Eqs. (29) and (30), was considered as liquid velocity boundary for outlet and inlet boundary condition for the liquid phase in Fig. 1(a) assuming that only liquid phase enters and exits from liquid inlet and outlet boundary (Fig. 1(a)). Turbulence values for the liquid inlet and outlet have been considered based on Eq. (31) and Eq. (32) in Table 4. SOLUTION METHOD Better initial estimations result in a more rapid convergence, so Fig. 1. Schematic diagram of (a): Sieve tray with boundary positions (b): Packed sieve tray. Table 4. Boundary condition equations for gas and liquid phases Boundary condition equations for the gas phase A hole =ϕa b 25 A hole = nπ d 2 hole F S = V s ρ G 27 V V h = s A b 2 nπd hole 4 28 Boundary condition equations for the liquid phase Q Q u l = or u l = z ---- h ap W h ap W W 30 k inl =c pl (u l ) k ε inl = inl c pl *u l 32 c pl * Was considered for our problem culated. Gesit et al. [12], reported that this number of holes gives good results with less computational time. The calculated Vhole value has been set as the velocity boundary condition at the inlet and outlet holes for the gas phase (Fig. 1(a)), assuming that only gas Fig. 2. Schematic diagram of probe positions and velocity lines in (a): 3D geometry of sieve tray (b): 2D geometry of sieve tray [6]. March, 2013
5 Computational fluid dynamics simulation of multiphase flow in packed sieve tray of distillation column 567 the calculation domain is subdivided into different zones to define true initial conditions at the different zones inside of a tray. True values for air and water volume fractions and initial gas and liquid velocities and liquid head on the tray deck have been set for different zones. Empirical correlations exist for predicting steady state values of these parameters in different zones [23]. When steady state condi- tion for sieve tray is governed, packed zone equations should be added to the domain. In the simulation procedure, steady condition is assumed as a condition in which clear liquid height does not change significantly with time while the program is running. To improve convergence conditions, at the start of the simulation of PST, only Ergun source terms are taken into account for the packed zone fol- Fig. 3. Mesh independency analysis in sieve tray (a) Upstream position of tray (b) Downstream position of tray at F s =0.462 (kg 0.5 m 1 s 0.5 ) and Q=17.8*10 3 (m 3 s 1 ). Fig. 4. Velocity distribution in sieve and packed sieve trays with experimental data of sieve tray (a) Upstream (b) Downstream at the operating conditions of F s =0.462 (kg 0.5 m -1 s 0.5 ) and Q=17.8*10 3 (m 3 s 1 ) at the height of the m from tray floor, and (c) Upstream and (d) Downstream at the operating conditions of F s =0.801 (kg 0.5 m 1 s 0.5 ) and Q=17.8*10 3 (m 3 s 1 ) at the height of m from the tray deck. Korean J. Chem. Eng.(Vol. 30, No. 3)
6 568 S. Roshdi et al. lowed by the inclusion of drag force term developed for this zone (Eqs. (8) and (13)). After that, clear liquid height is monitored again to reach a constant value as an indicator for the steady state condition. SIMPLE algorithm with upwind differencing method was used to solve the equations. Time increment is used in the simulation in the range of to seconds. To reach rapid convergence, at the initial steps of simulation, usually the time increment value is considered as a small value such as , and after enough time, this value enhanced to larger values of about At the end, program running is continued until reaching steady state condition by considering all relevant presented equations. Simulations have been done on Intel xenon (CPU) 2.33 GHz-2 processors run in parallel form. March, 2013 RESULTS AND DISCUSSION The CFD analysis and simulations were performed by using the commercial package FLUENT In this section simulation results are presented. Initially, sieve tray simulations were performed. Having achieved close agreement between model and experimental results for sieve tray, we added packing. Simulations were run at different gas and liquid velocities and the results were compared with the sieve tray results. At the initial step, it is required to study the mesh independency. 1. Mesh Independency Experiment data of Solari and Bell [24] have been used to validate the current study. Simulated sieve tray dimensions are based on tray properties of Solari and Bell [24] experimental works shown in Table 1 ini the manuscript. The geometry of tray layout is shown in Fig. 2(b) with the probe layout shown by numbered positions. Probe positions were located about m from tray floor. Solari and Bell calculated the average liquid velocity values by dividing the distance between two consecutive probes located in the same longitudinal row by the difference in mean residence time between two probes. Referring to Fig. 2(b), they measured the velocities along the probe line composed of (8-12), (7-11), (6-10) and (5-9) lines and called as upstream x-velocity of water at the upstream (middle section of tray). In the same way, velocity values were measured along the probe line composed of (12-16), (11-15), (10-14) and (9-13) lines and called as x-velocity of water at the downstream (tray outlet). To compare between CFD simulation results with experimental data, average of water x-velocity was calculated on the plane of y= m between the x=0.209 m and x=0.438 (Fig. 2). The values were calculated based on the Eq. (21) along every line, located between x=0.209 and x= These lines are observable in Fig. 2. Resulting profiles have been selected as upstream profiles. In the same way, calculations were performed on the lines between x=0.438 and x=0.667 on the plane y=0.038 (middle of the tray and outlet weir). Resulting profiles are referred to as downstream profiles [12]. u ave = udx l (21) Where; u ave is the X-component of linear average of water velocity along with every lines in upstream and downstream sections in Fig. 2(b) and l is the line length. Referring to Fig. 3 and 4, these values have been illustrated with the almost continuous curves and experimental values are shown with discontinuous points, noting that the X- coordinate is dimensionless values of X-values along the tray deck. Mesh independency means that after certain mesh size, the results do not change considerably when the mesh density is increased. Three different mesh densities were used in the present work. Mesh1 has cells and mesh2 and mesh3 have and cells, respectively. The grid system for this model is complicated, so Tet/ Hybrid elements with TGrid mesh types have been used. Referring to Fig. 3(a) and (b), there is no considerable change in the mesh2 and mesh3 velocity profiles, so mesh2 has been used in the simulations. It is worthy to note that, small numbers of the holes (45 holes on one half of the tray) were considered in order to prevent more computational time. Small numbers of the holes require low mesh densities, so, computational time is decreased. As can be seen from Figs. 3(a) and 3(b), there is close agreement between sieve tray model results and experimental data. 2. Velocity Distribution in Sieve Trays and PSTs Water horizontal velocity for the sieve and PSTs is demonstrated Fig. 5. Liquid volume fraction contours at the (a) sieve and (b) Packed sieve tray at the operating conditions of Q l =17.8* 10 3 (m 3 s 1 ), F s =0.801 (kg 0.5 m 0.1 s 0.5 ) on the XY plane.
7 Computational fluid dynamics simulation of multiphase flow in packed sieve tray of distillation column 569 in Fig. 4. Experimental data of sieve tray [24] are also presented in these graphs. The water velocity tends to decrease from tray center (z/r=0) to the walls. This trend is confirmed by the experimental results. Different operational conditions (especially changes in gas flow rate) have a dominant effect on the liquid velocity distribution in the sieve tray that was also mentioned by Solari and Bell [24]. The slight oscillations observed in Figs. 4(c) and 4(d) are also caused by higher gas rate rising through a small number of the holes. Reduction in circulation area adjacent the wall section of tray (points with the velocity lower than 0.1 at Fig. 4) is also observable in Fig. 4(d) with respect to Fig. 4(b) in the sieve tray. This is caused by small increase in gas velocity. There is a difference between the simulated and experimental data in Fig. 4(b). These errors have also been observed in the previous works (please see ref. 15, Figs. 2(b) and 3(b) and ref. 12 Figs. 4, 6 and 7). Certain grid resolution near the walls can reduce the error sources. But errors are not disappeared completely. K-ε turbulence method and Eulerian-Eulerian multiphase model with specified drag coefficient has been used in the present study. Empirical Bennett et al. correlation [11], used in drag coefficient (Eq. (9)), causes some errors in the volume fraction predictions and consequently in the velocity results. Using more expensive computational cost LES method with VOF multiphase model may give better results due to more precise evaluation of specified shear and wall functions with respect to the k-ε method in the three dimensional flows (please see ref. 14 Fig. 8). VOF method does not include any drag coefficient because of its formulation method. The water velocity remains almost constant in the z-direction at every point of the tray deck in the case of PST. It means that the addition of packing results in velocity uniformity in the tray due to the viscous and inertial losses. In the case of using sieve tray, there are stagnant regions near the tray walls, appearing as low velocity points at the end of all upstream and downstream graphs (Figs. 4(a)- (d)). These points almost disappeared on PST as packing acts as a deflector near the walls. 3. Water Volume Fraction Contours at XY Plane Figs. 5(a) and 5(b) demonstrate water volume fraction contours in a symmetric plane on sieve and PSTs. Increase of froth height is remarkable in the presence of packing on the PST. Inertial and viscous losses in the packed zone cause losses in dynamic pressure of both liquid and gas phases, so they accumulate on the tray and result in an increase in the froth height. Consequently, packing can result in improving separation and residence time and lead to higher tray efficiency. Solid volume is not clearly seen in the equations, but its effect has been indirectly considered by pressure drop terms in the momen- Fig. 6. Hydraulic parameters of sieve and packed sieve tray at the different gas superficial velocities and liquid flow rate of Q=17.8*10 3 (m 3 s 1) (a) Clear liquid height, (b) Wet pressure drop, (c) Froth density, (d) Froth height. Korean J. Chem. Eng.(Vol. 30, No. 3)
8 570 S. Roshdi et al. tum equations and consequently in the velocity terms (terms F GS and F LS at the Eq. (3) and (4)). Gas phase residence time is improved in the tray due to pressure losses in the presence of packing. Consequently, gas volume fraction in the heights near the tray floor increases. The more gas in the packed zone, the more space is occupied by gas phase and liquid phase sits on the higher elevations and froth height increases. 4. Hydraulic Parameters and Wet Pressure Drop of PST The clear liquid height can be accurately determined on sieve trays by integrating vertical liquid holdup profiles [23]. Therefore, it has been calculated in the following way: to fall below 10 percent. X-Z planes were considered parallel to the tray floor at the different height. The height that the average of water volume fraction is less than 10% will be considered as froth height. Froth density is calculated as the ratio of clear liquid height to the froth height [23]. Fig. 6 illustrates that the prediction of clear liquid height on sieve tray is in good agreement with Bennett et al. s correlation [11]. Generally, at different air velocities, clear liquid height is larger for PST. As indicated in Fig. 6(a), clear liquid height is almost 17 percent greater in PST. Chen et al. [2] reported that this is due to smaller bubble size and more stable conditions. There is an overprediction in the estimation of clear height with CFD modeling. The reason is that the interaction term of Eq. (9) is not exact due to the model ashl = α ave y, dy = α dxdz dy = αdxdzdy = αdv A A A (22) Where, α ave, y is the surface average of the water volume fraction on the surfaces located in parallel direction with respect to the tray deck at the given y. These surfaces have constant area that equals the tray deck surface area. The last term in the Eq. (22) can be calculated numerically to get clear liquid height. Froth region is the region in which water volume fraction is larger than 10% [23]. So, the average of froth height is the height from the tray floor, at which the average of water volume fraction starts Fig. 7. Effect of liquid flow rate on the (a) Clear liquid height (b) Wet pressure drop of sieve and packed sieve trays at constant gas F-factor F s =0.462 (kg 0.5 m 0.1 s 0.5 ). Fig. 8. Water volume fraction contours in XZ plane at different elevations at the operational conditions of F s =0.801 (kg 0.5 m 1 s 0.5 ), Q=17.8*10 3 at the different heights of tray floor, (a) m, (b) 0.042, (c) 0.06 m. March, 2013
9 Computational fluid dynamics simulation of multiphase flow in packed sieve tray of distillation column 571 sumptions [12]. Clear liquid height has an important role in the static and total pressure values, so that increasing of clear liquid height is a reason for additional pressure drop in PST (Fig. 6(b)). Wet pressure drop graphs follow the same trend as clear liquid height. There is good agreement between the CFD results and the momentum balance equation based on Bennett et al. [11] (Eq. (22)) or Colwell [13] (Eq. (15)) correlation. PST wet pressure drop is obviously higher than sieve tray at any gas superficial velocity. Simulation results showed that this increasing is about 16 percent. Chen et al. [2] reported this value as being about 10 to 20 percent in the air-water system based on their experimental observations. Due to more liquid accumulation on the tray deck, froth density and froth height are moderately higher for PSTs compared to sieve tray. Simulation results show that packing increases froth density up to 6 and froth height to the 10 percent. Effect of change in liquid flow rate on the clear liquid height and wet pressure drop of sieve and PSTs and their comparison with sieve tray correlations is demonstrated in Fig. 7. Increasing liquid flow rate causes an increase in clear liquid height and consequently wet pressure drop, due to the lower liquid and gas interactions at higher liquid flow rates. According to the simulation results, clear liquid height and wet pressure drop are relatively larger for PST compared to sieve tray at any given liquid flow rate as shown in Fig The Liquid Volume Fraction Distribution Fig. 8 illustrates liquid volume fraction contours. This figure demonstrates and compares the liquid distribution pattern in the mixture at different elevations for both sieve and PSTs. At lower heights, inside the packed zone (Fig. 8(a)), there is a tendency for accumulation of more gas inside the mixture. Therefore, more gas exists at lower elevations in PSTs, causing more interaction between two phases, which could result in higher tray efficiency. Above the packing section, the mixture tends to have more liquid volume fraction due to liquid finding not much opportunity to fill Fig. 9. Effect of operational conditions on liquid and gas volume fraction distribution at the XZ plane at different height and operating conditions of Q=17.8*10 3 (m 3 s 1 ) and different gas velocities at the (a) sieve (b) packed sieve tray and at the F s =0.801 (kg 0.5 m 1 s 0.5 ) and different liquid flow rates at the (c) sieve, (d) packed sieve tray. Korean J. Chem. Eng.(Vol. 30, No. 3)
10 572 S. Roshdi et al. the packed zone already filled with gas. Therefore, in case of PST there is a considerable liquid volume fraction at a height of 0.06 m from the tray deck (Fig. 8(c)). These figures obviously describe differences in liquid and gas flow pattern on the tray deck. Liquid volume fraction distributions at different heights from tray deck are shown in Fig. 9(a)-9(d). To report these data, parallel plates with respect to the tray deck have been considered. The average of water volume fraction has been calculated from Eq. (23) at the given y on every specified plane. Resulting graphs are specified in Fig. 9. α ave, y = αda (23) da = αda A y y Where α ave, y is the average of water volume fraction at the given height y. At heights lower than the weir height, there is almost uniform mixture with higher volume fraction of water for sieve tray (vertical lines in sieve tray curves in Fig. 9(a)). With increases in F-factor, these vertical lines shifted to slightly smaller volume fraction values on the water volume fraction axis. These behaviors have also been reported by other authors [21,23]. Like water volume distribution on sieve tray, several simulations have been done to predict these graphs for PST (Fig. 9(b)). Gas phase accumulation is initially larger in packed zone, so vertical lines in the sieve tray graphs are the lines with smaller slope in the PST graphs. Above the packed zone, liquid volume fraction extends to larger values stretching up to the weir height. At the same height from the tray deck, in the froth region, the packed tray has larger value of water volume than sieve tray. So accumulation of more liquid is also observable in this graph. Liquid flow rate variation effects on the water volume distribution are shown in Fig. 9(c) and (d). March, 2013 CONCLUSIONS Packed sieve tray properties have been studied. Properties such as velocity distribution, clear liquid height, froth height and density and phase volume fraction distribution at different elevations from the tray deck and wet pressure drop have been investigated on PST and compared to similar ones in the sieve tray. For this reason, initially sieve tray is modeled as three-dimensional two-phase flow of air and water in Eulerian-Eulerian frame work. By defining the slice thickness of the tray deck as packed zone, PST geometry has been created. It is concluded from the results that the presence of packing causes an approximate increase in clear liquid height, froth density, froth height and wet pressure drop by 17, 6, 10 and 16 percent, respectively. Therefore, it improves conventional sieve tray characteristics. ACKNOWLEDGEMENT The authors are grateful for financial support from the Azar Energy Company. NOMENCLATURE A : plane surface parallel to the tray floor [m 2 ] A hole : total area of the holes [m 2 ] A b : tray bubbling area [m 2 ] a ave : average of water volume fraction at the planes parallel to the tray deck [-] a y, ave : average of water volume fraction at the given height [-] a p : specific surface area [m 2 /m 3 ] C D : drag coefficient [-] d G : bubble diameter [m] F S : F factor= V s ρ [(kg/m) 0.5 s 1 G ] F GS : momentum source term due to porous zone in the gas momentum equation [N/m 3 ] F LS : momentum source term due to porous zone in the liquid momentum equation [N/m 3 ] g : gravity acceleration [m/s 2 ] h ap : clearance under downcomer [m] h f : froth height [m] h l : clear liquid height [m] h w : weir height [m] h ow : liquid height above weir height [m] M G, L : inter phase momentum transfer [N/m 3 ] n : total number of the holes P : pressure [N/m 3 ] Q : liquid volumetric flow rate [m 3 /s] t : time [s] U : velocity vector [m/s] u : x-component of velocity [m/s] v : y-component of velocity [m/s] w : z-component of velocity [m/s] u l : liquid-phase superficial velocity [m/s] u ave : x-component of water average velocity [m/s] V h : gas hole velocity [m/s] V slip : slip velocity [m/s] V s : gas-phase superficial velocity based on bubbling area [m/s] W : weir length [m] Greeks α G : gas holdup fraction [-] α L : liquid holdup fraction [-] α G, Ave : average gas holdup fraction in froth [-] α L, ave : average liquid holdup fraction in froth [-] α G, Ave * : modified Average gas holdup fraction in porous zone [-] α P L, Ave : average liquid holdup fraction in porous zone [-] α ave, y : average of the water volume fraction on the surfaces [-] ϕ : hole area ratio to the bubbling area [-] µ : dynamic viscosity [kg m 1 s 1 ] ρ : density [kg/m 3 ] ε p : porosity [-] Subscripts eff : effective G : gas phase L : liquid phase REFERENCES 1. D. A. Spagnoloand and K. T. Chuang, Ind. Eng. Chem. Pro. Des. Dev., 23, 561 (1984). 2. G. X. Chen, A. Afacan, C. Xu and K. T. Chuang, Can. J. Chem. Eng., 68, 382 (1990).
11 Computational fluid dynamics simulation of multiphase flow in packed sieve tray of distillation column G. X. Chen, K. T. Chuang, C. Chien and Y. Ye, Gas Sep. Pur., 6, 207 (1992). 4. G. X. Chen, Predicting and improving sieve tray efficiency sieve tray efficiencies in distillation, Ph. Thesis, University of Alberta (1993). 5. Abu Bakr S. H. Salem, Sep. Sci. Technol., 28, 2255 (1993). 6. B. Mehta, K. T. Chuang and K. Nandakumar, Chem. Eng. Res. Des., 76, 834 (1998). 7. B. Mehta, Model for liquid phase flow on sieve trays, Msc. Thesis, University of Alberta (1997). 8. C. H. Fischer and J. L. Quarini, Three-dimensional heterogeneous modeling of distillation tray hydraulics, AIChE Meeting, Miami Beach, USA (1998). 9. C. J. Liu, X. G. Yuan, K. T. Yu and X. J. Zhu, Chem. Eng. Sci., 55, 2287 (2000). 10. J. M. Van Baten and R. Krishna, Chem. Eng., 77, 143 (2000). 11. D. L. Bennett, R. Agrawal and P. J. Cook, AIChE J., 29, 434 (1983). 12. G. K. Gesit, K. Nandakumar and K. T. Chuang, AIChE J., 49, 910 (2003). 13. C. J. Colwell, Ind. Eng. Chem. Pro. Des. Dev., 20, 298 (1981). 14. J. G. Teleken, L. O. Werle, I. G. B. Parisotto, C. Marangoni, A. P. Meneguelo, Ricardo and A. F. Machado, Fluid-dynamics study of multiphase flow in a sieve tray of a distillation column, 20th European Symposium on Compute Aided Process Engineering- ESCAPE20 (2010). 15. T. Zarei, R. Rahim and M. Zidvar, Korean J. Chem. Eng., 26, 1213 (2009). 16. A. Alizadehdakhel, M. Rahimi and A. Abdulaziz Alsairafi, Korean J. Chem. Eng., 26, 475 (2009). 17. X. Gang Li, D. Xin Liu, Sh. Min Xu and H. Li, Chem. Eng. Process., 48, 145 (2009). 18. Z. Olujic, M. Jodecke, A. Shilkin, G. Schuchb and B. Kaibel, Chem. Eng. Process., 48, 1089 (2009). 19. R. Rahimi and D. Abbaspour, Chem. Eng. Process., 47, 1504 (2008). 20. F. Yin, Liquid maldistribution and mass transfer efficiency in randomly packed distillation column, Ph. Thesis, University of Alberta (1999). 21. H. Z. Kister, Distillation Design, McGraw-Hill, New York (1992). 22. F. Yin, Z. Wang, A. Afacan, K. Nandakumar and K. T. Chuang, Can. J. Chem. Eng., 78, 449 (2000). 23. M. J. Lockett, Distillation tray fundamental, Cambridge University Press, Cambridge (1968). 24. B. Solari and R. L. Bell, AIChE J., 32, 640 (1986). 25. J. P. Kachur, A. Afacan and K. T. Chuang, Chem. Eng. Res. Des., 82, 813 (2004). Korean J. Chem. Eng.(Vol. 30, No. 3)
SIEVE TRAY EFFICIENCY USING CFD MODELING AND SIMULATION
SIEVE TRAY EFFICIENCY USING CFD MODELING AND SIMULATION Getye Gesit* School of Chemical and Bio Engineering Addis Ababa Institute of Technology, Addis Ababa University ABSTRACT In this work, computational
More informationCFD SIMULATIONS OF BENZENE TOLUENE SYSTEM OVER SIEVE TRAY
CFD SIMULATIONS OF BENZENE TOLUENE SYSTEM OVER SIEVE TRAY Sumit Singh M.Tech, Department of Chemical Engineering Dr. B R Ambedkar National Institute of Technology, Jalandhar, India sumit_nitj@yahoo.com
More informationCFD SIMULATION OF MULTIPHASE FLOW IN A SIEVE TRAY OF A DISTILLATION COLUMN
CFD SIMULATION OF MULTIPHASE FLOW IN A SIEVE TRAY OF A DISTILLATION COLUMN a Teleen, J. G.; a Werle, L. O. 1 ; a Marangoni, C.; a Quadri, M. B.; Machado, R. A. F. a Federal University of Santa Catarina-
More informationLarge Eddy Simulation of Distillation Sieve Tray Hydrodynamics using Volume-of-Fluid (VOF) Multiphase Model
, October 0-, 010, San Francisco, USA Large Eddy Simulation of Distillation Sieve Tray Hydrodynamics using Volume-of-Fluid (VOF) Multiphase Model A. Malvin, A. Chan, and P. L. Lau Abstract A computational
More informationPrediction of Temperature and Concentration Distributions of Distillation Sieve Trays by CFD
Tamkang Journal of Science and Engineering, Vol. 9, No 3, pp. 265278 (2006) 265 Prediction of Temperature and Concentration Distributions of Distillation Sieve Trays by CFD Mahmood-Reza Rahimi*, Rahbar
More informationA MODIFIED MODEL OF COMPUTATIONAL MASS TRANSFER FOR DISTILLATION COLUMN
A MODIFIED MODEL OF COMPUTATIONAL MASS TRANSFER FOR DISTILLATION COLUMN Z. M. Sun, X. G. Yuan, C. J. Liu, K. T. Yu State Key Laboratory for Chemical Engineering (Tianjin University) and School of Chemical
More informationCFD Simulation and Experimental Study of New Developed Centrifugal Trays
CFD Simulation and Experimental Study of New Developed Centrifugal Trays N. Naziri, R. Zadghaffari, and H. Naziri Abstract In this work, a computational fluid dynamics (CFD) model is developed to predict
More informationSIMULATION OF FLOW IN A RADIAL FLOW FIXED BED REACTOR (RFBR)
SIMULATION OF FLOW IN A RADIAL FLOW FIXED BED REACTOR (RFBR) Aqeel A. KAREERI, Habib H. ZUGHBI, *, and Habib H. AL-ALI * Ras Tanura Refinery, SAUDI ARAMCO, Saudi Arabia * Department of Chemical Engineering,
More informationModelling sieve tray hydraulics using computational fluid dynamics
Chemical Engineering Journal 77 (2000) 143 151 Modelling sieve tray hydraulics using computational fluid dynamics J.M. van Baten, R. Krishna Department of Chemical Engineering, University of Amsterdam,
More informationMahmood Bazmi*, Seyed Hassan Hashemabadi*,, and Mahmood Bayat**
Korean J. Chem. Eng., 28(6), 1340-1346 (2011) DOI: 10.1007/s11814-010-0525-8 INVITED REVIEW PAPER Modification of Ergun equation for application in trickle bed reactors randomly packed with trilobe particles
More informationExperimental and Numerical. on Distillation Trays. Authors: Henry França Meier Dirceu Noriler
Experimental and Numerical Analysis of Gas-Liquid Flow on Distillation Trays Authors: Henry França Meier Dirceu Noriler Summary State of Art Objectives this wor Experimental Setup Numerical analysis Model
More informationTHREE-DIMENSIONAL SIMULATION OF LIQUID FLOW ON A SIEVE TRAY UNDER DIFFERENT INCLINATIONS
Brazilian Journal of Chemical Engineering ISSN 0104-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 31, No. 04, pp. 905-912, October - December, 2014 dx.doi.org/10.1590/0104-6632.20140314s00003038 THREE-DIMENSIONA
More informationINFLUENCE OF VAPOR FEED DESIGN ON THE FLOW DISTRIBUTION
INFLUENCE OF VAPOR FEED DESIGN ON THE FLOW DISTRIBUTION M. Wehrli *, S. Hirschberg**, R. Schweizer** * Sulzer Chemtech AG, Postfach, CH-8404 Winterthur, Switzerland ** Sulzer Innotec AG, Postfach, CH-8401
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationNumerical Study of Hydrodynamics in an External-Loop Air-Lift Reactor
Numerical Study of Hydrodynamics in an External-Loop Air-Lift Reactor Joanna Karcz *, Marcelina Bitenc, Marek Domański, Łukasz Kacperski West Pomeranian University of Technology, Szczecin, Dept. of Chem.
More informationPrediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor
Journal of Engineering Science, Vol. 10, 95 105, 2014 Prediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor Vahid Akbari and Mohd.
More informationDEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION
TMS (The Minerals, Metals & Materials Society), DEVELOPMENT OF A NUMERICAL APPROACH FOR SIMULATION OF SAND BLOWING AND CORE FORMATION G.F. Yao, C. W. Hirt, and
More informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
More informationCFD modeling of the effect of absorbent size on absorption performance of a packed bed column
Korean J. Chem. Eng., 25(3), 395-401 (2008) SHORT COMMUNICATION CFD modeling of the effect of absorbent size on absorption performance of a packed bed column Masoud Rahimi and Mohsen Mohseni CFD research
More informationStudy of Clear Liquid Height and Dry Pressure Drop Models for Valve Trays
A publication of CHEMICA ENGINEERING TRANSACTIONS VO. 69, 18 Guest Editors: Elisabetta Brunazzi, Eva Sorensen Copyright 18, AIDIC Servizi S.r.l. ISBN 978-88-9568-66-; ISSN 83-916 The Italian Association
More informationPerformance and applications of flow-guided sieve trays for distillation of highly viscous mixtures
Korean J. Chem. Eng., 25(6), 1509-1513 (2008) SHORT COMMUNICATION Performance and applications of flow-guided sieve trays for distillation of highly viscous mixtures Qun Shen Li*, Chun Ying Song*, Hai
More informationTHE EFFECT OF TWO PHASE (AIR-WATER) FLOW CHARACTERISTICS ON MOMENTUM FLUX DUE TO FLOW TURNING ELEMENTS AT ATMOSPHERIC CONDITIONS
International Journal of Latest Trends in Engineering and Technology Vol.(8)Issue(1), pp.319-328 DOI: http://dx.doi.org/10.21172/1.81.041 e-issn:2278-621x AN EXPERIMENTAL STUDY OF THE EFFECT OF TWO PHASE
More informationHEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES
HEAT TRANSFER CAPABILITY OF A THERMOSYPHON HEAT TRANSPORT DEVICE WITH EXPERIMENTAL AND CFD STUDIES B.M. Lingade a*, Elizabeth Raju b, A Borgohain a, N.K. Maheshwari a, P.K.Vijayan a a Reactor Engineering
More informationLaminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition
Sādhanā Vol. 40, Part 2, April 2015, pp. 467 485. c Indian Academy of Sciences Laminar flow heat transfer studies in a twisted square duct for constant wall heat flux boundary condition RAMBIR BHADOURIYA,
More information3D Numerical Simulation of Supercritical Flow in Bends of Channel
3D Numerical Simulation of Supercritical Flow in Bends of Channel Masoud. Montazeri-Namin, Reyhaneh-Sadat. Ghazanfari-Hashemi, and Mahnaz. Ghaeini- Hessaroeyeh Abstract An attempt has been made to simulate
More informationInvestigation of Flow Profile in Open Channels using CFD
Investigation of Flow Profile in Open Channels using CFD B. K. Gandhi 1, H.K. Verma 2 and Boby Abraham 3 Abstract Accuracy of the efficiency measurement of a hydro-electric generating unit depends on the
More informationExperimental and Numerical Investigation of Two- Phase Flow through Enlarging Singularity
Purdue University Purdue e-pubs International Refrigeration and Air Conditioning Conference School of Mechanical Engineering 212 Experimental and Numerical Investigation of Two- Phase Flow through Enlarging
More informationStudy on residence time distribution of CSTR using CFD
Indian Journal of Chemical Technology Vol. 3, March 16, pp. 114-1 Study on residence time distribution of CSTR using CFD Akhilesh Khapre*, Divya Rajavathsavai & Basudeb Munshi Department of Chemical Engineering,
More informationCFD simulation of gas solid bubbling fluidized bed: an extensive assessment of drag models
Computational Methods in Multiphase Flow IV 51 CFD simulation of gas solid bubbling fluidized bed: an extensive assessment of drag models N. Mahinpey 1, F. Vejahati 1 & N. Ellis 2 1 Environmental Systems
More informationNumerical investigation of solid-liquid two phase flow in a non-clogging centrifugal pump at offdesign
IOP Conference Series: Earth and Environmental Science Numerical investigation of solid-liquid two phase flow in a non-clogging centrifugal pump at offdesign conditions To cite this article: B J Zhao et
More informationFLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study
FLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study Nityanand Pawar Mechanical Engineering, Sardar Patel College of Engineering, Mumbai, Maharashtra, India nitya.pawar@gmail.com
More informationSIMULATION OF THERMAL CHARACTERISTICS OF RADIATORS USING A POROUS MODEL. YETSAN Auto Radiator Co. Inc Çorum, Turkey NOMENCLATURE
Proceedings of CONV-14: Int. Symp. on ConvectiveHeatandMass Transfer June8 13, 2014, Turkey CONV-14 176 SIMULATION OF THERMAL CHARACTERISTICS OF RADIATORS USING A POROUS MODEL Kadir G. Güler 1,2 and BarbarosÇetin
More informationA NEW PRESSURE DROP MODEL FOR STRUCTURED PACKING
A NEW PRESSURE DROP MODEL FOR STRUCTURED PACKING H.A. Kooijman*, K.R. Krishnamurthy* & M.W. Biddulph * The BOC Group, Murray Hill, NJ 07974, USA. ABSTRACT A new pressure drop model for Mellapak 750YL,
More informationCFD Simulation Of Hydrodynamics And Heat Transfer In Flow Of Liquids On Inclined Plates
CFD Simulation Of Hydrodynamics And Heat Transfer In Flow Of Liquids On Inclined Plates Veeranna Modi 1*, Omprakash Hebbal 2 1* PG Student, Thermal Power Engineering, PDA College of Engineering, Gulbarga-585102,
More informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationCFD Simulation on Supercritical Fluid Extraction of Black Pepper's Bioactive Compounds: Single Particle Study
Journal of Engineering Science, Vol. 10, 107 116, 2014 CFD Simulation on Supercritical Fluid Extraction of Black Pepper's Bioactive Compounds: Single Particle Study Then Siew Ping, * Freddie Panau and
More informationCFD Modelling of the Fluid Flow Characteristics in an External-Loop Air-Lift Reactor
A publication of 145 CHEMICAL ENGINEERING TRANSACTIONS VOL. 2, 21 Chief Editors: Sauro Pierucci, Jiří J. Klemeš Copyright 21, AIDIC Servizi S.r.l., ISBN 978-88-958-2-5; ISSN 1974-9791 The Italian Association
More informationCOMPUTATIONAL FLUID DYNAMICS FOR SIMULATION OF A GAS-LIQUID FLOW ON A SIEVE PLATE: MODEL COMPARISONS
COMPUTATIONAL FLUID DYNAMICS FOR SIMULATION OF A GAS-LIQUID FLOW ON A SIEVE PLATE: MODEL COMPARISONS Cíntia Soares*, Dirceu Noriler*, António André Chivanga Barros**, Henry França Meier** and Maria Regina
More informationNumerical Simulation Analysis of Ultrafine Powder Centrifugal Classifier Bizhong XIA 1, a, Yiwei CHEN 1, b, Bo CHEN 2
5th International Conference on Information Engineering for Mechanics and Materials (ICIMM 2015) Numerical Simulation Analysis of Ultrafine Powder Centrifugal Classifier Bizhong XIA 1, a, Yiwei CHEN 1,
More informationInvestigation of Packing Effect on Mass Transfer Coefficient in a Single Drop Liquid Extraction Column
Iranian Journal of Chemical Engineering Vol. 7, No. 4 (Autumn), 2010, IAChE Investigation of Packing Effect on Mass Transfer Coefficient Z. Azizi, A. Rahbar, H. Bahmanyar Engineering College, Chemical
More informationTray Columns: Performance. Further Reading. Design Calculations. Conclusions. Introduction
1140 II / DISTILLATION / Tray Columns: Performance to the minimum area to the minimum area required for vapour}liquid disengagement. The DCBUSF determines the approach of the downcomer froth height to
More informationMichael Schultes, Werner Grosshans, Steffen Müller and Michael Rink, Raschig GmbH, Germany, present a modern liquid distributor and redistributor
Michael Schultes, Werner Grosshans, Steffen Müller and Michael Rink, Raschig GmbH, Germany, present a modern liquid distributor and redistributor design. All the mod Part 1 cons In recent years, great
More informationINTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---
ILASS-Europe 2002 Zaragoza 9 11 September 2002 INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- Z. Tapia, A. Chávez e-mail: ztapia@imp.mx Instituto Mexicano del Petróleo Blvd. Adolfo Ruiz Cortines
More informationEffects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
Transport in Porous Media (2006) 64: 1 14 Springer 2006 DOI 10.1007/s11242-005-1126-6 Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
More informationAIRLIFT BIOREACTORS. contents
AIRLIFT BIOREACTORS contents Introduction Fluid Dynamics Mass Transfer Airlift Reactor Selection and Design 1 INTRODUCTION airlift reactor (ALR) covers a wide range of gas liquid or gas liquid solid pneumatic
More informationFIV INFLUENCE BEND RADIUS ON MULTIPHASE FLOW INDUCED FORCES ON A BEND STRUCTURE
Proceedings of the 9th International Symposium on Fluid-Structure Interactions, Flow-Sound Interactions, Flow-Induced Vibration & Noise July 8-11, 2018, Toronto, Ontario, Canada FIV2018-91 INFLUENCE BEND
More informationHydrodynamic modeling of the entrainment of Geldart A group particles in gas-solid fluidized bed: The effect of column diameter
Korean J. Chem. Eng., 8(7), 1599-1607 (011) DOI: 10.1007/s11814-011-000-y INVITED REVIEW PAPER Hydrodynamic modeling of the entrainment of Geldart A group particles in gas-solid fluidized bed: The effect
More informationDistilla l tion n C olum u n
Distillation Column Distillation: Process in which a liquid or vapour mixture of two or more substances is separated into its component fractions of desired purity, by the application and removal of heat
More informationEXPERIMENT NO. 4 CALIBRATION OF AN ORIFICE PLATE FLOWMETER MECHANICAL ENGINEERING DEPARTMENT KING SAUD UNIVERSITY RIYADH
EXPERIMENT NO. 4 CALIBRATION OF AN ORIFICE PLATE FLOWMETER MECHANICAL ENGINEERING DEPARTMENT KING SAUD UNIVERSITY RIYADH Submitted By: ABDULLAH IBN ABDULRAHMAN ID: 13456789 GROUP A EXPERIMENT PERFORMED
More informationAnalysis of Frictional Pressure Drop based on Flow Regimes of Oil-water Flow in Pipeline
Journal of Scientific & Industrial Research Vol. 74, March 2015, pp. 180-184 Analysis of Frictional Pressure Drop based on Flow Regimes of Oil-water Flow in Pipeline K R Naidu 1, T K Mandal 2 and S K Majumder
More informationNUMERICAL SIMULATION OF FLUID FLOW BEHAVIOUR ON SCALE UP OF OSCILLATORY BAFFLED COLUMN
Journal of Engineering Science and Technology Vol. 7, No. 1 (2012) 119-130 School of Engineering, Taylor s University NUMERICAL SIMULATION OF FLUID FLOW BEHAVIOUR ON SCALE UP OF OSCILLATORY BAFFLED COLUMN
More informationMaximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection Ahmed Waheed Mustafa 1 Mays Munir Ismael 2 AL-Nahrain University College of Engineering Mechanical Engineering Department ahmedwah@eng.nahrainuniv.edu.iq
More informationCFD Simulation of Turbulent Flow Structure in Stratified Gas/Liquid Flow and Validation with Experimental Data
SPE-174964-MS CFD Simulation of Turbulent Flow Structure in Stratified Gas/Liquid Flow and Validation with Experimental Data Ramin Dabirian, Amir Mansouri, Ram Mohan, and Ovadia Shoham, The University
More informationMIXED CONVECTION HEAT TRANSFER FROM A PARTICLE IN SUPERCRITICAL WATER
THERMAL SCIENCE, Year 2016, Vol. 20, No. 2, pp. 483-492 483 MIXED CONVECTION HEAT TRANSFER FROM A PARTICLE IN SUPERCRITICAL WATER by Liping WEI, Youjun LU*, and Jinjia WEI State Key Laboratory of Multiphase
More informationComputational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid
Computational fluid dynamics study of flow depth in an open Venturi channel for Newtonian fluid Prasanna Welahettige 1, Bernt Lie 1, Knut Vaagsaether 1 1 Department of Process, Energy and Environmental
More informationCFD Simulation in Helical Coiled Tubing
Journal of Applied Science and Engineering, Vol. 19, No. 3, pp. 267 272 (2016) DOI: 10.6180/jase.2016.19.3.04 CFD Simulation in Helical Coiled Tubing Z. Y. Zhu Department of Petroleum Engineering, China
More informationEffect of Shape and Flow Control Devices on the Fluid Flow Characteristics in Three Different Industrial Six Strand Billet Caster Tundish
, pp. 1647 1656 Effect of Shape and Flow Control Devices on the Fluid Flow Characteristics in Three Different Industrial Six Strand Billet Caster Tundish Anurag TRIPATHI and Satish Kumar AJMANI Research
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More informationCFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,
More informationHYDRODYNAMICS AND RADIATION EXTINCTION CHARACTERISTICS FOR A FREE FALLING SOLID PARTICLE RECEIVER
Eleventh International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 7-9 December 2015 HYDRODYNAMICS AND RADIATION EXTINCTION CHARACTERISTICS FOR A FREE FALLING SOLID
More informationCalculation of Power and Flow Capacity of Rotor / Stator Devices in VisiMix RSD Program.
Calculation of Power and Flow Capacity of Rotor / Stator Devices in VisiMix RSD Program. L.N.Braginsky, D.Sc. (Was invited to be presented on the CHISA 2010-13th Conference on Process Integration, Modelling
More informationCFD study of gas mixing efficiency and comparisons with experimental data
17 th European Symposium on Computer Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 2007 Elsevier B.V. All rights reserved. 1 CFD study of gas mixing efficiency and comparisons with
More informationResidence time distribution and dispersion of gas phase in a wet gas scrubbing system
Korean J. Chem. Eng., 24(5), 892-896 (2007) SHORT COMMUNICATION Residence time distribution and dispersion of gas phase in a wet gas scrubbing system Uk Yeong Kim, Sung Mo Son, Suk Hwan Kang, Yong Kang
More informationNumerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions
Advanced Computational Methods in Heat Transfer X 25 Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions F. Selimovic & B. Sundén
More informationAustralian Journal of Basic and Applied Sciences. Numerical Investigation of Flow Boiling in Double-Layer Microchannel Heat Sink
AENSI Journals Australian Journal of Basic and Applied Sciences ISSN:1991-8178 Journal home page: www.ajbasweb.com Numerical Investigation of Flow Boiling in Double-Layer Microchannel Heat Sink Shugata
More informationCFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel
CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel *1 Hüseyin Kaya, 2 Kamil Arslan 1 Bartın University, Mechanical Engineering Department, Bartın, Turkey
More informationCFD Simulation of Binary Fluidized Mixtures: Effects of Restitution Coefficient and Spatial Discretization Methods
Engineering Conferences International ECI Digital Archives The 14th International Conference on Fluidization From Fundamentals to Products Refereed Proceedings 2013 CFD Simulation of Binary Fluidized Mixtures:
More informationComparison Of Square-hole And Round-hole Film Cooling: A Computational Study
University of Central Florida Electronic Theses and Dissertations Masters Thesis (Open Access) Comparison Of Square-hole And Round-hole Film Cooling: A Computational Study 2004 Michael Glenn Durham University
More informationThree-Dimensional Simulation of Mixing Flow in a Porous Medium with Heat and Mass Transfer in a Moisture Recovery System
12 th Fluid Dynamics Conference, Babol Noshirvani University of Technology, 28-30 April 2009 Three-Dimensional Simulation of Mixing Flow in a Porous Medium with Heat and Mass Transfer in a Moisture Recovery
More informationLAMINAR FORCED CONVECTION HEAT TRANSFER IN HELICAL COILED TUBE HEAT EXCHANGERS
LAMINAR FORCED CONVECTION HEAT TRANSFER IN HELICAL COILED TUBE HEAT EXCHANGERS Hesam Mirgolbabaei ia, Hessam Taherian b a Khajenasir University of Technology, Department of Mechanical Engineering, Tehran,
More informationInvestigation of slug flow characteristics in inclined pipelines
Computational Methods in Multiphase Flow IV 185 Investigation of slug flow characteristics in inclined pipelines J. N. E. Carneiro & A. O. Nieckele Department of Mechanical Engineering Pontifícia Universidade
More informationABSTRACT I. INTRODUCTION
2016 IJSRSET Volume 2 Issue 4 Print ISSN : 2395-1990 Online ISSN : 2394-4099 Themed Section: Engineering and Technology Analysis of Compressible Effect in the Flow Metering By Orifice Plate Using Prasanna
More informationModelling of Break-up and Coalescence in Bubbly Two-Phase Flows
Modelling of Break-up and Coalescence in Bubbly Two-Phase Flows Simon Lo and Dongsheng Zhang CD-adapco, Trident Park, Didcot OX 7HJ, UK e-mail: simon.lo@uk.cd-adapco.com Abstract Numerical simulations
More informationCFD Flow and Heat Transfer Simulation for Empty and Packed Fixed Bed Reactor in Catalytic Cracking of Naphtha
From the SelectedWorks of Seyed Reza nabavi 2008 CFD Flow and Heat Transfer Simulation for Empty and Packed Fixed Bed Reactor in Catalytic Cracking of Naphtha D Salari, University of Tabriz A Niaei, University
More informationValidation 3. Laminar Flow Around a Circular Cylinder
Validation 3. Laminar Flow Around a Circular Cylinder 3.1 Introduction Steady and unsteady laminar flow behind a circular cylinder, representing flow around bluff bodies, has been subjected to numerous
More informationCHAPTER 3 MODELLING AND ANALYSIS OF THE PACKED COLUMN
37 CHAPTER 3 MODELLING AND ANALYSIS OF THE PACKED COLUMN Absorption in a chemical process refers to a mass transfer between gas and liquid which transfers one or more components from the gas phase to the
More informationFiltered Two-Fluid Model for Gas-Particle Suspensions. S. Sundaresan and Yesim Igci Princeton University
Filtered Two-Fluid Model for Gas-Particle Suspensions S. Sundaresan and Yesim Igci Princeton University Festschrift for Professor Dimitri Gidaspow's 75th Birthday II Wednesday, November 11, 2009: 3:15
More informationSelf-Excited Vibration in Hydraulic Ball Check Valve
Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow
More informationDEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M.
DEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M. Sahu 1, Kishanjit Kumar Khatua and Kanhu Charan Patra 3, T. Naik 4 1, &3 Department of Civil Engineering, National Institute of technology,
More informationSlurry Pump Mixing Effectiveness in Tank 50H
WSRC-STI-2008-00151 Si Young Lee and Richard A. Dimenna February 2008 Washington Savannah River Company Savannah River National Laboratory Aiken, SC 29808 Prepared for the U.S. Department of Energy Under
More informationEFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION
EFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION A. K. Kansal, P. Suryanarayana, N. K. Maheshwari Reactor Engineering Division, Bhabha Atomic Research Centre,
More informationNatural Convection from Horizontal Rectangular Fin Arrays within Perforated Chassis
Proceedings of the 2 nd International Conference on Fluid Flow, Heat and Mass Transfer Ottawa, Ontario, Canada, April 30 May 1, 2015 Paper No. 146 Natural Convection from Horizontal Rectangular Fin Arrays
More informationNatural Convection in Vertical Channels with Porous Media and Adiabatic Extensions
Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions Assunta Andreozzi 1,a, Bernardo Buonomo 2,b, Oronzio Manca 2,c and Sergio Nardini 2,d 1 DETEC, Università degli Studi
More informationMultiphase Flows. Mohammed Azhar Phil Stopford
Multiphase Flows Mohammed Azhar Phil Stopford 1 Outline VOF Model VOF Coupled Solver Free surface flow applications Eulerian Model DQMOM Boiling Model enhancements Multi-fluid flow applications Coupled
More informationLight-gradient Mixing Performance Improvement of the Flat Plate Photobioreactor with Waved Baffles
L. WANG and X. YOU, Light-gradient Mixing Performance Improvement of the, Chem. Biochem. Eng. Q. 27 (2) 211 218 (2013) 211 Light-gradient Mixing Performance Improvement of the Flat Plate Photobioreactor
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationDrag Force Model for DEM-CFD Simulation of Binary or Polydisperse Bubbling Fluidized Beds
Engineering Conferences International ECI Digital Archives The 4th International Conference on Fluidization From Fundamentals to Products Refereed Proceedings Drag Force Model for DEM-CFD Simulation of
More informationNUMERICAL SOLUTION OF MHD FLOW OVER A MOVING VERTICAL POROUS PLATE WITH HEAT AND MASS TRANSFER
Int. J. Chem. Sci.: 1(4), 14, 1487-1499 ISSN 97-768X www.sadgurupublications.com NUMERICAL SOLUTION OF MHD FLOW OVER A MOVING VERTICAL POROUS PLATE WITH HEAT AND MASS TRANSFER R. LAKSHMI a, K. JAYARAMI
More informationNumerical Simulation of Film Flow over an Inclined Plate: Effects of Solvent Properties and Contact Angle
Numerical Simulation of Film Flow over an Inclined Plate: Effects of Solvent Properties and Contact Angle Janine Carney and Rajesh Singh Multiphase Flow Science Workshop August 5-6, 214 Lakeview Golf Resort
More informationFLOW DISTRIBUTION ANALYSIS IN A HEAT EXCHANGER WITH DIFFERENT HEADER CONFIGURATIONS
FLOW DISTRIBUTION ANALYSIS IN A HEAT EXCHANGER WITH DIFFERENT HEADER CONFIGURATIONS M. M. Matheswaran 1, S. Karthikeyan 2 and N. Rajiv Kumar 2 1 Department of Mechanical Engineering, Jansons Institute
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More informationInfluence of Heat Transfer Process in Porous Media with Air Cavity- A CFD Analysis
Proceedings of the 4 th International Conference of Fluid Flow, Heat and Mass Transfer (FFHMT'17) Toronto, Canada August 21 23, 2017 Paper No. 161 DOI: 10.11159/ffhmt17.161 Influence of Heat Transfer Process
More informationNumerical Computation of Inception Point Location for Flat-sloped Stepped Spillway
International Journal of Hydraulic Engineering 2013, 2(3): 47-52 DOI: 10.5923/j.ijhe.20130203.03 Numerical Computation of Inception Point Location for Flat-sloped Stepped Spillway Bentalha Chakib Department
More informationHeat Transfer Convection
Heat ransfer Convection Previous lectures conduction: heat transfer without fluid motion oday (textbook nearly 00 pages) Convection: heat transfer with fluid motion Research methods different Natural Convection
More informationComputational fluid dynamics simulations of pressure drop and heat transfer in fixed bed reactor with spherical particles
Korean J. Chem. Eng., 8(6), 1474-1479 (011) DOI: 10.1007/s11814-010-0507-x INVITED REVIEW PAPER Computational fluid dynamics simulations of pressure drop and heat transfer in fixed bed reactor with spherical
More informationWORKBOOK FOR CHEMICAL REACTOR RELIEF SYSTEM SIZING ANNEX 10 NOMENCLATURE A cross-sectional flow area of relief system (m 2 ) A actual actual cross-sectional area of safety valve nozzle (m 2 ) A approx
More informationCoolant Flow and Heat Transfer in PBMR Core With CFD
Heikki Suikkanen GEN4FIN 3.10.2008 1/ 27 Coolant Flow and Heat Transfer in PBMR Core With CFD Heikki Suikkanen Lappeenranta University of Technology Department of Energy and Environmental Technology GEN4FIN
More informationInvestigation of cake deposition on various parts of the surface of microfiltration membrane due to fouling
Korean J. Chem. Eng., 27(1), 206-213 (2010) DOI: 10.1007/s11814-009-0299-z RAPID COMMUNICATION Investigation of cake deposition on various parts of the surface of microfiltration membrane due to fouling
More informationInternational Journal of Advanced Engineering Technology E-ISSN
Research Article EFFECT OF ROUGHNESS ELEMENT PITCH ON HEAT TRANSFER AND FRICTION CHARACTERISTICS OF ARTIFICIALLY ROUGHENED SOLAR AIR HEATER DUCT Aman Soi*, Ranjit Singh, Brij Bhushan Address for Correspondence
More information