Thin Film Equations
The mathematical modeling of the dynamics of thin liquid as well as thin solid films typically uses the small ratio of vertical to horizontal scales of the evolving morphologies to simplify the underlying free boundary problems to dimensionreduced partial differential equations for evolving free surface and/or interfaces of the thin films.
Methods of singular perturbation theory show that the resulting partial differential equations are typically of high order, such as lubrication type equations, or systems thereof, or convective CahnHilliard type equations.
To understand the complex nonlinear solution structure of these equations and to complement numerical simulation, matched asymptotic expansions are being derived and exponential asymptotics are extended for these higher order equations. Both, formal and rigorous asymptotic analysis lead to corresponding sharpinterface models that are used to systematically describe important properties, such as the stability of moving contact lines or the longtime behaviour of coarsening binary alloys, the longtime convergence to the selfsimilar approach towards finitetime blowup, describing the process of film rupture, or the discovery of novel stationary solutions describing the facetting of socalled quantum dots.
Sharp interface limits of generalized NavierStokesKorteweg systems
Sharp interface limits are considered for isothermal and nonisothermal Korteweg and NavierStokesKorteweg phase field models.
The sharp interface limit is obtained by matched asymptotic expansions of the phase fields in powers of the interface width ε. These expansions are considered in the interfacial region (inner expansion) and in the bulk (outer expansion), and are matched order by order. This results in partial differential equations for the diffuse fields and a series of boundary conditions at the interface.
For different scalings, solvability criteria for the inner equations are established. This leads to various settings of sharp interface models with different jump conditions at the sharp interface.
Sharp limits of regularized diffusion equations with mechanical coupling for applications in energy technology
Two futureoriented energy storage systems, i.e. the storage of electrical energy in lithiumion batteries and the storage of hydrogen in metal hydrides, are studied. In both applications foreign atoms, lithium respectively hydrogen, will reversibly be stored in crystals. During loading and unloading the foreign atoms form a two phase system with high and low concentration. The phases are separated by a moving interface. Moreover, the storage of foreign atoms in the crystals is accompanied by a volume change leading to mechanical stresses. The physical behavior of the storage process is described by a sharpinterface model.
The sharpinterface model consists of a system of parabolic and elliptic differential equations. The evolution of the interface is described by algebraic and ordinary differential equations, which are called jump conditions.
The regularization of the sharpinterface model by higher gradients and viscous terms leads to local phase field models. Here the jump of the concentration between the two phases is a smooth transition of thickness ε. Therefore the jump conditions at the interface are redundant. However, in the sharp limit ε to zero, the jump conditions of a sharpinterface model must be recovered.
By means of formal asymptotic methods several regularizations of the sharpinterface model are studied. One could show, that the regularization of the diffusion equation implies a regularization of the balance of momentum. Furthermore there are regularizations that lead to wrong jump conditions, so that they can not be used to describe the energy storage systems.
Publications
Articles in Refereed Journals

M. Heida, S. Nesenenko, Stochastic homogenization of ratedependent models of monotone type in plasticity, Asymptotic Analysis, 112 (2019), pp. 185212, DOI 10.3233/ASY181502 .
Abstract
In this work we deal with the stochastic homogenization of the initial boundary value problems of monotone type. The models of monotone type under consideration describe the deformation behaviour of inelastic materials with a microstructure which can be characterised by random measures. Based on the Fitzpatrick function concept we reduce the study of the asymptotic behaviour of monotone operators associated with our models to the problem of the stochastic homogenization of convex functionals within an ergodic and stationary setting. The concept of Fitzpatrick's function helps us to introduce and show the existence of the weak solutions for ratedependent systems. The derivations of the homogenization results presented in this work are based on the stochastic twoscale convergence in Sobolev spaces. For completeness, we also present some twoscale homogenization results for convex functionals, which are related to the classical Gammaconvergence theory. 
M. Heida, B. Schweizer, Stochastic homogenization of plasticity equations, ESAIM. Control, Optimisation and Calculus of Variations, 24 (2018), pp. 153176.
Abstract
In the context of infinitesimal strain plasticity with hardening, we derive a stochastic homogenization result. We assume that the coefficients of the equation are random functions: elasticity tensor, hardening parameter and flowrule function are given through a dynamical system on a probability space. A parameter ε > 0 denotes the typical length scale of oscillations. We derive effective equations that describe the behavior of solutions in the limit ε → 0. The homogenization procedure is based on the fact that stochastic coefficients “allow averaging”: For one representative volume element, a strain evolution [0,T]∋ t ↦ξ(t) ∈ ℝ_{s}^{dxd} induces a stress evolution [0,T]∋ t ↦Σ (ξ) (t)∈ℝ_{s}^{dxd}. Once the hysteretic evolution law Σ is justified for averages, we obtain that the macroscopic limit equation is given by ∇ ⋅ Σ(∇^{s}u)=f. 
M. Heida, On convergences of the squareroot approximation scheme to the FokkerPlanck operator, Mathematical Models & Methods in Applied Sciences, 28 (2018), pp. 25992635, DOI 10.1142/S0218202518500562 .
Abstract
We study the qualitative convergence properties of a finite volume scheme that recently was proposed by Lie, Fackeldey and Weber [SIAM Journal on Matrix Analysis and Applications 2013 (34/2)] in the context of conformation dynamics. The scheme was derived from physical principles and is called the squareroot approximation (SQRA) scheme. We show that solutions to the SQRA equation converge to solutions of the FokkerPlanck equation using a discrete notion of Gconvergence. Hence the squareroot approximation turns out to be a usefull approximation scheme to the FokkerPlanck equation in high dimensional spaces. As an example, in the special case of stationary Voronoi tessellations we use stochastic twoscale convergence to prove that this setting satisfies the Gconvergence property. In particular, the class of tessellations for which the Gconvergence result holds is not trivial. 
E. Meca Álvarez, A. Münch, B. Wagner, Localized instabilities and spinodal decomposition in driven systems in the presence of elasticity, Physical Review E. Statistical, Nonlinear, and Soft Matter Physics, 97 (2018), pp. 012801/1012801/12, DOI 10.1103/PhysRevE.97.012801 .
Abstract
We study numerically and analytically the instabilities associated with phase separation in a solid layer on which an external material flux is imposed. The first instability is localized within a boundary layer at the exposed free surface by a process akin to spinodal decomposition. In the limiting static case, when there is no material flux, the coherent spinodal decomposition is recovered. In the present problem stability analysis of the timedependent and nonuniform base states as well as numerical simulations of the full governing equations are used to establish the dependence of the wavelength and onset of the instability on parameter settings and its transient nature as the patterns eventually coarsen into a flat moving front. The second instability is related to the MullinsSekerka instability in the presence of elasticity and arises at the moving front between the two phases when the flux is reversed. Stability analyses of the full model and the corresponding sharpinterface model are carried out and compared. Our results demonstrate how interface and bulk instabilities can be analysed within the same framework which allows to identify and distinguish each of them clearly. The relevance for a detailed understanding of both instabilities and their interconnections in a realistic setting are demonstrated for a system of equations modelling the lithiation/delithiation processes within the context of Lithium ion batteries. 
T. Ahnert, A. Münch, B. Wagner, Models for the twophase flow of concentrated suspensions, European Journal of Applied Mathematics, 30 (2019), pp. 585617 (published online on 04.06.2018), DOI 10.1017/S095679251800030X .
Abstract
A new twophase model is derived that make use of a constitutive law combining nonBrownian suspension with granular rheology, that was recently proposed by Boyer et al. [PRL, 107(18),188301 (2011)]. It is shown that for the simple channel flow geometry, the stress model naturally exhibits a Bingham type flow property with an unyielded finitesize zone in the center of the channel. As the volume fraction of the solid phase is increased, the various transitions in the flow fields are discussed using phase space methods for a boundary value problem, that is derived from the full model. The predictions of this analysis is then compared to the direct finiteelement numerical solutions of the full model. 
L. Donati, M. Heida, M. Weber, B. Keller, Estimation of the infinitesimal generator by squareroot approximation, Journal of Physics: Condensed Matter, 30 (2018), pp. 425201/1425201/14, DOI 10.1088/1361648X/aadfc8 .
Abstract
For the analysis of molecular processes, the estimation of timescales, i.e., transition rates, is very important. Estimating the transition rates between molecular conformations is  from a mathematical point of view  an invariant subspace projection problem. A certain infinitesimal generator acting on function space is projected to a lowdimensional rate matrix. This projection can be performed in two steps. First, the infinitesimal generator is discretized, then the invariant subspace is approximated and used for the subspace projection. In our approach, the discretization will be based on a Voronoi tessellation of the conformational space. We will show that the discretized infinitesimal generator can simply be approximated by the geometric average of the Boltzmann weights of the Voronoi cells. Thus, there is a direct correlation between the potential energy surface of molecular structures and the transition rates of conformational changes. We present results for a 2ddiffusion process and Alanine dipeptide. 
W. Dreyer, C. Guhlke, R. Müller, Bulksurface electrothermodynamics and applications to electrochemistry, Entropy. An International and Interdisciplinary Journal of Entropy and Information Studies, 20 (2018), pp. 939/1939/44, DOI 10.3390/e20120939 .
Abstract
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. 
M. Liero, S. Reichelt, Homogenization of CahnHilliardtype equations via evolutionary Gammaconvergence, NoDEA. Nonlinear Differential Equations and Applications, 25 (2018), pp. 6/16/31, DOI 10.1007/s0003001804959 .
Abstract
In this paper we discuss two approaches to evolutionary Γconvergence of gradient systems in Hilbert spaces. The formulation of the gradient system is based on two functionals, namely the energy functional and the dissipation potential, which allows us to employ Γconvergence methods. In the first approach we consider families of uniformly convex energy functionals such that the limit passage of the timedependent problems can be based on the theory of evolutionary variational inequalities as developed by Daneri and Savaré 2010. The second approach uses the equivalent formulation of the gradient system via the energydissipation principle and follows the ideas of Sandier and Serfaty 2004. We apply both approaches to rigorously derive homogenization limits for CahnHilliardtype equations. Using the method of weak and strong twoscale convergence via periodic unfolding, we show that the energy and dissipation functionals Γconverge. In conclusion, we will give specific examples for the applicability of each of the two approaches. 
S. Bergmann, D.A. BarraganYani, E. Flegel, K. Albe, B. Wagner, Anisotropic solidliquid interface kinetics in silicon: An atomistically informed phasefield model, Modelling and Simulation in Materials Science and Engineering, 25 (2017), pp. 065015/1065015/20, DOI 10.1088/1361651X/aa7862 .
Abstract
We present an atomistically informed parametrization of a phasefield model for describing the anisotropic mobility of liquidsolid interfaces in silicon. The model is derived from a consistent set of atomistic data and thus allows to directly link molecular dynamics and phase field simulations. Expressions for the free energy density, the interfacial energy and the temperature and orientation dependent interface mobility are systematically fitted to data from molecular dynamics simulations based on the StillingerWeber interatomic potential. The temperaturedependent interface velocity follows a VogelFulcher type behavior and allows to properly account for the dynamics in the undercooled melt. 
M. Heida, A. Mielke, Averaging of timeperiodic dissipation potentials in rateindependent processes, Discrete and Continuous Dynamical Systems  Series S, 10 (2017), pp. 13031327.
Abstract
We study the existence and wellposedness of rateindependent systems (or hysteresis operators) with a dissipation potential that oscillates in time with period ε. In particular, for the case of quadratic energies in a Hilbert space, we study the averaging limit ε→0 and show that the effctive dissipation potential is given by the minimum of all friction thresholds in one period, more precisely as the intersection of all the characteristic domains. We show that the rates of the process do not converge weakly, hence our analysis uses the notion of energetic solutions and relies on a detailed estimates to obtain a suitable equicontinuity of the solutions in the limit ε→0. 
M. Heida, Stochastic homogenization of rateindependent systems, Continuum Mechanics and Thermodynamics, 29 (2017), pp. 853894, DOI 10.1007/s001610170564z .
Abstract
We study the stochastic and periodic homogenization 1homogeneous convex functionals. We proof some convergence results with respect to stochastic twoscale convergence, which are related to classical Gammaconvergence results. The main result is a general liminfestimate for a sequence of 1homogeneous functionals and a twoscale stability result for sequences of convex sets. We apply our results to the homogenization of rateindependent systems with 1homogeneous dissipation potentials and quadratic energies. In these applications, both the energy and the dissipation potential have an underlying stochastic microscopic structure. We study the particular homogenization problems of PrandltReuss plasticity, Coulomb friction on a macroscopic surface and Coulomb friction on microscopic fissure. 
M. Liero, A. Mielke, M.A. Peletier, D.R.M. Renger, On microscopic origins of generalized gradient structures, Discrete and Continuous Dynamical Systems  Series S, 10 (2017), pp. 135, DOI 10.3934/dcdss.2017001 .
Abstract
Classical gradient systems have a linear relation between rates and driving forces. In generalized gradient systems we allow for arbitrary relations derived from general nonquadratic dissipation potentials. This paper describes two natural origins for these structures. A first microscopic origin of generalized gradient structures is given by the theory of largedeviation principles. While Markovian diffusion processes lead to classical gradient structures, Poissonian jump processes give rise to coshtype dissipation potentials. A second origin arises via a new form of convergence, that we call EDPconvergence. Even when starting with classical gradient systems, where the dissipation potential is a quadratic functional of the rate, we may obtain a generalized gradient system in the evolutionary Gammalimit. As examples we treat (i) the limit of a diffusion equation having a thin layer of low diffusivity, which leads to a membrane model, and (ii) the limit of diffusion over a high barrier, which gives a reactiondiffusion system. 
M. Dziwnik, A. Münch, B. Wagner, An anisotropic phasefield model for solidstate dewetting and its sharpinterface limit, Nonlinearity, 30 (2017), pp. 14651496.
Abstract
We propose a phase field model for solid state dewetting in form of a CahnHilliard equation with weakly anisotropic surface energy and a degenerate mobility together with a free boundary condition at the filmsubstrate contact line. We derive the corresponding sharp interface limit via matched asymptotic analysis involving multiple inner layers. The resulting sharp interface model is consistent with the pure surface diffusion model. In addition, we show that the natural boundary conditions, as indicated from the first variation of the total free energy, imply a contact angle condition for the dewetting front, which, in the isotropic case, is consistent with the wellknown Young's equation. 
A. Münch, B. Wagner, L.P. Cook, R.R. Braun, Apparent slip for an upper convected Maxwell fluid, SIAM Journal on Applied Mathematics, 77 (2017), pp. 537564, DOI 10.1137/16M1056869 .
Abstract
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. 
A. Caiazzo, F. Caforio, G. Montecinos, L.O. Müller, P.J. Blanco, E.F. Toro, Assessment of reduced order Kalman filter for parameter identification in onedimensional blood flow models using experimental data, International Journal of Numerical Methods in Biomedical Engineering, 33 (2017), pp. e2843/1e2843/26, DOI 10.1002/cnm.2843 .
Abstract
This work presents a detailed investigation of a parameter estimation approach based on the reduced order unscented Kalman filter (ROUKF) in the context of onedimensional blood flow models. In particular, the main aims of this study are (i) to investigate the effect of using real measurements vs. synthetic data (i.e., numerical results of the same in silico model, perturbed with white noise) for the estimation and (ii) to identify potential difficulties and limitations of the approach in clinically realistic applications in order to assess the applicability of the filter to such setups. For these purposes, our numerical study is based on the in vitro model of the arterial network described by [Alastruey et al. 2011, J. Biomech. bf 44], for which experimental flow and pressure measurements are available at few selected locations. In order to mimic clinically relevant situations, we focus on the estimation of terminal resistances and arterial wall parameters related to vessel mechanics (Young's modulus and thickness) using few experimental observations (at most a single pressure or flow measurement per vessel). In all cases, we first perform a theoretical identifiability analysis based on the generalized sensitivity function, comparing then the results obtained with the ROUKF, using either synthetic or experimental data, to results obtained using reference parameters and to available measurements. 
W. Dreyer, C. Guhlke, M. Landstorfer, R. Müller, New insights on the interfacial tension of electrochemical interfaces and the Lippmann equation, European Journal of Applied Mathematics, 29 (2018), pp. 708753, DOI 10.1017/S0956792517000341 .
Abstract
The Lippmann equation is considered as universal relationship between interfacial tension, double layer charge, and cell potential. Based on the framework of continuum thermoelectrodynamics 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 mercuryelectrolyte 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 curvature radius is comparable to the Debye length. 
E. Meca Álvarez, A. Münch, B. Wagner, Sharpinterface formation during lithium intercalation into silicon, European Journal of Applied Mathematics, 29 (2018), pp. 118145, DOI 10.1017/S0956792517000067 .
Abstract
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. 
M. Heida, B. Schweizer, Nonperiodic homogenization of infinitesimal strain plasticity equations, ZAMM. Zeitschrift für Angewandte Mathematik und Mechanik, 96 (2016), pp. 523.

M. Khodayari, P. Reinsberg, A.A. AbdElLatif, Ch. Merdon, J. Fuhrmann, H. Baltruschat, Determining solubility and diffusivity by using a flow cell coupled to a mass spectrometer, ChemPhysChem, 17 (2016), pp. 16471655.

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, DOI 10.1039/C6CP04142F .
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. 
K. Disser, M. Liero, On gradient structures for Markov chains and the passage to Wasserstein gradient flows, Networks and Heterogeneous Media, 10 (2015), pp. 233253.
Abstract
We study the approximation of Wasserstein gradient structures by their finitedimensional analog. We show that simple finitevolume discretizations of the linear FokkerPlanck equation exhibit the recently established entropic gradientflow structure for reversible Markov chains. Then, we reprove the convergence of the discrete scheme in the limit of vanishing mesh size using only the involved gradientflow structures. In particular, we make no use of the linearity of the equations nor of the fact that the FokkerPlanck equation is of second order. 
D. Peschka, Thinfilm free boundary problems for partial wetting, Journal of Computational Physics, 295 (2015), pp. 770778.
Abstract
We present a novel framework to solve thinfilm equations with an explicit nonzero contact angle, where the support of the solution is treated as an unknown. The algorithm uses a finite element method based on a gradient formulation of the thinfilm equations coupled to an arbitrary LagrangianEulerian method for the motion of the support. Features of this algorithm are its simplicity and robustness. We apply this algorithm in 1D and 2D to problems with surface tension, contact angles and with gravity. 
P. Colli, G. Gilardi, E. Rocca, J. Sprekels, Vanishing viscosities and error estimate for a CahnHilliard type phase field system related to tumor growth, Nonlinear Analysis. Real World Applications. An International Multidisciplinary Journal, 26 (2015), pp. 93108.
Abstract
In this paper we perform an asymptotic analysis for two different vanishing viscosity coefficients occurring in a phase field system of CahnHilliard type that was recently introduced in order to approximate a tumor growth model. In particular, we extend some recent results obtained in [ColliGilardiHilhorst 2015], letting the two positive viscosity parameters tend to zero independently from each other and weakening the conditions on the initial data in such a way as to maintain the nonlinearities of the PDE system as general as possible. Finally, under proper growth conditions on the interaction potential, we prove an error estimate leading also to the uniqueness result for the limit system. 
M.G. Hennessy, V.M. Burlakov, A. Münch, B. Wagner, A. Goriely, Controlled topological transitions in thinfilm phase separation, SIAM Journal on Applied Mathematics, 75 (2015), pp. 3860.
Abstract
In this paper the evolution of a binary mixture in a thinfilm geometry with a wall at the top and bottom is considered. Bringing the mixture into its miscibility gap so that no spinodal decomposition occurs in the bulk, a slight energetic bias of the walls towards each one of the constituents ensures the nucleation of thin boundary layers that grow until the constituents have moved into one of the two layers. These layers are separated by an interfacial region, where the composition changes rapidly. Conditions that ensure the separation into two layers with a thin interfacial region are investigated based on a phasefield model and using matched asymptotic expansions a corresponding sharpinterface problem for the location of the interface is established. It is then argued that a thus created twolayer system is not at its energetic minimum but destabilizes into a controlled selfreplicating pattern of trapezoidal vertical stripes by minimizing the interfacial energy between the phases while conserving their area. A quantitative analysis of this mechanism is carried out via a new thinfilm model for the free interfaces, which is derived asymptotically from the sharpinterface model. 
R. Huth, S. Jachalski, G. Kitavtsev, D. Peschka, Gradient flow perspective on thinfilm bilayer flows, Journal of Engineering Mathematics, 94 (2015), pp. 4361.
Abstract
We study gradient flow formulations of thinfilm bilayer flows with triplejunctions between liquid/liquid/air. First we highlight the gradient structure in the Stokes freeboundary flow and identify its solutions with the well known PDE with boundary conditions. Next we propose a similar gradient formulation for the corresponding thinfilm model and formally identify solutions with those of the corresponding freeboundary problem. A robust numerical algorithm for the thinfilm gradient flow structure is then provided. Using this algorithm we compare the sharp triplejunction model with precursor models. For their stationary solutions a rigorous connection is established using Gammaconvergence. For timedependent solutions the comparison of numerical solutions shows a good agreement for small and moderate times. Finally we study spreading in the zerocontact angle case, where we compare numerical solutions with asymptotically exact sourcetype solutions. 
CH. Bayer, H.A. Hoel, A. Kadir, P. Plechac, M. Sandberg, A. Szepessy, Computational error estimates for BornOppenheimer molecular dynamics with nearly crossing potential surfaces, Applied Mathematics Research Express, 2015 (2015), pp. 329417.

A. Caiazzo, I. RamisConde, Multiscale modeling of palisade formation in glioblastoma multiforme, Journal of Theoretical Biology, 383 (2015), pp. 145156.
Abstract
Palisades are characteristic tissue aberrations that arise in glioblastomas. Observation of palisades is considered as a clinical indicator of the transition from a noninvasive to an invasive tumour. In this article we propose a computational model to study the influence of genotypic and phenotypic heterogeneity in palisade formation. For this we produced three dimensional realistic simulations, based on a multiscale hybrid model, coupling the evolution of tumour cells and the oxygen diffusion in tissue, that depict the shape of palisades during its formation. Our results can be summarized as the following: (1) we show that cell heterogeneity is a crucial factor in palisade formation and tumour growth; (2) we present results that can explain the observed fact that recursive tumours are more malignant than primary tumours; and (3) the presented simulations can provide to clinicians and biologists for a better understanding of palisades 3D structure as well as glioblastomas growth dynamics 
W. Dreyer, C. Guhlke, R. Müller, Modeling of electrochemical double layers in thermodynamic nonequilibrium, Physical Chemistry Chemical Physics, 17 (2015), pp. 2717627194, DOI 10.1039/C5CP03836G .
Abstract
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 metalelectrolyte interface, we derive a qualitative description of the double layer capacitance without the need to resolve space charge layers. 
H. Hanke, D. Knees, Homogenization of elliptic systems with nonperiodic, state dependent coefficients, Asymptotic Analysis, 92 (2015), pp. 203234.
Abstract
In this paper, a homogenization problem for an elliptic system with nonperiodic, state dependent coefficients representing microstructure is investigated. The state functions defining the tensor of coefficients are assumed to have an intrinsic length scale denoted by ε > 0. The aim is the derivation of an effective model by investigating the limit process ε → 0 of the state functions rigorously. The effective model is independent of the parameter ε > 0 but preserves the microscopic structure of the state functions (ε > 0), meaning that the effective tensor is given by a unit cell problem prescribed by a suitable microscopic tensor. Due to the nonperiodic structure of the state functions and the corresponding microstructure, the effective tensor turns out to vary from point to point (in contrast to a periodic microscopic model). In a forthcoming paper, these states will be solutions of an additional evolution law describing changes of the microstructure. Such changes could be the consequences of temperature changes, phase separation or damage progression, for instance. Here, in addition to the above and as a preparation for an application to timedependent damage models (discussed in a future paper), we provide a Γconvergence result of sequences of functionals being related to the previous microscopic models with state dependent coefficients. This requires a penalization term for piecewise constant state functions that allows us to extract from bounded sequences those sequences converging to a Sobolev function in some sense. The construction of the penalization term is inspired by techniques for Discontinuous Galerkin methods and is of own interest. A compactness and a density result are provided. 
S. Jachalski, G. Kitavtsev, R. Taranets, Weak solutions to lubrication systems describing the evolution of bilayer thin films, Communications in Mathematical Sciences, 12 (2014), pp. 527544.
Abstract
The existence of global nonnegative weak solutions is proved for coupled onedimen sional lubrication systems that describe the evolution of nanoscopic bilayer thin polymer films that take account of Navierslip or noslip conditions at both liquidliquid and liquidsolid interfaces. In addition, in the presence of attractive van der Waals and repulsive Born intermolecular interactions existence of positive smooth solutions is shown. 
G. Aki, W. Dreyer, J. Giesselmann, Ch. Kraus, A quasiincompressible diffuse interface model with phase transition, Mathematical Models & Methods in Applied Sciences, 24 (2014), pp. 827861.
Abstract
This work introduces a new thermodynamically consistent diffuse model for twocomponent flows of incompressible fluids. For the introduced diffuse interface model, we investigate physically admissible sharp interface limits by matched asymptotic techniques. To this end, we consider two scaling regimes where in one case we recover the Euler equations and in the other case the NavierStokes equations in the bulk phases equipped with admissible interfacial conditions. For the NavierStokes regime, we further assume the densities of the fluids are close to each other in the sense of a small parameter which is related to the interfacial thickness of the diffuse model. 
P. Colli, G. Gilardi, P. Krejčí, J. Sprekels, A vanishing diffusion limit in a nonstandard system of phase field equations, Evolution Equations and Control Theory, 3 (2014), pp. 257275.
Abstract
We are concerned with a nonstandard phase field model of CahnHilliard type. The model, which was introduced by PodioGuidugli (Ric. Mat. 2006), describes twospecies phase segregation and consists of a system of two highly nonlinearly coupled PDEs. It has been recently investigated by Colli, Gilardi, PodioGuidugli, and Sprekels in a series of papers: see, in particular, SIAM J. Appl. Math. 2011, and Boll. Unione Mat. Ital. 2012. In the latter contribution, the authors can treat the very general case in which the diffusivity coefficient of the parabolic PDE is allowed to depend nonlinearly on both variables. In the same framework, this paper investigates the asymptotic limit of the solutions to the initialboundary value problems as the diffusion coefficient $sigma$ in the equation governing the evolution of the order parameter tends to zero. We prove that such a limit actually exists and solves the limit problem, which couples a nonlinear PDE of parabolic type with an ODE accounting for the phase dynamics. In the case of a constant diffusivity, we are able to show uniqueness and to improve the regularity of the solution. 
A. Caiazzo, J. Mura, Multiscale modeling of weakly compressible elastic materials in harmonic regime and application to microscale structure estimation, Multiscale Modeling & Simulation. A SIAM Interdisciplinary Journal, 12 (2014), pp. 514537.
Abstract
This article is devoted to the modeling of elastic materials composed by an incompressible elastic matrix and small compressible gaseous inclusions, under a time harmonic excitation. In a biomedical context, this model describes the dynamics of a biological tissue (e.g. lung or liver) when wave analysis methods (such as Magnetic Resonance Elastography) are used to estimate tissue properties. Due to the multiscale nature of the problem, direct numerical simulations are prohibitive. We extend the homogenized model introduced in [Baffico, Grandmont, Maday, Osses, SIAM J. Mult. Mod. Sim., 7(1), 2008] to a time harmonic regime to describe the solidgas mixture from a macroscopic point of view in terms of an effective elasticity tensor. Furthermore, we derive and validate numerically analytical approximations for the effective elastic coefficients in terms of macroscopic parameters. This simplified description is used to to set up an efficient variational approach for the estimation of the tissue porosity, using the mechanical response to external harmonic excitations. 
W. Dreyer, J. Giesselmann, Ch. Kraus, A compressible mixture model with phase transition, Physica D. Nonlinear Phenomena, 273274 (2014), pp. 113.
Abstract
We introduce a new thermodynamically consistent diffuse interface model of AllenCahn/NavierStokes type for multicomponent flows with phase transitions and chemical reactions. For the introduced diffuse interface model, we investigate physically admissible sharp interface limits by matched asymptotic techniques. We consider two scaling regimes, i.e. a nondissipative and a dissipative regime, where we recover in the sharp interface limit a generalized AllenCahn/Euler system for mixtures with chemical reactions in the bulk phases equipped with admissible interfacial conditions. The interfacial conditions satify, for instance, a YoungLaplace and a Stefan type law. 
S. Neukamm, H. Olbermann, Homogenization of the nonlinear bending theory for plates, Calculus of Variations and Partial Differential Equations, (published online on Sept. 14, 2014), DOI 10.1007/s0052601407652 .
Abstract
We carry out the spatially periodic homogenization of nonlinear bending theory for plates. The derivation is rigorous in the sense of Gammaconvergence. In contrast to what one naturally would expect, our result shows that the limiting functional is not simply a quadratic functional of the second fundamental form of the deformed plate as it is the case in nonlinear plate theory. It turns out that the limiting functional discriminates between whether the deformed plate is locally shaped like a "cylinder" or not. For the derivation we investigate the oscillatory behavior of sequences of second fundamental forms associated with isometric immersions, using twoscale convergence. This is a nontrivial task, since one has to treat twoscale convergence in connection with a nonlinear differential constraint. 
P. Colli, G. Gilardi, P. PodioGuidugli, J. Sprekels, An asymptotic analysis for a nonstandard CahnHilliard system with viscosity, Discrete and Continuous Dynamical Systems  Series S, 6 (2013), pp. 353368.
Abstract
This paper is concerned with a diffusion model of phasefield type, consisting of a parabolic system of two partial differential equations, interpreted as balances of microforces and microenergy, for two unknowns: the problem's order parameter $rho$ and the chemical potential $mu$; each equation includes a viscosity term  respectively, $varepsilon,partial_tmu$ and $delta,partial_trho$  with $varepsilon$ and $delta$ two positive parameters; the field equations are complemented by Neumann homogeneous boundary conditions and suitable initial conditions. In a recent paper [5], we proved that this problem is wellposed and investigated the longtime behavior of its $(varepsilon,delta)$solutions. Here we discuss the asymptotic limit of the system as $eps$ tends to 0. We prove convergence of $(varepsilon,delta)$solutions to the corresponding solutions for the case $eps$ =0, whose longtime behavior we characterize; in the proofs, we employ compactness and monotonicity arguments. 
P. Colli, G. Gilardi, P. PodioGuidugli, J. Sprekels, Global existence and uniqueness for a singular/degenerate CahnHilliard system with viscosity, Journal of Differential Equations, 254 (2013), pp. 42174244.
Abstract
Existence and uniqueness are investigated for a nonlinear diffusion problem of phasefield type, consisting of a parabolic system of two partial differential equations, complemented by Neumann homogeneous boundary conditions and initial conditions. This system aims to model twospecies phase segregation on an atomic [19]; in the balance equations of microforces and microenergy, the two unknowns are the order parameter $rho$ and the chemical potential $mu$. A simpler version of the same system has recently been discussed in [8]. In this paper, a fairly more general phasefield equation for $rho$ is coupled with a genuinely nonlinear diffusion equation for $mu$. The existence of a globalintime solution is proved with the help of suitable a priori estimates. In the case of costant atom mobility, a new and rather unusual uniqueness proof is given, based on a suitable combination of variables. 
P. Colli, G. Gilardi, P. PodioGuidugli, J. Sprekels, Continuous dependence for a nonstandard CahnHilliard system with nonlinear atom mobility, Rendiconti del Seminario Matematico. Universita e Politecnico Torino, 70 (2012), pp. 2752.
Abstract
This note is concerned with a nonlinear diffusion problem of phasefield type, consisting of a parabolic system of two partial differential equations, complemented by Neumann homogeneous boundary conditions and initial conditions. The system arises from a model of twospecies phase segregation on an atomic lattice [PodioGuidugli 2006]; it consists of the balance equations of microforces and microenergy; the two unknowns are the order parameter $rho$ and the chemical potential $mu$. Some recent results obtained for this class of problems is reviewed and, in the case of a nonconstant and nonlinear atom mobility, uniqueness and continuous dependence on the initial data are shown with the help of a new line of argumentation developed in Colli/Gilardi/PodioGuidugli/Sprekels 2012. 
P. Colli, G. Gilardi, P. PodioGuidugli, J. Sprekels, Global existence for a strongly coupled CahnHilliard system with viscosity, Bollettino della Unione Matematica Italiana. Serie 9, 5 (2012), pp. 495513.
Abstract
An existence result is proved for a nonlinear diffusion problem of phasefield type, consisting of a parabolic system of two partial differential equations, complemented by Neumann homogeneous boundary conditions and initial conditions. This system is meant to model twospecies phase segregation on an atomic lattice under the presence of diffusion. A similar system has been recently introduced and analyzed in [CGPS11]. Both systems conform to the general theory developed in [Pod06]: two parabolic PDEs, interpreted as balances of microforces and microenergy, are to be solved for the order parameter $rho$ and the chemical potential $mu$. In the system studied in this note, a phasefield equation in $rho$ fairly more general than in [CGPS11] is coupled with a highly nonlinear diffusion equation for $mu$, in which the conductivity coefficient is allowed to depend nonlinearly on both variables. 
G. Aki, J. Daube, W. Dreyer, J. Giesselmann, M. Kränkel, Ch. Kraus, A diffuse interface model for quasiincompressible flows: Sharp interface limits and numerics, ESAIM Proceedings, 38 (2012), pp. 5477.
Abstract
In this contribution, we investigate a diffuse interface model for quasiincompressible flows. We determine corresponding sharp interface limits of two different scalings. The sharp interface limit is deduced by matched asymptotic expansions of the fields in powers of the interface. In particular, we study solutions of the derived system of inner equations and discuss the results within the general setting of jump conditions for sharp interface models. Furthermore, we treat, as a subproblem, the convective CahnHilliard equation numerically by a Local Discontinuous Galerkin scheme. 
W. Dreyer, J. Giesselmann, Ch. Kraus, Ch. Rohde, Asymptotic analysis for Korteweg models, Interfaces and Free Boundaries. Mathematical Modelling, Analysis and Computation, 14 (2012), pp. 105143.
Abstract
This paper deals with a sharp interface limit of the isothermal NavierStokesKorteweg system. The sharp interface limit is performed by matched asymptotic expansions of the fields in powers of the interface width. These expansions are considered in the interfacial region (inner expansions) and in the bulk (outer expansion) and are matched order by order. Particularly we consider the first orders of the corresponding inner equations obtained by a change of coordinates in an interfacial layer. For a specific scaling we establish solvability criteria for these inner equations and recover the results within the general setting of jump conditions for sharp interface models. 
K. Hermsdörfer, Ch. Kraus, D. Kröner, Interface conditions for limits of the NavierStokesKorteweg model, Interfaces and Free Boundaries. Mathematical Modelling, Analysis and Computation, 13 (2011), pp. 239254.
Abstract
In this contribution we will study the behaviour of the pressure across phase boundaries in liquidvapour flows. As mathematical model we will consider the static version of the NavierStokesKorteweg model which belongs to the class of diffuse interface models. From this static equation a formula for the pressure jump across the phase interface can be derived. If we perform then the sharp interface limit we see that the resulting interface condition for the pressure seems to be inconsistent with classical results of hydrodynamics. Therefore we will present two approaches to recover the results of hydrodynamics in the sharp interface limit at least for special situations. 
G. Kitavtsev, L. Recke, B. Wagner, Center manifold reduction approach for the lubrication equation, Nonlinearity, 24 (2011), pp. 23472369.
Abstract
The goal of this study is the reduction of the lubrication equation, modelling thin film dynamics, onto an approximate invariant manifold. The reduction is derived for the physical situation of the late phase evolution of a dewetting thin liquid film, where arrays of droplets connected by an ultrathin film of thickness $eps$ undergo a slowtime coarsening dynamics. With this situation in mind, we construct an asymptotic approximation of the corresponding invariant manifold, that is parametrized by a family of droplet pressures and positions, in the limit when $epsto 0$. The approach is inspired by the paper by Mielke and Zelik [Mem. Amer. Math. Soc., Vol. 198, 2009], where the center manifold reduction was carried out for a class of semilinear systems. In this study this approach is considered for quasilinear degenerate parabolic PDE's such as lubrication equations. While it has previously been shown by Glasner and Witelski [Phys. Rev. E, Vol. 67, 2003], that the system of ODEs governing the coarsening dynamics, can be obtained via formal asymptotic methods, the center manifold reduction approach presented here, pursues the rigorous justification of this asymptotic limit. 
A. Münch, B. Wagner, Impact of slippage on the morphology and stability of a dewetting rim, Journal of Physics: Condensed Matter, 23 (2011), pp. 184101/1184101/12.
Abstract
In this study lubrication theory is used to describe the stability and morphology of the rim that forms as a thin polymer film dewets from a hydrophobized silicon wafer. Thin film equations are derived from the governing hydrodynamic equations for the polymer to enable the systematic mathematical and numerical analysis of the properties of the solutions for different regimes of slippage and for a range of time scales. Dewetting rates and the cross sectional profiles of the evolving rims are derived for these models and compared to experimental results. Experiments also show that the rim is typically unstable in the spanwise direction and develops thicker and thinner parts that may grow into “fingers”. Linear stability analysis as well as nonlinear numerical solutions are presented to investigate shape and growth rate of the rim instability. It is demonstrated that the difference in morphology and the rate at which the instability develops can be directly attributed to the magnitude of slippage. Finally, a derivation is given for the dominant wavelength of the bulges along the unstable rim. 
A. Münch, C.P. Please, B. Wagner, Spin coating of an evaporating polymer solution, Physics of Fluids, 23 (2011), pp. 102101/1102101/12.
Abstract
We consider a mathematical model of spin coating of a single polymer blended in a solvent. The model describes the onedimensional development of the thin layer of the mixture as the layer thins due to flow created by a balance of viscous forces and centrifugal forces and due to evaporation of the solvent. In the model both the diffusivity of the solvent in the polymer and the viscosity of the mixture are very rapidly varying functions of the solvent volume fraction. Guided by numerical solutions an asymptotic analysis reveals a number of different possible behaviours of the thinning layer dependent on the nondimensional parameters describing the system. The main practical interest is in controlling the appearance and development of a “skin” on the polymer where the solvent concentration reduces rapidly on the outer surface leaving the bulk of the layer still with high concentrations of solvent. The critical parameters controlling this behaviour are found to be $eps$ the ratio of the diffusion to advection time scales, $delta$ the ratio of the evaporation to advection time scales and $exp(gamma)$, the ratio of the diffusivity of the initial mixture and the pure polymer. In particular, our analysis shows that for very small evaporation with $delta ll exp(3/(4gamma)) eps^3/4$ skin formation can be prevented. 
G. Kitavtsev, B. Wagner, Coarsening dynamics of slipping droplets, Journal of Engineering Mathematics, 66 (2010), pp. 271292.
Abstract
This paper studies the late phase dewetting process of nanoscopic thin polymer films on hydrophobized substrates using some recently derived lubrication models that take account of large slippage at the polymersubstrate interface. The late phase of this process is characterized by the slowtime coarsening dynamics of arrays of droplets that remain after rupture and the initial dewetting phases. For this situation a reduced system of ordinary differential equations is derived from the lubrication model for large slippage using asymptotic analysis. This extends known results for the noslip case. On the basis of the reduced model, the role of the slippage as a control parameter for droplet migration is analysed and several new qualitative effects for the coarsening process are identified. 
D. Peschka, A. Münch, B. Niethammer, Selfsimilar rupture of viscous thin films in the strongslip regime, Nonlinearity, 23 (2010), pp. 409427.

D. Peschka, A. Münch, B. Niethammer, Thin film rupture for large slip, Journal of Engineering Mathematics, 66 (2010), pp. 3351.
Abstract
This paper studies the rupture of thin liquid films on hydrophobic substrates, assuming large slip at the liquidsolid interface. Using a recently developed em strong slip lubrication model, it is shown that the rupture passes through up to three selfsimilar regimes with different dominant balances and different scaling exponents. For one of these regimes the similarity is of second kind, and the similarity exponent is determined by solving a boundary value problem for a nonlinear ODE. For this regime we also prove finitetime rupture. 
W. Dreyer, Ch. Kraus, On the van der WaalsCahnHilliard phasefield model and its equilibria conditions in the sharp interface limit, Proceedings of the Royal Society of Edinburgh. Section A. Mathematics, 140 A (2010), pp. 11611186.
Abstract
We study the equilibria of liquidvapor phase transitions of a single substance at constant temperature and relate the sharp interface model of classical thermodynamics to a phase field model that determines the equilibria by the stationary van der WaalsCahnHilliard theory.
For two reasons we reconsider this old problem. 1. Equilibria in a two phase system can be established either under fixed total volume of the system or under fixed external pressure. The latter case implies that the domain of the twophase system varies. However, in the mathematical literature rigorous sharp interface limits of phase transitions are usually considered under fixed volume. This brings the necessity to extend the existing tools for rigorous sharp interface limits to changing domains since in nature most processes involving phase transitions run at constant pressure. 2. Thermodynamics provides for a single substance two jump conditions at the sharp interface, viz. the continuity of the specific Gibbs free energies of the adjacent phases and the discontinuity of the corresponding pressures, which is balanced by the mean curvature. The existing estimates for rigorous sharp interface limits show only the first condition. We identify the cause of this phenomenon and develop a strategy that yields both conditions up to the first order.
The necessary information on the equilibrium conditions are achieved by an asymptotic expansion of the density which is valid for an arbitrary double well potential. We establish this expansion by means of local energy estimates, uniform convergence results of the density and estimates on the Laplacian of the density. 
D. Peschka, A. Münch, B. Niethammer, Selfsimilar rupture of viscous thin films in the strongslip regime, Nonlinearity, 23 (2010), pp. 409427.
Abstract
We consider rupture of thin viscous films in the strongslip regime with small Reynolds numbers. Numerical simulations indicate that near the rupture point viscosity and vanderWaals forces are dominant and that there are selfsimilar solutions of the second kind. For a corresponding simplified model we rigorously analyse selfsimilar behaviour. There exists a oneparameter family of selfsimilar solutions and we establish necessary and sufficient conditions for convergence to any selfsimilar solution in a certain parameter regime. We also present a conjecture on the domains of attraction of all selfsimilar solutions which is supported by numerical simulations. 
J.R. King, A. Münch, B. Wagner, Linear stability analysis of a sharpinterface model for dewetting thin films, Journal of Engineering Mathematics, 63 (2009), pp. 177195.
Abstract
The topic of this study concerns the stability of the threephase contactline of a dewetting thin liquid film on a hydrophobised substrate driven by van der Waals forces. The role of slippage in the emerging instability at the threephase contactline is studied by deriving a sharpinterface model for the dewetting thin film via matched asymptotic expansions. This allows for a derivation of travelling waves and their linear stability via eigenmode analysis. In contrast to the dispersion relations typically encountered for the fingerinstabilty, where the dependence of the growth rate on the wave number is quadratic, here it is linear. Using the separation of time scales of the slowly growing rim of the dewetting film and time scale on which the contact line destabilises, the sharpinterface results are compared to earlier results for the full lubrication model and good agreement for the most unstable modes is obtained. 
K. Afanasiev, A. Münch, B. Wagner, Thin film dynamics on a vertically rotating disk partially immersed in a liquid bath, Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems. Elsevier Science Inc., New York, NY. English, English abstracts., 32 (2008), pp. 18941911.
Abstract
The axisymmetric flow of a thin liquid film subject to surface tension, gravity and centrifugal forces is considered for the problem of a vertically rotating disk that is partially immersed in a liquid bath. This problem constitutes a generalization of the classic LandauLevich dragout problem to axisymmetric flow. A generalized lubrication model that includes the meniscus region connecting the thin film to the bath is derived. The resulting nonlinear fourthorder partial differential equation is solved numerically using a finite element scheme. For a range of parameters steady states are found. While the solutions for the height profile of the film near the dragout region show excellent agreement with the asymptotic solutions to the corresponding classic LandauLevich problem, they show novel patterns away from the meniscus region. The implications for possible industrial applications are discussed. 
A. Münch, B. Wagner, Galerkin method for feedback controlled RayleighBénard convection, Nonlinearity, 21 (2008), pp. 26252651.

M. Rauscher, R. Blossey, A. Münch, B. Wagner, Spinodal dewetting of thin films with large interfacial slip: Implications from the dispersion relation, Langmuir, 24 (2008), pp. 1229012294.

K. Afanasiev, A. Münch, B. Wagner, On the LandauLevich problem for nonNewtonian liquids, Physical Review E. Statistical, Nonlinear, and Soft Matter Physics, 76 (2007), pp. 036307/1036307/12.
Abstract
In this paper the dragout problem for shearthinning liquids at variable inclination angle is considered. For this free boundary problem dimensionreduced lubrication equations are derived for the most commonly used viscosity models, namely, the powerlaw, Ellis and Carreau model. For the resulting lubrication models a system of ordinary differential equation governing the steady state solutions is obtained. Phase plane analysis is used to characterize the type of possible steady state solutions and their dependence on the rheological parameters. 
P. Krejčí, J. Sprekels, Elasticideally plastic beams and PrandtlIshlinskii hysteresis operators, Mathematical Methods in the Applied Sciences, 30 (2007), pp. 23712393.

P. Evans, A. Münch, Interaction of advancing fronts and meniscus profiles formed by surfacetensionongradientdriven liquid films, SIAM Journal on Applied Mathematics, 66 (2006), pp. 16101631.

P.L. Evans, J.R. King, A. Münch, Intermediateasymptotic structure of a dewetting rim with strong slip, Applied Mathematics Research Express, (2006), pp. 25262/125262/25.
Abstract
When a thin viscous liquid film dewets, it typically forms a rim which spreads outwards, leaving behind a growing dry region. We consider the dewetting behaviour of a film, when there is strong slip at a liquidsubstrate interface. The film can be modelled by two coupled partial differential equations (PDEs) describing the film thickness and velocity. Using asymptotic methods, we describe the structure of the rim as it evolves in time, and the rate of dewetting, in the limit of large slip lengths. An inner region emerges, closest to the dewetted region, where surface tension is important; in an outer region, three subregions develop. This asymptotic description is compared with numerical solutions of the full system of PDEs. 
J.R. King, A. Münch, B. Wagner, Linear stability of a ridge, Nonlinearity, 19 (2006), pp. 28132831.
Abstract
We investigate the stability of the threephase contactline of a thin liquid ridge on a hydrophobic substrate for flow driven by surface tension and van der Waals forces. We study the role of slippage in the emerging instability at the threephase contactline by comparing the lubrication models for noslip and slipdominated conditions at the liquid/substrate interface. For both cases we derive a sharpinterface model via matched asymptotic expansions and derive the eigenvalues from a linear stability analysis of the respective reduced models. We compare our asymptotic results with the eigenvalues obtained numerically for the full lubrication models. 
A. Münch, B. Wagner, M. Rauscher, R. Blossey, A thin film model for corotational Jeffreys fluids under strong slip, The European Physical Journal. E. Soft Matter, 20 (2006), pp. 365368.
Abstract
We derive a thin film model for viscoelastic liquids under strong slip which obey the stress tensor dynamics of corotational Jeffreys fluids. 
B. Jin, A. Acrivos, A. Münch, The dragout problem in film coating, Physics of Fluids, 17 (2005), pp. 103603/1103603/12.

A. Münch, P.L. Evans, Marangonidriven liquid films rising out of a meniscus onto a nearlyhorizontal substrate, Physica D. Nonlinear Phenomena, 209 (2005), pp. 164177.

A. Münch, B. Wagner, Th.P. Witelski, Lubrication models for small to large slip lengths, Journal of Engineering Mathematics, 53 (2005), pp. 359383.

A. Münch, B. Wagner, Contactline instability for dewetting thin films, Physica D. Nonlinear Phenomena, 209 (2005), pp. 178190.
Preprints, Reports, Technical Reports

A. Münch, B. Wagner, Selfconsistent field theory for a polymer brush. Part II: The effective chemical potential, Preprint no. 2649, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2649 .
Abstract, PDF (318 kByte)
The most successful meanfield model to describe the collective behaviour of the large class of macromolecular polymers is the selfconsistent field theory (SCFT). Still, even for the simple system of a grafted dry polymer brush, the meanfield equations have to be solved numerically. As one of very few alternatives that offer some analytical tractability the strongstretching theory (SST) has led to explicit expressions for the effective chemical potential and consequently the free energy to promote an understanding of the underlying physics. Yet, a direct derivation of these analytical results from the SCFT model is still outstanding. In this study we present a systematic asymptotic theory based on matched asymtptotic expansions to obtain the effective chemical potential from the SCFT model for a dry polymer brush for large but finite stretching. 
A. Münch, B. Wagner, Selfconsistent field theory for a polymer brush. Part I: Asymptotic analysis in the strongstretching limit, Preprint no. 2648, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2648 .
Abstract, PDF (854 kByte)
In this study we consider the selfconsistent field theory for a dry, in compressible polymer brush, densely grafted on a substrate, describing the average segment density of a polymer in terms of an effective chemical potential for the interaction between the segments of the polymer chain. We present a systematic singular perturbation analysis of the selfconsistent field theory in the strongstretching limit, when the length scale of the ratio of the radius of gyration of the polymer chain to the extension of the brush from the substrate vanishes. Our analysis yields, for the first time, an approximation for the average segment density that is correct to leading order in the outer scaling and resolves the boundary layer singularity at the end of the polymer brush in the strongstretching limit. We also show that in this limit our analytical results agree increasingly well with our numerical solutions to the full model equations comprising the selfconsistent field theory. 
P. Colli, G. Gilardi, J. Sprekels, Asymptotic analysis of a tumor growth model with fractional operators, Preprint no. 2625, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2625 .
Abstract, PDF (341 kByte)
In this paper, we study a system of three evolutionary operator equations involving fractional powers of selfadjoint, monotone, unbounded, linear operators having compact resolvents. This system constitutes a generalized and relaxed version of a phase field system of CahnHilliard type modelling tumor growth that has originally been proposed in HawkinsDaarud et al. (Int. J. Numer. Math. Biomed. Eng. 28 (2012), 324). The original phase field system and certain relaxed versions thereof have been studied in recent papers coauthored by the present authors and E. Rocca. The model consists of a CahnHilliard equation for the tumor cell fraction φ, coupled to a reactiondiffusion equation for a function S representing the nutrientrich extracellular water volume fraction. Effects due to fluid motion are neglected. Motivated by the possibility that the diffusional regimes governing the evolution of the different constituents of the model may be of different (e.g., fractional) type, the present authors studied in a recent note a generalization of the systems investigated in the abovementioned works. Under rather general assumptions, wellposedness and regularity results have been shown. In particular, by writing the equation governing the evolution of the chemical potential in the form of a general variational inequality, also singular or nonsmooth contributions of logarithmic or of double obstacle type to the energy density could be admitted. In this note, we perform an asymptotic analysis of the governing system as two (small) relaxation parameters approach zero separately and simultaneously. Corresponding wellposedness and regularity results are established for the respective cases; in particular, we give a detailed discussion which assumptions on the admissible nonlinearities have to be postulated in each of the occurring cases. 
M. Heida, S. Neukamm, M. Varga, Stochastic homogenization of Lambdaconvex gradient flows, Preprint no. 2594, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2594 .
Abstract, PDF (429 kByte)
In this paper we present a stochastic homogenization result for a class of Hilbert space evolutionary gradient systems driven by a quadratic dissipation potential and a Λconvex energy functional featuring random and rapidly oscillating coefficients. Specific examples included in the result are AllenCahn type equations and evolutionary equations driven by the pLaplace operator with p ∈ in (1, ∞). The homogenization procedure we apply is based on a stochastic twoscale convergence approach. In particular, we define a stochastic unfolding operator which can be considered as a random counterpart of the wellestablished notion of periodic unfolding. The stochastic unfolding procedure grants a very convenient method for homogenization problems defined in terms of (Λ)convex functionals. 
A. Caiazzo, R. Maier, D. Peterseim, Reconstruction of quasilocal numerical effective models from lowresolution measurements, Preprint no. 2577, WIAS, Berlin, 2019, DOI 10.20347/WIAS.PREPRINT.2577 .
Abstract, PDF (1550 kByte)
We consider the inverse problem of reconstructing an effective model for a prototypical diffusion process in strongly heterogeneous media based on lowresolution measurements. We rely on recent quasilocal numerical effective models that, in contrast to conventional homogenized models, are provably reliable beyond periodicity assumptions and scale separation. The goal of this work is to show that the identification of the matrix representation of these effective models is possible. Algorithmic aspects of the inversion procedure and its performance are illustrated in a series of numerical experiments.
Talks, Poster

M. Heida, Stochastic homogenization of PDE on nonuniformly Lipschitz and percolating structures, DMVJahrestagung 2019, September 23  26, 2019, KIT  Karlsruher Institut für Technologie, September 24, 2019.

M. Heida, The SQRA operator: convergence behaviour and applications, Universität Wien, Austria, March 19, 2019.

B. Wagner, Free boundary problems of active and driven hydrogels, PIMSGermany Workshop on Modelling, Analysis and Numerical Analysis of PDEs for Applications, June 24  26, 2019, Universität Heidelberg, Interdisciplinary Center for Scientifiv Computing (IWR) and BIOQUANT Center, June 24, 2019.

B. Wagner, Illposedness of twophase flow models of concentrated suspensions, International Congress on Industrial and Applied Mathematics ICIAM2019, Minisymposium MS ME07 6 Recent advances in understanding suspensions and granular media flow  Part 2 of 2, July 15  19, 2019, Valencia, Spain, July 17, 2019.

B. Wagner, Mathematical modeling of real world processes, CERN Academic Training Programme 20182019, March 14  15, 2019, CERN, Genf, Switzerland.

M. Heida, The SQRA operator: convergence behaviour and applications, Polytechnico Milano, Italy, March 13, 2019.

A. Caiazzo, Homogenization methods for weakly compressible elastic materials forward and inverse problem, Workshop on Numerical Inverse and Stochastic Homogenization, February 13  17, 2017, Universität Bonn, Hausdorff Research Institute for Mathematics, February 17, 2017.

K. Disser, Econvergence to the quasisteadystate approximation in systems of chemical reactions, ERC Workshop on Modeling Materials and Fluids using Variational Methods, February 22  26, 2016, WIAS Berlin, Berlin, February 25, 2016.

S. Reichelt, Competing patterns in antisymmetrically coupled SwiftHohenberg equations, International Conference on Control of Complex Systems and Networks, SFB 910 ``Control of SelfOrganizing Nonlinear Systems: Theoretical Methods and Concepts of Application'', Heringsdorf/Usedom, September 4  8, 2016.

M. Heida, Homogenization of the random conductance model, 7th European Congress of Mathematics (ECM), session ``Probability, Statistics and Financial Mathematics'', July 18  22, 2016, Technische Universität Berlin, Berlin, July 20, 2016.

M. Heida, Homogenization of the random conductance model, Workshop of the GAMM Activity Group ``Analysis of Partial Differential Equations'', September 26  28, 2016, Technische Universität Dortmund, Fachbereich Mathematik, Dortmund, September 26, 2016.

M. Heida, Stochastic homogenization of 1homogeneous functionals, 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.

M. Heida, Stochastic homogenization of rateindependent systems, Berlin Dresden Prague Würzburg Workshop ``Homogenization and Related Topics'', Technische Universität Dresden, Fachbereich Mathematik, June 22, 2016.

M. Heida, Stochastic homogenization of rateindependent systems, Joint Annual Meeting of DMV and GAMM, Session ``Multiscales and Homogenization'', March 7  11, 2016, Technische Universität Braunschweig, Braunschweig, March 10, 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.

J. Fuhrmann, Ch. Merdon, A thermodynamically consistent numerical approach to NernstPlanckPoisson systems with volume constraints, The 67th Annual Meeting of the International Society of Electrochemistry, Den Haag, Netherlands, August 21  26, 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.

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

A. Caiazzo, Multiscale modeling of weakly compressible elastic materials in harmonic regime, Université de FrancheComté, Laboratoire de Mathématiques de Besançon, France, July 1, 2014.

S. Neukamm, Optimal decay estimate on the semigroup associated with a random walk among random conductances, Dirichlet Forms and Applications, GermanJapanese Meeting on Stochastic Analysis, September 9  13, 2013, Universität Leipzig, Mathematisches Institut, September 9, 2013.

H. Abels, J. Daube, Ch. Kraus, D. Kröner, Sharp interface limit for the NavierStokesKorteweg model, DIMO2013  Diffuse Interface Models, Levico Terme, Italy, September 10  13, 2013.

CH. Kraus, Damage and phase separation processes: Modeling and analysis of nonlinear PDE systems, DIMO2013  Diffuse Interface Models, September 10  13, 2013, Levico Terme, Italy, September 11, 2013.

CH. Kraus, Modeling and analysis of a nonlinear PDE system for phase separation and damage, Università di Pavia, Dipartimento di Matematica, Italy, January 22, 2013.

CH. Kraus, Sharp interface limit of a diffuse interface model of NavierStokesAllenCahn type for mixtures, Workshop ``Hyperbolic Techniques for Phase Dynamics'', June 10  14, 2013, Mathematisches Forschungsinstitut Oberwolfach, June 11, 2013.

CH. Kraus, A nonlinear PDE system for phase separation and damage, Universität Freiburg, Abteilung Angewandte Mathematik, November 13, 2012.

CH. Kraus, CahnLarché systems coupled with damage, Università degli Studi di Milano, Dipartimento di Matematica, Italy, November 28, 2012.

CH. Kraus, Phase field systems for phase separation and damage processes, 12th International Conference on Free Boundary Problems: Theory and Applications, June 11  15, 2012, Frauenchiemsee, June 12, 2012.

CH. Kraus, Phasenfeldsysteme für Entmischungs und Schädigungsprozesse, Mathematisches Kolloquium, Universität Stuttgart, Fachbereich Mathematik, May 15, 2012.

CH. Kraus, Diffuse interface systems for phase separation and damage, Seminar on Partial Differential Equations, Institute of Mathematics, Academy of Sciences of the Czech Republic, Prague, May 3, 2011.

CH. Kraus, Phase separation systems coupled with elasticity and damage, ICIAM 2011, July 18  22, 2011, Vancouver, Canada, July 18, 2011.

W. Dreyer, On a paradox within the phase field modeling of storage systems and its resolution, 8th AIMS International Conference on Dynamical Systems, Differential Equations and Applications, May 25  28, 2010, Technische Universität Dresden, May 26, 2010.

CH. Kraus, An inhomogeneous, anisotropic and elastically modified GibbsThomson law as singular limit of a diffuse interface model, 81st Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM), March 22  26, 2010, Karlsruhe, March 23, 2010.

CH. Kraus, Inhomogeneous and anisotropic phasefield quantities in the sharp interface limit, 6th Singular Days 2010, April 29  May 1, 2010, WIAS, Berlin, April 30, 2010.

G. Kitavtsev, L. Recke, B. Wagner, Derivation, analysis and numerics of reduced ODE models describing coarsening dynamics, 3textsuperscriptrd European Postgraduate Fluid Dynamics Conference, Nottingham, UK, July 13  16, 2009.

G. Kitavtsev, Derivation, analysis and numerics of reduced ODE models describing coarsening dynamics, 3$^rm rd$ European Postgraduate Fluid Dynamics Conference, July 13  16, 2009, University of Nottingham, UK, July 15, 2009.

G. Kitavtsev, Reduced ODE models describing coarsening dynamics of slipping droplets and a geometrical approach for their derivation, Oberseminar, Universität Bonn, Institut für Angewandte Mathematik, July 23, 2009.

W. Dreyer, On a paradox within the phase field modeling of storage systems and its resolution, PF09  2nd Symposium on PhaseField Modelling in Materials Science, August 30  September 2, 2009, Universität Aachen, Kerkrade, Netherlands, August 31, 2009.

W. Dreyer, Phase field models and their sharp limits in the context of hydrogen storage and lithiumion batteries, 1textsuperscriptst International Conference on Material Modelling (1textsuperscriptst ICMM), September 15  18, 2009, Technische Universität Dortmund, September 16, 2009.

CH. Kraus, A phasefield model with anisotropic surface tension in the sharp interface limit, Second GAMMSeminar on Multiscale Material Modelling, July 10  12, 2008, Universität Stuttgart, Institut für Mechanik (Bauwesen), July 12, 2008.

CH. Kraus, Ein Phasenfeldmodell vom CahnHilliardTyp im singulären Grenzwert, Oberseminar Analysis, Universität Regensburg, Fakultät für Mathematik, April 25, 2008.

CH. Kraus, Phase field models and corresponding GibbsThomson laws. Part II, SIMTECH Seminar Multiscale Modelling in Fluid Mechanics, Universität Stuttgart, Institut für Angewandte Analysis und Numerische Simulation, November 5, 2008.

CH. Kraus, On jump conditions at phase interfaces, Oberseminar über Angewandte Mathematik, December 10  15, 2007, Universität Freiburg, Abteilung für Angewandte Mathematik, December 11, 2007.

CH. Kraus, Equilibrium conditions for liquidvapor system in the sharp interface limit, Seminar Thermodynamische Modellierung und Analyse von Phasenübergängen, WIAS, Berlin, July 18, 2006.

CH. Kraus, Equilibria conditions in the sharp interface limit of the van der WaalsCahnHilliard phase model, Recent Advances in Free Boundary Problems and Related Topics (FBP2006), September 14  16, 2006, Levico, Italy, September 14, 2006.

CH. Kraus, The sharp interface limit of the van der WaalsCahnHilliard model, PolishGerman Workshop ``Modeling Structure Formation'', Interdisciplinary Centre for Mathematical and Computational Modelling, University of Warsaw, Poland, September 8, 2006.

CH. Kraus, On the sharp limit of the Van der WaalsCahnHilliard model, WIAS Workshop ``Dynamic of Phase Transitions'', November 30  December 3, 2005, Berlin, December 2, 2005.

CH. Kraus, On the sharp limit of the Van der WaalsCahnHilliard model, Workshop ``MicroMacro Modeling and Simulation of LiquidVapor Flows'', November 16  18, 2005, Universität Freiburg, Mathematisches Institut, Kirchzarten, November 17, 2005.

CH. Kraus, Maximale Konvergenz in höheren Dimensionen, Seminar Thermodynamische Modellierung und Analyse von Phasenübergängen, WIAS, Berlin, May 24, 2005.
Contributing Groups of WIAS
Projects/Grants
 Datadriven electronic structure calculation for nanostructures
 Effective models for materials and interfaces with multiple scales
 Electrochemical Double Layers in Solid Oxide Cells (EDLSOC)
 Fault networks and scaling properties of deformation accumulation
 Hybrid models for the electrothermal behavior of organic semiconductor devices
 Interface dynamics: Bridging stochastic and hydrodynamic descriptions
 Modellbasierte Abschätzung der Lebensdauer von gealterten LiBatterien für die 2ndLife Anwendung als stationärer Stromspeicher (MALLi2)
 Reliability of efficient approximation schemes for material discontinuities described by functions of bounded variation
 Stochastic spatial coagulation particle processes
 Verbundprojekt “Wege zu sekundären Mg/Ca  LuftBatterien”