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.

Application of composition gradient solid electrolytes in sensors and thermodynamic measurements

New compos~tiong radient solid electrolytes are developed which have application in high temperature solid state galvanic sensors and provide a new tool for thermodynamic measurements. The electrolyte consists oi a solid solution between two ionic conductors with a common mobile ion and spatial variation in composition of otber coxup nents. Incorporation of the composite electrolyte in sensors permits the use oi dissimilar gas electrodes. It is demonsuated, both experimentall y and theoretically, that the composition gradient of the relativeiy immobile species does not give rise to a diffusion potential.The emi of a cell is determined by the activity of the mobile species at the two eiectrodes. The thermodynamic properties of solid solutions can be measured using the gradient solid electrolyte. The experimental stuay is based on model systems A?(COj)x(S04)l-x (A=Na,K),where S \.aria across the electrolyte. The functionally gradient solid electrolytes used for activity measurements consist of pure carbonate at one ena and the solid solution under stuav at the other. The identical vaiues of activity, obtained h m t hree different modes of operation of the ceil. indicate unit transport number for the ddi metal ion in the graciient electrolyte. Tlle activities in the solid solutions exhibit moderate positive deviations from Raoult 's law.

Phase relations in the system Sr–Ir–O and thermodynamic measurements on SrIrO3, Sr2IrO4 and Sr4IrO6 using solid-state cells with buffer electrodes

The standard Gibbs energies of formation of SrIrO3, Sr2IrO4 and Sr4IrO6 have been determined in the temperature range from 975 to 1400 K using solid-state cells with (Y2O3) ZrO2 as the electrolyte and pure oxygen gas at a pressure of 0.1 MPa as the reference electrode. For the design of appropriate working electrodes, phase relations in the ternary system Sr–Ir–O were investigated at 1350 K. The only stable oxide detected along the binary Ir–O was IrO2. Three ternary oxides, SrIrO3, Sr2IrO4 and Sr4IrO6, compositions of which fall on the join SrO–IrO2, were found to be stable. Each of the oxides coexisted with pure metal Ir. Therefore, three working electrodes were prepared consisting of mixtures of Ir+SrO+Sr4IrO6, Ir+Sr4IrO6+Sr2IrO4, and Ir+Sr2IrO4+SrIrO3. These mixtures unambiguously define unique oxygen chemical potentials under isothermal and isobaric conditions. Used for the measurements was a novel apparatus, in which a buffer electrode was introduced between reference and working electrodes to absorb the electrochemical flux of oxygen through the solid electrolyte. The buffer electrode prevented polarization of the measuring electrode and ensured accurate data. The standard Gibbs energies of formation of the compounds, obtained from the emf of the cells, can be represented by the following equations: View the MathML sourcem View the MathML source View the MathML source where Δf (ox)Go represents the standard Gibbs energy of formation of the ternary compound from its component binary oxides SrO and IrO2. Based on the thermodynamic information, chemical potential diagrams for the system Sr–Ir–O were developed.

Solid-state cells with buffer electrodes for accurate thermodynamic measurements: system Nd-Ir-O

A new design for the solid-state cell incorporating a buffer electrode for high-temperature thermodynamic measurements is presented. The function of the buffer electrode, placed between the reference and working electrodes, is to absorb the electrochemical flux of the mobile species through the solid electrolyte caused by trace electronic conductivity. The buffer electrode prevents polarization of the measuring electrode and ensures accurate data. The application of this novel design and its advantages are demonstrated by measurement of the standard Gibbs energies of formation of Nd6Ir2O13 (low-temperature form) and Nd2Ir2O7 in the temperature range from 975 to 1450 K. Yttria-stabilized zirconia is used as the solid electrolyte and pure oxygen gas at a pressure of 0.1 MPa as the reference electrode. For the design of appropriate working electrodes, phase relations in the ternary system NdIrO were investigated at 1350 K. The two ternary oxides, Nd6Ir2O13 and Nd2Ir2O7, compositions of which fall on the join Nd2O3IrO2, were found to coexist with pure metal Ir. Therefore, two working electrodes were prepared consisting of mixtures of Ir+Nd2O3+Nd6Ir2O13 and Ir+Nd6Ir2O13+ Nd2Ir2O7. These mixtures unambiguously define unique oxygen chemical potentials under isothermal and isobaric conditions. The standard Gibbs energies of formation (ΔG°f (ox)) of the compounds from their component binary oxides Nd2O3 and IrO2, obtained from the emf of the cells, can be represented by the equations:View the MathML source View the MathML source Based on the thermodynamic information, chemical potential diagrams for the system NdIrO are developed.

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.

Electrical Breakdown of Gases Between Coaxial Cylindrical Electrodes in Crossed Electric and Magnetic Fields

Sparking potentials have been measured in nitrogen and dry air between coaxial cylindrical electrodes for values of n = R2/R1 = approximately 1 to 30 (R1 = inner electrode radius, R2 = outer electrode radius) in the presence of crossed uniform magnetic fields. The magnetic flux density was varied from 0 to 3000 Gauss. It has been shown that the minimum sparking potentials in the presence of the crossed magnetic field can be evaluated on the basis of the equivalent pressure concept when the secondary ionization coefficient does not vary appreciably with B/p (B = magnetic flux density, p = gas pressure). The values of secondary ionization coefficients �¿B in nitrogen in crossed fields calculated from measured values of sparking potentials and Townsend ionization coefficients taken from the literature, have been reported. The calculated values of collision frequencies in nitrogen from minimum sparking potentials in crossed fields are found to increase with increasing B/p at constant E/pe (pe = equivalent pressure). Studies on the similarity relationship in crossed fields has shown that the similarity theorem is obeyed in dry air for both polarities of the central electrode in crossed fields.

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.

Part I. Influence of membrane morphology on ion-exchange and charge propagation in composite polyelectrolyte electrode coatings. Part II. Dynamic surface tension measurements of polarization relaxation at mercury pool electrodes by the method of Wilhelmy.

Part I. Novel composite polyelectrolyte materials were developed that exhibit desirable charge propagation and ion-retention properties. The morphology of electrode coatings cast from these materials was shown to be more important for its electrochemical behavior than its chemical composition.

Part II. The Wilhelmy plate technique for measuring dynamic surface tension was extended to electrified liquid-liquid interphases. The dynamical response of the aqueous NaF-mercury electrified interphase was examined by concomitant measurement of surface tension, current, and applied electrostatic potential. Observations of the surface tension response to linear sweep voltammetry and to step function perturbations in the applied electrostatic potential (e.g., chronotensiometry) provided strong evidence that relaxation processes proceed for time-periods that are at least an order of magnitude longer than the time periods necessary to establish diffusion equilibrium. The dynamical response of the surface tension is analyzed within the context of non-equilibrium thermodynamics and a kinetic model that requires three simultaneous first order processes.

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.