Abstract: Multi-modal film boiling simulations are carried out on adaptive octree grids. The liquid-vapor interface is captured using the volume-of-fluid framework adjusted to account for exchanges of mass, momentum, and energy across the interface. Surface tension effects are included using a volumetric source term in the momentum equations. The phase change calculations are conducted based on the exact location and orientation of the interface; however, the source terms are calculated using the mixture variables to be consistent with the one field formulation used to represent the entire fluid domain. The numerical model on octree representation of the computational grid is first verified using test cases including advection tests in severely deforming velocity fields, gravity-based instabilities and bubble growth in uniformly superheated liquid under zero gravity. The model is then used to simulate both single and multi-modal film boiling simulations. The octree grid is dynamically adapted in order to maintain the highest grid resolution on the instability fronts using markers of interface location, volume fraction, and thermal gradients. The method thus provides an efficient platform to simulate fluid instabilities with or without phase change in the presence of body forces like gravity or shear layer instabilities.
Abstract: This paper investigates the dynamic behavior of laminar water flows inside superhydrophobic micro-pipes patterned with square micro-posts features under different operating conditions. It also investigates the effects of air fraction and Reynolds number on the frictional performance of these pipes. Rather than modeling the air-water interfaces of superhydrophobic as a flat inflexible surface, a transient, incompressible, three-dimensional, volume-of-fluid (VOF) methodology has been employed to continuously track the air–water interface shape inside micro-pipes. Also, the entrance effects on the flow field have been taken into consideration. The results revealed the strong dependency of the frictional performance on the air fractions and Reynolds number. The frictional resistance reduction becomes increasingly more significant at large air fractions and low Reynolds numbers. Increasing Reynolds number has an adverse effect on the frictional resistance reduction.
Abstract: Carbon Capture & Storage (CCS) is one of the various
methods that can be used to reduce the carbon footprint of the
energy sector. This paper focuses on the absorption of CO2 from
flue gas using packed columns, whose efficiency is highly dependent
on the structure of the liquid films within the column. To study the
characteristics of liquid films a CFD solver, OpenFOAM is utilised
to solve two-phase, isothermal film flow using the volume-of-fluid
(VOF) method. The model was validated using existing experimental
data and the Nusselt theory. It was found that smaller plate inclination
angles, with respect to the horizontal plane, resulted in larger wetted
areas on smooth plates. However, only a slight improvement in
the wetted area was observed. Simulations were also performed
using a ridged plate and it was observed that these surface textures
significantly increase the wetted area of the plate. This was mainly
attributed to the channelling effect of the ridges, which helped to
oppose the surface tension forces trying to minimise the surface area.
Rivulet formations on the ridged plate were also flattened out and
spread across a larger proportion of the plate width.
Abstract: In this article, various models of surface tension force (CSF, CSS and PCIL) for interfacial flows have been applied to dynamic case and the results were compared. We studied the Kelvin- Helmholtz instabilities, which are produced by shear at the interface between two fluids with different physical properties. The velocity inlet is defined as a sinusoidal perturbation. When gravity and surface tension are taking into account, we observe the development of the Instability for a critic value of the difference of velocity of the both fluids. The VOF Model enables to simulate Kelvin-Helmholtz Instability as dynamic case.
Abstract: The effect of a uniform magnetic field on the
formation of drops of specific size has been investigated numerically
in a T-shaped microchannel. Previous researches indicated that the
drop sizes of secondary stream decreases, with increasing main
stream flow rate and decreasing interfacial tension. In the present
study the effect of a uniform magnetic field on the main stream is
considered, and it is proposed that by increasing the Hartmann
number, the size of the drops of the secondary stream will be
decreased.
Abstract: The Navier–Stokes equations for unsteady, incompressible, viscous fluids in the axisymmetric coordinate system are solved using a control volume method. The volume-of-fluid (VOF) technique is used to track the free-surface of the liquid. Model predictions are in good agreement with experimental measurements. It is found that the dynamic processes after impact are sensitive to the initial droplet velocity and the liquid pool depth. The time evolution of the crown height and diameter are obtained by numerical simulation. The critical We number for splashing (Wecr) is studied for Oh (Ohnesorge) numbers in the range of 0.01~0.1; the results compares well with those of the experiments.
Abstract: Different numerical methods are employed and developed for simulating interfacial flows. A large range of applications belong to this group, e.g. two-phase flows of air bubbles in water or water drops in air. In such problems surface tension effects often play a dominant role. In this paper, various models of surface tension force for interfacial flows, the CSF, CSS, PCIL and SGIP models have been applied to simulate the motion of small air bubbles in water and the results were compared and reviewed. It has been pointed out that by using SGIP or PCIL models, we are able to simulate bubble rise and obtain results in close agreement with the experimental data.
Abstract: In this paper, the two-dimensional stagger grid
interface pressure (SGIP) model has been generalized and presented
into three-dimensional form. For this purpose, various models of
surface tension force for interfacial flows have been investigated and
compared with each other. The VOF method has been used for
tracking the interface. To show the ability of the SGIP model for
three-dimensional flows in comparison with other models, pressure
contours, maximum spurious velocities, norm spurious flow
velocities and pressure jump error for motionless drop of liquid and
bubble of gas are calculated using different models. It has been
pointed out that SGIP model in comparison with the CSF, CSS and
PCIL models produces the least maximum and norm spurious
velocities. Additionally, the new model produces more accurate
results in calculating the pressure jumps across the interface for
motionless drop of liquid and bubble of gas which is generated in
surface tension force.