GRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE

Size: px
Start display at page:

Download "GRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE"

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

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 information

DNS OF FALLING DROPLETS IN A VERTICAL CHANNEL

DNS 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 information

A STUDY ON BINARY COLLISION OF GNF DROPLETS USING A CONSERVATIVE LEVEL-SET METHOD

A 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 information

Portal del coneixement obert de la UPC

Portal 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:

, 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 information

Oblique Drop Impact on Deep and Shallow Liquid

Oblique 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 information

A 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 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 information

Height function interface reconstruction algorithm for the simulation of boiling flows

Height 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 information

Level 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 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 information

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 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 information

INTERACTION BETWEEN A PAIR OF DROPS ASCENDING IN A LINEARLY STRATIFIED FLUID

INTERACTION 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 information

Direct Numerical Simulation of Single Bubble Rising in Viscous Stagnant Liquid

Direct 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 information

LARGE-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 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 information

Numerical simulation of wave breaking in turbulent two-phase Couette flow

Numerical 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 information

Fluid-Structure Interaction of a Reed Type Valve

Fluid-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 information

Fluid Dynamics Exercises and questions for the course

Fluid 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 information

Drop 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 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 information

Flow Field and Oscillation Frequency of a Rotating Liquid Droplet

Flow 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 information

An OpenFOAM-based electro-hydrodynamical model

An 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 information

We 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. 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 information

Application 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 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 information

MULTI-SCALE VOLUME OF FLUID MODELLING OF DROPLET COALESCENCE

MULTI-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 information

ON THE FLOW PAST A CIRCULAR CYLINDER FROM

ON 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 information

Numerical 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 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 information

Computational Simulation of Marangoni Convection Under Microgravity Condition

Computational 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 information

Fundamentals of Fluid Dynamics: Elementary Viscous Flow

Fundamentals 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 information

IHMTC EULER-EULER TWO-FLUID MODEL BASED CODE DEVELOPMENT FOR TWO-PHASE FLOW SYSTEMS

IHMTC 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 information

Nonlinear oscillations and rotations of a liquid droplet

Nonlinear 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 information

Direct Numerical Simulations of Gas-Liquid Flows

Direct 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 information

Simulation 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, 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 information

Fluid Dynamics: Theory, Computation, and Numerical Simulation Second Edition

Fluid 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 information

A numerical study on the effects of cavitation on orifice flow

A 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 information

CCC Annual Report. UIUC, August 19, Argon Bubble Behavior in EMBr Field. Kai Jin. Department of Mechanical Science & Engineering

CCC 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 information

FEM-Level Set Techniques for Multiphase Flow --- Some recent results

FEM-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 information

Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles

Model 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 information

2. FLUID-FLOW EQUATIONS SPRING 2019

2. 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

Simulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework

Simulating 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 information

Experience 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 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 information

Reduction of parasitic currents in the DNS VOF code FS3D

Reduction 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 information

A numerical analysis of drop impact on liquid film by using a level set method

A 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 information

Modeling Complex Flows! Direct Numerical Simulations! Computational Fluid Dynamics!

Modeling 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 information

NON-OBERBECK-BOUSSINESQ NATURAL CONVECTION IN A TALL DIFFERENTIALLY HEATED CAVITY

NON-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 information

Contribution 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 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 information

A 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 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 information

On surface tension modelling using the level set method

On 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 information

Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition

Detailed 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 information

IMPINGEMENT OF A DROPLET ONTO A DRY WALL: A NUMERICAL INVESTIGATION

IMPINGEMENT 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 information

6.2 Governing Equations for Natural Convection

6.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 information

Simulation of floating bodies with lattice Boltzmann

Simulation 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 information

Implementation of a symmetry-preserving discretization in Gerris

Implementation 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 information

Experiments at the University of Minnesota (draft 2)

Experiments 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 information

Phase-Field simulation of small capillary-number two-phase flow in a microtube

Phase-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 information

A MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD

A 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 information

MODELING OF DIESEL- COMPRESSED NATURAL GAS BUBBLY FLOW UNDER INFLUENCING OF A MAGNETIC FIELD

MODELING 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 information

Topics in Other Lectures Droplet Groups and Array Instability of Injected Liquid Liquid Fuel-Films

Topics 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 information

Numerical study of bubble dynamics with the Boundary Element Method

Numerical 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 information

Application of the Immersed Boundary Method to particle-laden and bubbly flows

Application 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 information

NEAR-WALL TURBULENCE-BUBBLES INTERACTIONS IN A CHANNEL FLOW AT Re =400: A DNS/LES INVESTIGATION

NEAR-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 information

Numerical Simulation of the Hagemann Entrainment Experiments

Numerical 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 information

CFD in COMSOL Multiphysics

CFD 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 information

FUNDAMENTAL STUDY OF BINGHAM FLUID BY MEANS OF DAM-BREAK FLOW MODEL

FUNDAMENTAL 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 information

ABSTRACT NUMERICAL SIMULATION OF LIQUID - VAPOR INTERFACES IN THE SHARP INTERFACE LIMIT

ABSTRACT 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 information

FOUR-WAY COUPLED SIMULATIONS OF TURBULENT

FOUR-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 information

3D Compressible Simulation Of a Muffler With Pseudosound Prediction Levels

3D 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 information

The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets

The 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 information

The Lateral Migration of a Drop under Gravity between Two Parallel Plates at Finite Reynolds Numbers

The 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 information

Fluid Mechanics Theory I

Fluid 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 information

Numerical study of wall effects on buoyant gasbubble rise in a liquid-filled finite cylinder

Numerical 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 information

ScienceDirect. Evaporating falling drop

ScienceDirect. 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 information

Modeling of Humidification in Comsol Multiphysics 4.4

Modeling 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 information

PHYS 432 Physics of Fluids: Instabilities

PHYS 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 information

Numerical 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 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 information

Development and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations

Development 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 information

Numerical Study of Air Inside Refrigerating Compartment of Frost-free Domestic Refrigerators

Numerical 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 information

CFD modelling of multiphase flows

CFD 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 information

Contents. I Introduction 1. Preface. xiii

Contents. 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 information

A 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 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 information

Absorption of gas by a falling liquid film

Absorption 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 information

FLUID MECHANICS. Chapter 9 Flow over Immersed Bodies

FLUID 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 information

NUMERICAL SIMULATION OF GAS-LIQUID TWO-PHASE CONVECTIVE HEAT TRANSFER IN A MICRO TUBE

NUMERICAL 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 information

DIRECT 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 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 information

Numerical 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 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 information

Numerical Analysis of Metal Transfer in Gas Metal Arc Welding

Numerical 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 information

Tutorial for the supercritical pressure pipe with STAR-CCM+

Tutorial 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 information

International Journal of Multiphase Flow

International 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 information

DIRECT NUMERICAL SIMULATION OF LIQUID- SOLID FLOW

DIRECT 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 information

A formulation for fast computations of rigid particulate flows

A 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 information

1D-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 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 information

NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS

NONLINEAR 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 information

NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE

NUMERICAL 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 information

LINEAR 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 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 information

Validation 3. Laminar Flow Around a Circular Cylinder

Validation 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 information

Two-Dimensional and Axisymmetric Viscous. Flow in Apertures

Two-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 information

Study on residence time distribution of CSTR using CFD

Study 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 information

A 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 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 information

Modified DLM method for finite-volume simulation of particle flow

Modified 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 information

x j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2

x 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 information

The Role of Splatting Effect in High Schmidt Number Turbulent Mass Transfer Across an Air-Water Interface

The 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 information

Modelling multiphase flows in the Chemical and Process Industry

Modelling 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 information

NUMERICAL SIMULATIONS OF BUBBLE DYNAMICS. Dr. Nikos A. Pelekasis

NUMERICAL 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