Variational Methods and Planar Elliptic Growth

A nested family of growing or shrinking planar domains is called a Laplacian growth process if the normal velocity of each domain's boundary is proportional to the gradient of the domain's Green function with a fixed singularity on the interior. In this paper we review the Laplacian growth model and its key underlying assumptions, so that we may consider a generalization to so-called elliptic growth, wherein the Green function is replaced with that of a more general elliptic operator--this models, for example, inhomogeneities in the underlying plane. In this paper we continue the development of the underlying mathematics for elliptic growth, considering perturbations of the Green function due to those of the driving operator, deriving characterizations and examples of growth, developing a weak formulation of growth via balayage, and discussing of a couple of inverse problems in the spirit of Calder\'on. We conclude with a derivation of a more delicate, reregularized model for Hele-Shaw flow

Theory of voltammetry in charged porous media

We couple the Leaky Membrane Model, which describes the diffusion and electromigration of ions in a homogenized porous medium of fixed background charge, with Butler-Volmer reaction kinetics for flat electrodes separated by such a medium in a simple mathematical theory of voltammetry. The model is illustrated for the prototypical case of copper electro-deposition/dissolution in aqueous charged porous media. We first consider the steady state with three different experimentally relevant boundary conditions and derive analytical or semi-analytical expressions for concentration profiles, electric potential profiles, current-voltage relations and overlimiting conductances. Next, we perform nonlinear least squares fitting on experimental data, consider the transient response for linear sweep voltammetry and demonstrate good agreement of the model predictions with experimental data. The experimental datasets are for copper electrodeposition from copper(II) sulfate solutions in a variety of nanoporous media, such as anodic aluminum oxide, cellulose nitrate and polyethylene battery separators, whose internal surfaces are functionalized with positively and negatively charged polyelectrolyte polymers.Comment: 39 pages, 12 figures, 5 tables; clarified where other parameters are taken from and fixed typo

Electrochemical Impedance Imaging via the Distribution of Diffusion Times

We develop a mathematical framework to analyze electrochemical impedance spectra in terms of a distribution of diffusion times (DDT) for a parallel array of random finite-length Warburg (diffusion) or Gerischer (reaction-diffusion) circuit elements. A robust DDT inversion method is presented based on Complex Nonlinear Least Squares (CNLS) regression with Tikhonov regularization and illustrated for three cases of nanostructured electrodes for energy conversion: (i) a carbon nanotube supercapacitor, (ii) a silicon nanowire Li-ion battery, and (iii) a porous-carbon vanadium flow battery. The results demonstrate the feasibility of non-destructive "impedance imaging" to infer microstructural statistics of random, heterogeneous materials

Homogenization of the Poisson-Nernst-Planck Equations for Ion Transport in Charged Porous Media

Effective Poisson-Nernst-Planck (PNP) equations are derived for macroscopic ion transport in charged porous media under periodic fluid flow by an asymptotic multi-scale expansion with drift. The microscopic setting is a two-component periodic composite consisting of a dilute electrolyte continuum (described by standard PNP equations) and a continuous dielectric matrix, which is impermeable to the ions and carries a given surface charge. Four new features arise in the upscaled equations: (i) the effective ionic diffusivities and mobilities become tensors, related to the microstructure; (ii) the effective permittivity is also a tensor, depending on the electrolyte/matrix permittivity ratio and the ratio of the Debye screening length to the macroscopic length of the porous medium; (iii) the microscopic fluidic convection is replaced by a diffusion-dispersion correction in the effective diffusion tensor; and (iv) the surface charge per volume appears as a continuous "background charge density", as in classical membrane models. The coefficient tensors in the upscaled PNP equations can be calculated from periodic reference cell problems. For an insulating solid matrix, all gradients are corrected by the same tensor, and the Einstein relation holds at the macroscopic scale, which is not generally the case for a polarizable matrix, unless the permittivity and electric field are suitably defined. In the limit of thin double layers, Poisson's equation is replaced by macroscopic electroneutrality (balancing ionic and surface charges). The general form of the macroscopic PNP equations may also hold for concentrated solution theories, based on the local-density and mean-field approximations. These results have broad applicability to ion transport in porous electrodes, separators, membranes, ion-exchange resins, soils, porous rocks, and biological tissues

Multiphase Porous Electrode Theory

Porous electrode theory, pioneered by John Newman and collaborators, provides a useful macroscopic description of battery cycling behavior, rooted in microscopic physical models rather than empirical circuit approximations. The theory relies on a separation of length scales to describe transport in the electrode coupled to intercalation within small active material particles. Typically, the active materials are described as solid solution particles with transport and surface reactions driven by concentration fields, and the thermodynamics are incorporated through fitting of the open circuit potential. This approach has fundamental limitations, however, and does not apply to phase-separating materials, for which the voltage is an emergent property of inhomogeneous concentration profiles, even in equilibrium. Here, we present a general theoretical framework for "multiphase porous electrode theory" implemented in an open-source software package called "MPET", based on electrochemical nonequilibrium thermodynamics. Cahn-Hilliard-type phase field models are used to describe the solid active materials with suitably generalized models of interfacial reaction kinetics. Classical concentrated solution theory is implemented for the electrolyte phase, and Newman's porous electrode theory is recovered in the limit of solid-solution active materials with Butler-Volmer kinetics. More general, quantum-mechanical models of Faradaic reactions are also included, such as Marcus-Hush-Chidsey kinetics for electron transfer at metal electrodes, extended for concentrated solutions. The full equations and numerical algorithms are described, and a variety of example calculations are presented to illustrate the novel features of the software compared to existing battery models

Theory of Nucleation in Phase-separating Nanoparticles

The basic physics of nucleation in solid \hl{single-crystal} nanoparticles is revealed by a phase-field theory that includes surface energy, chemical reactions and coherency strain. In contrast to binary fluids, which form arbitrary contact angles at surfaces, complete "wetting" by one phase is favored at binary solid surfaces. Nucleation occurs when surface wetting becomes unstable, as the chemical energy gain (scaling with area) overcomes the elastic energy penalty (scaling with volume). The nucleation barrier thus decreases with the area-to-volume ratio and vanishes below a critical size, and nanoparticles tend to transform in order of increasing size, leaving the smallest particles homogeneous (in the phase of lowest surface energy). The model is used to simulate phase separation in realistic nanoparticle geometries for \ce{Li_XFePO4}, a popular cathode material for Li-ion batteries, and collapses disparate experimental data for the nucleation barrier, with no adjustable parameters. Beyond energy storage, the theory generally shows how to tailor the elastic and surface properties of a solid nanostructure to achieve desired phase behavior.Comment: 7 pages, 4 fig

Theory of Sorption Hysteresis in Nanoporous Solids: II. Molecular condensation

Motivated by the puzzle of sorption hysteresis in Portland cement concrete or cement paste, we develop in Part II of this study a general theory of vapor sorption and desorption from nanoporous solids, which attributes hysteresis to hindered molecular condensation with attractive lateral interactions. The classical mean-field theory of van der Waals is applied to predict the dependence of hysteresis on temperature and pore size, using the regular solution model and gradient energy of Cahn and Hilliard. A simple "hierarchical wetting" model for thin nanopores is developed to describe the case of strong wetting by the first monolayer, followed by condensation of nanodroplets and nanobubbles in the bulk. The model predicts a larger hysteresis critical temperature and enhanced hysteresis for molecular condensation across nanopores at high vapor pressure than within monolayers at low vapor pressure. For heterogeneous pores, the theory predicts sorption/desorption sequences similar to those seen in molecular dynamics simulations, where the interfacial energy (or gradient penalty) at nanopore junctions acts as a free energy barrier for snap-through instabilities. The model helps to quantitatively understand recent experimental data for concrete or cement paste wetting and drying cycles and suggests new experiments at different temperatures and humidity sweep rates.Comment: 26 pages, 10 fig