Browsing by Author "Müller, Rüdiger"

Now showing 1 - 18 of 18
  • Thumbnail Image
    Item
    Bulk-surface electro-thermodynamics and applications to electrochemistry
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Dreyer, Wolfgang; Guhlke, Clemens; Müller, Rüdiger
    We propose a modeling framework for magnetizable, polarizable, elastic, viscous, heat conducting, reactive mixtures in contact with interfaces. To this end we first introduce bulk and surface balance equations that contain several constitutive quantities. For further modeling the constitutive quantities, we formulate constitutive principles. They are based on an axiomatic introduction of the entropy principle and the postulation of Galilean symmetry. We apply the proposed formalism to derive constitutive relations in a rather abstract setting. For illustration of the developed procedure, we state an explicit isothermal material model for liquid electrolyte|metal electrode interfaces in terms of free energy densities in the bulk and on the surface. Finally we give a survey of recent advancements in the understanding of electrochemical interfaces that were based on this model.
  • Thumbnail Image
    Item
    Bulk-Surface Electrothermodynamics and Applications to Electrochemistry
    (Basel : MDPI, 2018) Dreyer, Wolfgang; Guhlke, Clemens; Müller, Rüdiger
    We propose a modeling framework for magnetizable, polarizable, elastic, viscous, heat conducting, reactive mixtures in contact with interfaces. To this end, we first introduce bulk and surface balance equations that contain several constitutive quantities. For further modeling of the constitutive quantities, we formulate constitutive principles. They are based on an axiomatic introduction of the entropy principle and the postulation of Galilean symmetry. We apply the proposed formalism to derive constitutive relations in a rather abstract setting. For illustration of the developed procedure, we state an explicit isothermal material model for liquid electrolyte|metal electrode interfaces in terms of free energy densities in the bulk and on the surface. Finally, we give a survey of recent advancements in the understanding of electrochemical interfaces that were based on this model.
  • Thumbnail Image
    Item
    A continuum model for yttria-stabilised zirconia incorporating triple phase boundary, lattice structure and immobile oxide ions
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2019) Vágner, Petr; Guhlke, Clemens; Miloš, Vojtěch; Müller, Rüdiger; Fuhrmann, Jürgen
    A continuum model for yttria-stabilised zirconia (YSZ) in the framework of non-equilibrium thermodynamics is developed. Particular attention is given to i) modeling of the YSZ-metal-gas triple phase boundary, ii) incorporation of the lattice structure and immobile oxide ions within the free energy model and iii) surface reactions. A finite volume discretization method based on modified Scharfetter-Gummel fluxes is derived in order to perform numerical simulations. The model is used to study the impact of yttria and immobile oxide ions on the structure of the charged boundary layer and the double layer capacitance. Cyclic voltammograms of an air-half cell are simulated to study the effect of parameter variations on surface reactions, adsorption and anion diffusion.
  • Thumbnail Image
    Item
    A continuum model for yttria-stabilized zirconia incorporating triple phase boundary, lattice structure and immobile oxide ions
    (Berlin ; Heidelberg ; New York : Springer, 2019) Vágner, Petr; Guhlke, Clemens; Miloš, Vojtěch; Müller, Rüdiger; Fuhrmann, Jürgen
    A continuum model for yttria-stabilized zirconia (YSZ) in the framework of non-equilibrium thermodynamics is developed. Particular attention is given to (i) modeling of the YSZ-metal-gas triple phase boundary, (ii) incorporation of the lattice structure and immobile oxide ions within the free energy model and (iii) surface reactions. A finite volume discretization method based on modified Scharfetter-Gummel fluxes is derived in order to perform numerical simulations. The model is used to study the impact of yttria and immobile oxide ions on the structure of the charged boundary layer and the double layer capacitance. Cyclic voltammograms of an air-half cell are simulated to study the effect of parameter variations on surface reactions, adsorption and anion diffusion. © 2019, The Author(s).
  • Thumbnail Image
    Item
    Error control for the approximation of Allen-Cahn and Cahn-Hilliard equations with a logarithmic potential
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2010) Bartels, Sören; Müller, Rüdiger
    A fully computable upper bound for the finite element approximation error of Allen-Cahn and Cahn-Hilliard equations with logarithmic potentials is derived. Numerical experiments show that for the sharp interface limit this bound is robust past topological changes. Modifications of the abstract results to derive quasi-optimal error estimates in different norms for lowest order finite element methods are discussed and lead to weaker conditions on the residuals under which the conditional error estimates hold.
  • Thumbnail Image
    Item
    Galilean Bulk-Surface Electrothermodynamics and Applications to Electrochemistry
    (Basel : MDPI, 2023) Müller, Rüdiger; Landstorfer, Manuel
    In this work, the balance equations of non-equilibrium thermodynamics are coupled to Galilean limit systems of the Maxwell equations, i.e., either to (i) the quasi-electrostatic limit or (ii) the quasi-magnetostatic limit. We explicitly consider a volume (Formula presented.), which is divided into (Formula presented.) and (Formula presented.) by a possibly moving singular surface S, where a charged reacting mixture of a viscous medium can be present on each geometrical entity (Formula presented.). By the restriction to the Galilean limits of the Maxwell equations, we achieve that only subsystems of equations for matter and electromagnetic fields are coupled that share identical transformation properties with respect to observer transformations. Moreover, the application of an entropy principle becomes more straightforward and finally helps estimate the limitations of the more general approach based the full set of Maxwell equations. Constitutive relations are provided based on an entropy principle, and particular care is taken in the analysis of the stress tensor and the momentum balance in the general case of non-constant scalar susceptibility. Finally, we summarise the application of the derived model framework to an electrochemical system with surface reactions.
  • Thumbnail Image
    Item
    The impact of solvation and dissociation on the transport parameters of liquid electrolytes: Continuum modeling and numerical study
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Dreyer, Wolfgang; Guhlke, Clemens; Müller, Rüdiger
    Electro-thermodynamics provides a consistent framework to derive continuum models for electrochemical systems. For the application to a specific experimental system, the general model must be equipped with two additional ingredients: a free energy model to calculate the chemical potentials and a kinetic model for the kinetic coefficients. Suitable free energy models for liquid electrolytes incorporating ion-solvent interaction, finite ion sizes and solvation already exist and have been validated against experimental measurements. In this work, we focus on the modeling of the mobility coefficients based on MaxwellStefan setting and incorporate them into the general electro-thermodynamic framework. Moreover, we discuss the impact of model parameter on conductivity, transference numbers and salt diffusion coefficient. In particular, the focus is set on the solvation of ions and incomplete dissociation of a non-dilute electrolyte.
  • Thumbnail Image
    Item
    Modeling Polycrystalline Electrode-electrolyte Interfaces: The Differential Capacitance
    (Bristol : IOP Publishing, 2020) Müller, Rüdiger; Fuhrmann, Jürgen; Landstorfer, Manuel
    We present and analyze a model for polycrystalline electrode surfaces based on an improved continuum model that takes finite ion size and solvation into account. The numerical simulation of finite size facet patterns allows to study two limiting cases: While for facet size diameter dfacet →0 we get the typical capacitance of a spatially homogeneous but possible amorphous or liquid surface, in the limit 1[nm] < dfacet, an ensemble of non-interacting single crystal surfaces is approached. Already for moderate size of the facet diameters, the capacitance is remarkably well approximated by the classical approach of adding the single crystal capacities of the contributing facets weighted by their respective surface fraction. As a consequence, the potential of zero charge is not necessarily attained at a local minimum of capacitance, but might be located at a local capacitance maximum instead. Moreover, the results show that surface roughness can be accurately taken into account by multiplication of the ideally flat polycrystalline surface capacitance with a single factor. In particular, we find that the influence of the actual geometry of the facet pattern in negligible and our theory opens the way to a stochastic description of complex real polycrystal surfaces. © 2020 The Author(s). Published on behalf of The Electrochemical Society by IOP Publishing Limited.
  • Thumbnail Image
    Item
    Modeling polycrystalline electrode-electrolyte interfaces: The differential capacitance
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2019) Fuhrmann, Jürgen; Landstorfer, Manuel; Müller, Rüdiger
    We present and analyze a model for polycrystalline electrode surfaces based on an improved continuum model that takes finite ion size and solvation into account. The numerical simulation of finite size facet patterns allows to study two limiting cases: While for facet size diameter $d^facet to 0$ we get the typical capacitance of a spatially homogeneous but possible amorphous or liquid surface, in the limit $L^Debye << d^facet$ , an ensemble of non-interacting single crystal surfaces is approached. Already for moderate size of the facet diameters, the capacitance is remarkably well approximated by the classical approach of adding the single crystal capacities of the contributing facets weighted by their respective surface fraction. As a consequence, the potential of zero charge is not necessarily attained at a local minimum of capacitance, but might be located at a local capacitance maximum instead. Moreover, the results show that surface roughness can be accurately taken into account by multiplication of the ideally flat polycrystalline surface capacitance with a single factor. In particular, we find that the influence of the actual geometry of the facet pattern in negligible and our theory opens the way to a stochastic description of complex real polycrystal surfaces.
  • Thumbnail Image
    Item
    New insights on the interfacial tension of electrochemical interfaces and the Lippmann equation
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Dreyer, Wolfgang; Guhlke, Clemens; Landstorfer, Manuel; Neumann, Johannes; Müller, Rüdiger
    The Lippmann equation is considered as universal relationship between interfacial tension, double layer charge, and cell potential. Based on the framework of continuum thermo-electrodynamics we provide some crucial new insights to this relation. In a previous work we have derived a general thermodynamic consistent model for electrochemical interfaces, which showed a remarkable agreement to single crystal experimental data. Here we apply the model to a curved liquid metal electrode. If the electrode radius is large compared to the Debye length, we apply asymptotic analysis methods and obtain the Lippmann equation. We give precise definitions of the involved quantities and show that the interfacial tension of the Lippmann equation is composed of the surface tension of our general model, and contributions arising from the adjacent space charge layers. This finding is confirmed by a comparison of our model to experimental data of several mercury-electrolyte interfaces. We obtain qualitative and quantitative agreement in the 2V potential range for various salt concentrations. We also discuss the validity of our asymptotic model when the electrode radius is comparable to the Debye length.
  • Thumbnail Image
    Item
    A new perspective on the electron transfer: Recovering the Butler-Volmer equation in non-equilibrium thermodynamics
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Dreyer, Wolfgang; Guhlke, Clemens; Müller, Rüdiger
    Understanding and correct mathematical description of electron transfer reaction is a central question in electrochemistry. Typically the electron transfer reactions are described by the Butler-Volmer equation which has its origin in kinetic theories. The Butler-Volmer equation relates interfacial reaction rates to bulk quantities like the electrostatic potential and electrolyte concentrations. Since in the classical form, the validity of the Butler-Volmer equation is limited to some simple electrochemical systems, many attempts have been made to generalize the Butler-Volmer equation. Based on non-equilibrium thermodynamics we have recently derived a reduced model for the electrode-electrolyte interface. This reduced model includes surface reactions and adsorption but does not resolve the charge layer at the interface. Instead it is locally electroneutral and consistently incorporates all features of the double layer into a set of interface conditions. In the context of this reduced model we are able to derive a general Butler-Volmer equation. We discuss the application of the new Butler-Volmer equations to different scenarios like electron transfer reactions at metal electrodes, the intercalation process in lithium-iron-phosphate electrodes and adsorption processes. We illustrate the theory by an example of electroplating.
  • Thumbnail Image
    Item
    Optimal and robust a posteriori error estimates in L∞(L2) for the approximation of Allen-Cahn equations past singularities
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Bartels, Sören; Müller, Rüdiger
    Optimal a posteriori error estimates in L∞(L2) are derived for the finite element approximation of Allen-Cahn equations. The estimates depend on the inverse of a small parameter only in a low order polynomial and are valid past topological changes of the evolving interface. The error analysis employs an elliptic reconstruction of the approximate solution and applies to a large class of conforming, nonconforming, mixed, and discontinuous Galerkin methods. Numerical experiments illustrate the theoretical results.
  • Thumbnail Image
    Item
    Overcoming the shortcomings of the Nernst-Planck model
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2012) Dreyer, Wolfgang; Guhlke, Clemens; Müller, Rüdiger
    This is a study on electrolytes that takes a thermodynamically consistent coupling between mechanics and diffusion into account. It removes some inherent deficiencies of the popular Nernst-Planck model. A boundary problem for equilibrium processes is used to illustrate the new features of our model.
  • Thumbnail Image
    Item
    A posteriori error control for stationary coupled bulk-surface equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2015) Eigel, Martin; Müller, Rüdiger
    We consider a system of two coupled elliptic equations, one defined on a bulk domain and the other one on the boundary surface. Problems of this kind of problem are relevant for applications in engineering, chemistry and in biology like e.g. biological signal transduction. For the a posteriori error control of the coupled system, a residual error estimator is derived which takes into account the approximation errors due to the finite element discretisation in space as well as the polyhedral approximation of the surface. An adaptive refinement algorithm controls the overall error. Numerical experiments illustrate the performance of the a posteriori error estimator and the proposed adaptive algorithm with several benchmark examples.
  • Thumbnail Image
    Item
    Rational modeling of electrochemical double layers in thermodynamic non-equilibrium
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Dreyer, Wolfgang; Guhlke, Clemens; Müller, Rüdiger
    We consider the contact between an electrolyte and a solid electrode. At first we formulate a thermodynamic consistent model that resolves boundary layers at interfaces. The model includes charge transport, diffusion, chemical reactions, viscosity, elasticity and polarization under isothermal conditions. There is a coupling between these phenomena that particularly involves the local pressure in the electrolyte. Therefore the momentum balance is of major importance for the correct description of the layers. The width of the boundary layers is typically very small compared to the macroscopic dimensions of the system. In a second step we thus apply the method of asymptotic analysis to derive a simpler reduced model that does not resolve the boundary layers but instead incorporates the electrochemical properties of the layers into a set of new boundary conditions. For a metal-electrolyte interface, we derive a qualitative description of the double layer capacitance without the need to resolve space charge layers.
  • Thumbnail Image
    Item
    Rational modeling of electrochemical double-layers and derivation of Butler-Volmer equations
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Dreyer, Wolfgang; Guhlke, Clemens; Müller, Rüdiger
    We derive the boundary conditions for the contact between an electrolyte and a solid electrode. At first we revisit the thermodynamic consistent complete model that resolves the actual electrodeelectrolyte interface and its adjacent boundary layers. The width of these layers is controlled by the Debye length that is typically very small, leading to strongly different length scales in the system. We apply the method of asymptotic analysis to derive a simpler reduced model that does not resolve the boundary layers but instead incorporates the electrochemical properties of the layers into a set of new boundary conditions. This approach fully determines the relation of bulk quantities to the boundary conditions of the reduced model. In particular, the Butler-Volmer equations for electrochemical reactions, which are still under discussion in the literature, are rational consequences of our approach. For illustration and to compare with the literature, we consider a simple generic reaction.
  • Thumbnail Image
    Item
    Stable Crank-Nicolson discretisation for incompressible miscible displacement problems of low regularity
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Jensen, Max; Müller, Rüdiger
    In this article we study the numerical approximation of incompressible miscible displacement problems with a linearised Crank-Nicolson time discretisation, combined with a mixed finite element and discontinuous Galerkin method. At the heart of the analysis is the proof of convergence under low regularity requirements. Numerical experiments demonstrate that the proposed method exhibits second-order convergence for smooth and robustness for rough problems.
  • Thumbnail Image
    Item
    Thermodynamic models for a concentration and electric field dependent susceptibility in liquid electrolytes
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2021) Landstorfer, Manuel; Müller, Rüdiger
    The dielectric susceptibility $chi$ is an elementary quantity of the electrochemical double layer and the associated Poisson equation. While most often $chi$ is treated as a material constant, its dependency on the salt concentration in liquid electrolytes is demonstrated by various bulk electrolyte experiments. This is usually referred to as dielectric decrement. Further, it is theoretically well accepted that the susceptibility declines for large electric fields. This effect is frequently termed dielectric saturation. We analyze the impact of a variable susceptibility in terms of species concentrations and electric fields based on non-equilibrium thermodynamics. This reveals some non-obvious generalizations compared to the case of a constant susceptibility. In particular the consistent coupling of the Poisson equation, the momentum balance and the chemical potentials functions are of ultimate importance. In a numerical study, we systematically analyze the effects of a concentration and field dependent susceptibility on the double layer of a planar electrode electrolyte interface. We compute the differential capacitance and the spatial structure of the electric potential, solvent concentration and ionic distribution for various non-constant models of $chi$.