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

Direct Numerical Simulation of Complex Multi-Fluid Flows Using a Combined Volume of Fluid and Immersed Boundary Method

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of complex multi-fluid flows in which simultaneously (moving) deformable (drops or bubbles) and non-deformable (moving) elements (particles) are present, possibly with the additional presence of free surfaces. Our model combines the VOF model developed by van Sint Annaland et al. (2005) and the Immersed Boundary (IB) model developed by van der Hoef et al. (2006). The Volume of Fluid (VOF) part features i) an interface reconstruction technique based on piecewise linear interface representation ii) a three-dimensional version of the CSF model of Brackbill et al. (1992). The Immersed Boundary (IB) part incorporates both particle-fluid and particle-particle interaction via a Direct Forcing Method (DFM) and a hard sphere Discrete Particle (DP) approach. In our model a fixed (Eulerian) grid is utilized to solve the Navier-Stokes equations for the entire computational domain. The no-slip condition at the surface of the moving particles is enforced via a momentum source term which only acts in the vicinity of the particle surface. For the enforcement of the no-slip condition Lagrangian force points are used which are distributed evenly over the surface of the particle. Dissipative particle-particle and/or particle-wall collisions are accounted via a hard sphere DP approach (Hoomans et al., 1996) using a three-parameter particle-particle interaction model accounting for normal and tangential restitution and tangential friction. The capabilities of the hybrid VOF-IB model are demonstrated with a number of examples in which complex topological changes in the interface are encountered

Direct numerical simulation of the drag force in bubble swarms

This paper studies the swarm effect on the drag force in bubbly flows. From literature it is well-known that for relatively small bubbles, the drag force increases with the bubble hold-up due to additional hindrance experienced by the bubbles caused by the modified flow field. Very large (spherical cap) bubbles on the other hand may rise cooperatively. The unique capabilities of a 3D Front Tracking model were used to investigate the influence of important parameters like the gas fraction, Reynolds number and the bubble size independently. It was found that the relative drag force increases for bubbles in the range of 2 to 5 mm when the gas fraction is increased up to 13%, while the bubbles become more spherical. Also the influence of the Reynolds number and the bubble aspect ratio on the increased drag force has been determined. It can be concluded that there is only a very weak effect over several decades of the Reynolds number, while there is a strong effect of the bubble aspect ratio.\ud This also helps explaining why the increase in drag is smaller for larger bubbles: when the gas fraction is increased deformable bubbles become more spherical, thus reducing the drag force