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

Timing and chemistry of fluid-flow events in the Lawn Hill Platform, Northern Australia

Mudrocks and carbonates of the Isa superbasin in the Lawn Hill platform in northern Australia host major base metal sulfide mineralization, including the giant strata-bound Century Zn-Pb deposit. Mineral paragenesis, stable isotope, and K-Ar dating studies demonstrate that long-lived structures such as the Termite Range fault acted as hot fluid conduits several times during the Paleoproterozoic and Mesoproterozoic in response to major tectonic events. Illite and chlorite crystallinity studies suggest the southern part of the platform has experienced higher temperatures (up to 300 degrees C) than similar stratigraphic horizons in the north. The irregular downhole variation of illite crystallinity values provides further information oil the thermal regime in the basin and shows that clay formation was controlled not only by temperature increase with depth but also by high water/rock ratios along relatively permeable zones. K-Ar dating of illite, in combination with other data, may indicate three major thermal events in the central and northern Lawn Hill platform Lit 1500, 1440 to 1400, and 1250 to 1150 Ma. This study did not detect the earlier Century base metal mineralizing event at 1575 Ma. 1500 Ma ages are recorded only in the south and correspond to the age of the Late Isan orogeny and deposition of the Lower Roper superbasin. They may reflect exhumation of a provenance region. The 1440 to 1300 Ma ages are related to fault reactivation and a thermal pulse at similar to 1440 to 1400 Ma possibly accompanied by fluid flow, with subsequent enhanced cooling possibly due to thermal relaxation or further crustal exhumation. The youngest thermal and/or fluid-flow event at 1250 to 1150 Ma is recorded mainly to the cast of the Tern-lite Range fault and may be related to the assembly of the Rodinian supercontinent. Fluids in equilibrium with illite that formed over a range of temperatures, at different times in different parts of the platform. have relatively uniform oxygen isotope compositions and more variable hydrogen isotope compositions (delta O-18 = 3.5-9.7 parts per thousand V-SMOW; delta D = -94 to -36 parts per thousand V-SMOW). The extent of the 180 enrichment and the variably depleted hydrogen isotope compositions suggest the illite interacted with deep-basin hypersaline brines that were composed of evaporated seawater and/or highly evolved meteoric water. Siderite is the most abundant iron-rich gangue phase in the Century Zn-Pb deposit, which is surrounded by all extensive ferroan carbonate alteration halo. Modeling suggests that the ore siderite formed at temperatures of 120 degrees to 150 degrees C, whereas siderite and ankerite in the alteration halo formed at temperatures of 150 degrees to 180 degrees C. The calculated isotopic compositions of the fluids are consistent with O-18-rich basinal brines and mixed inorganic and organic carbon Sources (6180 = 3-10 parts per thousand V-SMOW, delta C-13 = -7 to -3 parts per thousand V-PDB). in the northeast Lawn Hill platform carbonate-rich rocks preserve marine to early diagenetic carbon and oxygen isotope compositions, whereas ferroan carbonate cements in siltstones and shales in the Desert Creek borehole are O-18 and C-13 depleted relative to the sedimentary carbonates. The good agreement between temperature estimates from illite crystallinity and organic reflectance (160 degrees-270 degrees C) and inverse correlation with carbonate delta O-18 values indicates that organic maturation and carbonate precipitation in the northeast Lawn Hill platform resulted from interaction with the 1250 to 1150 Ma fluids. The calculated isotopic compositions of the fluid are consistent with evolved basinal brine (delta O-18 = 5.1-9.4 parts per thousand V-SMOW; delta C-13 = -13.2 to -3.7 parts per thousand V-PDB) that contained a variable organic carbon component from the oxidation and/or hydrolysis of organic matter in the host sequence. The occurrence of extensive O-18- and C-13-depleted ankerite and siderite alteration in Desert Creek is related to the high temperature of the 1250 to 1150 Ma fluid-flow event in the northeast Lawn Hill platform, in contrast to the lower temperature fluids associated with the earlier Century Zn-Pb deposit in the central Lawn Hill platform.

Zinc deposits and related mineralization of the Burketown mineral field, including the world-class Century Deposit, Northern Australia: Fluid inclusion and stable isotope evidence for basin fluid sources

The stratiform Century Zn-Pb deposit and the discordant Zn-Pb lode deposits of the Burketown mineral field, northern Australia, host ore and gangue minerals with primary fluid inclusions that have not been affected by the Isan orogeny, thus providing a unique opportunity to investigate the nature of the ore-forming brines. All of the deposits are hosted in shales and siltstones belonging to the Isa superbasin and comprise sphalerite, pyrite, carbonate, quartz, galena, minor chalcopyrite, and minor illite. According to Pb model ages, the main ore stage of mineralization at Century formed at I575 Ma, some 20 m.y. after deposition of the host shale sequence. Microthermometry on undeformed, primary fluid inclusions hosted in porous sphalerite shows that the Zn at Century was transported to the deposit by a homogeneous, Ca2+- and Na+-bearing brine with a salinity of 21.6 wt percent NaCl equiv. delta D-fluid of the fluid inclusion water ranges from -89 to -83 per mil, consistent with a basinal brine that evolved from meteoric water. Fluid inclusion homogenization temperatures range between 74 degrees and 125 degrees C, which are lower than the 120 degrees to 160 degrees C range calculated from vitrinite reflectance and illite crystallinity data from the deposit. This discrepancy indicates that mineralization likely formed at 50 to 85 Mpa, corresponding to a depth of 1,900 to 3,100 m. Transgressive galena-sphalerite veins that cut stratiform mineralization at Century and breccia-filled quartz-dolomite-sphalerite-galena veins in the discordant Zn-Pb lodes have Pb model ages between 1575 and 1485 Ma. Raman spectroscopy and microthermometry reveal that the primary fluid inclusions in these veins contain Ca2+, Na+. but they have lower salinities between 23 and 10 wt percent NaCl equiv and higher delta D-fluid values ranging from -89 to -61 per mil than fluid inclusions in porous sphalerite from Century. Fluid inclusion water from sphalerite in one of the lode deposits has delta O-18(fluid) values of 1.6 and 2.4 per mil, indistinguishable from delta O-18(fluid) values between -0.3 to +7.4 per mil calculated from the isotopic composition of coexisting quartz, dolomite, and illite. The trend toward lower salinities and higher delta D-fluid values relative to the earlier mineralizing fluids is attributed to mixing between the fluid that formed Century and a seawater-derived fluid from a different source. Based on seismic data from the Lawn Hill platform and paragenetic and geochemical results from the Leichhardt River fault trough to the south, diagenetic aquifers in the Underlying Calvert superbasin appear to have been the most likely sources for the fluids that formed Century and the discordant lode deposits. Paragenetically late sphalerite and calcite cut sphalerite, quartz, and dolomite in the lode deposits and contain Na+-dominated fluid inclusions with much lower salinities than their older counterparts. The isotopic composition of calcite also indicates delta O-18(fluid) from 3.3 to 10.7 per mil, which is larger than the range obtained from synmineralization minerals, supporting the idea that a unique fluid source was involved. The absolute timing of this event is unclear, but a plethora of Pb model, K-Ar, and Ar-40/Ar-39 ages between 1440 and 1300 Ma indicate that a significant volume of fluid was mobilized at this time. The deposition of the Roper superbasin from ca. 1492 +/- 4 Ma suggests that these late veins formed from fluids that may have been derived from aquifers in overlying sediments of the Roper superbasin. Clear, buck, and drusy quartz in veins unrelated to any form of Pb-Zn mineralization record the last major fluid event in the Burketown mineral field and form distinct outcrops and ridges in the district. Fluid inclusions in these veins indicate formation from a low-salinity, 300 degrees +/- 80 degrees C fluid. Temperatures approaching 300 degrees C recorded in organic matter adjacent to faults and at sequence boundaries correspond to K-Ar ages spanning 1300 to 1100 Ma, which coincides with regional hydrothermal activity in the northern Lawn Hill platform and the emplacement of the Lakeview Dolerite at the time of assemblage of the Rodinia supercontinent.