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: An experimental study of saturated pool boiling on a single artificial nucleation site without and with the application of an electric field on the boiling surface has been conducted. N-pentane is boiling on a copper surface and is recorded with a high speed camera providing high quality pictures and movies. The accuracy of the visualization allowed establishing an experimental bubble growth law from a large number of experiments. This law shows that the evaporation rate is decreasing during the bubble growth, and underlines the importance of liquid motion induced by the preceding bubble. Bubble rise is therefore studied: once detached, bubbles accelerate vertically until reaching a maximum velocity in good agreement with a correlation from literature. The bubbles then turn to another direction. The effect of applying an electric field on the boiling surface in finally studied. In addition to changes of the bubble shape, changes are also shown in the liquid plume and the convective structures above the surface. Lower maximum rising velocities were measured in the presence of electric fields, especially with a negative polarity.
Abstract: Boiling process is characterized by the rapid
formation of vapour bubbles at the solid–liquid interface (nucleate
boiling) with pre-existing vapour or gas pockets. Computational fluid
dynamics (CFD) is an important tool to study bubble dynamics. In
the present study, CFD simulation has been carried out to determine
the bubble detachment diameter and its terminal velocity. Volume of
fluid method is used to model the bubble and the surrounding by
solving single set of momentum equations and tracking the volume
fraction of each of the fluids throughout the domain. In the
simulation, bubble is generated by allowing water-vapour to enter a
cylinder filled with liquid water through an inlet at the bottom. After
the bubble is fully formed, the bubble detaches from the surface and
rises up during which the bubble accelerates due to the net balance
between buoyancy force and viscous drag. Finally when these forces
exactly balance each other, it attains a constant terminal velocity. The
bubble detachment diameter and the terminal velocity of the bubble
are captured by the monitor function provided in FLUENT. The
detachment diameter and the terminal velocity obtained are compared
with the established results based on the shape of the bubble. A good
agreement is obtained between the results obtained from simulation
and the equations in comparison with the established results.
Abstract: A bubbly flow in a vertical miniature tube is analyzed theoretically. The liquid and gas phase are co-current flowing upward. The gas phase is injected via a nozzle whose inner diameter is 0.11mm and it is placed on the axis of the tube. A force balance is applied on the bubble at its detachment. The set of governing equations are solved by use of Mathematica software. The bubble diameter and the bubble generation frequency are determined for various inlet phase velocities represented by the inlet mass quality. The results show different behavior of bubble growth and detachment depending on the tube size.
Abstract: In this research, the changes in bubbles diameter and
number that may occur due to the change in heat flux of pure water
during pool boiling process. For this purpose, test equipment was
designed and developed to collect test data. The bubbles were graded
using Caliper Screen software. To calculate the growth and
nucleation rates of bubbles under different fluxes, population balance
model was employed. The results show that the increase in heat flux
from q=20 kw/m2 to q= 102 kw/m2 raised the growth and nucleation
rates of bubbles.
Abstract: With the long-term objective of Critical Heat Flux (CHF) prediction, a Direct Numerical Simulation (DNS) framework for simulation of subcooled and saturated nucleate pool boiling is developed. One case of saturated, and three cases of subcooled boiling at different subcooling levels are simulated. Grid refinement study is also reported. Both boiling and condensation phenomena can be computed simultaneously in the proposed numerical framework. Computed bubble detachment diameters of the saturated nucleate pool boiling cases agree well with the experiment. The flow structures around the growing bubble are presented and the accompanying physics is described. The relation between heat flux evolution from the heated wall and the bubble growth is studied, along with investigations of temperature distribution and flow field evolutions.