Abstract:In this study, the direct numerical simulation of a fullydeveloped turbulent channel flow with deformed bubbles were conducted by means of the MARS method, turbulent Reynolds number 150, and Particle Reynolds number 120. As the results, large scale wake motions round bubbles were confirmed. At the bubble related region, mean velocity was degreased and turbulent intensites and Reynolds shear stress were incresed by the effects of the large scale wake motions round bubbles. On the other hands, near wall region, effects of bubbly flow might act on flow laminarlize and drag reduction. Two types of drag coefficient of bubble were estimated from the accelerated velocity of bubble and experimental equation as a function of Particle Reynolds number. Drag coefficient estimated by the accelerated velocity of bubble was below the results estimated by the experimental equation. INTRODUCTION Turbulent flows with bubbles are very often found in engineering devices such as a light-water nuclear reactor and chemical reactor. Therefore, both experimental and numerical investigations have been conducted, extensively. In the view points of understanding coherent structures and turbulent statistics behaviours, the DNS (Direct Numerical Simulation) is expected to have advantages of the experimental approaches. Recently, DNS of the turbulent shear flows with deformability bubbles have been conducted by Kawamura and Kodama, 2002, Lu et. al., 2005 and Lu and Tryggvason, 2008, based on the fronttracking method. These DNS focused on interaction between turbulent structures near wall and bubbles to investigate friction drag reduction effects on wall. In this study, DNS of a fully-developed turbulent channel flow with bubbles were conducted to investigate forces acted on deformability bubbles such as drag coefficient. NUMERICAL METHOD Numerical procedure was based on the MARS method (Kunugi, 1997, Kunugi, 2001); the governing equations are consisted of Navier-Stokes equations with the CSF (Continuum Surface Force) model (Brackbill, 1992) , continuity equation and transport equations of a volume fraction function (VOF) F. To refine the suface tention acted on bubbles, a kind of a level-set function estimated by using VOF function, was used for calculation of the surface tention. As for the discretization of the governing equations on the Cartesian coordinate system, the second-order scheme for the spatial differencing terms is used on the staggered grid system. Time integration method was lowstorage 3rd-order Runge-Kutta scheme for convection terms, Crank-Nicolson scheme for viscous terms and Euler Implicit scheme for Pressure terms, respectively. The physical problem treated here is the motion of two Newtonian incompressible fluids allowed the interface deformation between them. NUMERICAL CONDITION Figure 1 shows flow geometry and coordinate systems and Numerical condition was tabled in Table 1. In this study, minimum flow units (Jimenez and Moin, 1991) was selected as numerical domain. Thermal properties at normal airwater system were adapted. A fully developed single-phase turbulent channel flow data was used as the initial velocity fields and 4 bubbles which diameter was 1 mm were entrained in channel center domain. Sixth International Symposium on Turbulence and Shear Flow Phenomena Seoul, Korea, 22-24 June 2009

Numerical investigation on the drag force of a single bubble and bubble swarm

Experiment and Numerical Simulation of the Motion of Bubble Swarm in a Bubble Column

A numerical study on the drag law of a gas bubble using dynamic body force method

Numerical Modelling of Bubbling Flow Dynamic Behaviour in a Bubble Column

Numerical Investigation of Microbubble Drag Reduction on an Axisymmetric Body Based on Eulerian Multiphase Model

Direct numerical simulation of deformable rising bubbles at low Reynolds numbers

Experimental And Numerical Study Of Bubble Drag Reduction On A Flat Plate

Drag Force of Bubble Swarms and Numerical Simulations of a Bubble Column with a CFD-PBM Coupled Model

Numerical Study on the Deformation of a Gas Bubble in Uniform Flow

Analysis of Drag and Lift Coefficient Expressions of Bubbly Flow System for Low to Medium Reynolds Number

Stream-Wise Distribution of Skin-Friction Drag Reduction on A Flat Plate with Bubble Injection

Numerical Studies on Effects of Bubbles Regular Array on the Liquid-Phase Turbulence

Spatiotemporal evolutions of forces and vortices of flow past ellipsoidal bubbles: Direct numerical simulation based on a Cartesian grid scheme

DRAG FLUCTUATIONS ON DEFORMED BUBBLES IN A TURBULENT SHEAR FLOW BY MEANS OF DNS

Analyses and Modified Models for Bubble Shape and Drag Coefficient Covering a Wide Range of Working Conditions

Euler–Lagrange Study of Bubble Breakup and Coalescence in a Turbulent Boundary Layer for Bubble Drag Reduction

Correlation of Drag Coefficient Between Rising Bubbles in Chain

The Characteristics and Mechanism of Microbubble Drag Reduction on the Axisymmetric Body

Numerical Analysis of Single Pulsating Bubble Interaction with a Free Surface in a Viscous Compressible Liquid

Euler-Lagrange Study of Bubble Drag Reduction in Turbulent Channel Flow and Boundary Layer Flow

An Improved Computational Fluid Dynamics-Discrete Bubble Model for Multiphase Bubble Flow Simulation with a Wide Grid/bubble Size Ratio