14 results
Search Results
Now showing 1 - 10 of 14
- ItemOptimal 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üdigerOptimal 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.
- ItemError 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üdigerA 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.
- ItemOvercoming the shortcomings of the Nernst-Planck model(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2012) Dreyer, Wolfgang; Guhlke, Clemens; Müller, RüdigerThis 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.
- ItemStable 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üdigerIn 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.
- ItemA 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üdigerUnderstanding 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.
- ItemThermodynamic 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üdigerThe 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$.
- ItemModeling 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üdigerWe 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.
- ItemThe 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üdigerElectro-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.
- ItemRational 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üdigerWe 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.
- ItemRational 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üdigerWe 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.