A numerical model of a liquid-feed solid polymer electrolyte DMFC and its experimental validation

A one-dimensional, biphasic, multicomponent steady-state model based on phenomenological transport equations for the catalyst layer, diffusion layer, and polymeric electrolyte membrane has been developed for a liquid-feed solid polymer electrolyte direct methanol fuel cell (SPE- DMFC). The model employs three important requisites: (i) implementation of analytical treatment of nonlinear terms to obtain a faster numerical solution as also to render the iterative scheme easier to converge, (ii) an appropriate description of two-phase transport phenomena in the diffusive region of the cell to account for flooding and water condensation/evaporation effects, and (iii) treatment of polarization effects due to methanol crossover. An improved numerical solution has been achieved by coupling analytical integration of kinetics and transport equations in the reaction layer, which explicitly include the effect of concentration and pressure gradient on cell polarization within the bulk catalyst layer. In particular, the integrated kinetic treatment explicitly accounts for the nonhomogeneous porous structure of the catalyst layer and the diffusion of reactants within and between the pores in the cathode. At the anode, the analytical integration of electrode kinetics has been obtained within the assumption of macrohomogeneous electrode porous structure, because methanol transport in a liquid-feed SPE- DMFC is essentially a single-phase process because of the high miscibility of methanol with water and its higher concentration in relation to gaseous reactants. A simple empirical model accounts for the effect of capillary forces on liquid-phase saturation in the diffusion layer. Consequently, diffusive and convective flow equations, comprising Nernst-Plank relation for solutes, Darcy law for liquid water, and Stefan-Maxwell equation for gaseous species, have been modified to include the capillary flow contribution to transport. To understand fully the role of model parameters in simulating the performance of the DMCF, we have carried out its parametric study. An experimental validation of model has also been carried out. (C) 2003 The Electrochemical Society.

Study of gas-fines in decreasing gas velocity in a packed bed

It has been well recognized now that the blast furnace aerodynamics can be represented more accurately under the decreasing gas velocity condition. Therefore, gas-fines study has been carried out in a packed bed under the decreasing gas velocity condition. Gas and fines flow equations have been developed and solved, for two-dimensional case using finite volume method. To take into account the turbulence, k-e turbulent flow model has also been developed in two-dimension. The model's predictions have been validated against the published experimental data for the increasing gas velocity case, as no experimental data are available in open literature for the decreasing gas velocity. This study shows the difference in the results for increasing and decreasing gas velocity cases under various conditions which have been reported here. Implication of the results to the blast furnace condition has also been discussed.

Modeling of nonreacting and reacting turbulent spray jets using a fully stochastic separated flow approach

Numerical modeling of several turbulent nonreacting and reacting spray jets is carried out using a fully stochastic separated flow (FSSF) approach. As is widely used, the carrier-phase is considered in an Eulerian framework, while the dispersed phase is tracked in a Lagrangian framework following the stochastic separated flow (SSF) model. Various interactions between the two phases are taken into account by means of two-way coupling. Spray evaporation is described using a thermal model with an infinite conductivity in the liquid phase. The gas-phase turbulence terms are closed using the k-epsilon model. A novel mixture fraction based approach is used to stochastically model the fluctuating temperature and composition in the gas phase and these are then used to refine the estimates of the heat and mass transfer rates between the droplets and the surrounding gas-phase. In classical SSF (CSSF) methods, stochastic fluctuations of only the gas-phase velocity are modeled. Successful implementation of the FSSF approach to turbulent nonreacting and reacting spray jets is demonstrated. Results are compared against experimental measurements as well as with predictions using the CSSF approach for both nonreacting and reacting spray jets. The FSSF approach shows little difference from the CSSF predictions for nonreacting spray jets but differences are significant for reacting spray jets. In general, the FSSF approach gives good predictions of the flame length and structure but further improvements in modeling may be needed to improve the accuracy of some details of the Predictions. (C) 2011 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

CFD Based Design Optimizations of the Diffuser of a Gas Dynamically Driver Laser Cavity

Based on the an earlier CFD analysis of the performance of the gas-dynamically controlled laser cavity [1]it was found that there is possibility of optimizing the geometry of the diffuser that can bring about reductions in both size and cost of the system by examining the critical dimensional requirements of the diffuser. Consequently,an extensive CFD analysis has been carried out for a range of diffuser configurations by simulating the supersonic flow through the arrangement including the laser cavity driven by a bank of converging – diverging nozzles and the diffuser. The numerical investigations with 3D-RANS code are carried out to capture the flow patterns through diffusers past the cavity that has multiple supersonic jet interactions with shocks leading to complex flow pattern. Varying length of the diffuser plates is made to be the basic parameter of the study. The analysis reveals that the pressure recovery pattern during the flow through the diffuser from the simulation, being critical for the performance of the laser device shows its dependence on the diffuser length is weaker beyond a critical lower limit and this evaluation of this limit would provide a design guideline for a more efficient system configuration.The observation based on the parametric study shows that the pressure recovery transients in the near vicinity of the cavity is not affected for the reduction in the length of the diffuser plates up to its 10% of the initial size, indicating the design in the first configuration that was tested experimentally has a large factor of margin. The flow stability in the laser cavity is found to be unaffected since a strong and stable shock is located at the leading edge of the diffuser plates while the downstream shock and flow patterns are changed, as one would expect. Results of the study for the different lengths of diffusers in the range of 10% to its full length are presented, keeping the experimentally tested configuration used in the earlier study [1] as the reference length. The conclusions drawn from the analysis is found to be of significance since it provides new design considerations based on the understanding of the intricacies of the flow, allowing for a hardware optimization that can lead to substantial size reduction of the device with no loss of performance.

Dynamic contact angle beating from drops impacting onto solid surfaces exhibiting anisotropic wetting

This paper reports an experimental investigation of low Weber number water drops impacting onto solid surfaces exhibiting anisotropic wetting. The wetting anisotropy is created by patterning the solid surfaces with unidirectional parallel grooves. Temporal measurements of impacting drop parameters such as drop base contact diameter, apparent contact angle of drop, and drop height at the center are obtained from high-speed video recordings of drop impacts. The study shows that the impact of low Weber number water drops on the grooved surface exhibits beating phenomenon in the temporal variations of the dynamic contact angle anisotropy and drop height at the center of the impacting drop. It is observed that the beating phenomenon of impacting drop parameters is caused by the frequency difference between the dynamic contact angle oscillations of impacting drop liquid oriented perpendicular and parallel to the direction of grooves on the grooved surface. The primary trigger for the phenomenon is the existence of non-axisymmetric drop flow on the grooved surface featuring pinned and free motions of drop liquid in the directions perpendicular and parallel to the grooves, respectively. The beat frequency is almost independent of the impact drop Weber number. Further experimental measurements with solid surfaces of different groove textures show that the grooved surface with larger wetting anisotropy may be expected to show a dominant beating phenomenon. The phenomenon is gradually damped out with time and is fully unrecognizable at higher drop impact Weber numbers. (C) 2011 Elsevier B.V. All rights reserved.

Ru(4+) ion in CeO(2) (Ce(0.95)Ru(0.05)O(2-delta)): A non-deactivating, non-platinum catalyst for water gas shift reaction

Hydrogen is a clean energy carrier and highest energy density fuel. Water gas shift (WGS) reaction is an important reaction to generate hydrogen from steam reforming of CO. A new WGS catalyst, Ce(1-x)Ru(x)O(2-delta) (0 <= x <= 0.1) was prepared by hydrothermal method using melamine as a complexing agent. The Catalyst does not require any pre-treatment. Among the several compositions prepared and tested, Ce(0.95)Ru(0.05)O(2-delta) (5% Ru(4+) ion substituted in CeO(2)) showed very high WGS activity in terms of high conversion rate (20.5 mu mol.g(-1).s(-1) at 275 degrees C) and low activation energy (12.1 kcal/mol). Over 99% conversion of CO to CO(2) by H(2)O is observed with 100% H(2) selectivity at >= 275 degrees C. In presence of externally fed CO(2) and H(2) also, complete conversion of CO to CO(2) was observed with 100% H(2) selectivity in the temperature range of 305-385 degrees C. Catalyst does not deactivate in long duration on/off WGS reaction cycle due to absence of surface carbon and carbonate formation and sintering of Ru. Due to highly acidic nature of Ru(4+) ion, surface carbonate formation is also inhibited. Sintering of noble metal (Ru) is avoided in this catalyst because Ru remains in Ru(4+) ionic state in the Ce(1-x)Ru(x)O(2-delta) catalyst.

Ce(0.67)Fe(0.33)O(2-delta) and Ce(0.65)Fe(0.33)Pt(0.02)O(2-delta): New water gas shift (WGS) catalysts

Ce(0.65)Fe(0.33)Pt(0.02)O(2-delta) and Ce(0.67)Fe(0.33)O(2-delta) have been synthesized by a new low temperature sonochemical method using diethylenetriamine as a complexing agent. Due to the substitution of Fe and Pt ions in CeO(2), lattice oxygen is activated in Ce(0.67)Fe(0.33)O(2-delta) and Ce(0.65)Fe(0.33)Pt(0.02)O(2-delta). Hydrogen uptake studies show strong reduction peaks at 125 C in Ce(0.65)Fe(0.33)Pt(0.02)O(2-delta) against a hydrogen uptake peak at 420 degrees C in Ce(0.67)Fe(0.33)O(2-delta). Fe substituted ceria, Ce(0.67)Fe(0.33)O(2-delta) itself acts as a catalyst for CO oxidation and water gas shift (WGS) reactions at moderate temperatures. The rate of CO conversion in WGS with Pt free Ce(0.65)Fe(0.33)O(2-delta) is 2.8 mu mol g(-1) s(-1) at 450 C and with Pt substituted Ce(0.65)Fe(0.33)Pt(0.02)O(2-delta) is 4.05 mu mol g(-1) s(-1) at 275 degrees C. Due to the synergistic interaction of the Pt ion with Ce and Fe ions in Ce(0.65)Fe(0.33)Pt(0.02)O(2-delta), the catalyst showed much higher activity for CO oxidation and WGS reactions compared to Ce(0.67)Fe(0.33)O(2-delta). A reverse WGS reaction does not occur over Ce(0.65)Fe(0.33)Pt(0.02)O(2-delta). The catalyst also does not deactivate even when operated for a long time. Nearly 100% conversion of CO to CO(2) with 100% H(2) selectivity is observed in WGS reactions even up to 550 degrees C. (C) 2011 Elsevier B.V. All rights reserved.

Chemical potential of germanium in its solid solutions with Cu, Ag and Au

Thin foils of copper, silver and gold were equilibrated with tetragonal GeO2 under controlled View the MathML source gas streams at 1000 K. The equilibrium concentration of germanium in the foils was determined by the X-ray fluorescence technique. The standard free energy of formation of tetragonal GeO2 was measured by a solid oxide galvanic cell. The chemical potential of germanium calculated from the experimental data and the free energies of formation of carbon monoxide and carbon dioxide was found to decrease in the sequence Ag + Ge > Au + Ge > Cu + Ge. The more negative value for the chemical potential of germanium in solid copper, compared to that in solid gold, cannot be explained in terms of the strain energy factor, electro-negativity differences or the vaporization energies of the solvent, and suggests that the d band and its hybridization with s electrons are an important factor in determining the absolute values for the chemical potential in dilute solutions. However, the variation of the chemical potential with solute concentration can be correlated to the concentration of s and p electrons in the outer shell.

Ion exchange equilibria between (Mn,Co)O solid solution and (Mn,Co)Cr2O4 and (Mn,Co)Al2O4 spinel solid solutions at 1100oC

The compositions of the (Mn,Co)O solid solution with rock salt structure in equilibrium with (Mn,Co)Cr2O4 and (Mn,Co)Al2O4 spinel solid solutions have been determined by X-ray diffraction measurements at 1100° C and an oxygen partial pressure of 10–10 atm. The ion exchange equilibria are quantitatively analysed, using values for activities in the (Mn,Co)O solid solution available in the literature, in order to obtain activities in the spinel solid solutions. The MnAl2O4-CoAl2O4 solid solution exhibits negative deviations from Raoult's law, consistent with the estimated cation disorder in the solid solution, while the MnCr2O4-CoCr2O4 solid solution shows slightly positive deviations. The difference in the Gibbs free energy of formation of the two pure chromites and aluminates derived from the results of this study are in good agreement with recent results obtained from solid oxide galvanic cells and gas-equilibrium techniques.

A solid-state probe for SO2/SO3 based on Na2SO4 I electrolyte

Conductivity measurements as a function of temperature and partial pressures of SOs, SO2, and O2, and transference experiments indicate that the transport number of Na + ions is unity in Na2SO4-I. A concentration cell based on this electrolyte Pt, O2' + SO2' + SOs'/Na2SO4-I/SOa" + SO~" + O~", Pt produces emf's that are in agreement with those calculated from the Nernst equation when equilibrium is assumed between the gas species at the electrodes. The cell can be used for monitoring the SO#SOs pollution in air, and in combination with an oxygen probe can be used for the determination of SO=/SOs concentrations in coal combustion reactors, for the evaluation of the partial pressure of $2 in coal gasification systems, and for emission control in nonferrous smelters using sulfide ores. The probe is similar to that developed recently by Gauthier et aL (4, 5) using K=SO4 as the electrolyte, but can operate at higher pressures of SO3. Because of the greater polarizing power of the Na+ ion compared to the K + ion, Na2S207 is less stable and can be formed only at a considerably higher pressure of S03 than that required for K~20~.

Thermodynamic properties of Fe3O4 FeAl2O4 spinel solid solution

The Gibbs energy of mixing for the system Fe3O4-FeAl2O4 was determined at 1573 K using a gas-metal-oxide equilibration technique. Oxide solid solution samples were equilibrated with Pt foils under controlled CO+CO2 gas streams. The equilibrium iron concentration in the foil was determined by chemical analysis. The cation distribution between tetrahedral and octahedral sites in the spinel crystal can be calculated from site-preference energies and used as an alternate method of determining some thermodynamic properties, including the Gibbs energy of mixing. The solvus occurring at low temperatures in the system Fe3C4-FeAl2C4 was used to derive the effect of lattice distortion due to cation size difference on the enthalpy of mixing and to obtain a better approximation to the measured thermodynamic quantities.

Thermodynamic properties of Fe3O4 FeV2O4 and Fe3O4 FeCr2O4 spinel solid solutions

Solid solutions of Fe304-FeV204 and Fe304-FeCr204 were prepared and equilibrated with Pt under controlled streams of CO/CO, gas mixtures at 1673 K. The concentration of Fe in Pt was used to determine the activity of Fe304 in the solid solutions. The activity of the second component was calculated by Gibbshhem integration. From these data, the Gibbs energy of mixing was derived for both systems. The experimental results and theoretical values which are determined from calculated cation distribution compare favorably in the case of vanadite solid solutions but not in the case of chromite solid solutions. The difference is attributed to a heat term arising from lattice distortion due to cation size difference. The positive heat of mixing will give rise to a miscibility gap in the system Fe304-FeCr204 at lower temperatures.

Thermodynamic properties of Cr Mo solid alloys

The activity of Cr in solid Cr-Mo alloys has been measured at 1873 K using a metal-oxide-gas equilibrium technique. Thin foils of Mo were equilibrated with solid Cr203 under flowing gas mixtures of argon, hydrogen and watervapourof known composition. The equilibrium concentration of Cr in Mo was determined by chemical analysis. These measurements indicate positive deviations from Raoult's law. The activity data obtained in the study at 1873 K are combined with free energy of mixing at 1471 K, calorimetric enthalpy of mixing at 1673 K, and experimental evidence of phase separation at lower temperatures, reported in the literature, to obtain an optimised set of thermodynamic parameters for the Cr-Mosystem in the solid state.

Applications of solid electrolytes in galvanic sensors

The concepts and theoretical origins of conduction domains for solid electrolytes and electrode polarization are outlined briefly. The point electrode made of the ' solid electrolyte material is useful for deflecting the semipermeability flux away from the electrode. The emf of galvanic sensors consisting of two solid electrolytes in intimate contact with each other and in which transport occurs by a common ion is reviewed. The voltage of such cells depends on the chemical potential of the active species at the interface between the two electrolytes, which can be evaluated from the transport properties of electrolytes using a numerical procedure. The factors governing the speed of response of solid electrolyte gas sensors are analyzed. A rigorous expression for the emf of non-isothermal galvanic sensors and the criterion for the design of temperature compensated reference electrodes for nonisothermal galvanic sensors are outlined. Non-isothermal sensors are useful for the continuous monitoring of concentrations or chemical potentials in reactive systems at high temperatures. The principles of operation of galvanic sensors for oxygen, sulphur, oxides of sulphur (SOx,x=2,3), carbon, oxides uf carbon (COx,x= 1,2), oxides of nitrogen (NOx,x= 1,2) and silicon are discussed. The use of auxiliary electrodes in galvanic sensors to expand the detection capability of known solid electrolytes to a large number of species is explained with reference to sensors for sulphur and oxides of sulphur (SOx,x=2,3).Finally the cause of the common errors in galvanic measurements and test for the correct functioning of galvanic sensors is given.