Gas-liquid nucleation at large metastability: unusual features and a new formalism

Nucleation at large metastability is still largely an unsolved problem, even though it is a problem of tremendous current interest, with wide-ranging practical value, from atmospheric research to materials science. It is now well accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation-free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this paper, we present an alternative formalism to treat nucleation at large supersaturation by introducing an extended set of order parameters in terms of the kth largest liquid-like clusters, where k = 1 is the largest cluster in the system, k = 2 is the second largest cluster and so on. At low supersaturation, the size of the largest liquid-like cluster acts as a suitable order parameter. At large supersaturation, the free energy barrier for the largest liquid-like cluster disappears. We identify this supersaturation as the one at the onset of kinetic spinodal. The kinetic spinodal is system-size-dependent. Beyond kinetic spinodal many clusters grow simultaneously and competitively and hence the nucleation and growth become collective. In order to describe collective growth, we need to consider the full set of order parameters. We derive an analytic expression for the free energy of formation of the kth largest cluster. The expression predicts that, at large metastability (beyond kinetic spinodal), the barrier of growth for several largest liquid-like clusters disappears, and all these clusters grow simultaneously. The approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts weaker supersaturation dependence than what is predicted by CNT at large metastability. Such a crossover behavior has indeed been observed in recent experiments (but eluded an explanation till now). In order to understand the large numerical discrepancy between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interactions of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for the Lennard-Jones potential, both in two and three dimensions. The value of surface tension and also the free energy difference between the gas and the liquid phase increase significantly and converge only when the range of interaction is extended beyond 6-7 molecular diameters. We find, with the full range of interaction potential, that the surface tension shows only a weak dependence on supersaturation, so the reason for the breakdown of CNT (with simulated values of surface tension and free energy gap) cannot be attributed to the supersaturation dependence of surface tension. This remains an unsettled issue at present because of the use of the value of surface tension obtained at coexistence.

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Production of biosurfactant “Biosur-Pm” by Pseudomonas maltophila CSV89: characterization and role in hydrocarbon uptake

Pseudomonas maltophilia CSV89, a soil bacterium, produces an extracellular biosurfactant, ''Biosur-Pm''. The partially purified product is nondialyzable and chemically composed of 50% protein and 12-15% sugar, which indicates the complex nature of Biosur-Pm. It reduces the surface tension of water from 73 to 53 x 10(-3) N m(-1) and has a critical micellar concentration of 80 mg/l. Compared to aliphatic hydrocarbons, Biosur-Pm shows good activity against aromatic hydrocarbons. The emulsion formed is stable and does not require any metal ions for emulsification. The kinetics of Biosur-Pm production suggest that its synthesis isa growth-associated and pH-dependent process. At pH 7.0, cells produced more Biosur-Pm with less cell surface hydrophobicity. At pH 8.0, however, the cells produced less Biosur-Pm with more cell surface hydrophobicity and showed a twofold higher affinity for aromatic hydrocarbons compared to the cells grown at pH 7.0. The Biosur-Pm showed a pH-dependent release, stimulated growth of the producer strain on mineral salts medium with 1-naphthoic acid when added externally, and facilitated the conversion of salicylate to catechol. All these results suggest that Biosur-Pm is probably a cell-wall component and helps in hydrocarbon assimilation/uptake.

Levitation of a drop over a film flow

A vertical jet of water impinging on a horizontal surface produces a radial film flow followed by a circular hydraulic jump. We report a phenomenon where fairly large (1 mi) drops of liquid levitate just upstream of the jump on a thin air layer between the drop and the film flow. We explain the phenomenon using lubrication theory. Bearing action both in the air film and the water film seems to be necessary to support large drops. Horizontal support is given to the drop by the hydraulic jump. A variety of drop shapes is observed depending on the volume of the drop and liquid properties. We show that interaction of the forces due to gravity, surface tension, viscosity and inertia produces these various shapes.

Pronounced hydrogel formation by the self-assembled aggregates of N-alkyl disaccharide amphiphiles

Six disaccharide amphiphiles were synthesized and their hydrogel-forming behavior was extensively studied. These amphiphiles were based on maltose and lactose. Since the gels formed from some of these systems showed the ability to "trap" water molecules upon gelation, these gels were described as "hydrogels". When these gels were heated to similar to 70 degrees C, the samples turned into clear, isotropic fluids, and upon gradual cooling, the hydrogels could be reproduced. Thus these systems were also "thermoreversible". The low molecular mass (MW 565) of the gelators compared to that of a typical polymeric gelator forming substance implies pronounced aggregation of the disaccharide amphiphiles into larger microstructures during gelation. To discern the aggregate textures and morphologies, the specimen hydrogel samples were examined by high-resolution scanning electron microscopy (SEM). A possible reason for the exceptionally high water gelating capacities (>6000 molecules of water per gelator molecule) exhibited by these N-alkyl disaccharide amphiphiles is the presence of large interlamellar spaces into which the water molecules get entrapped due to surface tension. In contrast to their single-chain counterparts, the double-chain lactosyl and maltosylamine amphiphiles upon solubilization in EtOH-H2O afforded hydrogels with reduced mechanical strengths. Interestingly, the corresponding microstructures were found to be quite different from the corresponding hydrogels of their single-chain counterparts. Rheological studies provided further insights into the behavior of these hydrogels. Varying the chain length of the alcohol cosolvent could modulate the gelation capacities, melting temperatures, and the mechanical properties of these hydrogels. To explain the possible reasons of gelation, the results of molecular modeling and energy minimization studies were also included.

Thermodynamics of surfaces and adsorption in the Fe-S, Fe-N, and Fe-S-N systems at 1823 K

The present investigation calculates the surface tension and adsorption functions of the Fe-S, Fe-N, and Fe-S-N melts at 1823 K using the modified form of Butler's equations and the derived values of the surface interaction parameters of the systems. The calculated values are found to be in good agreement with those of the experimental data of the systems. The present analysis indicates similar adsorption behavior of sulfur for the Fe-S and Fe-S-N systems at 1823 K. Although a linear adsorption behavior is observed in the Fe-N system, an inverse relationship in the adsorption behavior exists between sulur and nitrogen in the Fe-S-N system.

An inverse problem for the Helmholtz equation involving two semi-infinite fluids

An analytical method is developed for solving an inverse problem for Helmholtz's equation associated with two semi-infinite incompressible fluids of different variable refractive indices, separated by a plane interface. The unknowns of the inverse problem are: (i) the refractive indices of the two fluids, (ii) the ratio of the densities of the two fluids, and (iii) the strength of an acoustic source assumed to be situated at the interface of the two fluids. These are determined from the pressure on the interface produced by the acoustic source. The effect of the surface tension force at the interface is taken into account in this paper. The application of the proposed analytical method to solve the inverse problem is also illustrated with several examples. In particular, exact solutions of two direct problems are first derived using standard classical methods which are then used in our proposed inverse method to recover the unknowns of the corresponding inverse problems. The results are found to be in excellent agreement.

Thermal instability in a three-dimensional rigid container with prescribed heat flux at lower boundary

The Bénard–Marangoni convection is studied in a three-dimensional container with thermally insulated lateral walls and prescribed heat flux at lower boundary. The upper surface of the incompressible, viscous fluid is assumed to be flat with temperature dependent surface tension. A Galerkin–Tau method with odd and even trial functions satisfying all the essential boundary conditions except the natural boundary conditions at the free surface has been used to solve the problem. The critical Marangoni and Rayleigh numbers are determined for the onset of steady convection as a function of aspect ratios x0 and y0 for the cases of Bénard–Marangoni, pure Marangoni and pure Bénard convections. It is observed that critical parameters are decreasing with an increase in aspect ratios. The flow structures corresponding to the values of the critical parameters are presented in all the cases. It is observed that the critical parameters are higher for case with heat flux prescribed than those corresponding to the case with prescribed temperature. The critical Marangoni number for pure Marangoni convection is higher than critical Rayleigh number corresponding to pure Bénard convection for a given aspect ratio whereas the reverse was observed for two-dimensional infinite layer.

Instabilities of soft elastic microtubes filled with viscous fluids: Pearls, wrinkles, and sausage strings

A linear stability analysis is presented to study the self-organized instabilities of a highly compliant elastic cylindrical shell filled with a viscous liquid and submerged in another viscous medium. The prototype closely mimics many components of micro-or nanofluidic devices and biological processes such as the budding of a string of pearls inside cells and sausage-string formation of blood vessels. The cylindrical shell is considered to be a soft linear elastic solid with small storage modulus. When the destabilizing capillary force derived from the cross-sectional curvature overcomes the stabilizing elastic and in-plane capillary forces, the microtube can spontaneously self-organize into one of several possible configurations; namely, pearling, in which the viscous fluid in the core of the elastic shell breaks up into droplets; sausage strings, in which the outer interface of the mircrotube deforms more than the inner interface; and wrinkles, in which both interfaces of the thin-walled mircrotube deform in phase with small amplitudes. This study identifies the conditions for the existence of these modes and demonstrates that the ratios of the interfacial tensions at the interfaces, the viscosities, and the thickness of the microtube play crucial roles in the mode selection and the relative amplitudes of deformations at the two interfaces. The analysis also shows asymptotically that an elastic fiber submerged in a viscous liquid is unstable for Y = gamma/(G(e)R) > 6 and an elastic microchannel filled with a viscous liquid should rupture to form spherical cavities (pearling) for Y > 2, where gamma, G(e), and R are the surface tension, elastic shear modulus, and radius, respectively, of the fiber or microchannel.

Thermally induced secondary atomization of droplet in an acoustic field

We study the thermal effects that lead to instability and break up in acoustically levitated vaporizing fuel droplets. For selective liquids, atomization occurs at the droplet equator under external heating. Short wavelength [Kelvin-Helmholtz (KH)] instability for diesel and bio-diesel droplets triggers this secondary atomization. Vapor pressure, latent heat, and specific heat govern the vaporization rate and temperature history, which affect the surface tension gradient and gas phase density, ultimately dictating the onset of KH instability. We develop a criterion based on Weber number to define a condition for the inception of secondary atomization. (C) 2012 American Institute of Physics. [doi:10.1063/1.3680257]

A k-epsilon model for turbulent weld pool convection in gas metal arc welding process

In this paper, we present a modified k - epsilon model capable of addressing turbulent weld-pool convection in a GMAW process, taking into account the morphology of the phase change interface during a Gas Metal Arc Welding (GMAW) process. A three-dimensional turbulence mathematical model has been developed to study the heat transfer and fluid flow within the weld pool by considering the combined effect of three driving forces, viz., buoyancy, Lorentz force and surface tension (Marangoni convection). Mass and energy transports by the droplets are considered through the thermal analysis of the electrode. The falling droplet's heat addition to the molten pool is considered to be a volumetric heat source distributed in an imaginary cylindrical cavity ("cavity model") within the weld pool. This nature of heat source distribution takes into account the momentum and the thermal, energy of the falling droplets. The numerically predicted weld pool dimensions both from turbulence and laminar models are then compared with the experimental post-weld results sectioned across the weld axis. The above comparison enables us to analyze the overall effects of turbulent convection on the nature of heat and fluid flow and hence on the weld pool shape/size during the arc welding processes.

Sensitivity of nucleation phenomena on range of interaction potential

Theoretical and computational investigations of nucleation have been plagued by the sensitivity of the phase diagram to the range of the interaction potential. As the surface tension depends strongly on the range of interaction potential and as the classical nucleation theory (CNT) predicts the free energy barrier to be directly proportional to the cube of the surface tension, one expects a strong sensitivity of nucleation barrier to the range of the potential; however, CNT leaves many aspects unexplored. We find for gas-liquid nucleation in Lennard-Jones system that on increasing the range of interaction the kinetic spinodal (KS) (where the mechanism of nucleation changes from activated to barrierless) shifts deeper into the metastable region. Therefore the system remains metastable for larger value of supersaturation and this allows one to explore the high metastable region without encountering the KS. On increasing the range of interaction, both the critical cluster size and pre-critical minima in the free energy surface of kth largest cluster, at respective kinetic spinodals, shift towards smaller cluster size. In order to separate surface tension contribution to the increase in the barrier from other non-trivial factors, we introduce a new scaling form for surface tension and use it to capture both the temperature and the interaction range dependence of surface tension. Surprisingly, we find only a weak non-trivial contribution from other factors to the free energy barrier of nucleation. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.3685835]

Exciting forces due to interaction of water waves with a submerged sphere in an ice-covered two-layer fluid of finite depth

Scattering of water waves by a sphere in a two-layer fluid, where the upper layer has an ice-cover modelled as an elastic plate of very small thickness, while the lower one has a rigid horizontal bottom surface, is investigated within the framework of linearized water wave theory. The effects of surface tension at the surface of separation is neglected. There exist two modes of time-harmonic waves - the one with lower wave number propagating along the ice-cover and the one with higher wave number along the interface. Method of multipole expansions is used to find the particular solution for the problem of wave scattering by a submerged sphere placed in either of the layers. The exciting forces for vertical and horizontal directions are derived and plotted against different values of the wave number for different submersion depths of the sphere and flexural rigidity of the ice-cover. When the flexural rigidity and the density of the ice-cover are taken to be zero, the numerical results for the exciting forces for the problem with free surface are recovered as particular cases. (C) 2011 Elsevier Ltd. All rights reserved.

Surfactant solubility and micellization in ternary eutectic melt (acetamide plus urea plus ammonium nitrate)

Cationic surfactants such as dodecyltrimethylammonium bromide (DTAB), tetradecyltrimehtylammonium bromide (TTAB) and hexadecyltrimethylammonium bromide (HTAB); and anionic surfactants such as sodium decyl sulphate (SDeS), sodium dodecyl sulphate (SDS) and sodium tetradecyl sulphate (STDS) have been used to determine their solubility and micellization in ternary eutectic melt (acetamide + urea + ammonium nitrate) at 50 degrees C. We employed the electrical conductivity and the surface tension measurement techniques to determine the critical micelle concentration (CMC). The deviation in the slope of the specific conductance/surface tension against surfactant concentration plots indicated the aggregations of surfactants and hence, their CMC. CMC decreases with increase of alkyl chain length due to the increased van der Waals forces. The calculated increment in Gibb's energy per methylene group for cationic and anionic surfactants is about -6 kJ mol(-1) and -4 kJ mol(-1) respectively. It is found that, the CMCs of the surfactants in the ternary melt are higher than the CMCs of same surfactants in water (similar to 25 degrees C). (C) 2012 Elsevier B.V. All rights reserved.

Directional motion of impacting drops on dual-textured surfaces

In this work, we analyze the directional movement of impacting liquid drops on dual-textured solid surfaces comprising two different surface morphologies: a textured surface and a smooth surface. The dynamics of liquid drops impacting onto the junction line between the two parts of the dual-textured surfaces is studied experimentally for varying drop impact velocity. The dual-textured surfaces used here featured a variation in their textures' geometrical parameters as well as their surface chemistry. Two types of liquid drop differing in their surface tension were used. The impact process develops a net horizontal drop velocity towards the higher-wettability surface portion and results in a bulk movement of the impacting drop liquid. The final distance moved by the impacting drop from the junction line decreases with increasing impacting drop Weber number We. A fully theoretical model, employing a balance of forces acting at the drop contact line as well as energy conservation, is formulated to determine the variation, with We, of net horizontal drop velocity and subsequent movement of the impacting drop on the dual-textured surfaces.