Computational study on cell structure evolution of random liquid and metal-melt foams

A method for geometrical and topological modeling the evolution of close-cell metallic foams based on the Voronoi tessellation in three-dimensional space is presented. Numerical computations were carried out to examine the evolution of the bubble size distribution and topological and geometric properties of aluminum foams in the liquid state, which were implemented by using McPherson’s new theory on coarsening of microstructures as well as the topological transition rules (T1 and T2 processes) in 3D foams, accounting for remarkable effects of both the gas diffusion and surface tension. Computational results show that the bubble size distributions of metallic foams are strongly coupled to the evolution of the cellular structure and dependent on the gas diffusivity and surface tension. The way of foam coarsening can be expressed as RR ₃₂=−mt ²+1 approximately, and gas diffusion between bubbles dominates the evolution of bubble sizes and foam structures. It is found that the average number of faces per bubble is 〈f〉=13.8, which is in good agreement with the values reported in the literature.     

Physicochemical approach to the theory of foam drainage

We have investigated theoretically the effect of surface viscoelasticity on the drainage of an aqueous foam. Former theories consider that the flow in Plateau borders is either Poiseuille flow or plug-flow. In the last case, the dissipation is attributed to flow in the nodes connecting Plateau borders. Although we do not include this dissipation in our model, we obtain a drainage equation that includes terms equivalent to those of the earlier models. We show that when the water solubility of the surfactant stabilizing the foam is low, the control parameter M for the transition between the two regimes is the ratio , where μ is the bulk viscosity, D _s the surface diffusion coefficient, R the radius of curvature of the Plateau border and ɛ the surface elasticity. When the surfactant is more soluble M is rather related to the bulk diffusion coefficient. Within the frame of this approach and in view of the estimated M values, we show that the flow in Plateau borders is Poiseuille-like. Received 26 June 2001     

Viscosity effects in foam drainage: Newtonian and non-newtonian foaming fluids

We have studied the drainage of foams made from Newtonian and non-Newtonian solutions of different viscosities. Forced-drainage experiments first show that the behavior of Newtonian solutions and of shear-thinning ones (foaming solutions containing either Carbopol or Xanthan) are identical, provided one considers the actual viscosity corresponding to the shear rate found inside the foam. Second, for these fluids, a drainage regime transition occurs as the bulk viscosity is increased, illustrating a coupling between surface and bulk flow in the channels between bubbles. The properties of this transition appear different from the ones observed in previous works in which the interfacial viscoelasticity was varied. Finally, we show that foams made of solutions containing long flexible PolyEthylene Oxide (PEO) molecules counter-intuitively drain faster than foams made with Newtonian solutions of the same viscosity. Complementary experiments made with fluids having all the same viscosity but different responses to elongational stresses (PEO-based Boger fluids) suggest an important role of the elastic properties of the PEO solutions on the faster drainage.     

Structural properties of silver nanoparticle agglomerates based on transmission electron microscopy: relationship to particle mobility analysis

In this work, the structural properties of silver nanoparticle agglomerates generated using condensation and evaporation method in an electric tube furnace followed by a coagulation process are analyzed using Transmission Electron Microscopy (TEM). Agglomerates with mobility diameters of 80, 120, and 150 nm are sampled using the electrostatic method and then imaged by TEM. The primary particle diameter of silver agglomerates was 13.8 nm with a standard deviation of 2.5 nm. We obtained the relationship between the projected area equivalent diameter (d _pa) and the mobility diameter (d _m), i.e., d _pa = 0.92 ± 0.03 d _m for particles from 80 to 150 nm. We obtained fractal dimensions of silver agglomerates using three different methods: (1) D _f = 1.84 ± 0.03, 1.75 ± 0.06, and 1.74 ± 0.03 for d _m = 80, 120, and 150 nm, respectively from projected TEM images using a box counting algorithm; (2) fractal dimension (D _fL) = 1.47 based on maximum projected length from projected TEM images using an empirical equation proposed by Koylu et al. (1995) Combust Flame 100:621–633; and (3) mass fractal-like dimension (D _fm) = 1.71 theoretically derived from the mobility analysis proposed by Lall and Friedlander (2006) J Aerosol Sci 37:260–271. We also compared the number of primary particles in agglomerate and found that the number of primary particles obtained from the projected surface area using an empirical equation proposed by Koylu et al. (1995) Combust Flame 100:621–633 is larger than that from using the relationship, d _pa = 0.92 ± 0.03 d _m or from using the mobility analysis.     

Annihilation rates of ortho—positronium in organic liquids: correlation with liquid compressibility and polarization ratio

The bubble model of ortho—positronium annihilation rates customarily employs the value of the planar-surface tension in the estimation of the (surface) energy required to expand a normal interstitial cavity to form the bubble. However, bubble radii are estimated to be only about 0.4 nm; the surface of such a small bubble is expected to have a tension different from that of the planar surface. We have replaced the surface tension parameter by the isothermal compressibility, which is thought not to have a nanostructure limitation. The model has been further improved by introducing the polarization ratio (n ² _D ? 1)/(n ² _D + 2) as a parameter, n _D being the refractive index of the organic liquid. Calculations are presented for 20 liquids.     

The foam/emulsion analogy in structure and drainage

The often quoted analogy between foams and emulsions is experimentally tested by studying properties after settling and under forced drainage of oil-in-water emulsions of drop size similar as for bubbles generally used in foam experiments. Observations with regard to structure, water fraction and drainage wave properties confirm the expected similarity in the low flow rate range. However, while for foams a convective circulation on the scale of the container sets in for values of water fraction exceeding about 0.2, no such convection is found in emulsions. Here instabilities are only encountered at water fractions of about 0.4, close to the void fraction of random packings of spheres. These take on the form of descending pulses of increased water fraction and lead to the transition from a frozen to a locally agitated structure.Received: 12 December 2003, Published online: 24 August 2004PACS: 82.70.-y Disperse systems; complex fluids - 47.20.-k Hydrodynamic stability - 47.55.Dz Drops and bubbles     

A physical model for dimensional reduction and its effects on the observable parameters of the universe

We construct a physical model to study the effects of dimensional reduction that might have taken place during the inflationary phase of the universe. The model we propose is a (1 + D)-dimensional (D > 3), nonsingular, spatially homogeneous and isotropic Friedmann model. We consider dimensional reduction to take place in a stepwise manner and interpret each step as a phase transition. Independent of the details of the process of dimensional reduction, we impose suitable boundary conditions across the transitions and trace the effects of dimensional reduction to the currently observable parameters of the universe. In order to exhibit the cosmological features of the proposed model, we construct a (1 + 4)-dimensional toy model for both closed and open cases of Friedmann geometries. It is shown that in these models the universe makes transition into the lower dimension when the critical length parameter l _4,3, which signals dimensional reduction, reaches the Planck length in D = 3. The numerical models we present in this paper have the capability of making definite predictions about the cosmological parameters of the universe such as the Hubble parameter, age and density.     

一维液态泡沫渗流实验研究及表面能和粘性耗散分析

液态泡沫由大量气泡密集堆积在微量表面活性剂溶液中形成，是远离平衡态的软物质. 泡沫强制渗流在微观上是指以恒定流率输入的液体在气泡间隙内的微流动过程，是影响泡沫稳定的主要因素之一. 采用在表面活性剂溶液中添加微量色素以显示泡沫中液体流动的方法，确定了透射率与液体分率的对应关系，测量得到了一维液态泡沫强制渗流中渗流波传播规律以及液体分率的演变规律；理论推导了泡沫基本单元，即开尔文单元结构(Kelvin cell)的粘性耗散能表达式，并依据Surface Evolver软件计算得到了不同液体分率时开尔文单元结构对应的的表面能，并计算出了与实验系统对应的开尔文单元结构的表面能和粘性耗散. 基于开尔文单元结构内液体分率演变的准静态假设，分析了表面能和粘性耗散的演变规律.     

Acoustic emission measurement of foam evolution in H2O–C2H5OH–air systems with content of detergent triton X-100

The acoustic emission (AE) of short-lived static foams from non-ionic detergent Triton X-100 in ethanol–water solutions has been investigated. The AE of the foam drainage process has shown two different stages of behaviour. The first step (up to 80 s) was without the appearance of AE signals. During this stage, bubble rearrangement has been observed. The next, main stage of drainage process proceeded with a high intensity AE generation. The results suggest that the AE technique gives the direct possibility of foam stability measurement.     

Quasi-one-dimensional foam drainage

Foam drainage is considered in a froth flotation cell. Air flow through the foam is described by a simple two-dimensional deceleration flow, modelling the foam spilling over a weir. Foam microstructure is given in terms of the number of channels (Plateau borders) per unit area, which scales as the inverse square of bubble size. The Plateau border number density decreases with height in the foam, and also decreases horizontally as the weir is approached. Foam drainage equations, applicable in the dry foam limit, are described. These can be used to determine the average cross-sectional area of a Plateau border, denoted A, as a function of position in the foam. Quasi-one-dimensional solutions are available in which A only varies vertically, in spite of the two-dimensional nature of the air flow and Plateau border number density fields. For such situations the liquid drainage relative to the air flow is purely vertical. The parametric behaviour of the system is investigated with respect to a number of dimensionless parameters: K (the strength of capillary suction relative to gravity), α (the deceleration of the air flow), and n and h (respectively, the horizontal and vertical variations of the Plateau border number density). The parameter K is small, implying the existence of boundary layer solutions: capillary suction is negligible except in thin layers near the bottom boundary. The boundary layer thickness (when converted back to dimensional variables) is independent of the height of the foam. The deceleration parameter α affects the Plateau border area on the top boundary: weaker decelerations give larger Plateau border areas at the surface. For weak decelerations, there is rapid convergence of the boundary layer solutions at the bottom onto ones with negligible capillary suction higher up. For strong decelerations, two branches of solutions for A are possible in the K = 0 limit: one is smooth, and the other has a distinct kink. The full system, with small but non-zero capillary suction, lies relatively close to the kinked solution branch, but convergence from the lower boundary layer onto this branch is distinctly slow. Variations in the Plateau border number density (non-zero n and h) increase individual Plateau border areas relative to the case of uniformly sized bubbles. For strong decelerations and negligible capillarity, solutions closely follow the kinked solution branch if bubble sizes are only slightly non-uniform. As the extent of non-uniformity increases, the Plateau border area reaches a maximum corresponding to no net upward velocity of foam liquid. In the case of vertical variation of number density, liquid content profiles and Plateau border area profiles cease to be simply proportional to one another. Plateau border areas match at the top of the foam independent of h, implying a considerable difference in liquid content for foams which exhibit different number density profiles. Received 3 July 2001     

Optical characterization of methanol adsorption on the bare and oxygen precovered Cu(1 1 0) surface

We have studied the adsorption and reaction of methanol on the bare and oxygen precovered Cu(1 1 0) surface at 200 K using reflectance difference spectroscopy (RDS). On the bare and fully oxygen covered surface, the sticking coefficient is close to zero. In contrast, on the partially oxygen covered surface, a sticking coefficient close to unity is obtained. This observation suggests a high mobility of methanol on both bare and oxygen covered Cu(1 1 0) and of methoxy on Cu(1 1 0). Two reaction regimes, an oxygen supply limited and an adsorption site limited regime are identified. The transition between these two regimes occurs for an oxygen coverage of about 0.2.     

Nonlinear non-spherical oscillations of gas bubble under a periodic variation of ambient liquid pressure

A mathematical model is constructed for the bubble dynamics, in which the interphase surface variation is presented in the form of a series in spherical harmonics, and the equations are written with the accuracy up to the squared amplitude of the distortion of the spherical shape of the bubble. In the oscillation regimes close to periodic sonoluminescence of a single bubble in a standing acoustic wave, the character of air bubble oscillations in water was studied depending on the bubble initial radius and the amplitude of the liquid pressure variation. It was found that non-spherical oscillations of bounded amplitude can take place outside the region of linearly stable spherical oscillations. Both the oscillations with a period equal to one or several periods of the liquid pressure variation and aperiodic oscillations are observed. It is shown that neglecting the distortions in the form of spherical harmonics with large numbers (i > 3) may lead to a change of oscillation regimes. The influence of distortions on the bubble surface shape for the harmonics with i > 8 is insignificant.     

Laser image measurement of twin bubbles formation in shear-thinning fluids

A laser image system for investigating twin bubbles formation in shear-thinning fluid was established. The process of twin-bubble formation could be directly visualized and real-time recorded through computer by means of He–Ne laser as light source using the beam expanding and light amplification technology. The shape and size of bubbles generating in carboxymethylcellulose (CMC) aqueous solutions were studied experimentally at orifice diameter 1 mm, 1.6 mm and 2.4 mm, the orifices interval 1D_o, 2D_o and 3D_o (D_o: orifice diameter) and the gas flow rate from 0.1 to 1.0 ml/s, respectively. The effects of solution mass concentration, orifice diameter and orifice interval on bubble detachment volume were investigated. The results reveals that twin bubbles gradually touch each other and then deviate from the vertical axis crossing the middle point of the line joining the two orifice during the formation process. However compared with the perfect teardrop terminal shapes in glycerol solution, the bubbles formed in CMC solutions are stretched vertically due to the shear-thinning effect of fluids. The bubble detachment volume increases with the solution mass concentration, whereas decreases with orifice diameter. The detachment volume generated at twin orifices is less affected by orifices interval, but still smaller than that at single orifice.     

Dynamics of laser-induced bubble collapse visualized by time-resolved optical shadowgraph

Abstract  

Following the first shock wave generation and the successive single bubble expansion after the breakdown by the Nd:YAG laser pulse with 35 mJ and 10 ns in distilled water, the strong secondary shock wave is generated at the instant of the bubble collapse. The single bubble expands up to 0.59 mm in radius, and then closes up by the pressure difference between the ambient liquid pressure at 10² kPa and the vapor pressure inside the bubble at 2 kPa. The maximum pressure up to 3 GPa is attained without the strong rebounding surface motion at about 93 μs after the laser shedding. We present time-resolved velocity measurements for estimating the extreme peak pressures of the first and second shock waves with the Rankine–Hugoniot analysis.     

Structure of random foam

The Surface Evolver was used to compute the equilibrium microstructure of dry soap foams with random structure and a wide range of cell-size distributions. Topological and geometric properties of foams and individual cells were evaluated. The theory for isotropic Plateau polyhedra describes the dependence of cell geometric properties on their volume and number of faces. The surface area of all cells is about 10% greater than a sphere of equal volume; this leads to a simple but accurate theory for the surface free energy density of foam. A novel parameter based on the surface-volume mean bubble radius R32 is used to characterize foam polydispersity. The foam energy, total cell edge length, and average number of faces per cell all decrease with increasing polydispersity. Pentagonal faces are the most common in monodisperse foam but quadrilaterals take over in highly polydisperse structures.     

The escape transition of a polymer: A unique case of non-equivalence between statistical ensembles

A flexible polymer chain under good solvent conditions, end-grafted on a flat repulsive substrate surface and compressed by a piston of circular cross-section with radius L may undergo the so-called “escape transition” when the height of the piston D above the substrate and the chain length N are in a suitable range. In this transition, the chain conformation changes from a quasi-two-dimensional self-avoiding walk of “blobs” of diameter D to an inhomogeneous “flower” state, consisting of a “stem” (stretched string of blobs extending from the grafting site to the piston border) and a “crown” outside of the confining piston. The theory of this transition is developed using a Landau free-energy approach, based on a suitably defined (global) order parameter and taking also effects due to the finite chain length N into account. The parameters of the theory are determined in terms of known properties of limiting cases (unconfined mushroom, chain confined between infinite parallel walls). Due to the non-existence of a local order parameter density, the transition has very unconventional properties (negative compressibility in equilibrium, non-equivalence between statistical ensembles in the thermodynamic limit, etc.). The reasons for this very unusual behavior are discussed in detail. Using Molecular Dynamics (MD) simulation for a simple bead-spring model, with N in the range 50 N 300 , a comprehensive study of both static and dynamic properties of the polymer chain was performed. Even though for the considered rather short chains the escape transition is still strongly rounded, the order parameter distribution does reveal the emerging transition clearly. Time autocorrelation functions of the order parameter and first passage times and their distribution indicate clearly the strong slowing down associated with the chain escape. The theory developed here is in good agreement with all these simulation results.     

Imaging properties of photon sieve with a large aperture

We report the optimization design and experimental results for the imaging properties of a photon sieve, which is formed on a layer of metal film supported by a thin glass substrate. As an example, we considered a micro-optical element with parameters of diameter D=50 mm, 3,564,290 hole number, and 10 μm minimum micro-hole diameter, which was designed and fabricated by means of surface machining technique in the lab. To evaluate its imaging performance, both on-axis and off-axis imaging experiments were carried out using the element. Compared to a Fresnel zone plate lens with the same feature size, the photon sieve has super imaging performance. Some quantitative analyses and initial qualitative explanations were given for the imaging characteristics.     

Statics and dynamics of adhesion between two soap bubbles

An original set-up is used to study the adhesive properties of two hemispherical soap bubbles put into contact. The contact angle at the line connecting the three films is extracted by image analysis of the bubbles profiles. After the initial contact, the angle rapidly reaches a static value slightly larger than the standard 120^° angle expected from Plateau rule. This deviation is consistent with previous experimental and theoretical studies: it can be quantitatively predicted by taking into account the finite size of the Plateau border (the liquid volume trapped at the vertex) in the free energy minimization. The visco-elastic adhesion properties of the bubbles are further explored by measuring the deviation Δθ_d(t) of the contact angle from the static value as the distance between the two bubbles supports is sinusoidally modulated. It is found to linearly increase with Δr _c/r _c , where r_c is the radius of the central film and Δr _c the amplitude of modulation of this length induced by the displacement of the supports. The in-phase and out-of-phase components of Δθ_d(t) with the imposed modulation frequency are systematically probed, which reveals a transition from a viscous to an elastic response of the system with a crossover pulsation of the order 1rad · s^-1. Independent interfacial rheological measurements, obtained from an oscillating bubble experiment, allow us to develop a model of dynamic adhesion which is confronted to our experimental results. The relevance of such adhesive dynamic properties to the rheology of foams is briefly discussed using a perturbative approach to the Princen 2D model of foams.     

Vortex dynamics of a trapezoidal bluff body placed inside a circular pipe

The influence of Reynolds number and blockage ratio on the vortex dynamics of a trapezoidal bluff body placed inside a circular pipe is studied experimentally and numerically. Low aspect ratio, high blockage ratio, curved end conditions (junction of pipe and bluff body), axisymmetric upstream flow with shear and turbulence are some of the intrinsic features of this class of bluff body flows which have been scarcely addressed in the literature. A large range (200:200,000) of Reynolds number (Re_D) is covered in this study, encompassing all the three pipe flow regimes (laminar, transition and turbulent). Four different flow regimes are defined based on the distinct features of Strouhal number (St)–Re_D relation: steady, laminar irregular, transition and turbulent. The wake in the steady regime is stationary with no oscillations in the shear layer. The laminar regime is termed as irregular owing to irregular vortex shedding. The vortex shedding in this regime is observed to be symmetric. The emergence of separation bubble downstream of the bluff body on either side is another interesting feature of this regime, which is further observed to be symmetric. Two pairs of mean streamwise vortices are noticed in the near-wake regime, which are termed as reverse dipole-type wake topology. Beyond the irregular laminar regime, the Strouhal number falls gradually and vortex shedding becomes more periodic. This regime is named transition and occurs close to the Reynolds number at which transition to turbulence takes place in a fully developed pipe. The turbulent regime is characterised by a nearly constant Strouhal number. Typical Karman-type vortex shedding is noticed in this regime. The convection velocity, wake width formation length and irrecoverable pressure loss are quantified to highlight the influence of blockage ratio. These results will be useful to develop basic understanding of vortex dynamics of confined bluff body flow for several practical applications.     

Control of the dynamics of a boiling vapour bubble using pressure-modulated high intensity focused ultrasound without the shock scattering effect: A first proof-of-concept study

Boiling histotripsy is a promising High-Intensity Focused Ultrasound (HIFU) technique that can be used to induce mechanical tissue fractionation at the HIFU focus via cavitation. Two different types of cavitation produced during boiling histotripsy exposure can contribute towards mechanical tissue destruction: (1) a boiling vapour bubble at the HIFU focus and (2) cavitation clouds in between the boiling bubble and the HIFU source. Control of the extent and degree of mechanical damage produced by boiling histotripsy is necessary when treating a solid tumour adjacent to normal tissue or major blood vessels. This is, however, difficult to achieve with boiling histotripsy due to the stochastic formation of the shock scattering-induced inertial cavitation clouds. In the present study, a new histotripsy method termed pressure-modulated shockwave histotripsy is proposed as an alternative to or in addition to boiling histotripsy without inducing the shock scattering effect. The proposed concept is (a) to generate a boiling vapour bubble via localised shockwave heating and (b) subsequently control its extent and lifetime through manipulating peak pressure magnitudes and a HIFU pulse length. To demonstrate the feasibility of the proposed method, bubble dynamics induced at the HIFU focus in an optically transparent liver tissue phantom were investigated using a high speed camera and a passive cavitation detection systems under a single 10, 50 or 100 ms-long 2, 3.5 or 5 MHz pressure-modulated HIFU pulse with varying peak positive and negative pressure amplitudes from 5 to 89 MPa and −3.7 to −14.6 MPa at the focus. Furthermore, a numerical simulation of 2D nonlinear wave propagation with the presence of a boiling bubble at the focus of a HIFU field was conducted by numerically solving the generalised Westervelt equation. The high speed camera experimental results showed that, with the proposed pressure-modulated shockwave histotripsy, boiling bubbles generated by shockwave heating merged together, forming a larger bubble (of the order of a few hundred micron) at the HIFU focus. This coalesced boiling bubble then persisted and maintained within the HIFU focal zone until the end of the exposure (10, 50, or 100 ms). Furthermore, and most importantly, no violent cavitation clouds which typically appear in boiling histotripsy occurred during the proposed histotripsy excitation (i.e. no shock scattering effect). This was likely because that the peak negative pressure magnitude of the backscattered acoustic field by the boiling bubble was below the cavitation cloud intrinsic threshold. The size of the coalesced boiling bubble gradually increased with the peak pressure magnitudes. In addition, with the proposed method, an oval shaped lesion with a length of 0.6 mm and a width of 0.1 mm appeared at the HIFU focus in the tissue phantom, whereas a larger lesion in the form of a tadpole (length: 2.7 mm, width: 0.3 mm) was produced by boiling histotripsy. Taken together, these results suggest that the proposed pressure-modulated shockwave histotripsy could potentially be used to induce a more spatially localised tissue destruction with a desired degree of mechanical damage through controlling the size and lifetime of a boiling bubble without the shock scattering effect.