Numerical Simulation of Heat Transport in Dispersed Gas-Liquid Two-Phase Flow using a Front Tracking Approach

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of heat transport in dispersed gas-liquid two-phase flow using the Front Tracking (FT) approach. Our model extends the FT model developed by van Sint Annaland et al. (2006) to non-isothermal conditions. In FT an unstructured dynamic mesh is used to represent and track the interface explicitly by a number of interconnected marker points. The Lagrangian representation of the interface avoids the necessity to reconstruct the interface from the local distribution of the fractions of the phases and, moreover, allows a direct and accurate calculation of the surface tension force circumventing the (problematic) computation of the interface curvature. The extended model is applied to predict the heat exchange rate between the liquid and a hot wall kept at a fixed temperature. It is found that the wall-to-liquid heat transfer coefficient exhibits a maximum in the vicinity of the bubble that can be attributed to the locally decreased thickness of the thermal boundary layer

Direct numerical simulation of heat transport in dispersed gas-liquid two-phase flow using a front tracking approach

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of heat transport in dispersed gas-liquid two-phase flow using the Front Tracking (FT) approach. Our model extends the FT model developed by van Sint Annaland et al. (2006) to non-isothermal conditions. In FT an unstructured dynamic mesh is used to represent and track the interface explicitly by a number of interconnected marker points. The Lagrangian representation of the interface avoids the necessity to reconstruct the interface from the local distribution of the fractions of the phases and, moreover, allows a direct and accurate calculation of the surface tension force circumventing the (problematic) computation of the interface curvature. The extended model is applied to predict the heat exchange rate between the liquid and a hot wall kept at a fixed temperature. It is found that the wall-to-liquid heat transfer coefficient exhibits a maximum in the vicinity of the bubble that can be attributed to the locally decreased thickness of the thermal boundary layer

Dynamic simulation of dispersed gas-liquid two-phase flow using a discrete bubble model.

In this paper a detailed hydrodynamic model for gas-liquid two-phase flow will be presented. The model is based on a mixed Eulerian-Lagrangian approach and describes the time-dependent two-dimensional motion of small, spherical gas bubbles in a bubble column operating in the homogeneous regime. The motion of these bubbles is calculated from a force balance for each individual bubble, accounting for all relevant forces acting on them. Contributions from liquid-phase pressure gradient, drag, virtual mass, liquid-phase vorticity and gravity are considered, whereas direct bubble-bubble interactions are accounted for via an interaction model resembling the collision model developed by Hoomans et al. (1996) to model gas-fluidized beds. The liquid-phase hydrodynamics are described using the volume-averaged, unsteady, Navier-Stokes equations. A preliminary model validation has been performed by comparing the computational results with experimental observations published previously in literature by various authors. The model is shown to predict correctly the motion of a bubble plume in a pseudo-two-dimensional bubble column operated at different superficial gas velocities, provided that a detailed description of the bubble dynamics is incorporated in the model. The effect of bubble column aspect ratio on the hydrodynamic behaviour of the column has also been investigated. Our model predicts the effect of aspect ratio on the flow structure in the bubble column. The importance of the various forces acting on the bubbles will also be discussed and it will be shown that the added mass force and the lift force cannot be neglected in bubble column simulation. Finally, the model has been used to study the start-up behaviour of a two-dimensional bubble column. It will be shown that the history of the gas-liquid two-phase flow significantly affects the flow structure ultimately obtained in a bubble column. This finding has, to our knowledge, not been reported before in literature

Preferential accumulation of bubbles in Couette-Taylor flow patterns

We investigate the migration of bubbles in several flow patterns occurring within the gap between a rotating inner cylinder and a concentric fixed outer cylinder. The time-dependent evolution of the two-phase flow is predicted through three-dimensional Euler-Lagrange simulations. Lagrangian tracking of spherical bubbles is coupled with direct numerical simulation of the Navier-Stokes equations. We assume that bubbles do not influence the background flow (one-way coupling simulations). The force balance on each bubble takes into account buoyancy, added-mass, viscous drag and shear-induced lift forces. For increasing velocities of the rotating inner cylinder, the flow in the fluid gap evolves from the purely azimuthal steady Couette flow to Taylor toroidal vortices and eventually a wavy vortex flow. The migration of bubbles is highly dependent on the balance between buoyancy and centripetal forces (mostly due to the centripetal pressure gradient) directed toward the inner cylinder and the vortex cores. Depending on the rotation rate of the inner cylinder, bubbles tend to accumulate alternatively along the inner wall, inside the core of Taylor vortices or at particular locations within the wavy vortices. A stability analysis of the fixed points associated with bubble trajectories provides a clear understanding of their migration and preferential accumulation. The location of the accumulation points is parameterized by two dimensionless parameters expressing the balance of buoyancy, centripetal attraction toward the inner rotating cylinder, and entrapment in Taylor vortices. A complete phase diagram summarizing the various regimes of bubble migration is built. Several experimental conditions considered by Djéridi et al.1 are reproduced; the numerical results reveal a very good agreement with the experiments. When the rotation rate is further increased, the numerical results indicate the formation of oscillating bubble strings, as observed experimentally by Djéridi et al.2. After a transient state, bubbles collect at the crests or troughs of the wavy vortices. An analysis of the flow characteristics clearly indicates that bubbles accumulate in the low-pressure regions of the flow field

Bubbly and Buoyant Particle-Laden Turbulent Flows

Fluid turbulence is commonly associated with stronger drag, greater heat transfer, and more efficient mixing than in laminar flows. In many natural and industrial settings, turbulent liquid flows contain suspensions of dispersed bubbles and light particles. Recently, much attention has been devoted to understanding the behavior and underlying physics of such flows by use of both experiments and high-resolution direct numerical simulations. This review summarizes our present understanding of various phenomenological aspects of bubbly and buoyant particle-laden turbulent flows. We begin by discussing different dynamical regimes, including those of crossing trajectories and wake-induced oscillations of rising particles, and regimes in which bubbles and particles preferentially accumulate near walls or within vortical structures. We then address how certain paradigmatic turbulent flows, such as homogeneous isotropic turbulence, channel flow, Taylor-Couette turbulence, and thermally driven turbulence, are modified by the presence of these dispersed bubbles and buoyant particles. We end with a list of summary points and future research questions.Comment: 29 pages, 14 figure

Eulerian-Lagrangian method for simulation of cloud cavitation

We present a coupled Eulerian-Lagrangian method to simulate cloud cavitation in a compressible liquid. The method is designed to capture the strong, volumetric oscillations of each bubble and the bubble-scattered acoustics. The dynamics of the bubbly mixture is formulated using volume-averaged equations of motion. The continuous phase is discretized on an Eulerian grid and integrated using a high-order, finite-volume weighted essentially non-oscillatory (WENO) scheme, while the gas phase is modeled as spherical, Lagrangian point-bubbles at the sub-grid scale, each of whose radial evolution is tracked by solving the Keller-Miksis equation. The volume of bubbles is mapped onto the Eulerian grid as the void fraction by using a regularization (smearing) kernel. In the most general case, where the bubble distribution is arbitrary, three-dimensional Cartesian grids are used for spatial discretization. In order to reduce the computational cost for problems possessing translational or rotational homogeneities, we spatially average the governing equations along the direction of symmetry and discretize the continuous phase on two-dimensional or axi-symmetric grids, respectively. We specify a regularization kernel that maps the three-dimensional distribution of bubbles onto the field of an averaged two-dimensional or axi-symmetric void fraction. A closure is developed to model the pressure fluctuations at the sub-grid scale as synthetic noise. For the examples considered here, modeling the sub-grid pressure fluctuations as white noise agrees a priori with computed distributions from three-dimensional simulations, and suffices, a posteriori, to accurately reproduce the statistics of the bubble dynamics. The numerical method and its verification are described by considering test cases of the dynamics of a single bubble and cloud cavitaiton induced by ultrasound fields.Comment: 28 pages, 16 figure

Deformable ellipsoidal bubbles in Taylor-Couette flow with enhanced Euler-Lagrange tracking

In this work we present numerical simulations of $10^5$ sub-Kolmogorov deformable bubbles dispersed in Taylor-Couette flow (a wall-bounded shear system) with rotating inner cylinder and outer cylinder at rest. We study the effect of deformability of the bubbles on the overall drag induced by the carrier fluid in the two-phase system. We find that an increase in deformability of the bubbles results in enhanced drag reduction due to a more pronounced accumulation of the deformed bubbles near the driving inner wall. This preferential accumulation is induced by an increase in the resistance on the motion of the bubbles in the wall-normal direction. The increased resistance is linked to the strong deformation of the bubbles near the wall which makes them prolate (stretched along one axes) and orient along the stream-wise direction. A larger concentration of the bubbles near the driving wall implies that they are more effective in weakening the plume ejections which results in stronger drag reduction effects. These simulations which are practically impossible with fully resolved techniques are made possible by coupling a sub-grid deformation model with two-way coupled Euler-Lagrangian tracking of sub-Kolmogorov bubbles dispersed in a turbulent flow field which is solved through direct numerical simulations. The bubbles are considered to be ellipsoidal in shape and their deformation is governed by an evolution equation which depends on the local flow conditions and their surface tension