This report documents the author's efforts in the deterministic modeling of coppersulfidation corrosion on non-planar substrates such as diodes and electrical connectors. A new framework based on Goma was developed for multi-dimensional modeling of atmospheric copper-sulfidation corrosion on non-planar substrates. In this framework, the moving sulfidation front is explicitly tracked by treating the finite-element mesh as a pseudo solid with an arbitrary Lagrangian-Eulerian formulation and repeatedly performing re-meshing using CUBIT and re-mapping using MAPVAR. Three one-dimensional studies were performed for verifying the framework in asymptotic regimes. Limited model validation was also carried out by comparing computed copper-sulfide thickness with experimental data. The framework was first demonstrated in modeling one-dimensional copper sulfidation with charge separation. It was found that both the thickness of the space-charge layers and the electrical potential at the sulfidation surface decrease rapidly as the Cu 2 S layer thickens initially but eventually reach equilibrium values as Cu 2 S layer becomes sufficiently thick; it was also found that electroneutrality is a reasonable approximation and that the electro-migration flux may be estimated by using the equilibrium potential difference between the sulfidation and annihilation surfaces when the Cu 2 S layer is sufficiently thick. The framework was then employed to model copper sulfidation in the solid-state-diffusion controlled regime (i.e. stage II sulfidation) on a prototypical diode until a continuous Cu 2 S film was formed on the diode surface. The framework was also applied to model copper sulfidation on an intermittent electrical contact between a gold-plated copper pin and gold-plated copper pad; the presence of Cu 2 S was found to raise the effective electrical resistance drastically. Lastly, future research needs in modeling atmospheric copper sulfidation are discussed.

In this report, we document the accomplishments in our Laboratory Directed Research and Development project in which we employed a technical approach of combining experiments with computational modeling and analyses to elucidate the performance of hydrogen-fed proton exchange membrane fuel cells (PEMFCs). In the first part of this report, we document our focused efforts on understanding water transport in and removal from a hydrogen-fed PEMFC. Using a transparent cell, we directly visualized the evolution and growth of liquid-water droplets at the gas diffusion layer (GDL)/gas flow channel (GFC) interface. We further carried out a detailed experimental study to observe, via direct visualization, the formation, growth, and instability of water droplets at the GDL/GFC interface using a speciallydesigned apparatus, which simulates the cathode operation of a PEMFC. We developed a simplified model, based on our experimental observation and data, for predicting the onset of water-droplet instability at the GDL/GFC interface. Using a state-of-the-art neutron imaging instrument available at NIST (National Institute of Standard and Technology), we probed liquid-water distribution inside an operating PEMFC under a variety of operating conditions and investigated effects of evaporation due to local heating by waste heat on water removal. Moreover, we developed computational models for analyzing the effects of micro-porous layer on net water transport across the membrane and GDL * Current address: Connecticut Global Fuel Cell Center, and Department of Mechanical Engineering, University of Connecticut, Storrs, CT 06269-5233.4 anisotropy on the temperature and water distributions in the cathode of a PEMFC. We further developed a two-phase model based on the multiphase mixture formulation for predicting the liquid saturation, pressure drop, and flow maldistribution across the PEMFC cathode channels.In the second part of this report, we document our efforts on modeling the electrochemical performance of PEMFCs. We developed a constitutive model for predicting proton conductivity in polymer electrolyte membranes and compared model prediction with experimental data obtained in our laboratory and from literature. Moreover, we developed a one-dimensional analytical model for predicting electrochemical performance of an idealized PEMFC with small surface over-potentials. Furthermore, we developed a multi-dimensional computer model, which is based on the finite-element method and a fully-coupled implicit solution scheme via Newton's technique, for simulating the performance of PEMFCs. We demonstrated utility of our finite-element model by comparing the computed current density distribution and overall polarization with those measured using a segmented cell. In the last part of this report, we document an exploratory experimental study on MEA (membrane electrode assembly) degradation. AcknowledgmentsThis work was funded by the Laboratory Directed Research and Development (LDRD) program at Sandia National Laboratories. We would...

No abstract

This document contains a summary of the work performed under the LDRD project entitled "Interface Physics in Microporous Media". The presence of fluid-fluid interfaces, which can carry non-zero stresses, distinguishes multiphase flows from more readily understood single-phase flows. In this work the physics active at these interfaces has been examined via a combined experimental and computational approach. One of the major difficulties of examining true microporous systems of the type found in filters, membranes, geologic media, etc. is the geometric uncertainty. To help facilitate the examination of transport at the pore-scale without this 4 complication, a significant effort has been made in the area of fabrication of both two-dimensional and three-dimensional micromodels. Using these micromodels, multiphase flow experiments have been performed for liquid-liquid and liquid-gas systems. Laser scanning confocal microscopy has been utilized to provide high resolution, three-dimensional reconstructions as well as time resolved, twodimensional reconstructions. Computational work has focused on extending lattice Boltzmann (LB) and finite element methods for probing the interface physics at the pore scale. A new LB technique has been developed that provides over 100x speed up for steady flows in complex geometries. A new LB model has been developed that allows for arbitrary density ratios, which has been a significant obstacle in applying LB to air-water flows. A new reduced order model has been developed and implemented in finite element code for examining non-equilibrium wetting in microchannel systems.These advances will enhance Sandia's ability to quantitatively probe the rich interfacial physics present in microporous systems. 5 ACKNOWLEDGMENTSWe would like to acknowledge the contributions of Joanne Fredrich who proposed this project and stewarded it during her time at Sandia. We would also like to acknowledge the work of both Mike Hickner and David Holdych who contributed to this work while at Sandia. We would also like to thank Raymond O. Cote for fabricating the cartridge for the micromodels and Brian Wroblewski and David J. Saiz for performing the photolithography work. 6 7

Two of the most daunting problems facing humankind in the twenty-first century are energy security and climate change. This report summarizes work accomplished towards addressing these problems through the execution of a Grand Challenge LDRD project (FY09-11). The vision of Sunshine to Petrol is captured in one deceptively simple chemical equation: Solar Energy + xCO 2 + (x+1) H 2 O → C x H 2x+2 (liquid fuel) + (1.5x+.5) O 2