Computational modeling of intermediate temperature proton exchange membrane (PEM) fuel cells

A two-phase three-dimensional computational model of an intermediate temperature (120--190°C) proton exchange membrane (PEM) fuel cell is presented. This represents the first attempt to model PEM fuel cells employing intermediate temperature membranes, in this case, phosphoric acid doped polybenzimidazole (PBI). To date, mathematical modeling of PEM fuel cells has been restricted to low temperature operation, especially to those employing Nafion ® membranes; while research on PBI as an intermediate temperature membrane has been solely at the experimental level. This work is an advancement in the state of the art of both these fields of research. With a growing trend toward higher temperature operation of PEM fuel cells, mathematical modeling of such systems is necessary to help hasten the development of the technology and highlight areas where research should be focused.^ This mathematical model accounted for all the major transport and polarization processes occurring inside the fuel cell, including the two phase phenomenon of gas dissolution in the polymer electrolyte. Results were presented for polarization performance, flux distributions, concentration variations in both the gaseous and aqueous phases, and temperature variations for various heat management strategies. The model predictions matched well with published experimental data, and were self-consistent.^ The major finding of this research was that, due to the transport limitations imposed by the use of phosphoric acid as a doping agent, namely low solubility and diffusivity of dissolved gases and anion adsorption onto catalyst sites, the catalyst utilization is very low (∼1--2%). Significant cost savings were predicted with the use of advanced catalyst deposition techniques that would greatly reduce the eventual thickness of the catalyst layer, and subsequently improve catalyst utilization. The model also predicted that an increase in power output in the order of 50% is expected if alternative doping agents to phosphoric acid can be found, which afford better transport properties of dissolved gases, reduced anion adsorption onto catalyst sites, and which maintain stability and conductive properties at elevated temperatures.^

Weak coupling of acoustic-like phonons and magnon dynamics in incommensurate spin ladder compound Sr14Cu24O41

Intriguing lattice dynamics has been predicted for aperiodic crystals that contain incommensurate substructures. Here we report inelastic neutron scattering measurements of phonon and magnon dispersions in Sr14Cu24O41, which contains incommensurate one-dimensional (1D) chain and two-dimensional (2D) ladder substructures. Two distinct acoustic phonon-like modes, corresponding to the sliding motion of one sublattice against the other, are observed for atomic motions polarized along the incommensurate axis. In the long wavelength limit, it is found that the sliding mode shows a remarkably small energy gap of 1.7-1.9 meV, indicating very weak interactions between the two incommensurate sublattices. The measurements also reveal a gapped and steep linear magnon dispersion of the ladder sublattice. The high group velocity of this magnon branch and weak coupling with acoustic phonons can explain the large magnon thermal conductivity in Sr14Cu24O41 crystals. In addition, the magnon specific heat is determined from the measured total specific heat and phonon density of states, and exhibits a Schottky anomaly due to gapped magnon modes of the spin chains. These findings offer new insights into the phonon and magnon dynamics and thermal transport properties of incommensurate magnetic crystals that contain low-dimensional substructures.

Application of Numerical Methods to Study Arrangement and Fracture of Lithium-Ion Microstructure

The focus of this work is to develop and employ numerical methods that provide characterization of granular microstructures, dynamic fragmentation of brittle materials, and dynamic fracture of three-dimensional bodies.

We first propose the fabric tensor formalism to describe the structure and evolution of lithium-ion electrode microstructure during the calendaring process. Fabric tensors are directional measures of particulate assemblies based on inter-particle connectivity, relating to the structural and transport properties of the electrode. Applying this technique to X-ray computed tomography of cathode microstructure, we show that fabric tensors capture the evolution of the inter-particle contact distribution and are therefore good measures for the internal state of and electronic transport within the electrode.

We then shift focus to the development and analysis of fracture models within finite element simulations. A difficult problem to characterize in the realm of fracture modeling is that of fragmentation, wherein brittle materials subjected to a uniform tensile loading break apart into a large number of smaller pieces. We explore the effect of numerical precision in the results of dynamic fragmentation simulations using the cohesive element approach on a one-dimensional domain. By introducing random and non-random field variations, we discern that round-off error plays a significant role in establishing a mesh-convergent solution for uniform fragmentation problems. Further, by using differing magnitudes of randomized material properties and mesh discretizations, we find that employing randomness can improve convergence behavior and provide a computational savings.

The Thick Level-Set model is implemented to describe brittle media undergoing dynamic fragmentation as an alternative to the cohesive element approach. This non-local damage model features a level-set function that defines the extent and severity of degradation and uses a length scale to limit the damage gradient. In terms of energy dissipated by fracture and mean fragment size, we find that the proposed model reproduces the rate-dependent observations of analytical approaches, cohesive element simulations, and experimental studies.

Lastly, the Thick Level-Set model is implemented in three dimensions to describe the dynamic failure of brittle media, such as the active material particles in the battery cathode during manufacturing. The proposed model matches expected behavior from physical experiments, analytical approaches, and numerical models, and mesh convergence is established. We find that the use of an asymmetrical damage model to represent tensile damage is important to producing the expected results for brittle fracture problems.

The impact of this work is that designers of lithium-ion battery components can employ the numerical methods presented herein to analyze the evolving electrode microstructure during manufacturing, operational, and extraordinary loadings. This allows for enhanced designs and manufacturing methods that advance the state of battery technology. Further, these numerical tools have applicability in a broad range of fields, from geotechnical analysis to ice-sheet modeling to armor design to hydraulic fracturing.

Decolourisation of Acid orange 7 in a microbial fuel cell with a laccase-based biocathode: Influence of mitigating pH changes in the cathode chamber

Biocathodes may be a suitable replacement of platinum in microbial fuel cells (MFCs) if the cost of MFCs is to be reduced. However, the use of enzymes as bio-cathodes is fraught with loss of activity as time progresses. A possible cause of this loss in activity might be pH increase in the cathode as pH gradients in MFCs are well known. This pH increase is however, accompanied by simultaneous increase in salinity; therefore salinity may be a confounding variable. This study investigated various ways of mitigating pH changes in the cathode of MFCs and their effect on laccase activity and decolourisation of a model azo dye Acid orange 7 in the anode chamber. Experiments were run with catholyte pH automatically controlled via feedback control or by using acetate buffers (pH 4.5) of various strength (100 mM and 200 mM), with CMI7000 as the cation exchange membrane. A comparison was also made between use of CMI7000 and Nafion 117 as the transport properties of cations for both membranes (hence their potential effects on pH changes in the cathode) are different.

Electrochemical Investigations into Blended Electrolytes Containing Ionic Liquids and a Glyme Solvent for Li-O­2 batteries

The Li-O2 battery may theoretically possess practical gravimetric energy densities several times greater than the current state-of-the-art Li-ion batteries.1 This magnitude of development is a requisite for true realization of electric vehicles capable of competing with the traditional combustion engine. However, significant challenges must be addressed before practical application may be considered. These include low efficiencies, low rate capabilities and the parasitic decomposition reactions of electrolyte/electrode materials resulting in very poor rechargeability.2-4 Ionic liquids, ILs, typically display several properties, extremely low vapor pressure and high electrochemical and thermal stability, which make them particularly interesting for Li-O2 battery electrolytes. However, the typically sluggish transport properties generally inhibit rate performance and cells suffer similar inefficiencies during cycling.5,6

In addition to the design of new ILs with tailored properties, formulating blended electrolytes using molecular solvents with ILs has been considered to improve their performance.7,8 In this work, we will discuss the physical properties vs. the electrochemical performance of a range of formulated electrolytes based on tetraglyme, a benchmark Li-O2 battery electrolyte solvent, and several ILs. The selected ILs are based on the bis{(trifluoromethyl)sulfonyl}imide anion and alkyl/ether functionalized cyclic alkylammonium cations, which exhibit very good stability and moderate viscosity.9 O2 electrochemistry will be investigated in these media using macro and microdisk voltammetry and O2 solubility/diffusivity is quantified as a function of the electrolyte formulation. Furthermore, galvanostatic cycling of selected electrolytes in Li-O2 cells will be discussed to probe their practical electrochemical performance. Finally, the physical characterization of the blended electrolytes will be reported in parallel to further determine structure (or formulation) vs. property relationships and to, therefore, assess the importance of certain electrolyte properties (viscosity, O2supply capability, donor number) on their performance.

This work was funded by the EPSRC (EP/L505262/1) and Innovate UK for the Practical Lithium-Air Batteries project (project number: 101577).

1. P. G. Bruce, S. A. Freunberger, L. J. Hardwick and J.-M. Tarascon, Nat. Mater., 11, 19 (2012).

2. S. A. Freunberger, Y. Chen, N. E. Drewett, L. J. Hardwick, F. Barde and P. G. Bruce, Angew. Chem., Int. Ed., 50, 8609 (2011).

3. B. D. McCloskey, A. Speidel, R. Scheffler, D. C. Miller, V. Viswanathan, J. S. Hummelshøj, J. K. Nørskov and A. C. Luntz, J. Phys. Chem. Lett., 3, 997 (2012).

4. D. G. Kwabi, T. P. Batcho, C. V. Amanchukwu, N. Ortiz-Vitoriano, P. Hammond, C. V. Thompson and Y. Shao-Horn, J. Phys. Chem. Lett., 5, 2850 (2014).

5. Z. H. Cui, W. G. Fan and X. X. Guo, J. Power Sources, 235, 251 (2013).

6. F. Soavi, S. Monaco and M. Mastragostino, J. Power Sources, 224, 115 (2013).

7. L. Cecchetto, M. Salomon, B. Scrosati and F. Croce, J. Power Sources, 213, 233 (2012).

8. A. Khan and C. Zhao, Electrochem. Commun., 49, 1 (2014).

9. Z. J. Chen, T. Xue and J.-M. Lee, RSC Adv., 2, 10564 (2012).

An ether-functionalised cyclic sulfonium based ionic liquid as an electrolyte for electrochemical double layer capacitors

A novel cyclic sulfonium cation-based ionic liquid (IL) with an ether-group appendage and the bis{(trifluoromethyl)sulfonyl}imide anion was synthesised and developed for electrochemical double layer capacitor (EDLC) testing. The synthesis and chemical-physical characterisation of the ether-group containing IL is reported in parallel with a similarly sized alkyl-functionalised sulfonium IL. Results of the chemical-physical measurements demonstrate how important transport properties, i.e. viscosity and conductivity, can be promoted through the introduction of the ether-functionality without impeding thermal, chemical or electrochemical stability of the IL. Although the apparent transport properties are improved relative to the alkyl-functionalised analogue, the ether-functionalised sulfonium cation-based IL exhibits moderately high viscosity, and poorer conductivity, when compared to traditional EDLC electrolytes based on organic solvents (propylene carbonate and acetonitrile). Electrochemical testing of the ether-functionalised sulfonium IL was conducted using activated carbon composite electrodes to inspect the performance of the IL as a solvent-free electrolyte for EDLC application. Good cycling stability was achieved over the studied range and the performance was comparable to other solvent free,
IL-based EDLC systems. Nevertheless, limitations of the attainable performance are primarily the result of sluggish transport properties and a restricted operative voltage of the IL, thus highlighting key aspects of this field which require further attention.

On the statistics of transverse permeability of randomly distributed fibers

An RVE–based stochastic numerical model is used to calculate the permeability of randomly generated porous media at different values of the fiber volume fraction for the case of transverse flow in a unidirectional ply. Analysis of the numerical results shows that the permeability is not normally distributed. With the aim of proposing a new understanding on this particular topic, permeability data are fitted using both a mixture model and a unimodal distribution. Our findings suggest that permeability can be fitted well using a mixture model based on the lognormal and power law distributions. In case of a unimodal distribution, it is found, using the maximum-likelihood estimation method (MLE), that the generalized extreme value (GEV) distribution represents the best fit. Finally, an expression of the permeability as a function of the fiber volume fraction based on the GEV distribution is discussed in light of the previous results.

The use of binary mixtures of 1-butyl-1-methylpyrrolidinium bis{(trifluoromethyl)sulfonyl}imide and aliphatic nitrile solvents as electrolyte for supercapacitors

Abstract The development of high voltage electrolytes is one of the key aspects for increasing both energy and power density of electrochemical double layer capacitors (EDLCs). The usage of blends of ionic liquids and organic solvents has been considered as a feasible strategy since these electrolytes combine high usable voltages and good transport properties at the same time. In this work, the ionic liquid 1-butyl-1-methylpyrrolidinium bis{(trifluoromethyl)sulfonyl}imide ([Pyrr14][TFSI]) was mixed with two nitrile-based organic solvents, namely butyronitrile and adiponitrile, and the resulting blends were investigated regarding their usage in electrochemical double layer capacitors. Both blends have a high electrochemical stability, which was confirmed by prolonged float tests at 3.2 V, as well as, good transport properties. In fact, the butyronitrile blend reaches a conductivity of 17.14 mS·cm−1 and a viscosity of 2.46 mPa·s at 20 °C, which is better than the state-of-the-art electrolyte (1 mol·dm−3 of tetraethylammonium tetrafluoroborate in propylene carbonate).

Epitaxially Grown Strained Pentacene Thin Film on Graphene Membrane

Organic-graphene system has emerged as a new platform for various applications such as flexible organic photovoltaics and organic light emitting diodes. Due to its important implication in charge transport, the study and reliable control of molecular packing structures at the graphene-molecule interface are of great importance for successful incorporation of graphene in related organic devices. Here, an ideal membrane of suspended graphene as a molecular assembly template is utilized to investigate thin-film epitaxial behaviors. Using transmission electron microscopy, two distinct molecular packing structures of pentacene on graphene are found. One observed packing structure is similar to the well-known bulk-phase, which adapts a face-on molecular orientation on graphene substrate. On the other hand, a rare polymorph of pentacene crystal, which shows significant strain along the c-axis, is identified. In particular, the strained film exhibits a specific molecular orientation and a strong azimuthal correlation with underlying graphene. Through ab initio electronic structure calculations, including van der Waals interactions, the unusual polymorph is attributed to the strong graphene-pentacene interaction. The observed strained organic film growth on graphene demonstrates the possibility to tune molecular packing via graphene-molecule interactions.

Graphene-MoS2 Hybrid Technology for Large-Scale Two-Dimensional Electronics

Two-dimensional (2D) materials have generated great interest in the last few years as a new toolbox for electronics. This family of materials includes, among others, metallic graphene, semiconducting transition metal dichalcogenides (such as MoS2) and insulating Boron Nitride. These materials and their heterostructures offer excellent mechanical flexibility, optical transparency and favorable transport properties for realizing electronic, sensing and optical systems on arbitrary surfaces. In this work, we develop several etch stop layer technologies that allow the fabrication of complex 2D devices and present for the first time the large scale integration of graphene with molybdenum disulfide (MoS2) , both grown using the fully scalable CVD technique. Transistor devices and logic circuits with MoS2 channel and graphene as contacts and interconnects are constructed and show high performances. In addition, the graphene/MoS2 heterojunction contact has been systematically compared with MoS2-metal junctions experimentally and studied using density functional theory. The tunability of the graphene work function significantly improves the ohmic contact to MoS2. These high-performance large-scale devices and circuits based on 2D heterostructure pave the way for practical flexible transparent electronics in the future. The authors acknowledge financial support from the Office of Naval Research (ONR) Young Investigator Program, the ONR GATE MURI program, and the Army Research Laboratory. This research has made use of the MI.

Stability of a Nanofiltration Membrane After Contact with a Low-Level Liquid Radioactive Waste

This study investigated the treatment of a liquid radioactive waste containing uranium (235U + 238U) using nanofiltration membranes. The membranes were immersed in the waste for 24–5000 h, and their transport properties were evaluated before and after the immersion. Surface of the membranes changed after immersion in the waste. The SW5000 h specimen lost its coating layer of polyvinyl alcohol, and its rejection of sulfate ions and uranium decreased by about 35% and 30%, respectively. After immersion in the waste, the polyamide selective layer of the membranes became less thermally stable than that before immersion.

BETA-GA2O3: A TRANSPARENT CONDUCTIVE OXIDE FOR POTENTIAL RESISTIVE SWITCHING APPLICATIONS

Thesis (Ph.D.)--University of Washington, 2016-08

Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

Effect of experimental parameters on alginate/chitosan microparticles for BCG encapsulation

The aim of the present study was to develop novel Mycobacterium bovis bacille Calmette-Guérin (BCG)-loaded polymeric microparticles with optimized particle surface characteristics and biocompatibility, so that whole live attenuated bacteria could be further used for pre-exposure vaccination against Mycobacterium tuberculosis by the intranasal route. BCG was encapsulated in chitosan and alginate microparticles through three different polyionic complexation methods by high speed stirring. For comparison purposes, similar formulations were prepared with high shear homogenization and sonication. Additional optimization studies were conducted with polymers of different quality specifications in a wide range of pH values, and with three different cryoprotectors. Particle morphology, size distribution, encapsulation efficiency, surface charge, physicochemical properties and biocompatibility were assessed. Particles exhibited a micrometer size and a spherical morphology. Chitosan addition to BCG shifted the bacilli surface charge from negative zeta potential values to strongly positive ones. Chitosan of low molecular weight produced particle suspensions of lower size distribution and higher stability, allowing efficient BCG encapsulation and biocompatibility. Particle formulation consistency was improved when the availability of functional groups from alginate and chitosan was close to stoichiometric proportion. Thus, the herein described microparticulate system constitutes a promising strategy to deliver BCG vaccine by the intranasal route.

Fluctuations in and out of equilibrium: Thermalization, quantum measurements and Coulomb disorder

This purely theoretical thesis covers aspects of two contemporary research fields: the non-equilibrium dynamics in quantum systems and the electronic properties of three-dimensional topological insulators. In the first part we investigate the non-equilibrium dynamics in closed quantum systems. Thanks to recent technologies, especially from the field of ultracold quantum gases, it is possible to realize such systems in the laboratory. The focus is on the influence of hydrodynamic slow modes on the thermalization process. Generic systems in equilibrium, either classical or quantum, in equilibrium are described by thermodynamics. This is characterized by an ensemble of maximal entropy, but constrained by macroscopically conserved quantities. We will show that these conservation laws slow down thermalization and the final equilibrium state can be approached only algebraically in time. When the conservation laws are violated thermalization takes place exponential in time. In a different study we calculate probability distributions of projective quantum measurements. Newly developed quantum microscopes provide the opportunity to realize new measurement protocols which go far beyond the conventional measurements of correlation functions. The second part of this thesis is dedicated to a new class of materials known as three-dimensional topological insulators. Also here new experimental techniques have made it possible to fabricate these materials to a high enough quality that their topological nature is revealed. However, their transport properties are not fully understood yet. Motivated by unusual experimental results in the optical conductivity we have investigated the formation and thermal destruction of spatially localized electron- and hole-doped regions. These are caused by charged impurities which are introduced into the material in order to make the bulk insulating. Our theoretical results are in agreement with the experiment and can explain the results semi-quantitatively. Furthermore, we study emergent lengthscales in the bulk as well as close to the conducting surface.