Publications
Some of the current members of the research group "Thermodynamic Modeling and Analysis of Phase Transitions" were members of a former group or of RG 1 respectively. Therefore, the corresponding publications can found on the web pages of these groups:
 former Young Scientists' Group "Modeling of Damage Processes"
 former Leibniz Group "Mathematical Models for LithiumIon Batteries"
 Research group 1 "Partial Differential Equations"

M. Dimian, P. Gurevich, O. Klein, D. Knees, D. Rachinskii, S. Tikhomirov, eds., MURPHYSHSFS2014: 7th International Workshop on MUltiRate Processes and HYSteresis (MURPHYS) & 2nd International Workshop on Hysteresis and SlowFast Systems (HSFS), 727 of Journal of Physics: Conference Series, IOP Publishing, 2016, 252 pages, (Collection Published).

E. Meca Álvarez, A. Münch, B. Wagner, Thinfilm electrodes for highcapacity lithiumion batteries: Influence of phase transformations on stress, Proceedings of The Royal Society of London. Series A. Mathematical, Physical and Engineering Sciences, 472 (2016) pp. 20160093/120160093/15.

M. Korzec, A. Münch, E. Süli, B. Wagner, Anisotropy in wavelet based phase field models, Discrete and Continuous Dynamical Systems. Series B. A Journal Bridging Mathematics and Sciences, 21 (2016) pp. 11671187.
Abstract
Anisotropy is an essential feature of phasefield models, in particular when describing the evolution of microstructures in solids. The symmetries of the crystalline phases are reflected in the interfacial energy by introducing corresponding directional dependencies in the gradient energy coefficients, which multiply the highest order derivative in the phasefield model. This paper instead considers an alternative approach, where the anisotropic gradient energy terms are replaced by a wavelet analogue that is intrinsically anisotropic and linear. In our studies we focus on the classical coupled temperature  GinzburgLandau type phasefield model for dendritic growth. For the resulting derivativefree wavelet analogue existence, uniqueness and continuous dependence on initial data for weak solutions is proved. The ability to capture dendritic growth similar to the results obtained from classical models is investigated numerically. 
W. Dreyer, C. Guhlke, M. Landstorfer, Theory and structure of the metal/electrolyte interface incorporating adsorption and solvation effects, Electrochimica Acta, 201 (2016) pp. 187219.
Abstract
In this work we present a continuum theory for the metal/electrolyte interface which explicitly takes into account adsorption and partial solvation on the metal surface. It is based on a general theory of coupled thermoelectrodynamics for volumes and surfaces, utilized here in equilibrium and a 1D approximation. We provide explicit free energy models for the volumetric metal and electrolyte phases and derive a surface free energy for the species present on the metal surface. This surface mixture theory explicitly takes into account the very different amount of sites an adsorbate requires, originating from solvation effects on the surface. Additionally we account for electron transfer reactions on the surface and the associated stripping of the solvation shell. Based on our overall surface free energy we thus provide explicit expressions of the surface chemical potentials of all constituents. The equilibrium representations of the coverages and the overall charge are briefly summarized.
Our model is then used to describe two examples: (i) a silver single crystal electrode with (100) face in contact to a (0.01M NaF + 0.01M KPF6) aqueous solution, and (ii) a general metal surface in contact to some electrolytic solution AC for which an electron transfer reaction occurs in the potential range of interest. We reflect the actual modeling procedure for these examples and discuss the respective model parameters. Due to the representations of the coverages in terms of the applied potential we provide an adsorption map and introduce adsorption potentials. Finally we investigate the structure of the space charge layer at the metal/surface/electrolyte interface by means of numerical solutions of the coupled Poissonmomentum equation system for various applied potentials. It turns out that various layers selfconsistently form within the overall space charge region, which are compared to historic and recent pictures of the double layer. Based on this we present new interpretations of what is known as inner and outer Helmholtzplanes and finally provide a thermodynamic consistent picture of the metal/electrolyte interface structure. 
W. Dreyer, C. Guhlke, R. Müller, A new perspective on the electron transfer: Recovering the ButlerVolmer equation in nonequilibrium thermodynamics, Physical Chemistry Chemical Physics, 18 (2016) pp. 2496624983.
Abstract
Understanding and correct mathematical description of electron transfer reaction is a central question in electrochemistry. Typically the electron transfer reactions are described by the ButlerVolmer equation which has its origin in kinetic theories. The ButlerVolmer equation relates interfacial reaction rates to bulk quantities like the electrostatic potential and electrolyte concentrations. Since in the classical form, the validity of the ButlerVolmer equation is limited to some simple electrochemical systems, many attempts have been made to generalize the ButlerVolmer equation. Based on nonequilibrium thermodynamics we have recently derived a reduced model for the electrodeelectrolyte interface. This reduced model includes surface reactions 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 ButlerVolmer equation. We discuss the application of the new ButlerVolmer equations to different scenarios like electron transfer reactions at metal electrodes, the intercalation process in lithiumironphosphate electrodes and adsorption processes. We illustrate the theory by an example of electroplating. 
O. Klein, A representation result for rateindependent systems, Phys. B, 486 (2016) pp. 8183.

O. Klein, V. Recupero, Hausdorff metric BV discontinuity of sweeping processes, in: MURPHYSHSFS2014: 7th International Workshop on MUltiRate Processes & HYSteresis (MURPHYS) & the 2nd International Workshop on Hysteresis and SlowFast Systems (HSFS), O. Klein, M. Dimian, P. Gurevich, D. Knees, D. Rachinskii, S. Tikhomirov, eds., 727 of Journal of Physics: Conference Series, 2016, pp. 012006/1012006/12, DOI 10.1088/17426596/727/1/012006 .
Abstract
Sweeping processes are a class of evolution differential inclusions arising in elastoplasticity and were introduced by J.J. Moreau in the early seventies. The solution operator of the sweeping processes represents a relevant example of emphrate independent operator containing as a particular case the so called emphplay operator which is widely used in hysteresis. The continuity properties of these operators were studied in several works. In this note we address the continuity with respect to the strict metric in the space of functions of bounded variation with values in the metric space of closed convex subsets of a Hilbert space. We provide a counterexample showing that the solution operator of the sweeping process is not continuous when its domain is endowed with the strict topology of $BV$ and its codomain is endowed with the $L^1$topology. This is at variance with the case of the play operator which instead is continuous in this sense. 
S. Bommer, S. Jachalski, D. Peschka, R. Seemann, B. Wagner, Structure formation in thin liquidliquid films, Preprint no. 2380, WIAS, Berlin, 2017, DOI 10.20347/WIAS.PREPRINT.2380 .
Abstract, PDF (7275 kByte)
We revisit the problem of a liquid polymer that dewets from another liquid polymer substrate with the focus on the direct comparison of results from mathematical modeling, rigorous analysis, numerical simulation and experimental investigations of rupture, dewetting dynamics and equilibrium patterns of a thin liquidliquid system. The experimental system uses as a model system a thin polystyrene (PS) / polymethylmethacrylate (PMMA) bilayer of a few hundred nm. The polymer systems allow for in situ observation of the dewetting process by atomic force microscopy (AFM) and for a precise ex situ imaging of the liquidliquid interface. In the present study, the molecular chain length of the used polymers is chosen such that the polymers can be considered as Newtonian liquids. However, by increasing the chain length, the rheological properties of the polymers can be also tuned to a viscoelastic flow behavior. The experimental results are compared with the predictions based on the thin film models. The system parameters like contact angle and surface tensions are determined from the experiments and used for a quantitative comparison. We obtain excellent agreement for transient drop shapes on their way towards equilibrium, as well as dewetting rim profiles and dewetting dynamics. 
M. Landstorfer, Boundary conditions for electrochemical interfaces, Preprint no. 2361, WIAS, Berlin, 2016, DOI 10.20347/WIAS.PREPRINT.2361 .
Abstract, PDF (1954 kByte)
Consistent boundary conditions for electrochemical interfaces, which cover double layer charging, pseudocapacitive effects and transfer reactions, are of high demand in electrochemistry and adjacent disciplines. Mathematical modeling and optimization of electrochemical systems is a strongly emerging approach to reduce cost and increase efficiency of supercapacitors, batteries, fuel cells, and electrocatalysis. However, many mathematical models which are used to describe such systems lack a real predictive value. Origin of this shortcoming is the usage of oversimplified boundary conditions. In this work we derive the boundary conditions for some general electrodeelectrolyte interface based on nonequilibrium thermodynamics for volumes and surfaces. The resulting equations are widely applicable and cover also tangential transport. The general framework is then applied to a specific material model which allows the deduction of a currentvoltage relation and thus a comparison to experimental data. Some simplified 1D examples show the range of applicability of the new approach. 
E. Meca Álvarez, V.B. Shenoy, J. Lowengrub, H2dependent attachment kinetics and shape evolution in chemical vapor deposition graphene growth, Preprint no. 2358, WIAS, Berlin, 2016, DOI 10.20347/WIAS.PREPRINT.2358 .
Abstract, PDF (1763 kByte)
Experiments on graphene growth through chemical vapor deposition (CVD) involving methane (CH4) and hydrogen (H2) gases reveal a complex shape evolution and a nonmonotonic dependence on the partial pressure of H2 (pH2). To explain these intriguing observations, we develop a microkinetic model for the stepwise decomposition of CH4 into mobile radicals and consider two possible mechanisms of attachment to graphene crystals: CH radicals to hydrogendecorated edges of the crystals and C radicals to bare crystal edges. We derive an effective mass flux and an effective kinetic coefficient, both of which depend on pH2, and incorporate these into a phase field model. The model reproduces both the nonmonotonic dependence on pH2 and the characteristic shapes of graphene crystals observed in experiments. At small pH2, growth is limited by the kinetics of attachment while at large pH2 growth is limited because the effective mass flux is small. We also derive a simple analytical model that captures the nonmonotone behavior, enables the two mechanisms of attachment to be distinguished and provides guidelines for CVD growth of defectfree 2D crystals. 
S. Bommer, R. Seemann, S. Jachalski, D. Peschka, B. Wagner, Liquidliquid dewetting: Morphologies and rates, Preprint no. 2346, WIAS, Berlin, 2016.
Abstract, PDF (1107 kByte)
The dependence of the dissipation on the local details of the flow field of a liquid polymer film dewetting from a liquid polymer substrate is shown, solving the free boundary problem for a twolayer liquid system. As a key result we show that the dewetting rates of such a liquid bilayer system can not be described by a single power law but shows transient behaviour of the rates, changing from increasing to decreasing behaviour. The theoretical predictions on the evolution of morphology and rates of the free surfaces and free interfaces are compared to measurements of the evolution of the polystyrene(PS)air, the polymethyl methacrylate (PMMA)air and the PSPMMA interfaces using in situ atomic force microscopy (AFM), and they show excellent agreement. 
M. Landstorfer, The partial molar volume and area of solvated ions and some aspects of partial charge transfer, Preprint no. 2337, WIAS, Berlin, 2016, DOI 10.20347/WIAS.PREPRINT.2337 .
Abstract, PDF (10181 kByte)
The double layer capacity is one of the central quantities in theoretical and experimental electrochemistry of metal/electrolyte interfaces. It turns out that the capacity is related to two central thermodynamic quantities, i.e. the partial molar volume of an ionic constituent and the partial molar area of the respective adsorbate. Since ions in solution (or on the surface) accumulated solvent molecules in their solvation shell, the partial molar volume and area are effected by this phenomena. In this work we discuss several aspects of the relationship between the molar volume and area of an ion, the solvation number and the charge number. In addition, we account for partial charge transfer on the metal surface which explains naturally the difference of the capacity maxima between ceF and ceClO4 on silver. We provide simple yet validated analytical expressions for the partial molar volume and area of multivalent ions and parameter values for aqueous solutions. 
M. Dziwnik, S. Jachalski, Existence of solutions to an anisotropic degenerate CahnHilliardtype equation, Preprint no. 2332, WIAS, Berlin, 2016.
Abstract, PDF (278 kByte)
We prove existence of solutions to an anisotropic CahnHilliardtype equation with degenerate diffusional mobility. In particular, the mobility vanishes at the pure phases, which is typically used to model motion by surface diffusion. The main difficulty of the present existence result is the strong nonlinearity given by the fourthorder anisotropic operator. Imposing particular assumptions on the domain and assuming that the strength of the anisotropy is sufficiently small enables to establish appropriate auxiliary results which play an essential part in the present existence proof. In addition to the existence we show that the absolute value of the corresponding solutions is bounded by 1. 
W. Dreyer, P. Friz, P. Gajewski, C. Guhlke, M. Maurelli, Stochastic model for LFPelectrodes, Preprint no. 2329, WIAS, Berlin, 2016.
Abstract, PDF (1531 kByte)
In the framework of nonequilibrium thermodynamics we derive a new model for porous electrodes. The model is applied to LiFePO4 (LFP) electrodes consisting of many LFP particles of nanometer size. The phase transition from a lithiumpoor to a lithiumrich phase within LFP electrodes is controlled by surface fluctuations leading to a system of stochastic differential equations. The model is capable to derive an explicit relation between battery voltage and current that is controlled by thermodynamic state variables. This voltagecurrent relation reveals that in thin LFP electrodes lithium intercalation from the particle surfaces into the LFP particles is the principal rate limiting process. There are only two constant kinetic parameters in the model describing the intercalation rate and the fluctuation strength, respectively. The model correctly predicts several features of LFP electrodes, viz. the phase transition, the observed voltage plateaus, hysteresis and the rate limiting capacity. Moreover we study the impact of both the particle size distribution and the active surface area on the voltagecharge characteristics of the electrode. Finally we carefully discuss the phase transition for varying charging/discharging rates. 
W. Dreyer, P.É. Druet, P. Gajewski, C. Guhlke, Existence of weak solutions for improved NernstPlanckPoisson models of compressible reacting electrolytes, Preprint no. 2291, WIAS, Berlin, 2016.
Abstract, PDF (638 kByte)
We consider an improved NernstPlanckPoisson model for compressible electrolytes first proposed by Dreyer et al. in 2013. The model takes into account the elastic deformation of the medium. In particular, large pressure contributions near electrochemical interfaces induce an inherent coupling of mass and momentum transport. The model consists of convectiondiffusionreaction equations for the constituents of the mixture, of the NavierStokes equation for the barycentric velocity and the Poisson equation for the electrical potential. Crossdiffusion phenomena occur due to the principle of mass conservation. Moreover, the diffusion matrix (mobility matrix) has a zero eigenvalue, meaning that the system is degenerate parabolic. In this paper we establish the existence of a globalin time weak solution for the full model, allowing for crossdiffusion and an arbitrary number of chemical reactions in the bulk and on the active boundary. 
E. Meca Álvarez, A. Münch, B. Wagner, Sharpinterface formation during lithium intercalation into silicon, Preprint no. 2257, WIAS, Berlin, 2016.
Abstract, PDF (518 kByte)
In this study we present a phasefield model that describes the process of intercalation of Li ions into a layer of an amorphous solid such as aSi. The governing equations couple a viscous CahnHilliardReaction model with elasticity in the framework of the CahnLarché system. We discuss the parameter settings and flux conditions at the free boundary that lead to the formation of phase boundaries having a sharp gradient in ion concentration between the initial state of the solid layer and the intercalated region. We carry out a matched asymptotic analysis to derive the corresponding sharpinterface model that also takes into account the dynamics of triple points where the sharp interface in the bulk of the layer intersects the free boundary. We numerically compare the interface motion predicted by the sharpinterface model with the longtime dynamics of the phasefield model. 
A. Münch, B. Wagner, L.P. Cook, R.R. Braun, Apparent slip for an upper convected Maxwell fluid, Preprint no. 2255, WIAS, Berlin, 2016.
Abstract, PDF (367 kByte)
In this study the flow field of a nonlocal, diffusive upper convected Maxwell (UCM) fluid with a polymer in a solvent undergoing shearing motion is investigated for pressure driven planar channel flow and the free boundary problem of a liquid layer on a solid substrate. For large ratios of the zero shear polymer viscosity to the solvent viscosity, it is shown that channel flows exhibit boundary layers at the channel walls. In addition, for increasing stress diffusion the flow field away from the boundary layers undergoes a transition from a parabolic to a plug flow. Using experimental data for the wormlike micelle solutions CTAB/NaSal and CPyCl/NaSal, it is shown that the analytic solution of the governing equations predicts these signatures of the velocity profiles. Corresponding flow structures and transitions are found for the free boundary problem of a thin layer sheared along a solid substrate. Matched asymptotic expansions are used to first derive sharpinterface models describing the bulk flow with expressions for an em apparent slip for the boundary conditions, obtained by matching to the flow in the boundary layers. For a thin film geometry several asymptotic regimes are identified in terms of the order of magnitude of the stress diffusion, and corresponding new thin film models with a slip boundary condition are derived. 
O. Klein, Uncertainty quantification for hysteresis operators and a model for magnetomechanical hysteresis, Preprint no. 2246, WIAS, Berlin, 2016.
Abstract, PDF (708 kByte)
Many models for magnetomechanical components involve hysteresis operators. The parameter within these operators have to be identified from measurements and are therefore subject to uncertainties. To quantify the influence of these uncertainties, the parameter in the hysteresis operator are considered as functions of random variables. Combining this with the hysteresis operator, we get new random variables and we can compute stochastic properties of the output of the model. For two hysteresis operators corresponding numerical results are presented in this paper. Moreover, the influence of the variation of the parameters in a model for a magnetomechanical component is investigated. 
S. Bergmann, Phasefield modeling of Si thinfilm growth via liquidphase crystallization, Winter Workshop on Microstructure Characterization and Modeling for Solar Cells, February 12  16, 2017, SchlierseeSpitzingsee, February 13, 2017.

B. Wagner, Phasefield modeling of Si thinfilm growth via liquidphase crystallization, Winter Workshop on Microstructure Characterization and Modeling for Solar Cells, February 13  16, 2017, SchlierseeSpitzingsee.

W. Dreyer, J. Fuhrmann, C. Guhlke, M. Landstorfer, M. Maurelli, R. Müller, Stochastic model for LiFePO4electrodes, ModVal14  14th Symposium on Fuel Cell and Battery Modeling and Experimental Validation, Karlsruhe, March 2  3, 2017.

S. Jachalski, Structure formation in thin liquidliquid films, SPP 1506 Kolloquium, May 9  10, 2016, Aachen, May 9, 2016.

M. Landstorfer, Towards a model based understanding of cyclic voltammetry, Electrochemistry 2016, Gesellschaft Deutscher Chemiker e.V. (GDCh), Goslar, September 26  28, 2016.

M. Maurelli , P. Gajewski, Stochastic methods for lithiumion batteries, Matheon Center Days, April 11, 2016.

E. Meca Álvarez, Si electrodes for Liion batteries: Phasefield modeling, University of Oxford, Mathematical Institute, UK, February 18, 2016.

E. Meca Alvarez, Thinfilm electrodes for highcapacity lithiumion batteries: Influence of phase transformations on stress, ECMI 2016  The 19th European Conference on Mathematics for Industry, CT04 Numerical simulation of batteries, June 13  17, 2016, Santiago de Compostela, Spain, June 17, 2016.

R. Müller, W. Dreyer, J. Fuhrmann, C. Guhlke, New insights into ButlerVolmer kinetics from thermodynamic modeling, The 67th Annual Meeting of the International Society of Electrochemistry, Den Haag, Netherlands, August 21  26, 2016.

R. Müller, Aposteriori error control for stationary coupled bulksurface equations, Joint Annual Meeting GAMM und DMV, March 7  11, 2016, TU Braunschweig, March 10, 2016.

R. Müller, Adaptive FEM and a posteriori error control for stationary coupled bulksurface PDEs, 2nd Applied Mathematical Symposium Münster, Workshop: Numerical Schemes for Surface PDEs, February 22  24, 2016, Westfälische WilhelmsUniversität Münster, February 22, 2016.

R. Müller, The Lippmann equation for liquid metal electrodes, 11th DFGCNRS Workshop ``MicroMacro Modeling and Simulation of LiquidVapor Flows'', March 2  4, 2016, Université Pierre et Marie Curie, Laboratoire JacquesLouis Lions, Paris, France, March 3, 2016.

C. Guhlke, Modeling of LiFePO4Electrode, VII European Congress on Computational Methods in Applied Sciences and Engineering  ECCOMAS 2016, Minisymposium ``Advanced Computational Modeling of Batteries and Fuel Cells'', June 5  10, 2016, Kreta, Greece, June 7, 2016.

B. Wagner, Concentrated suspensions: A multiscale analysis, University of Limerick, Mathematics Applications Consortium for Science and Industry, Ireland, April 1, 2016.

B. Wagner, Dewetting patterns of viscoelastic liquid bilayers, Workshop ``Nonequilibrium dynamics of thin films  solids, liquids and bioactiv materials'', September 20  23, 2016, CECAMHQEPFL, Lausanne, Switzerland, September 20, 2016.

B. Wagner, Models for the twophase flow of concentrated suspensions, British Applied Mathematics Colloquium, April 5  8, 2016, University of Oxford, Mathematical Institute, UK, April 6, 2016.

B. Wagner, Thin films and complex liquids for industrial application, LeibnizKolleg for Young Researchers: Challenges and Chances of Interdisciplinary Research, Thematic Workshop ``Models and Modelling'', LeibnizGemeinschaft, Berlin, November 9, 2016.

W. Dreyer, Models of lithiumionbatteries: A paradigm for consistent modelling, 1st Leibniz MMS Days, January 27  29, 2016, WIAS, Berlin, January 28, 2016.

J. Fuhrmann, W. Dreyer, C. Guhlke, M. Landstorfer, R. Müller, A. Linke, Ch. Merdon, Modeling and numerics for electrochemical systems, Micro Battery and Capacitive Energy Harvesting Materials  Results of the MatFlexEnd Project, Universität Wien, Austria, September 19, 2016.

J. Fuhrmann, A. Linke, Ch. Merdon, W. Dreyer, C. Guhle, M. Landstorfer, R. Müller, Numerical methods for electrochemical systems, 2nd Graz Battery Days, Graz, Austria, September 27  28, 2016.

C. Guhlke, W. Dreyer, R. Müller, M. Landstorfer, J. Fuhrmann, Beyond Newman's battery model, 2nd Graz Battery Days, Graz, Austria, September 27  28, 2016.

C. Guhlke, W. Dreyer, R. Müller, M. Landstorfer , Models of lithiumion batteries: A paradigm for consistent modelling, 1st Leibniz MMS Days, WIAS Berlin, January 27  29, 2016.

C. Guhlke, J. Fuhrmann, W. Dreyer, R. Müller, M. Landstorfer, Modeling of batteries, Batterieforum Deutschland 2016, Berlin, April 6  8, 2016.

C. Guhlke, Continuum thermodynamics in the immediate vicinity of an electrodeelectrolyte interface, Soft Matter at Interfaces 2016, February 28  March 2, 2016, Schloss Ringberg, Kreuth, February 28, 2016.

C. Guhlke, Impact of stress on the electrical double layer, 2016 MRS Spring Meeting & Exhibit, Symposium EE7 ``Mechanics of Energy Storage and Conversion  Batteries, Thermoelectrics and Fuel Cells'', March 28  April 1, 2016, Phoenix, USA, March 31, 2016.

O. Klein, On uncertainty quantification for hysteresis operators, Silesian University, Mathematical Institute, Opava, Czech Republic, December 7, 2016.

O. Klein, Uncertainty quantification for hysteresis operators, 7th European Congress of Mathematics (7ECM), Minisymposium 29 ``Nonsmooth PDEs in the Modeling Damage, Delamination, and Fracture'', July 18  22, 2016, Technische Universität Berlin, July 22, 2016.

O. Klein, Uncertainty quantification for hysteresis operators and models for magnetomechanical hysteresis, Conference ``Advances in Magnetics'' (AIM) 2016, March 14  16, 2016, Bormio, Italy, March 14, 2016.
Monographs
Articles in Refereed Journals
Contributions to Collected Editions
Preprints, Reports, Technical Reports
Talks, Poster
Research Groups
 Partial Differential Equations
 Laser Dynamics
 Numerical Mathematics and Scientific Computing
 Nonlinear Optimization and Inverse Problems
 Interacting Random Systems
 Stochastic Algorithms and Nonparametric Statistics
 Thermodynamic Modeling and Analysis of Phase Transitions
 Nonsmooth Variational Problems and Operator Equations