GRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE
|
|
- Neal McKinney
- 5 years ago
- Views:
Transcription
1 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 Rigola, Assensi Oliva Heat and Mass Transfer Technological Center (CTTC) Universidad Politécnica de Cataluña ETSEIAT, Colom 11, Terrassa (Barcelona), Spain cttc@upc.edu, web page: Termo Fluids, S.L. Avda Jacquard 97 1-E, Terrassa (Barcelona), Spain termofluids@termofluids.com Key words: Two-phase flow, Level-set method, Bubbles, Droplets, Parallel CFD applications Abstract. In this work, Direct Numerical Simulations of gravity-driven motion of a swarm of bubbles in a vertical pipe are presented. These simulations were carried out in the contex of a multiple marker level-set methodology, which is integrated in a finitevolume framework on collocated unstructured grids. Each fluid particle is described by a separate level-set function, thus, different interfaces can be solved at the same control volume, avoiding the coalescence of the bubbles. Present simulations have been performed in a periodic vertical domain divided in hexahedral control volumes, distributed in 2048 CPUs. Collective and individual behaviour of the bubbles are characterized and compared with previous results reported in the literature for other flow configurations. 1 INTRODUCTION Gravity-driven bubbly flows play an important role in many natural and industrial processes. Steam generators in nuclear plants, unit operations in chemical engineering such as distillation, absorption, extraction, heterogeneous catalysis and bubble reactors are only a few among a multitude of applications that involve the motion of bubble swarms [9]. These applications have motivated a large number of numerical and experimental investigations of bubble dynamics, however, despite those efforts, the current understanding of such flows and their predictive models are far from satisfactory, and many challenging problems still remain as pointed out in recent reviews [9, 13]. The development of parallel computation has promoted Direct Numerical Simulation (DNS) of the Navier-Stokes equations as another means of performing controlled experiments [14, 13], providing a good way to non-invasive measure of droplet and bubble flows, although computationally demanding. There are multiple methods for DNS of bubbly flows, for instance: the front tracking (FT) method [14], level set (LS) methods 1
2 [11, 10, 1, 2, 3], volume-of-fluid (VOF) methods [7], and hybrid VOF/LS methods [12]. In these methods, two-phase flow is treated as a single flow with the density and viscosity varying smoothly across the moving interface which is captured in an Eulerian framework (VOF, LS, VOF/LS) or in a Lagrangian framework (FT), moreover, although the idea behind these methods is similar, their numerical implementation may differ greatly. Forthepresent DNSstudy, amultiplemarker level-set methoddeployed in[2]isusedto perform simulations of gravity-driven motion of a swarm of bubbles in a periodic vertical pipe. The main advantage of this method is its capacity to avoid the numerical coalescence of the bubbles, which is an artifact of interface capturing methods such as VOF and LS approaches. Hence, bubbles are able to approach each other closely within the size of one grid cell, and can even collide, while the bubble volume is kept constant throughout the simulation [2]. To the best of the authors knowledge, there are not previous DNS studies of gravity-driven bubbly flows in vertical pipes, therefore this is the main motivation of the present research. This paper is organized in the following order: The governing equations and numerical methods are presented in Section 2. Section 3 is devoted to the discussion of the numerical results. Finally, concluding remarks are given in Section 4. 2 MATHEMATICAL FORMULATION AND NUMERICAL METHODS The Navier-Stokes equations for the dispersed fluid in Ω d and continuous fluid in Ω c can be combined into a set of equations in a global domain Ω = Ω d Ω c, with a singular source term for the surface tension force at the interface Γ: (ρv) t + (ρvv) = p+ µ ( v+( v) T) +(ρ ρ 0 )g+f σ v = 0 (1) ρ = ρ d H d +ρ c (1 H d ) µ = µ d H d +µ c (1 H d ) (2) wherevandpdenotethefluidvelocityandpressurefieldrespectively, ρisthefluiddensity, µ is the dynamic viscosity, g is the gravitational acceleration, f σ is the surface tension force, subscripts d and c are used for the dispersed and continuous fluids respectively, while H d the Heaviside step function that is one in Ω d and zero elsewhere. Because a periodic domain is used in the y axis direction, a force ρ 0 g is added to the Navier-Stokes equations [3], with ρ 0 = V 1 Ω (ρ Ω dφ d +ρ c (1 φ d ))dv. The conservative level-set method (CLS) deployed in [1] has been selected for interface capturing on unstructured meshes. Moreover, multiple level-set functions are used in order to avoid the coalescence of the bubbles, according to the work reported by [2], therefore, the interface of the i th fluid particle is defined as the 0.5 iso-surface of a regularized indicator function φ i, where i = 1,2,...,n d and n d is the total number of fluid particles contained by the dispersed phase. The i th interface transport equation can be written in conservative form provided the velocity field is solenoidal, v = 0, namely, φ i t + φ iv = 0 (3) 2
3 Furthermore, an additional re-initialization equation is introduced in order to keep a sharp and constant interface profile φ i τ + φ i(1 φ i )n i = ε φ i (4) This equation is advanced in pseudo-time τ up to steady state. It consists of a compressive term, φ i (1 φ i )n i τ=0, which forces the level-set function to be compressed onto the interface along the normal vector n i, and of a diffusion term ε φ i that ensure the profile remains of characteristic thickness ε. Geometrical information on the interface Γ i, such as normal vector n i or curvature κ i, is obtained through: n i (φ i ) = φ i φ i κ i (φ i ) = n i (5) Surface tension forces are calculated by the continuous surface force model [5], which has been extended to include multiple markers in the same grid cell [2], f σ = σ i κ i (φ i ) φ i (6) Finally, in order to avoid numerical instabilities at the interface, fluid properties in Eq. 2 are regularized by using a global level-set function H d = φ d, defined as φ d (x,t) = max{φ 1 (x,t),...,φ nd 1(x,t),φ nd (x,t)} (7) The governing equations have been discretized on a collocated unstructured grid arrangement by means of the finite-volume method, according to [1]. In order to avoid unphysical oscillations of the level-set function, a TVD Superbee limiter [1] is used to discretize the convective term of advection Eq. 3, while the compressive term of the reinitialization Eq. 4 is discretized by using a central difference scheme (CD). CD scheme is also used to discretize the convective term of momentum Eq. 1. A distance-weighted linear interpolation is used to find the cell-face values, while gradients are computed at cell centroids by using the least-squares method [1]. The velocity-pressure coupling has been solved by means of a classical fractional step projection method. The reader is referred to [1, 2, 3] for technical details on the discretization of the Navier-Stokes and Level-set equations on unstructured meshes that are beyond the scope of this paper. 3 NUMERICAL RESULTS AND DISCUSSION The multiple marker level-set method has been validated and verified in previous publications by performing simulations of single bubbles rising in an initially quiescent liquid [1, 3], drop collision with a fluid-fluid interface [2] and binary droplet collision with bouncing outcome [2]. For present simulations, the computational set-up is defined as a vertical cylinder bounded by a rigid wall, with gravity in the y direction. The size of the domain is (D Ω,H Ω ) = (5d,4d), where d is the initial bubble diameter, D Ω is the cylinder diameter 3
4 and H Ω is the cylinder height, as shown in Fig. 3a. The domain Ω is divided in hexahedral control volumes distributed on 2048 CPUs, whereas the grid size is given by h = d/60. As initial condition, both bubbles and liquid are quiescent. Imposed boundary conditions are non-slip at the rigid wall and periodic on the streamwise (y-direction). In this way bubbles go out of the domain on the top side, and they come back in the domain again from the opposite side. Bubbles and droplets rising or falling freely in infinite media can be characterized by four dimensionless numbers: M gµ4 c ρ ρ 2 c σ3 Eo gd2 ρ c σ η ρ ρ c ρ d η µ µ c µ d (8) where, η ρ andη µ arerespectively thedensity andviscosity ratio; M isthemortonnumber; Eo is the Eötvös number; ρ = ρ c ρ d, specifies the density difference between the continuous and dispersed fluid phases, and d is the spherical volume equivalent diameter of the droplet, while the subscript d denotes the dispersed fluid phase and c the continuous fluid phase. Simulations presented in this work were performed at Eo = 3, M = 10 6, η ρ = 10 and η µ = 10, which correspond to deformable bubbles. Finally, in order to report the numerical results, the velocity of the bubble centroid, v i, and the Reynolds number, Re i, are calculated as follows: v i (t) = vφ Ω i(x,t)dv φ Ω i(x,t)dv Re i (t) = ρ cd v i e y µ c for i = 1,..,n d (9) where e y is a unit vector parallel to +y direction. In addition the average Reynolds number of the bubble swarm is calculated as Re d = 1 n d Re i (10) n d Once the velocity of each bubble has been calculated by Eq. 9, the trajectory of the bubble can be determined by: x i (t) = x 0 i + i=1 t 0 v i (t)dt (11) where x 0 i is the initial position of the i th bubble centroid. First, we consider the motion of two bubbles in an initially quiescent fluid where the bubbles are initially released with different configuration angles and initial centroidcentroid distance of 2d in a periodic vertical pipe. As the bubbles move upward due to the buoyancy effect, their interaction leads to the formation of complex trajectories, as is illustrated in Fig. 1. Figure 2a shows the dimensionless separation distance between the bubbles, s/d, versus the dimensionless time for the different initial configuration angles, θ 0 = {0 o,45 o,60 o,75 o,90 o }. It can be seen that two deformable bubbles repel each other, except for θ 0 = 90 o where the bubbles experience the drafting-kissing-tumbling (DKT) 4
5 Figure 1: Interaction of two deformable bubbles in a vertical pipe at Eo = 3, M = 10 6, η ρ = 10 and η µ = 10, for different configuration angles Figure 2: Interaction of two deformable bubbles at Eo = 3, M = 10 6, η ρ = 10 and η µ = 10 (a) The centroid-centroid distance between two bubbles versus dimensionless time. (b) The angle between two bubbles versus dimensionless time. phenomenon. Moreover, the bubbles tend to migrate to the wall, however a repulsion force avoids the bubble-wall collision. Figure 2b illustrates the time evolution of θ for different initial configuration angles θ 0. A decrease of θ is observed as the time advances, 5
6 Figure 3: Eo = 3, M = 10 6, η ρ = 10 and η µ = 10. (a) Bubble distribution and velocity vectors on the symmetry plane x y for t = {2,20} (b) e z v on the plane x y, t = 20. (c) Bubble trajectories integrated from t = 0 up to t = 28. therefore it indicates a torque action on the bubbles that tends to align them side-by-side. We now study a free bubble array of 16 fluid particles which are initially placed in the periodic vertical pipe following a random pattern. This corresponds to a dilute bubbly flow, with an overall volumetric fraction of α Ω = 5.33%, moreover, since both fluids are assumed to be incompressible and coalescence of the bubbles is not allowed, α is constant throughout the simulation, so that statistics can be obtained for a constant number of bubbles. Fig. 3a shows snapshots of the swarm of bubbles and their periodic images at the dimensionless time t = t(g/d) 1/2 = {2,20}, and the velocity vectors on the symmetry plane of the cylinder. Fig. 3b illustrates the voticity field e y v generated by both the wall and bubbles on the x y plane at z = 0 and t = 20. Close to the wall, the vorticity generated by bubbles and wall have opposite signs, moreover, its interaction produces a velocity reduction, which leads to a pressure increment due to the Bernoulli effect, and a repulsion force directed from the wall to the bubbles. Fig. 3c shows the bubble trajectories integrated from t = 0 up to t = 28, where it is observed that bubbles move not only in streamwise direction, both also in radial direction. Fig. 4a shows the Reynolds number versus time, defined for each bubble as Re i = ρ c d(v i e y )µ 1 c with i = 1,2,..,n d. Even though the individual bubble motion show a transient behaviour due the wake interaction, the whole swarm of bubbles reaches a steady state at approximately t = 4. Fig. 4b shows the radial position for each bubble centroid, which lead to the formation of a bubble curtain between the wall and the symmetry axis, produced by the alignment of the bubbles at approximately constant distance from the wall. When 6
7 Figure 4: (a) Time evolution of Re for each bubble, and averaged Re. (b) Time evolution of radial position for each bubble centroid. all bubbles are at approximately the same vertical position, the probability of collision increases and these interactions are addressed by the DKT behaviour, which lead to the formation of horizontal clusters, as shown Fig. 4a. Present results are consistent with simulations reported by [6] using the front-tracking method [13]. 4 CONCLUSIONS Direct Numerical Simulations of gravity-driven motion of a swarm of bubbles in a vertical channel have been performed. These numerical experiments demonstrate the ability of the multiple marker level-set method [2] to simulate bubbly flows without numerical coalescence. Regarding the interaction of two-bubbles in a vertical pipe, a repulsion force was observed except for the in-line configuration where the bubbles follows the DKT phenomenon. For the present flow conditions, it was observed that deformable bubbles do not collide with the wall, moreover bubble-bubble interactions in the swarm follows also the DKT behaviour which leads to the formation of horizontal clusters. 5 ACKNOWLEDGMENTS This work has been financially supported by the Ministerio de Economía y Competitividad, Secretaría de Estado de Investigación, Desarrollo e Innovación, Spain (ENE , ENE R). The simulations reported in this work were carried out using computer time provided by PRACE (project ) through the MareNostrum III supercomputer based in Barcelona, Spain. REFERENCES [1] Balcázar, N., Jofre, L., Lehmkhul, O., Castro, J., Rigola, J., A finitevolume/level-set method for simulating two-phase flows on unstructured grids. In- 7
8 ternational Journal of Multiphase Flow 64 : [2] Balcázar, N., Lehmkhul, O., Rigola, J., Oliva, A., A multiple marker level-set method for simulation of deformable fluid particles. International Journal of Multiphase Flow 74 : [3] Balcázar, N., Lehmkhul, O., Rigola, J., Oliva, A Level-set simulations of buoyancy-driven motion of single and multiple bubbles. International Journal of Heat and Fluid Flow 56 : [4] Balcázar, N., Jofre, L., Lehmkuhl, O.,Castro, J., Oliva, A., A multiple marker levelset method for simulation of bubbly flows, Proceedings of the 6th European Conference on Computational Fluid Dynamics (ECFD VI) (2014) [5] Brackbill, J.U., Kothe, D.B., Zemach, C., A Continuum Method for Modeling Surface Tension, J. Comput. Phys. (1992) 100 : [6] Bunner, B., Tryggvason, G., Effect of bubble deformation on the properties of bubbly flows. J. Fluid Mech. (2003) 495 : [7] Hirt, C., Nichols, B., Volume of fluid (VOF) method for the dynamics of free boundary. J. Comput. Phys. (1981) 39 : [8] Lehmkuhl, O., Perez-Segarra, C.D., Soria, M., Oliva, A., A new Parallel unstructured CFD code for the simulation of turbulent industrial problems on low cost PC cluster, Proceedings of the Parallel CFD 2007 Conference (2007), pp.1 8. [9] Mudde, R., Gravity-Driven bubbly flows, Annu. Rev. Fluid Mech. (2005), pp [10] Olsson, E., Kreiss, G., A conservative level set method for two phase flow. J. Comput. Phys. (2005) 210 : [11] Sussman, M., Smereka, P., Osher, S., A Level Set Approach for Computing Solutions to Incompressible Two-Phase Flow, J. Comput. Phys. (1994), 144 : [12] Sussman, M., Smereka, P., Osher, S., A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows, J. Comput. Phys. (2000), 162 : [13] Tryggvason, G., Dabiri, S., Abouhasanzadeh, B., Jaicai, L., Multiscale considerations in direct numerical simulations of multiphase flows. Phys. Fluids (2013) [14] Tryggvason, G., Bunner, B., Esmaeeli, A., Juric, D., Al-Rawahi, N., Tauber, W., Han, J., Nas, S., Jan, Y-J., A Front-Tracking Method for the Computations of Multiphase Flow. J. Comput. Phys. (2001) 169 :
Portal del coneixement obert de la UPC
UPCommons Portal del coneixement obert de la UPC http://upcommons.upc.edu/e-prints 2017. Aquesta versió està disponible sota la llicència CC-BY-NC-ND 4.0 http://creativecommons.org/licenses/by-nc-nd/4.0/
More informationDNS OF FALLING DROPLETS IN A VERTICAL CHANNEL
N. Balcázar, et al., Int. J. Comp. Meth. and Exp. Meas., Vol. 6, No. 2 (2018) 398 410 DNS OF FALLING DROPLETS IN A VERTICAL CHANNEL NÉSTOR BALCÁZAR 1,2, JESÚS CASTRO 1, JORGE CHIVA 1 & ASSENSI OLIVA 1
More informationA STUDY ON BINARY COLLISION OF GNF DROPLETS USING A CONSERVATIVE LEVEL-SET METHOD
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK A STUDY ON BINARY COLLISION OF GNF DROPLETS USING
More informationPortal del coneixement obert de la UPC
UPCommons Portal del coneixement obert de la UPC http://upcommons.upc.edu/e-prints Publicat sota llicència de la Ltd. a Journal of Physics: Conference Series El contingut d'aquest treball pot ser utilitzat
More information, where the -function is equal to:
Paper ID ILASS08-000 ILASS08-9-4 ILASS 2008 Sep. 8-10, 2008, Como Lake, Italy BINARY COLLISION BETWEEN UNEQUAL SIZED DROPLETS. A NUMERICAL INVESTIGATION. N. Nikolopoulos 1, A. Theodorakakos 2 and G. Bergeles
More informationOblique Drop Impact on Deep and Shallow Liquid
1 2 3 4 5 6 7 8 9 11 10 12 13 14 15 16 17 18 19 20 21 22 23 24 25 Commun. Comput. Phys. doi: 10.4208/cicp.XXX.XXX Oblique Drop Impact on Deep and Shallow Liquid B. Ray 1, G. Biswas 1,2, and A. Sharma 3
More informationA Two-Phase Flow Interface Tracking Algorithm Using a Fully Coupled Pressure-Based Finite Volume Method
A Two-Phase Flow Interface Tracking Algorithm Using a Fully Coupled Pressure-Based Finite Volume Method Shidvash Vakilipour, Scott Ormiston, Masoud Mohammadi, Rouzbeh Riazi, Kimia Amiri, Sahar Barati Abstract
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 informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationDetailed 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 informationINTERACTION BETWEEN A PAIR OF DROPS ASCENDING IN A LINEARLY STRATIFIED FLUID
Proceedings of the ASME 2013 Fluids Engineering Division Summer Meeting FEDSM 2013 July 7-11, 2013, Incline Village, Nevada, USA FEDSM 2013-16046 INTERACTION BETWEEN A PAIR OF DROPS ASCENDING IN A LINEARLY
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 informationLARGE-EDDY SIMULATIONS OF FLUID FLOW AND HEAT TRANSFER AROUND A PARABOLIC TROUGH SOLAR COLLECTOR
LARGE-EDDY SIMULATIONS OF FLUID FLOW AND HEAT TRANSFER AROUND A PARABOLIC TROUGH SOLAR COLLECTOR Ahmed Amine Hachicha, Ivette Rodríguez and Assensi Oliva 1 1 Heat and Mass Transfer Technological Center
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationFluid-Structure Interaction of a Reed Type Valve
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2016 Fluid-Structure Interaction of a Reed Type Valve Ignacio Gonzalez Heat and Mass Transfer
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
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 informationFlow Field and Oscillation Frequency of a Rotating Liquid Droplet
Flow Field and Oscillation Frequency of a Rotating Liquid Droplet TADASHI WATANABE Center for Computational Science and e-systems Japan Atomic Energy Agency (JAEA) Tokai-mura, Naka-gun, Ibaraki-ken, 319-1195
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 informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,9 116, 12M Open access books available International authors and editors Downloads Our authors
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationMULTI-SCALE VOLUME OF FLUID MODELLING OF DROPLET COALESCENCE
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 MULTI-SCALE VOLUME OF FLUID MODELLING OF DROPLET COALESCENCE Lachlan R. MASON,
More informationON THE FLOW PAST A CIRCULAR CYLINDER FROM
ON THE FLOW PAST A CIRCULAR CYLINDER FROM CRITICAL TO SUPER-CRITICAL REYNOLDS NUMBERS: WAKE TOPOLOGY AND VORTEX SHEDDING I. Rodríguez 1, O. Lehmkuhl 1,2, J. Chiva 1, R. Borrell 2 and A. Oliva 1 Heat and
More informationNumerical simulation of heat transfer and fluid flow in a flat plate solar collector with TIM and ventilation channel
Numerical simulation of heat transfer and fluid flow in a flat plate solar collector with TIM and ventilation channel Hamdi Kessentini 1, Roser Capdevila 1, Oriol Lehmkuhl 1,2, Jesus Castro 1 and Assensi
More informationComputational Simulation of Marangoni Convection Under Microgravity Condition
Transaction B: Mechanical Engineering Vol. 16, No. 6, pp. 513{524 c Sharif University of Technology, December 2009 Computational Simulation of Marangoni Convection Under Microgravity Condition Abstract.
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationIHMTC EULER-EULER TWO-FLUID MODEL BASED CODE DEVELOPMENT FOR TWO-PHASE FLOW SYSTEMS
Proceedings of the 24th National and 2nd International ISHMT-ASTFE Heat and Mass Transfer Conference (IHMTC-2017), December 27-30, 2017, BITS-Pilani, Hyderabad, India IHMTC2017-13-0160 EULER-EULER TWO-FLUID
More informationNonlinear oscillations and rotations of a liquid droplet
Nonlinear oscillations and rotations of a liquid droplet Tadashi Watanabe watanabe.tadashi66@jaea.go.jp Abstract Nonlinear oscillations and rotations of a liquid droplet are simulated numerically by solving
More informationDirect Numerical Simulations of Gas-Liquid Flows
Direct Numerical Simulations of Gas-Liquid Flows 1 Gretar Tryggvason*; 1 Jiacai Lu; 2 Ming Ma 1 Johns Hopkins University, Baltimore, MD, USA; 2 University of Notre Dame, Notre Dame, IN, USA Introduction
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationFluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition
Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition C. Pozrikidis m Springer Contents Preface v 1 Introduction to Kinematics 1 1.1 Fluids and solids 1 1.2 Fluid parcels and flow
More informationA numerical study on the effects of cavitation on orifice flow
PHSICS OF FLUIDS, A numerical study on the effects of cavitation on orifice flow S. Dabiri, W. A. Sirignano, and D. D. Joseph, University of California, Irvine, California 9697, USA University of Minnesota,
More informationCCC Annual Report. UIUC, August 19, Argon Bubble Behavior in EMBr Field. Kai Jin. Department of Mechanical Science & Engineering
CCC Annual Report UIUC, August 19, 2015 Argon Bubble Behavior in EMBr Field Kai Jin Department of Mechanical Science & Engineering University of Illinois at Urbana-Champaign Introduction Argon bubbles
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
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 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 informationExperience with DNS of particulate flow using a variant of the immersed boundary method
Experience with DNS of particulate flow using a variant of the immersed boundary method Markus Uhlmann Numerical Simulation and Modeling Unit CIEMAT Madrid, Spain ECCOMAS CFD 2006 Motivation wide range
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 informationA numerical analysis of drop impact on liquid film by using a level set method
Journal of Mechanical Science and Technology 25 (10) (2011) 2567~2572 wwwspringerlinkcom/content/1738-494x DOI 101007/s12206-011-0613-7 A numerical analysis of drop impact on liquid film by using a level
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 informationNON-OBERBECK-BOUSSINESQ NATURAL CONVECTION IN A TALL DIFFERENTIALLY HEATED CAVITY
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 21 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal,14-17 June 21 NON-OBERBECK-BOUSSINESQ NATURAL CONVECTION IN A TALL DIFFERENTIALLY
More informationContribution of inter-particle collisions on kinetic energy modification in a turbulent channel flow
Contribution of inter-particle collisions on kinetic energy modification in a turbulent channel flow Valentina Lavezzo a, Alfredo Soldati a,b a Dipartimento di Energetica e Macchine and b Centro Interdipartimentale
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 informationOn surface tension modelling using the level set method
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2008) Published online in Wiley InterScience (www.interscience.wiley.com)..1804 On surface tension modelling using the
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationIMPINGEMENT OF A DROPLET ONTO A DRY WALL: A NUMERICAL INVESTIGATION
IMPINGEMENT OF A DROPLET ONTO A DRY WALL:. Introduction A NUMERICAL INVESTIGATION N. Nikolopoulos, and G. Bergeles Department Mechanical Engineering Nat. Technical University of Athens 57 Zografos, Greece
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationSimulation of floating bodies with lattice Boltzmann
Simulation of floating bodies with lattice Boltzmann by Simon Bogner, 17.11.2011, Lehrstuhl für Systemsimulation, Friedrich-Alexander Universität Erlangen 1 Simulation of floating bodies with lattice Boltzmann
More informationImplementation of a symmetry-preserving discretization in Gerris
Implementation of a symmetry-preserving discretization in Gerris Daniel Fuster Cols: Pierre Sagaut, Stephane Popinet Université Pierre et Marie Curie, Institut Jean Le Rond D Alembert Introduction 10/11:
More informationExperiments at the University of Minnesota (draft 2)
Experiments at the University of Minnesota (draft 2) September 17, 2001 Studies of migration and lift and of the orientation of particles in shear flows Experiments to determine positions of spherical
More informationPhase-Field simulation of small capillary-number two-phase flow in a microtube
Fluid Dynamics Research 40 (008) 497 509 Phase-Field simulation of small capillary-number two-phase flow in a microtube Qunwu He, Nobuhide Kasagi Department of Mechanical Engineering, The University of
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 informationMODELING OF DIESEL- COMPRESSED NATURAL GAS BUBBLY FLOW UNDER INFLUENCING OF A MAGNETIC FIELD
Journal of Engineering Science and Technology Vol. 12, No. 7 (2017) 1930-1938 School of Engineering, Taylor s University MODELING OF DIESEL- COMPRESSED NATURAL GAS BUBBLY FLOW UNDER INFLUENCING OF A MAGNETIC
More informationTopics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid Fuel-Films
Lecture Topics Transient Droplet Vaporization Convective Vaporization Liquid Circulation Transcritical Thermodynamics Droplet Drag and Motion Spray Computations Turbulence Effects Topics in Other Lectures
More informationNumerical study of bubble dynamics with the Boundary Element Method
Numerical study of bubble dynamics with the Boundary Element Method N. Méndez and R. González-Cinca Departament de Física Aplicada, Barcelona Tech-Universitat Politècnica de Catalunya Esteve Terradas 5,
More informationApplication of the Immersed Boundary Method to particle-laden and bubbly flows
Application of the Immersed Boundary Method to particle-laden and bubbly flows, B. Vowinckel, S. Schwarz, C. Santarelli, J. Fröhlich Institute of Fluid Mechanics TU Dresden, Germany EUROMECH Colloquium
More informationNEAR-WALL TURBULENCE-BUBBLES INTERACTIONS IN A CHANNEL FLOW AT Re =400: A DNS/LES INVESTIGATION
ABSTRACT NEAR-WALL TURBULENCE-BUBBLES INTERACTIONS IN A CHANNEL FLOW AT Re =400: A DNS/LES INVESTIGATION D. Métrailler, S. Reboux and D. Lakehal ASCOMP GmbH Zurich, Technoparkstr. 1, Switzerland Metrailler@ascomp.ch;
More informationNumerical Simulation of the Hagemann Entrainment Experiments
CCC Annual Report UIUC, August 14, 2013 Numerical Simulation of the Hagemann Entrainment Experiments Kenneth Swartz (BSME Student) Lance C. Hibbeler (Ph.D. Student) Department of Mechanical Science & Engineering
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationFUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL
Annual Journal of Hydraulic Engineering, JSCE, Vol.54, 2010, February FUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL How Tion PUAY1 and Takashi HOSODA2 1 Member of JSCE, Phd Student,
More informationABSTRACT NUMERICAL SIMULATION OF LIQUID - VAPOR INTERFACES IN THE SHARP INTERFACE LIMIT
ABSTRACT Title of Thesis: NUMERICAL SIMULATION OF LIQUID - VAPOR INTERFACES IN THE SHARP INTERFACE LIMIT Vamsee K. Yerramilli, Master of Science, 2005 Thesis directed by: Dr. Elias Balaras Department of
More informationFOUR-WAY COUPLED SIMULATIONS OF TURBULENT
FOUR-WAY COUPLED SIMULATIONS OF TURBULENT FLOWS WITH NON-SPHERICAL PARTICLES Berend van Wachem Thermofluids Division, Department of Mechanical Engineering Imperial College London Exhibition Road, London,
More information3D Compressible Simulation Of a Muffler With Pseudosound Prediction Levels
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2016 3D Compressible Simulation Of a Muffler With Pseudosound Prediction Levels Jesus Ruano
More informationThe effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets
ILASS Americas, 26 th Annual Conference on Liquid Atomization and Spray Systems, Portland, OR, May 204 The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization
More informationThe Lateral Migration of a Drop under Gravity between Two Parallel Plates at Finite Reynolds Numbers
Journal of Applied Fluid Mechanics, Vol. 5, No. 1, pp. 11-21, 2012. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. The Lateral Migration of a Drop under Gravity between Two Parallel
More informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationNumerical study of wall effects on buoyant gasbubble rise in a liquid-filled finite cylinder
University of Pennsylvania ScholarlyCommons Departmental Papers (MEAM) Department of Mechanical Engineering & Applied Mechanics September 7 Numerical study of wall effects on buoyant gasbubble rise in
More informationScienceDirect. Evaporating falling drop
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 15 (2015 ) 201 206 IUTAM Symposium on Multiphase flows with phase change: challenges and opportunities, Hyderabad, India (December
More informationModeling of Humidification in Comsol Multiphysics 4.4
Modeling of Humidification in Comsol Multiphysics 4.4 Indrajit Wadgaonkar *1 and Suresh Arikapudi 1 1 Tata Motors Ltd. Pimpri, Pune, India, 411018. *Corresponding author: Indrajit Wadgaonkar, Tata Motors
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More informationNumerical study of liquid metal film flows in a varying spanwise. magnetic field
Numerical study of liquid metal film flows in a varying spanwise magnetic field D. Gao, N.B. Morley, V. Dhir Department of Mechanical and Aerospace Engineering, University of California, Los Angeles, CA
More informationDevelopment and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Development and analysis of a Lagrange-Remap sharp interface solver
More informationNumerical Study of Air Inside Refrigerating Compartment of Frost-free Domestic Refrigerators
Purdue University Purdue e-pubs International Refrigeration and Air Conditioning Conference School of Mechanical Engineering 2010 Numerical Study of Air Inside Refrigerating Compartment of Frost-free Domestic
More informationCFD modelling of multiphase flows
1 Lecture CFD-3 CFD modelling of multiphase flows Simon Lo CD-adapco Trident House, Basil Hill Road Didcot, OX11 7HJ, UK simon.lo@cd-adapco.com 2 VOF Free surface flows LMP Droplet flows Liquid film DEM
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
More informationA fundamental study of the flow past a circular cylinder using Abaqus/CFD
A fundamental study of the flow past a circular cylinder using Abaqus/CFD Masami Sato, and Takaya Kobayashi Mechanical Design & Analysis Corporation Abstract: The latest release of Abaqus version 6.10
More informationAbsorption of gas by a falling liquid film
Absorption of gas by a falling liquid film Christoph Albert Dieter Bothe Mathematical Modeling and Analysis Center of Smart Interfaces/ IRTG 1529 Darmstadt University of Technology 4th Japanese-German
More informationFLUID MECHANICS. Chapter 9 Flow over Immersed Bodies
FLUID MECHANICS Chapter 9 Flow over Immersed Bodies CHAP 9. FLOW OVER IMMERSED BODIES CONTENTS 9.1 General External Flow Characteristics 9.3 Drag 9.4 Lift 9.1 General External Flow Characteristics 9.1.1
More informationNUMERICAL SIMULATION OF GAS-LIQUID TWO-PHASE CONVECTIVE HEAT TRANSFER IN A MICRO TUBE
ECI International Conference on Heat Transfer and Fluid Flow in Microscale Castelvecchio Pascoli, 5-30 September 005 NUMERICAL SIMULATION OF GAS-LIQUID TWO-PHASE CONVECTIVE HEAT TRANSFER IN A MICRO TUBE
More informationDIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING
DIRECT NUMERICAL SIMULATION OF SPATIALLY DEVELOPING TURBULENT BOUNDARY LAYER FOR SKIN FRICTION DRAG REDUCTION BY WALL SURFACE-HEATING OR COOLING Yukinori Kametani Department of mechanical engineering Keio
More informationNumerical Simulation of Gas-Liquid-Reactors with Bubbly Flows using a Hybrid Multiphase-CFD Approach
Numerical Simulation of Gas-Liquid-Reactors with Bubbly Flows using a Hybrid Multiphase-CFD Approach TFM Hybrid Interface Resolving Two-Fluid Model (HIRES-TFM) by Coupling of the Volume-of-Fluid (VOF)
More informationNumerical Analysis of Metal Transfer in Gas Metal Arc Welding
Numerical Analysis of Metal Transfer in Gas Metal Arc Welding G. WANG, P.G. HUANG, and Y.M. ZHANG The present article describes a numerical procedure to simulate metal transfer and the model will be used
More informationTutorial for the supercritical pressure pipe with STAR-CCM+
Tutorial for the supercritical pressure pipe with STAR-CCM+ For performing this tutorial, it is necessary to have already studied the tutorial on the upward bend. In fact, after getting abilities with
More informationInternational Journal of Multiphase Flow
International Journal of Multiphase Flow 52 (2013) 22 34 Contents lists available at SciVerse ScienceDirect International Journal of Multiphase Flow journal homepage: www.elsevier.com/locate/ijmulflow
More informationDIRECT NUMERICAL SIMULATION OF LIQUID- SOLID FLOW
DIRECT NUMERICAL SIMULATION OF LIQUID- SOLID FLOW http://www.aem.umn.edu/solid-liquid_flows Sponsored by NSF-Grand Challenge Grant Fluid Mechanics & CFD Computer Scientists D.D. Joseph Y. Saad R. Glowinski
More informationA formulation for fast computations of rigid particulate flows
Center for Turbulence Research Annual Research Briefs 2001 185 A formulation for fast computations of rigid particulate flows By N. A. Patankar 1. Introduction A formulation is presented for the direct
More information1D-3D COUPLED SIMULATION OF THE FUEL INJECTION INSIDE A HIGH PERFORMANCE ENGINE FOR MOTORSPORT APPLICATION: SPRAY TARGETING AND INJECTION TIMING
1D-3D COUPLED SIMULATION OF THE FUEL INJECTION INSIDE A HIGH PERFORMANCE ENGINE FOR MOTORSPORT APPLICATION: SPRAY TARGETING AND INJECTION TIMING M. Fiocco, D. Borghesi- Mahindra Racing S.P.A. Outline Introduction
More informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
More informationNUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00786 NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A
More informationLINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION
LINEAR STABILITY ANALYSIS AND DIRECT NUMERICAL SIMUALATION OF DOUBLE-LAYER RAYLEIGH-BÉNARD CONVECTION É. FONTANA 1, E. MANCUSI 1,2, A. A. ULSON DE SOUZA 1, S. M. A. GUELLI U. SOUZA 1 1 Universidade Federal
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 informationTwo-Dimensional and Axisymmetric Viscous. Flow in Apertures
Under consideration for publication in J. Fluid Mech. Two-Dimensional and Axisymmetric Viscous Flow in Apertures By S A D E G H D A B I R I, W I L L I A M A. S I R I G N A N O A N D D A N I E L D. J O
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 informationA Coupled VOF-Eulerian Multiphase CFD Model To Simulate Breaking Wave Impacts On Offshore Structures
A Coupled VOF-Eulerian Multiphase CFD Model To Simulate Breaking Wave Impacts On Offshore Structures Pietro Danilo Tomaselli Ph.d. student Section for Fluid Mechanics, Coastal and Maritime Engineering
More informationModified DLM method for finite-volume simulation of particle flow
Modified DLM method for finite-volume simulation of particle flow A. M. Ardekani, S. Dabiri, and R. H. Rangel Department of Mechanical and Aerospace Engineering, University of California, Irvine, CA 92697-3975,
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 informationThe Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface
Turbulence, Heat and Mass Transfer 4 K. Hanjalic, Y. Nagano and M. Tummers (Editors) 3 Begell House, Inc. The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water
More informationModelling multiphase flows in the Chemical and Process Industry
Modelling multiphase flows in the Chemical and Process Industry Simon Lo 9/11/09 Contents Breakup and coalescence in bubbly flows Particle flows with the Discrete Element Modelling approach Multiphase
More informationNUMERICAL SIMULATIONS OF BUBBLE DYNAMICS. Dr. Nikos A. Pelekasis
NUMERICAL SIMULATIONS OF BUBBLE DYNAMICS Dr. Nikos A. Pelekasis Assistant Professor of Computational Fluid Dynamics, Dept. Mechanical & Industrial Engineering, University of Thessaly, Volos 38334, GREECE
More information