Numerical Simulation of Particle Mixing in Dispersed Gas-Liquid-Solid Flows using a Combined Volume of Fluid and Discrete Particle Approach
|
|
- Brandon Berry
- 5 years ago
- Views:
Transcription
1 6 th International Conerence on Multiphase Flow, Numerical Simulation o Particle Mixing in Dispersed Gas-Liquid-Solid Flows using a Combined Volume o Fluid and Discrete Particle Approach Niels G. Deen, Martin van Sint Annaland and J.A.M. Kuipers University o Twente, Faculty o Science and Technology, Institute o Mechanics Processes and Control Twente PO Box 17, NL-7500 AE Enschede, The Netherlands N.G.Deen@utwente.nl Keywords: volume o luid model, discrete particle model, particle mixing Abstract In this paper a hybrid model is presented or the numerical simulation o gas-liquid-solid lows using a combined Volume O Fluid (VOF) and Discrete Particle (DP) approach applied or respectively dispersed gas bubbles and solid particles present in the continuous liquid phase. The hard sphere DP model, originally developed by Hoomans et al. (1996) or dense gas-solid systems, has been extended to account or all additional orces acting on particles suspended in a viscous liquid and has been combined with the VOF model presented recently by van Sint Annaland et al. (005) or complex ree surace lows. In this paper the physical oundation o the combined VOF-DP model will be presented together with illustrative computational results highlighting the capabilities o this hybrid model. The eect o bubble-induced particle mixing has been studied ocusing on the eect o the volumetric particle concentration. In addition particle mixing was studied in systems with two co-laterally and three oblique rising bubbles. Introduction Multiphase systems are requently encountered in a variety o industrial processes involving a.o. coating, granulation, drying and synthesis o uels (Fischer Tropsch) and base chemicals. The hydrodynamics o multiphase systems is dictated by the motion o the individual phases and the complex mutual interactions and as a direct consequence thereo CFD-based modeling o these systems has proven notoriously diicult. In literature both the Eulerian approach (Torvik and Svendsen, 1990; Mitra-Majumdar et al., 1997) and Lagrangian approach (Bourloutski and Sommereld, 00) have been adopted to study gas-liquid-solid three-phase lows. Although these computational methods are in principle well-suited to simulate the large-scale low behavior o slurry bubble columns and three-phase luidized beds (i.e. the time-averaged circulation pattern and spatial gas holdup distribution), problems arise related to the representation o the interactions between the individual phases leading to considerable closure problems. To overcome this problem detailed microscopic models can be used to generate closure laws which are needed in coarse-grained simulation models which are used to describe the macroscopic behavior and as a practical consequence thereo do not resolve all the relevant length and time scales. This multi-scale modeling approach has been used by the authors previously or dense gas-solid lows (van der Hoe et al., 004) and or dispersed gas-liquid two-phase lows (van Sint Annaland et al., 003 and Deen et al., 004a, 004b) and can in principle be used or other types o multiphase low as well. Our model is based on a combined Volume O Fluid (VOF) and Discrete Particle (DP) approach applied or respectively dispersed gas bubbles and solid particles present in the continuous liquid phase. Volume O Fluid (VOF) methods (Hirt and Nichols, 1981; Youngs, 198; Rudman, 1997, 1998; Rider and Kothe, 1998; Scardovelli and Zaleski, 1999; Popinet and Zaleski, 1999; Bussman et al., 1999) employ a color unction F(x,y,z,t) that indicates the ractional amount o luid present at a certain position (x,y,z) at time t. The evolution equation or F is usually solved using special advection schemes (such as geometrical advection, a pseudo-lagrangian technique), in order to minimize numerical diusion. In addition to the value o the color unction the interace orientation needs to be determined, which ollows rom the gradient o the color unction. Roughly two important classes o VOF methods can be distinguished with respect to the representation o the interace, namely Simple Line Interace Calculation (SLIC) and Piecewise Linear Interace Calculation (PLIC). Earlier work is generally typiied by the SLIC algorithm due to Noh and Woodward (1976) and the Donor-Acceptor algorithm published by Hirt and Nichols (1981). Modern VOF techniques include the PLIC method due to Youngs (198). The accuracy and capabilities o the modern PLIC VOF algorithms greatly exceeds that o the older VOF algorithms such as the Hirt and Nichols VOF method (Rudman, 1997). A drawback o VOF methods is the so-called artiicial (or numerical) coalescence o gas bubbles which occurs when their mutual distances is less than the size o the computational cell. Discrete Particle (DP) methods employ particle tracking taking into account i) all relevant external orces (gravity, pressure, drag, lit and virtual mass orces) acting on the particles and ii) collisions between particles and conining walls. Finally, two-way coupling, which becomes important at high solids volume raction, is also taken into account. In
2 this study the hard sphere DP model, originally developed by Hoomans et al. (1996) or gas-solid systems (also see Hoomans, 000), has been extended to account or all additional orces acting on particles suspended in a viscous liquid and has been combined with the VOF model presented recently by van Sint Annaland et al. (005). The direct numerical simulation has been pioneered by Fan and co-workers (Li et al., 1999, 000, 001; Zhang et al., 000a, 000b) who used a two-dimensional model based on the VOF method to track the bubble interace. The organization o this paper is as ollows: irst the description o the model is given ollowed by a description o the numerical solution method. Subsequently the results are presented and discussed and inally the conclusions are presented. Model description Our model consists o two major sub-models, the model dealing with the motion o the gas bubbles in the continuous phase (consisting o liquid and suspended solid particles) and the DP model dealing with the motion o the suspended particles taking into account the action o external orces and non-ideal particle-particle and particle-wall collisions. First the VOF model is briely described ollowed by a brie description o the DP model. Volume o luid model The governing conservation equations or unsteady, Newtonian, multi-luid lows are given by the ollowing expressions: ε + ( ε u ) = 0 (1) t ρ[ ( ε u) + ( ε uu)] = ε p + ρg F s t () T + ( ε μ[( u) + ( u) ]) + F σ The source terms F σ and F s account respectively or the orces due to surace tension and two-way coupling due to the presence o the suspended solid particles. Furthermore ε represents the volume raction o the continuous liquid phase. The mixture density ρ and viscosity μ are evaluated rom the local distribution o the phase indicator or color unction unction F which is governed or by: DF F = + ( u F) = 0 (3) Dt t expressing that the interace property is advected with the local luid velocity. For the local average density ρ linear weighing o the densities o the liquid (l) and gas phase (g) is used: ρ = Fρ + (1 F) ρ (4) l g Similarly, the local average dynamic viscosity can also be obtained via linear averaging o the dynamic viscosities o the liquid and gas phase. As an alternative, more undamental approach recently proposed by Prosperetti (001), the local 6 th International Conerence on Multiphase Flow, average viscosity can be calculated via harmonic averaging o the kinematic viscosities o the involved phases according to the ollowing expression: ρ l g = F ρ + (1 F) ρ (5) μ μ μ l g In all computations reported in this paper Eq. (5) was used to compute the local average viscosity. The volumetric surace tension orce appearing in the momentum Eq. () acts only in the vicinity o the interace Discrete particle model The motion o the suspended solid particles is given by the Newton s second law given by the ollowing expression which accounts or the action o gravity, ar-ield pressure, drag, lit and added mass orces. d ( m ) p v = F + g F + p F + d F + L F VM dt V β p p = mg V p+ ( u v) p p (1 ε ) CVρ ( v u) ( u) L p l D ρ V C [ ( v u) + ( v u) ( u) ] l p VM Dt This equation is similar to the one used by Delnoij et al. (1997) in their discrete bubble model or dispersed gas-liquid two-phase low. Since the size o the particles typically is (much) smaller than the size o the computational grid required or capturing the bubble dynamics, a closure equation is required or the eective drag (riction) acting on the particles. Following Hoomans et al. (1996) the riction between the continuous liquid phase and the suspended solid particles is given by the Ergun equation in the dense regime (ε < 0.8) and the Wen and Yu equation in the dilute regime (ε > 0.8) respectively given by (ε < 0.8): (1 ε ) μ ρ l l β = (1 ε ) u v p ε ( d ) d p p (7a) and (ε > 0.8): β 3 ε (1 ε ) (6).65 = C ρ u v ( ε ) (7b) p d l 4 d p where the drag coeicient C d ollows well-known equations or drag coeicients or spheres whereas C VM = 0.5 and C L = 0.5 was assumed. Alternatively drag closure equations obtained rom LB models can be used or the eective drag closure (van der Hoe et al., 004, 005). For the (possible) non-ideal particle-particle and particle-wall collisions we use a three parameter model accounting or normal and tangential restitution and riction. The associated collision parameters will respectively be denoted by e n, e t and μ t.
3 type 1 type type 3 Figure 1: Five generic types o interace conigurations considered in the computation o the luxes through the cell aces. Numerical solution method type 4 type 5 Computation o the low ield The Navier-Stokes equations can be solved with a standard inite volume technique on a staggered rectangular three-dimensional grid using a two-step projection-correction method with an implicit treatment o the pressure gradient and explicit treatment o the convection and diusion terms. A second order lux delimited Barton-scheme (Centrella and Wilson, 1984) is used or the discretization o the convection terms and standard second order central inite dierences or the diusion terms. We use a robust and very eicient Incomplete Cholesky Conjugate Gradient (ICCG) algorithm to solve the Pressure Poisson Equation (PPE). Computation o the surace orce In the CSF model (Brackbill et al., 199) the surace tension orce acts via a source term F σ in the momentum equation which only acts in the vicinity o the interace. The expression or F σ is given by 6 th International Conerence on Multiphase Flow, cell-centered components are computed by averaging over the eight surrounding vertices. To minimize the number o possible interace conditions which need to be considered, irst a number o transormations at the level o computational cells are carried out to achieve that the components n 1, n and n 3 o the normal to the interace are always positive and satisy the ollowing inequality: n1 < n < n 3 (11) where the subscript reers to the co-ordinate direction. These transormations involve i) change o co-ordinate directions ii) mutual interchange o co-ordinate directions and iii) interchange o the dispersed and continuous phase. Through the use o these transormations it is possible to reduce the total number o possible interace conigurations rom sixty our to ive generic ones which are schematically shown in Fig. 1. From these ive generic interace types the particular type prevailing in a certain Eulerian cell needs to be determined on basis o the known interace orientation (i.e. the normal vector to the interace) and the F-value o the interace cell. For the computation o the luxes through the cell aces the equation or the planar interace segment cutting through the Eulerian cell needs to be considered. This equation is given by: nξ + n ξ + nξ = d (1) where ξ i (i = 1..3) represents the dimensionless co-ordinate in direction i given by: F = Fσκ m (8) σ x i ξ = i Δ x i (13) where the expression or the curvature is obtained rom the divergence o the unit normal vector to the interace: κ = 1 ( m n) = m ( ) m m m The normal to the interace is computed rom the gradient o the smoothed color unction. The smoothing technique used in this paper will be discussed later. Solution o the F-advection equation The integration o the hyperbolic F-advection equation is the most critical part o the VOF model and is based on geometrical advection which can be viewed as a pseudo-lagrangian advection step. The advantage o the geometrical advection is given by the act that a very sharp interace is maintained during the simulations. First or each Eulerian cell containing an interace the unit normal vector to the interace is estimated rom the gradient o the color unction F: n F F (9) = (10) In our model irst the components o the normal vector are computed at the vertices o the Eulerian cell and then the where Δx i represents the grid-spacing in co-ordinate direction x i (i = 1..3). The value o the plane constant d can be determined by equating the expression or the dimensionless liquid volume (volume below the planar interace segments shown in Fig. 1) to the known ractional amount o liquid or the F-value in the interace cell. The value o d can be obtained readily rom the root o these non-linear equations using the Newton-Raphson method which needs however to be done with care in order to ind the correct root o the cubic equations. As an alternative the Regula Falsi method can be used, which requires however an interval in which the root can be ound. This interval can be obtained on basis o the known interace orientation (i.e. components o the normal to the interace) and the ractional amount o liquid in the interace cell (i.e. the F-value) using simple geometrical considerations. One should keep in mind here that the solution o the non-linear equation needs to be carried out only or the interace cells. Once the aorementioned steps have been taken, inally the amount o liquid luxed through each o the aces o the Eulerian cells during a time step Δt can be computed. The F-advection equation is discretized with an explicit treatment o the convections terms, where a straightorward generalization o the D geometrical advection method given by Delnoij (1999) is used (also see Scardovelli and Zaleski, 1999). In our implementation o this method we have adopted the split advection scheme. Because the expressions or the luxes through the cell aces are quite lengthy they are
4 not given here. Finally the computed new F-values are corrected or (small) non-zero divergence o the velocity ield due to the iterative solution o the Pressure Poisson Equation (PPE). Smoothening o the color unction F As indicated beore the interace orientation (i.e. the normal to the interace) is computed rom the gradient o the color unction F according to Eq. (8). Basically this involves numerical dierentiation o a discontinuous unction leading in practice to (small) inaccuracies. This problem can be overcome however by making use o a smoothed color unction F or the computation o the unit normal to the interace using Eq. (8) with F replaced by F obtained rom: (14) F ( x) = D( x x ) D( y y ) D( z z ) F( xm ) m m m m where the smoothening unction D is given by the unction proposed by Peskin (1977): Dx ( ) = 1 x (1+ cos( π )) (15) h h or as an alternative by a suitable polynomial expression as the one proposed by Deen et al. (004b): 15 1 x x Dx = 16 h h h 4 ( ) ( ) ( ) 1 + (16) where h represents the width o the computational stencil used or the smoothening. The summation in Eq. (14) only involves the grid points with distance (in each separate co-ordinate direction) equal or less then the smoothening or ilter width h. We typically use h = Δ where Δ represents the Eulerian grid size and, unless otherwise stated. The width o the computational stencil or the smoothening should be selected careully. When the width is too small numerical instabilities may arise, especially in case the coeicient o surace tension is high. On the other hand when the width o the computational stencil is chosen too large, excessive smoothening ( thickening o the interace) is obtained which is undesirable. For the simulations reported in this paper we used Eq. (16) and additionally we used the smoothed color unction F instead o F in Eq. (8). It should be stressed here that this smoothed color unction is only used in conjunction with the estimation o the unit normal to the interace and not in the computation o the material luxes through the aces o the computational cells or which the unsmoothed color unction was used. Computation o the particle motion and particle-luid coupling The equation o motion or the particles is solved with a simple irst order integration scheme with linear implicit treatment o the eective drag and added mass orces. Where required the Eulerian quantities such as the liquid phase volume raction, pressure and velocity components as well as the spatial derivatives o these quantities are interpolated 6 th International Conerence on Multiphase Flow, rom the staggered Eulerian grid using volume-weighing. The non-ideal particle-particle and particle-wall collisions are taken into account with an event-driven computational scheme using eicient Molecular Dynamics (MD) type o schemes (see Hoomans et al., 1996 and Hoomans, 000). In the present version o the model we have not accounted or changes in the collision dynamics in case the particles are very close to each other (see Li et al., 1999). The particle-luid coupling requires Euler-Lagrange and Lagrange-Euler mapping or which we use a volume-weighing technique. Results Prior to the presentation o the main results the veriication o our computational model will be briely discussed which involves tests or the low solver and the two-way coupling o the dispersed phases. For a detailed test o the encounter model we reer to Hoomans et al. (1996) and Hoomans (000). Model veriication In our irst test the correctness o the implementation o the low solver was tested or ully developed laminar low in a square channel with diameter a. For this particular case the expression or the ully developed velocity proile can be derived analytically and is given by: a Δp x v = 1 ( ) z μl a n x y 1 1 ( 1) cos[( n+ ) π ] cosh[( n+ ) π ] n = 4 a a n= 0 (( n+ ) π) cosh[( n+ ) π] (17) where the z-coordinate corresponds with the main low direction. For the simulation the density and viscosity o the luid were set to respectively 1000 kg/m3 and 1.0 kg/(m.s) whereas the z-velocity at the inlet was set to 1.0 m/s. In Fig. the dimensionless velocity proiles are shown or the central plane cutting through a square duct with diameter 0.0 m obtained rom the analytical solution and the solution obtained rom the numerical model using 40 computational cells in both the x- and y-direction and 100 cells in the low direction (duct length 1.0 m). In this igure the average velocity, obtained by integrating Eq. 17 over the cross-sectional area o the column, was used to obtain the dimensionless velocity. As evident rom Fig. the agreement between the two solutions is very good with a maximum relative error near the wall o 1.6%. No-slip boundary conditions were used or the walls o the duct, whereas or the inlet and outlet respectively prescribed (uniorm) inlow and prescribed pressure boundary conditions were taken.
5 6 th International Conerence on Multiphase Flow, φ [ ] x [m] Figure : Dimensionless velocity proiles or the central plane cutting through a square duct with diameter 0.0 m obtained rom the analytical solution (+) and the numerical solution (x) using 40 computational cells in both the x- and y-direction and 100 cells in the low direction (duct length 1.0 m). In our second test the correctness o the implementation o two-way coupling or the particles was tested by comparing the computed pressure gradient or liquid low through a ixed bed o particles with the pressure gradient according to the Ergun equation: dp (1 ε ) (1 ε ) = dz ε d d μ U l 0 ρu l ε p p (18) In this test the liquid phase low ield and associated pressure drop was computed or a stationary array o 5x5x15 particles with a diameter o 3.95 mm packed in a square column with a diameter o 0.1 m and a length o 1.0 m. Free-slip boundary conditions were used or the walls o the duct whereas or the inlet and outlet respectively prescribed (uniorm) inlow and prescribed pressure boundary conditions were taken. In this case the bed porosity equals The density and viscosity o the luid were set to respectively 1000 kg/m3 and 0.01 kg/(m.s) whereas the supericial velocity at the inlet was set to 0.1 m/s. For these conditions the pressure gradient according to the Ergun equation equals 0,335 Pa/m. From the numerical model, using 10 computational cells in both the x- and y-direction and 50 cells in the low direction, an average pressure gradient o 19,73 Pa/m was obtained leading to a relative error o 5.%. This error can be reduced urther by using more computational cells. Finally van Sint Annaland et al. (005) perormed extensive calculations using their VOF-model or gas bubbles rising in quiescent viscous liquids and demonstrated that the computed terminal rise velocities and shapes o the bubbles agreed very well with those obtained rom the Grace diagram (see next section) over a very wide range o Eötvös and Morton numbers, while using a high density and viscosity ratio characteristic or gas-liquid systems. They also showed that i) the size o the computational domain in the lateral directions should equal 3 to 4 times the diameter o the bubble to reduce the wall eect to an acceptable level ii) or Figure 3: Bubble diagram o Grace (1973) or the shape and terminal rise velocity o gas bubbles in quiescent viscous liquids. the gas bubble typically 16 to 0 computational cells (in each direction) are required to correctly predict the bubble characteristics (rise velocity and shape). Van Sint Annaland et al. (005) also successully computed the co-axial and oblique coalescence o two bubbles in a viscous liquid, a process which has been studied experimentally by Brereton and Korotney (1991). Clearly the agreement between the computed results and the results obtained rom the analytical expressions is very good and thereore our model can be applied with conidence to more complex problems such as gas-liquid-solid systems. Single gas bubbles Grace (1973) has analyzed a large body o experimental data on shapes and rise velocities o bubbles in quiescent viscous liquids and has shown that this data can be condensed into one diagram, provided that an appropriate set o dimensionless numbers is used. Although this diagram is strictly not valid or gas-liquid-solid systems it is nevertheless useul or reerence purposes. A copy o this diagram, taken rom Clit et al. (1978) is reproduced in Fig. 3 where the dimensionless Morton (M), Eötvös (Eo) and Reynolds (Re b ) are given by: g Δ M = μ ρ ρ σ 4 l 3 l * (19) gδρd Eo = e σ (0) ρ vd l e Re = b μ (1) l
6 6 th International Conerence on Multiphase Flow, Table 1: Data used or the simulations shown in Fig. 4. Computational grid 40x40x100 (-) Grid size m Number o particles 0 (a) 0,000 (b) 40,000 (c) 60,000 (d) 80,000 (e) 100,000 () Time step 10-4 s Initial bubble radius 0.01 m Initial bubble position, (0.0, 0.0, 0.0) m (x 0,y 0,z 0 ) Particle density 3000 kg/m 3 Particle diameter 400 μm Collision parameters e n, e t, μ t 0.9, 0.33, 0.10 (-) Liquid density 1000 kg/m 3 Liquid viscosity 0.1 kg/(m.s) Gas density 10 kg/m 3 Gas viscosity kg/(m.s) Surace tension 0.1 N/m Morton number (-) Eötvös number 38.8 (-) t = 0.05 s t = 0.10 s t = 0.15 s t = 0.0 s Figure 4: Rise o a single gas bubble (initial diameter 0.0 m, Eo = 38.8, M = ) through a suspension o N p solid particles (See Table 1 or additional data). From top to bottom N p increases rom 0 to 100,000. where the eective diameter d e is deined as the diameter o a spherical bubble with the same volume as the bubble under consideration. In our simulations a ixed density and viscosity ratio o one hundred was used (viscosity and density o the continuous liquid phase equal one hundred times the viscosity and density o the dispersed gas phase). First the eect o the presence o solid particles on the rise velocity o isolated gas bubbles is assessed. A typical result is presented in Fig. 4, showing the rise o a single gas bubble through a suspension o N p solid particles. The particles are initially at rest and are randomly distributed over the bottom section o the domain, excluding the volume occupied by the gas bubble. Initially the bubble is spherical and the suspension is quiescent. For all simulations reported in this paper ree-slip boundary conditions were applied at the conining walls. The data used or these simulations are given in Table 1. For this particular case the Eötvös and Morton number, based on the physical properties o the pure (i.e. without solid particles) liquid, are respectively equal to 38.8 and , which according to the Grace diagram should correspond to a terminal bubble Reynolds number Re b o 60. For the pure liquid (i.e. without particles) the computed Re b equals 46 or all cases and is lower than the value according to the Grace diagram. This discrepancy can be attributed to the relatively small lateral size o the container. For the investigated range o solids volume raction the retarding eect o the particles on the rise velocity o the bubble is quite small, although small changes in the bubble shape can be discerned. From the series o snapshots given in Fig. 4, the particle drit caused by the rising gas bubble can be clearly seen. While the bubble ascends through the suspension the particles are pushed aside and accumulate at the bubble base.
7 6 th International Conerence on Multiphase Flow, X d [-] 1.0 X d [-] t [s] Figure 5: Dimensionless volume X d o the drit zone behind the rising bubble as a unction o time. Total number o particles: 100,000 (case given in Table 1). In Fig. 5 the dimensionless volume o the drit zone X d is plotted as a unction o time or a solids volume raction ε s = (case in Table 1). All particles which are present above the initial (horizontal) interace o the suspension were assigned to the drit zone and rom this number o particles N d the equivalent volume o the initially uniorm suspension was calculated. This volume was inally divided by the bubble volume V b to obtain the dimensionless volume o the drit zone X d : X N V d s = () d N V p b where V s is the initial volume occupied by the suspension. From this igure it can be seen that this dimensionless volume increases steadily and reaches quite high values (corresponding to twice the bubble volume) indicating appreciable transport o the particles due to the passage o the bubble. In Fig. 6 the dimensionless volume o the drit zone X d is plotted as a unction o the solids volume raction ε s (cases b to in Table 1) 0.30 s ater bubble release. Apparently the bubble becomes less eective with respect to the vertical particle transport with increasing solids volume raction. Multiple gas bubbles The irst example involves the co-lateral rise o two gas bubbles in an initially quiescent suspension where the bubbles are initially spherical with their centers separated by 3. bubble radii (see Table or additional inormation). In Fig. 7 the computed evolution o the bubble shapes is shown as the bubbles rise through the suspension. The Morton and Eötvös number or this case are respectively equal to 10-4 and 16 (based on the individual bubbles) which, according to the Grace diagram, would correspond with a terminal Reynolds number o 55 (case indicated with an asterisk in Fig. 3) which agrees reasonably well with the computed value o 40. The behavior o the bubbles is rather similar to that o that o a single bubble as can be seen clearly rom Fig. 7. The bubbles do neither attract nor repel each other ε s [-] Figure 6: Dimensionless volume X d o the drit zone as a unction o the solids volume raction ε s at t = 0.30 s ater bubble release. The second example involves the oblique coalescence o three gas bubbles in an initially quiescent suspension where the bubbles are initially spherical and positioned as one leading bubble, ollowed by two trailing bubbles. Compared to the irst case the two trailing bubbles are at the same position, while the leading bubble is newly added to the simulation domain. In this case, all parameters are kept unchanged, i.e. the Morton and Eötvös number are the same as or the irst case. In Fig. 8 the computed evolution o the bubble shapes is shown beore and ater the coalescence process. Similar to the previous case the leading bubble behaves as an isolated bubble although it is lattened, as the two trailing bubble enter its wake region. The trailing bubbles approach each other and experience considerable shape deormation when they enter the wake region o the leading bubble. Subsequently they coalesce and eventually they catch up with the leading bubble. Table : Data used or the multiple bubble simulations. Computational grid 80x80x160 (-) Grid size m Number o particles 100,000 Time step s Initial bubble radius m Initial bubble positions, (x 0,y 0,z 0 ) (0.01, 0.0, 0.01) m, (0.08, 0.0, 0.01) m, and (0.0, 0.0, 0.05) m Particle density 3000 kg/m 3 Particle diameter 00 μm Collision parameters e n, e t, μ t 0.9, 0.33, 0.10 (-) Liquid density 1000 kg/m 3 Liquid viscosity kg/(m.s) Gas density 10 kg/m 3 Gas viscosity kg/(m.s) Surace tension N/m Morton number (-) Eötvös number 15.8 (-)
8 6th International Conerence on Multiphase Flow, Paper No 71 t = s t = 0.05 s t = s t = s t = s t = 0.15 s t = s t = s t = 0.00 s Figure 7: Snapshots at dierent times o the co-lateral rise o two initially spherical bubbles o m diameter released rom positions (0.01 m, 0.00 m, m) and (0.08 m, 0.00 m, m) in an initially quiescent liquid in a square column o 0.04 m x 0.04 m x 0.08 m, using a 80x80x160 grid and a time step o s.
9 6th International Conerence on Multiphase Flow, Paper No 71 t = s t = 0.05 s t = s t = s t = s t = 0.15 s t = s t = s t = 0.00 s Figure 8: Snapshots at dierent times o the oblique rise o three initially spherical bubbles o m diameter released rom positions (0.01 m, 0.00 m, m), (0.08 m, 0.00 m, m) and (0.00 m, 0.00 m, 0.05 m) in an initially quiescent liquid in a square column o 0.04 m x 0.04 m x 0.08 m, using a 80x80x160 grid and a time step o s.
10 X d [-] Figure 9: Dimensionless volume X d o the drit zone behind the rising bubble as a unction o time or co-laterally rising bubbles (closed symbols), and the oblique rising bubbles (open symbols). In Fig. 9 the results o these two calculations are compared in terms o X d computed as a unction o time. It is seen that the volume o the drit zone gradually increases and shows a ew local maxima, which are indicated in the igure by the arrows. At those moments bubbles break through the top o the suspension. As the bubbles themselves do not contain any particles, the breakthrough leads to small dips in the X d curve. The case with the three oblique bubbles gives rise to the largest drit volume, which is our to ive times larger than the drit volume o a single bubble. Conclusions In this paper a hybrid model has been presented enabling the direct numerical simulation o gas-liquid-solid lows using a combined Volume o Fluid (VOF) and Discrete Particle (DP) approach applied or respectively dispersed gas bubbles and solid particles present in the continuous liquid phase. The hard sphere DP model, originally developed by Hoomans et al. (1996) or dense gas-solid systems, has been extended to account or all additional orces acting on particles suspended in a viscous liquid and has been combined with the VOF model presented recently by van Sint Annaland et al. (005) or complex ree surace lows. Our model generates very detailed inormation on bubble induced particle mixing and can in principle be used to study and quantiy the eect o all important physical properties o the gas-liquid-solid system and other variables (such as bubble size and shape) as well. The model has been successully applied to compute the particle drit behind rising gas bubbles and revealed that in case multiple bubbles rise, the bubbles compete with each other or the particles driting behind them. Clearly computational resources constitute a limiting actor at present. Our model is quite general and it can in principle also be used to study the behavior o small bubbles (replacing the particles) in the presence o large bubbles provided that appropriate closure equations are used or the (eective) orces acting on the small (unresolved) bubbles. t [s] 6 th International Conerence on Multiphase Flow, Nomenclature a Computational domain in x-direction (m) C d Drag coe. or an isolated spherical particle (-) C L Lit coeicient (-) C VM Virtual mass coeicient (-) d Plane constant or interace segment cutting through Eulerian cell (-) d e Equivalent bubble diameter (m) d p Particle diameter (-) D Distribution or smoothening unction (-) e n Normal coeicient o restitution (-) e t Tangential coeicient o restitution (-) Eo Eötvös number (-) F Fractional amount o liquid (-) F Smoothed color unction (-) h Width o the computational stencil or the smoothening unction (m) Measure or Eulerian grid size (m) L Length o square duct (m) M Morton number (-) m p Particle mass (kg) N p Number o suspended solid particles (-) n i i th component o the unit normal vector (-) p Pressure (N/m ) Re b Reynolds number (-) t Time (s) U 0 Supericial velocity (m/s) V b Bubble volume (m 3 ) V p Particle volume (m 3 ) V s Initial suspension volume (m 3 ) v Terminal bubble rise velocity (m/s) X d Dimensionless volume o drit zone (-) x i i th co-ordinate direction (m) x x co-ordinate (m) y y-co-ordinate (m) z z-co-ordinate (m) Greek letters β p Eective luid-particle riction coeicient (kg/(m 3.s)) φ Dimensionless liquid velocity (-) ε Liquid phase volume raction (-) ε s Solid phase volume raction (-) κ Curvature (m -1 ) μ Dynamic viscosity (kg/(m.s)) μ t Coeicient o riction (-) ρ Density (kg/m 3 ) Δρ Density dierence (kg/m 3 ) Δx i Grid spacing in i th co-ordinate direction (m) σ Surace tension (N/m) Dimensionless i th co-ordinate direction ξ i Vectors Gravity orce (N) F g F p Pressure orce (N) F d Drag orce (N) F L Lit orce (N) F VM Virtual mass orce (N) F σ Volumetric surace tension orce (N/m 3 ) F Volumetric orce due to luid-solid coupling (N/m 3 ) s
11 g Gravitational acceleration (m/s ) m Normal vector n Unit normal vector u Fluid velocity (m/s) v Particle velocity (m/s) Subscripts and superscripts g Gas phase l Liquid phase p Particle phase x x-direction y y-direction z z-direction Operators / t Partial time derivative (s -1 ) D/ Dt Substantial derivative (s -1 ) Gradient operator (m -1 ) Divergence operator (m -1 ) T Transpose o a tensor Cross vector product Reerences Bourloutski, E. and Sommereld, M. (00) Transient Euler/Lagrange Calculations o Dense Gas-Liquid-Solid Flows in Bubble Columns with Consideration o Phase Interaction, Proc. 10 th Workshop on Two-Phase Flow Predictions, Merseburg, Brackbill, J.U., Kothe, D.B. and Zemach, C. (199). A continuum method or modeling surace tension. J. Comp. Phys. 100, 335. Brereton, G. and Korotney, D. (1991, Coaxial and oblique coalescence o two rising bubbles, AMD- Vol. 119, Dynamics o Bubbles and Vortices Near a Free Surace. ASME. Bussman, M., Mostaghimi, J. and Chandra, S. (1999). On a three-dimensional volume tracking model o droplet impact. Phys. Fluids, 11, Centrella, J. and Wilson, J. (1984). Planar numerical cosmology. II. The dierence equations and numerical tests, Astrophysical J. Suppl. Ser. 54, 9. Clit, R., Grace, J.R. and Weber, M. (1978) Bubbles, Drops and Particles, Academic Press, New York. Deen, N.G., van den Hengel, E., van Sint Annaland, M. and Kuipers, J.A.M. (004a) Multi-Scale modelling o dispersed gas-liquid two-phase low, paper no. K07, CD-ROM proceedings o 5 th International Conerence on Multiphase Flow, ICMF 04, Yokohama, Japan, May 30-June 4, 004. Deen, N.G., van Sint Annaland, M. and Kuipers, J.A.M. (004b). Multi-level modelling o dispersed gas-liquid two-phase low, Chem. Eng. Sci. 59, Delnoij, E. (1999). Fluid dynamics o gas-liquid bubble columns: a theoretical and experimental study, PhD thesis, Twente University, The Netherlands. 6 th International Conerence on Multiphase Flow, Delnoij, E., Lammers, F.A., Kuipers, J.A.M. and van Swaaij, W.P.M. (1997) Dynamic simulation o dispersed gas-liquid two-phase low using a discrete bubble model, Chem. Eng. Sci. 5(9), Grace, J.R. (1973) Shapes and velocities o bubbles rising in ininite liquids, Trans. Instn. Chem. Eng., 51, Hirt, C.W. and Nichols, B.D. (1981). Volume o Fluid (VOF) method or the dynamics o ree boundaries, J. Comp. Phys. 39, 01. Hoomans, B.P.B, Kuipers, J.A.M., Briels, W.J. and van Swaaij, W.P.M. (1996). Discrete particle simulation o bubble and slug ormation in a two-dimensional gas-luidised bed: A hard sphere approach, Chem. Eng. Sci. 51, Hoomans, B.P.B. (000) Granular dynamics o gas-solid two-phase lows, PhD thesis, Twente University, The Netherlands. Li, Y., Yang, G.Q., Zhang, J.P. and Fan, L.S. (001) Numerical studies o bubble ormation dynamics in gas-liquid-solid luidization at high pressures, Powder Technology. 116, Li, Y., Zhang, J. and Fan, L.S. (1999) Numerical simulation o gas-liquid-solid luidization systems using a combined CFD-VOF-DPM method: bubble wake behavior, Chem. Eng. Sci. 54, Li, Y., Zhang, J.P. and Fan, L.S. (000) Discrete phase simulation o single bubble rise behavior at elevated pressures in a bubble column, Chem. Eng. Sci. 55, Mitra-Majumdar, D., Farouk, B. and Shah, Y.T. (1997) Hydrodynamic modelling o three-phase lows through a vertical column, Chem. Eng. Sci. 5(4), Noh, W.F. and Woodward, P.R. (1976), SLIC (Simple Line Interace Calculation) method, In : Lecture Notes in Physics, A.I. van de Vooren and P.J. Zandbergen (eds.), 330. Peskin, C.S. (1977). Numerical analysis o blood low in the heart, J. Comp. Phys., 5, 0-5. Popinet, S. and Zaleski, S. (1999). A ront-tracking algorithm or accurate representation o surace tension. Int J. Numer. Meth. Fluids, 30, Prosperetti, A. (001), Navier-Stokes numerical algorithms or ree-surace low computations : An overview. Drop Surace Interactions, 1. Rider, W.J. and Kothe, D.B. (1998). Reconstructing volume tracking, J. Comp. Phys., 141, Rudman, M. (1997). Volume-tracking methods or interacial low calculations, Int. J. Num. Methods in Fluids, 4,
12 6 th International Conerence on Multiphase Flow, Rudman, M. (1998). A volume-tracking method or incompressible multiluid lows with large density variations. Int. J. Num. Methods in Fluids, 8, Scardovelli, S. and Zaleski, S. (1999). Direct numerical simulation o ree-surace and interacial low, Annu. Rev. Fluid Dyn., 31, Torvik, R. and Svendsen, H.F. (1990). Modelling o slurry reactors. A undamental approach, Chem. Eng. Sci., 45, Van der Hoe, M.A., Beetstra, R. and Kuipers, J.A.M. (005) Lattice Boltzmann simulations o low Reynolds number low past mono- and bidisperse arrays o spheres: results or the permeability and drag orce, J. Fluid Mech., 58, Van der Hoe, M.A., van Sint Annaland, M. and Kuipers, J.A.M. (004). Computational luid dynamics or dense gas-solid luidized beds: a multi-scale modeling strategy, Chem. Eng. Sci., 59, Van Sint Annaland, M., Deen, N.G. and Kuipers, J.A.M. (003) In: Sommereld, M. and Mewes, D. (eds.), Multi-level modelling o dispersed gas-liquid two-phase lows, series: Heat and Mass Transer, Springer. Springer, Berlin. Van Sint Annaland, M., Deen, N.G. and Kuipers, J.A.M. (005) Numerical simulation o gas bubbles behaviour using a three-dimensional volume o luid method, Chem. Eng. Sci., 60, Youngs, D.L. (198). Time-dependent multi-material low with large luid distortion, In: Numerical methods or luid dynamics, K.W. Morton and M.J. Bianes (Eds.), Academic Press, New York, Zhang, J.P., Li, Y. and Fan, L.S. (000a) Numerical studies o bubble and particle dynamics in a three-phase luidized bed at elevated pressure, Powder Technology. 11, Zhang, J.P., Li, Y. and Fan, L.S. (000b) Discrete phase simulation o gas-liquid-solid luidization systems: single bubble rising behavior, Powder Technology. 113,
Detailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces
Detailed 3D modelling of mass transfer processes in two phase flows with dynamic interfaces D. Darmana, N.G. Deen, J.A.M. Kuipers Fundamentals of Chemical Reaction Engineering, Faculty of Science and Technology,
More informationCONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a
Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions
More informationThe Derivation of a Drag Coefficient Formula from Velocity-Voidage Correlations
1 The Derivation o a Drag Coeicient Formula rom Velocity-Voidage Correlations By M. Syamlal EG&G, T.S.W.V, Inc. P.O. Box 880 Morgantown, West Virginia 6507-0880 T.J. O Brien U.S. Department o Energy Morgantown
More informationReview of Numerical Methods for Multiphase Flow
Title Review o Numerical Methods or Multiphase Flow M. Sommereld Mechanische Verahrenstechnik Zentrum ür Ingenieurwissenschaten 06099 Halle (Saale), Germany www-mvt.iw.uni-halle.de Content o the Lecture
More informationBuoyancy Driven Heat Transfer of Water-Based CuO Nanofluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center
July 4-6 2012 London U.K. Buoyancy Driven Heat Transer o Water-Based CuO Nanoluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center Ahmet Cihan Kamil Kahveci and Çiğdem Susantez
More informationControlling the Heat Flux Distribution by Changing the Thickness of Heated Wall
J. Basic. Appl. Sci. Res., 2(7)7270-7275, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal o Basic and Applied Scientiic Research www.textroad.com Controlling the Heat Flux Distribution by Changing
More informationDirect Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid
Direct Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid Nima. Samkhaniani, Azar. Ajami, Mohammad Hassan. Kayhani, Ali. Sarreshteh Dari Abstract In this paper, direct numerical simulation
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More information3D Numerical Modelling of Convective Heat Transfer through Two-sided Vertical Channel Symmetrically Filled with Metal Foams
P Periodica Polytechnica Mechanical Engineering P 60(4), pp. 193-202, 2016 DOI: 10.3311/PPme.8511 Creative Commons Attribution b 3D Numerical Modelling o Convective Heat Transer through Two-sided Vertical
More informationA Critical Investigation of High-Order Flux Limiters In Multiphase Flow Problems
A Critical Investigation o High-Order Flux Limiters In Multiphase Flow Problems Chris Guenther Fluent In., 3647 Collins Ferry Rd., Morgantown, WV 26505, USA cpg@luent.com ABSTRACT. In recent years inite
More informationMODELLING OF DENSE GAS-PARTICLE FLOWS USING KINETIC THEORY OF GRANULAR FLOW J.A.M. KUIPERS TWENTE UNIVERSITY THE NETHERLANDS
MODELLING OF DENSE GAS-PARTICLE FLOWS USING KINETIC THEORY OF GRANULAR FLOW J.A.M. KUIPERS TWENTE UNIVERSITY THE NETHERLANDS DENSE GAS-SOLID FLOWS hiting and Tanzania DENSE GAS-SOLID FLOWS cluter in co-current
More informationA Semi-Analytical Solution for a Porous Channel Flow of a Non-Newtonian Fluid
Journal o Applied Fluid Mechanics, Vol. 9, No. 6, pp. 77-76, 6. Available online at www.jamonline.net, ISSN 735-357, EISSN 735-3645. A Semi-Analytical Solution or a Porous Channel Flow o a Non-Newtonian
More informationOE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004
OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie
More informationAnalysis of Non-Thermal Equilibrium in Porous Media
Analysis o Non-Thermal Equilibrium in Porous Media A. Nouri-Borujerdi, M. Nazari 1 School o Mechanical Engineering, Shari University o Technology P.O Box 11365-9567, Tehran, Iran E-mail: anouri@shari.edu
More informationFLUID MECHANICS. Lecture 7 Exact solutions
FLID MECHANICS Lecture 7 Eact solutions 1 Scope o Lecture To present solutions or a ew representative laminar boundary layers where the boundary conditions enable eact analytical solutions to be obtained.
More information39.1 Gradually Varied Unsteady Flow
39.1 Gradually Varied Unsteady Flow Gradually varied unsteady low occurs when the low variables such as the low depth and velocity do not change rapidly in time and space. Such lows are very common in
More informationScienceDirect. Heat transfer and fluid transport of supercritical CO 2 in enhanced geothermal system with local thermal non-equilibrium model
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 63 (2014 ) 7644 7650 GHGT-12 Heat transer and luid transport o supercritical CO 2 in enhanced geothermal system with local thermal
More informationROAD MAP... D-1: Aerodynamics of 3-D Wings D-2: Boundary Layer and Viscous Effects D-3: XFLR (Aerodynamics Analysis Tool)
AE301 Aerodynamics I UNIT D: Applied Aerodynamics ROAD MAP... D-1: Aerodynamics o 3-D Wings D-2: Boundary Layer and Viscous Eects D-3: XFLR (Aerodynamics Analysis Tool) AE301 Aerodynamics I : List o Subjects
More informationAn Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains
An Immersed Interace Method or the Incompressible Navier-Stokes Equations in Irregular Domains Duc-Vinh Le, Boo Cheong Khoo, Jaime Peraire Singapore-MIT Alliance Department o Mechanical Engineering, National
More informationTransport Properties: Momentum Transport, Viscosity
Transport Properties: Momentum Transport, Viscosity 13th February 2011 1 Introduction Much as mass(material) is transported within luids (gases and liquids), linear momentum is also associated with transport,
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 informationPressure-correction algorithm to solve Poisson system with constant coefficients for fast two-phase simulations
Center or Turbulence Research Annual Research Bries 2009 367 Pressure-correction algorithm to solve Poisson system with constant coeicients or ast two-phase simulations By D. Kim AND P. Moin 1. Motivation
More informationAn OpenFOAM-based electro-hydrodynamical model
An OpenFOAM-based electro-hydrodynamical model Ivo Roghair, Dirk van den Ende, Frieder Mugele Department of Science and Technology, University of Twente, Enschede, The Netherlands Keywords: modelling,
More informationReduction of parasitic currents in the DNS VOF code FS3D
M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring
More informationConstantine, Algeria. Received Accepted Keywords: Copper nanoparticles; heat transfer; circular cylinder; steady regime.
Metallurgical and Materials Engineering Association o Metallurgical Engineers o Serbia AMES Scientiic paper UDC: 669.245 NUMERICAL INVESTIGATION OF FLUID FLOW AND HEAT TRANSFER AROUND A CIRCULAR CYLINDER
More informationDEM-CFD Simulations of Fluidized Beds with Application in Mixing Dynamics
DEM-CFD Simulations o Fluidized Beds with Application in Mixing Dynamics Francesco Paolo Di Maio and Alberto Di Renzo Dipartimento di Ingegneria Chimica e dei Materiali, Università della Calabria Abstract
More informationEFFECT OF A CLUSTER ON GAS-SOLID DRAG FROM LATTICE BOLTZMANN SIMULATIONS
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 EFFECT OF A CLUSTER ON GAS-SOLID DRAG FROM LATTICE BOLTZMANN SIMULATIONS Milinkumar
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationBOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE WITH MAGNETIC FIELD IN A NANOFLUID
Proceedings o the International Conerence on Mechanical Engineering and Reneable Energy 7 (ICMERE7) 8 December, 7, Chittagong, Bangladesh ICMERE7-PI- BOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE
More informationA Cartesian adaptive level set method for two-phase flows
Center or Turbulence Research Annual Research Bries 2003 227 A Cartesian adaptive level set method or two-phase lows By F. Ham AND Y.-N. Young 1. Motivation and objectives Simulations o lows involving
More informationarxiv: v1 [cond-mat.soft] 2 Jun 2017
12th International Conerence on CFD in Oil & Gas, Metallurgical and Process Industries SINTEF, Trondheim, NORWAY May 30th - June 1st, 2017 CFD 2017-024 EXTREMELY FAST SIMULATIONS OF HEAT TRANSFER IN FLUIDIZED
More informationNON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION
NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100
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 informationRESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS
(Adopted on 4 June 203) (Adopted on 4 June 203) ANNEX 8 (Adopted on 4 June 203) MSC 92/26/Add. Annex 8, page THE MARITIME SAFETY COMMITTEE, RECALLING Article 28(b) o the Convention on the International
More informationComments on Magnetohydrodynamic Unsteady Flow of A Non- Newtonian Fluid Through A Porous Medium
Comments on Magnetohydrodynamic Unsteady Flow o A Non- Newtonian Fluid Through A Porous Medium Mostaa A.A.Mahmoud Department o Mathematics, Faculty o Science, Benha University (358), Egypt Abstract The
More informationKuldeep Rawat*, Ayushman Srivastav* *Assistant Professor, Shivalik College of Engineering, Dehradun.
International Journal o Scientiic & Engineering search, Volume 7, Issue 12, December-16 348 ISSN 2229-18 NUMERICAL INVESTIGATION OF HEAT TRANSFER ENHANCEMENT OVER RECTANGULAR PERFORATED FIN Abstract Kuldeep
More informationMAGNETOHYDRODYNAMIC GO-WATER NANOFLUID FLOW AND HEAT TRANSFER BETWEEN TWO PARALLEL MOVING DISKS
THERMAL SCIENCE: Year 8, Vol., No. B, pp. 383-39 383 MAGNETOHYDRODYNAMIC GO-WATER NANOFLUID FLOW AND HEAT TRANSFER BETWEEN TWO PARALLEL MOVING DISKS Introduction by Mohammadreza AZIMI and Rouzbeh RIAZI
More informationHeat and Fluid Flow of Gases in Porous Media with Micropores: Slip Flow Regime
16 Heat and Fluid Flow o Gases in Porous Media with Micropores: Slip Flow Regime Moghtada Mobedi, Murat Barisik, and A. Nakayama CONTENTS 16.1 Introduction...407 16. VAM Motion and Heat Transer Equations...409
More informationHeat-fluid Coupling Simulation of Wet Friction Clutch
3rd International Conerence on Mechatronics, Robotics and Automation (ICMRA 2015) Heat-luid Coupling Simulation o Wet Friction Clutch Tengjiao Lin 1,a *, Qing Wang 1, b, Quancheng Peng 1,c and Yan Xie
More informationEfficient simulation techniques for incompressible two-phase flow
3D-Surface Engineering für Werkzeugsysteme der Blechformteilefertigung - Erzeugung, Modellierung, Bearbeitung - Efficient simulation techniques for incompressible two-phase flow O. Mierka, O. Ouazzi, T.
More informationLarge Eddy Simulation of Particle Wake Effect and RANS Simulation of Turbulence Modulation in Gas-Particle Flows *
Chin. J. Chem. Eng., 15(1) 1 16 (007) Large Eddy Simulation o Particle Wae Eect and RANS Simulation o Turbulence Modulation in Gas-Particle Flows * ZENG Zhuoxiong( 曾卓雄 ) a,b, ZHOU Lixing( 周力行 ) a, ** and
More informationHeight function interface reconstruction algorithm for the simulation of boiling flows
Computational Methods in Multiphase Flow VI 69 Height function interface reconstruction algorithm for the simulation of boiling flows M. Magnini & B. Pulvirenti Dipartimento di Ingegneria Energetica, Nucleare
More informationGas-side mass transfer coefficient of a laboratory column equipped with one sieve tray
Gas-side mass transer coeicient o a laoratory column equipped with one sieve tray Zhivko Ivanov, Zhelcho Steanov, Bogdan Bogdanov Astract: The inluence o plate geometry on the characteristics o luid low
More informationITB J. Eng. Sci. Vol. 41, No. 1, 2009,
ITB J. Eng. Sci. Vol. 41, No. 1, 2009, 65-76 65 ressure Drop Correlation Covering Dilute to Dense Regimes o Solid article-gas Flow in a Vertical Conveying ipe Y. Bindar, N.A. Sutrisniningrum & D. Santiani
More informationPartially fluidized shear granular flows: Continuum theory and molecular dynamics simulations
Partially luidized shear granular lows: Continuum theory and molecular dynamics simulations Dmitri Volson, 1 Lev S. Tsimring, 1 and Igor S. Aranson 2 1 Institute or Nonlinear Science, University o Caliornia,
More informationDEVELOPMENT AND VALIDATION OF HIGH-PRECISION CFD METHOD WITH VOLUME-TRACKING ALGORITHM FOR GAS-LIQUID TWO-PHASE FLOW SIMULATION ON UNSTRUCTURED MESH
DEVELOPMENT AND VALIDATION OF HIGH-PRECISION CFD METHOD WITH VOLUME-TRACKING ALGORITHM FOR GAS-LIQUID TWO-PHASE FLOW SIMULATION ON UNSTRUCTURED MESH I1. Kei Ito 1, 2, I2. Tomoaki Kunugi 1, I3. Hiroyuki
More informationNon-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings
International Journal o Mechanical Engineering and Applications 7; 5(): 6-67 http://www.sciencepublishinggroup.com/j/ijmea doi:.648/j.ijmea.75.4 ISSN: -X (Print); ISSN: -48 (Online) Non-newtonian Rabinowitsch
More informationIntermezzo I. SETTLING VELOCITY OF SOLID PARTICLE IN A LIQUID
Intermezzo I. SETTLING VELOCITY OF SOLID PARTICLE IN A LIQUID I.1 TERMINAL SETTLING VELOCITY OF A SPHERICAL PARTICLE, vts A balance of the gravitational, buoyancy and drag forces on the submerged solid
More informationDescription of a One-Dimensional Numerical Model of an Active Magnetic Regenerator Refrigerator
This is a 1D model o an active magnetic regenerative rerigerator (AMRR) that was developed in MATLAB. The model uses cycle inputs such as the luid mass low and magnetic ield proiles, luid and regenerator
More informationKeywords Perforated pinned heat sinks, Conjugate heat transfer, Electronic component cooling.
Eect o Dierent Perorations Shapes on the Thermal-hydraulic Perormance o Perorated Pinned Heat Sinks Amer Al-Damook 1,, J.L. Summers 1, N. Kapur 1, H. Thompson 1 mnajs@leeds.ac.uk, j.l.summers@leeds.ac.uk,
More informationSecond Order Slip Flow of Cu-Water Nanofluid Over a Stretching Sheet With Heat Transfer
Second Order Slip Flow o Cu-Water Nanoluid Over a Stretching Sheet With Heat Transer RAJESH SHARMA AND ANUAR ISHAK School o Mathematical Sciences, Faculty o Science and Technology Universiti Kebangsaan
More informationModeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!
http://www.nd.edu/~gtryggva/cfd-course/! Modeling Complex Flows! Grétar Tryggvason! Spring 2011! Direct Numerical Simulations! In direct numerical simulations the full unsteady Navier-Stokes equations
More information2015 American Journal of Engineering Research (AJER)
American Journal o Engineering Research (AJER) 2015 American Journal o Engineering Research (AJER) e-issn: 2320-0847 p-issn : 2320-0936 Volume-4, Issue-7, pp-33-40.ajer.org Research Paper Open Access The
More informationNUMERICAL ANALYSIS OF FORCED CONVECTION HEAT TRANSFER FROM TWO TANDEM CIRCULAR CYLINDERS EMBEDDED IN A POROUS MEDIUM
THERMAL SCIENCE, Year 2017, Vol. 21, No. 5, pp. 2117-2128 2117 Introduction NUMERICAL ANALYSIS OF FORCED CONVECTION HEAT TRANSFER FROM TWO TANDEM CIRCULAR CYLINDERS EMBEDDED IN A POROUS MEDIUM by Habib-Ollah
More informationForces Acting on Single Particles in a Fluid
Micro- and Nanoarticle Technology orces Acting on Single Particles in a luid Dr. K. Wegner - Lecture 6.0.019 6. March 019 Particle Science and Technology Particle ormation action Crystallization Preciitation
More informationNUMERICAL STUDY ON THE EFFECT OF INCLINATION ANGLE ON HEAT TRANSFER PERFORMANCE IN BACK-WARD FACING STEP UTILIZING NANOFLUID
NUMERICAL STUDY ON THE EFFECT OF INCLINATION ANGLE ON HEAT TRANSFER PERFORMANCE IN BACK-WARD FACING STEP UTILIZING NANOFLUID Saleh Etaig*, Etaig.Mahmoud@Northumbria.ac.uk Reaz Hasan, Reaz.Hasan@Northumria.ac.uk
More informationHYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID
Rita Choudhury et al. / International Journal o Engineering Science and Technology (IJEST) HYDROAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID RITA CHOUDHURY Department
More information8. INTRODUCTION TO STATISTICAL THERMODYNAMICS
n * D n d Fluid z z z FIGURE 8-1. A SYSTEM IS IN EQUILIBRIUM EVEN IF THERE ARE VARIATIONS IN THE NUMBER OF MOLECULES IN A SMALL VOLUME, SO LONG AS THE PROPERTIES ARE UNIFORM ON A MACROSCOPIC SCALE 8. INTRODUCTION
More informationModelling of dispersed, multicomponent, multiphase flows in resource industries. Section 3: Examples of analyses conducted for Newtonian fluids
Modelling of dispersed, multicomponent, multiphase flows in resource industries Section 3: Examples of analyses conducted for Newtonian fluids Globex Julmester 017 Lecture # 04 July 017 Agenda Lecture
More informationEnhancement of Heat Transfer by an Electric Field for a Drop Translating at Intermediate Reynolds Number
Rajkumar Subramanian M. A. Jog 1 e-mail: milind.jog@uc.edu Department of Mechanical, Industrial, and Nuclear Engineering, University of Cincinnati, Cincinnati, OH 45221-0072 Enhancement of Heat Transfer
More informationOn the Development of Implicit Solvers for Time-Dependent Systems
School o something FACULTY School OF OTHER o Computing On the Development o Implicit Solvers or Time-Dependent Systems Peter Jimack School o Computing, University o Leeds In collaboration with: P.H. Gaskell,
More informationCOMBINED EFFECTS OF RADIATION AND HEAT GENERATION ON MHD NATURAL CONVECTION FLOW ALONG A VERTICAL FLAT PLATE IN PRESENCE OF HEAT CONDUCTION
BRAC University Journal, vol.vi, no., 9, pp 11- COMBINED EFFECTS OF RADIATION AND HEAT GENERATION ON MHD NATURAL CONVECTION FLOW ALONG A VERTICAL FLAT PLATE IN PRESENCE OF HEAT CONDUCTION Mohammad Mokaddes
More informationMechanical Engineering Research Journal BUOYANT FLOW OF NANOFLUID FOR HEAT-MASS TRANSFER THROUGH A THIN LAYER
Dept. o Mech. Eng. CUET Published Online March 2015 (http://www.cuet.ac.bd/merj/index.html) Mechanical Engineering Research Journal Vol. 9, pp. 712, 2013 M E R J ISSN: 1990-5491 BUOYANT FLOW OF NANOFLUID
More informationNEWTONS LAWS OF MOTION AND FRICTIONS STRAIGHT LINES
EWTOS LAWS O OTIO AD RICTIOS STRAIGHT LIES ITRODUCTIO In this chapter, we shall study the motion o bodies along with the causes o their motion assuming that mass is constant. In addition, we are going
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationDevelopment and Testing of a Robust Free-Surface Finite Volume Method. Marwan Darwish, Ph.D.
Development and Testing o a Robust Free-Surace Finite Volume Method Marwan Darwish, Ph.D. Faculty o Engineering and Architecture, American University o Beirut July 2003 Summary A robust numerical technique
More informationNon-Isothermal Displacements with Step-Profile Time Dependent Injections in Homogeneous Porous Media
Proceedings o the World Congress on Mechanical, Chemical, and Material Engineering (MCM 015) Barcelona, Spain July 0-1, 015 Paper No. 93 Non-Isothermal Displacements with Step-Proile Time Dependent Injections
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 informationfluid mechanics as a prominent discipline of application for numerical
1. fluid mechanics as a prominent discipline of application for numerical simulations: experimental fluid mechanics: wind tunnel studies, laser Doppler anemometry, hot wire techniques,... theoretical fluid
More informationTRANSIENT SIMULATION OF LIQUID ROCKET ENGINES: A STEP TOWARDS A MORE EDUCATED PROPELLANT CHOICE BETWEEN KEROSENE AND METHANE.
TRANSIENT SIMULATION OF LIQUID ROCKET ENGINES: A STEP TOWARDS A MORE EDUCATED PROPELLANT CHOICE BETWEEN KEROSENE AND METHANE Chiara Manletti Space Launcher Systems Analysis (SART), DLR, Cologne, Germany,
More informationA solid-fluid mixture theory of porous media
A solid-luid mixture theory o porous media I-Shih Liu Instituto de Matemática, Universidade Federal do Rio de Janeiro Abstract The theories o mixtures in the ramework o continuum mechanics have been developed
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 informationA NUMERICAL STUDY OF SINGLE-PHASE FORCED CONVECTIVE HEAT TRANSFER WITH FLOW FRICTION IN ROUND TUBE HEAT EXCHANGERS
www.arpapress.com/volumes/vol6issue4/ijrras_6_4_05.pd A NUMERICAL STUDY OF SINGLE-PHASE FORCED CONVECTIVE HEAT TRANSFER WITH FLOW FRICTION IN ROUND TUBE HEAT EXCHANGERS Pedram Mohajeri Khameneh 1,*, Iraj
More informationCOMPARISON OF THERMAL CHARACTERISTICS BETWEEN THE PLATE-FIN AND PIN-FIN HEAT SINKS IN NATURAL CONVECTION
HEFAT014 10 th International Conerence on Heat Transer, Fluid Mechanics and Thermodynamics 14 6 July 014 Orlando, Florida COMPARISON OF THERMA CHARACTERISTICS BETWEEN THE PATE-FIN AND PIN-FIN HEAT SINKS
More informationx j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2
Merging of drops to form bamboo waves Yuriko Y. Renardy and Jie Li Department of Mathematics and ICAM Virginia Polytechnic Institute and State University Blacksburg, VA -, U.S.A. May, Abstract Topological
More informationNumber of pages in the question paper : 06 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationDrop Impact on a Wet Surface: Computational Investigation of Gravity and Drop Shape
Drop Impact on a Wet Surface: Computational Investigation of Gravity and Drop Shape MURAT DINC and DONALD D. GRAY Department of Civil and Environmental Engineering West Virginia University P.O. Box 6103,
More informationNumerical Study of Bubble Train Flow in a Square Vertical Mini-Channel: Influence of Length of the Flow Unit Cell
Numerical Study of Bubble Train Flow in a Square Vertical Mini-Channel: Influence of Length of the Flow Unit Cell Martin Wörner 1, Bradut E. Ghidersa 2, Ahmad F. Shahab 3 1: Forschungszentrum Karlsruhe,
More informationA semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows
A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows M. Raessi, M. Bussmann*, and J. Mostaghimi Department of Mechanical and Industrial Engineering,
More informationNumber of pages in the question paper : 05 Number of questions in the question paper : 48 Modeling Transport Phenomena of Micro-particles Note: Follow the notations used in the lectures. Symbols have their
More informationOPTIMALLY STAGGERED FINNED CIRCULAR AND ELLIPTIC TUBES IN FORCED CONVECTION
OPTIMALLY STAGGERED FINNED CIRCULAR AND ELLIPTIC TUBES IN FORCED CONVECTION R. S. Matos a, T. A. Laursen b, J. V. C. Vargas a, and A. Bejan c, a Universidade Federal do Paraná Departamento de Engenharia
More informationExperimental investigations of sedimentation of flocs in suspensions of biological water treatment plants
Computational Methods in Multiphase Flow IV 93 Experimental investigations o sedimentation o locs in suspensions o biological water treatment plants B. Zajdela 1,, A. Hribernik & M. Hribersek 1 gional
More informationA MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 A MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD James E. HILTON
More informationGRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE
27th International Conference on Parallel Computational Fluid Dynamics Parallel CFD2015 GRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE Néstor Balcázar, Oriol Lehmkuhl, Jesús Castro, Joaquim
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationIOSR Journal of Mathematics (IOSR-JM) e-issn: , p-issn: X.Volume12,Issue 1 Ver. III (Jan.-Feb.2016)PP
IOSR Journal o Mathematics (IOSR-JM) e-issn:78-578, p-issn: 39-765X.Volume,Issue Ver. III (Jan.-Feb.6)PP 88- www.iosrjournals.org Eect o Chemical Reaction on MHD Boundary Layer Flow o Williamson Nanoluid
More informationChapter 3 Water Flow in Pipes
The Islamic University o Gaza Faculty o Engineering Civil Engineering Department Hydraulics - ECI 33 Chapter 3 Water Flow in Pipes 3. Description o A Pipe Flow Water pipes in our homes and the distribution
More informationOne-Dimensional Motion Review IMPORTANT QUANTITIES Name Symbol Units Basic Equation Name Symbol Units Basic Equation Time t Seconds Velocity v m/s
One-Dimensional Motion Review IMPORTANT QUANTITIES Name Symbol Units Basic Equation Name Symbol Units Basic Equation Time t Seconds Velocity v m/s v x t Position x Meters Speed v m/s v t Length l Meters
More informationFLOW DYNAMICS OF COMPLEX FLUIDS USING NUMERICAL MODELS
DEPARTMENT OF PHYSICS UNIVERSITY OF JYVÄSKYLÄ RESEARCH REPORT No. 8/3 FLOW DYNAMICS OF COMPLEX FLUIDS USING NUMERICAL MODELS BY Amir Shakib-Manesh Academic Dissertation or the Degree o Doctor o Philosophy
More informationAn Improved Expression for a Classical Type of Explicit Approximation of the Colebrook White Equation with Only One Internal Iteration
International Journal o Hydraulic Engineering 06, 5(): 9-3 DOI: 0.593/j.ijhe.06050.03 An Improved Expression or a Classical Type o Explicit Approximation o the Colebrook White Equation with Only One Internal
More informationFLOW CHARACTERISTICS OF HFC-134a IN AN ADIABATIC HELICAL CAPILLARY TUBE
E HEFAT7 5 th International Conerence on Heat Transer, Fluid Mechanics and Thermodynamics Sun City, South Arica Paper number: KM1 FLOW CHARACTERISTICS OF HFC-1a IN AN ADIABATIC HELICAL CAPILLARY TUBE Khan
More informationThe Ascent Trajectory Optimization of Two-Stage-To-Orbit Aerospace Plane Based on Pseudospectral Method
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (014) 000 000 www.elsevier.com/locate/procedia APISAT014, 014 Asia-Paciic International Symposium on Aerospace Technology,
More informationIterative Methods for Stokes/Darcy Coupling
Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationChapter 3: Newtonian Fluid Mechanics QUICK START W
Chapter 3: Newtonian Fluid Mechanics TWO GOAL Derive governing equations (mass and momentum balances olve governing equations or velocity and stress ields QUICK TART W First, beore we get deep into derivation,
More informationarxiv: v1 [math.na] 29 Oct 2018
A Multirate Approach or Fluid-Structure Interaction Computation with Decoupled Methods Lian Zhang 1, Mingchao Cai 2 Mo Mu 3, arxiv:1810.11933v1 [math.na] 29 Oct 2018 Abstract We investigate a multirate
More informationD.R. Rector, M.L. Stewart and A.P. Poloski Pacific Northwest National Laboratory P.O. Box 999, Richland, WA
Modeling of Sediment Bed Behavior for Critical Velocity in Horizontal Piping 9263 D.R. Rector, M.L. Stewart and A.P. Poloski Pacific Northwest National Laboratory P.O. Box 999, Richland, WA ABSTRACT A
More informationFILM STACKING IMPREGNATION MODEL FOR THERMOPLASTIC COMPOSITES APPLIED TO A NOVEL NET-SHAPE PREFORMING PROCESS
The 8 th International Conerence on Flow Processes in Composite Materials (FPCM8) FILM STACKING IMPREGNATION MODEL FOR THERMOPLASTIC COMPOSITES APPLIED TO A NOVEL NET-SHAPE PREFORMING PROCESS S.T.Jespersen
More informationComputers & Fluids. Improving lattice Boltzmann simulation of moving particles in a viscous flow using local grid refinement
Computers & Fluids 00 (204) 2 Computers & Fluids Improving lattice Boltzmann simulation o moving particles in a viscous low using local grid reinement Songying Chen a,, Cheng Peng b, Yihua Teng c, Lian-Ping
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More information