2011, Dreyer, Wolfgang, Hantke, Maren, Warnecke, Gerald

We consider the isothermal Euler equations with phase transition between a liquid and a vapor phase. The mass transfer is modeled by a kinetic relation. We prove existence and uniqueness results. Further, we construct the exact solution for Riemann problems. We derive analogous results for the cases of initially one phase with resulting condensation by compression or evaporation by expansion. Further we present numerical results for these cases. We compare the results to similar problems without phase transition.

2014, Dziwnik, Marion, Münch, Andreas, Wagner, Barbara

We propose a phase field model for solid state dewetting in form of a Cahn-Hilliard equation with weakly anisotropic surface energy and a degenerate mobility together with a free boundary condition at the film-substrate 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 well-known Young's equation

2019, Hintermüller, Michael, Keil, Tobias

Optimal control problems for geometric (evolutionary) partial differential inclusions are considered. The focus is on problems which, in addition to the nonlinearity due to geometric evolution, contain optimization theoretic challenges because of non-smoothness. The latter might stem from energies containing non-smooth constituents such as obstacle-type potentials or terms modeling, e.g., pinning phenomena in microfluidics. Several techniques to remedy the resulting constraint degeneracy when deriving stationarity conditions are presented. A particular focus is on Yosida-type mollifications approximating the original degenerate problem by a sequence of nondegenerate nonconvex optimal control problems. This technique is also the starting point for the development of numerical solution schemes. In this context, also dual-weighted residual based error estimates are addressed to facilitate an adaptive mesh refinement. Concerning the underlying state model, sharp and diffuse interface formulations are discussed. While the former always allows for accurately tracing interfacial motion, the latter model may be dictated by the underlying physical phenomenon, where near the interface mixed phases may exist, but it may also be used as an approximate model for (sharp) interface motion. In view of the latter, (sharp interface) limits of diffuse interface models are addressed. For the sake of presentation, this exposition confines itself to phase field type diffuse interface models and, moreover, develops the optimal control of either of the two interface models along model applications. More precisely, electro-wetting on dielectric is used in the sharp interface context, and the control of multiphase fluids involving spinodal decomposition highlights the phase field technique. Mathematically, the former leads to a Hele-Shaw flow with geometric boundary conditions involving a complementarity system due to contact line pinning, and the latter gives rise to a Cahn-Hilliard Navier-Stokes model including a non-smooth obstacle type potential leading to a variational inequality constraint.

2014, Korzec, Maciek, Münch, Andreas, Süli, Endre, Wagner, Barbara

Anisotropy is an essential feature of phase-field 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 phase-field 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