A comparative study of two-phase flow models relevant to bubble column dynamics
Minev, P. D.; Lange, U.; Nandakumar, K.
1999-09-01
Multiphase flow modelling is still a major challenge in fluid dynamics and, although many different models have been derived, there is no clear evidence of their relevance to certain flow situations. That is particularly valid for bubbly flows, because most of the studies have considered the case of fluidized beds. In the present study we give a general formulation to five existing models and study their relevance to bubbly flows. The results of the linear analysis of those models clearly show that only two of them are applicable to that case. They both show a very similar qualitative linear stability behaviour. In the subsequent asymptotic analysis we derive an equation hierarchy which describes the weakly nonlinear stability of the models. Their qualitative behaviour up to first order with respect to the small parameter is again identical. A permanent-wave solution of the first two equations of the hierarchy is found. It is shown, however, that the permanent-wave (soliton) solution is very unlikely to occur for the most common case of gas bubbles in water. The reason is that the weakly nonlinear equations are unstable due to the low magnitude of the bulk modulus of elasticity. Physically relevant stabilization can eventually be achieved using some available experimental data. Finally, a necessary condition for existence of a fully nonlinear soliton is derived.
Dholakia, Nikhilesh; Turcan, Romeo V.
2013-01-01
A goal of our ongoing research stream is to develop a multidisciplinary metatheory of bubbles. In this viewpoint paper we put forward a typology of bubbles by comparing four types of assets – entertainment, commodities, financial securities (stocks), and housing properties – where bubbles could and...... do form occasionally. Cutting across and comparing such varied asset types provides some rich insights into the nature of bubbles – and offers an inductive way to arrive at the typology of bubbles....
Cox, D. P.
Modeling the Local Bubble is one of those activities fraught with danger. It is very easy to be too naive, to fail to consider the dependence of the model on assumptions about the nearby ambient state, or the likelihood of such a structure. It is similarly easy to become so caught up in the details of the vicinity that it is unclear where to begin a necessarily idealized modeling effort. And finally, it is important to remember that the data we have may in some cases be lying to us, and that we have not yet learned to read their facial expressions quite carefully enough. That said, I've tried in this paper to be helpful to those who may wish to take the risks. I surveyed the very most basic stories that the data seem to tell, and pointed out the standard coincidences that may be telling us a lot about what is happening, but may turn out once again to have been just coincidences. I've described 5 distinct conceptions that in one flavor or another pretty well survey the collection of mental images that have so far been carried by those who've attempted models. One may be right, or something entirely different may be more appropriate. It's at least vital to realize that a conception comes first, followed by a simplified model of details. I've also included a long list of questions directed at observers. Some have partial answers, some one wouldn't know today quite how to approach. But it is a list that students of the soft x-ray background, interstellar absorption lines, possible instrumentation, and the heliosphere may wish to review from time to time, just to see whether they can figure out how to be more helpful. There is another list for modelers, things the models must address, however-so-flimsily if necessary, because there are strong observational constraints (and stronger ones coming) on what can and cannot be present in the local ISM. To that I've added a few remarks concerning x-ray emission coming from beyond the Local Bubble, and another few on how x
Large-Scale Clustering in Bubble Models
Borgani, S
1993-01-01
We analyze the statistical properties of bubble models for the large-scale distribution of galaxies. To this aim, we realize static simulations, in which galaxies are mostly randomly arranged in the regions surrounding bubbles. As a first test, we realize simulations of the Lick map, by suitably projecting the three-dimensional simulations. In this way, we are able to safely compare the angular correlation function implied by a bubbly geometry to that of the APM sample. We find that several bubble models provide an adequate amount of large-scale correlation, which nicely fits that of APM galaxies. Further, we apply the statistics of the count-in-cell moments to the three-dimensional distribution and compare them with available observational data on variance, skewness and kurtosis. Based on our purely geometrical constructions, we find that a well defined hierarchical scaling of higher order moments up to scales $\\sim 70\\hm$. The overall emerging picture is that the bubbly geometry is well suited to reproduce ...
Argonne Bubble Experiment Thermal Model Development II
Buechler, Cynthia Eileen [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2016-07-01
This report describes the continuation of the work reported in “Argonne Bubble Experiment Thermal Model Development”. The experiment was performed at Argonne National Laboratory (ANL) in 2014. A rastered 35 MeV electron beam deposited power in a solution of uranyl sulfate, generating heat and radiolytic gas bubbles. Irradiations were performed at three beam power levels, 6, 12 and 15 kW. Solution temperatures were measured by thermocouples, and gas bubble behavior was observed. This report will describe the Computational Fluid Dynamics (CFD) model that was developed to calculate the temperatures and gas volume fractions in the solution vessel during the irradiations. The previous report described an initial analysis performed on a geometry that had not been updated to reflect the as-built solution vessel. Here, the as-built geometry is used. Monte-Carlo N-Particle (MCNP) calculations were performed on the updated geometry, and these results were used to define the power deposition profile for the CFD analyses, which were performed using Fluent, Ver. 16.2. CFD analyses were performed for the 12 and 15 kW irradiations, and further improvements to the model were incorporated, including the consideration of power deposition in nearby vessel components, gas mixture composition, and bubble size distribution. The temperature results of the CFD calculations are compared to experimental measurements.
Modeling of bubble dynamics in relation to medical applications
In various pulsed-laser medical applications, strong stress transients can be generated in advance of vapor bubble formation. To better understand the evolution of stress transients and subsequent formation of vapor bubbles, two-dimensional simulations are presented in channel or cylindrical geometry with the LATIS (LAser TISsue) computer code. Differences with one-dimensional modeling are explored, and simulated experimental conditions for vapor bubble generation are presented and compared with data. 22 refs., 8 figs
Modeling bubble collapse aggressiveness in traveling bubble cavitation using bubble breakup model
Zima, Patrik; Sedlář, M.
Singapore : Research Publishing Services, 2012 - (Ohl, C.; Klaseboer, E.; Ohl, S.; Gong, S.; Khoo, B.), s. 182-186 ISBN 978-981-07-2826-7. [International Symposium on Cavitation /8./ CAV 2012. Singapur (SG), 13.08.2012-16.08.2012] R&D Projects: GA ČR GAP101/10/1428 Institutional research plan: CEZ:AV0Z20760514 Keywords : bubble collapse * traveling bubble cavitation * numerical modeling Subject RIV: BK - Fluid Dynamics http://rpsonline.com.sg/proceedings/9789810728267/html/209.xml51.xml
Bubbles in Non-Newtonian Fluids: A Multiscale Modeling
Frank X.
2013-06-01
Full Text Available In this paper, the concept of a multiscale modeling approach is highlighted with which physical phenomena at different scales can be studied. The work reports a multiscale approach to describe the dynamics of a chain of bubbles rising in non-Newtonian fluids. By means of the Particle Image Velocimetry (PIV and the Lattice Boltzmann (LB simulation, a deep understanding of the complex flow pattern around a single bubble is gained at microscale. The interactions and coalescences between bubbles rising in non-Newtonian fluids are experimentally investigated by the PIV measurements, birefringence and rheological characterization for both an isolated bubble and a chain of bubbles formed from a submerged orifice. Two aspects are identified as central to interactions and coalescence: the stress creation by the passage of bubbles and their relaxation due to the fluid’s memory. This competition between the creation and relaxation of stresses displays non-linear complex dynamics. Along with the detailed knowledge around a single bubble, these fundamental mechanisms governing bubbles’ collective behavior in a train of bubbles at mesoscale lead to a cognitive modeling based on behavioral rules. By simulating bubbles as adaptive agents with the surround fluid via residual stresses, model predictions for consecutive coalescence between a great number of bubbles compare very satisfactorily with the experimental investigation at macroscale. Obviously this new approach captures important quantitative and qualitative features of the collective behaviors of bubbles at macroscale level which are predicted by the mesoscopic cognitive modeling approach of the interactions rules which are deduced from the understanding of the microscopic mechanism of the flow around a single bubble.
Bubble models, data acquisition and model applicability
Jebavá, Marcela; Kloužek, Jaroslav; Němec, Lubomír
Vsetín : GLASS SERVICE ,INC, 2005, s. 182-191. ISBN 80-239-4687-0. [International Seminar on Mathematical Modeling and Advanced Numerical Methods in Furnace Design and Operation /8./. Velké Karlovice (CZ), 19.05.2005-20.05.2005] Institutional research plan: CEZ:AV0Z40320502 Keywords : bubble models Subject RIV: CA - Inorganic Chemistry
IMPROVEMENT OF BUBBLE MODEL FOR CAVITATING FLOW SIMULATIONS
TAMURA Y.; MATSUMOTO Y.
2009-01-01
In the present research,a bubble dynamics based model for cavitating flow simulations is extended to higher void fraction region for wider range of applications.The present bubble model is based on the so-called Rayleigh-Plesset equation that calculates a temporal bubble radius with the surrounding liquid pressure and is considered to be valid in an area below a certain void fraction.The solution algorithm is modified so that the Rayleigh-Plesset equation is no more solved once the bubble radius(or void fraction)reaches at a certain value till the liquid pressure recovers above the vapor pressure in order to overcome this problem.This procedure is expected to stabilize the numerical calculation.The results of simple two-dimensional flow field are presented compared with the existing bubble model.
Modeling bubbly-cap flows using two-group average bubble number density
The basic concept of two-group average bubble number density equations along with three-fluid model has been demonstrated for vertical gas-liquid flow. Specifically, the current study focused on: (i)classification of bubble interaction between spherical bubbles (Group-1) and cap bubbles (Group- 2), (ii) preliminary consideration of source and sink terms in the averaged bubble number density equations via the model of Hibiki and Ishii [1] and (iii) assessment by means of experimental data sets at bubbly-to-cap flow transition. Reasonable agreement was achieved between measured and predicted distributions of void fraction, interfacial area concentration (IAC) and volume equivalent bubble diameter. (author)
The parton bubble model compared to central Au Au collisions (0% to 5%) at $\\sqrt{s_{NN}}$=200 GeV
Longacre, R S
2010-01-01
In an earlier paper we developed a Parton Bubble Model (PBM) for RHIC, high-energy heavy-ion collisions. PBM was based on a substructure of a ring of localized bubbles (gluonic hot spots) which initially contain 3-4 partons composed of almost entirely gluons. The bubbles emitted correlated particles at kinetic freezeout, leading to a lumpy fireball surface. For a selection of charged particles (0.8 GeV/c $<$ $p_t$ $<$ 4.0 GeV/c), the PBM reasonably quantitatively (within a few percent) explained high precision RHIC experimental correlation analyses in a manner which was consistent with the small observed HBT source size in this transverse momentum range. We demonstrated that surface emission from a distributed set of surface sources (as in the PBM) was necessary to obtain this consistency. In this paper we give a review of the above comparison to central Au Au collisions. The bubble formation can be associated with gluonic objects predicted by a Glasma Flux Tube Model (GFTM) that formed longitudinal flu...
Modeling and simulation of bubbles and particles
Dorgan, Andrew James
negligible (in terms of particle concentration predictions) when terminal velocity was oriented in the wall-normal direction. The history force was shown to damp particle diffusion and have some minor impacts on particle concentration. This effect was augmented by using the creeping flow Basset expression and shows that the creeping flow expression should not be used in finite-Reynolds number conditions. The effects due to the finite-size extensions are also considered as are effects due to spatial reconstruction of the fluid properties. In general, little effect of the finite-size model or choice of spatial interpolation was observed in terms of particle concentration. However, Lagrangian statistics show some interesting sensitivities. Finally, the new equation of motion was applied to air bubbles and sand particles of several different diameters. Particle-fluid interactions observed through flow-visualization, particle concentration, particle-wall interactions, and Lagrangian statistics were all considered. These results were interpreted and compared to heavy-particle results where appropriate. Particle deposition was found to be well-described by the heavy-particle model of Young & Leeming and root-mean-square relative velocities were found to also agree with previous heavy-particle work. A model for the latter is suggested for heavy-particles and found to work similarly well for low-density particles. Non-tracer behavior was observed for bubbles with small Stokes numbers, a result not expected based on heavy-particle expectations. Local clustering of particles was observed in certain fluid structures which may indicate the importance of modeling particle collisions in future studies.
Modeling of Wall Effects on Drag and Lift Acting on Bubbles in Bubbly Flow
The two-fluid model based on Eulerian-Eulerian approach has been widely used for simulating two-phase flows in many industrial applications. However, the two-fluid approach needs accurate modeling for interfacial momentum exchange such as drag, shear induced lift, and wall-induced lift. In particular, it is important to accurately model the wall effect in order to predict 'wall peaking' or 'core peaking' phenomena observed in bubbly pipe flow. Those phenomena characterized by the radial distribution of void fraction are mainly determined by the balance between shear induced lift and wall-induced lift usually called 'wall lubrication force'. However, the wall effect is not fully understood yet and the wall force coefficient in previous studies has a wide range of values, probably tuned to the experimental results. Therefore, we propose a new model considering the wall effect on drag and lift forces and evaluate its accuracy by simulating laminar bubbly flows with available data for comparison. In this study, we proposed a new model for interfacial momentum exchange for wall-bounded bubbly flow. In particular, to accurately consider the wall effects on drag and lift, separate simulations were performed for the flow around a moving sphere near the wall. The present model was verified by solving the laminar bubbly flow in a vertical pipe and comparing the results with previous ones. The present void fraction and water velocity profiles showed good results
K. Ekambara
2012-01-01
Full Text Available Modelling of gas-liquid bubbly flows is achieved by coupling a population balance equation with the three-dimensional, two-fluid, hydrodynamic model. For gas-liquid bubbly flows, an average bubble number density transport equation has been incorporated in the CFD code CFX 5.7 to describe the temporal and spatial evolution of the gas bubbles population. The coalescence and breakage effects of the gas bubbles are modeled. The coalescence by the random collision driven by turbulence and wake entrainment is considered, while for bubble breakage, the impact of turbulent eddies is considered. Local spatial variations of the gas volume fraction, interfacial area concentration, Sauter mean bubble diameter, and liquid velocity are compared against experimental data in a horizontal pipe, covering a range of gas (0.25 to 1.34 m/s and liquid (3.74 to 5.1 m/s superficial velocities and average volume fractions (4% to 21%. The predicted local variations are in good agreement with the experimental measurements reported in the literature. Furthermore, the development of the flow pattern was examined at three different axial locations of L/D = 25, 148, and 253. The first location is close to the entrance region where the flow is still developing, while the second and the third represent nearly fully developed bubbly flow patterns.
Hydrodynamic models for slurry bubble column reactors
Gidaspow, D. [IIT Center, Chicago, IL (United States)
1995-12-31
The objective of this investigation is to convert a {open_quotes}learning gas-solid-liquid{close_quotes} fluidization model into a predictive design model. This model is capable of predicting local gas, liquid and solids hold-ups and the basic flow regimes: the uniform bubbling, the industrially practical churn-turbulent (bubble coalescence) and the slugging regimes. Current reactor models incorrectly assume that the gas and the particle hold-ups (volume fractions) are uniform in the reactor. They must be given in terms of empirical correlations determined under conditions that radically differ from reactor operation. In the proposed hydrodynamic approach these hold-ups are computed from separate phase momentum balances. Furthermore, the kinetic theory approach computes the high slurry viscosities from collisions of the catalyst particles. Thus particle rheology is not an input into the model.
无
2001-01-01
A model for a bubble column slurry reactor is developed based on the experiment of Rhenpreussen Koppers demonstration plant for slurry phase Fischer-Tropsch synthesis reported by Koelble et al. This model is applicable to the operation in the churn-turbulent regime and incorporates the information on the bubble size. The axial dispersion model is adopted to describe the flow characteristics of the Fischer-Tropsch slurry reactor. With the model developed, simulations are performed to identify the steady state behavior of a Fischer-Tropsch slurry reactor of commercial size. Predictions of the two-bubble class model is compared with that of the conventional single- bubble class model. The results show that under a variety of conditions, the two-bubble class model gives results different from those for the single-bubble class model.
Modelling of boiling bubbly flows using a polydisperse approach
The objective of this work was to improve the modelling of boiling bubbly flows.We focused on the modelling of the polydisperse aspect of a bubble population, i.e. the fact that bubbles have different sizes and different velocities. The multi-size aspect of a bubble population can originate from various mechanisms. For the bubbly flows we are interested in, bubble coalescence, bubble break-up, phase change kinematics and/or gas compressibility inside the bubbles can be mentioned. Since, bubble velocity depends on bubble size, the bubble size spectrum also leads to a bubble velocity spectrum. An averaged model especially dedicated to dispersed flows is introduced in this thesis. Closure of averaged interphase transfer terms are written in a polydisperse framework, i.e. using a distribution function of the bubble sizes and velocities. A quadratic law and a cubic law are here proposed for the modelling of the size distribution function, whose evolution in space and time is then obtained with the use of the moment method. Our averaged model has been implemented in the NEPTUNE-CFD computation code in order to simulate the DEBORA experiment. The ability of our model to deal with sub-cooled boiling flows has therefore been evaluated. (author)
Free Surface Lattice Boltzmann with Enhanced Bubble Model
Anderl, Daniela; Rauh, Cornelia; Rüde, Ulrich; Delgado, Antonio
2016-01-01
This paper presents an enhancement to the free surface lattice Boltzmann method (FSLBM) for the simulation of bubbly flows including rupture and breakup of bubbles. The FSLBM uses a volume of fluid approach to reduce the problem of a liquid-gas two-phase flow to a single-phase free surface simulation. In bubbly flows compression effects leading to an increase or decrease of pressure in the suspended bubbles cannot be neglected. Therefore, the free surface simulation is augmented by a bubble model that supplies the missing information by tracking the topological changes of the free surface in the flow. The new model presented here is capable of handling the effects of bubble breakup and coalesce without causing a significant computational overhead. Thus, the enhanced bubble model extends the applicability of the FSLBM to a new range of practically relevant problems, like bubble formation and development in chemical reactors or foaming processes.
A computational model of gas bubble evolution in liquid filled straight tubes
Himm, Jeff; Halpern, David
1996-11-01
Deep sea divers suffer from decompression sickness (DCS) when their rate of ascent to the surface is too quick. When the ambient pressure drops, inert gas bubbles are usually formed in blood vessels and tissues of divers. It is believed that the existence of gas bubbles is the cause of DCS that manifests itself as itching, joint pain, and neurological abnormalities. While models of gas bubbles in tissues are relatively well developed, the mechanism of bubble growth in the circulation is far less well understood. The existence of gas bubbles may affect gas exchange in small blood vessels by blocking the flow of blood. Gas bubble evolution in the circulation is investigated using an analytical method for small bubbles and the boundary element method for bubbles whose effective radius is close to the tube radius. The concentration field for the dissolved gas surrounding the bubble is solved numerically using finite differences. The bubble volume is adjusted over time according to the mass flux at the surface. It is shown that the effect of increasing the flow rate is to enhance bubble evolution, up to a factor of two compared with the evolution in tissue where there is no flow. This work was supported by the Naval Medical Research and Development Command work unit 62233N.MM33P30.0041509.
Photon Number Conserving Models of HII Bubbles during Reionization
Paranjape, Aseem; Choudhury, T. Roy; Padmanabhan, Hamsa
2016-05-01
Traditional excursion set based models of HII bubble growth during the epoch of reionization are known to violate photon number conservation, in the sense that the mass fraction in ionized bubbles in these models does not equal the ratio of the number of ionizing photons produced by sources and the number of hydrogen atoms in the intergalactic medium. E.g., for a Planck13 cosmology with electron scattering optical depth τ ≃ 0.066, the discrepancy is ˜15 per cent for xHII = 0.1 and ˜5 per cent for xHII = 0.5. We demonstrate that this problem arises from a fundamental conceptual shortcoming of the excursion set approach (already recognised in the literature on this formalism) which only tracks average mass fractions instead of the exact, stochastic source counts. With this insight, we build an approximately photon number conserving Monte Carlo model of bubble growth based on partitioning regions of dark matter into halos. Our model, which is formally valid for white noise initial conditions (ICs), shows dramatic improvements in photon number conservation, as well as substantial differences in the bubble size distribution, as compared to traditional models. We explore the trends obtained on applying our algorithm to more realistic ICs, finding that these improvements are robust to changes in the ICs. Since currently popular semi-numerical schemes of bubble growth also violate photon number conservation, we argue that it will be worthwhile to pursue new, explicitly photon number conserving approaches. Along the way, we clarify some misconceptions regarding this problem that have appeared in the literature.
Studies on modelling of bubble driven flows in chemical reactors
Grevskott, Sverre
1997-12-31
Multiphase reactors are widely used in the process industry, especially in the petrochemical industry. They very often are characterized by very good thermal control and high heat transfer coefficients against heating and cooling surfaces. This thesis first reviews recent advances in bubble column modelling, focusing on the fundamental flow equations, drag forces, transversal forces and added mass forces. The mathematical equations for the bubble column reactor are developed, using an Eulerian description for the continuous and dispersed phase in tensor notation. Conservation equations for mass, momentum, energy and chemical species are given, and the k-{epsilon} and Rice-Geary models for turbulence are described. The different algebraic solvers used in the model are described, as are relaxation procedures. Simulation results are presented and compared with experimental values. Attention is focused on the modelling of void fractions and gas velocities in the column. The energy conservation equation has been included in the bubble column model in order to model temperature distributions in a heated reactor. The conservation equation of chemical species has been included to simulate absorption of CO{sub 2}. Simulated axial and radial mass fraction profiles for CO{sub 2} in the gas phase are compared with measured values. Simulations of the dynamic behaviour of the column are also presented. 189 refs., 124 figs., 1 tab.
Modelling for three dimensional coalescence of two bubbles
Han, R.; Li, S.; Zhang, A. M.; Wang, Q. X.
2016-06-01
This paper is concerned with the three dimensional (3D) interaction and coalescence of two bubbles subject to buoyancy and the dynamics of the subsequent joined bubble using the boundary integral method (BIM). An improved density potential method is implemented to control the mesh quality. It helps to avoid the numerical instabilities, which occur after coalescence. Numerical convergence tests are conducted in terms of mesh sizes and time steps. The 3D numerical model agrees well with an axisymmetric BIM model for axisymmetric cases as well as experimental results captured by high-speed camera. The bubble jetting, interaction, and coalescence of the two bubbles depend on the maximum bubble radii, the centre distance between two bubbles at inception, and the angle β between the centre line and the direction of buoyancy. We investigate coalescence of two bubbles for β = 0, π/4, and π/2, respectively, and at various centre distances at inception. Numerical results presented include the bubble and jet shapes, the velocity, and pressure fields surrounding the bubbles, as well as the time histories of bubble volumes, jet velocities, and positions of centroid of the bubble system.
A mechanistic model of separation bubble
Krechetnikov, R; Nagib, H M
2007-01-01
This work uncovers the low-dimensional nature the complex dynamics of actuated separated flows. Namely, motivated by the problem of model-based predictive control of separated flows, we identify the requirements on a model-based observer and the key variables and propose a prototype model in the case of thick airfoils as motivated by practical applications. The approach in this paper differs fundamentally from the logic behind known models, which are either linear or based on POD-truncations and are unable to reflect even the crucial bifurcation and hysteresis inherent in separation phenomena. This new look at the problem naturally leads to several important implications, such as, firstly, uncovering the physical mechanisms for hysteresis, secondly, predicting a finite amplitude instability of the bubble, and thirdly to new issues to be studied theoretically and tested experimentally. More importantly, by employing systematic reasoning, the low-dimensional nature of these complex phenomena at the coarse level...
Force Balance Model for Bubble Rise, Impact, and Bounce from Solid Surfaces.
Manica, Rogerio; Klaseboer, Evert; Chan, Derek Y C
2015-06-23
A force balance model for the rise and impact of air bubbles in a liquid against rigid horizontal surfaces that takes into account effects of buoyancy and hydrodynamic drag forces, bubble deformation, inertia of the fluid via an added mass force, and a film force between the bubble and the rigid surface is proposed. Numerical solution of the governing equations for the position and velocity of the center of mass of the bubbles is compared against experimental data taken with ultraclean water. The boundary condition at the air-water interface is taken to be stress free, which is consistent for bubbles in clean water systems. Features that are compared include bubble terminal velocity, bubbles accelerating from rest to terminal speed, and bubbles impacting and bouncing off different solid surfaces for bubbles that have already or are yet to attain terminal speed. Excellent agreement between theory and experiments indicates that the forces included in the model constitute the main physical ingredients to describe the bouncing phenomenon. PMID:26035016
A new pressure formulation for gas-compressibility dampening in bubble dynamics models.
Gadi Man, Yezaz Ahmed; Trujillo, Francisco J
2016-09-01
We formulated a pressure equation for bubbles performing nonlinear radial oscillations under ultrasonic high pressure amplitudes. The proposed equation corrects the gas pressure at the gas-liquid interface on inertial bubbles. This pressure formulation, expressed in terms of gas-Mach number, accounts for dampening due to gas compressibility during the violent collapse of cavitation bubbles and during subsequent rebounds. We refer to this as inhomogeneous pressure, where the gas pressure at the gas-liquid interface can differ to the pressure at the centre of the bubble, in contrast to homogenous pressure formulations that consider that pressure inside the bubble is spatially uniform from the wall to the centre. The pressure correction was applied to two bubble dynamic models: the incompressible Rayleigh-Plesset equation and the compressible Keller and Miksis equation. This improved the predictions of the nonlinear radial motion of the bubble vs time obtained with both models. Those simulations were also compared with other bubble dynamics models that account for liquid and gas compressibility effects. It was found that our corrected models are in closer agreement with experimental data than alternative models. It was concluded that the Rayleigh-Plesset family of equations improve accuracy by using our proposed pressure correction. PMID:27150768
THE SEMIEMPIRICAL MODEL OF THE MULTICOMPONENT BUBBLE BEHAVIOUR IN GLASS MELTS
LUBOMÍR NĚMEC
2012-12-01
Full Text Available A semi-empirical model of the bubble growth and dissolution in glasses with a fining agent has been derived. This model applies the experimental data from bubble observation at melting and fining temperatures. The experimental data needed for the model involved the temperature dependences of the average growth rate of the bubble radius and the average concentration of the fining gas in the bubbles. Both sets of values were measured in the laboratory in the glass of the float type and applied in the model. The measurements of the solubilities and diffusion coefficients of the gases present in the glass – needed for the analytical model of multicomponent bubbles – were thus avoided. The course of the partial bubble absorption with the temperature decreasing was simulated by means of two factors modifying the experimental values of the bubble growth rates at constant temperature. The temperature dependence of the resulting bubble growth rate qualitatively corresponded to the experimental observations in the soda-lime-silica glass, but a more detailed experimental and comparative study has yet to be performed. Such a study is being prepared.
Extension of the inhomogeneous MUSIG model for bubble condensation
Highlights: ► The inhomogenous MUSIG model allows 3D simulations for poly-dispersed bubbly flows. ► The model is now extended to consider flows with phase transfer. ► Experimental data for the condensation of steam bubbles in sub-cooled vertical pipe flow are used for validation. ► There is a good agreement between experimental data and CFD simulations with the ANSYS-CFX code. - Abstract: Bubble condensation plays an important role, e.g. in sub-cooled boiling or steam injection into pools. Since the condensation rate is proportional to the interfacial area density, bubble size distributions have to be considered in an adequate modeling of the condensation process. The effect of bubble sizes was clearly shown in experimental investigations done previously at the TOPFLOW facility of FZD. Steam bubbles were injected into a sub-cooled upward pipe flow via orifices in the pipe wall located at different distances from measuring plane. 1 mm and 4 mm injection orifices were used to vary the initial bubble size distribution. Measurements were done using a wire-mesh sensor. Condensation is clearly faster in case of the injection via the smaller orifices, i.e. in case of smaller bubble sizes. Recently the Inhomogeneous MUSIG model was implemented into the CFD code CFX from ANSYS enabling the simulation of poly-dispersed flows including the effects of separation of small and large bubbles due to bubble size dependent lift force inversion. It allows to divide the dispersed phase into size classes regarding the mass as well as regarding the momentum balance. Up to now transfers between the classes in the mass balance can be considered only by bubble coalescence and breakup (population balance). Here an extension of the model is proposed to include the effects due to phase transfer. The paper focuses on the derivation of equations for the extension of the Inhomogeneous MUSIG model and presents some first results for verification and validation.
The effect of high viscosity on compressible and incompressible Rayleigh–Plesset-type bubble models
Highlights: • Spherical gas bubbles in glycerol have been examined numerically and experimentally. • The bubble was generated using a Q-switched Nd:YAG laser. • The radius was measured with a novel shadowing technique of a He–Ne laser beam. • The measurements were compared with a compressible and an incompressible bubble models. • The validity domain of the incompressible assumption has been given. -- Abstract: Free oscillations of a single spherical gas bubble in glycerol have been examined numerically and experimentally at different ambient temperatures and pressures. The bubble was generated using a Q-switched Nd:YAG laser and the unsteady radius measurement was based on a shadowing technique of a He–Ne laser beam. The measurements were compared to computations obtained from two models, first taking into consideration the liquid compressibility and then assuming an incompressible liquid domain, respectively. In both cases the temperature fields inside and outside the bubble were computed by solving the energy equation in both phases as the thermodynamic processes have great importance to the bubble behavior. For high amplitude oscillations the incompressible model provides poor agreement with the measurements and the modeling of the liquid compressibility becomes necessary. In contrast to the standard method, a practical region of applicability for the incompressible approach was determined as a function of the instantaneous Mach and Reynolds numbers, rather than specifying a simple threshold Mach number
In previous theoretical studies of the behaviour of the fission gases in nuclear fuel, the Nelson single knock-on model of the fission induced re-solution of gas atoms from fission gas bubbles has been employed. In the present investigation, predictions from this model are compared with those from a complete bubble destruction model of the re-solution process. The main conclusions of the study are that the complete bubble destruction model predicts more gas release after a particular irradiation time than the single knock-on model, for the same choice of the model parameters, and that parameter sets chosen to give the same gas release predict significantly different bubble size distribution functions. (author)
Li, Linmin; Li, Baokuan
2016-03-01
In ladle metallurgy, bubble-liquid interaction leads to complex phase structures. Gas bubble behavior, as well as the induced slag layer behavior, plays a significant role in the refining process and the steel quality. In the present work, a mathematical model using the large eddy simulation (LES) is developed to investigate the bubble transport and slag layer behavior in a water model of an argon-stirred ladle. The Eulerian volume of fluid model is adopted to track the liquid steel-slag-air free surfaces while the Lagrangian discrete phase model is used for tracking and handling the dynamics of discrete bubbles. The bubble coalescence is considered using O'Rourke's algorithm to solve the bubble diameter redistribution and bubbles are removed after leaving the air-liquid interface. The turbulent liquid flow that is induced by bubble-liquid interaction is solved by LES. The slag layer fluactuation, slag droplet entrainment and spout eye open-close phenomenon are well revealed. The bubble diameter distribution and the spout eye size are compared with the experiment. The results show that the hybrid Eulerian-Lagrangian-LES model provides a valid modeling framework to predict the unsteady gas bubble-slag layer coupled behaviors.
Under decompression, bubbles can form in the human body, and these can be found both within the body tissues and the bloodstream. Mathematical models for the growth of both types of bubbles have previously been presented, but they have not been coupled together. This work thus explores the interaction between the growth of tissue and blood-borne bubbles under decompression, specifically looking at the extent to which they compete for the common resource of inert gas held in solution in the tissues. The influence of tissue bubbles is found to be significant for densities as low as 10 ml-1 for tissues which are poorly perfused. However, the effects of formation of bubbles in the blood are not found until the density of bubble production sites reaches 106 ml-1. From comparison of the model predictions with experimental evidence for bubbles produced in animals and man under decompression, it is concluded that the density of tissue bubbles is likely to have a significant effect on the number of bubbles produced in the blood. However, the density of nucleation sites in the blood is unlikely to be sufficiently high in humans for the formation of bubbles in the blood to have a significant impact on the growth of the bubbles in the tissue
Bubble Size Models for the Prediction of Bubbly Flow with CMFD Code
Bak, Jin-yeong; Yun, Byong-jo; Jeong, Jae-jun [Pusan National Univ., Busan (Korea, Republic of)
2015-05-15
In recent years, the use of computational multi-fluid dynamics (CMFD) codes has been extended to the analysis of multi-dimensional two-phase flow for the operation and safety analysis of nuclear power plants (NPP). In these applications, an accurate prediction of bubble behaviors is one of major concerns. Yao and Morel and Yeoh and Tu respectively applied interfacial area concentration transport (IACT) equation and bubble number density transport equation into CMFD code. Recently Lo and Zhang tried to apply the generalized S{sub γ} model to the predictions of not only droplet size in the oil-water flow but also bubble size in the air-water flow. In this paper, three-dimensional numerical simulations for the gas-liquid two-phase flow were conducted to validate and confirm the performance of S{sub γ} bubble size model for the further application to the narrow rectangular boiling channel for the research reactor core, using the commercial CFD code STAR CCM''+ ver. 9.06. For this, S{sub γ} model was evaluated against air-water data of DEDALE and Hibiki et al.'s experiment. These experimental data were obtained in a vertically arranged pipe under upwards air-water flow condition. Detailed descriptions on the S{sub γ} with its breakup and coalescence model are presented in the present manuscript.
A 3D Bubble Merger Model for RTI Mixing
Cheng, Baolian
2015-11-01
In this work we present a model for the merger processes of bubbles at the edge of an unstable acceleration driven mixing layer. Steady acceleration defines a self-similar mixing process, with a time-dependent inverse cascade of structures of increasing size. The time evolution is itself a renormalization group evolution. The model predicts the growth rate of a Rayleigh-Taylor chaotic fluid-mixing layer. The 3-D model differs from the 2-D merger model in several important ways. Beyond the extension of the model to three dimensions, the model contains one phenomenological parameter, the variance of the bubble radii at fixed time. The model also predicts several experimental numbers: the bubble mixing rate, the mean bubble radius, and the bubble height separation at the time of merger. From these we also obtain the bubble height to the radius aspect ratio, which is in good agreement with experiments. Applications to recent NIF and Omega experiments will be discussed. This work was performed under the auspices of the U.S. Department of Energy by the Los Alamos National Laboratory under Contract No. W-7405-ENG-36.
Time-Dependent Stochastic Acceleration Model for the Fermi Bubbles
Sasaki, Kento; Terasawa, Toshio
2015-01-01
We study stochastic acceleration models for the Fermi bubbles. Turbulence is excited just behind the shock front via Kelvin-Helmholtz, Rayleigh-Taylor or Richtmyer-Meshkov instabilities, and plasma particles are continuously accelerated by the interaction with the turbulence. The turbulence gradually decays as it goes away from the shock fronts. Adopting a phenomenological model for the stochastic acceleration, we explicitly solve the temporal evolution of the particle energy distribution in the turbulence. Our results show that the spatial distribution of high-energy particles is different from those for a steady solution. We also show that the contribution of electrons escaped from the acceleration regions significantly softens the photon spectrum. The photon spectrum and surface brightness profile are reproduced by our models. If the escape efficiency is very high, the radio flux from the escaped low-energy electrons can be comparable to that of the WMAP haze. We also demonstrate hadronic models with the s...
Modeling the dynamics of single-bubble sonoluminescence
Vignoli, Lucas L; Thomé, Roberto C A; Nogueira, A L M A; Paschoal, Ricardo C; Rodrigues, Hilario
2014-01-01
Sonoluminescence (SL) is the phenomenon in which acoustic energy is (partially) transformed into light. It may occur by means of many or just one bubble of gas inside a liquid medium, giving rise to the terms multi-bubble- and single-bubble sonoluminescence (MBSL and SBSL). In the last years some models have been proposed to explain this phenomenon, but there is still no complete theory for the light emission mechanism (especially in the case of SBSL). In this work, we will not address this more complicated particular issue, but only present a simple model describing the dynamical behaviour of the sonoluminescent bubble, in the SBSL case. Using simple numerical techniques within the software Matlab, we discuss solutions considering various possibilities for some of the parameters involved: liquid compressibility, superficial tension, viscosity, and type of gas. The model may be used as an introductory study of sonoluminescence in physics courses at undergraduate or graduate levels, as well as a quite clarifyi...
Experimental investigation and mechanistic modelling of dilute bubbly bulk boiling
During evaporation the geometric shape of the vapour is not described using thermodynamics. In bubbly flows the bubble shape is considered spheric with small diameters and changing into various shapes upon growth. The heat and mass transfer happens at the interfacial area. The forces acting on the bubbles depend on the bubble diameter and shape. In this work the prediction of the bubble diameter and/or bubble number density in bulk boiling was considered outside the vicinity of the heat input area. Thus the boiling effects that happened inside the nearly saturated bulk were under investigation. This situation is relevant for nuclear safety analysis concerning a stagnant coolant in the spent fuel pool. In this research project a new experimental set-up to investigate was built. The experimental set-up consists of an instrumented, partly transparent, high and slender boiling container for visual observation. The direct visual observation of the boiling phenomena is necessary for the identification of basic mechanisms, which should be incorporated in the simulation model. The boiling process has been recorded by means of video images and subsequently was evaluated by digital image processing methods, and by that data concerning the characteristics of the boiling process were generated for the model development and validation. Mechanistic modelling is based on the derivation of relevant mechanisms concluded from observation, which is in line with physical knowledge. In this context two mechanisms were identified; the growth/-shrink mechanism (GSM) of the vapour bubbles and sudden increases of the bubble number density. The GSM was implemented into the CFD-Code ANSYS-CFX using the CFX Expression Language (CEL) by calculation of the internal bubble pressure using the Young-Laplace-Equation. This way a hysteresis is realised as smaller bubbles have an increased internal pressure. The sudden increases of the bubble number density are explainable by liquid super
Experimental investigation and mechanistic modelling of dilute bubbly bulk boiling
Kutnjak, Josip
2013-06-27
During evaporation the geometric shape of the vapour is not described using thermodynamics. In bubbly flows the bubble shape is considered spheric with small diameters and changing into various shapes upon growth. The heat and mass transfer happens at the interfacial area. The forces acting on the bubbles depend on the bubble diameter and shape. In this work the prediction of the bubble diameter and/or bubble number density in bulk boiling was considered outside the vicinity of the heat input area. Thus the boiling effects that happened inside the nearly saturated bulk were under investigation. This situation is relevant for nuclear safety analysis concerning a stagnant coolant in the spent fuel pool. In this research project a new experimental set-up to investigate was built. The experimental set-up consists of an instrumented, partly transparent, high and slender boiling container for visual observation. The direct visual observation of the boiling phenomena is necessary for the identification of basic mechanisms, which should be incorporated in the simulation model. The boiling process has been recorded by means of video images and subsequently was evaluated by digital image processing methods, and by that data concerning the characteristics of the boiling process were generated for the model development and validation. Mechanistic modelling is based on the derivation of relevant mechanisms concluded from observation, which is in line with physical knowledge. In this context two mechanisms were identified; the growth/-shrink mechanism (GSM) of the vapour bubbles and sudden increases of the bubble number density. The GSM was implemented into the CFD-Code ANSYS-CFX using the CFX Expression Language (CEL) by calculation of the internal bubble pressure using the Young-Laplace-Equation. This way a hysteresis is realised as smaller bubbles have an increased internal pressure. The sudden increases of the bubble number density are explainable by liquid super
ONE-DIMENSIONAL DYNAMICAL MODELS OF THE CARINA NEBULA BUBBLE
We have tested the two main theoretical models of bubbles around massive star clusters, Castor et al. and Chevalier and Clegg, against observations of the well-studied Carina nebula. The Castor et al. theory overpredicts the X-ray luminosity in the Carina bubble by a factor of 60 and expands too rapidly, by a factor of 4; if the correct radius and age are used, the predicted X-ray luminosity is even larger. In contrast, the Chevalier and Clegg model underpredicts the X-ray luminosity by a factor of 10. We modify the Castor et al. theory to take into account lower stellar wind mass-loss rates, radiation pressure, gravity, and escape of or energy loss from the hot shocked gas. We argue that energy is advected rather than radiated from the bubble. We undertake a parameter study for reduced stellar mass-loss rates and for various leakage rates and are able to find viable models. The X-ray surface brightness in Carina is highest close to the bubble wall, which is consistent with conductive evaporation from cold clouds. The picture that emerges is one in which the hot gas pressure is far below that found by dividing the time-integrated wind luminosity by the bubble volume; rather, the pressure in the hot gas is set by pressure equilibrium with the photoionized gas at T = 104 K. It follows that the shocked stellar winds are not dynamically important in forming the bubbles.
Numerical simulations of bubbly flows using an averaged equations' model
The paper discusses a numerical method for solving a two phase flow model based on the interpenetrating continua hypothesis. The model incorporates terms to account for the effects of virtual mass force, different pressures for the two phases and the viscous dissipation. Our numerical scheme extends the incremental projection scheme for the incompressible Navier-Stokes equation toward the multiphase flows. An optimal stability is obtained by slightly modifying the Galerkin formulation. The stabilized Galerkin technique we used is based on a two-level hierarchical decomposition of the approximation space. Numerical simulations of the three-dimensional bubbly flows in a periodic domain are presented. These simulations are compared with experiments. The stability of this flow with respect to 3D perturbations is studied numerically and a discussion of the results is presented. (author)
Fluid dynamic modelling of bubble column reactors
Khan, Khurram Imran
2014-01-01
Numerical simulations of rectangular shape bubble column reactors (BCR) are validated starting from preliminary simulations aimed at identifying proper simulation parameters for a given system and resulting up to the numerical simulation with mass transfer and chemical reactions. The transient, three dimensional simulations are carried out using FLUENT software and the results obtained for a system with low gas flow rate (48 L/h) indicated that we need enough fine mesh grid and appropriate cl...
Algebraic model for bubble tracking in horizontal gas-liquid flow
Freitas, Felipe G.C. de; Tisserant, Hendy R. [Universidade Tecnologica Federal do Parana (UTFPR), Curitiba, PR (Brazil); Morales, Rigoberto E.M. [Universidade Tecnologica Federal do Parana (UTFPR), Curitiba, PR (Brazil). Programa de Pos-Graduacao em Engenharia Mecanica e de Materiais; Mazza, Ricardo A.; Rosa, Eugenio S. [Universidade Estadual de Campinas (UNICAMP), SP (Brazil). Fac. de Engenharia Mecanica
2008-07-01
The current work extends the concept of unit-cell applied in gas-liquid slug flow models to predict the evolution of the gas and liquid flow properties along a horizontal pipe. The motivation of this model is its simplicity, easiness of application and low computational cost. It is a useful tool of reference data generation in order to check the consistency of numerical slug tracking models. The potential of the model is accessed by comparing the gas bubbles and liquid slug sizes, the translational bubble velocity and the pressure drop against experimental data. (author)
Using a dynamic point-source percolation model to simulate bubble growth
Accurate modeling of nucleation, growth and clustering of helium bubbles within metal tritide alloys is of high scientific and technological importance. Of interest is the ability to predict both the distribution of these bubbles and the manner in which these bubbles interact at a critical concentration of helium-to-metal atoms to produce an accelerated release of helium gas. One technique that has been used in the past to model these materials, and again revisited in this research, is percolation theory. Previous efforts have used classical percolation theory to qualitatively and quantitatively model the behavior of interstitial helium atoms in a metal tritide lattice; however, higher fidelity models are needed to predict the distribution of helium bubbles and include features that capture the underlying physical mechanisms present in these materials. In this work, we enhance classical percolation theory by developing the dynamic point-source percolation model. This model alters the traditionally binary character of site occupation probabilities by enabling them to vary depending on proximity to existing occupied sites, i.e. nucleated bubbles. This revised model produces characteristics for one and two dimensional systems that are extremely comparable with measurements from three dimensional physical samples. Future directions for continued development of the dynamic model are also outlined
Using a dynamic point-source percolation model to simulate bubble growth.
Zimmerman, Jonathan A.; Zeigler, David A.; Cowgill, Donald F.
2004-05-01
Accurate modeling of nucleation, growth and clustering of helium bubbles within metal tritide alloys is of high scientific and technological importance. Of interest is the ability to predict both the distribution of these bubbles and the manner in which these bubbles interact at a critical concentration of helium-to-metal atoms to produce an accelerated release of helium gas. One technique that has been used in the past to model these materials, and again revisited in this research, is percolation theory. Previous efforts have used classical percolation theory to qualitatively and quantitatively model the behavior of interstitial helium atoms in a metal tritide lattice; however, higher fidelity models are needed to predict the distribution of helium bubbles and include features that capture the underlying physical mechanisms present in these materials. In this work, we enhance classical percolation theory by developing the dynamic point-source percolation model. This model alters the traditionally binary character of site occupation probabilities by enabling them to vary depending on proximity to existing occupied sites, i.e. nucleated bubbles. This revised model produces characteristics for one and two dimensional systems that are extremely comparable with measurements from three dimensional physical samples. Future directions for continued development of the dynamic model are also outlined.
Argonne Bubble Experiment Thermal Model Development
Buechler, Cynthia Eileen [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2015-12-03
This report will describe the Computational Fluid Dynamics (CFD) model that was developed to calculate the temperatures and gas volume fractions in the solution vessel during the irradiation. It is based on the model used to calculate temperatures and volume fractions in an annular vessel containing an aqueous solution of uranium . The experiment was repeated at several electron beam power levels, but the CFD analysis was performed only for the 12 kW irradiation, because this experiment came the closest to reaching a steady-state condition. The aim of the study is to compare results of the calculation with experimental measurements to determine the validity of the CFD model.
Numerical modeling of bubble dynamics in viscoelastic media with relaxation
Warnez, M. T.; Johnsen, E.
2015-06-01
Cavitation occurs in a variety of non-Newtonian fluids and viscoelastic materials. The large-amplitude volumetric oscillations of cavitation bubbles give rise to high temperatures and pressures at collapse, as well as induce large and rapid deformation of the surroundings. In this work, we develop a comprehensive numerical framework for spherical bubble dynamics in isotropic media obeying a wide range of viscoelastic constitutive relationships. Our numerical approach solves the compressible Keller-Miksis equation with full thermal effects (inside and outside the bubble) when coupled to a highly generalized constitutive relationship (which allows Newtonian, Kelvin-Voigt, Zener, linear Maxwell, upper-convected Maxwell, Jeffreys, Oldroyd-B, Giesekus, and Phan-Thien-Tanner models). For the latter two models, partial differential equations (PDEs) must be solved in the surrounding medium; for the remaining models, we show that the PDEs can be reduced to ordinary differential equations. To solve the general constitutive PDEs, we present a Chebyshev spectral collocation method, which is robust even for violent collapse. Combining this numerical approach with theoretical analysis, we simulate bubble dynamics in various viscoelastic media to determine the impact of relaxation time, a constitutive parameter, on the associated physics. Relaxation time is found to increase bubble growth and permit rebounds driven purely by residual stresses in the surroundings. Different regimes of oscillations occur depending on the relaxation time.
Using DNS and Statistical Learning to Model Bubbly Channel Flow
Ma, Ming; Lu, Jiacai; Tryggvason, Gretar
2015-11-01
The transient evolution of laminar bubbly flow in a vertical channel is examined by direct numerical simulation (DNS). Nearly spherical bubbles, initially distributed evenly in a fully developed parabolic flow, are driven relatively quickly to the walls, where they increase the drag and reduce the flow rate on a longer time scale. Once the flow rate has been decreased significantly, some of the bubbles move back into the channel interior and the void fraction there approaches the value needed to balance the weight of the mixture and the imposed pressure gradient. A database generated by averaging the DNS results is used to model the closure terms in a simple model of the average flow. Those terms relate the averaged lateral flux of the bubbles, the velocity fluctuations and the averaged surface tension force to the fluid shear, the void fraction and its gradient, as well as the distance to the nearest wall. An aggregated neural network is used for the statistically leaning of unknown closures, and closure relationships are tested by following the evolution of bubbly channel flow with different initial conditions. It is found that the model predictions are in reasonably good agreement with DNS results. Supported by NSF.
CFD Approaches for Modelling Bubble Entrainment by an Impinging Jet
Martin Schmidtke
2009-01-01
Full Text Available This contribution presents different approaches for the modeling of gas entrainment under water by a plunging jet. Since the generation of bubbles happens on a scale which is smaller than the bubbles, this process cannot be resolved in meso-scale simulations, which include the full length of the jet and its environment. This is why the gas entrainment has to be modeled in meso-scale simulations. In the frame of a Euler-Euler simulation, the local morphology of the phases has to be considered in the drag model. For example, the gas is a continuous phase above the water level but bubbly below the water level. Various drag models are tested and their influence on the gas void fraction below the water level is discussed. The algebraic interface area density (AIAD model applies a drag coefficient for bubbles and a different drag coefficient for the free surface. If the AIAD model is used for the simulation of impinging jets, the gas entrainment depends on the free parameters included in this model. The calculated gas entrainment can be adapted via these parameters. Therefore, an advanced AIAD approach could be used in future for the implementation of models (e.g., correlations for the gas entrainment.
Cogné, C; Labouret, S; Peczalski, R; Louisnard, O; Baillon, F; Espitalier, F
2016-03-01
This paper deals with the inertial cavitation of a single gas bubble in a liquid submitted to an ultrasonic wave. The aim was to calculate accurately the pressure and temperature at the bubble wall and in the liquid adjacent to the wall just before and just after the collapse. Two different approaches were proposed for modeling the heat transfer between the ambient liquid and the gas: the simplified approach (A) with liquid acting as perfect heat sink, the rigorous approach (B) with liquid acting as a normal heat conducting medium. The time profiles of the bubble radius, gas temperature, interface temperature and pressure corresponding to the above models were compared and important differences were observed excepted for the bubble size. The exact pressure and temperature distributions in the liquid corresponding to the second model (B) were also presented. These profiles are necessary for the prediction of any physical phenomena occurring around the cavitation bubble, with possible applications to sono-crystallization. PMID:26044460
Modeling DNA bubble formation at the atomic scale
We describe the fluctuations of double stranded DNA molecules using a minimalist Go model over a wide range of temperatures. Minimalist models allow us to describe, at the atomic level, the opening and formation of bubbles in DNA double helices. This model includes all the geometrical constraints in helix melting imposed by the 3D structure of the molecule. The DNA forms melted bubbles within double helices. These bubbles form and break as a function of time. The equilibrium average number of broken base pairs shows a sharp change as a function of T. We observe a temperature profile of sequence dependent bubble formation similar to those measured by Zeng et al. Long nuclei acid molecules melt partially through the formations of bubbles. It is known that CG rich sequences melt at higher temperatures than AT rich sequences. The melting temperature, however, is not solely determined by the CG content, but by the sequence through base stacking and solvent interactions. Recently, models that incorporate the sequence and nonlinear dynamics of DNA double strands have shown that DNA exhibits a very rich dynamics. Recent extensions of the Bishop-Peyrard model show that fluctuations in the DNA structure lead to opening in localized regions, and that these regions in the DNA are associated with transcription initiation sites. 1D and 2D models of DNA may contain enough information about stacking and base pairing interactions, but lack the coupling between twisting, bending and base pair opening imposed by the double helical structure of DNA that all atom models easily describe. However, the complexity of the energy function used in all atom simulations (including solvent, ions, etc) does not allow for the description of DNA folding/unfolding events that occur in the microsecond time scale.
Time-dependent Stochastic Acceleration Model for Fermi Bubbles
Sasaki, Kento; Asano, Katsuaki; Terasawa, Toshio
2015-12-01
We study stochastic acceleration models for the Fermi bubbles. Turbulence is excited just behind the shock front via Kelvin-Helmholtz, Rayleigh-Taylor, or Richtmyer-Meshkov instabilities, and plasma particles are continuously accelerated by the interaction with the turbulence. The turbulence gradually decays as it goes away from the shock fronts. Adopting a phenomenological model for the stochastic acceleration, we explicitly solve the temporal evolution of the particle energy distribution in the turbulence. Our results show that the spatial distribution of high-energy particles is different from those for a steady solution. We also show that the contribution of electrons that escaped from the acceleration regions significantly softens the photon spectrum. The photon spectrum and surface brightness profile are reproduced by our models. If the escape efficiency is very high, the radio flux from the escaped low-energy electrons can be comparable to that of the WMAP haze. We also demonstrate hadronic models with the stochastic acceleration, but they are unlikely in the viewpoint of the energy budget.
Development and validation of models for bubble coalescence and breakup
A generalized model for bubble coalescence and breakup has been developed, which is based on a comprehensive survey of existing theories and models. One important feature of the model is that all important mechanisms leading to bubble coalescence and breakup in a turbulent gas-liquid flow are considered. The new model is tested extensively in a 1D Test Solver and a 3D CFD code ANSYS CFX for the case of vertical gas-liquid pipe flow under adiabatic conditions, respectively. Two kinds of extensions of the standard multi-fluid model, i.e. the discrete population model and the inhomogeneous MUSIG (multiple-size group) model, are available in the two solvers, respectively. These extensions with suitable closure models such as those for coalescence and breakup are able to predict the evolution of bubble size distribution in dispersed flows and to overcome the mono-dispersed flow limitation of the standard multi-fluid model. For the validation of the model the high quality database of the TOPFLOW L12 experiments for air-water flow in a vertical pipe was employed. A wide range of test points, which cover the bubbly flow, turbulent-churn flow as well as the transition regime, is involved in the simulations. The comparison between the simulated results such as bubble size distribution, gas velocity and volume fraction and the measured ones indicates a generally good agreement for all selected test points. As the superficial gas velocity increases, bubble size distribution evolves via coalescence dominant regimes first, then breakup-dominant regimes and finally turns into a bimodal distribution. The tendency of the evolution is well reproduced by the model. However, the tendency is almost always overestimated, i.e. too much coalescence in the coalescence dominant case while too much breakup in breakup dominant ones. The reason of this problem is discussed by studying the contribution of each coalescence and breakup mechanism at different test points. The redistribution of the
Modeling Air Bubble Transport in Hydraulic Jump Flows using Population Balance Approach
Min Xiang
2016-01-01
Full Text Available This paper proposed a numerical model aiming at coupling the MUltiple-SIze-Group (MUSIG with the semiempirical air entrainment model based on the Euler-Euler two-fluid framework to handle the bubble transport in hydraulic jump flows. The internal flow structure including the recirculation region, the shear layer region and the jet region was accurately predicted. The flow parameters such as the water velocity and void fraction distributions were examined and compared with the experimental data, validating the effectiveness of the numerical model. Prediction of the Sauter mean bubble diameter distributions by the population balance approach at different axial locations confirmed the dominance of breakage due to the high turbulent intensity in the shear layer region which led to the generation of small gas bubbles at high void fraction. Comparison between different cases indicates that high Froude number not only give rise to longer recirculation region and higher void fraction due to larger air entrainment rate, but also generate larger bubble number density and smaller bubble size because of the stronger turbulence intensity in the same axial position.
Evolutionary thinking in microeconomic models: prestige bias and market bubbles.
Bell, Adrian Viliami
2013-01-01
Evolutionary models broadly support a number of social learning strategies likely important in economic behavior. Using a simple model of price dynamics, I show how prestige bias, or copying of famed (and likely successful) individuals, influences price equilibria and investor disposition in a way that exacerbates or creates market bubbles. I discuss how integrating the social learning and demographic forces important in cultural evolution with economic models provides a fruitful line of inquiry into real-world behavior. PMID:23544100
A modelling and experimental study of the bubble trajectory in a non-Newtonian crystal suspension
Hassan, N M S [Process Engineering and Light Metals (PELM) Centre, Faculty of Sciences, Engineering and Health, CQUniversity, Rockhampton, QLD 4702 (Australia); Khan, M M K; Rasul, M G, E-mail: m.rasul@cqu.edu.a [School of Engineering and Built Environment, Faculty of Sciences, Engineering and Health, CQUniversity, Rockhampton, QLD 4702 (Australia)
2010-12-15
This paper presents an experimental and computational study of air bubbles rising in a massecuite-equivalent non-Newtonian crystal suspension. The bubble trajectory inside the stagnant liquid of a 0.05% xanthan gum crystal suspension was investigated and modelled using the computational fluid dynamics (CFD) model to gain an insight into the bubble flow characteristics. The CFD code FLUENT was used for numerical simulation, and the bubble trajectory calculations were performed through a volume of fluid (VOF) model. The influences of the Reynolds number (Re), the Weber number (We) and the bubble aspect ratio (E) on the bubble trajectory are discussed. The conditions for the bubbles' path oscillations are identified. The experimental results showed that the path instability for the crystal suspension was less rapid than in water. The trajectory analysis indicated that 5.76 mm diameter bubbles followed a zigzag motion in the crystal suspension. Conversely, the smaller bubbles (5.76 mm) followed a path of least horizontal movement and larger bubbles (21.21 mm) produced more spiral motion within the crystal suspension. Path instability occurred for bubbles of 15.63 and 21.21 mm diameter, and they induced both zigzag and spiral trajectories within the crystal suspension. At low Re and We, smaller bubbles (5.76 mm) produced a zigzag trajectory, whereas larger bubbles (15.63 and 21.21 mm) showed both zigzag and spiral trajectories at intermediate and moderately high Re and We in the crystal suspension. The simulation results illustrated that a repeating pattern of swirling vortices was created for smaller bubbles due to the unstable wake and unsteady flow of these bubbles. This is the cause of the smaller bubbles moving in a zigzag way. Larger bubbles showed two counter-rotating trailing vortices at the back of the bubble. These vortices induced a velocity component to the gas-liquid interface and caused a deformation. Hence, the larger bubbles produced a path
A modelling and experimental study of the bubble trajectory in a non-Newtonian crystal suspension
This paper presents an experimental and computational study of air bubbles rising in a massecuite-equivalent non-Newtonian crystal suspension. The bubble trajectory inside the stagnant liquid of a 0.05% xanthan gum crystal suspension was investigated and modelled using the computational fluid dynamics (CFD) model to gain an insight into the bubble flow characteristics. The CFD code FLUENT was used for numerical simulation, and the bubble trajectory calculations were performed through a volume of fluid (VOF) model. The influences of the Reynolds number (Re), the Weber number (We) and the bubble aspect ratio (E) on the bubble trajectory are discussed. The conditions for the bubbles' path oscillations are identified. The experimental results showed that the path instability for the crystal suspension was less rapid than in water. The trajectory analysis indicated that 5.76 mm diameter bubbles followed a zigzag motion in the crystal suspension. Conversely, the smaller bubbles (5.76 mm) followed a path of least horizontal movement and larger bubbles (21.21 mm) produced more spiral motion within the crystal suspension. Path instability occurred for bubbles of 15.63 and 21.21 mm diameter, and they induced both zigzag and spiral trajectories within the crystal suspension. At low Re and We, smaller bubbles (5.76 mm) produced a zigzag trajectory, whereas larger bubbles (15.63 and 21.21 mm) showed both zigzag and spiral trajectories at intermediate and moderately high Re and We in the crystal suspension. The simulation results illustrated that a repeating pattern of swirling vortices was created for smaller bubbles due to the unstable wake and unsteady flow of these bubbles. This is the cause of the smaller bubbles moving in a zigzag way. Larger bubbles showed two counter-rotating trailing vortices at the back of the bubble. These vortices induced a velocity component to the gas-liquid interface and caused a deformation. Hence, the larger bubbles produced a path transition.
Development of bubble-induced turbulence model for advanced two-fluid model
diameter pipes. As a result, the analysis results of void fraction distribution were improved by using the suggested model. However, the analytical velocity distribution was flat compared with the experiments. As the main cause of this, we thought that turbulent viscosity in two-phase flow was overestimated in these analyses. To express the velocity distribution, we improved the two-phase k-ε turbulent model in this study. Especially, turbulent dissipation rate transport equation was introduced in order to consider the influence of bubble-induced turbulence. These models were incorporated to the advanced two-fluid model code ACE-3D, and numerical simulations for air-water two-phase flow experiment in 200 mm-, 60 mm- and 38 mm-diameter vertical pipe were performed. From these results, the qualitative phenomena could be expressed and the dependency of the suggested model was confirmed. (author)
Kálal Zbyněk
2014-09-01
Full Text Available The main topic of this study is the mathematical modelling of bubble size distributions in an aerated stirred tank using the population balance method. The air-water system consisted of a fully baffled vessel with a diameter of 0.29 m, which was equipped with a six-bladed Rushton turbine. The secondary phase was introduced through a ring sparger situated under the impeller. Calculations were performed with the CFD software CFX 14.5. The turbulent quantities were predicted using the standard k-ε turbulence model. Coalescence and breakup of bubbles were modelled using the MUSIG method with 24 bubble size groups. For the bubble size distribution modelling, the breakup model by Luo and Svendsen (1996 typically has been used in the past. However, this breakup model was thoroughly reviewed and its practical applicability was questioned. Therefore, three different breakup models by Martínez-Bazán et al. (1999a, b, Lehr et al. (2002 and Alopaeus et al. (2002 were implemented in the CFD solver and applied to the system. The resulting Sauter mean diameters and local bubble size distributions were compared with experimental data.
Highlights: • Flow behaviors with H and Vg were simulated using CFD coupled with PBM. • Bubble size and gas holdup are mainly determined by vortical flow. • Interfacial area enhanced by increasing Vg and H, but the enhancement effect is not obvious while the Vg is too fast. - Abstract: Using bubble column to extract tritium from lead lithium (Pb–17Li) eutectic is an effective way in the process of tritium extraction in liquid blanket system, where the hydrodynamic characteristics of the gas–liquid two-phase flow in the columns play a very important role. In order to understand the two-phase flow details and investigate the influence factors on the hydrodynamic performance, in this paper the flow behaviors in the cylindrical bubble columns with different heights and purge gas inlet velocities using computational fluid dynamics coupled with population balance model were simulated. Liquid flow field, bubble Sauter mean diameter, time-averaged gas holdup and two-phase interfacial area for the different cases were obtained and compared. The simulation results showed good agreement with previous studies, and which indicated that bubble size and gas holdup formation are mainly determined by vortical flow. In addition, interfacial area can be enhanced by increasing the purge gas inlet velocity and column height. However, the enhancement effect will trail away when the gas inlet velocity is too fast, and the contribution of column height is relatively small
Numerical Modeling and Prediction of Bubbling Fluidized Beds
England, Jonas Andrew
2011-01-01
Numerical modeling and prediction techniques are used to determine pressure drop, minimum fluidization velocity and segregation for bubbling fluidized beds. The computational fluid dynamics (CFD) code Multiphase Flow with Interphase eXchange (MFIX) is used to study a two-stage reactor geometry with a binary mixture. MFIX is demonstrated to accurately predict pressure drop versus inlet gas velocity for binary mixtures. A new method is developed to predict the pressure drop versus inlet gas v...
Horizontal bubbly flow with elbow restrictions: Interfacial area transport modeling
The present study develops an interfacial area transport equation applicable to an air-water horizontal bubbly flow, along which two types of horizontal elbows are installed as flow restrictions. Two sets of experiments are performed in a round glass tube of 50.3 mm inner diameter. Along the test section, a 90-degree elbow is installed at L/D = 206.6 from the two-phase mixture inlet and then a 45-degree elbow is installed at L/D = 353.5. In total, 15 different flow conditions in the bubbly flow regime for each of the two flow restriction experiments are studied. Detailed local two-phase flow parameters are acquired by a double-sensor conductivity probe at four different axial locations in the 90-degree experiment and three different axial locations in the 45-degree experiment. The effect of the elbows is found to be evident in the distribution of local parameters as well as in the development of interfacial structures. It is clear that the elbows make an effect on the bubble interactions resulting in significant changes to both the void fraction and interfacial area concentration. In the present analysis, the interfacial area transport equation is developed in one-dimensional form via area averaging. In the averaging process, characteristic non-uniform distributions of the flow parameters in horizontal two-phase flow are treated mathematically through a distribution parameter. The mechanistic models for the major bubble interaction phenomena developed in vertical flow analysis are employed in the present study. Furthermore, the change in pressure due to the minor loss of an elbow is taken into consideration by using a newly developed correlation analogous to Lockhart and Martinelli's. In total, 105 area-averaged data points are employed to benchmark the present model. The present model predicts the data relatively well with an average percent difference of approximately ±20%.
Bubbling phenomenon in a discrete economic model for the interaction of demand and supply
Yang Xiaozhong [School of Mathematics and Statistics, North China Electric Power University, Beijing 102 206 (China)], E-mail: yxiaozh@ncepu.edu.cn; Peng Mingshu [School of Science, Beijing Jiao Tong University, Beijing 100 044 (China)], E-mail: mshpeng@bjtu.edu.cn; Hu Jiping [School of Civil Engineering, Beijing Jiao Tong University, Beijing 100 044 (China)], E-mail: jphu@bjtu.edu.cn; Jiang Xiaoxia [School of Science, Beijing Jiao Tong University, Beijing 100 044 (China)
2009-11-15
In this paper, we study rich dynamics of a nonlinear economic model. Chaotic and bubbling phenomena are shown by numerical simulation, which clearly parallels with phenomena from technology bubbling during 1999-2000, or economic bubbling such as in the global industry of real state/stock market/production (oil) market since 2006.
Liang, Mingchao; Wei, Junhong; Han, Hongmei; Fu, Chengguo; Liu, Jianjun
2015-09-01
The capillary pressure is one of the crucial parameters in many science and engineering applications such as composite materials, interface science, chemical engineering, oil exploration, etc. The drop/bubble formation and its mechanisms that affect the permeability of porous media have steadily attracted much attention in the past. When a drop/bubble moves from a larger capillary to a smaller one, it is often obstructed by an additional pressure difference caused by the capillary force. In this paper, a comprehensive model is derived for the capillary pressure difference when a drop/bubble flows through a constricted capillary, i.e. a geometrically constricted passage with an abrupt change in radius. The proposed model is expressed as a function of the smaller capillary radius, pore-throat ratio, contact angle, surface tension and length of the drop/bubble in the smaller capillary. The model predictions are compared with the available experimental data, and good agreement is found between them.
An application of the method of moments to the modeling of bubbly flow
Bubbly flows are relevant in nuclear reactors thermalhydraulics and safety analysis. Regularly, empirical constitutive laws are required to close the two-fluid equations, particularly in relating the interfacial area and the bubble number densities to the local void fraction. In this article, starting from a generalized Boltzmann transport equation for the bubble size spectrum, a convection equation for the bubble number density is derived using the method of moments. The equation is analyzed for a vertical bubbly flow in stagnated liquid, showing excellent agreement with experimental data. The model is useful as a mean to provide conservation-based correlations to complement the existing two-fluid models
Bubbly flows with fixed polydispersity: Validation of a baseline closure model
Rzehak, Roland, E-mail: R.Rzehak@hzdr.de; Krepper, Eckhard, E-mail: E.Krepper@hzdr.de
2015-06-15
Highlights: • Consideration of regime with significant but non-varying polydispersity. • Facilitates qualification of closure models. • Including bubble forces and bubble-induced turbulence. • Validation of the models with an experimental database for developing flow. • Shows reversing direction of shear lift force with bubble size. - Abstract: For practical applications the Euler–Euler two-fluid model relies on suitable closure relations describing interfacial exchange processes. In dispersed gas–liquid multiphase flow, closures are needed for bubble forces, bubble-induced turbulence, as well as bubble-coalescence and -breakup. The quest for models with a broad range of applicability allowing predictive simulations is an ongoing venture. Reasonable success has been achieved so far for flows that are amenable to a monodisperse approximation for the bubble size which limits the latter to no more than a few mm. In the present work we extend the validation to flow in which bubbles with a broad distribution of sizes up to ∼10 mm are present, but the shape of the distribution remains unchanged during the flow development. The existence of such conditions, which we term “fixed polydispersity”, is deduced from the experimental data. For this kind of situation the complexity of the closure problem is reduced since a balance between bubble-coalescence and -breakup prevails that allows to neglect these processes and simply impose a fixed bubble size distribution. Conclusions towards best practice guidelines for modeling bubbly flows are drawn and needs for further research identified.
Bubbly flows with fixed polydispersity: Validation of a baseline closure model
Highlights: • Consideration of regime with significant but non-varying polydispersity. • Facilitates qualification of closure models. • Including bubble forces and bubble-induced turbulence. • Validation of the models with an experimental database for developing flow. • Shows reversing direction of shear lift force with bubble size. - Abstract: For practical applications the Euler–Euler two-fluid model relies on suitable closure relations describing interfacial exchange processes. In dispersed gas–liquid multiphase flow, closures are needed for bubble forces, bubble-induced turbulence, as well as bubble-coalescence and -breakup. The quest for models with a broad range of applicability allowing predictive simulations is an ongoing venture. Reasonable success has been achieved so far for flows that are amenable to a monodisperse approximation for the bubble size which limits the latter to no more than a few mm. In the present work we extend the validation to flow in which bubbles with a broad distribution of sizes up to ∼10 mm are present, but the shape of the distribution remains unchanged during the flow development. The existence of such conditions, which we term “fixed polydispersity”, is deduced from the experimental data. For this kind of situation the complexity of the closure problem is reduced since a balance between bubble-coalescence and -breakup prevails that allows to neglect these processes and simply impose a fixed bubble size distribution. Conclusions towards best practice guidelines for modeling bubbly flows are drawn and needs for further research identified
Hydrodynamics of Bubble Bouncing on a Wall, Experiment and Modelling
Zedníková, Mária; Vejražka, Jiří; Růžička, Marek; Drahoš, Jiří
Prague: Institute of Hydrodynamics ASCR, v.v.i, 2008 - (Chára, Z.; Klaboch, L.), s. 20-27 ISBN 978-80-87117-04-0. [XXII. Symposium on Anemometry. Holany-Litice (CZ), 03.06.2008-04.06.2008] R&D Projects: GA ČR GA104/07/1110; GA AV ČR(CZ) IAA200720801; GA AV ČR(CZ) KJB200720801 Institutional research plan: CEZ:AV0Z40720504 Keywords : hydrodynamics of bubble * experiments * modelling Subject RIV: CI - Industrial Chemistry, Chemical Engineering
The Accretion Wind Model of Fermi Bubbles. II. Radiation
Mou, Guobin; Yuan, Feng; Gan, Zhaoming; Sun, Mouyuan
2015-09-01
In a previous work, we have shown that the formation of Fermi bubbles can be due to the interaction between winds launched from the hot accretion flow in Sgr A* and the interstellar medium (ISM). In that work, we focus only on the morphology. In this paper we continue our study by calculating the gamma-ray radiation. Some cosmic-ray protons (CRp) and electrons (CRe) must be contained in the winds, which are likely formed by physical processes such as magnetic reconnection. We have performed MHD simulations to study the spatial distribution of CRp, considering the advection and diffusion of CRp in the presence of magnetic field. We find that a permeated zone is formed just outside of the contact discontinuity between winds and the ISM, where the collisions between CRp and thermal nuclei mainly occur. The decay of neutral pions generated in the collisions, combined with the inverse Compton scattering of background soft photons by the secondary leptons generated in the collisions and primary CRe, can well explain the observed gamma-ray spectral energy distribution. Other features such as the uniform surface brightness along the latitude and the boundary width of the bubbles are also explained. The advantage of this “accretion wind” model is that the adopted wind properties come from the detailed small-scale MHD numerical simulation of accretion flows and the value of mass accretion rate has independent observational evidences. The success of the model suggests that we may seriously consider the possibility that cavities and bubbles observed in other contexts such as galaxy clusters may be formed by winds rather than jets.
Modeling high-energy gamma-rays from the Fermi Bubbles
Splettstoesser, Megan
2015-09-17
In 2010, the Fermi Bubbles were discovered at the galactic center of the Milky Way. These giant gamma-ray structures, extending 55° in galactic latitude and 20°-30° in galactic longitude, were not predicted. We wish to develop a model for the gamma-ray emission of the Fermi Bubbles. To do so, we assume that second order Fermi acceleration requires charged particles and irregular magnetic fields- both of which are present in the disk of the Milky Way galaxy. By solving the steady-state case of the transport equation, I compute the proton spectrum due to second order Fermi acceleration. I compare the analytical solutions of the proton spectrum to a numerical solution. I find that the numerical solution to the transport equation converges to the analytical solution in all cases. The gamma-ray spectrum due to proton-proton interaction is compared to Fermi Bubble data (from Ackermann et al. 2014), and I find that second order Fermi acceleration is a good fit for the gamma-ray spectrum of the Fermi Bubbles at low energies with an injection source term of S = 1.5 x 10⁻¹⁰ GeV⁻¹cm⁻³yr⁻¹. I find that a non-steady-state solution to the gamma-ray spectrum with an injection source term of S = 2 x 10⁻¹⁰ GeV⁻¹cm⁻³yr⁻¹ matches the bubble data at high energies.
Modeling Space-Time Dependent Helium Bubble Evolution in Tungsten Armor under IFE Conditions
The High Average Power Laser (HAPL) program is a coordinated effort to develop Laser Inertial Fusion Energy. The implosion of the D-T target produces a spectrum of neutrons, X-rays, and charged particles, which arrive at the first wall (FW) at different times within about 2.5 μs at a frequency of 5 to 10 Hz. Helium is one of several high-energy charged particle constituents impinging on the candidate tungsten armored low activation ferritic steel First Wall. The spread of the implanted debris and burn helium energies results in a unique space-time dependent implantation profile that spans about 10 μm in tungsten. Co-implantation of X-rays and other ions results in spatially dependent damage profiles and rapid space-time dependent temperature spikes and gradients. The rate of helium transport and helium bubble formation will vary significantly throughout the implanted region. Furthermore, helium will also be transported via the migration of helium bubbles and non-equilibrium helium-vacancy clusters. The HEROS code was developed at UCLA to model the spatial and time-dependent helium bubble nucleation, growth, coalescence, and migration under transient damage rates and transient temperature gradients. The HEROS code is based on kinetic rate theory, which includes clustering of helium and vacancies, helium mobility, helium-vacancy cluster stability, cavity nucleation and growth and other microstructural features such as interstitial loop evolution, grain boundaries, and precipitates. The HEROS code is based on space-time discretization of reaction-diffusion type equations to account for migration of mobile species between neighboring bins as single atoms, clusters, or bubbles. HAPL chamber FW implantation conditions are used to model helium bubble evolution in the implanted tungsten. Helium recycling rate predictions are compared with experimental results of helium ion implantation experiments. (author)
Mathematical model of diffusion-limited evolution of multiple gas bubbles in tissue
Srinivasan, R. Srini; Gerth, Wayne A.; Powell, Michael R.
2003-01-01
Models of gas bubble dynamics employed in probabilistic analyses of decompression sickness incidence in man must be theoretically consistent and simple, if they are to yield useful results without requiring excessive computations. They are generally formulated in terms of ordinary differential equations that describe diffusion-limited gas exchange between a gas bubble and the extravascular tissue surrounding it. In our previous model (Ann. Biomed. Eng. 30: 232-246, 2002), we showed that with appropriate representation of sink pressures to account for gas loss or gain due to heterogeneous blood perfusion in the unstirred diffusion region around the bubble, diffusion-limited bubble growth in a tissue of finite volume can be simulated without postulating a boundary layer across which gas flux is discontinuous. However, interactions between two or more bubbles caused by competition for available gas cannot be considered in this model, because the diffusion region has a fixed volume with zero gas flux at its outer boundary. The present work extends the previous model to accommodate interactions among multiple bubbles by allowing the diffusion region volume of each bubble to vary during bubble evolution. For given decompression and tissue volume, bubble growth is sustained only if the bubble number density is below a certain maximum.
Pseudopotential MRT lattice Boltzmann model for cavitation bubble collapse with high density ratio
Shan, Ming-Lei; Yao, Cheng; Yin, Cheng; Jiang, Xiao-Yan
2016-01-01
The dynamics of the cavitation bubble collapse is a fundamental issue for the bubble collapse application and prevention. In present work, the modified forcing scheme for the pseudopotential multi-relaxation-time lattice Boltzmann model developed by Li Q. et al. is adopted to develop a cavitation bubble collapse model. In the respects of coexistence curves and Laplace law verification, the improved pseudopotential multi-relaxation-time lattice Boltzmann model is investigated. The independence between the kinematic viscosity and the thermodynamic consistency, surface tension is founded. By homogeneous and heterogeneous cavitation simulation, the capability of the present model to describe the cavitation bubble development as well as the cavitation inception is verified. The bubble collapse between two parallel walls is simulated. The dynamic process of collapsing bubble is consistent with the results from experiments and simulations by other numerical method. It is demonstrated that the present pseudopotential...
Multiple bifurcations and periodic 'bubbling' in a delay population model
In this paper, the flip bifurcation and periodic doubling bifurcations of a discrete population model without delay influence is firstly studied and the phenomenon of Feigenbaum's cascade of periodic doublings is also observed. Secondly, we explored the Neimark-Sacker bifurcation in the delay population model (two-dimension discrete dynamical systems) and the unique stable closed invariant curve which bifurcates from the nontrivial fixed point. Finally, a computer-assisted study for the delay population model is also delved into. Our computer simulation shows that the introduction of delay effect in a nonlinear difference equation derived from the logistic map leads to much richer dynamic behavior, such as stable node → stable focus → an lower-dimensional closed invariant curve (quasi-periodic solution, limit cycle) or/and stable periodic solutions → chaotic attractor by cascading bubbles (the combination of potential period doubling and reverse period-doubling) and the sudden change between two different attractors, etc
Hu, Shenyang Y.; Burkes, Douglas; Lavender, Curt A.; Senor, David J.; Setyawan, Wahyu; Xu, Zhijie
2016-07-08
Nano-gas bubble superlattices are often observed in irradiated UMo nuclear fuels. However, the for- mation mechanism of gas bubble superlattices is not well understood. A number of physical processes may affect the gas bubble nucleation and growth; hence, the morphology of gas bubble microstructures including size and spatial distributions. In this work, a phase-field model integrating a first-passage Monte Carlo method to investigate the formation mechanism of gas bubble superlattices was devel- oped. Six physical processes are taken into account in the model: 1) heterogeneous generation of gas atoms, vacancies, and interstitials informed from atomistic simulations; 2) one-dimensional (1-D) migration of interstitials; 3) irradiation-induced dissolution of gas atoms; 4) recombination between vacancies and interstitials; 5) elastic interaction; and 6) heterogeneous nucleation of gas bubbles. We found that the elastic interaction doesn’t cause the gas bubble alignment, and fast 1-D migration of interstitials along 〈110〉 directions in the body-centered cubic U matrix causes the gas bubble alignment along 〈110〉 directions. It implies that 1-D interstitial migration along [110] direction should be the primary mechanism of a fcc gas bubble superlattice which is observed in bcc UMo alloys. Simulations also show that fission rates, saturated gas concentration, and elastic interaction all affect the morphology of gas bubble microstructures.
A model of the interaction of bubbles and solid particles under acoustic excitation
Hay, Todd Allen
The Lagrangian formalism utilized by Ilinskii, Hamilton and Zabolotskaya [J. Acoust. Soc. Am. 121, 786-795 (2007)] to derive equations for the radial and translational motion of interacting bubbles is extended here to obtain a model for the dynamics of interacting bubbles and elastic particles. The bubbles and particles are assumed to be spherical but are otherwise free to pulsate and translate. The model is accurate to fifth order in terms of a nondimensional expansion parameter R/d, where R is a characteristic radius and d is a characteristic distance between neighboring bubbles or particles. The bubbles and particles may be of nonuniform size, the particles elastic or rigid, and external acoustic sources are included to an order consistent with the accuracy of the model. Although the liquid is assumed initially to be incompressible, corrections accounting for finite liquid compressibility are developed to first order in the acoustic Mach number for a cluster of bubbles and particles, and to second order in the acoustic Mach number for a single bubble. For a bubble-particle pair consideration is also given to truncation of the model at fifth order in R/d via automated derivation of the model equations to arbitrary order. Numerical simulation results are presented to demonstrate the effects of key parameters such as particle density and size, liquid compressibility, particle elasticity and model order on the dynamics of single bubbles, pairs of bubbles, bubble-particle pairs and clusters of bubbles and particles under both free response conditions and sinusoidal or shock wave excitation.
Pellacani, Filippo
2012-12-04
A local mechanistic model for bubble coalescence and breakup for the one-group interfacial area transport equation has been developed, in agreement and within the limits of the current understanding, based on an exhaustive survey of the theory and of the state of the art models for bubble dynamics simulation. The new model has been tested using the commercial 3D CFD code ANSYS CFX. Upward adiabatic turbulent air-water bubbly flow has been simulated and the results have been compared with the data obtained in the experimental facility PUMA. The range of the experimental data available spans between 0.5 to 2 m/s liquid velocity and 5 to 15 % volume fraction. For the implementation of the models, both the monodispersed and the interfacial area transport equation approaches have been used. The first one to perform a detailed analysis of the forces and models to reproduce the dynamic of the dispersed phase adequately and to be used in the next phases of the work. Also two different bubble induced turbulence models have been tested to consider the effect of the presence of the gas phase on the turbulence of the liquid phase. The interfacial area transport equation has been successfully implemented into the CFD code and the state of the art breakup and coalescence models have been used for simulation. The limitations of the actual theory have been shown and a new bubble interactions model has been developed. The simulations showed that a considerable improvement is achieved if compared to the state of the art closure models. Limits in the implementation derive from the actual understanding and formulation of the bubbly dynamics. A strong dependency on the interfacial non-drag force models and coefficients have been shown. More experimental and theory work needs to be done in this field to increase the prediction capability of the simulation tools regarding the distribution of the phases along the pipe radius.
A local mechanistic model for bubble coalescence and breakup for the one-group interfacial area transport equation has been developed, in agreement and within the limits of the current understanding, based on an exhaustive survey of the theory and of the state of the art models for bubble dynamics simulation. The new model has been tested using the commercial 3D CFD code ANSYS CFX. Upward adiabatic turbulent air-water bubbly flow has been simulated and the results have been compared with the data obtained in the experimental facility PUMA. The range of the experimental data available spans between 0.5 to 2 m/s liquid velocity and 5 to 15 % volume fraction. For the implementation of the models, both the monodispersed and the interfacial area transport equation approaches have been used. The first one to perform a detailed analysis of the forces and models to reproduce the dynamic of the dispersed phase adequately and to be used in the next phases of the work. Also two different bubble induced turbulence models have been tested to consider the effect of the presence of the gas phase on the turbulence of the liquid phase. The interfacial area transport equation has been successfully implemented into the CFD code and the state of the art breakup and coalescence models have been used for simulation. The limitations of the actual theory have been shown and a new bubble interactions model has been developed. The simulations showed that a considerable improvement is achieved if compared to the state of the art closure models. Limits in the implementation derive from the actual understanding and formulation of the bubbly dynamics. A strong dependency on the interfacial non-drag force models and coefficients have been shown. More experimental and theory work needs to be done in this field to increase the prediction capability of the simulation tools regarding the distribution of the phases along the pipe radius.
A mathematical model and numerical simulation of pressure wave in horizontal gas-liquid bubbly flow
HUANG Fei; BAI Bofeng; GUO Liejin
2004-01-01
By using an ensemble-averaged two-fluid model,with valid closure conditions of interfacial momentum exchange due to virtual mass force,viscous shear stress and drag force,a model for pressure wave propagation in a horizontal gas-liquid bubbly flow is proposed.According to the small perturbation theory and solvable condition of one-order linear uniform equations,a dispersion equation of pressure wave is induced.The pressure wave speed calculated from the model is compared and in good agreement with existing data.According to the dispersion equation,the propagation and attenuation of pressure wave are investigated systemically.The factors affecting pressure wave,such as void fraction,pressure,wall shear stress,perturbation frequency,virtual mass force and drag force,are analyzed.The result shows that the decrease in system pressure,the increase in void fraction and the existence of wall shear stress,will cause a decrease in pressure wave speed and an increase in the attenuation coefficient in the horizontal gas-liquid bubbly flow.The effects of perturbation frequency,virtual mass and drag force on pressure wave in the horizontal gas-liquid bubbly flow at low perturbation frequency are different from that at high perturbation frequency.
Bubble Economics and Structural Change: The Cases of Spain and France Compared
Agnese, Pablo; Hromcová, Jana
2015-01-01
This paper delves into the recent events that led to the formation of the housing bubble in Spain and the resulting structural change that is arguably needed to put the economy back into the right track. For this purpose we calibrate a model with different equilibria descriptive of the labor markets in Spain and France, where the unemployment rates went from the same initial spot to very different levels. In addition to this, we run two counterfactual analyses that throw some more light on th...
COMPUTATIONAL AND EXPERIMENTAL MODELING OF SLURRY BUBBLE COLUMN REACTORS
Paul C.K. Lam; Isaac K. Gamwo; Dimitri Gidaspow
2002-05-01
The objective of this study was to develop a predictive experimentally verified computational fluid dynamics (CFD) model for gas-liquid-solid flow. A three dimensional transient computer code for the coupled Navier-Stokes equations for each phase was developed and is appended in this report. The principal input into the model is the viscosity of the particulate phase which was determined from a measurement of the random kinetic energy of the 800 micron glass beads and a Brookfield viscometer. The details are presented in the attached paper titled ''CFD Simulation of Flow and Turbulence in a Slurry Bubble Column''. This phase of the work is in press in a referred journal (AIChE Journal, 2002) and was presented at the Fourth International Conference on Multiphase Flow (ICMF 2001) in New Orleans, May 27-June 1, 2001 (Paper No. 909). The computed time averaged particle velocities and concentrations agree with Particle Image Velocimetry (PIV) measurements of velocities and concentrations, obtained using a combination of gamma-ray and X-ray densitometers, in a slurry bubble column, operated in the bubbly-coalesced fluidization regime with continuous flow of water. Both the experiment and the simulation show a down-flow of particles in the center of the column and up-flow near the walls and nearly uniform particle concentration. Normal and shear Reynolds stresses were constructed from the computed instantaneous particle velocities. The PIV measurement and the simulation produced instantaneous particle velocities. The PIV measurement and the simulation produced similar nearly flat horizontal profiles of turbulent kinetic energy of particles. To better understand turbulence we studied fluidization in a liquid-solid bed. This work was also presented at the Fourth International Conference on Multiphase Flow (ICMF 2001, Paper No. 910). To understand turbulence in risers, measurements were done in the IIT riser with 530 micron glass beads using a PIV
An Axial Dispersion Model for Evaporating Bubble Column Reactor
谢刚; 李希
2004-01-01
Evaporating bubble column reactor (EBCR) is a kind of aerated reactor in which the reaction heat is removed by the evaporation of volatile reaction mixture. In this paper, a mathematical model that accounts for the gas-liquid exothermic reaction and axial dispersions of both gas and liquid phase is employed to study the performance of EBCR for the process of p-xylene(PX) oxidation. The computational results show that there are remarkable concentration and temperature gradients in EBCR for high ratio of height to diameter (H/DT). The temperature is lower at the bottom of column and higher at the top, due to rapid evaporation induced by the feed gas near the bottom. The concentration profiles in the gas phase are more nonuniform than those (except PX) in the liquid phase, which causes more solvent burning consumption at high H/DT ratio. For p-xylene oxidation, theo ptimal H/DT is around 5.
Cheung, S.C.P.; Deju, L. [School of Aerospace, Mechanical and Manufacturing Engineering, RMIT University, Victoria 3083 (Australia); Yeoh, G.H. [Australian Nuclear Science and Technology Organisation, Locked Bag 2001, Kirrawee DC, New South Wales 2232 (Australia); School of Mechanical and Manufacturing Engineering, University of New South Wales, New South Wales 2052 (Australia); Tu, J.Y., E-mail: jiyuan.tu@rmit.edu.au [School of Aerospace, Mechanical and Manufacturing Engineering, RMIT University, Victoria 3083 (Australia)
2013-12-15
Highlights: • Analyze gas–liquid flow in medium and large pipe under various flow regimes. • Identifying relative merits and capabilities of DQMOM. • Comparison of DQMOM with ABND and homogeneous MUSIG model. • Numerical results validation against gas–liquid flow experiments. - Abstract: Gas–liquid flows are commonly encountered in many industrial flow systems. In many cases, the evolution of bubble size distribution is a crucial factor governing the momentum, heat and mass transfer between phases within the system. Aiming to evaluate the capability of existing models, numerical assessment of three different population balance approaches – direct quadrature method of moments (DQMOM), average bubble number density (ABND) model and homogeneous MUlti-SIze-Group (MUSIG) models – is presented in this paper. Model predictions were validated against experimental measurements from medium and large scale bubble columns where bubble sizes within the system were found to be dominant by coalescence and break-up mechanism, respectively. In result of the bubble size change, transitions of phase distribution from wall peak to core peak profile were also found in both experiments. In general, predictions of the three models were in satisfactory agreement with experiment measurements clearly demonstrating its applicability for large scale practical systems. Encouraging results have also been obtained in capturing the evolution of bubble size distribution. Nevertheless, noticeable errors were also found in predictions of the MUSIG and DQMOM model indicating some potential deficiencies of the model. To evaluate the numerical efficiency of the three models, computational requirements of each model were also compared.
Highlights: • Analyze gas–liquid flow in medium and large pipe under various flow regimes. • Identifying relative merits and capabilities of DQMOM. • Comparison of DQMOM with ABND and homogeneous MUSIG model. • Numerical results validation against gas–liquid flow experiments. - Abstract: Gas–liquid flows are commonly encountered in many industrial flow systems. In many cases, the evolution of bubble size distribution is a crucial factor governing the momentum, heat and mass transfer between phases within the system. Aiming to evaluate the capability of existing models, numerical assessment of three different population balance approaches – direct quadrature method of moments (DQMOM), average bubble number density (ABND) model and homogeneous MUlti-SIze-Group (MUSIG) models – is presented in this paper. Model predictions were validated against experimental measurements from medium and large scale bubble columns where bubble sizes within the system were found to be dominant by coalescence and break-up mechanism, respectively. In result of the bubble size change, transitions of phase distribution from wall peak to core peak profile were also found in both experiments. In general, predictions of the three models were in satisfactory agreement with experiment measurements clearly demonstrating its applicability for large scale practical systems. Encouraging results have also been obtained in capturing the evolution of bubble size distribution. Nevertheless, noticeable errors were also found in predictions of the MUSIG and DQMOM model indicating some potential deficiencies of the model. To evaluate the numerical efficiency of the three models, computational requirements of each model were also compared
COMPUTATIONAL AND EXPERIMENTAL MODELING OF SLURRY BUBBLE COLUMN REACTORS
Paul Lam; Dimitri Gidaspow
2000-09-01
The objective if this study was to develop a predictive experimentally verified computational fluid dynamics (CFD) model for gas-liquid-solid flow. A three dimensional transient computer code for the coupled Navier-Stokes equations for each phase was developed. The principal input into the model is the viscosity of the particulate phase which was determined from a measurement of the random kinetic energy of the 800 micron glass beads and a Brookfield viscometer. The computed time averaged particle velocities and concentrations agree with PIV measurements of velocities and concentrations, obtained using a combination of gamma-ray and X-ray densitometers, in a slurry bubble column, operated in the bubbly-coalesced fluidization regime with continuous flow of water. Both the experiment and the simulation show a down-flow of particles in the center of the column and up-flow near the walls and nearly uniform particle concentration. Normal and shear Reynolds stresses were constructed from the computed instantaneous particle velocities. The PIV measurement and the simulation produced instantaneous particle velocities. The PIV measurement and the simulation produced similar nearly flat horizontal profiles of turbulent kinetic energy of particles. This phase of the work was presented at the Chemical Reaction Engineering VIII: Computational Fluid Dynamics, August 6-11, 2000 in Quebec City, Canada. To understand turbulence in risers, measurements were done in the IIT riser with 530 micron glass beads using a PIV technique. The results together with simulations will be presented at the annual meeting of AIChE in November 2000.
GISAXS modelling of helium-induced nano-bubble formation in tungsten and comparison with TEM
Thompson, Matt; Sakamoto, Ryuichi; Bernard, Elodie; Kirby, Nigel; Kluth, Patrick; Riley, Daniel; Corr, Cormac
2016-05-01
Grazing-incidence small angle x-ray scattering (GISAXS) is a powerful non-destructive technique for the measurement of nano-bubble formation in tungsten under helium plasma exposure. Here, we present a comparative study between transmission electron microscopy (TEM) and GISAXS measurements of nano-bubble formation in tungsten exposed to helium plasma in the Large Helical Device (LHD) fusion experiment. Both techniques are in excellent agreement, suggesting that nano-bubbles range from spheroidal to ellipsoidal, displaying exponential diameter distributions with mean diameters μ=0.68 ± 0.04 nm and μ=0.6 ± 0.1 nm measured by TEM and GISAXS respectively. Depth distributions were also computed, with calculated exponential depth distributions with mean depths of 8.4 ± 0.5 nm and 9.1 ± 0.4 nm for TEM and GISAXS. In GISAXS modelling, spheroidal particles were fitted with an aspect ratio ε=0.7 ± 0.1. The GISAXS model used is described in detail.
Development of Bubble Lift-off Diameter Model for Subcooled Boiling Flows
A lot of models and correlations for predicting the bubble departure/lift-off diameter are available in the literature. Most of them were developed based on a hydrodynamic principle, which balances forces acting on a bubble at the departure/lift-off point. One difficulty of these models is lack of essential information, such as bubble front velocity, liquid velocity, or relative velocity, to estimate the active force elements. Hence, the lift-off bubble diameter predicted by these hydrodynamic-controlled models may be suffered a large uncertainty. In contract to the hydrodynamic approach, there are few models developed based on the heat transfer aspect. By balancing the heat conducted through a microlayer underneath a bubble with the heat taken away by condensation at the upper part of the bubble, Unal derived a heat-controlled model of the bubble lift-off diameter. This model did not consider the role of superheat liquid layer surrounding the bubble as well as the effect of liquid properties on the heat transfer process. Beside these two approaches, several empirical correlations have been proposed based on dimensionless analyses for measured experimental databases. The application of these correlations to different experiments conditions is, of course, questionable because of the lack of physical bases. Regarding the heat transfer accompanied by a vapor bubble, four involved heat transfer regions surrounding this bubble can be defined as in Fig. 1. These are dry region, microlayer, superheated liquid layer (SpLL) and subcooled liquid layer (SbLL). The existing of the microlayer is confirmed by experiments, and it is considered to be very effective in the heat transfer. Sernas and Hoper defined five types of the microlayer and indicated that the microlayer acting as a very thick liquid layer gives a best prediction for the bubble growth. However, beside the microlayer, the SpLL might play an important role in the heat transfer if its effective heat transfer area
Modeling and Measurements of Multiphase Flow and Bubble Entrapment in Steel Continuous Casting
Jin, Kai; Thomas, Brian G.; Ruan, Xiaoming
2016-02-01
In steel continuous casting, argon gas is usually injected to prevent clogging, but the bubbles also affect the flow pattern, and may become entrapped to form defects in the final product. To investigate this behavior, plant measurements were conducted, and a computational model was applied to simulate turbulent flow of the molten steel and the transport and capture of argon gas bubbles into the solidifying shell in a continuous slab caster. First, the flow field was solved with an Eulerian k- ɛ model of the steel, which was two-way coupled with a Lagrangian model of the large bubbles using a discrete random walk method to simulate their turbulent dispersion. The flow predicted on the top surface agreed well with nailboard measurements and indicated strong cross flow caused by biased flow of Ar gas due to the slide-gate orientation. Then, the trajectories and capture of over two million bubbles (25 μm to 5 mm diameter range) were simulated using two different capture criteria (simple and advanced). Results with the advanced capture criterion agreed well with measurements of the number, locations, and sizes of captured bubbles, especially for larger bubbles. The relative capture fraction of 0.3 pct was close to the measured 0.4 pct for 1 mm bubbles and occurred mainly near the top surface. About 85 pct of smaller bubbles were captured, mostly deeper down in the caster. Due to the biased flow, more bubbles were captured on the inner radius, especially near the nozzle. On the outer radius, more bubbles were captured near to narrow face. The model presented here is an efficient tool to study the capture of bubbles and inclusion particles in solidification processes.
To qualify CFD codes for two-phase flows, they have to be equipped with constitutive laws describing the interaction between the gaseous and the liquid phases. In the case of bubble flow this particularly concerns the forces acting on the bubbles and bubble coalescence and break-up. To obtain detailed experimental data, an electrode wire-mesh sensor was used, which enables the measurement of the phase distribution with a very high resolution in space and in time. Air-water flow at ambient conditions in a vertical pipe (51.2 mm inner diameter) is investigated to have well defined boundary conditions. Local bubble size distributions are calculated from the data. The measurements were done in different distances from the gas injection device. As a result the development of bubble size distributions as well as the development of the radial gas fraction profiles can be studied. It was found, that the bubble size distribution as well as local effects determine the transition from bubble flow to slug flow. The data are used for the development of a model, which predicts the development of the bubble size distribution and the transition from bubble flow to slug flow in case of stationary flow in a vertical pipe. (orig.)
Mathematical models for tritium permeation analysis in liquid metal flows with helium bubbles
Highlights: • Tritium transport in wall-attached He bubbles is studied. • Tritium mainly by-passes the bubble and its partial pressure in the bubble is homogeneous. • The bubble boundary layer depth is about twice the bubble radius. • Simplified tritium transport models are developed and validated. • A preliminary tritium wall function is proposed for wall-attached He bubbles flows. - Abstract: In LIBRETTO-2 test, evidence was obtained that helium bubbles nucleated and grew in the neutron irradiated PbLi probes. If such phenomenon occurs inside liquid metal (LM) breeding blanket channels, the study of its effect on tritium permeation and heat transfer in the near wall region will acquire utmost importance. The T4F research group has developed in the past a nucleation, growth and transport model for helium bubbles in LM flows, as well as a tritium transport model in such a multi-fluid system. In the present study, we are focused on the near-wall region analysis in order to obtain a wall function that allow reproducing the tritium permeation with coarse meshes and, hence, reduce the computational time. First, we perform some detailed CFD simulations of the near-wall region where bubbles might be attached. In these simulations, tritium diffusion processes as well as tritium recombination and dissociation are modelled. The analysis of such simulations allows us to further understand the complex phenomena and justify the use of simplified models. As a result, a new model for tritium transport across a LM–solid interface partially covered by helium bubbles is developed, implemented and validated. This simplified model can be seen as a wall function for the CFD simulation which substantially reduces computational time
Modeling High-Energy Gamma-Rays from the Fermi Bubbles - Oral Presentation
Splettstoesser, Megan [SLAC National Accelerator Lab., Menlo Park, CA (United States)
2015-08-25
In 2010, the Fermi Bubbles were discovered at the galactic center of the Milky Way. These giant gamma-ray structures, extending 50 degrees in galactic latitude and 20-30 degrees in galactic longitude, were not predicted. We wish to develop a model for the gamma-ray emission of the Fermi Bubbles. To do so, we assume that second order Fermi acceleration is responsible for the high-energy emission of the bubbles. Second order Fermi acceleration requires charged particles and irregular magnetic fields—both of which are present in the disk of the Milky Way galaxy. I use the assumption of second order Fermi acceleration in the transport equation, which describes the diffusion of particles. By solving the steady-state case of the transport equation, I compute the proton spectrum due to Fermi second order acceleration and compare this analytical solution to a numerical solution provided by Dr. P. Mertsch. Analytical solutions to the transport equation are taken from Becker, Le, & Dermer and are used to further test the numerical solution. I find that the numerical solution converges to the analytical solution in all cases. Thus, we know the numerical solution accurately calculates the proton spectrum. The gamma-ray spectrum follows the proton spectrum, and will be computed in the future.
Numerical and physical modelling of bubbly flow phenomena
Sangani, A. S.
1991-01-01
The objective of this study is to develop theoretical tools -- analytical as well as numerical -- for understanding how the flows of bubbly liquids are affected by its microstructure, i.e., the detailed spatial, size, and velocity distribution of bubbles, and how the microstructure, in turn, is affected by the flow. This report describes the progress made to date on the several problems that are being studied. The first problem is concerned with the molecular-dynamics type simulations of monodispersed bubbly liquids under equilibrium and homogeneous conditions and their application to slightly inhomogeneous flows. The Reynolds number is large and the Weber and Froud numbers are small in these simulations. The second problem is concerned with the simulations of flows of bubbly liquids undergoing small amplitude oscillatory motion. Both the cases of bubbles with rigid (due to impurities) and stress-free interfaces are examined. The results are related to the added mass, Basset, and viscous drag coefficients. The third problem is concerned with the acoustic wave propagation in bubbly liquids at frequencies above natural frequency of the bubbles. The second problem is completed as of this writing. Work on the other two problems is currently in progress. A summary of the work to be carried out during the period 1/91 to 6/92 is given in the last section.
Micro-bubble drag reduction on a high speed vessel model
Yanuar; Gunawan; Sunaryo; Jamaluddin, A.
2012-09-01
Ship hull form of the underwater area strongly influences the resistance of the ship. The major factor in ship resistance is skin friction resistance. Bulbous bows, polymer paint, water repellent paint (highly water-repellent wall), air injection, and specific roughness have been used by researchers as an attempt to obtain the resistance reduction and operation efficiency of ships. Micro-bubble injection is a promising technique for lowering frictional resistance. The injected air bubbles are supposed to somehow modify the energy inside the turbulent boundary layer and thereby lower the skin friction. The purpose of this study was to identify the effect of injected micro bubbles on a navy fast patrol boat (FPB) 57 m type model with the following main dimensions: L=2 450 mm, B=400 mm, and T=190 mm. The influence of the location of micro bubble injection and bubble velocity was also investigated. The ship model was pulled by an electric motor whose speed could be varied and adjusted. The ship model resistance was precisely measured by a load cell transducer. Comparison of ship resistance with and without micro-bubble injection was shown on a graph as a function of the drag coefficient and Froude number. It was shown that micro bubble injection behind the mid-ship is the best location to achieve the most effective drag reduction, and the drag reduction caused by the micro-bubbles can reach 6%-9%.
Corrigan, Jackie
2004-01-01
, a computational model developed at Glenn, that simulates the cavitational collapse of a single bubble in a liquid (water) and the subsequent combustion of the gaseous contents inside the bubble. The model solves the time-dependent, compressible Navier-Stokes equations in one-dimension with finite-rate chemical kinetics using the CHEMKIN package. Specifically, parameters such as frequency, pressure, bubble radius, and the equivalence ratio were varied while examining their effect on the maximum temperature, radius, and chemical species. These studies indicate that the radius of the bubble is perhaps the most critical parameter governing bubble combustion dynamics and its efficiency. Based on the results of the parametric studies, we plan on conducting experiments to study the effect of ultrasonic perturbations on the bubble generation process with respect to the bubble radius and size distribution.
To predict the behavior of air-water two-phase flows in a centrifugal pump impeller, a three-dimensional numerical method is proposed based on a bubbly flow model. If it is assumed that the mixtures are homogeneous bubbly flow containing fine bubbles compared with the characteristic length of the impeller channel, then the equations of motion of the mixtures are represented by those of liquid phase and its velocity is expressed as a potential for the quasi-harmonic equation. The equations are solved by use of the finite element method to obtain the velocities and pressures, and the equation of motion of an air bubble is integrated numerically on this flow field to obtain the void fraction. These calculations are repeated until the solutions converge. The results obtained show good agreement with experiments within the range of bubbly flow regime. (author)
A model of phase distribution in bubble flow
The knowledge of phase distribution in a channel with bubble flow is important for several problems in nuclear reactor technology. It is possible to assume the random bubble motion as a Markov diffusion process caused by turbulent liquid velocity fluctuations. A diffusion coefficient can be evaluated in this case. Diffusion equations and their boundary conditions are given to calculate bubble distribution in some cases of practical interest. The solution is possible for the most problems on numerical way only. Analytical solutions are shown in the case of a small bubble source near the wall, for instance a water to sodium microleak in a sodium-heated steam generator. Results of a two-dimensional calculation are discussed from the point of view of microleak detection. (author)
Modelling of bubble trajectories in a pump impeller
Dupoiron, Marine; Linden, Paul
2015-11-01
A vertical rotating flow in an annulus gap with an increasing diameter is used to approximate the flow in a pump impeller. We study a spherical gas bubble released at the flow inlet, subject to turbulent drag and added mass forces. Bubbles trajectories have been computed for different geometries, rotation speeds and bubble size, showing a deviation from the liquid streamlines in the angular and radial directions. This effect is related to the pump performance in multiphase conditions: the velocity difference between the gas and the liquid phases changes the final pressure rise produced by the impeller. In some extreme cases, the centrifugal force can be large enough to prevent bubbles from exiting the impeller at all, leading to an unwanted gas accumulation and the blockage of the pump. We eventually quantify the effects of geometrical and operational parameters on the pump behaviour. Work done in collaboration with Schlumberger Gould Research, Cambridge.
Comment on "Simple improvements to classical bubble nucleation models"
Schmelzer, Jürn W. P.; Baidakov, Vladimir G.
2016-08-01
A critical analysis of several statements concerning experimental studies, molecular dynamics simulations, and the theoretical interpretation of bubble nucleation processes is performed. In particular, it is shown that the Tolman equation does not supply us, in general, with a satisfactory theoretically founded description of the curvature dependence of the surface tension and the dependence of the steady-state nucleation rate of bubbles and droplets on supersaturation in the framework of classical nucleation theory.
Inhomogeneous MUSIG Model - a population balance approach for polydispersed bubbly flows
Many flow regimes in Nuclear Reactor Safety (NRS) Research are characterized by multiphase flows, with one phase being a continuous liquid and the other phase consisting of gas or vapour of the liquid phase. In the range of low to intermediate volume fraction of the gaseous phase the multiphase flow under consideration is a bubbly or slug flow, where the disperse phase is characterized by an evolving bubble size distribution due to bubble breakup and coalescence processes. The paper presents a generalized inhomogeneous Multiple Size Group (MUSIG) Model. Within this model the disperse gaseous phase is divided into N inhomogeneous velocity groups (phases) and each of these groups is subdivided into M bubble size classes. Bubble breakup and coalescence processes between all bubble size classes are taken into account by appropriate models. The derived inhomogeneous MUSIG model has been validated against experimental data from the TOPFLOW test facility at the Research Center Rossendorf (FZR). Comparisons of gas volume fraction and velocity profiles with TOPFLOW-074 test case data are provided, showing the applicability and accuracy of the model for polydispersed bubbly flow in large diameter vertical pipe flow. (author)
Time and Space Dependent Stochastic Acceleration Model for the Fermi Bubbles
Sasaki, K; Terasawa, T
2015-01-01
Fermi-LAT reveals two huge gamma-ray bubbles existing in the Galactic Center, called 'Fermi Bubbles'. The existence of two microwave bubbles at the same region are also reported by the observation by WMAP, dubbed 'WMAP haze'. In order to explain these components, It has been argued that the gamma-rays arise from Inverse-Compton scattering of relativistic electrons accelerated by plasma turbulence, and the microwaves are radiated by synchrotron radiation. But no previous research reproduces both the Fermi Bubbles and WMAP haze under typical magnetic fields in the galaxy. We assume that shocks present in the bubbles and the efficiency of the acceleration by plasma turbulence, 'stochastic acceleration', changes with the distance from the shock front. The distance from the shock front increases with time, accordingly the efficiency of the acceleration changes with time. We also consider the time development of the electrons escape from the turbulence by diffusive loss. Our model succeed to reproduce both the obse...
A heat transfer model for evaporating micro-channel coalescing bubble flow
Consolini, L.; Thome, J.R. [Ecole Polytechnique Federale de Lausanne (Switzerland). Lab. de Transfert de Chaleur et de Masse], e-mail: lorenzo.consolini@epfl.ch, E-mail: john.thome@epfl.ch
2009-07-01
The current study presents a one-dimensional model of confined coalescing bubble flow for the prediction of micro-channel convective boiling heat transfer. Coalescing bubble flow has recently been identified as one of the characteristic flow patterns to be found in micro-scale systems, occurring at intermediate vapor qualities between the isolated bubble and the fully annular regimes. As two or more bubbles bond under the action of inertia and surface tension, the passage frequency of the bubble liquid slug pair declines, with a redistribution of liquid among the remaining flow structures. Assuming heat transfer to occur only by conduction through the thin evaporating liquid film surrounding individual bubbles, the present model includes a simplified description of the dynamics of the thin film evaporation process that takes into account the added mass transfer by breakup of the bridging liquid slugs. The new model has been confronted against experimental data taken within the coalescing bubble flow mode that have been identified by a diabatic micro-scale flow pattern map. The comparisons for three different fluids (R-134a, R-236fa and R-245fa) gave encouraging results with 83% of the database predicted within a {+-} 30% error band. (author)
A new condensation model for the FLUENT 6 code was developed. This model is based on a chemical species mixing model. In the areas where condensation takes place, the model modifies the terms in the transfer equations. The model allows condensation in the bulk and condensation on walls to be calculated. Condensation during steam passage over the area with water is modelled in a simplified manner by means of a porous body. Condensation during steam passage through the lattice can be modelled in a simplified manner by means of the 'porous jump' area. The condensation model was tested on a model of the EREC facility and on a model of the VVER-440/213 vacuum bubble condenser in LB LOCA circumstances. The calculated pressure levels in the EREC facility were compared with the observed data, the pressure levels in the bubble condensation tower was compared with COCOSYS data. The fit was reasonable in both cases. The results point to and unbalanced load of the caps in the EREC facility: the condenser right side is loaded more than the left side. In the bubble condenser tower, the lower plates are stressed more than the upper plates
A Revised Nuclear Star Formation Driven, Hadronic Model for the Fermi Bubbles
Crocker, Roland M; Carretti, Ettore; Hill, Alex S; Sutherland, Ralph S
2013-01-01
The Fermi Bubbles are enigmatic \\gamma-ray features of the Galactic bulge. Both putative activity (within $\\sim$ few $\\times$ Myr) connected to the Galactic center super-massive black hole and, alternatively, nuclear star formation have been claimed as the energising source of the Bubbles. Likewise, both inverse-Compton emission by non-thermal electrons (`leptonic' models) and collisions between non-thermal protons and gas (`hadronic' models) have been advanced as the process supplying the Bubbles' \\gamma -ray emission. An issue for any steady state hadronic model is that the very low density of the Bubbles' plasma seems to require that they accumulate protons over a multi-Gyr timescale, much longer than other natural timescales occurring in the problem. Here we present a hadronic model where the timescale for generating the Bubbles' hadronic \\gamma -ray emission is $\\sim$ few $\\times 10^8$ years. Our model invokes collapse of the Bubbles' thermally-unstable plasma, leading to an accumulation of cosmic rays a...
Widera, Paweł
2011-01-01
The process of comparison of computer generated protein structural models is an important element of protein structure prediction. It has many uses including model quality evaluation, selection of the final models from a large set of candidates or optimisation of parameters of energy functions used in template free modelling and refinement. Although many protein comparison methods are available online on numerous web servers, their ability to handle a large scale model comparison is often very limited. Most of the servers offer only a single pairwise structural comparison, and they usually do not provide a model-specific comparison with a fixed alignment between the models. To bridge the gap between the protein and model structure comparison we have developed the Protein Models Comparator (pm-cmp). To be able to deliver the scalability on demand and handle large comparison experiments the pm-cmp was implemented "in the cloud". Protein Models Comparator is a scalable web application for a fast distributed comp...
Modeling the impediment of methane ebullition bubbles by seasonal lake ice
S. Greene
2014-07-01
Full Text Available Microbial methane (CH4 ebullition (bubbling from anoxic lake sediments comprises a globally significant flux to the atmosphere, but ebullition bubbles in temperate and polar lakes can be trapped by winter ice cover and later released during spring thaw. This "ice-bubble storage" (IBS constitutes a novel mode of CH4 emission. Before bubbles are encapsulated by downward-growing ice, some of their CH4 dissolves into the lake water, where it may be subject to oxidation. We present field characterization and a model of the annual CH4 cycle in Goldstream Lake, a thermokarst (thaw lake in interior Alaska. We find that summertime ebullition dominates annual CH4 emissions to the atmosphere. Eighty percent of CH4 in bubbles trapped by ice dissolves into the lake water column in winter, and about half of that is oxidized. The ice growth rate and the magnitude of the CH4 ebullition flux are important controlling factors of bubble dissolution. Seven percent of annual ebullition CH4 is trapped as IBS and later emitted as ice melts. In a future warmer climate, there will likely be less seasonal ice cover, less IBS, less CH4 dissolution from trapped bubbles, and greater CH4 emissions from northern lakes.
Measuring and modeling the bubble population produced by an underwater explosion.
Holt, Fred D; Lee Culver, R
2011-11-01
Underwater explosions have been studied intensively in the United States since 1941 [e.g., R. H. Cole, Underwater Explosions (Princeton University Press, Princeton, NJ, 1945), pp. 3-13]. Research to date has primarily focused on the initial shock and subsequent pressure waves caused by the oscillations of the "gas-globe" resulting from charge detonation. These phenomena have relatively short timescales (typically less than 2 s). However, after the gas-globe rises through the water column and breaks the surface, there remains behind a cloud of bubbles and perhaps debris from the explosion container which has been markedly less studied. A recent experiment measured the spatial and temporal acoustic response of the bubble cloud resulting from a 13.6 kg PBXN-111 charge detonated at 15.2 m (50 ft) depth. A directional projector was used to propagate linear frequency-modulated (5-65 kHz) and 40 kHz tonal pulses through the bubble cloud. Two hydrophone arrays were positioned so as to measure the energy lost in propagating through the bubble cloud. Three methods have been utilized to invert measurements and estimate the bubble population. The bubble population estimates have been used to develop a model for the bubble population resulting from an underwater explosion. PMID:22088003
A multi-size model for sub-cooled boiling bubbly flows
In this paper we present a model which takes into account the multi-size aspect of a bubble population in a two-phase flow. Several methods exist in the literature to deal with such problematic. For its relative simplicity, the method of moments is here used in addition to a classical two-fluid model. It consists in determining the evolution of the topological and statistical quantities of the two-phase flow by solving additional transport equations on the bubble size distribution moment densities. These so-called geometrical equations involve all phenomena responsible for the changes in the bubble size distribution, namely coalescence and break-up but also dispersed phase compressibility and phase change. To close the bubble size distribution appearing in different source terms of the geometrical equations, a two-parameter quadratic expression is here proposed. The simplicity of this quadratic law allows to derived quite easily the geometrical equation source terms and implies to solve only two additional equations. The resulting multi-size model is therefore easy to implement and numerically efficient. Numerical simulation of a DEBORA test case (upward bubbly boiling flow in a vertical tube) has been performed with the NEPTUNE-CFD code to test the ability of the multi-size model to represent sub-cooled boiling bubbly flows. (authors)
An equation of motion for bubble growth
A mathematical model is developed which describes asymmetric bubble growth, either during boiling or bubble injection from submerged orifices. The model is developed using the integral form of the continuity and momentum equations, resulting in a general expression for the acceleration of the bubble's centre of gravity. The proposed model highlights the need to include acceleration due to an asymmetric gain or loss of mass in order to accurately predict bubble motion. Some scenarios are posed by which the growth of bubbles, particularly idealized bubbles that remain a section of a sphere, must include the fact that bubble growth can be asymmetric. In particular, for approximately hemispherical bubble growth the sum of the forces acting on the bubble is negligible compared with the asymmetric term. Further, for bubble injection from a submerged needle this component in the equation of motion is very significant during the initial rapid growth phase as the bubble issues from the nozzle changing from a near hemisphere to truncated sphere geometry. (author)
Prediction of bubble departure in forced convection boiling: a mechanistic model
Colombo, M; Fairweather, M.
2015-01-01
In the context of computational fluid dynamic simulations of boiling flows using time-averaged Eulerian multi-phase approaches, the many sub-models required to describe such a complex phenomena are of particular importance. Of interest here, wall boiling requires calculation of the contribution of evaporation to global heat transfer, which in turn relies on determination of the active nucleation site density, bubble departure diameter and frequency of bubble departure. In this paper, an impro...
Computational Model for Determination of Vapor Bubble Growth Speed Maximimum in Superheated Liquids
Turlajs, D; Grivcovs, V; Jaundālders, S
2008-01-01
The article deals with problems of investigation of boiling of liquids - the most intensive heat transfer phenomena. There are a lot of research reports of boiling dynamics and heat transfer studies in micro- and nanoscopic levels, carried out with aim to understund phenomena in whole. Simulation model for vapor bubble growth speed maximum calculations in the region of bubble critical radius is elaborated. Numerical calculations were carried out for water in the region of system p...
Zhang, L. C.; Zhu, X. L.; Huang, Y. F.; Liu, Z.; Yan, K.
2016-07-01
In this paper, a simple model was developed to predict the dynamics of a spark-induced bubble under different ambient pressures. This work helps in developing a deep-towed plasma sparker, as the model can predict the dynamics of bubbles subjected to very high ambient pressures (about 20 MPa) which normally are difficult to obtain experimentally. Experimental results indicate that the maximum bubble radius for a fixed discharge energy decreases as a power-law function of the ambient pressure up to 1.0 MPa; the bubble period also decreases quickly with increasing ambient pressure. For a constant value of the ratio of bubble energy to discharge energy, the modeling results for both maximum radius and bubble period are in good agreement with the experimental results. Both sets of results indicate that the bubble period is proportional to the maximum radius under different ambient pressures.
Micro-bubble Drag Reduction on a High Speed Vessel Model
Yanuar; Gunawan; Sunaryo; A. Jamaluddin
2012-01-01
Ship hull form of the underwater area strongly influences the resistance of the ship.The major factor in ship resistance is skin friction resistance.Bulbous bows,polymer paint,water repellent paint (highly water-repellent wall),air injection,and specific roughness have been used by researchers as an attempt to obtain the resistance reduction and operation efficiency of ships.Micro-bubble injection is a promising technique for lowering frictional resistance.The injected air bubbles are supposed to somehow modify the energy inside the turbulent boundary layer and thereby lower the skin friction.The purpose of this study was to identify the effect of injected micro bubbles on a navy fast patrol boat (FPB) 57 m type model with the following main dimensions:L=2 450 mm,B=400 mm,and T=190 mm.The influence of the location of micro bubble injection and bubble velocity was also investigated.The ship model was pulled by an electric motor whose speed could be varied and adjusted.The ship model resistance was precisely measured by a load cell transducer.Comparison of ship resistance with and without micro-bubble injection was shown on a graph as a function of the drag coefficient and Froude number.It was shown that micro bubble injection behind the mid-ship is the best location to achieve the most effective drag reduction,and the drag reduction caused by the micro-bubbles can reach 6％-9％.
Barnaveli, A T; Barnaveli, Andro; Gogberashvili, Merab
1995-01-01
We investigate the gravitational behavior of spherical domain walls (bubbles) arising during the phase transitions in the early Universe. In the thin-wall approximation we show the existence of the new solution of Einstein equations with negative gravitational mass of bubbles and the reversed direction of time flow on the shell. This walls exhibit gravitational repulsion just as the planar walls are assumed to do. The equilibrium radius and critical mass of such objects are found for realistic models.
Mario Forni; Luca Gambetti; Marco Lippi; Luca Sala
2014-01-01
We introduce noisy information into a standard present value stock price model. Agents receive a noisy signal about the structural shock driving future dividend variations. The resulting equilibrium stock price includes a transitory component — the "noise bubble" — which can be responsible for boom and bust episodes unrelated to economic fundamentals. We propose a non-standard VAR procedure to estimate the structural shock and the "noise" shock, their impulse response functions and the bubble...
The universal applicability of the bubble model is examined and attention is drawn to a number of situations in which it is inadequate. Some possible areas of improvements are considered. The effects of van der Waals forces on the bubble radius are estimated to be insignificant, but the problem of the potential well depth and shape is more intractable. We find that when a linear combination of finite potential components is used, the Ps pressure can become independent of well depth when the infinite component becomes large, as may be the case in certain liquids. We have developed the idea of Tao on the description of the wavefunction overlap occurring in a thin skin in the inside of the bubble surface. By considering the contribution from the protrusion of hydrogen atoms from a hard core bubble surface, we calculate that the pick-off rate can be significantly altered. (orig.)
Numerical Modeling of the Photothermal Processing for Bubble Forming around Nanowire in a Liquid
Anis Chaari
2014-01-01
Full Text Available An accurate computation of the temperature is an important factor in determining the shape of a bubble around a nanowire immersed in a liquid. The study of the physical phenomenon consists in solving a photothermic coupled problem between light and nanowire. The numerical multiphysic model is used to study the variations of the temperature and the shape of the created bubble by illumination of the nanowire. The optimization process, including an adaptive remeshing scheme, is used to solve the problem through a finite element method. The study of the shape evolution of the bubble is made taking into account the physical and geometrical parameters of the nanowire. The relation between the sizes and shapes of the bubble and nanowire is deduced.
PERTURBANCE ANALYSIS AND FISSION MODEL OF COLLAPSING BUBBLE
Petrík, P.; Zima, Patrik; Ridha, A.
Vol. 1. Liberec: Technical University of Liberec, 2010 - (Vít, T.; Dančová, P.), s. 507-514 ISBN 978-80-7372-670-6. [International Conference Experimental Fluid Mechanics 2010. Liberec (CZ), 24.11.2010-26.11.2010] R&D Projects: GA ČR GAP101/10/1428 Institutional research plan: CEZ:AV0Z20760514 Keywords : cavitation bubble * collapse * fission Subject RIV: BK - Fluid Dynamics
Modelling of Methane Bubble Dissolution in Freshwater Lakes
Štiková, L.; Stanovský, Petr; Orvalho, Sandra; Kordač, M.
Bratislava: Slovak Society of Chemical Engineering, 2015 - (Markoš, J.), s. 161 ISBN 978-80-89475-14-8. [International Conference of Slovak Society of Chemical Engineering /42./. Tatranské Matliare (SK), 25.05.2015-29.05.2015] R&D Projects: GA ČR GAP504/12/1186 Institutional support: RVO:67985858 Keywords : methane ebullition * mass transfer * single bubbles Subject RIV: EH - Ecology, Behaviour
Some Misconceptions in Bubble and Drops Hydrodynamics Models
Wichterle, Kamil; Smutná, K.; Večeř, M.; Růžička, Marek
Bratislava: Slovak University of Technology, 2010 - (Markoš, J.), s. 288 ISBN 978-80-227-3290-1. [International Conference of Slovak Society of Chemical Engineering /37./. Tatranské Matliare (SK), 24.05.2010-28.05.2010] R&D Projects: GA ČR GA104/07/1110 Institutional research plan: CEZ:AV0Z40720504 Keywords : bubbles * surface tension * viscosity Subject RIV: CI - Industrial Chemistry, Chemical Engineering
袁德文; 潘良明; 陈德奇
2009-01-01
The process of bubble growth on heating wall in subcooled boiling includes the micro-layer evaporation on heating wall and the bubble top coagulation when the bubbles grow to a certain size and emerge into the subcooled mainstream fluid. Based on this consideration, a model for the single bubble growth of subcooled flow boiling in vertical narrow rectangular channel was proposed. Compared with experimental results, the error of the simulation results using the proposed model is less than ?5%. The simulation results indicated that as the wall superheat increases, the bubble growth gets faster, with the subcooled degree of mainstream increases, the bubble growth in later stage would be slowed, with the contact angle increases, the contact radius of the bubble bottom and the wall tension would be strengthened, resulting in faster bubble growth to make the bubble to be flat and more easily exposed to the mainstream. The velocity of mainstream has no significant effects on bubble growth rate.%@@ 引言 沸腾具有极高的换热效率的原因就在于沸腾过程中大量气相对周围液体的扰动.对汽泡形成、生长、脱离等过程机理的了解是理解和优化沸腾换热过程的关键.
Kochetkov, I. I.; Pinaev, A. V.
2013-03-01
Strong shock and detonation waves in inert and chemically active bubble media, which are generated by a wire explosion initiated by a capacitor with a stored energy W_0 =12.3-1,600 J, is experimentally studied. The measurements are performed near the wire and far from the wire in a vertical shock tube 4.5 m long with a volume fraction of the gas in the medium β _0 =1-4 %. It is shown that in inert bubble medium, a short intensely decaying shock wave (SW) with intense pressure oscillations is formed in the vicinity of wire explosion point; near the explosion point at β _0 le 2 % the SW propagates with the velocity of sound in a liquid. In chemically active bubble medium, an unsteady detonation wave generated by a wire explosion is formed. The pressure amplitude and the velocity of this wave are greater and the length is smaller than those of SW in an inert bubble medium in the same range of explosion energy. It is found that in the interval of low energy explosion from {˜ }12 to 64 J, the formation of the bubble detonation wave occurs faster than that at high energies (3× 102-103 J).
A mathematical model of particles entrainment in the freeboard of bubbling fluidized bed
A two-component flow with a low concentration occurs in the freeboard after the eruption of the bubbles at the free surface of the bubbling fluidized beds. A three-dimensional mathematical model of gas-particle mixture turbulent flow is developed in this work. The gas turbulence is modeled using standard k-ε turbulence model. The dispersed phase is treated by the Lagrangian approach. Coupling between the gas phase and the dispersed particles is modeled by adding, source term in the momentum equation for gas phase. By using the equation for determining the drag coefficient, the particle's shape is involved in the model, as well. Experimental investigations for determining the conditions at the bed surface, the origin of erupting bubbles and their erupting mechanism and for selecting dominant influencing parameters are done. Comparisons of the experimentally obtained results with the results of other authors, which are in reasonable agreement, are used for determining the dominant influencing parameter related to the phenomena of bursting bubbles. On a base of observation, the dynamic eruption mechanism of single bursting bubble is proposed. This bubble eruption mechanism is used like a boundary condition in the numerical experimentation for investigation of the entrainment of particles and its separation along the freeboard height. The numerical solving of the developed mathematical model is accomplished by using the CFD technology. For graphical design and mesh generation of the flow domain and for numerical solving of the equations of the developed mathematical model, the software packages Gambit and FLUENT are used, respectively. The testing and verification of the proposed erupting bubble mechanism and the developed mathematical model for two-component flow in the freeboard, is made by numerical experimentation in 3D cylindrical flow domain, in the following conditions: eruption of isolated central bursting bubble; determining of particles terminal velocity
Proussevitch, Alexander
2014-05-01
Parameterization of volcanic ash transport and dispersion (VATD) models strongly depends on particle morphology and their internal properties. Shape of ash particles affects terminal fall velocities (TFV) and, mostly, dispersion. Internal density combined with particle size has a very strong impact on TFV and ultimately on the rate of ash cloud thinning and particle sedimentation on the ground. Unlike other parameters, internal particle density cannot be measured directly because of the micron scale sizes of fine ash particles, but we demonstrate that it varies greatly depending on the particle size. Small simple type ash particles (fragments of bubble walls, 5-20 micron size) do not contain whole large magmatic bubbles inside and their internal density is almost the same as that of volcanic glass matrix. On the other side, the larger compound type ash particles (>40 microns for silicic fine ashes) always contain some bubbles or the whole spectra of bubble size distribution (BSD), i.e. bubbles of all sizes, bringing their internal density down as compared to simple ash. So, density of the larger ash particles is a function of the void fraction inside them (magmatic bubbles) which, in turn, is controlled by BSD. Volcanic ash is a product of the fragmentation of magmatic foam formed by pre-eruptive bubble population and characterized by BSD. The latter can now be measured from bubble imprints on ash particle surfaces using stereo-scanning electron microscopy (SSEM) and BubbleMaker software developed at UNH, or using traditional high-resolution X-Ray tomography. In this work we present the mathematical and statistical formulation for this problem connecting internal ash density with particle size and BSD, and demonstrate how the TFV of the ash population is affected by variation of particle density.
Gleiser, Marcello
1994-01-01
The evolution of spherically symmetric unstable scalar field configurations (``bubbles'') is examined for both symmetric (SDWP) and asymmetric (ADWP) double-well potentials. Bubbles with initial static energies $E_0\\la E_{{\\rm crit}}$, where $E_{{\\rm crit}}$ is some critical value, shrink in a time scale determined by their linear dimension, or ``radius''. Bubbles with $E_0\\ga E_{{\\rm crit}}$ evolve into time-dependent, localized configurations which are {\\it very} long-lived compared to characteristic time-scales in the models examined. The stability of these configurations is investigated and possible applications are briefly discussed.tic time-scales in the models examined. The stability of these configurations is investigated and possible applications are briefly discussed.
Hydrodynamics of Bubble Bouncing on a Wall, Experiment and Modelling
Zedníková, Mária; Vejražka, Jiří; Růžička, Marek; Drahoš, Jiří
Praha: Process Engineering Publisher, 2008, s. 878. ISBN 978-80-02-02050-9. [18th International Congress of Chemical and Process Engineering CHISA 2008. Praha (CZ), 24.08.2008-28.08.2008] R&D Projects: GA ČR GA104/07/1110; GA AV ČR(CZ) KJB200720801; GA AV ČR(CZ) IAA200720801 Institutional research plan: CEZ:AV0Z40720504 Keywords : bubble-wall collision * impact velocity * contact time Subject RIV: CI - Industrial Chemistry, Chemical Engineering
Bubble burst as jamming phase transition
Nishinari, Katsuhiro; Saito, Yukiko Umeno; Watanabe, Tsutomu
2010-01-01
Recently research on bubble and its burst attract much interest of researchers in various field such as economics and physics. Economists have been regarding bubble as a disorder in prices. However, this research strategy has overlooked an importance of the volume of transactions. In this paper, we have proposed a bubble burst model by focusing the transactions incorporating a traffic model that represents spontaneous traffic jam. We find that the phenomenon of bubble burst shares many similar properties with traffic jam formation by comparing data taken from US housing market. Our result suggests that the transaction could be a driving force of bursting phenomenon.
A force balance model for the motion, impact, and bounce of bubbles
Klaseboer, Evert; Manica, Rogerio; Hendrix, Maurice H. W.; Ohl, Claus-Dieter; Chan, Derek Y. C.
2014-09-01
A force balance model has been developed to predict the terminal velocity of a sub-millimetric bubble as its rises in water under buoyancy. The dynamics of repeated collisions and rebounds of the bubble against a horizontal solid surface is modeled quantitatively by including forces due to buoyancy, added mass, drag, and hydrodynamic lubrication—the last arises from the drainage of water trapped in the thin film between the solid surface and the surface of the deformable bubble. The result is a self-contained, parameter-free model that is capable of giving quantitative agreement with measured trajectories and observed collisions and rebounds against a solid surface as well as the spatio-temporal evolution of the thin film during collision as measured by interferometry.
Bubble nucleation in $\\phi^{4}$ models at all temperatures
Ferrera, A
1995-01-01
One possible way in which phase transitions in the early universe may have ocurred is via nucleation of bubbles of the new phase (true vacuum) in the old phase (false vacuum). The technique most widely used to compute the probability of bubble nucleation is based on instanton methods in the context of the semiclassical approximation. At zero temperature in (3+1) dimensions the nucleation rate is dominated by the O(4) symmetric instanton, a sphere of radius R, while at temperatures T >> R^{-1}, the decay is dominated by a ``cylindrical'' (static) instanton wtih O(3) invariance. There has been discussion in the literature as to whether the transition between these two regimens would be first order (discontinuity in the first derivative of the nucleation rate at the transition temperature T_{c}), or second order (continuity of the first derivative, but discontinuity of the second derivative at T_{c}). In this paper we obtain the finite temperature solutions corresponding to the quantum and the thermal regimes, a...
Li, Yulan; Hu, Shenyang Y.; Montgomery, Robert O.; Gao, Fei; Sun, Xin
2012-05-30
Experiments show that inter-granular and intra-granular gas bubbles have different growth kinetics which results in heterogeneous gas bubble microstructures in irradiated nuclear fuels. A science-based model predicting the heterogeneous microstructure evolution kinetics is desired, which enables one to study the effect of thermodynamic and kinetic properties of the system on gas bubble microstructure evolution kinetics and morphology, improve the understanding of the formation mechanisms of heterogeneous gas bubble microstructure, and provide the microstructure to macroscale approaches to study their impact on thermo-mechanical properties such as thermo-conductivity, gas release, volume swelling, and cracking. In our previous report 'Mesoscale Benchmark Demonstration, Problem 1: Mesoscale Simulations of Intra-granular Fission Gas Bubbles in UO2 under Post-irradiation Thermal Annealing', we developed a phase-field model to simulate the intra-granular gas bubble evolution in a single crystal during post-irradiation thermal annealing. In this work, we enhanced the model by incorporating thermodynamic and kinetic properties at grain boundaries, which can be obtained from atomistic simulations, to simulate fission gas bubble growth kinetics in polycrystalline UO2 fuels. The model takes into account of gas atom and vacancy diffusion, vacancy trapping and emission at defects, gas atom absorption and resolution at gas bubbles, internal pressure in gas bubbles, elastic interaction between defects and gas bubbles, and the difference of thermodynamic and kinetic properties in matrix and grain boundaries. We applied the model to simulate gas atom segregation at grain boundaries and the effect of interfacial energy and gas mobility on gas bubble morphology and growth kinetics in a bi-crystal UO2 during post-irradiation thermal annealing. The preliminary results demonstrate that the model can produce the equilibrium thermodynamic properties and the morphology of gas
Modeling of reaction kinetics in bubbling fluidized bed biomass gasification reactor
Thapa, R.K.; Halvorsen, B.M. [Telemark University College, Kjolnes ring 56, P.O. Box 203, 3901 Porsgrunn (Norway); Pfeifer, C. [University of Natural Resources and Life Sciences, Vienna (Austria)
2013-07-01
Bubbling fluidized beds are widely used as biomass gasification reactors as at the biomass gasification plant in Gussing, Austria. The reactor in the plant is a dual circulating bubbling fluidized bed gasification reactor. The plant produces 2MW electricity and 4.5MW heat from the gasification of biomass. Wood chips as biomass and olivine particles as hot bed materials are fluidized with high temperature steam in the reactor. As a result, biomass undergoes endothermic chemical reaction to produce a mixture of combustible gases in addition to some carbon-dioxide (CO2). The combustible gases are mainly hydrogen (H2), carbon monoxide (CO) and methane (CH4). The gas is used to produce electricity and heat via utilization in a gas engine. Alternatively, the gas is further processed for gaseous or liquid fuels, but still on the process of development level. Composition and quality of the gas determine the efficiency of the reactor. A computational model has been developed for the study of reaction kinetics in the gasification rector. The simulation is performed using commercial software Barracuda virtual reactor, VR15. Eulerian-Lagrangian approach in coupling of gas-solid flow has been implemented. Fluid phase is treated with an Eulerian formulation. Discrete phase is treated with a Lagrangian formulation. Particle-particle and particle-wall interactions and inter-phase heat and mass transfer have been taken into account. Series of simulations have been performed to study model prediction of the gas composition. The composition is compared with data from the gasifier at the CHP plant in Güssing, Austria. The model prediction of the composition of gases has good agreements with the result of the operating plant.
Bubble and kink solitons in the φ6-model of nonlinear field theory
We have studied the φ6-model in the parameter domain A>1, with A being the relevant parameter of the model. For this case we have found localized soliton-like solutions: kinks and bubbles. The investigation of waves propagating through a stable vacuum shows that the sound velocity provides a rigid constraint for these oscillations to be stable or not. (orig.)
Mechanism of the Dark Matter and Condensed Bubble Objects Formation in the Model of Extended Space
Andreev, V A
2004-01-01
Within the framework of Extended Space Model (ESM) the processes connected to birth of photons in a gravitational field are studied. These photons have a nonzero mass. It can be both positive, and negative, and photon's energy and strength of the gravitational field determine its absolute value. It is shown that in ESM model formation of bubble gravitational structures is possible. In the frame of ESM one can obtain the follow physical picture. Bubble gravitational objects have a halo formed by dark matter generated by photons with a positive mass. The photons with a negative mass are throw away in free deep space and create there antigravitating vacuum with negative pressure. The comparison ESM bubble structures with similar objects of a type "gravastar", existing in a General Theory of Relativity (GR) is made.
Modeling of Liquid Level and Bubble Behavior in Vacuum Chamber of RH Process
Yi-hong LI; Yan-ping BAO; Rui WANG; Min WANG; Qing-xue HUANG; Yu-gui LI
2016-01-01
In the Ruhrstahl-Heraeus (RH)refining process,liquid steel flow pattern in a ladle is controlled by the fluid flow behavior in the vacuum chamber.Potassium chloride solution and NaOH solution saturated with CO 2 were respectively used as a tracer to investigate the liquid and gas flow behaviors in the vacuum chamber.Principal compo-nent and comparative analysis were made to show the factors controlling mixing and circulation flow rate.The liquid level and bubble behavior in the vacuum chamber greatly affect fluid flow in RH process.Experiments were per-formed to investigate the effects of liquid steel level,gas flow rate,bubble residence time,and gas injection mode on mixing,decarburization,and void fraction.The results indicate that the mixing process can be divided into three re-gions:the flow rate-affected zone,the concentration gradient-affected zone,and their combination.The liquid steel level in the vacuum chamber of 300 mm is a critical point in the decarburization transition.For liquid level lower than 300 mm,liquid steel circulation controls decarburization,while for liquid level higher than 300 mm,bubble behavior is the main controlling factor.During the RH process,it is recommended to use the concentrated bubble injection mode for low gas flow rates and the uniform bubble injection mode for high gas flow rates.
Modelling chemical reactions in dc plasma inside oxygen bubbles in water
Takeuchi, N.; Ishii, Y.; Yasuoka, K.
2012-02-01
Plasmas generated inside oxygen bubbles in water have been developed for water purification. Zero-dimensional numerical simulations were used to investigate the chemical reactions in plasmas driven by dc voltage. The numerical and experimental results of the concentrations of hydrogen peroxide and ozone in the solution were compared with a discharge current between 1 and 7 mA. Upon increasing the water vapour concentration inside bubbles, we saw from the numerical results that the concentration of hydrogen peroxide increased with discharge current, whereas the concentration of ozone decreased. This finding agreed with the experimental results. With an increase in the discharge current, the heat flux from the plasma to the solution increased, and a large amount of water was probably vaporized into the bubbles.
Bubble spreading during the boiling crisis: modelling and experimenting in microgravity
Nikolayev, Vadim; Garrabos, Y; Lecoutre, C; Chatain, D
2016-01-01
Boiling is a very efficient way to transfer heat from a heater to the liquid carrier. We discuss the boiling crisis, a transition between two regimes of boiling: nucleate and film boiling. The boiling crisis results in a sharp decrease in the heat transfer rate, which can cause a major accident in industrial heat exchangers. In this communication, we present a physical model of the boiling crisis based on the vapor recoil effect. Under the action of the vapor recoil the gas bubbles begin to spread over the heater thus forming a germ for the vapor film. The vapor recoil force not only causes its spreading, it also creates a strong adhesion to the heater that prevents the bubble departure, thus favoring the further spreading. Near the liquid-gas critical point, the bubble growth is very slow and allows the kinetics of the bubble spreading to be observed. Since the surface tension is very small in this regime, only microgravity conditions can preserve a convex bubble shape. In the experiments both in the Mir spa...
A vapour bubble collapse model to describe the fragmentation of low-melting materials
By means of a model, the fragmentation of a hot melt of metal in consequence of collapsing vapour-bubbles is investigated. In particular the paper deals with the development of the physical model-ideas for calculation of the temperature of contact that adjusts between the temperature of the melt and the coolant, of the waiting-time until bubble-nucleation occurs and of the maximal obtainable vapour-bubble-radius in dependence of the coolant-temperature. After that follows the description of the computing-program belonging to this model and of the results of an extensive parameter-study. The study examined the influence of the temperature of melt and coolant, the melted mass, the nucleation-site-density, the average maximum bubble-radius, the duration of film-breakdown and the coefficient of heat-transition. The calculation of the process of fragmentation turns out to be according to expectation, whereas the duration of this process seems to be somewhat too long. The dependence of the surface-enlargement on the subcooling of the water-bath and the initial temperature of the melt is not yet reproduced satisfactorily by the model. The reasons for this are the temperature-increase of the water-bath as well as the fact that the coupling of heat-flux-density and nucleation-site-density are not taken into consideration. Further improvement of the model is necessary and may improve the results in the sense of the experimental observations. (orig.)
Mahood Hameed B.
2016-01-01
Full Text Available An analytical model for the convective heat transfer coefficient and the two-phase bubble size of a three-phase direct contact heat exchanger was developed. Until the present, there has only been a theoretical model available that deals with a single two-phase bubble and a bubble train condensation in an immiscible liquid. However, to understand the actual heat transfer process within the three-phase direct contact condenser, characteristic models are required. A quasi - steady energy equation in a spherical coordinate system with a potential flow assumption and a cell model configuration has been simplified and solved analytically. The convective heat transfer in terms of Nu number has been derived, and it was found to be a function to Pe number and a system void fraction. In addition, the two-phase bubble size relates to the system void fraction and has been developed by solving a simple energy balance equation and using the derived convective heat transfer coefficient expression. Furthermore, the model correlates well with previous experimental data and theoretical results.
A gas bubble-based parallel micro manipulator: conceptual design and kinematics model
The parallel mechanism has become an alternative solution when micro manipulators are demanded in the fields of micro manipulation and micro assembly. In this technical note, a three-degree-of-freedom (3-DOF) parallel micro manipulator is presented, which is directly driven by three micro gas bubbles. Since the micro gas bubbles are generated and maintained due to the surface tension between the gas and liquid media, the proposed novel system can be used in the liquid environment which allows for rotation about the X and Y axes and translation along the Z axis. In this technical note, the conceptual design of micro gas bubble-based parallel manipulator is introduced and the input/output characteristic of the actuator is analyzed in detail. The kinematics model of the parallel micro manipulator is also established, based on which the workspace and the system motion resolution are analyzed as a criterion and reference for future prototype development. (technical note)
A MODEL OF MIRA'S COMETARY HEAD/TAIL ENTERING THE LOCAL BUBBLE
We model the cometary structure around Mira as the interaction of an asymptotic giant branch stellar wind from Mira A with a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment, Mira entered the Local Bubble (low-density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well collimated for a >100 kyr timescale. The result is a broad-head/narrow-tail structure that resembles the observed morphology of Mira's comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.
A model of Mira's cometary head/tail entering the Local Bubble
Esquivel, A; Canto, J; Rodriguez-Gonzalez, A; Lopez-Camara, D; Velazquez, P F; De Colle, F
2010-01-01
We model the cometary structure around Mira as the interaction of an AGB wind from Mira A, and a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment Mira entered the Local Bubble (low density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well collimated for a 100 kyr timescale. The result is a broad-head/narrow-tail structure that resembles the observed morphology of Mira's comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.
Javed Bin Kamal
2012-09-01
Full Text Available The paper aims at constructing an optimal portfolio by applying Sharpe’s single index model of capital asset pricing in different scenarios, one is ex ante stock price bubble scenario and stock price bubble and bubble burst is second scenario. Here we considered beginning of year 2010 as rise of stock price bubble in Dhaka Stock Exchange. Hence period from 2005 -2009 is considered as ex ante stock price bubble period. Using DSI (All share price index in Dhaka Stock Exchange as market index and considering daily indices for the March 2005 to December 2009 period, the proposed method formulates a unique cut off point (cut off rate of return and selects stocks having excess of their expected return over risk-free rate of return surpassing this cut-off point. Here, risk free rate considered to be 8.5% per annum (Treasury bill rate in 2009. Percentage of an investment in each of the selected stocks is then decided on the basis of respective weights assigned to each stock depending on respective ‘β’ value, stock movement variance representing unsystematic risk, return on stock and risk free return vis-à-vis the cut off rate of return. Interestingly, most of the stocks selected turned out to be bank stocks. Again we went for single index model applied to same stocks those made to the optimum portfolio in ex ante stock price bubble scenario considering data for the period of January 2010 to June 2012. We found that all stocks failed to make the pass Single Index Model criteria i.e. excess return over beta must be higher than the risk free rate. Here for the period of 2010 to 2012, the risk free rate considered to be 11.5 % per annum (Treasury bill rate during 2012.
The dynamic behaviour of SonoVue microbubbles, a new generation ultrasound contrast agent, is investigated in real time with light scattering method. Highly diluted SonoVue microbubbles are injected into a diluted gel made of xanthan gum and water. The responses of individual SonoVue bubbles to driven ultrasound pulses are measured. Both linear and nonlinear bubble oscillations are observed and the results suggest that SonoVue microbubbles can generate strong nonlinear responses. By fitting the experimental data of individual bubble responses with Sarkar's model, the shell coating parameter of the bubbles and dilatational viscosity is estimated to be 7.0 nm·s·Pa
Numerical Modeling of Air-Water Flows in Bubble Columns and Airlift Reactors
Studley, Allison F
2010-01-01
Bubble columns and airlift reactors were modeled numerically to better understand the hydrodynamics and analyze the mixing characteristics for each configuration. An Eulerian-Eulerian approach was used to model air as the dispersed phase within a continuous phase of water using the commercial software FLUENT. The Schiller-Naumann drag model was employed along with virtual mass and the standard k-e turbulence model. The equations were discretized using the QUICK scheme and solved with the SIMP...
Effect of Water Vapour to Temperature Inside Sonoluminescing Bubble
安宇; 谢崇国; 应崇福
2003-01-01
Using the model based on the homo-pressure approximation, we explain why the maximum temperature is sensitive to the ambient temperature in the single bubble sonoluminescence. The numerical simulation shows that the maximum temperature inside a sonoluminescing bubble depends on how much water vapour evaporates or coagulates at the bubble wall during the bubble shrinking to its minimum size. While the amount of water vapour inside the bubble at the initial and the final state of the compression depends on the saturated water vapour pressure which is sensitive to the ambient temperature. The lower the saturated vapour pressure is, the higher the maximum temperature is. This may lead to more general conclusion that those liquids with lower saturated vapour pressure are more favourable for the single bubble sonoluminescence. We also compare those bubbles with different noble gases, the result shows that the maximum temperatures in the different gas bubbles are almost the same for those with the same ambient temperature.
Bubble generation in a twisted and bent DNA-like model
Larsen, Peter Ulrik Vingaard; Christiansen, Peter Leth; Bang, Ole;
2004-01-01
The DNA molecule is modeled by a parabola embedded chain with long-range interactions between twisted base pair dipoles. A mechanism for bubble generation is presented and investigated in two different configurations. Using random normally distributed initial conditions to simulate thermal fluctu...
Hong, Ban Zhen; Keong, Lau Kok; Shariff, Azmi Mohd
2016-05-01
The employment of different mathematical models to address specifically for the bubble nucleation rates of water vapour and dissolved air molecules is essential as the physics for them to form bubble nuclei is different. The available methods to calculate bubble nucleation rate in binary mixture such as density functional theory are complicated to be coupled along with computational fluid dynamics (CFD) approach. In addition, effect of dissolved gas concentration was neglected in most study for the prediction of bubble nucleation rates. The most probable bubble nucleation rate for the water vapour and dissolved air mixture in a 2D quasi-stable flow across a cavitating nozzle in current work was estimated via the statistical mean of all possible bubble nucleation rates of the mixture (different mole fractions of water vapour and dissolved air) and the corresponding number of molecules in critical cluster. Theoretically, the bubble nucleation rate is greatly dependent on components' mole fraction in a critical cluster. Hence, the dissolved gas concentration effect was included in current work. Besides, the possible bubble nucleation rates were predicted based on the calculated number of molecules required to form a critical cluster. The estimation of components' mole fraction in critical cluster for water vapour and dissolved air mixture was obtained by coupling the enhanced classical nucleation theory and CFD approach. In addition, the distribution of bubble nuclei of water vapour and dissolved air mixture could be predicted via the utilisation of population balance model.
This paper describes the updating of the sub-cooled boiling model with the more recent and better sub-models. The improved sub-models include: Hibiki and Ishii [1] correlation for nucleation site density, Kocamustafaogullari [2] correlation for bubble departure diameter and the S-gamma model of Lo and Rao [3] for bubble size distribution in the flow. The new model has been tested against measured data from Debora [4] and Bartolomei [5]. The results show that improvement in the bubble size prediction has the most significant impact on the accuracy of the model. (author)
Dynamics of micro-bubble sonication inside a phantom vessel
Qamar, Adnan
2013-01-10
A model for sonicated micro-bubble oscillations inside a phantom vessel is proposed. The model is not a variant of conventional Rayleigh-Plesset equation and is obtained from reduced Navier-Stokes equations. The model relates the micro-bubble oscillation dynamics with geometric and acoustic parameters in a consistent manner. It predicts micro-bubble oscillation dynamics as well as micro-bubble fragmentation when compared to the experimental data. For large micro-bubble radius to vessel diameter ratios, predictions are damped, suggesting breakdown of inherent modeling assumptions for these cases. Micro-bubble response with acoustic parameters is consistent with experiments and provides physical insight to the micro-bubble oscillation dynamics.
How do the models to measure bubbles in the stock market
Ding, Yun
2007-01-01
Bubble economy is a phenomenon frequently occurs in the process of finance opening and an essential factor contributing to financial crisis. Once bubbles burst, financial crisis is quite possible. Stock market bubble is a key component of a countrys economic bubble, thus it is important to analyze the developing process of stock market bubble and to find out the method to eliminate harmful bubble component. The thesis organizes a static index system to measure the bubble component in a st...
The Accretion Wind Model of the Fermi Bubbles (II): Radiation
Mou, Guobin; Gan, Zhaoming; Sun, Mouyuan
2015-01-01
In a previous work, we have shown that the formation of the Fermi bubbles can be due to the interaction between winds launched from the hot accretion flow in Sgr A* and the interstellar medium (ISM). In that work, we focus only on the morphology. In this paper we continue our study by calculating the gamma-ray radiation. Some cosmic ray protons (CRp) and electrons must be contained in the winds, which are likely formed by physical processes such as magnetic reconnection. We have performed MHD simulations to study the spatial distribution of CRp, considering the advection and diffusion of CRp in the presence of magnetic field. We find that a permeated zone is formed just outside of the contact discontinuity between winds and ISM, where the collisions between CRp and thermal nuclei mainly occur. The decay of neutral pions generated in the collisions, combined with the inverse Compton scattering of background soft photons by the secondary leptons generated in the collisions and primary CR electrons can well expl...
Bubble merger model for the nonlinear Rayleigh-Taylor instability driven by a strong blast wave
A bubble merger model is presented for the nonlinear evolution of the Rayleigh-Taylor instability driven by a strong blast wave. Single bubble motion is determined by an extension of previous buoyancy-drag models extended to the blast-wave-driven case, and a simple bubble merger law in the spirit of the Sharp-Wheeler model allows for the generation of larger scales. The blast-wave-driven case differs in several respects from the classical case of incompressible fluids in a uniform gravitational field. Because of material decompression in the rarefaction behind the blast front, the asymptotic bubble velocity and the merger time depend on time as well as the transverse scale and the drive. For planar blast waves, this precludes the emergence of a self-similar regime independent of the initial conditions. With higher-dimensional blast waves, divergence restores the properties necessary for the establishment of the self-similar state, but its establishment requires a very high initial characteristic mode number and a high Mach number for the incident blast wave
A hierarchy of simple hyperbolic two-fluid models for bubbly flows
Drui, Florence; Kokh, Samuel; Massot, Marc
2016-01-01
With the objective of modeling both separate and disperse two-phase flows, we use in this paper a methodology for deriving two-fluid models that do not assume any flow topology. This methodology is based on a variational principle and on entropy dissipation requirement. Some of the models that are such derived and studied are already known in the contexts of the description of separate-or disperse-phase flows. However, we here propose an arrangement of these models into a hierarchy based on their links through relaxation parameters. Moreover, the models are shown to be compatible with the description of a monodisperse bubbly flow and, within this frame, the relaxation parameters can be identified. This identification is finally verified and discussed through comparisons with experimental measures of sound dispersion and with dispersion relations of a reference model for bubbly media.
Prediction of a subcooled boiling flow with mechanistic wall boiling and bubble size models
Subcooled boiling is one of the crucial phenomena for the design, operation and safety analysis of a nuclear power plant. In recent years, developers of multiphase CFD (Computational Fluid Dynamics) codes focused their development activity on the mechanistic prediction of DNB (Departure from Nucleate Boiling) in PWR. Wall boiling model is one of the key parameters for this purpose. In order to enhance prediction capability of the subcooled boiling flow, an advanced wall boiling model consisting of a mechanistic bubble departure model (Klausner et al., 1993), Hibiki et al.'s (2009) active nucleate site model and Cole's bubble departure frequency model was explored for the CFD code. To ensure a wide range applicability of the advanced wall boiling model, each model was evaluated separately according to the flow conditions such as pressure, temperature and flow rate. Finally, the advanced wall boiling model was implemented into the STAR-CD as a form of user FORTRAN file. One of the other important parameters for an accurate prediction of the subcooled boiling flow is bubble size which governs interfacial transfer terms between two phases. In this study, the S-gamma model, which was developed for the STAR-CD (Lo, 2006), was applied as a bubble size model. For the validation of the present wall boiling and bubble size models, benchmark calculations were carried out against SUBO and DEBORA subcooled boiling flow data. Working fluid of SUBO test is steam/water and its pressure condition is about 2 bars. In contrast to this, working fluid of DEBORA test is Refrigerant-12 (R-12) and phasic density ratio of the tests is equivalent to that of steam/water around 90 to 170 bars. Therefore, present benchmark calculation covers wide range pressure condition of steam/water. The calculation results confirms that the new mechanistic wall boiling and bubble size models follow well the tendency on the change of flow conditions and they can be applicable to the wide range of flow
Multifluid modeling of the desulfurization process within a bubbling fluidized bed coal gasifier
Armstrong, L-M.; Gu, Sai; Luo, K. H.; P Mahanta
2013-01-01
The desulfurization process to a two-dimensional (2-D) and three-dimensional (3-D) Eulerian-Eulerian computational fluid dynamic (CFD) model of a coal bubbling fluidized gasifier is introduced. The desulfurization process is important for the reduction of harmful SOx emissions; therefore, the development of a CFD model capable of predicting chemical reactions involving desulfurization is key to the optimization of reactor designs and operating conditions. To model the process, one gaseous pha...
To enhance the multi-dimensional analysis capability for a subcooled boiling two-phase flow, the one-group interfacial area transport equation was improved with a source term for the bubble lift-off. It included the bubble lift-off diameter model and the lift-off frequency reduction factor model. The bubble lift-off diameter model took into account the bubble's sliding on a heated wall after its departure from a nucleate site, and the lift-off frequency reduction factor was derived by considering the coalescences of the sliding bubbles. To implement the model, EAGLE (elaborated analysis of gas-liquid evolution) code was developed for a multi-dimensional analysis of two-phase flow. The developed model and EAGLE code were validated with the experimental data of SUBO (subcooled boiling) and SNU (Seoul National University) test, where the subcooled boiling phenomena in a vertical annulus channel were observed. Locally measured two-phase flow parameters included a void fraction, interfacial area concentration, and bubble velocity. The results of the computational analysis revealed that the interfacial area transport equation with the bubble lift-off model showed a good agreement with the experimental results of SUBO and SNU. It demonstrates that the source term for the wall nucleation by considering a bubble sliding and lift-off mechanism enhanced the prediction capability for the multi-dimensional behavior of void fraction or interfacial area concentration in the subcooled boiling flow. From the point of view of the bubble velocity, the modeling of an increased turbulence induced by boiling bubbles at the heated wall enhanced the prediction capability of the code.
Thomas, Johannes
2014-01-01
An extended analytical model for particle dynamics in fields of a highly-nonlinear plasma wake field (the bubble or blow out regime) is derived. A recently proposed piecewise model (Kostyukov et al., New J. Phys., {\\bf 12}, 045009 (2010)) is generalized to include a time dependent bubble radius and full field solution in the acceleration direction. Incorporation of the cavity dynamics in the model is required to simulate the particle trapping properly. On the other hand, it is shown that the previously reported piecewise model does not reproduce the formation of a mono energetic peak in the particle spectrum. The mono energetic electron beams are recovered only when the full longitudinal field gradient is included in the model.
On the modelling of bubble entrainment by impinging jets in CFD-simulations
This contribution presents different approaches for the modeling of air entrainment under water by plunging jets in CFD codes. In simulations which include the full length of the jet and its environment, the process of bubble generation cannot be resolved due to computational limitations. This is why the air entrainment has to be modeled in meso-scale simulations. In the frame of an Euler- Euler simulation, the local morphology of the phases has to be considered in the drag model. In the impinging jet configuration, the air is a continuous phase above the water level but bubbly below the water level. Various drag models are implemented in the CFD solver CFX11 and their influence on the gas void fraction below the water level is discussed. The algebraic interface area density (AIAD) model applies a drag coefficient for bubbles and a different drag coefficient for the free surface. If the AIAD model is used for the simulation of impinging jets, the gas entrainment depends on the free parameters included in this model. The calculated gas entrainment can be adapted via these parameters. Therefore, an advanced AIAD approach could be used in future for the implementation of models (e.g. correlations) for the gas entrainment. (authors)
Predictions of bubbly flows in vertical pipes using two-fluid models in CFDS-FLOW3D code
Banas, A.O.; Carver, M.B. [Chalk River Laboratories (Canada); Unrau, D. [Univ. of Toronto (Canada)
1995-09-01
This paper reports the results of a preliminary study exploring the performance of two sets of two-fluid closure relationships applied to the simulation of turbulent air-water bubbly upflows through vertical pipes. Predictions obtained with the default CFDS-FLOW3D model for dispersed flows were compared with the predictions of a new model (based on the work of Lee), and with the experimental data of Liu. The new model, implemented in the CFDS-FLOW3D code, included additional source terms in the {open_quotes}standard{close_quotes} {kappa}-{epsilon} transport equations for the liquid phase, as well as modified model coefficients and wall functions. All simulations were carried out in a 2-D axisymmetric format, collapsing the general multifluid framework of CFDS-FLOW3D to the two-fluid (air-water) case. The newly implemented model consistently improved predictions of radial-velocity profiles of both phases, but failed to accurately reproduce the experimental phase-distribution data. This shortcoming was traced to the neglect of anisotropic effects in the modelling of liquid-phase turbulence. In this sense, the present investigation should be considered as the first step toward the ultimate goal of developing a theoretically sound and universal CFD-type two-fluid model for bubbly flows in channels.
PROGRESS TOWARDS MODELING OF FISCHER TROPSCH SYNTHESIS IN A SLURRY BUBBLE COLUMN REACTOR
Donna Post Guillen; Tami Grimmett; Anastasia M. Gandrik; Steven P. Antal
2010-11-01
The Hybrid Energy Systems Testing (HYTEST) Laboratory is being established at the Idaho National Laboratory to develop and test hybrid energy systems with the principal objective to safeguard U.S. Energy Security by reducing dependence on foreign petroleum. A central component of the HYTEST is the slurry bubble column reactor (SBCR) in which the gas-to-liquid reactions will be performed to synthesize transportation fuels using the Fischer Tropsch (FT) process. SBCRs are cylindrical vessels in which gaseous reactants (for example, synthesis gas or syngas) is sparged into a slurry of liquid reaction products and finely dispersed catalyst particles. The catalyst particles are suspended in the slurry by the rising gas bubbles and serve to promote the chemical reaction that converts syngas to a spectrum of longer chain hydrocarbon products, which can be upgraded to gasoline, diesel or jet fuel. These SBCRs operate in the churn-turbulent flow regime which is characterized by complex hydrodynamics, coupled with reacting flow chemistry and heat transfer, that effect reactor performance. The purpose of this work is to develop a computational multiphase fluid dynamic (CMFD) model to aid in understanding the physico-chemical processes occurring in the SBCR. Our team is developing a robust methodology to couple reaction kinetics and mass transfer into a four-field model (consisting of the bulk liquid, small bubbles, large bubbles and solid catalyst particles) that includes twelve species: (1) CO reactant, (2) H2 reactant, (3) hydrocarbon product, and (4) H2O product in small bubbles, large bubbles, and the bulk fluid. Properties of the hydrocarbon product were specified by vapor liquid equilibrium calculations. The absorption and kinetic models, specifically changes in species concentrations, have been incorporated into the mass continuity equation. The reaction rate is determined based on the macrokinetic model for a cobalt catalyst developed by Yates and Satterfield [1]. The
In this paper, a computational simulation study is presented on the prediction of helium bubble evolution during repair welding of irradiated 304 stainless steel. Realistic spatial and temporal temperature and stress evolution during welding were obtained from simulation of the repair welding operation using the finite element model approach. The helium bubble evolution model by Kawano et al. was adopted as a user subroutine in the finite element model to predict the spatial distribution and temporal evolution of the helium bubble size and density in the heat-affected zone (HAZ) of partial penetration welds. Comparisons with experimental results available in open literature show that the predicted average helium bubble sizes were consistent with those observed experimentally under similar conditions. In addition, the computer simulation revealed strong spatial variation of helium bubble size due to the differences in combined thermal and stress conditions experienced in different locations in the HAZ. The predicted location of the maximum helium bubble agreed well with the observed helium-induced cracking site. The effect of welding heat input and welding speed was also investigated numerically. The modeling approach adopted in this study could be used as a cost-effective tool to quantitatively correlate the welding condition, radiation damage, and the likelihood of cracking, under the influence of welding-induced thermal and stress cycles. The model will also be useful in studying the degradation of properties from helium bubble formation of post-welded structures, even if a successful weld is made. (authors)
Bubble Growth in Lunar Basalts
Zhang, Y.
2009-05-01
Although Moon is usually said to be volatile-"free", lunar basalts are often vesicular with mm-size bubbles. The vesicular nature of the lunar basalts suggests that they contained some initial gas concentration. A recent publication estimated volatile concentrations in lunar basalts (Saal et al. 2008). This report investigates bubble growth on Moon and compares with that on Earth. Under conditions relevant to lunar basalts, bubble growth in a finite melt shell (i.e., growth of multiple regularly-spaced bubbles) is calculated following Proussevitch and Sahagian (1998) and Liu and Zhang (2000). Initial H2O content of 700 ppm (Saal et al. 2008) or lower is used and the effect of other volatiles (such as carbon dioxide, halogens, and sulfur) is ignored. H2O solubility at low pressures (Liu et al. 2005), concentration-dependent diffusivity in basalt (Zhang and Stolper 1991), and lunar basalt viscosity (Murase and McBirney 1970) are used. Because lunar atmospheric pressure is essentially zero, the confining pressure on bubbles is completely supplied by the overlying magma. Due to low H2O content in lunar basaltic melt (700 ppm H2O corresponds to a saturation pressure of 75 kPa), H2O bubbles only grow in the upper 16 m of a basalt flow or lake. A depth of 20 mm corresponds to a confining pressure of 100 Pa. Hence, vesicular lunar rocks come from very shallow depth. Some findings from the modeling are as follows. (a) Due to low confining pressure as well as low viscosity, even though volatile concentration is very low, bubble growth rate is extremely high, much higher than typical bubble growth rates in terrestrial melts. Hence, mm-size bubbles in lunar basalts are not strange. (b) Because the pertinent pressures are so low, bubble pressure due to surface tension plays a main role in lunar bubble growth, contrary to terrestrial cases. (c) Time scale to reach equilibrium bubble size increases as the confining pressure increases. References: (1) Liu Y, Zhang YX (2000) Earth
Implementation of a new bubbly-slug interphase drag model in RELAP5/MOD2
The implementation of a new bubbly-slug interphase drag model in the RELAP5/MOD2 code is described. The model is based on the determination of an effective interphase drag coefficient from a set of best-estimate void fraction correlations covering the full range of geometries and flow conditions encountered in PWR safety analysis. Calculations are reported which show that the new model leads to a much better prediction of void fraction profile for low flows in rod bundles than the standard model. Further work is necessary to derive a model formulation which can be guaranteed to produce physical drag coefficients in all flow situations. (author)
Modeling of Bubble Column Slurry Reactor for Dimethyl Ether Synthesis from Syngas
张海涛; 应卫勇; 房鼎业
2005-01-01
A mathematical model for a bubble column slurry reactor is presented for dimethyl ether synthesis from syngas. Methanol synthesis from carbon monoxide and carbon dioxide by hydrogenation and the methanol dehydration are considered as independent reactions, in which methanol, dimethyl ether and carbon dioxide are the key components. In this model, the gas phase is considered to be in plug flow and the liquid phase to be in partly back mixing with axial distribution of solid catalyst. The simulation results show that the axial dispersion of solid catalysts, the operational height of the slurry phase in the bubble column slurry reactor, and the reaction results are influenced by the reaction temperature and pressure, which are the basic data for the scale-up of reactor.
Prediction of bubble detachment diameter in flow boiling based on force analysis
Highlights: ► All the forces acting on the growing bubbles are taken into account in the model. ► The bubble contact diameter has significant effect on bubble detachment. ► Bubble growth force and surface tension are more significant in narrow channel. ► A good agreement between the predicted and the measured results is achieved. - Abstract: Bubble detachment diameter is one of the key parameters in the study of bubble dynamics and boiling heat transfer, and it is hard to be measured in a boiling system. In order to predict the bubble detachment diameter, a theoretical model is proposed based on forces analysis in this paper. All the forces acting on a bubble are taken into account to establish a model for different flow boiling configurations, including narrow and conventional channels, upward, downward and horizontal flows. A correlation of bubble contact circle diameter is adopted in this study, and it is found that the bubble contact circle diameter has significant effect on bubble detachment. A new correlation taking the bubble contact circle diameter into account for the evaluation of bubble growth force is proposed in this study, and it is found that the bubble growth force and surface tension force are more significant in narrow channel when comparing with that in conventional channel. A visual experiment was carried out in order to verify present model; and the experimental data from published literature are used also. A good agreement between predicted and measured results is achieved.
Soap bubble hadronic states in a QCD-motivated Nambu-Jona-Lasinio model
Kutnii, Sergii
2015-01-01
Inhomogeneous solutions of the gap equation in the mean field approach to Nambu-Jona-Lasinio model are studied. An approximate Ginzburg-Landau-like gap equation is obtained and the domain wall solution is found. Binding of fermions to the domain wall is demonstrated. Compact domain wall with bound fermions is studied and stabilisation by fermion pressure is demonstrated which opens a possibility for existence of "soap bubble" hadronic states.
On the Dynamics of a Fluid-Particle Interaction Model: The Bubbling Regime
Carrillo, Jose' A.; Karper, Trygve; Trivisa, Konstantina
2010-01-01
This article deals with the issues of global-in-time existence and asymptotic analysis of a fluid-particle interaction model in the so-called bubbling regime. The mixture occupies the physical space $\\Omega \\subset \\mathbb{R}^3$ which may be unbounded. The system under investigation describes the evolution of particles dispersed in a viscous compressible fluid and is expressed by the conservation of fluid mass, the balance of momentum and the balance of particle density often referred as the ...
Modeling and simulation of the bubble-induced flow in wine fermentation vessels
Schmidt Dominik; Velten Kai
2015-01-01
Detailed flow pattern analyses regarding wine fermentations conducted without mechanical agitation are limited to lab-scale investigations, as industrial size measurements are expensive and difficult to realize. Computational fluid dynamic (CFD) methods can offer an alternative and more flexible approach to gain insight into such bubble induced fluid flows. Therefore, the aim of this study was to transfer the findings of existing research onto a CFD model capable of capturing the three- dimen...
Implementation of two-phase tritium models for helium bubbles in HCLL breeding blanket modules
Fradera, Jordi; Sedano, L.A.; Mas de les Valls Ortiz, Elisabet; Batet Miracle, Lluís
2011-01-01
Tritium self-sufficiency requirement of future DT fusion reactors involves large helium production rates in the breeding blankets; this might impact on the conceptual design of diverse fusion power reactor units, such as Liquid Metal (LM) blankets. Low solubility, long residence-times and high production rates create the conditions for Helium nucleation, which could mean effective T sinks in LM channels. A model for helium nano-bubble formation and tritium conjugate transport phen...
Modelling of bubble-mediated gas transfer: Fundamental principles and a laboratory test
Woolf, D.K.; Leifer, I.S.; Nightingale, P.D.; Rhee, T.S.; Bowyer, P.; Caulliez, G.; Leeuw, G. de; Larsen, S.E.; Liddicoat, M.; Baker, J.; Andreae, M.O.
2007-01-01
The air-water exchange of gases can be substantially enhanced by wave breaking and specifically by bubble-mediated transfer. A feature of bubble-mediated transfer is the additional pressure on bubbles resulting from the hydrostatic forces on a submerged bubble and from surface tension and curvature.
Debasish Majumdar; Kamakshya Prasad Modak; Subhendu Rakshit
2016-02-01
We propose a two-component dark matter (DM) model, each component of which is a real singlet scalar, to explain results from both direct and indirect detection experiments. We put the constraints on the model parameters from theoretical bounds, PLANCK relic density results and direct DM experiments. The -ray flux is computed from DM annihilation in this framework and is then compared with the Fermi-LAT observations from galactic centre region and Fermi bubble.
Passive Underwater Noise Attenuation Using Large Encapsulated Air Bubbles.
Lee, Kevin M; Wochner, Mark S; Wilson, Preston S
2016-01-01
Measurements demonstrating low-frequency underwater sound attenuation using arrays of large, tethered, stationary encapsulated bubbles to surround a sound source were compared with various effective medium models for the acoustic dispersion relationship in bubbly liquids. Good agreement was observed between measurements for the large bubbles (on the order of 10 cm) at frequencies below 1 kHz and a model originally intended to describe the acoustic behavior of ultrasound contrast agents. The primary goal is to use the model for designing encapsulated-bubble-based underwater noise abatement systems and to reduce uncertainty in system performance. PMID:26611010
Hendricks, Vincent Fella; Pedersen, David Budtz
2013-01-01
Much like the trade and trait sof bubbles in financial markets,similar bubbles appear on the science market. When economic bubbles burst, the drop in prices causes the crash of unsustainable investments leading to an investor confidence crisis possibly followed by a financial panic. But when...... bubbles appear in science, truth and reliability are the first victims. This paper explores how fashions in research funding and research management may turn science into something like a bubble economy....
Sonochemical effects on single-bubble sonoluminescence
Yuan, L
2005-01-01
A refined hydro-chemical model for single-bubble sonoluminescence is presented. The processes of water vapor evaporation and condensation, mass diffusion, and chemical reactions are taken into account. Numerical simulations of Xe-, Ar- and He-filled bubbles are carried out. The results show that the trapped water vapor in conjunction with its endothermic chemical reactions significantly reduces the temperature within the bubble so that the degrees of ionization are generally very low. The chemical radicals generated from water vapor are shown to play an increasingly important role in the light emission from Xe to He bubbles. Light spectra and pulses computed from an optically thin model and from an essentially blackbody model are compared with recent experimental results. It is found that the results of the blackbody model generally match better with the experiment ones than those of the optically thin model. Suggestions on how to reconcile the conflict are given.
Bubbles tomorrow and bubbles yesterday, but never bubbles today?
Williams, John C.
2013-01-01
Standard asset price models have generally failed to detect bubbles, with enormous costs to the economy. Economists are now creating promising new models that account for bubbles by relaxing the assumption of rational expectations and allowing people’s decisions to be driven by their perceptions of what the future may hold. ; This letter is adapted from a presentation by the president and CEO of the Federal Reserve Bank of San Francisco to the National Association for Business Economics in Sa...
Shape of isolated bubble in downwardly inclined gas-liquid two-phase flows
Experimental and theoretical studies on the shape of a single bubble in downwardly inclined gas-liquid two-phase flows are carried out. Measurements of the shape are made using conductance probes and visual observation in this paper. The experimental results show that bubble turning phenomenon of bubble shape happens at low Froude number, which is defined by gas/liquid mixture velocity. The turning phenomenon significantly affects the bubble velocity. It is also found that the shape of bubble is mainly determined by the Froude number rather than the bubble length. A hydrodynamics model for the shape of bubble tail is derived based on a one-dimension two-fluid model. The model well predicts the characteristics of the bubble tail when compared with the experimental results. The theoretical analysis shows that the critical Froude number, at which the turning phenomenon happens, increases as the downward inclination and pipe diameter increase. (authors)
Lindenbaum, S J
2008-01-01
In an earlier paper we developed a QCD inspired theoretical parton bubble model (PBM) for RHIC/LHC. The PBM quantitatively agreed with the strong charged particle pair correlations observed by the STAR collaboration at RHIC in the highest energy Au + Au central collisions, and also agreed with the Hanbury Brown and Twiss (HBT) observed small final state source size approximately 2f radii in the transverse momentum range above 0.8 GeV/c. The model assumed a substructure of a ring of localized adjoining 2f radius bubbles(gluonic hot spots) perpendicular to the collider beam direction, centered on the beam, at mid-rapidity and located on the expanding fireball surface of the Au + Au collisions. In this paper we extend the model (PBME) to include the changing development of bubbles with centrality from the most central region where bubbles are very important to the most peripheral where the bubbles are gone. Energy density is found to be related to bubble formation and as centrality decreases the maximum energy d...
TU Juan; GUAN J.F.; MATULA T.J.; Crum L.A.; WEI Rong-jue
2008-01-01
The dynamic behaviour of SonoVue microbubbles a new generation ultrasound contrast agent is investigated in real time with light scattering method.Highly diluted SonoVue microbubbles are injected into a diluted gel made of xanthan gum and water.The responses of individual SonoVue bubbles to driven ultrasound pulses are measured.Both linear and nonlinear bubble oscillations are observed and the results suggest that SonoVue microbubbles can generate strong nonlinear responses.By fitting the experimental data of individual bubble responses with Sarkar's model,the shell coating parameter of the bubbles and dilatational viscosity is estimated to be 7.0 nm·s·Pa.
Models for viscosity and shear localization in bubble-rich magmas
Vona, Alessandro; Ryan, Amy G.; Russell, James K.; Romano, Claudia
2016-09-01
Bubble content influences magma rheology and, thus, styles of volcanic eruption. Increasing magma vesicularity affects the bulk viscosity of the bubble-melt suspension and has the potential to promote non-Newtonian behavior in the form of shear localization or brittle failure. Here, we present a series of high temperature uniaxial deformation experiments designed to investigate the effect of bubbles on the magma bulk viscosity. The starting materials are cores of natural rhyolitic obsidian synthesized to have variable vesicularity (ϕ = 0- 66%). The foamed cores were deformed isothermally (T = 750 °C) at atmospheric conditions using a high-temperature uniaxial press under constant displacement rates (strain rates between 0.5- 1 ×10-4 s-1) and to total strains of 10-40%. The viscosity of the bubble-free melt (η0) was measured by micropenetration and parallel plate methods to establish a baseline for experiments on the vesicle rich cores. At the experimental conditions, rising vesicle content produces a marked decrease in bulk viscosity that is best described by a two-parameter empirical equation: log10 ηBulk =log10 η0 - 1.47[ ϕ / (1 - ϕ) ] 0.48. Our parameterization of the bubble-melt rheology is combined with Maxwell relaxation theory to map the potential onset of non-Newtonian behavior (shear localization) in magmas as a function of melt viscosity, vesicularity, and strain rate. For low degrees of strain (i.e. as in our study), the rheological properties of vesicular magmas under different flow types (pure vs. simple shear) are indistinguishable. For high strain or strain rates where simple and pure shear viscosity values may diverge, our model represents a maximum boundary condition. Vesicular magmas can behave as non-Newtonian fluids at lower strain rates than unvesiculated melts, thereby, promoting shear localization and (explosive or non-explosive) magma fragmentation. The extent of shear localization in magma influences outgassing efficiency
Models for viscosity and shear localization in bubble-rich magmas
Vona, Alessandro; Ryan, Amy G.; Russell, James K.; Romano, Claudia
2016-09-01
Bubble content influences magma rheology and, thus, styles of volcanic eruption. Increasing magma vesicularity affects the bulk viscosity of the bubble-melt suspension and has the potential to promote non-Newtonian behavior in the form of shear localization or brittle failure. Here, we present a series of high temperature uniaxial deformation experiments designed to investigate the effect of bubbles on the magma bulk viscosity. The starting materials are cores of natural rhyolitic obsidian synthesized to have variable vesicularity (ϕ = 0- 66%). The foamed cores were deformed isothermally (T = 750 °C) at atmospheric conditions using a high-temperature uniaxial press under constant displacement rates (strain rates between 0.5- 1 ×10-4 s-1) and to total strains of 10-40%. The viscosity of the bubble-free melt (η0) was measured by micropenetration and parallel plate methods to establish a baseline for experiments on the vesicle rich cores. At the experimental conditions, rising vesicle content produces a marked decrease in bulk viscosity that is best described by a two-parameter empirical equation: log10 ηBulk =log10 η0 - 1.47[ ϕ / (1 - ϕ) ] 0.48. Our parameterization of the bubble-melt rheology is combined with Maxwell relaxation theory to map the potential onset of non-Newtonian behavior (shear localization) in magmas as a function of melt viscosity, vesicularity, and strain rate. For low degrees of strain (i.e. as in our study), the rheological properties of vesicular magmas under different flow types (pure vs. simple shear) are indistinguishable. For high strain or strain rates where simple and pure shear viscosity values may diverge, our model represents a maximum boundary condition. Vesicular magmas can behave as non-Newtonian fluids at lower strain rates than unvesiculated melts, thereby, promoting shear localization and (explosive or non-explosive) magma fragmentation. The extent of shear localization in magma influences outgassing efficiency
The Ising model and bubbles in the quark-gluon plasma
Svetitsky, B
1997-01-01
I review evidence for the stability of bubbles in the quark-gluon plasma near the confinement phase transition. In analogy with the much-studied oil-water emulsions, this raises the possibility that there are many phases between the pure plasma and the pure hadron gas, characterized by spontaneous inhomogeneity and modulation. In studying emulsions, statistical physicists have reproduced many of their phases with microscopic models based on Ising-like theories with competing interactions. Hence we seek an effective Ising Hamiltonian for the SU(3) gauge theory near its transition.
Highlights: • RANS simulation of laminar separation bubbles. • Spalart–Allamaras unexpected capability. • Straightforward implementation of our SA modifications. • Applications of a high order DG incompressible solver. - Abstract: The present paper deals with the Reynolds Averaged Navier–Stokes (RANS) simulation of Laminar Separation Bubble (LSB). This phenomenon is of large interest in several engineering fields, such as the study of wind turbines, unmanned aerial vehicles (UAV) and micro-air vehicles (MAV) characterized by a low operating Reynolds number. In such contexts a laminar boundary layer separation followed by a turbulent transition and afterwards by a turbulent reattachment may appear in the flow-field. The main novelty of this work is that an almost standard Spalart–Allmaras (SA) model, without additional equations for transition modeling, was successfully employed. The result achieved is very surprising being the model not developed for this purpose, but for fully-turbulent flows or for cases with imposed transition location. This result is of large interest, since the SA model is widely used in commercial, open-source and research codes. However, our approach cannot be advocated to predict natural transition within an attached boundary layer, indeed it is only able to deal with transitions triggered by a separated flow. The reliability and accuracy of our approach are here proved computing, by means of a high-order Discontinuous Galerkin (DG) incompressible solver, the flow-field over two airfoils at different flow regimes showing the formation of a LSB
Insights in hydrodynamics of bubbling fluidized beds at elevated pressure by DEM-CFD approach
Zahra Mansourpour; Sedighe Karimi; Reza Zarghami; Navid Mostoufi; Rahmat Sotudeh-Gharebagh
2010-01-01
A numerical simulation was conducted to study the effect of pressure on bubble dynamics in a gas-solid fluidized bed. The gas flow was modeled using the continuum theory and the solid phase, by the dis-crete element method (DEM). To validate the simulation results, calculated local pressure fluctuations were compared with corresponding experimental data of 1-mm polyethylene particles. It was shown that the model successfully predicts the hydrodynamic features of the fluidized bed as observed in the experiments. Influence of pressure on bubble rise characteristics such as bubble rise path, bubble sta-bility, average bubbles diameter and bubble velocity through the bed was investigated. The simulation results are in conformity with current hydrodynamic theories and concepts for fluidized beds at high pressures. The results show further that elevated pressure reduces bubble growth, velocity and stability and enhances bubble gyration through the bed, leading to change in bed flow structure.
Measurements of sub-surface bubble populations and the modelling of air-sea gas flux
Coles, David Geoffrey Hallstaff
2010-01-01
Bubbles, formed by breaking waves, play an important role in the transfer of gases between the Earth’s oceans and atmosphere and have been shown to increase the flux of gases during periods of heightened sea state. Having been formed, these bubble clouds evolve through the effects of buoyancy, gas exsolution and dissolution, and the fragmentation and coalescence of bubbles. A number of experimenters have successfully measured sub-surface bubble clouds using a variety of acoustic and opt...
Primary Particles from different bubble generation techniques
Butcher, A. C.; King, S. M.; Rosenoern, T.; Nilsson, E. D.; Bilde, M.
2011-12-01
Sea spray aerosols (SSA) are of major interest to global climate models due to large uncertainty in their emissions and ability to form Cloud Condensation Nuclei (CCN). In general, SSA are produced from wind breaking waves that entrain air and cause bubble bursting on the ocean surface. Preliminary results are presented for bubble generation, bubble size distribution, and CCN activity for laboratory generated SSA. In this study, the major processes of bubble formation are examined with respect to particle emissions. It has been suggested that a plunging jet closely resembles breaking wave bubble entrainment processes and subsequent bubble size distributions (Fuentes, Coe et al. 2010). Figure 1 shows the different particle size distributions obtained from the various bubble generation techniques. In general, frits produce a higher concentration of particles with a stronger bimodal particle size distribution than the various jet configurations used. The experiments consist of a stainless steel cylinder closed at both ends with fittings for aerosol sampling, flow connections for the recirculating jet, and air supply. Bubble generation included a recirculating jet with 16 mm or 4 mm nozzles, a stainless steel frit, or a ceramic frit. The chemical composition of the particles produced via bubble bursting processes has been probed using particle CCN activity. The CCN activity of sodium chloride, artificial sea salt purchased from Tropic Marin, and laboratory grade artificial sea salt (Kester, Duedall et al. 1967) has been compared. Considering the the limits of the shape factor as rough error bars for sodium chloride and bubbled sea salt, the CCN activity of artificial sea salt, Tropic Marin sea salt, and sodium chloride are not significantly different. This work has been supported by the Carlsberg Foundation.
Jones, T. B.; Bliss, G. W.
1977-01-01
The theoretical principles related to bubble dielectrophoresis are examined, taking into account the polarization force, aspects of bubble deformation, the electrostatic bubble levitation theorem, and the equation of motion. The measurement of the dielectrophoretic force on static and dynamic bubbles represents a convenient experimental method for the study of the general problem of dielectrophoresis. The experiments reported include static-force measurements, static-levitation experiments, and dynamic-force measurements.
Sonochemistry and bubble dynamics.
Mettin, Robert; Cairós, Carlos; Troia, Adriano
2015-07-01
The details of bubble behaviour in chemically active cavitation are still not sufficiently well understood. Here we report on experimental high-speed observations of acoustically driven single-bubble and few-bubble systems with the aim of clarification of the connection of their dynamics with chemical activity. Our experiment realises the sonochemical isomerization reaction of maleic acid to fumaric acid, mediated by bromine radicals, in a bubble trap set-up. The main result is that the reaction product can only be observed in a parameter regime where a small bubble cluster occurs, while a single trapped bubble stays passive. Evaluations of individual bubble dynamics for both cases are given in form of radius-time data and numerical fits to a bubble model. A conclusion is that a sufficiently strong collapse has to be accompanied by non-spherical bubble dynamics for the reaction to occur, and that the reason appears to be an efficient mixing of liquid and gas phase. This finding corroborates previous observations and literature reports on high liquid phase sonochemical activity under distinct parameter conditions than strong sonoluminescence emissions. PMID:25194210
Long scale evolution of a nonlinear stochastic dynamic system for modeling market price bubbles
Kiselev, S. A.; Phillips, Andy; Gabitov, I.
2000-07-01
This Letter investigates the stochastic dynamics of a simplified agent-based microscopic model describing stock market evolution. Our mathematical model includes a stochastic market and a sealed-bid double auction. The dynamics of the model are determined by the game of two types of traders: (i) `intelligent' traders whose strategy is based on nonlinear technical data analysis 1 and (ii) `random' traders that act without a consistent strategy. We demonstrate the effect of time-scale separations on the market dynamics. We study the characteristics of the market relaxation in response to perturbations caused by large cash flows generated between these two groups of traders. We also demonstrate that our model exhibits the formation of a price bubble 2 and the subsequent transition to a bear market 3. Bear market - a macroscopically long stage of a market evolution when the stock price declines significantly, 15% or more.
Lu, Yurong; Wang, Zhongtong; Yong, Huadong; Zhou, Youhe
2016-07-01
Due to the larger current-carrying property, Bi2Sr2CaCu2Ox (Bi2212) superconductors have a great potential application in high field magnet. Bi2212 superconducting material can be fabricated as an isotropic round wire. However, there is 30% void space in the wire, such as gas bubbles. The void space has a larger influence on the property of the wire. In this paper, we will study the effect of gas bubble on the fracture behavior. Based on the double cantilever beam model and critical state theory, the mechanical behavior of Bi2212 wire is studied for decreasing field. Two different damage mechanisms are discussed using the strain energy release rate and strain of bridge. The results show that the large gas bubble can increase the strain of bridge. The central filaments with gas bubble are easier to be damaged than the edge filaments with gas bubble.
Liquid-bubble Interaction under Surf Zone Breaking Waves
Derakhti, M.; Kirby, J. T., Jr.
2014-12-01
Liquid-bubble interaction, especially in complex two-phase bubbly flow under breaking waves, is still poorly understood. Derakhti and Kirby (2014a,b) have recently studied bubble entrainment and turbulence modulation by dispersed bubbles under isolated unsteady breaking waves along with extensive model verifications and convergence tests. In this presentation, we continue this examination with attention turned to the simulation of periodic surf zone breaking waves. In addition, the relative importance of preferential accumulation of dispersed bubbles in coherent vortex cores is investigated. Heavier-than-liquid particles, i.e. sediment, tend to accumulate in regions of high strain rate and avoid regions of intense vorticity. In contrast, lighter-than-liquid particles such as bubbles tend to congregate in vortical regions. We perform a three dimensional (3D) large-eddy simulation (LES) using a Navier-Stokes solver extended to incorporate entrained bubble populations, using an Eulerian-Eulerian formulation for the polydisperse bubble phase. The volume of fluid (VOF) method is used for free surface tracking. The model accounts for momentum exchange between dispersed bubbles and liquid phase as well as bubble-induced dissipation. We investigate the formation and evolution of breaking-induced turbulent coherent structures (BTCS) under both plunging and spilling periodic breaking waves as well as BTCS's role on the intermittent 3D distributions of bubble void fraction in the surf zone. We particularly examine the correlation between bubble void fractions and Q-criterion values to quantify this interaction. Also, the vertical transport of dispersed bubbles by downburst type coherent structures in the transition region is compared to that by obliquely descending eddies. All the results are summarized at different zones from outer to inner surf zone.
Simulation of hydrogen bubble growth in tungsten by a hybrid model
A two dimensional hybrid code (HIIPC-MC) joining rate-theory and Monte Carlo (MC) methods is developed in this work. We evaluate the cascade-coalescence mechanism contribution to the bubble growth by MC. First, effects of the starting radius and solute deuterium concentration on the bubble growth are studied; then the impacts of the wall temperature and implantation ion flux on the bubble growth are assessed. The simulation indicates that the migration-coalescence of the bubbles and the high pressure inside the bubbles are the main driving forces for the bubble growth, and that neglect of the migration and coalescence would lead to an underestimation of the bubble growth or blistering
Potential and energy of the monoenergetic electrons in an alternative ellipsoid bubble model
The electron acceleration in the bubble regime is considered during the intense laser-plasma interaction. The presented ellipsoid cavity model is more consistent than the previous spherical model, and it explains the monoenergetic electron trajectory more accurately. At the relativistic region, the maximum energy of electrons in the ellipsoid model is about 24% more than the spherical model, and this is confirmed by PIC and the measured experimental results reported here. The electron energy spectrum is also calculated, and it is found that the energy distribution ratio of electrons ΔE/E for the ellipsoid model in the here reported condition is about 11% which is less than the one third that of the spherical model. It is in good agreement with the experimentally measured value in the same condition. In this regime, the parameters of the quasi-monoenergetic electrons output beam can be described more appropriately. In this work, 10 TW from 16.6 TW, 500 mJ, and 30-fs laser pulse was focused on the best matched point above a 2-mm-diameter pulsed He gas jet to obtain a stable ellipsoid bubble. Laser intensity of 1.42x1019 W cm-2 corresponding to a normalized vector potential of a0=2.6 focused in a 100-μm2 spot at the focal point and 1 mm above the edge of the gas jet with an electron density of 1x1019 cm-3 accelerates electrons to the relativistic velocities. The obtained monoenergetic electron energy spectrum is properly explained by the ellipsoid model introduced here.
Bubble rupture in bubble electrospinning
Chen Rouxi
2015-01-01
Full Text Available As the distinctive properties and different applications of nanofibers, the demand of nanofibers increased sharply in recently years. Bubble electrospinning is one of the most effective and industrialized methods for nanofiber production. To optimize the set-up of bubble electrospinning and improve its mass production, the dynamic properties of un-charged and charged bubbles are studied experimentally, the growth and rupture process of a bubble are also discussed in this paper.
Measurements of local interfacial area concentration in two-phase bubbly flow
The local interfacial characteristics of an air-water, bubbly upflow in a circular pipe has been investigated experimentally cased on the measurements of a miniature two-sensor resistivity probe. Radial profiles of interfacial area concentration (IAC), together with other structural parameters (void fraction, bubble frequency, bubble velocity and Sauter mean bubble diameter) were simultaneously measured at four axial positions with entrance length (L)-to-internal diameter (D) ratios of 30, 60, 90 and 120. The experiments were carried out under various fixed gas and liquid fluxes, with only the bubble size being changed at the flow entrance. It is found that the phase distribution ana IAC are very sensitive to the variation oi the bubble size and the bubble coalescence effects during the development of bubbly flow. This paper summarizes the experimental results on the associated local interfacial parameters and compares the present data with the existing models used in predicting the IAC. (author)
Extended two-fluid model applied to analysis of bubbly flow in multiphase rotodynamic pump impeller
Zhiyi YU; Guoyu WANG; Shuliang CAO
2009-01-01
This paper presents an extended two-fluid model based on the Navier-Stokes equations and the standard k-e turbulence model, to simulate the three-dimensional air-water bubbly flow in turbo machinery. In the governing equations, the drag force and added mass force are added and the additional source terms arising from fluctuations of gas volume fraction are considered. The discrete equations are solved using a developed two-phase semi-implicit method for pressure-linked equations, consistent (SIMPLEC) algorithm in body-fitted coordi-nates with a staggered grid system. Simulation is then carried out for the pure liquid flow and air-water two-phase flow with the inlet gas volume fraction being 15% in a multiphase rotodynamic pump impeller and the pump head performance is predicted. Comparison with experimental results shows the reliability and commonality of the numerical model.
A study on models of the inhomogeneity of the decompression in a bubble chamber
Before building a hydrogen bubble chamber with liquid decompression the 'Saturne' cyclotron department wished to study for this chamber a shape leading to a homogeneous decompression as far as possible, without the production of vortices even after prolonged operation. The 'Office National d'Etudes et de Recherches Aeronautiques' (ONERA) were ready to carry out experiments on a model by strioscopy. The model was filled with air but an attempt was made to simulate the actual conditions as far as possible by varying the speed of the piston. The model was placed at one end of a tunnel, at the other end of which were produced alternatively compression and decompression waves. The study made it possible to conclude that it was necessary to make the base of the chamber round and that, in the space between the decompression cylinder and the body of the chamber it was advantageous to use 5 fins instead of 3. (author)
Full text of publication follows: A computational code has been developed for inert gas behavior in the primary system of sodium cooled fast reactor (SFR). The inert gas exists in the circulating coolant sodium as resolved gas or free gas bubbles. The primary coolant system of the SFR has free surfaces which are covered by argon gas. The cover gas system is slightly pressurized above the atmospheric pressure. Therefore, the argon cover gas would dissolve in the liquid sodium and dispersed in the primary coolant system by diffusion and convection. In addition, helium gas is generated in the control rods of the reactor core and is emitted as small bubbles from the B4C pellet. Another source of gas bubbles is gas entrainment at the free surface of the reactor vessel if the sodium flow velocity at the free surface is large enough. From the viewpoint of design and safety of the SFR, the gas in the primary system may cause disturbance in reactivity, nucleation site for boiling and cavitation, flow instability, and/or influence on heat transfer. In the present method, the concentration distribution of the dissolved gas and free gas bubble are evaluated for steady state and transient conditions according to the mass and heat transfer in the cooling system. The bubble diameter is assumed to follow the lognormal distribution and discretized into numbers of representative diameters. Physical modeling of the gas behavior is described in the following. In the intermediate heat exchanger, nucleation of bubbles caused by the temperature decrease is evaluated at the heat transfer tube surface. The bubbles on the tube surface are detached from the wall according to the balance of drag, surface tension and buoyancy forces. At the gas liquid interface (free surface and bubble surface), evaporation, dissolution and diffusion are calculated. Bubbles traveling the primary system would break up at turbulent region such as core fuel assembly and mechanical coolant pump. At the free
McDeavitt, Sean [Texas A & M Univ., College Station, TX (United States); Shao, Lin [Texas A & M Univ., College Station, TX (United States); Tsvetkov, Pavel [Texas A & M Univ., College Station, TX (United States); Wirth, Brian [Univ. of Tennessee, Knoxville, TN (United States); Kennedy, Rory [Idaho National Lab. (INL), Idaho Falls, ID (United States)
2014-04-07
Advanced fast reactor systems being developed under the DOE's Advanced Fuel Cycle Initiative are designed to destroy TRU isotopes generated in existing and future nuclear energy systems. Over the past 40 years, multiple experiments and demonstrations have been completed using U-Zr, U-Pu-Zr, U-Mo and other metal alloys. As a result, multiple empirical and semi-empirical relationships have been established to develop empirical performance modeling codes. Many mechanistic questions about fission as mobility, bubble coalescience, and gas release have been answered through industrial experience, research, and empirical understanding. The advent of modern computational materials science, however, opens new doors of development such that physics-based multi-scale models may be developed to enable a new generation of predictive fuel performance codes that are not limited by empiricism.
Advanced fast reactor systems being developed under the DOE's Advanced Fuel Cycle Initiative are designed to destroy TRU isotopes generated in existing and future nuclear energy systems. Over the past 40 years, multiple experiments and demonstrations have been completed using U-Zr, U-Pu-Zr, U-Mo and other metal alloys. As a result, multiple empirical and semi-empirical relationships have been established to develop empirical performance modeling codes. Many mechanistic questions about fission as mobility, bubble coalescience, and gas release have been answered through industrial experience, research, and empirical understanding. The advent of modern computational materials science, however, opens new doors of development such that physics-based multi-scale models may be developed to enable a new generation of predictive fuel performance codes that are not limited by empiricism.
Laser-Generated Shocks and Bubbles as Laboratory-Scale Models of Underwater Explosions
Theodore G. Jones
2003-01-01
Full Text Available Underwater shocks and bubbles were generated using a high energy pulsed laser system. The advantages of this experimental approach are: (1 precisely controlled and measured experimental conditions; (2 improved diagnostics, including extensive imaging capabilities; (3 unique experiments, including a simultaneously detonated line charge; and (4 the ability to provide validation quality data for hydrodynamic simulation codes. Bubble sensitivity to variation of several experimental parameters was examined. Numerical simulations were performed corresponding to the experimental shots, showing that empirical bubble theory, experimental bubble data, and simulations were all in good agreement.
Atomistic modeling of growth and coalescence of helium nano-bubbles in tungsten
The mechanisms of growth and coalescence of helium nano-bubbles in tungsten are investigated using molecular dynamics simulations. It is shown that crystal symmetries and governed by them properties of dislocations, generated by the growing nano-bubbles, are responsible for main nano-bubble features revealed, including non-spherical shape and anisotropy of surrounding stress field. The transport of helium atoms in non-uniform stress field is simulated at different temperatures and the transport coefficients are determined. The implications of the considered dislocation and helium dynamics on nucleation and growth of bubbles in tungsten with implanted helium are discussed
Atomistic modeling of growth and coalescence of helium nano-bubbles in tungsten
Smirnov, R. D.; Krasheninnikov, S. I.; Guterl, J.
2015-08-01
The mechanisms of growth and coalescence of helium nano-bubbles in tungsten are investigated using molecular dynamics simulations. It is shown that crystal symmetries and governed by them properties of dislocations, generated by the growing nano-bubbles, are responsible for main nano-bubble features revealed, including non-spherical shape and anisotropy of surrounding stress field. The transport of helium atoms in non-uniform stress field is simulated at different temperatures and the transport coefficients are determined. The implications of the considered dislocation and helium dynamics on nucleation and growth of bubbles in tungsten with implanted helium are discussed.
Interfacial area concentration in gas–liquid bubbly to churn flow regimes in large diameter pipes
Highlights: • A systematic method to predict interfacial area concentration (IAC) is presented. • A correlation for group 1 bubble void fraction is proposed. • Correlations of IAC and bubble diameter are developed for group 1 bubbles. • Correlations of IAC and bubble diameter are developed for group 2 bubbles. • The newly-developed two-group IAC model compares well with collected databases. - Abstract: This study performed a survey on existing correlations for interfacial area concentration (IAC) prediction and collected an IAC experimental database of two-phase flows taken under various flow conditions in large diameter pipes. Although some of these existing correlations were developed by partly using the IAC databases taken in the low-void-fraction two-phase flows in large diameter pipes, no correlation can satisfactorily predict the IAC in the two-phase flows changing from bubbly, cap bubbly to churn flow in the collected database of large diameter pipes. So this study presented a systematic way to predict the IAC for the bubbly-to-churn flows in large diameter pipes by categorizing bubbles into two groups (group 1: spherical or distorted bubble, group 2: cap bubble). A correlation was developed to predict the group 1 void fraction by using the void fraction for all bubble. The group 1 bubble IAC and bubble diameter were modeled by using the key parameters such as group 1 void fraction and bubble Reynolds number based on the analysis of Hibiki and Ishii (2001, 2002) using one-dimensional bubble number density and interfacial area transport equations. The correlations of IAC and bubble diameter for group 2 cap bubbles were developed by taking into account the characteristics of the representative bubbles among the group 2 bubbles and the comparison between a newly-derived drift velocity correlation for large diameter pipes and the existing drift velocity correlation of Kataoka and Ishii (1987) for large diameter pipes. The predictions from the newly
Avdeev, Alexander A
2016-01-01
This monograph presents a systematic analysis of bubble system mathematics, using the mechanics of two-phase systems in non-equilibrium as the scope of analysis. The author introduces the thermodynamic foundations of bubble systems, ranging from the fundamental starting points to current research challenges. This book addresses a range of topics, including description methods of multi-phase systems, boundary and initial conditions as well as coupling requirements at the phase boundary. Moreover, it presents a detailed study of the basic problems of bubble dynamics in a liquid mass: growth (dynamically and thermally controlled), collapse, bubble pulsations, bubble rise and breakup. Special emphasis is placed on bubble dynamics in turbulent flows. The analysis results are used to write integral equations governing the rate of vapor generation (condensation) in non-equilibrium flows, thus creating a basis for solving a number of practical problems. This book is the first to present a comprehensive theory of boil...
Mathematical modeling of cold cap: Effect of bubbling on melting rate
Pokorny, Richard; Kruger, Albert A.; Hrma, Pavel R.
2014-12-31
The rate of melting is a primary concern in the vitrification of radioactive wastes because it directly influences the life cycle of nuclear waste cleanup efforts. To increase glass melting performance, experimental and industrial all-electric waste glass melters employ various melt-rate enhancement techniques, the most prominent being the application of bubblers submerged into molten glass. This study investigates various ways in which bubbling affects melting rate in a waste glass melter. Using the recently developed cold cap model, we suggest that forced convection of molten glass, which increases the cold cap bottom temperature, is the main factor. Other effects, such as stirring the feed into molten glass or reducing the insulating effect of foaming, also play a role.
Identifying bubble collapse in a hydrothermal system using hiddden Markov models
Dawson, Phillip B.; Benitez, M.C.; Lowenstern, Jacob B.; Chouet, Bernard A.
2012-01-01
Beginning in July 2003 and lasting through September 2003, the Norris Geyser Basin in Yellowstone National Park exhibited an unusual increase in ground temperature and hydrothermal activity. Using hidden Markov model theory, we identify over five million high-frequency (>15 Hz) seismic events observed at a temporary seismic station deployed in the basin in response to the increase in hydrothermal activity. The source of these seismic events is constrained to within ~100 m of the station, and produced ~3500–5500 events per hour with mean durations of ~0.35–0.45 s. The seismic event rate, air temperature, hydrologic temperatures, and surficial water flow of the geyser basin exhibited a marked diurnal pattern that was closely associated with solar thermal radiance. We interpret the source of the seismicity to be due to the collapse of small steam bubbles in the hydrothermal system, with the rate of collapse being controlled by surficial temperatures and daytime evaporation rates.
Nur-E- Mostafa
2016-01-01
Full Text Available This paper presents a numerical study with pressure-based finite volume method for prediction of non-cavitating and time dependent cavitating flow on hydrofoil. The phenomenon of cavitation is modeled through a mixture model. For the numerical simulation of cavitating flow, a bubble dynamics cavitation model is used to investigate the unsteady behavior of cavitating flow and describe the generation and evaporation of vapor phase. The non-cavitating study focuses on choosing mesh size and the influence of the turbulence model. Three turbulence models such as Spalart-Allmaras, Shear Stress Turbulence (SST k-ω model and Re-Normalization Group (RNG k-ε model with enhanced wall treatment are used to capture the turbulent boundary layer on the hydrofoil surface. The cavitating study presents an unsteady behavior of the partial cavity attached to the foil at different time steps for σ=0.8. Moreover, this study focuses on cavitation inception, the shape and general behavior of sheet cavitation, lift and drag forces for different cavitation numbers. Finally, the flow pattern and hydrodynamic characteristics are also studied at different angles of attack.
Numerical analyses of boiling two-phase bubbly flows in pipes with particle tracking method
The present study deals with the numerical modeling of two-phase bubbly flows with subcooled boiling in a flow channel using Particle Tracking Method. It is an on-going project whose final aim is to develop a computer code to more accurately predict, and better understand, the Departure from Nuclear Boiling (DNB) under PWR accidental conditions. At the first stage of the present study, this presentation reports the development of bubble dynamics models such as the bubble coalesces/break models and bubble heat transfer models. The validation of the models by comparing simulation results with experimental data in literature is also reported. (author)
Magnetic Bubble Expansion as an Experimental Model for Extra-Galactic Radio Lobes
Lynn, Alan; Zhang, Yue; Hsu, Scott
2010-11-01
The Plasma Bubble Expansion Experiment (PBEX) is conducting laboratory experiments to address outstanding nonlinear plasma physics issues related to how magnetic energy and helicity carried by extra-galactic jets interacts with the intergalactic medium to form radio lobe structures. Experiments are being conducted in the 4 meter long, 50 cm diameter HELCAT linear plasma device at UNM. A pulsed magnetized coaxial gun (˜10 kV, ˜100 kA, ˜2 mWb) forms and injects magnetized plasma bubbles perpendicularly into a lower pressure weakly magnetized background plasma formed by a helicon and/or hot cathode source in HELCAT. Ideal MHD simulations show that an MHD shock develops ahead of the bubble as it propagates, and that the bubble develops asymmetries due to the background field [1]. Experimental data from plasma bubble injection into a background plasma, particularly magnetic probe measurements, will be discussed. [4pt] [1] W. Liu et al., Phys. Plasmas 15, 072905 (2008).
Bubbles, shocks and elementary technical trading strategies
Fry, John
2014-01-01
In this paper we provide a unifying framework for a set of seemingly disparate models for bubbles, shocks and elementary technical trading strategies in financial markets. Markets operate by balancing intrinsic levels of risk and return. This seemingly simple observation is commonly over-looked by academics and practitioners alike. Our model shares its origins in statistical physics with others. However, under our approach, changes in market regime can be explicitly shown to represent a phase transition from random to deterministic behaviour in prices. This structure leads to an improved physical and econometric model. We develop models for bubbles, shocks and elementary technical trading strategies. The list of empirical applications is both interesting and topical and includes real-estate bubbles and the on-going Eurozone crisis. We close by comparing the results of our model with purely qualitative findings from the finance literature.
Bubble Dynamics and Shock Waves
2013-01-01
This volume of the Shock Wave Science and Technology Reference Library is concerned with the interplay between bubble dynamics and shock waves. It is divided into four parts containing twelve chapters written by eminent scientists. Topics discussed include shock wave emission by laser generated bubbles (W Lauterborn, A Vogel), pulsating bubbles near boundaries (DM Leppinen, QX Wang, JR Blake), interaction of shock waves with bubble clouds (CD Ohl, SW Ohl), shock propagation in polydispersed bubbly liquids by model equations (K Ando, T Colonius, CE Brennen. T Yano, T Kanagawa, M Watanabe, S Fujikawa) and by DNS (G Tryggvason, S Dabiri), shocks in cavitating flows (NA Adams, SJ Schmidt, CF Delale, GH Schnerr, S Pasinlioglu) together with applications involving encapsulated bubble dynamics in imaging (AA Doinikov, A Novell, JM Escoffre, A Bouakaz), shock wave lithotripsy (P Zhong), sterilization of ships’ ballast water (A Abe, H Mimura) and bubbly flow model of volcano eruptions ((VK Kedrinskii, K Takayama...
Lu Zhang
2015-07-01
Full Text Available The RNA polymerase II (Pol II is a eukaryotic enzyme that catalyzes the synthesis of the messenger RNA using a DNA template. Despite numerous biochemical and biophysical studies, it remains elusive whether the "secondary channel" is the only route for NTP to reach the active site of the enzyme or if the "main channel" could be an alternative. On this regard, crystallographic structures of Pol II have been extremely useful to understand the structural basis of transcription, however, the conformation of the unpaired non-template DNA part of the full transcription bubble (TB is still unknown. Since diffusion routes of the nucleoside triphosphate (NTP substrate through the main channel might overlap with the TB region, gaining structural information of the full TB is critical for a complete understanding of Pol II transcription process. In this study, we have built a structural model of Pol II with a complete transcription bubble based on multiple sources of existing structural data and used Molecular Dynamics (MD simulations together with structural analysis to shed light on NTP entry pathways. Interestingly, we found that although both channels have enough space to allow NTP loading, the percentage of MD conformations containing enough space for NTP loading through the secondary channel is twice higher than that of the main channel. Further energetic study based on MD simulations with NTP loaded in the channels has revealed that the diffusion of the NTP through the main channel is greatly disfavored by electrostatic repulsion between the NTP and the highly negatively charged backbones of nucleotides in the non-template DNA strand. Taken together, our results suggest that the secondary channel is the major route for NTP entry during Pol II transcription.
Lohse, Detlef
2003-01-01
With their ubiquitous occurrence in a multitude of fluid systems, bubbles occupy an important place in contemporary science and technology. One can readily cite several examples: the production and transport of oil, in which bubbles are purposely injected to help lift heavy oil to the surface; energ
Proposals for improving interphase drag modelling for the bubbly and slug regimes in RELAP5
The proposal is put forward that the effective interphase drag coefficient for the bubbly and slug regimes in RELAP5 should be calculated using best-estimate void fraction correlations. It is argued that this will lead to improvements in the code's modelling of interphase drag and evidence is given to corroborate this. The need for such improvements has been prompted by the poor performance of the current models in simulating rod bundle experiments. There is also concern that the models do not account for profile slip effects, which could be important in a variety of geometries, and that the slug flow equations may not be appropriate for large diameter vertical pipes. To support the proposal, a set of void fraction correlations is identified which is believed to cover the full range of geometries and flow conditions encountered in PWR safety analysis including the analysis of small-scale experimental facilities. This set is selected from a detailed appraisal of the most appropriate correlations found in the literature which takes account of comparisons with experimental data and physical considerations. This Report forms part of the UK's commitment to the ICAP Code Improvement Plan. The recommendations will now be implemented in a development version of RELAP5/MOD3 and a preliminary assessment made. The interphase drag models used in the annular-mist regime will also be examined and, if necessary, appropriate improvements will be proposed. (author)
De Backer, A., E-mail: andree.debacker@ccfe.ac.uk [UMET, UMR 8207, Université Lille 1, Villeneuve d’Ascq (France); CCFE, Culham Centre for Fusion Energy, Abingdon, Oxon (United Kingdom); Adjanor, G.; Domain, C.; Lescoat, M.L. [EDF R& D, MMC Centre des Renardières, Moret-sur-Loing (France); Jublot-Leclerc, S.; Fortuna, F.; Gentils, A. [CSNSM, Univ Paris-Sud, CNRS/IN2P3, Orsay (France); Ortiz, C.J. [CIEMAT, Laboratorio Nacional de Fusión por Confinamiento Magnético, Madrid (Spain); Souidi, A. [Université Dr. Tahar Moulay de Saida, Saida (Algeria); Becquart, C.S. [UMET, UMR 8207, Université Lille 1, Villeneuve d’Ascq (France)
2015-06-01
Implantation of 10 keV helium in 316L steel thin foils was performed in JANNuS-Orsay facility and modeled using a multiscale approach. Density Functional Theory (DFT) atomistic calculations [1] were used to obtain the properties of He and He-vacancy clusters, and the Binary Collision Approximation based code MARLOWE was applied to determine the damage and He-ion depth profiles as in [2,3]. The processes involved in the homogeneous He bubble nucleation and growth were defined and implemented in the Object Kinetic Monte Carlo code LAKIMOCA [4]. In particular as the He to dpa ratio was high, self-trapping of He clusters and the trap mutation of He-vacancy clusters had to be taken into account. With this multiscale approach, the formation of bubbles was modeled up to nanometer-scale size, where bubbles can be observed by Transmission Electron Microscopy. Their densities and sizes were studied as functions of fluence (up to 5 × 10{sup 19} He/m{sup 2}) at two temperatures (473 and 723 K) and for different sample thicknesses (25–250 nm). It appears that the damage is not only due to the collision cascades but is also strongly controlled by the He accumulation in pressurized bubbles. Comparison with experimental data is discussed and sensible agreement is achieved.
Louisnard, Olivier
2013-01-01
In a companion paper, a reduced model for propagation of acoustic waves in a cloud of inertial cavitation bubbles was proposed. The wave attenuation was calculated directly from the energy dissipated by a single bubble, the latter being estimated directly from the fully nonlinear radial dynamics. The use of this model in a mono-dimensional configuration has shown that the attenuation near the vibrating emitter was much higher than predictions obtained from linear theory, and that this strong attenuation creates a large traveling wave contribution, even for closed domain where standing waves are normally expected. In this paper, we show that, owing to the appearance of traveling waves, the primary Bjerknes force near the emitter becomes very large and tends to expel the bubbles up to a stagnation point. Two-dimensional axi-symmetric computations of the acoustic field created by a large area immersed sonotrode are also performed, and the paths of the bubbles in the resulting Bjerknes force field are sketched. C...
were compared with several experimental data sets of steam-water two-phase flow in large-diameter pipes. One-dimensional and axi-symmetric two-dimensional two-fluid models were used for this purpose. Since the present interfacial drag force model would be applicable only to bubbly two-phase flow, the experimental data in which the mean void fraction is less than 0.3 were used for this purpose. As a result, it was demonstrated that the new model is effective for accurate prediction in these experimental conditions. (author)
A One-Dimensional (1-D) Three-Region Model for a Bubbling Fluidized-Bed Adsorber
Lee, Andrew; Miller, David C.
2012-01-01
A general one-dimensional (1-D), three-region model for a bubbling fluidized-bed adsorber with internal heat exchangers has been developed. The model can predict the hydrodynamics of the bed and provides axial profiles for all temperatures, concentrations, and velocities. The model is computationally fast and flexible and allows for any system of adsorption and desorption reactions to be modeled, making the model applicable to any adsorption process. The model has been implemented in both gPROMS and Aspen Custom Modeler, and the behavior of the model has been verified.
Evaluation of two bubble-detachment models for two-phase flow
Two-phase flows with attached wall bubbles are considered as flows in roughened ducts with roughness elements equal in magnitude to the bubble detachment diameters determined by a force balance analysis. Two methods for determining the drag component of the forces acting on a bubble give distinctly different mixing-length theory expressions for velocity profile, friction factor, and void fraction. Available data are shown to be compatible with that set of expressions in which drag is due primarily to skin friction rather than form drag
Bubble, Bubble, Toil and Trouble.
Journal of Chemical Education, 2001
2001-01-01
Bubbles are a fun way to introduce the concepts of surface tension, intermolecular forces, and the use of surfactants. Presents two activities in which students add chemicals to liquid dishwashing detergent with water in order to create longer lasting bubbles. (ASK)
Dirk G Baur; Kristoffer Glover
2012-01-01
In this paper we use a test developed by Phillips et al. (2011) to identify a bubble in the gold market. We find that the price of gold followed an explosive price process between 2002 and 2012 interrupted only briefly by the subprime crisis in 2008. We also provide a theoretical foundation for such bubble tests based on a behavioural model of heterogeneous agents and demonstrate that periods of explosive price behaviour are consistent with increased chartist activity in the gold market. The ...
A Hadronic-Leptonic Model for the Fermi Bubbles: Cosmic-Rays in the Galactic Halo and Radio Emission
Fujita, Yutaka; Yamazaki, Ryo
2014-01-01
We investigate non-thermal emission from the Fermi bubbles on a hadronic model. Cosmic-ray (CR) protons are accelerated at the forward shock of the bubbles. They interact with the background gas in the Galactic halo and create $\\pi^0$-decay gamma-rays and secondary electrons through proton-proton interaction. We follow the evolution of the CR protons and electrons by calculating their distribution functions. We find that the spectrum and the intensity profile of $\\pi^0$-decay gamma-rays are consistent with observations. We predict that the shock front is located far ahead of the gamma-ray boundary of the Fermi bubbles. This naturally explains the fact that a clear temperature jump of thermal gas was not discovered at the gamma-ray boundary in recent Suzaku observations. We also consider re-acceleration of the background CRs in the Galactic halo at the shock front. We find that it can significantly affect the gamma-rays from the Fermi bubbles, unless the density of the background CRs is $\\lesssim 10$\\% of that...
Powering of cool filaments in cluster cores by buoyant bubbles - I. Qualitative model
Churazov, E.; Ruszkowski, M.; Schekochihin, A.
2013-11-01
Cool-core clusters (e.g. Perseus or M87) often possess a network of bright gaseous filaments, observed in radio, infrared, optical and X-ray bands. We propose that these filaments are powered by the reconnection of the magnetic field in the wakes of buoyant bubbles. Active galactic nucleus (AGN)-inflated bubbles of relativistic plasma rise buoyantly in the cluster atmosphere, stretching and amplifying the field in the wake to values of β = 8πPgas/B2 ˜ 1. The field lines in the wake have opposite directions and are forced together as the bubble motion stretches the filament. This setup bears strong similarity to the coronal loops on the Sun or to the Earth's magnetotail. The reconnection process naturally explains both the required level of local dissipation rate in filaments and the overall luminosity of filaments. The original source of power for the filaments is the potential energy of buoyant bubbles, inflated by the central AGN.
Powering of cool filaments in cluster cores by buoyant bubbles. I. Qualitative model
Churazov, E; Schekochihin, A
2013-01-01
Cool-core clusters (e.g., Perseus or M87) often possess a network of bright gaseous filaments, observed in radio, IR, optical and X-ray bands. We propose that these filaments are powered by the reconnection of the magnetic field in the wakes of buoyant bubbles. AGN-inflated bubbles of relativistic plasma rise buoyantly in the cluster atmosphere, stretching and amplifying the field in the wake to values of $\\beta =8\\pi P_{gas}/B^2\\sim 1$. The field lines in the wake have opposite directions and are forced together as the bubble motion stretches the filament. This setup bears strong similarity to the coronal loops on the Sun or the Earth magneto-tail. The reconnection process naturally explains both the required level of local dissipation rate in filaments and the overall luminosity of filaments. The original source of power for the filaments is the potential energy of buoyant bubbles, inflated by the central AGN.
Modeling and simulation of the bubble-induced flow in wine fermentation vessels
Schmidt Dominik
2015-01-01
Full Text Available Detailed flow pattern analyses regarding wine fermentations conducted without mechanical agitation are limited to lab-scale investigations, as industrial size measurements are expensive and difficult to realize. Computational fluid dynamic (CFD methods can offer an alternative and more flexible approach to gain insight into such bubble induced fluid flows. Therefore, the aim of this study was to transfer the findings of existing research onto a CFD model capable of capturing the three- dimensional flow pattern in industrial scale wine fermentation vessels. First results were obtained by using an extended version of the OpenFOAM® (v.2.2.x solver multiphaseEulerFoam for modeling the gas-liquid two phase system. With parameters from the most vigorous phase of wine fermentation a fully developed, unsteady flow regime could be established after approx. 120 s of real time. Thereby the groundwork for further evaluations of e.g. mixing efficiency or cooling equipment optimizations with CFD methods is laid.
A gas bubble-based parallel micro manipulator : conceptual design and kinematics model.
Dong, Wei; Gauthier, Michaël; Lenders, Cyrille; Lambert, Pierre
2012-01-01
The parallel mechanism has become an alternative solution when micro manipulators are demanded in the fields of micro manipulation and micro assembly. In this paper, a three-Degree-Of-Freedom (3-DOF) parallel micro manipulator is presented, which is directly driven by three micro gas bubbles. Since the micro gas bubbles are generated and maintained due to the surface tension between the gas and liquid media, the proposed novel system can be used in the liquid environment which allows for rota...
Characteristics of a bubble jet near a vertical wall
ZHANG A-man; YAO Xiong-liang; LI Jia
2008-01-01
A numerical model of a coupled bubble jet and wall was built on the assumption of potential flow and calculated by the boundary integral method. A three-dimensional computing program was then developed. Starting with the basic phenomenon of the interaction between a bubble and a wall, the dynamics of bubbles near rigid walls were studied systematically with the program. Calculated results agreed well with experimental results. The relationship between the Bjerknes effect of a wall and characteristic parameters was then studied and the calculated results of various cases were compared and discussed with the Blake criterion based on the Kelvin-impulse theory. Our analyses show that the angle of the jet's direction and the pressure on the rigid wall have a close relationship with collapse force and the bubble's characteristic parameters. From this, the application range of Blake criterion can be determined. This paper aims to provide a basis for future research on the dynamics of bubbles near a wall.
A dry-spot model for the prediction of critical heat flux in water boiling in bubbly flow regime
Ha, Sang Jun; No, Hee Cheon [Korea Advanced Institute of Science and Technology, Taejon (Korea, Republic of)
1997-12-31
This paper presents a prediction of critical heat flux (CHF) in bubbly flow regime using dry-spot model proposed recently by authors for pool and flow boiling CHF and existing correlations for forced convective heat transfer coefficient, active site density and bubble departure diameter in nucleate boiling region. Without any empirical constants always present in earlier models, comparisons of the model predictions with experimental data for upward flow of water in vertical, uniformly-heated round tubes are performed and show a good agreement. The parametric trends of CHF have been explored with respect to variations in pressure, tube diameter and length, mass flux and inlet subcooling. 16 refs., 6 figs., 1 tab. (Author)
The aim of this thesis is the 3-D modeling and numerical simulation of liquid/gas (water/vapor or water/air) two-phase flows in cooling circuits of nuclear power plants during normal and accidental situations. The development of a multidimensional dual-fluid model encounters two problems: the statistical effects of turbulence and the interface mass, momentum and energy transfers. The models developed in this study were introduced in the 3-D module of the CATHARE code developed by the CEA and the results were compared to experimental results available in the literature. The first chapter describes the equations of the local dual-fluid model for the 3-D description of two-phase flows. Closing relations adapted to dispersed flows with isothermal bubbles and without phase transformation are proposed and focus on the momentum transfer at the interfaces. The theoretical study of turbulence in the liquid phase of a bubble flow is modelled in chapter 2. Chapter 3 deals with the voluminal interface area used in the interface mass, momentum and energy transfers, and chapters 4 and 5 concern the application of the developed models to concrete situations. Chapter 4 describes in details the 3-D module of the CATHARE code while chapter 5 gives a comparison of numerical results obtained using the CATHARE code with other experimental results obtained at EdF. (J.S.)
Donna Post Guillen; Tami Grimmett; Anastasia M. Gribik; Steven P. Antal
2010-09-01
The Hybrid Energy Systems Testing (HYTEST) Laboratory is being established at the Idaho National Laboratory to develop and test hybrid energy systems with the principal objective to safeguard U.S. Energy Security by reducing dependence on foreign petroleum. A central component of the HYTEST is the slurry bubble column reactor (SBCR) in which the gas-to-liquid reactions will be performed to synthesize transportation fuels using the Fischer Tropsch (FT) process. SBCRs are cylindrical vessels in which gaseous reactants (for example, synthesis gas or syngas) is sparged into a slurry of liquid reaction products and finely dispersed catalyst particles. The catalyst particles are suspended in the slurry by the rising gas bubbles and serve to promote the chemical reaction that converts syngas to a spectrum of longer chain hydrocarbon products, which can be upgraded to gasoline, diesel or jet fuel. These SBCRs operate in the churn-turbulent flow regime which is characterized by complex hydrodynamics, coupled with reacting flow chemistry and heat transfer, that effect reactor performance. The purpose of this work is to develop a computational multiphase fluid dynamic (CMFD) model to aid in understanding the physico-chemical processes occurring in the SBCR. Our team is developing a robust methodology to couple reaction kinetics and mass transfer into a four-field model (consisting of the bulk liquid, small bubbles, large bubbles and solid catalyst particles) that includes twelve species: (1) CO reactant, (2) H2 reactant, (3) hydrocarbon product, and (4) H2O product in small bubbles, large bubbles, and the bulk fluid. Properties of the hydrocarbon product were specified by vapor liquid equilibrium calculations. The absorption and kinetic models, specifically changes in species concentrations, have been incorporated into the mass continuity equation. The reaction rate is determined based on the macrokinetic model for a cobalt catalyst developed by Yates and Satterfield [1]. The
Veshchunov, M.S. [Russian Academy of Sciences, Nuclear Safety Institute (IBRAE), 52, B. Tulskaya, Moscow 115191 (Russian Federation)]. E-mail: vms@ibrae.ac.ru
2005-11-15
A new mechanism of the lenticular grain face bubble migration which controls the bubble mobility and determines the drag force exerted on the grain boundary, is developed. It is shown that besides a more complicated (so called 'lenticular') shape of grain face bubbles, the migration mechanism of these bubbles might be essentially different from the intragranular bubbles, owing to their specific location on and interaction with a grain boundary. The model is validated against tests on grain growth kinetics during steady irradiation exposure and during post-irradiation annealing of UO{sub 2} fuel samples, and allows explanation of a strong retarding effect of irradiation on the grain growth observed in these tests.
To predict the behavior of gas-liquid two-phase flows in a centrifugal pump impeller, a three-dimensional numerical method is proposed on the basis of a bubbly flow model. Under the assumption of homogeneous bubbly flow entraining fine bubbles, the equation of motion of the mixture is represented by that of liquid-phase and the liquid velocity is expressed as a potential for a quasi-harmonic equation. This equation is solved with a finite element method to obtain the velocities, and the equation of motion of an air bubble is integrated numerically in the flow field to obtain the void fraction. These calculations are iterated to obtain a converged solution. The method has been applied to a radial-flow pump, and the results obtained have been confirmed by experiments within the range of bubbly flow regime
O'Geary, Melissa A.
Bubbles provide an enjoyable and festive medium through which to teach many concepts within the science topics of light, color, chemistry, force, air pressure, electricity, buoyancy, floating, density, among many others. In order to determine the nature of children's engagement within a museum setting and the learning opportunities of playing with bubbles, I went to a children's interactive museum located in a metropolitan city in the Northeastern part of the United States.
Modeling and simulation of bubbling hot well deaerator in condensers of ships
Highlights: • The model was based on the specific structure and heat transfer process. • The solubility of Oxygen could be calculated. • The simulation results are consistent with the theoretical analysis. - Abstract: Application of deaerators in ship has been restricted due to space limitation. Thus far, the function of deaerator has been integrated into the hot well of the condenser. The condensate water would be heated to saturation by extracting steam from turbine, which would make the solubility of Oxygen in condensate water fall to zero according to Henry’s law and Dalton’s law. The purpose of this paper was to build a mathematical model of bubbling hot well deaerator in the condensers of ships. In this paper, the heat exchange rate was calculated by empirical formulas which took the specific structure and process of heat exchange into account. When the operating conditions were in the application ranges of the empirical formulas, the simulation model would be performed by utilizing them; otherwise, calculations would be done by the conservation of energy, which assured the simulation model could be used at any operating condition. Different from previous works, the solubility of Oxygen in heated condensate water could be calculated by an empirical formula. The simulation results showed that the structure and heat exchange process considered could be highly accurate at the steady-state operations, and the main parameters trend curves during dynamic-state operations were consistent with theoretical analysis. The solubility of Oxygen could be calculated and the simulation results at the steady-state operations were verified against the practical situation, the trend curves during dynamic-state operations were consistent with theoretical analysis
A non-equilibrium ionization model of the Local Bubble (I)
de Avillez, Miguel A
2012-01-01
Aims. We present the first high-resolution non-equilibrium ionization simulation of the joint evolution of the Local Bubble (LB) and Loop I superbubbles in the turbulent supernova-driven interstellar medium (ISM). The time variation and spatial distribution of the Li-like ions Civ, Nv, and Ovi inside the LB are studied in detail. Methods. This work uses the parallel adaptive mesh refinement code EAF-PAMR coupled to the newly developed atomic and molecular plasma emission module E(A+M)PEC, featuring the time-dependent calculation of the ionization structure of H through Fe, using the latest revision of solar abundances. The finest AMR resolution is 1 pc within a grid that covers a representative patch of the Galactic disk (with an area of 1 kpc^2 in the midplane) and halo (extending up to 10 kpc above and below the midplane). Results. The evolution age of the LB is derived by the match between the simulated and observed absorption features of the Li-like ions Civ, Nv, and Ovi . The modeled LB current evolution...
Issues for Achieving an Experimental Model Concerning Bubble Deck Concrete Slab with Spherical Gaps
Sergiu Călin
2010-01-01
Full Text Available After realizing numerous constructions in the world, which use Bubble Deck concrete slabs with spherical gaps, valuable information were gathered, allowing a rigorous processing and systematization, with the purpose of realizing an experimental and documentary study. The paper presents some experimental programs which refer to concrete slabs with spherical gaps, existing in similar execution and loading conditions as those from a real construction; this implies the realization of a monolithic slab element at a scale of 1:1, which will be subjected to static gravitational loadings in order to determine the deformation (deflection, cracking and failing characteristics. The resultant conclusions will be used in defining the failing mechanisms, very useful in the formulation of an adequate mathematical model. The research proposed in the project offers an answer to the major objectives of the development of calculus methods and existent prescriptions of the concrete slabs with spherical gaps. The realization of the proposed objectives involves documentation activities, theoretical study, collaboration with different other partners, gathering and processing of the results obtained in laboratory and even in situ.
Identifying bubble collapse in a hydrothermal system using hidden Markov models
Dawson, P.B.; Benitez, M.C.; Lowenstern, J. B.; Chouet, B.A.
2012-01-01
Beginning in July 2003 and lasting through September 2003, the Norris Geyser Basin in Yellowstone National Park exhibited an unusual increase in ground temperature and hydrothermal activity. Using hidden Markov model theory, we identify over five million high-frequency (>15Hz) seismic events observed at a temporary seismic station deployed in the basin in response to the increase in hydrothermal activity. The source of these seismic events is constrained to within ???100 m of the station, and produced ???3500-5500 events per hour with mean durations of ???0.35-0.45s. The seismic event rate, air temperature, hydrologic temperatures, and surficial water flow of the geyser basin exhibited a marked diurnal pattern that was closely associated with solar thermal radiance. We interpret the source of the seismicity to be due to the collapse of small steam bubbles in the hydrothermal system, with the rate of collapse being controlled by surficial temperatures and daytime evaporation rates. copyright 2012 by the American Geophysical Union.
Bubble chambers may have almost vanished from the front line of physics research, but the vivid memory of their intricate and sometimes beautiful patterns of particle tracks lives on, and has greatly influenced the computer graphics of track reconstruction in today's big experiments. 'Seeing' an interaction makes it more understandable. Bubble chambers, with their big collaborations of physicists from many widely scattered research institutes, started another ball rolling. The groups formed are even now only surpassed in size by the big collaborations working on today's major detectors at colliding beam machines. From 14-16 July, about 130 physicists gathered at CERN to commemorate the 40th anniversary of the invention of the bubble chamber by Donald Glaser. The meeting, organized by Derek C. Colley from Birmingham, gave a comprehensive overview of bubble chamber contributions to physics, their challenging technology, and the usefulness of bubble chamber photographs in education, both for physics and the public at large. After opening remarks by CERN Director Carlo Rubbia, Donald Glaser began with a brief review of the work which led to his invention - there was much more to it than idly watching beer bubbles rise up the wall of the glass - before turning to his present line of research, biophysics, also very visually oriented
Movahed, Pooya; Kreider, Wayne; Maxwell, Adam D; Hutchens, Shelby B; Freund, Jonathan B
2016-08-01
A generalized Rayleigh-Plesset-type bubble dynamics model with a damage mechanism is developed for cavitation and damage of soft materials by focused ultrasound bursts. This study is linked to recent experimental observations in tissue-mimicking polyacrylamide and agar gel phantoms subjected to bursts of a kind being considered specifically for lithotripsy. These show bubble activation at multiple sites during the initial pulses. More cavities appear continuously through the course of the observations, similar to what is deduced in pig kidney tissues in shock-wave lithotripsy. Two different material models are used to represent the distinct properties of the two gel materials. The polyacrylamide gel is represented with a neo-Hookean elastic model and damaged based upon a maximum-strain criterion; the agar gel is represented with a strain-hardening Fung model and damaged according to the strain-energy-based Griffith's fracture criterion. Estimates based upon independently determined elasticity and viscosity of the two gel materials suggest that bubble confinement should be sufficient to prevent damage in the gels, and presumably injury in some tissues. Damage accumulation is therefore proposed to occur via a material fatigue, which is shown to be consistent with observed delays in widespread cavitation activity. PMID:27586763
Bubble point pressures of the selected model system for CatLiq® bio-oil process
Toor, Saqib Sohail; Rosendahl, Lasse; Baig, Muhammad Noman; Robbins, Phil; Santos, Regina; Nørgaard, Kristian
2010-01-01
work, the bubble point pressures of a selected model mixture (CO2 + H2O + Ethanol + Acetic acid + Octanoic acid) were measured to investigate the phase boundaries of the CatLiq® process. The bubble points were measured in the JEFRI-DBR high pressure PVT phase behavior system. The experimental results......The CatLiq® process is a second generation catalytic liquefaction process for the production of bio-oil from WDGS (Wet Distillers Grains with Solubles) at subcritical conditions (280-350 oC and 225-250 bar) in the presence of a homogeneous alkaline and a heterogeneous Zirconia catalyst. In this...... predicted data is 8.7% in the selected model mixture....
Investigation of two-phase flow around obstacles is important for thermal-hydraulic design and safety analyses of various equipments in nuclear power plants. The examples of obstacles are spacers in nuclear fuel rod bundles or in steam generator tube bundles, as well as perforated plates submerged in two-phase mixture in order to equalize void fraction and swell level spatial distribution in large volumes. In this paper the numerical simulation of a bubble column flow in a rectangular water pool is performed. The gas is injected in form of bubbles near the pool bottom through the short perforated pipe. A horizontal plate with a slit is introduced at the half-length of the pool and it forms an obstacle to the bubble column flow. The simulation is performed with the multidimensional transient two-fluid model. The model is based on the mass and momentum balance equations for the liquid and gas phase. The model is numerically solved with the application of the SIMPLE type algorithm developed for the conditions of two-phase flow. Obtained results show the ability of the applied modelling approach to simulate complex liquid and gas flow conditions around obstacles. The bubble column rise along the pool and passing through the plate slit is shown, together with the recirculation of water below and above the horizontal plate. The problem of achieving symmetric velocity and void fraction fields for symmetric rectangular pool geometry and boundary conditions is encountered. The simulations without and with the lateral force are performed, showing the importance of lateral force inclusion in the momentum balance equations. (author)
A Phase-Field Approach to Modeling Hydrate Formation on Methane Gas Bubbles in a Water Column
Fu, X.; Cueto-Felgueroso, L.; Waite, W. F.; Ruppel, C. D.; Juanes, R.
2014-12-01
Methane hydrates are water-based crystalline solids, where gas molecules are trapped inside the lattice structure formed by water. Most commonly found in deep ocean floors where low temperature and high pressure are primal conditions for hydrate to form, gas hydrates contain most of the world's mobile carbon and yet it remains an important and open question how methane leakage from gas hydrate impacts ocean and the atmosphere. While current work focus on the breakdown of gas hydrate in marine environment and the the release of methane from seafloor, few studies explore the fate of a single or a plume of methane bubbles when entering the water column after the release. We propose to study the fate of an individual and a series of methane bubbles through mathematical modeling, specifically using a phase-field approach. Phase-field modeling is a mathematical framework that describes systems that are out of thermodynamic equilibrium. First introduced in the context of solidification process and phase transitions, it has since been adopted in the field of multiphase flow. In this work, we present a new phase-field formulation for multiphase/multicomponent flows that allows us to model the fate of methane bubbles in the water system as a nonequilibrium process.
Mechanisms of gas bubble retention
Retention and episodic release of flammable gases are critical safety concerns regarding double-shell tanks (DSTs) containing waste slurries. Previous investigations have concluded that gas bubbles are retained by the slurry that has settled at the bottom of the DST. However, the mechanisms responsible for the retention of these bubbles are not well understood. In addition, the presence of retained gas bubbles is expected to affect the physical properties of the sludge, but essentially no literature data are available to assess the effect of these bubbles. The rheological behavior of the waste, particularly of the settled sludge, is critical to characterizing the tendency of the waste to retain gas bubbles. The objectives of this study are to elucidate the mechanisms contributing to gas bubble retention and release from sludge such as is in Tank 241-SY-101, understand how the bubbles affect the physical properties of the sludge, develop correlations of these physical properties to include in computer models, and collect experimental data on the physical properties of simulated sludges with bubbles. This report presents a theory and experimental observations of bubble retention in simulated sludge and gives correlations and new data on the effect of gas bubbles on sludge yield strength
Rational Bubbles in Stock Prices?
Behzad T. Diba; Grossman, Herschel I.
1985-01-01
This paper reports empirical tests for the existence of rational bubbles in stock prices. The analysis focuses on a familiar model that defines market fundamentals to be the expected present value of dividends, discounted at a constantrate, and defines a rational bubble to be a self-confirming divergence of stock prices from market fundamentals in response to extraneous variables. The tests are based on the theoretical result that, if rational bubbles exist, time series obtained by differenci...
Rational Asset Pricing Bubbles Revisited
Jan Werner
2012-01-01
Price bubble arises when the price of an asset exceeds the asset's fundamental value, that is, the present value of future dividend payments. The important result of Santos and Woodford (1997) says that price bubbles cannot exist in equilibrium in the standard dynamic asset pricing model with rational agents as long as assets are in strictly positive supply and the present value of total future resources is finite. This paper explores the possibility of asset price bubbles when either one of ...
A criterion for the occurrence of bubble fission and its modeling
Delale, C. F.; Pasinlioğlu, S.; Zima, Patrik
Singapore: Research Publishing Services, 2012 - (Ohl, C.; Klaseboer, E.; Ohl, S.; Gong, S.; Khoo, B.), s. 78-82 ISBN 978-981-07-2826-7. [International Symposium on Cavitation /8./ CAV 2012. Singapur (SG), 13.08.2012-16.08.2012] R&D Projects: GA ČR GAP101/10/1428 Institutional research plan: CEZ:AV0Z20760514 Keywords : bubble dynamics * bubble fission * cavitation Subject RIV: BK - Fluid Dynamics http://rpsonline.com.sg/proceedings/9789810728267/html/192.xml59.xml
In the channels of liquid metal (ML) regenerating sheaths of a fusion reactor, the possibility of bubbles of helium is not remote. Bubbles adhering to the wall of the ML channels would affect heat transfer and the permeation of tritium. Detailed analysis has been conducted (fine mesh), using OpenFOAM, from the environment of a bubble attached to the wall and has developed a model for permeation of tritium through a partially covered with bubbles of helium surface of contact Ml-solid. The model developed has implemented as wall function in OpenFOAM, has validated and has been applied to a case study, using a relatively thick mesh. The developed model substantially reduces the need for computing on the detailed calculation power.
Solares, H A Ayala; Hüntemeyer, P
2015-01-01
The Fermi Bubbles, which comprise two large and homogeneous regions of spectrally hard gamma-ray emission extending up to $55^{o}$ above and below the Galactic Center, were first noticed in GeV gamma-ray data from the Fermi Telescope in 2010. The mechanism or mechanisms which produce the observed hard spectrum are not understood. Although both hadronic and lep- tonic models can describe the spectrum of the bubbles, the leptonic model can also explain similar structures observed in microwave data from the WMAP and Planck satellites. Recent publications show that the spectrum of the Fermi Bubbles is well described by a power law with an exponential cutoff in the energy range of 100MeV to 500GeV. Observing the Fermi Bubbles at higher gamma-ray energies will help constrain the origin of the bubbles. A steeper cutoff will favor a leptonic model. The High Altitude Water Cherenkov (HAWC) Observatory, located 4100m above sea level in Mexico, is designed to measure high-energy gamma rays between 100GeV to 100TeV. With...
Muñoz-Cobo, José; Chiva, Sergio; El Aziz Essa, Mohamed; Mendes, Santos
2012-08-01
Two phase flow experiments with different superficial velocities of gas and water were performed in a vertical upward isothermal cocurrent air-water flow column with conditions ranging from bubbly flow, with very low void fraction, to transition flow with some cap and slug bubbles and void fractions around 25%. The superficial velocities of the liquid and the gas phases were varied from 0.5 to 3 m/s and from 0 to 0.6 m/s, respectively. Also to check the effect of changing the surface tension on the previous experiments small amounts of 1-butanol were added to the water. These amounts range from 9 to 75 ppm and change the surface tension. This study is interesting because in real cases the surface tension of the water diminishes with temperature, and with this kind of experiments we can study indirectly the effect of changing the temperature on the void fraction distribution. The following axial and radial distributions were measured in all these experiments: void fraction, interfacial area concentration, interfacial velocity, Sauter mean diameter and turbulence intensity. The range of values of the gas superficial velocities in these experiments covered the range from bubbly flow to the transition to cap/slug flow. Also with transition flow conditions we distinguish two groups of bubbles in the experiments, the small spherical bubbles and the cap/slug bubbles. Special interest was devoted to the transition region from bubbly to cap/slug flow; the goal was to understand the physical phenomena that take place during this transition A set of numerical simulations of some of these experiments for bubbly flow conditions has been performed by coupling a Lagrangian code, that tracks the three dimensional motion of the individual bubbles in cylindrical coordinates inside the field of the carrier liquid, to an Eulerian model that computes the magnitudes of continuous phase and to a 3D random walk model that takes on account the fluctuation in the velocity field of the
Highlights: • We simulate convective nucleate pool boiling with a novel phase-change model. • We simulate four cases at different sub-cooling and wall superheat levels. • We investigate the flow structures around the growing bubble and analyze the accompanying physics. • We accurately simulate bubble shape elongation and enhanced wall cooling due to the sliding and slanting motions of bubbles. • Bubble cycle durations show good agreement with experimental observations. - Abstract: With the long-term objective of Critical Heat Flux (CHF) prediction, bubble dynamics in convective nucleate boiling flows has been studied using a Direct Numerical Simulation (DNS). A sharp-interface phase change model which was originally developed for pool boiling flows is extended to convective boiling flows. For physical scales smaller than the smallest flow scales (smaller than the grid size), a micro-scale model was used. After a grid dependency study and a parametric study for the contact angle, four cases of simulation were carried out with different wall superheat and degree of subcooling. The flow structures around the growing bubble were investigated together with the accompanying physics. The relation between the heat flux evolution and the bubble growth was studied, along with investigations of bubble diameter and bubble base diameter evolutions across the four cases. As a validation, the evolutions of bubble diameter and bubble base diameter were compared to experimental observations. The bubble departure period and the bubble shapes show good agreement between the experiment and the simulation, although the Reynolds number of the simulation cases is relatively low
无
2011-01-01
The shift of China’s monetary policy stance from "moderately loose" to "prudent" in 2011 indicates curbing inflation and asset bubbles have become the Central Government’s top priority. But is China’s bubble problem short-term or long-term? Is it only monetary or related to economic structure? Is it the cause of China’s economic imbalance or the result? And what kind of deep-rooted problems in the macro economy does it reflect? All these questions call for deep thought,said Zhang Monan,a
CFD Model for Pneumatic Mixing with Bubble Chains: Application to Glass Melts
Šimčík, Miroslav; Růžička, Marek
2015-01-01
Roč. 127, MAY 4 (2015), s. 344-361. ISSN 0009-2509 R&D Projects: GA MŠk(CZ) LD13018 Institutional support: RVO:67985858 Keywords : bubble chain * viscous liquid * flow simulation Subject RIV: CI - Industrial Chemistry, Chemical Engineering Impact factor: 2.337, year: 2014