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.

Design of solid electrolytes for use with dissimilar gas electrodes

The design of a solid electrolyte that permits the use of dissimilar gas electrodes in an electrochemical cell is presented. It consists of a functionally gradient material with spatial variation in composition. The activity of the conducting ion is fixed at each electrode using different gas species. The system chosen for demonstrating the concept consists of a solid solution between K2CO3 and K2SO4. The composition of the solid solution varies from pure K2CO3 in contact with a CO2 + O2 gas mixture at one electrode to pure K2SO4 exposed to a mixture of SO3 + SO2 + O2 at the other. Two types of composition profiles are studied, one with monotonic variation in composition and the other with extrema. The e.m.f. of the cells is studied as a function of temperature and composition of the gas mixture at each electrode. The results indicate that the e.m.f. is determined primarily by the difference in the chemical potential of potassium at the two electrodes. The diffusion potential caused by ionic concentration gradients in the electrolyte appears to be negligible when the corresponding ionic transport numbers are insignificant. Studies on the response characteristics of the cell based on the gradient electrolyte indicate that the nature of the variation in composition of the electrolyte has only a minor effect on the time evolution of e.m.f. The gradient solid electrolytes have potential application in multielement galvanic sensors at high temperatures.

Gas-phase controlled mass transfer in a foam-bed reactor

Gas-phase controlled absorption of ammonia in foams made of solutions of sulphuric acid has been studied experimentally. Effects of gas-phase concentration of ammonia and type of surfactant on the performance of the foam-bed reactor are investigated. Gas-phase controlled absorption from a spherical bubble is anaylzed using the asymptotic value of Sherwood number (Sh = 6.58), for both negligible as well as significant changes in the volume of the bubble. The experimental data are shown to be in good agreement with the single-stage model of the foam-bed reactor using these asymptotic sub-models, as well as the diffusion-in-sphere analysis available in literature. Influence of effective diffusivity on the time dependence of fractional gas absorption has been found to be unimportant for foam columns with large times of contact. The asymptotic sub-models have been compared and use of the rigid-sphere asymptotic sub-model is recommended for foam columns of practical relevence.

Modelling, simulation, and graphical user interface for industrial gas carburising process

A mathematical model has been developed for the gas carburising (diffusion) process using finite volume method. The computer simulation has been carried out for an industrial gas carburising process. The model's predictions are in good agreement with industrial experimental data and with data collected from the literature. A study of various mass transfer and diffusion coefficients has been carried out in order to suggest which correlations should be used for the gas carburising process. The model has been interfaced in a Windows environment using a graphical user interface. In this way, the model is extremely user friendly. The sensitivity analysis of various parameters such as initial carbon concentration in the specimen, carbon potential of the atmosphere, temperature of the process, etc. has been carried out using the model.

Tribology of ethylene-air diffusion flame soot under dry and lubricated contact conditions

Soot particles are generated in a flame caused by burning ethylene gas. The particles are collected thermophoretically at different locations of the flame. The particles are used to lubricate a steel/steel ball on flat reciprocating sliding contact, as a dry solid lubricant and also as suspended in hexadecane. Reciprocating contact is shown to establish a protective and low friction tribo-film. The friction correlates with the level of graphitic order of the soot, which is highest in the soot extracted from the mid-flame region and is low in the soot extracted from the flame root and flame tip regions. Micro-Raman spectroscopy of the tribo-film shows that the a priori graphitic order, the molecular carbon content of the soot and the graphitization of the film as brought about by tribology distinguish between the frictions of soot extracted from different regions of the flame, and differentiate the friction associated with dry tribology from that recorded under lubricated tribology.

The generalized Onsager model for a binary gas mixture

The Onsager model for the secondary flow field in a high-speed rotating cylinder is extended to incorporate the difference in mass of the two species in a binary gas mixture. The base flow is an isothermal solid-body rotation in which there is a balance between the radial pressure gradient and the centrifugal force density for each species. Explicit expressions for the radial variation of the pressure, mass/mole fractions, and from these the radial variation of the viscosity, thermal conductivity and diffusion coefficient, are derived, and these are used in the computation of the secondary flow. For the secondary flow, the mass, momentum and energy equations in axisymmetric coordinates are expanded in an asymptotic series in a parameter epsilon = (Delta m/m(av)), where Delta m is the difference in the molecular masses of the two species, and the average molecular mass m(av) is defined as m(av) = (rho(w1)m(1) + rho(w2)m(2))/rho(w), where rho(w1) and rho(w2) are the mass densities of the two species at the wall, and rho(w) = rho(w1) + rho(w2). The equation for the master potential and the boundary conditions are derived correct to O(epsilon(2)). The leading-order equation for the master potential contains a self-adjoint sixth-order operator in the radial direction, which is different from the generalized Onsager model (Pradhan & Kumaran, J. Fluid Mech., vol. 686, 2011, pp. 109-159), since the species mass difference is included in the computation of the density, viscosity and thermal conductivity in the base state. This is solved, subject to boundary conditions, to obtain the leading approximation for the secondary flow, followed by a solution of the diffusion equation for the leading correction to the species mole fractions. The O(epsilon) and O(epsilon(2)) equations contain inhomogeneous terms that depend on the lower-order solutions, and these are solved in a hierarchical manner to obtain the O(epsilon) and O(epsilon(2)) corrections to the master potential. A similar hierarchical procedure is used for the Carrier-Maslen model for the end-cap secondary flow. The results of the Onsager hierarchy, up to O(epsilon(2)), are compared with the results of direct simulation Monte Carlo simulations for a binary hard-sphere gas mixture for secondary flow due to a wall temperature gradient, inflow/outflow of gas along the axis, as well as mass and momentum sources in the flow. There is excellent agreement between the solutions for the secondary flow correct to O(epsilon(2)) and the simulations, to within 15 %, even at a Reynolds number as low as 100, and length/diameter ratio as low as 2, for a low stratification parameter A of 0.707, and when the secondary flow velocity is as high as 0.2 times the maximum base flow velocity, and the ratio 2 Delta m/(m(1) + m(2)) is as high as 0.5. Here, the Reynolds number Re = rho(w)Omega R-2/mu, the stratification parameter A = root m Omega R-2(2)/(2k(B)T), R and Omega are the cylinder radius and angular velocity, m is the molecular mass, rho(w) is the wall density, mu is the viscosity and T is the temperature. The leading-order solutions do capture the qualitative trends, but are not in quantitative agreement.

Simultaneous measurement of temperature, density and velocity in gas-flows by modulated photoluminescence

A novel possibility to determine the temperature, density and velocity simultaneously in gas flows by measuring the average value, amplitude of modulation and phase shift of the photoluminescence excited by a temporally or spatially modulated light source is investigated. Time-dependent equations taking the flow, diffusion, excitation and decay into account are solved analytically. Different experimental arrangements are proposed. Measurements of velocity with two components, and temporal and spatial resolutions in the measurements are investigated. Numerical examples are given for N z with biacetyl as the seed gas. Practical considerations for the measurements and the relation between this method and some existing methods of lifetime measurement are discussed.

Gain saturation in gas flow and chemical lasers

A theoretical model for gain saturation in gas flow and chemical lasers is presented. The theory is applicable to all possible numerical values of τ/τc, where τ is the characteristie flow time for the flowing gas to move across the laser action region and τc is the characteristic collision relaxation time. The saturation effects of the convection and the "source flow" of the inverted population are revealed. A general relation of gain coefficient and some new gain saturation laws are obtained. For the special case of τ/τc1, the present theoretical results agree with the experimental results on the "anomalous" saturation phenomena in the supersonic diffusion HF chemical laser determined recently by Gross and Coffer[8]. The theory also agrees with the measured results of saturation intensity varying with τ/τc in gas flow CO2 lasers[7]. For the special case of τ/τc1, the present theory is consistent with both the standard theory[1] for gas lasers where the gas has no macroscopic motion and the known gain saturation theory[2-5] for gas flow and chemical lasers.

Diffusion in glassy polymers

Fluid diffusion in glassy polymers proceeds in ways that are not explained by the standard diffusion model. Although the reasons for the anomalous effects are not known, much of the observed behavior is attributed to the long times that polymers below their glass transition temperature take to adjust to changes in their condition. The slow internal relaxations of the polymer chains ensure that the material properties are history-dependent, and also allow both local inhomogeneities and differential swelling to occur. Two models are developed in this thesis with the intent of accounting for these effects in the diffusion process.

In Part I, a model is developed to account for both the history dependence of the glassy polymer, and the dual sorption which occurs when gas molecules are immobilized by the local heterogeneities. A preliminary study of a special case of this model is conducted, showing the existence of travelling wave solutions and using perturbation techniques to investigate the effect of generalized diffusion mechanisms on their form. An integral averaging method is used to estimate the penetrant front position.

In Part II, a model is developed for particle diffusion along with displacements in isotropic viscoelastic materials. The nonlinear dependence of the materials on the fluid concentration is taken into account, while pure displacements are assumed to remain in the range of linear viscoelasticity. A fairly general model is obtained for three-dimensional irrotational movements, with the development of the model being based on the assumptions of irreversible thermodynamics. With the help of some dimensional analysis, this model is simplified to a version which is proposed to be studied for Case II behavior.

New insights into the formation and modification of carbonate-bearing minerals and methane gas in geological systems using multiply substituted isotopologues

This thesis describes the use of multiply-substituted stable isotopologues of carbonate minerals and methane gas to better understand how these environmentally significant minerals and gases form and are modified throughout their geological histories. Stable isotopes have a long tradition in earth science as a tool for providing quantitative constraints on how molecules, in or on the earth, formed in both the present and past. Nearly all studies, until recently, have only measured the bulk concentrations of stable isotopes in a phase or species. However, the abundance of various isotopologues within a phase, for example the concentration of isotopologues with multiple rare isotopes (multiply substituted or 'clumped' isotopologues) also carries potentially useful information. Specifically, the abundances of clumped isotopologues in an equilibrated system are a function of temperature and thus knowledge of their abundances can be used to calculate a sample’s formation temperature. In this thesis, measurements of clumped isotopologues are made on both carbonate-bearing minerals and methane gas in order to better constrain the environmental and geological histories of various samples.

Clumped-isotope-based measurements of ancient carbonate-bearing minerals, including apatites, have opened up paleotemperature reconstructions to a variety of systems and time periods. However, a critical issue when using clumped-isotope based measurements to reconstruct ancient mineral formation temperatures is whether the samples being measured have faithfully recorded their original internal isotopic distributions. These original distributions can be altered, for example, by diffusion of atoms in the mineral lattice or through diagenetic reactions. Understanding these processes quantitatively is critical for the use of clumped isotopes to reconstruct past temperatures, quantify diagenesis, and calculate time-temperature burial histories of carbonate minerals. In order to help orient this part of the thesis, Chapter 2 provides a broad overview and history of clumped-isotope based measurements in carbonate minerals.

In Chapter 3, the effects of elevated temperatures on a sample’s clumped-isotope composition are probed in both natural and experimental apatites (which contain structural carbonate groups) and calcites. A quantitative model is created that is calibrated by the experiments and consistent with the natural samples. The model allows for calculations of the change in a sample’s clumped isotope abundances as a function of any time-temperature history.

In Chapter 4, the effects of diagenesis on the stable isotopic compositions of apatites are explored on samples from a variety of sedimentary phosphorite deposits. Clumped isotope temperatures and bulk isotopic measurements from carbonate and phosphate groups are compared for all samples. These results demonstrate that samples have experienced isotopic exchange of oxygen atoms in both the carbonate and phosphate groups. A kinetic model is developed that allows for the calculation of the amount of diagenesis each sample has experienced and yields insight into the physical and chemical processes of diagenesis.

The thesis then switches gear and turns its attention to clumped isotope measurements of methane. Methane is critical greenhouse gas, energy resource, and microbial metabolic product and substrate. Despite its importance both environmentally and economically, much about methane’s formational mechanisms and the relative sources of methane to various environments remains poorly constrained. In order to add new constraints to our understanding of the formation of methane in nature, I describe the development and application of methane clumped isotope measurements to environmental deposits of methane. To help orient the reader, a brief overview of the formation of methane in both high and low temperature settings is given in Chapter 5.

In Chapter 6, a method for the measurement of methane clumped isotopologues via mass spectrometry is described. This chapter demonstrates that the measurement is precise and accurate. Additionally, the measurement is calibrated experimentally such that measurements of methane clumped isotope abundances can be converted into equivalent formational temperatures. This study represents the first time that methane clumped isotope abundances have been measured at useful precisions.

In Chapter 7, the methane clumped isotope method is applied to natural samples from a variety of settings. These settings include thermogenic gases formed and reservoired in shales, migrated thermogenic gases, biogenic gases, mixed biogenic and thermogenic gas deposits, and experimentally generated gases. In all cases, calculated clumped isotope temperatures make geological sense as formation temperatures or mixtures of high and low temperature gases. Based on these observations, we propose that the clumped isotope temperature of an unmixed gas represents its formation temperature — this was neither an obvious nor expected result and has important implications for how methane forms in nature. Additionally, these results demonstrate that methane-clumped isotope compositions provided valuable additional constraints to studying natural methane deposits.

Test particle motion in a Lorentz gas

This is a two-part thesis concerning the motion of a test particle in a bath. In part one we use an expansion of the operator PLe^it(1-P)LLP to shape the Zwanzig equation into a generalized Fokker-Planck equation which involves a diffusion tensor depending on the test particle's momentum and the time.

In part two the resultant equation is studied in some detail for the case of test particle motion in a weakly coupled Lorentz Gas. The diffusion tensor for this system is considered. Some of its properties are calculated; it is computed explicitly for the case of a Gaussian potential of interaction.

The equation for the test particle distribution function can be put into the form of an inhomogeneous Schroedinger equation. The term corresponding to the potential energy in the Schroedinger equation is considered. Its structure is studied, and some of its simplest features are used to find the Green's function in the limiting situations of low density and long time.

Diffusion and chemical reaction in the porous structures of solid oxide fuel cells

The Rolls-Royce Integrated-Planar Solid Oxide Fuel Cell (IP-SOFC) consists of ceramic modules which have electrochemical cells printed on the outer surfaces. The cathodes are the outermost layer of each cell and are supplied with oxygen from air flowing over the outside of the module. The anodes are in direct contact with the ceramic structure and are supplied with fuel from internal gas channels. Natural gas is reformed into hydrogen for use by the fuel cells in a separate reformer module of similar design except that the fuel cells are replaced by a reforming catalyst layer. The performance of the modules is intrinsically linked to the behaviour of the gas flows within their porous structures. Because the porous layers are very thin, a one-dimensional flow model provides a good representation of the flow property variations between fuel channel and fuel cell or reforming catalyst. The multi-component convective-diffusive flows are simulated using a new theory of flow in porous material, the Cylindrical Pore Interpolation Model. The effects of the catalysed methane reforming and water-gas shift chemical reactions are also considered using appropriate kinetic models. It is found that the shift reaction, which is catalysed by the anode material, has certain beneficial effects on the fuel cell module performance. In the reformer module it was found that the flow resistance of the porous support structure makes it difficult to sustain a high methane conversion rate. Although the analysis is based on IP-SOFC geometry, the modelling approach and general conclusions are applicable to other types of SOFC.

Validation of a lattice Boltzmann model for gas-solid reactions with experiments

A lattice Boltzmann method is used to model gas-solid reactions where the composition of both the gas and solid phase changes with time, while the boundary between phases remains fixed. The flow of the bulk gas phase is treated using a multiple relaxation time MRT D3Q19 model; the dilute reactant is treated as a passive scalar using a single relaxation time BGK D3Q7 model with distinct inter- and intraparticle diffusivities. A first-order reaction is incorporated by modifying the method of Sullivan et al. [13] to include the conversion of a solid reactant. The detailed computational model is able to capture the multiscale physics encountered in reactor systems. Specifically, the model reproduced steady state analytical solutions for the reaction of a porous catalyst sphere (pore scale) and empirical solutions for mass transfer to the surface of a sphere at Re=10 (particle scale). Excellent quantitative agreement between the model and experiments for the transient reduction of a single, porous sphere of Fe 2O 3 to Fe 3O 4 in CO at 1023K and 10 5Pa is demonstrated. Model solutions for the reduction of a packed bed of Fe 2O 3 (reactor scale) at identical conditions approached those of experiments after 25 s, but required prohibitively long processor times. The presented lattice Boltzmann model resolved successfully mass transport at the pore, particle and reactor scales and highlights the relevance of LB methods for modelling convection, diffusion and reaction physics. © 2012 Elsevier Inc.

Steady-state modelling of the universal exhaust gas oxygen (UEGO) sensor

The universal exhaust gas oxygen (UEGO) sensor is a well-established device which was developed for the measurement of relative air fuel ratio in internal combustion engines. There is, however, little information available which allows for the prediction of the UEGO's behaviour when exposed to arbitrary gas mixtures, pressures and temperatures. Here we present a steady-state model for the sensor, based on a solution of the Stefan-Maxwell equation, and which includes a momentum balance. The response of the sensor is dominated by a diffusion barrier, which controls the rate of diffusion of gas species between the exhaust and a cavity. Determination of the diffusion barrier characteristics, especially the mean pore size, porosity and tortuosity, is essential for the purposes of modelling, and a measurement technique based on identification of the sensor pressure giving zero temperature sensitivity is shown to be a convenient method of achieving this. The model, suitably calibrated, is shown to make good predictions of sensor behaviour for large variations of pressure, temperature and gas composition. © 2012 IOP Publishing Ltd.