Modeling, shape analysis and computation of the equilibrium pore shape near a PEM–PEM intersection

Abstract In this paper we study the equilibrium shape of an interface that represents the lateral boundary of a pore channel embedded in an elastomer. The model consists of a system of PDEs, comprising a linear elasticity equation for displacements within the elastomer and a nonlinear Poisson equation for the electric potential within the channel (filled with protons and water). To determine the equilibrium interface, a variational approach is employed. We analyze: (i) the existence and uniqueness of the electrical potential, (ii) the shape derivatives of state variables and (iii) the shape differentiability of the corresponding energy and the corresponding Euler–Lagrange equation. The latter leads to a modified Young–Laplace equation on the interface. This modified equation is compared with the classical Young–Laplace equation by computing several equilibrium shapes, using a fixed point algorithm.

[1]  D. Kinderlehrer,et al.  An introduction to variational inequalities and their applications , 1980 .

[2]  Multiple existence of solutions for a semilinear elliptic problem with Neumann boundary condition , 2006 .

[3]  J. Simon.  Differentiation with Respect to the Domain in Boundary Value Problems , 1980 .

[4]  Robert B. Moore,et al.  State of understanding of nafion. , 2004, Chemical reviews.

[5]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

[6]  Jean-François Remacle,et al.  A stabilized finite element method using a discontinuous level set approach for the computation of bubble dynamics , 2007, J. Comput. Phys..

[7]  Haim Brezis,et al.  Nonlinear elliptic equations with measures revisited , 2007, 1312.6495.

[8]  H. Brezis.  Functional Analysis, Sobolev Spaces and Partial Differential Equations , 2010 .

[9]  Sven-Joachim Kimmerle,et al.  Effects of radially dependent parameters on proton transport in polymer electrolyte membrane nanopores. , 2011, Journal of Chemical Physics.

[10]  A. Quarteroni,et al.  Fluid–structure algorithms based on Steklov–Poincaré operators , 2006 .

[11]  L. Orsina,et al.  Semilinear equations with exponential nonlinearity and measure data , 2005 .

[12]  P. Berg,et al.  Exact solution of an electro-osmotic flow problem in a cylindrical channel of polymer electrolyte membranes , 2009, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[13]  Peter Berg,et al.  An Electrohydrodynamic Equilibrium Shape Problem for Polymer Electrolyte Membranes in Fuel Cells , 2011, System Modelling and Optimization.

[14]  Haim Brezis,et al.  Semi-linear second-order elliptic equations in L 1 , 1973 .

[15]  Philippe G. Ciarlet,et al.  Lectures on Three-Dimensional Elasticity , 1983 .

[16]  K. S. Dhathathreyan,et al.  Polymer Electrolyte Membrane Fuel Cell , 2007 .

[17]  Hans-Jürgen Butt,et al.  Physics and Chemistry of Interfaces , 2003 .

[18]  S. Paddison,et al.  Transport in proton conductors for fuel-cell applications: simulations, elementary reactions, and phenomenology. , 2004, Chemical reviews.

[19]  Qiang Chen,et al.  Parallel cylindrical water nanochannels in Nafion fuel-cell membranes. , 2008, Nature materials.

[20]  S. Agmon,et al.  Estimates near the boundary for solutions of elliptic partial differential equations satisfying general boundary conditions. I , 1959 .

[21]  Annegret Glitzky,et al.  On energy estimates for electro-diffusion equations arising in semiconductor technology , 2005 .