Prediction of bubble growth rates and departure volumes in nucleate boiling at isolated sites

A model has been developed to predict heat transfer rates and sizes of bubbles generated during nucleate pool boiling. This model assumes conduction and a natural convective heat transfer mechanism through the liquid layer under the bubble and transient conduction from the bulk liquid. The temperature of the bulk liquid in the vicinity of the bubble is obtained by assuming a turbulent natural convection process from the hot plate to the liquid bulk. The shape of the bubble is obtained by equilibrium analysis. The bubble departure condition is predicted by a force balance equation. Good agreement has been found between the bubble radii predicted by the present theory and the ones obtained experimentally.

Bubble growth by injection of gas into viscous liquids in cylindrical and conical tubes.

The laminar low Mach number flow of a gas in a tube is analyzed for very small and very large values of the inlet-to-wall temperature ratio. When this ratio tends to zero, pressure forces confine the cold gas to a thin core around the axis of the tube. This core is neatly bounded by an ablation front that consumes it at a finite distance from the tube inlet. When the temperature ratio tends to infinity, the temperature of the gas increases smoothly from the wall to the axis of the tube and the shear stress and heat flux are positive at the wall despite the fact that the viscosity and thermal conductivity of the gas scaled with their inlet values tend to zero at the wall.

Bubble Growth From First Principles

We consider the simplest relevant problem in the foaming of molten plastics, the growth of a single bubble in a sea of highly viscous Newtonian fluid, and without interference from other bubbles. This simplest problem has defied accurate solution from first principles. Despite plenty of research on foaming, classical approaches from first principles have neglected the temperature rise in the surrounding fluid, and we find that this oversimplification greatly accelerates bubble growth prediction. We use a transport phenomena approach to analyze the growth of a solitary bubble, expanding under its own pressure. We consider a bubble of ideal gas growing without the accelerating contribution from mass transfer into the bubble. We explore the roles of viscous forces, fluid inertia, and viscous dissipation. We find that bubble growth depends upon the nucleus radius and nucleus pressure. We begin with a detailed examination of the classical approaches (thermodynamics without viscous heating). Our failure to fit experimental data with these classical approaches, sets up the second part of our paper, a novel exploration of the essential decelerating role of viscous heating. We explore both isothermal and adiabatic bubble expansion, and also the decelerating role of surface tension. The adiabatic analysis accounts for the slight deceleration due to the cooling of the expanding gas, which depends on gas polyatomicity. We also explore the pressure profile, and the components of the extra stress tensor, in the fluid surrounding the growing bubble. These stresses can eventually be frozen into foamed plastics. We find that our new theory compares well with measured bubble behavior.

Bubble Growth From First Principals

We consider the simplest relevant problem in the foaming of molten plastics, the growth of a single bubble in a sea of highly viscous Newtonian fluid, and without interference from other bubbles. This simplest problem has defied accurate solution from first principles. Despite plenty of research on foaming, classical approaches from first principles have neglected the temperature rise in the surrounding fluid, and we find that this oversimplification greatly accelerates bubble growth prediction. We use a transport phenomena approach to analyze the growth of a solitary bubble, expanding under its own pressure. We consider a bubble of ideal gas growing without the accelerating contribution from mass transfer into the bubble. We explore the roles of viscous forces, fluid inertia, and viscous dissipation. We find that bubble growth depends upon the nucleus radius and nucleus pressure. We begin with a detailed examination of the classical approaches (thermodynamics without viscous heating). Our failure to fit experimental data with these classical approaches, sets up the second part of our paper, a novel exploration of the essential decelerating role of viscous heating. We explore both isothermal and adiabatic bubble expansion, and also the decelerating role of surface tension. The adiabatic analysis accounts for the slight deceleration due to the cooling of the expanding gas, which depends on gas polyatomicity. We also explore the pressure profile, and the components of the extra stress tensor, in the fluid surrounding the growing bubble. These stresses can eventually be frozen into foamed plastics. We find that our new theory compares well with measured bubble behavior.

Effect of pulse heating parameters on the microscale bubble dynamics at a microheater surface

We study the effects of pulse heating parameters on the micro bubble behavior of a platinum microheater (100 mu m x 20 mu m) immersed in a methanol pool. The experiment covers the heat fluxes of 10-37 MW/m(2) and pulse frequencies of 25-500 Hz. The boiling incipience is initiated at the superheat limit of methanol, corresponding to the homogeneous nucleation. Three types of micro boiling patterns are identified. The first type is named as the bubble explosion and regrowth, consisting of a violent explosive boiling and shrinking, followed by a slower bubble regrowth and subsequent shrinking, occurring at lower heat fluxes. The second type, named as the bubble breakup and attraction, consists of the violent explosive boiling, bubble breakup and emission, bubble attraction and coalescence process, occurring at higher heat fluxes than those of the first type. The third type, named as the bubble size oscillation and large bubble formation, involves the initial explosive boiling, followed by a short periodic bubble growth and shrinking. Then the bubble continues to increase its size, until a constant bubble size is reached which is larger than the microheater length. 

Liquid-vapour phase change : nucleate boiling of pure fluid and nanofluid under different gravity levels

This research was a step towards the comprehension of the nano-particles interaction with bubbles created during boiling. It was aimed at solving the controversies of whether the heat transfer is enhanced or deteriorated during the boiling of the nanofluid. Experiments were conducted in normal gravity and reduced gravity environments on-board the European Space Agency Parabolic Flight Program. The local modification of the thermo-physical properties of the fluid and moreover the modification experienced in the liquid microlayer under the growing vapour bubble were the dominant factors in explaining the mechanisms of the boiling behaviour of the nanofluid.

Nucleation dynamics and pool boiling characteristics of high pressure refrigerant using thermochromic liquid crystals

This paper reports an experimental investigation of oscillating temperature field beneath a single isolated nucleation site using a non-invasive TLC (thermochromic liquid crystal) based thermography technique. Empirical correlations are presented to demonstrate the influence of system pressure and wall heat flux on different ebullition characteristics in the nucleate pool boiling regime of refrigerant R-134a. TLC transient response and two-phase flow structure are captured using synchronized, high resolution imaging. It is observed that the area of influence of nucleation site exhibits a two-part distinct transient behavior during the bubble growth period and broadens to a maximum of 1.57 times the bubble diameter at the instant of bubble departure. This is accompanied by a sharp fall of 2.5 degrees C in the local excess temperature at the nucleation site, which results in momentary augmentation (similar to 40%) in the local heat transfer coefficient at the nucleation origin. The enhanced heat transfer rate observed during the bubble peel-off event is primarily due to transient micro-convection in the wake of the retreating bubble. Further, the results indicate that a slight increase in system pressure from 813.6 to 882.5 kPa has no considerable effect on either the wall superheat or the overall heat transfer coefficient and ebullition frequency. In addition, correlations have been obtained for bubble Reynolds number, Jackob number and the dimensionless bubble generation frequency in terms of modified boiling number.

The non-linear acousto-mechanical and thermodynamic response of spherical absorbers to laser pulses of various temporal profiles

This study is to theoretically investigate shockwave and microbubble formation due to laser absorption by microparticles and nanoparticles. The initial motivation for this research was to understand the underlying physical mechanisms responsible for laser damage to the retina, as well as the predict threshold levels for damage for laser pulses with of progressively shorter durations. The strongest absorbers in the retina are micron size melanosomes, and their absorption of laser light causes them to accrue very high energy density. I theoretically investigate how this absorbed energy is transferred to the surrounding medium. For a wide range of conditions I calculate shockwave generation and bubble growth as a function of the three parameters; fluence, pulse duration and pulse shape. In order to develop a rigorous physical treatment, the governing equations for the behavior of an absorber and for the surrounding medium are derived. Shockwave theory is investigated and the conclusion is that a shock pressure explanation is likely to be the underlying physical cause of retinal damage at threshold fluences for sub-nanosecond pulses. The same effects are also expected for non-biological micro and nano absorbers. ^

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Effect of boiling on deposition of metallic colloids

In this study the relationship between heterogeneous nucleate boiling surfaces and deposition of suspended metallic colloidal particles, popularly known as crud or corrosion products in process industries, on those heterogeneous sites is investigated. Various researchers have reported that hematite is a major constituent of crud which makes it the primary material of interest; however the models developed in this work are irrespective of material choice. Qualitative hypotheses on the deposition process under boiling as proposed by previous researchers have been tested, which fail to provide explanations for several physical mechanisms observed and analyzed. In this study a quantitative model of deposition rate has been developed on the basis of bubble dynamics and colloid-surface interaction potential. Boiling from a heating surface aids in aggregation of the metallic particulates viz. nano-particles, crud particulate, etc. suspended in a liquid, which helps in transporting them to heating surfaces. Consequently, clusters of particles deposit onto the heating surfaces due to various interactive forces, resulting in formation of porous or impervious layers. The deposit layer grows or recedes depending upon variations in interparticle and surface forces, fluid shear, fluid chemistry, etc. This deposit layer in turn affects the rate of bubble generation, formation of porous chimneys, critical heat flux (CHF) of surfaces, activation and deactivation of nucleation sites on the heating surfaces. Several problems are posed due to the effect of boiling on colloidal deposition, which range from research initiatives involving nano-fluids as a heat transfer medium to industrial applications such as light water nuclear reactors. In this study, it is attempted to integrate colloid and surface science with vapor bubble dynamics, boiling heat transfer and evaporation rate. Pool boiling experiments with dilute metallic colloids have been conducted to investigate several parameters impacting the system. The experimental data available in the literature is obtained by flow experiments, which do not help in correlating boiling mechanism with the deposition amount or structure. With the help of experimental evidences and analysis, previously proposed hypothesis for particle transport to the contact line due to hydrophobicity has been challenged. The experimental observations suggest that deposition occurs around the bubble surface contact line and extends underneath area of the bubble microlayer as well. During the evaporation the concentration gradient of a non-volatile species is created, which induces osmotic pressure. The osmotic pressure developed inside the microlayer draws more particles inside the microlayer region or towards contact line. The colloidal escape time is slower than the evaporation time, which leads to the aggregation of particles in the evaporating micro-layer. These aggregated particles deposit onto or are removed from the heating surface, depending upon their total interaction potential. Interaction potential has been computed with the help of surface charge and van der Waals potential for the materials in aqueous solutions. Based upon the interaction-force boundary layer thickness, which is governed by debye radius (or ionic concentration and pH), a simplified quantitative model for the attachment kinetics is proposed. This attachment kinetics model gives reasonable results in predicting attachment rate against data reported by previous researchers. The attachment kinetics study has been done for different pH levels and particle sizes for hematite particles. Quantification of colloidal transport under boiling scenarios is done with the help of overall average evaporation rates because generally waiting times for bubbles at the same position is much larger than growth times. In other words, from a larger measurable scale perspective, frequency of bubbles dictates the rate of collection of particles rather than evaporation rate during micro-layer evaporation of one bubble. The combination of attachment kinetics and colloidal transport kinetics has been used to make a consolidated model for prediction of the amount of deposition and is validated with the help of high fidelity experimental data. In an attempt to understand and explain boiling characteristics, high speed visualization of bubble dynamics from a single artificial large cavity and multiple naturally occurring cavities is conducted. A bubble growth and departure dynamics model is developed for artificial active sites and is validated with the experimental data. The variation of bubble departure diameter with wall temperature is analyzed with experimental results and shows coherence with earlier studies. However, deposit traces after boiling experiments show that bubble contact diameter is essential to predict bubble departure dynamics, which has been ignored previously by various researchers. The relationship between porosity of colloid deposits and bubbles under the influence of Jakob number, sub-cooling and particle size has been developed. This also can be further utilized in variational wettability of the surface. Designing porous surfaces can having vast range of applications varying from high wettability, such as high critical heat flux boilers, to low wettability, such as efficient condensers.

Three-dimensional bubble clusters : Shape, packing, and growth rate

Cox, S.J., and Graner, F. (2004) Three-dimensional bubble clusters: shape, packing and growth-rate. Physical review. E, Statistical, nonlinear, and soft matter physics . 69:031409.

Economic growth in Sweden, 2000-2010 : The dot-com bubble and the financial crisis

    Economic growth is the increase in the inflation-adjusted market value of the goods and services produced by an economy over time. The total output is the quantity of goods or servicesproduced in a given time period within a country. Sweden was affected by two crises during the period 2000-2010: a dot-com bubble and a financial crisis. How did these two crises affect the economic growth?     The changes of domestic output can be separated into four parts: changes in intermediate demand, final domestic demand, export demand and import substitution. The main purpose of this article is to analyze the economic growth during the period 2000-2010, with focus on the dot-com bubble in the beginning of the period 2000-2005, and the financial crisis at the end of the period 2005-2010. The methodology to be used is the structural decomposition method.     This investigation shows that the main contributions to the Swedish total domestic output increase in both the period 2000-2005 and the period 2005-2010 were the effect of domestic demand. In the period 2005-2010, financial crisis weakened the effect of export. The output of the primary sector went from a negative change into a positive, explained mainly by strong export expansion. In the secondary sector, export had most effect in the period 2000-2005. Nevertheless, domestic demand and import ratio had more effect during the financial crisis period. Lastly, in the tertiary sector, domestic demand can mainly explain the output growth in the whole period 2000-2010.

The Brazilian real state market in 2012: robust growth or speculative bubble?

A forte alta dos imóveis no Brasil nos últimos anos iniciou um debate sobre a possível existência de uma bolha especulativa. Dada a recente crise do crédito nos Estados Unidos, é factível questionar se a situação atual no Brasil pode ser comparada à crise americana. Considerando argumentos quantitativos e fundamentais, examina-se o contexto imobiliário brasileiro e questiona-se a sustentabilidade em um futuro próximo. Primeiramente, analisou-se a taxa de aluguel e o nível de acesso aos imóveis e também utilizou-se um modelo do custo real para ver se o mercado está em equilíbrio o não. Depois examinou-se alguns fatores fundamentais que afetam o preço dos imóveis – oferta e demanda, crédito e regulação, fatores culturais – para encontrar evidências que justificam o aumento dos preços dos imóveis. A partir dessas observações tentou-se chegar a uma conclusão sobre a evolução dos preços no mercado imobiliário brasileiro. Enquanto os dados sugerem que os preços dos imóveis estão supervalorizados em comparação ao preço dos aluguéis, há evidências de uma legítima demanda por novos imóveis na emergente classe média brasileira. Um risco maior pode estar no mercado de crédito, altamente alavancado em relação ao consumidor brasileiro. No entanto, não se encontrou evidências que sugerem mais do que uma temporária estabilização ou correção no preço dos imóveis.

The organic fraction of bubble-generated, accumulation mode Sea Spray Aerosol (SSA)

Recent studies have detected a dominant accumulation mode (~100 nm) in the Sea Spray Aerosol (SSA) number distribution. There is evidence to suggest that particles in this mode are composed primarily of organics. To investigate this hypothesis we conducted experiments on NaCl, artificial SSA and natural SSA particles with a Volatility-Hygroscopicity-Tandem-Differential-Mobility-Analyser (VH-TDMA). NaCl particles were atomiser generated and a bubble generator was constructed to produce artificial and natural SSA particles. Natural seawater samples for use in the bubble generator were collected from biologically active, terrestrially-affected coastal water in Moreton Bay, Australia. Differences in the VH-TDMA-measured volatility curves of artificial and natural SSA particles were used to investigate and quantify the organic fraction of natural SSA particles. Hygroscopic Growth Factor (HGF) data, also obtained by the VH-TDMA, were used to confirm the conclusions drawn from the volatility data. Both datasets indicated that the organic fraction of our natural SSA particles evaporated in the VH-TDMA over the temperature range 170–200°C. The organic volume fraction for 71–77 nm natural SSA particles was 8±6%. Organic volume fraction did not vary significantly with varying water residence time (40 secs to 24 hrs) in the bubble generator or SSA particle diameter in the range 38–173 nm. At room temperature we measured shape- and Kelvin-corrected HGF at 90% RH of 2.46±0.02 for NaCl, 2.35±0.02 for artifical SSA and 2.26±0.02 for natural SSA particles. Overall, these results suggest that the natural accumulation mode SSA particles produced in these experiments contained only a minor organic fraction, which had little effect on hygroscopic growth. Our measurement of 8±6% is an order of magnitude below two previous measurements of the organic fraction in SSA particles of comparable sizes. We stress that our results were obtained using coastal seawater and they can’t necessarily be applied on a regional or global ocean scale. Nevertheless, considering the order of magnitude discrepancy between this and previous studies, further research with independent measurement techniques and a variety of different seawaters is required to better quantify how much organic material is present in accumulation mode SSA.