Digital simulation of galvanostatic current-potential data for gas-diffusion electrodes and estimation of electrodekinetic parameters

A computerized non-linear-least-squares regression procedure to analyse the galvanostatic current-potential data for kinetically hindered reactions on porous gas-diffusion electrodes is reported. The simulated data fit well with the corresponding measured values. The analytical estimates of electrode-kinetic parameters and uncompensated resistance are found to be in good agreement with their respective values obtained from Tafel plots and the current-interrupter method. The procedure circumvents the need to collect the data in the limiting-current region where the polarization values are usually prone to errors. The polarization data for two typical cases, namely, methanol oxidation on a carbon-supported platinum-tin electrode and oxygen reduction on a Nafion-coated platinized carbon electrode, are successfully analysed.

Effect of diffusion-layer porosity on the performance of polymer electrolyte fuel cells

The gas-diffusion layer (GDL) influences the performance of electrodes employed with polymer electrolyte fuel cells (PEFCs). A simple and effective method for incorporating a porous structure in the electrode GDL using sucrose as the pore former is reported. Optimal (50 w/o) incorporation of a pore former in the electrode GDL facilitates the access of the gaseous reactants to the catalyst sites and improves the fuel cell performance. Data obtained from permeability and porosity measurements, single-cell performance, and impedance spectroscopy suggest that an optimal porosity helps mitigating mass-polarization losses in the fuel cell resulting in a substantially enhanced performance.

Experimental Analysis of Water Vapor Diffusion Through Porous Membranes in a Proton Exchange Membrane Fuel Cell

We report the diffusion characteristics of water vapor through two different porous media, viz., membrane electrode assembly (MEA) and gas diffusion layer (GDL) in a nonoperational fuel cell. Tunable diode laser absorption spectroscopy (TDLAS) was employed for measuring water vapor concentration in the test channel. Effects of the membrane pore size and the inlet humidity on the water vapor transport are quantified through mass flux and diffusion coefficient. Water vapor transport rate is found to be higher for GDL than for MEA. The flexibility and wide range of application of TDLAS in a fuel cell setup is demonstrated through experiments with a stagnant flow field on the dry side.

An electrochemical method for the synthesis of few layer graphene sheets for high temperature applications

We demonstrate an electrochemical technique for the large scale synthesis of high quality few layer graphene sheets (FLGS) directly from graphite using oxalic acid (a weak acid) as the electrolyte. One of the interesting observations is that the FLGS are stable at least up to 800 degrees C and hence have potential application in solid oxide fuel cells as a gas diffusion layer.

An Analysis of a Boundary Layer Diffusion Flame Over a Porous Flat Plate in a Confined Flow

A numerical analysis of the gas dynamic structure of a two-dimensional laminar boundary layer diffusion flame over a porous flat plate in a confined flow is made on the basis of the familiar boundary layer and flame sheet approximations neglecting buoyancy effects. The governing equations of aerothermochemistry with the appropriate boundary conditions are solved using the Patankar-Spalding method. The analysis predicts the flame shape, profiles of temperature, concentrations of variousspecies, and the density of the mixture across the boundary layer. In addition, it also predicts the pressure gradient in the flow direction arising from the confinement ofthe flow and the consequent velocity overshoot near the flame surface. The results of thecomputation performed for an n-pentane-air system are compared with experimental data andthe agreement is found to be satisfactory.

The ratio of diffusion coefficient to mobility for slow electrons in dry air

Using Huxley's solution of the diffusion equation for electron-attaching gases, the ratio of diffusion coefficient D to mobility μ for electrons in dry air was measured over the range 3·06 × 10-17 gas density (cm -3). The results agreed with those of earlier workers who used the diffusion equation for non-attaching gases.

Synthesis of FeCu Nanopowder by Levitational Gas Condensation Process

Condensation from the vapor state is an important technique for the preparation of nanopowders. Levitational gas condensation is one such technique that has a unique ability of attaining steady state. Here, we present the results of applying this technique to an iron-copper alloy (96Fe-4Cu). A qualitative model of the process is proposed to understand the process and the characteristics of resultant powder. A phase diagram of the alloy system in the liquid-vapor region was calculated to help understand the course of condensation, especially partitioning and coring during processing. The phase diagram could not explain coring in view of the simultaneous occurrence of solidification and the fast homogenization through diffusion in the nanoparticles; however, it could predict the very low levels of copper observed in the levitated drop. The enrichment of copper observed near the surface of the powder was considered to be a manifestation of the lower surface energy of copper compared with that of iron. Heat transfer calculations indicated that most condensed particles can undergo solidification even when they are still in the proximity of the levitated drop. It helped us to predict the temperature and the cooling rate of the powder particles as they move away from the levitated drop. The particles formed by the process seem to be single domain, single crystals that are magnetic in nature. They, thus, can agglomerate by forming a chain-like structure, which manifests as a three-dimensional network enclosing a large unoccupied space, as noticed in scanning electron microscopy and transmission electron microscopy studies. This also explains the observed low packing density of the nanopowders.

Gradient solid-electrolyte for use with dissimilar gas electrodes

The EMF of a solid-state cell, incorporating a composite solid-electrolyte with gradual variation in composition, and dissimilar gas electrodes, has been studied as a function of temperature and partial pressures at the electrodes. The cell with the configuration: Pt, CO2' + O2' parallel-to Na2CO3\Na(SO4)x(CO3)1-x\Na2SO4 parallel-to SO3'' + SO2'' + O2'', Pt x=0 x=1 was investigated in the temperature range 973 to 1079 K. The solid-electrolyte surface exposed to SO3 + SO2 + O2 gas mixture was doped-Na2SO4, whereas the CO2 + O2 gas mixture was in contact with pure Na2CO3. The composition of the solid solution between the carbonate and sulfate, with hexagonal structure, was varied gradually between the boundary values. It has been found that the EMF of the cell is close to that calculated from thermodynamic data, assuming unit transport number for Na+ ions. The gradient in the concentration of sulfate and carbonate ions in the electrolyte does not give rise to a significant diffusion potential.

Crossover dynamics at large metastability in gas-liquid nucleation

We have developed an alternate description of dynamics of nucleation in terms of an extended set of order parameters. The order parameters consist of an ordered set of kth largest clusters, ordered such that k = 1 is the largest cluster in the system, k = 2 is the second largest cluster, and so on. We have derived an analytic expression for the free energy for the kth largest cluster, which is in excellent agreement with the simulated results. At large supersaturation, the free energy barrier for the growth of the kth largest cluster disappears and the nucleation becomes barrierless. The major success of this extended theoretical formalism is that it can clearly explain the observed change in mechanism at large metastability P. Bhimalapuram et al., Phys. Rev. Lett. 98, 206104 (2007)] and the associated dynamical crossover. The classical nucleation theory cannot explain this crossover. The crossover from activated to barrierless nucleation is found to occur at a supersaturation where multiple clusters cross the critical size. We attribute the crossover as the onset of the kinetic spinodal. We have derived an expression for the rate of nucleation in the barrierless regime by modeling growth as diffusion on the free energy surface of the largest cluster. The model reproduces the slower increase in the rate of growth as a function of supersaturation, as observed in experiments.

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.