Pendant drop formation of shear-thinning and yield stress fluids

A volume-of-fluid numerical method is used to predict the dynamics of shear-thinning liquid drop formation in air from a circular orifice. The validity of the numerical calculation is confirmed for a Newtonian liquid by comparison with experimental measurements. For particular values of Weber number and Froude number, predictions show a more rapid pinch-off, and a reduced number of secondary droplets, with increasing shear-thinning. Also a minimum in the limiting drop length occurs for the smallest Weber number as the zero-shear viscosity is varied. At the highest viscosity, the drop length is reduced due to shear-thinning, whereas at lower viscosities there is little effect of shear-thinning. The evolution of predicted drop shape, drop thickness and length, and the configuration at pinch-off are discussed for shear-thinning drops. The evolution of a drop of Bingham yield stress liquid is also considered as a limiting case. In contrast to the shear-thinning cases, it exhibits a plug flow prior to necking, an almost step-change approach to pinch-off of a torpedo shaped drop following the onset of necking, and a much smaller neck length; no secondary drops are formed. The results demonstrate the potential of the numerical model as a design tool in tailoring the fluid rheology for controlling drop formation behaviour. (c) 2006 Elsevier Inc. All rights reserved.

The effects of energy sites on adsorption of Lennard–Jones fluids and phase transition in carbon slit pore of finite length a computer simulation study

A Monte Carlo simulation method is Used 10 study the effects of adsorption strength and topology of sites on adsorption of simple Lennard-Jones fluids in a carbon slit pore of finite length. Argon is used as a model adsorbate, while the adsorbent is modeled as a finite carbon slit pore whose two walls composed of three graphene layers with carbon atoms arranged in a hexagonal pattern. Impurities having well depth of interaction greater than that of carbon atom are assumed to be grafted onto the surface. Different topologies of the impurities; corner, centre, shelf and random topologies are studied. Adsorption isotherms of argon at 87.3 K are obtained for pore having widths of 1, 1.5 and 3 11111 using a Grand Canonical Monte Carlo simulation (GCMC). These results are compared with isotherms obtained for infinite pores. It is shown that the Surface heterogeneity affects significantly the overall adsorption isotherm, particularly the phase transition. Basically it shifts the onset of adsorption to lower pressure and the adsorption isotherms for these four impurity models are generally greater than that for finite pore. The positions of impurities on solid Surface also affect the shape of the adsorption isotherm and the phase transition. We have found that the impurities allocated at the centre of pore walls provide the greatest isotherm at low pressures. However when the pressure increases the impurities allocated along the edges of the graphene layers show the most significant effect on the adsorption isotherm. We have investigated the effect of surface heterogeneity on adsorption hysteresis loops of three models of impurity topology, it shows that the adsorption branches of these isotherms are different, while the desorption branches are quite close to each other. This suggests that the desorption branch is either the thermodynamic equilibrium branch or closer to it than the adsorption branch. (c) 2005 Elsevier Inc. All rights reserved.

Structure of saccharose-based carbon and transport of confined fluids: hybrid reverse Monte Carlo reconstruction and simulation studies

We present results of the reconstruction of a saccharose-based activated carbon (CS1000a) using hybrid reverse Monte Carlo (HRMC) simulation, recently proposed by Opletal et al. [1]. Interaction between carbon atoms in the simulation is modeled by an environment dependent interaction potential (EDIP) [2,3]. The reconstructed structure shows predominance of sp(2) over sp bonding, while a significant proportion of sp(3) hybrid bonding is also observed. We also calculated a ring distribution and geometrical pore size distribution of the model developed. The latter is compared with that obtained from argon adsorption at 87 K using our recently proposed characterization procedure [4], the finite wall thickness (FWT) model. Further, we determine self-diffusivities of argon and nitrogen in the constructed carbon as functions of loading. It is found that while there is a maximum in the diffusivity with respect to loading, as previously observed by Pikunic et al. [5], diffusivities in the present work are 10 times larger than those obtained in the prior work, consistent with the larger pore size as well as higher porosity of the activated saccharose carbon studied here.

Transport of simple fluids in nanopores: Theory and simulation

A theory is discussed of single-component transport in nanopores, recently developed by Bhatia and coworkers. The theory considers the oscillatory motion of molecules between diffuse wall collisions, arising from the fluid-wall interaction, along with superimposed viscous flow due to fluid-fluid interaction. The theory is tested against molecular dynamics simulations for hydrogen, methane, and carbon tetrafluoride flow in cylindrical nanopores in silica. Although exact at low densities, the theory performs well even at high densities, with the density dependency of the transport coefficient arising from viscous effects. Such viscous effects are reduced at high densities because of the large increase in viscosity, which explains the maximum in the transport coefficient with increase in density. Further, it is seen that in narrow pore sizes of less than two molecular diameters, where a complete monolayer cannot form on the surface, the mutual interference of molecules on opposite sides of the cross section can reduce the transport coefficient, and lead to a maximum in the transport coefficient with increasing density. The theory is also tested for the case of partially diffuse reflection and shows the viscous contribution to be negligible when the reflection is nearly specular. (c) 2005 American Institute of Chemical Engineers AIChE J, 52: 29-38, 2006.

Effects of graphene layer size on the adsorption of fluids on graphitized thermal carbon black. A computer simulation study

The adsorption of Lennard-Jones fluids (argon and nitrogen) onto a graphitized thermal carbon black surface was studied with a Grand Canonical Monte Carlo Simulation (GCMC). The surface was assumed to be finite in length and composed of three graphene layers. When the GCMC simulation was used to describe adsorption on a graphite surface, an over-prediction of the isotherm was consistently observed in the pressure regions where the first and second layers are formed. To remove this over-prediction, surface mediation was accounted for to reduce the fluid-fluid interaction. Do and co-workers have introduced the so-called surface-mediation damping factor to correct the over-prediction for the case of a graphite surface of infinite extent, and this approach has yielded a good description of the adsorption isotherm. In this paper, the effects of the finite size of the graphene layer on the adsorption isotherm and how these would affect the extent of the surface mediation were studied. It was found that this finite-surface model provides a better description of the experimental data for graphitized thermal carbon black of high surface area (i.e. small crystallite size) while the infinite- surface model describes data for carbon black of very low surface area (i.e. large crystallite size).

Using urea dilution to standardise cellular and non-cellular components of pleural and bronchoalveolar lavage (BAL) fluids in the cat

A technique to standardise the analysis of cellular and non-cellular components in epithelial lining fluid (ELF) collected during saline lavage of pulmonary and pleural cavities was developed using the urea dilution method. Bronchoalveolar lavage (BAL) and pleural lavage (PL) fluids were collected from 12 clinically healthy cats. Total and differential cell counts in BAL fluid were within normal ranges for the cat, while cell Counts in PL fluid were assumed to be normal based on clinical health during examination, auscultation and lactate dehydrogenase (LDH) activities being comparable with other species. The major clinical implication of this study was that nucleated cell counts within feline ELF could not be predicted from analysis of lavage fluid which suggests that calculation of the proportion of ELF in lavage fluid by the urea dilution method may be necessary to avoid misdiagnosis of health or disease in pulmonary or pleural cavities. (C) 2005 ESFM and AAFP. Published by Elsevier Ltd. All rights reserved.

Comparison of surfactant lipids between pleural and pulmonary lining fluids

Saturated phospholipids (PCs), particularly dipalmitoylphosphatidylcholine (DPPC), predominate in surfactant lining the alveoli, although little is known about the relationship between saturated and unsaturated PCs on the outer surface of the lung, the pleura. Seven healthy cats were anesthetized and a bronchoalveolar lavage (BAL) was performed, immediately followed by a pleural lavage (PL). Lipid was extracted from lavage fluid and then analyzed for saturated, primarily dipalmitoylphosphatidylcholine (DPPC), and unsaturated PC species using high-performance liquid chromatography (HPLC) with combined fluorescence and ultraviolet detection. Dilution of epithelial lining fluid (ELF) in lavage fluids was corrected for using the urea method. The concentration of DPPC in BAL fluid (85.3 +/- 15.7 mu g/mL) was significantly higher (P=0.021) than unsaturated PCs (similar to 40 mu g/mL). However, unsaturated PCs (similar to 34 mu g/mL), particularly stearoyl-linoleoyl-phosphatidylcholine (SLPC; 17.4 +/- 6.8), were significantly higher (P = 0.021) than DPPC (4.3 +/- 1.8 mu g/mL) in PL fluid. These results show that unsaturated PCs appear functionally more important in the pleural cavity, which may have implications for surfactant replenishment following pleural disease or thoracic surgery. (c) 2005 Published by Elsevier Ltd.

Drop formation and breakup of low viscosity elastic fluids: Effects of molecular weight and concentration

The dynamics of drop formation and pinch-off have been investigated for a series of low viscosity elastic fluids possessing similar shear viscosities, but differing substantially in elastic properties. On initial approach to the pinch region, the viscoelastic fluids all exhibit the same global necking behavior that is observed for a Newtonian fluid of equivalent shear viscosity. For these low viscosity dilute polymer solutions, inertial and capillary forces form the dominant balance in this potential flow regime, with the viscous force being negligible. The approach to the pinch point, which corresponds to the point of rupture for a Newtonian fluid, is extremely rapid in such solutions, with the sudden increase in curvature producing very large extension rates at this location. In this region the polymer molecules are significantly extended, causing a localized increase in the elastic stresses, which grow to balance the capillary pressure. This prevents the necked fluid from breaking off, as would occur in the equivalent Newtonian fluid. Alternatively, a cylindrical filament forms in which elastic stresses and capillary pressure balance, and the radius decreases exponentially with time. A (0+1)-dimensional finitely extensible nonlinear elastic dumbbell theory incorporating inertial, capillary, and elastic stresses is able to capture the basic features of the experimental observations. Before the critical "pinch time" t(p), an inertial-capillary balance leads to the expected 2/3-power scaling of the minimum radius with time: R-min similar to(t(p)-t)(2/3). However, the diverging deformation rate results in large molecular deformations and rapid crossover to an elastocapillary balance for times t>t(p). In this region, the filament radius decreases exponentially with time R-min similar to exp[(t(p)-t)/lambda(1)], where lambda(1) is the characteristic time constant of the polymer molecules. Measurements of the relaxation times of polyethylene oxide solutions of varying concentrations and molecular weights obtained from high speed imaging of the rate of change of filament radius are significantly higher than the relaxation times estimated from Rouse-Zimm theory, even though the solutions are within the dilute concentration region as determined using intrinsic viscosity measurements. The effective relaxation times exhibit the expected scaling with molecular weight but with an additional dependence on the concentration of the polymer in solution. This is consistent with the expectation that the polymer molecules are in fact highly extended during the approach to the pinch region (i.e., prior to the elastocapillary filament thinning regime) and subsequently as the filament is formed they are further extended by filament stretching at a constant rate until full extension of the polymer coil is achieved. In this highly extended state, intermolecular interactions become significant, producing relaxation times far above theoretical predictions for dilute polymer solutions under equilibrium conditions. (C) 2006 American Institute of Physics