-
Well-posedness of an evaporation model for a spherical droplet exposed to an air flow
Authors:
Eberhard Bänsch,
Martin Doß,
Carsten Gräser,
Nadja Ray
Abstract:
In this paper, we address the well-posedness of an evaporation model for a spherical liquid droplet taking into account the convective impact of an air flow in the ambient gas phase. From a mathematical perspective, we are dealing with a coupled ODE-PDE system for the droplet radius, the temperature distribution, and the vapor concentration. The nonlinear coupling arises from the evaporation rate…
▽ More
In this paper, we address the well-posedness of an evaporation model for a spherical liquid droplet taking into account the convective impact of an air flow in the ambient gas phase. From a mathematical perspective, we are dealing with a coupled ODE-PDE system for the droplet radius, the temperature distribution, and the vapor concentration. The nonlinear coupling arises from the evaporation rate modeled by the Hertz-Knudsen equation. Under physically meaningful assumptions, we prove existence and uniqueness of a weak solution until the droplet has evaporated completely. Numerical simulations are performed to illustrate how different air flows affect the evaporation process.
△ Less
Submitted 11 December, 2023;
originally announced December 2023.
-
Numerical simulation of multiscale fault systems with rate- and state-dependent friction
Authors:
Carsten Gräser,
Ralf Kornhuber,
Joscha Podlesny
Abstract:
We consider the deformation of a geological structure with non-intersecting faults that can be represented by a layered system of viscoelastic bodies satisfying rate- and state-depending friction conditions along the common interfaces. We derive a mathematical model that contains classical Dieterich- and Ruina-type friction as special cases and accounts for possibly large tangential displacements.…
▽ More
We consider the deformation of a geological structure with non-intersecting faults that can be represented by a layered system of viscoelastic bodies satisfying rate- and state-depending friction conditions along the common interfaces. We derive a mathematical model that contains classical Dieterich- and Ruina-type friction as special cases and accounts for possibly large tangential displacements. Semi-discretization in time by a Newmark scheme leads to a coupled system of non-smooth, convex minimization problems for rate and state to be solved in each time step. Additional spatial discretization by a mortar method and piecewise constant finite elements allows for the decoupling of rate and state by a fixed point iteration and efficient algebraic solution of the rate problem by truncated non-smooth Newton methods. Numerical experiments with a spring slider and a layered multiscale system illustrate the behavior of our model as well as the efficiency and reliability of the numerical solver.
△ Less
Submitted 14 January, 2022; v1 submitted 27 October, 2021;
originally announced October 2021.
-
An evolving space framework for Oseen equations on a moving domain
Authors:
Ana Djurdjevac,
Carsten Gräser,
Philip J. Herbert
Abstract:
This article considers non-stationary incompressible linear fluid equations in a moving domain. We demonstrate the existence and uniqueness of an appropriate weak formulation of the problem by making use of the theory of time-dependent Bochner spaces. It is not possible to directly apply established evolving Hilbert space theory due to the incompressibility constraint. After we have established th…
▽ More
This article considers non-stationary incompressible linear fluid equations in a moving domain. We demonstrate the existence and uniqueness of an appropriate weak formulation of the problem by making use of the theory of time-dependent Bochner spaces. It is not possible to directly apply established evolving Hilbert space theory due to the incompressibility constraint. After we have established the well-posedness, we derive and analyse a first order time discretisation of the system.
△ Less
Submitted 8 September, 2023; v1 submitted 15 June, 2021;
originally announced June 2021.
-
Free energy computation of particles with membrane-mediated interactions via Langevin dynamics
Authors:
Tobias Kies,
Carsten Gräser,
Luigi Delle Site,
Ralf Kornhuber
Abstract:
We apply well-established concepts of Langevin sampling to derive a new class of algorithms for the efficient computation of free energy differences of fluctuating particles embedded in a 'fast' membrane, i.e., a membrane that instantaneously adapts to varying particle positions. A geometric potential accounting for membrane-mediated particle interaction is derived in the framework of variational…
▽ More
We apply well-established concepts of Langevin sampling to derive a new class of algorithms for the efficient computation of free energy differences of fluctuating particles embedded in a 'fast' membrane, i.e., a membrane that instantaneously adapts to varying particle positions. A geometric potential accounting for membrane-mediated particle interaction is derived in the framework of variational hybrid models for particles in membranes. Recent explicit representations of the gradient of the geometric interaction potential allows to apply well-known gradient based Markov Chain Monte-Carlo (MCDC) methods such as Langevin-based sampling.
△ Less
Submitted 30 September, 2020;
originally announced September 2020.
-
Truncated Nonsmooth Newton Multigrid for phase-field brittle-fracture problems, with analysis
Authors:
Carsten Gräser,
Daniel Kienle,
Oliver Sander
Abstract:
We propose the Truncated Nonsmooth Newton Multigrid Method (TNNMG) as a solver for the spatial problems of the small-strain brittle-fracture phase-field equations. TNNMG is a nonsmooth multigrid method that can solve biconvex, block-separably nonsmooth minimization problems with linear time complexity. It exploits the variational structure inherent in the problem, and handles the pointwise irrever…
▽ More
We propose the Truncated Nonsmooth Newton Multigrid Method (TNNMG) as a solver for the spatial problems of the small-strain brittle-fracture phase-field equations. TNNMG is a nonsmooth multigrid method that can solve biconvex, block-separably nonsmooth minimization problems with linear time complexity. It exploits the variational structure inherent in the problem, and handles the pointwise irreversibility constraint on the damage variable directly, without regularization or the introduction of a local history field. In the paper we introduce the method and show how it can be applied to several established models of phase-field brittle fracture. We then prove convergence of the solver to a solution of the nonsmooth Euler--Lagrange equations of the spatial problem for any load and initial iterate. On the way, we show several crucial convexity and regularity properties of the models considered here. Numerical comparisons to an operator-splitting algorithm show a considerable speed increase, without loss of robustness.
△ Less
Submitted 10 March, 2023; v1 submitted 23 July, 2020;
originally announced July 2020.
-
Error bounds for PDE-regularized learning
Authors:
Carsten Gräser,
Prem Anand Alathur Srinivasan
Abstract:
In this work we consider the regularization of a supervised learning problem by partial differential equations (PDEs) and derive error bounds for the obtained approximation in terms of a PDE error term and a data error term. Assuming that the target function satisfies an unknown PDE, the PDE error term quantifies how well this PDE is approximated by the auxiliary PDE used for regularization. It is…
▽ More
In this work we consider the regularization of a supervised learning problem by partial differential equations (PDEs) and derive error bounds for the obtained approximation in terms of a PDE error term and a data error term. Assuming that the target function satisfies an unknown PDE, the PDE error term quantifies how well this PDE is approximated by the auxiliary PDE used for regularization. It is shown that this error term decreases if more data is provided. The data error term quantifies the accuracy of the given data. Furthermore, the PDE-regularized learning problem is discretized by generalized Galerkin discretizations solving the associated minimization problem in subsets of the infinite dimensional functions space, which are not necessarily subspaces. For such discretizations an error bound in terms of the PDE error, the data error, and a best approximation error is derived.
△ Less
Submitted 13 March, 2020;
originally announced March 2020.
-
The DUNE Framework: Basic Concepts and Recent Developments
Authors:
Peter Bastian,
Markus Blatt,
Andreas Dedner,
Nils-Arne Dreier,
Christian Engwer,
René Fritze,
Carsten Gräser,
Christoph Grüninger,
Dominic Kempf,
Robert Klöfkorn,
Mario Ohlberger,
Oliver Sander
Abstract:
This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state [1, 2]. This discussion is accompanied with a description of various advanced features, such as coupling of domains and…
▽ More
This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state [1, 2]. This discussion is accompanied with a description of various advanced features, such as coupling of domains and cut cells, grid modifications such as adaptation and moving domains, high order discretizations and node level performance, non-smooth multigrid methods, and multiscale methods. A brief discussion on current and future development directions of the framework concludes the paper.
△ Less
Submitted 22 June, 2020; v1 submitted 30 September, 2019;
originally announced September 2019.
-
Function space bases in the dune-functions module
Authors:
Christian Engwer,
Carsten Gräser,
Steffen Müthing,
Oliver Sander
Abstract:
The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for p…
▽ More
The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for product spaces as trees of simpler bases. From this description, many different numberings of degrees of freedom by multi-indices can be derived in a natural way. We describe the abstract concepts, document the programmer interface, and give a complete example program that solves the stationary Stokes equation using Taylor-Hood elements.
△ Less
Submitted 25 June, 2018;
originally announced June 2018.
-
CRC 1114 - Report Membrane Deformation by N-BAR Proteins: Extraction of membrane geometry and protein diffusion characteristics from MD simulations
Authors:
Jan Henning Peters,
Carsten Gräser,
Rupert Klein
Abstract:
We describe simulations of Proteins and artificial pseudo-molecules interacting and shaping lipid bilayer membranes. We extract protein diffusion Parameters, membrane deformation profiles and the elastic properties of the used membrane models in preparation of calculations based on a large scale continuum model.
We describe simulations of Proteins and artificial pseudo-molecules interacting and shaping lipid bilayer membranes. We extract protein diffusion Parameters, membrane deformation profiles and the elastic properties of the used membrane models in preparation of calculations based on a large scale continuum model.
△ Less
Submitted 7 December, 2017;
originally announced December 2017.
-
On differentiability of the membrane-mediated mechanical interaction energy of discrete-continuum membrane-particle models
Authors:
Tobias Kies,
Carsten Gräser
Abstract:
We consider a discrete-continuum model of a biomembrane with embedded particles. While the membrane is represented by a continuous surface, embedded particles are described by rigid discrete objects which are free to move and rotate in lateral direction. For the membrane we consider a linearized Canham-Helfrich energy functional and height and slope boundary conditions imposed on the particle boun…
▽ More
We consider a discrete-continuum model of a biomembrane with embedded particles. While the membrane is represented by a continuous surface, embedded particles are described by rigid discrete objects which are free to move and rotate in lateral direction. For the membrane we consider a linearized Canham-Helfrich energy functional and height and slope boundary conditions imposed on the particle boundaries resulting in a coupled minimization problem for the membrane shape and particle positions.
When considering the energetically optimal membrane shape for each particle position we obtain a reduced energy functional that models the implicitly given interaction potential for the membrane-mediated mechanical particle-particle interactions. We show that this interaction potential is differentiable with respect to the particle positions and orientations. Furthermore we derive a fully practical representation of the derivative only in terms of well defined derivatives of the membrane. This opens the door for the application of minimization algorithms for the computation of minimizers of the coupled system and for further investigation of the interaction potential of membrane-mediated mechanical particle--particle interaction.
The results are illustrated with numerical examples comparing the explicit derivative formula with difference quotient approximations. We furthermore demonstrate the application of the derived formula to implement a gradient flow for the approximation of optimal particle configurations.
△ Less
Submitted 28 April, 2021; v1 submitted 29 November, 2017;
originally announced November 2017.
-
Truncated Nonsmooth Newton Multigrid Methods for Block-Separable Minimization Problems
Authors:
Carsten Gräser,
Oliver Sander
Abstract:
The Truncated Nonsmooth Newton Multigrid (TNNMG) method is a robust and efficient solution method for a wide range of block-separable convex minimization problems, typically stemming from discretizations of nonlinear and nonsmooth partial differential equations. This paper proves global convergence of the method under weak conditions both on the objective functional, and on the local inexact subpr…
▽ More
The Truncated Nonsmooth Newton Multigrid (TNNMG) method is a robust and efficient solution method for a wide range of block-separable convex minimization problems, typically stemming from discretizations of nonlinear and nonsmooth partial differential equations. This paper proves global convergence of the method under weak conditions both on the objective functional, and on the local inexact subproblem solvers that are part of the method. It also discusses a range of algorithmic choices that allows to customize the algorithm for many specific problems. Numerical examples are deliberately omitted, because many such examples have already been published elsewhere.
△ Less
Submitted 25 October, 2017; v1 submitted 14 September, 2017;
originally announced September 2017.
-
Discretization error estimates for penalty formulations of a linearized Canham-Helfrich type energy
Authors:
Carsten Gräser,
Tobias Kies
Abstract:
This paper is concerned with minimization of a fourth-order linearized Canham-Helfrich energy subject to Dirichlet boundary conditions on curves inside the domain. Such problems arise in the modeling of the mechanical interaction of biomembranes with embedded particles. There, the curve conditions result from the imposed particle--membrane coupling. We prove almost-$H^{\frac{5}{2}}$ regularity of…
▽ More
This paper is concerned with minimization of a fourth-order linearized Canham-Helfrich energy subject to Dirichlet boundary conditions on curves inside the domain. Such problems arise in the modeling of the mechanical interaction of biomembranes with embedded particles. There, the curve conditions result from the imposed particle--membrane coupling. We prove almost-$H^{\frac{5}{2}}$ regularity of the solution and then consider two possible penalty formulations. For the combination of these penalty formulations with a Bogner-Fox-Schmit finite element discretization we prove discretization error estimates which are optimal in view of the solution's reduced regularity. The error estimates are based on a general estimate for linear penalty problems in Hilbert spaces. Finally, we illustrate the theoretical results by numerical computations. An important feature of the presented discretization is that it does not require to resolve the particle boundary. This is crucial in order to avoid re-meshing if the presented problem arises as subproblem in a model where particles are allowed to move or rotate.
△ Less
Submitted 26 September, 2017; v1 submitted 20 March, 2017;
originally announced March 2017.
-
Numerical approximation of multi-phase Penrose-Fife systems
Authors:
Carsten Gräser,
Max Kahnt,
Ralf Kornhuber
Abstract:
We consider a non-isothermal multi-phase field model. We subsequently discretize implicitly in time and with linear finite elements. The arising algebraic problem is formulated in two variables where one is the multi-phase field, and the other contains the inverse temperature field. We solve this saddle point problem numerically by a non-smooth Schur-Newton approach using truncated non-smooth Newt…
▽ More
We consider a non-isothermal multi-phase field model. We subsequently discretize implicitly in time and with linear finite elements. The arising algebraic problem is formulated in two variables where one is the multi-phase field, and the other contains the inverse temperature field. We solve this saddle point problem numerically by a non-smooth Schur-Newton approach using truncated non-smooth Newton multigrid methods. An application in grain growth as occurring in liquid phase crystallization of silicon is considered.
△ Less
Submitted 12 May, 2016;
originally announced May 2016.
-
The interface for functions in the dune-functions module
Authors:
Christian Engwer,
Carsten Gräser,
Steffen Müthing,
Oliver Sander
Abstract:
The dune-functions dune module introduces a new programmer interface for discrete and non-discrete functions. Unlike the previous interfaces considered in the existing dune modules, it is based on overloading operator(), and returning values by-value. This makes user code much more readable, and allows the incorporation of newer C++ features such as lambda expressions. Run-time polymorphism is imp…
▽ More
The dune-functions dune module introduces a new programmer interface for discrete and non-discrete functions. Unlike the previous interfaces considered in the existing dune modules, it is based on overloading operator(), and returning values by-value. This makes user code much more readable, and allows the incorporation of newer C++ features such as lambda expressions. Run-time polymorphism is implemented not by inheritance, but by type erasure, generalizing the ideas of the std::function class from the C++11 standard library. We describe the new interface, show its possibilities, and measure the performance impact of type erasure and return-by-value.
△ Less
Submitted 18 December, 2015;
originally announced December 2015.
-
A note on Poincaré- and Friedrichs-type inequalities
Authors:
Carsten Gräser
Abstract:
We introduce a simple criterion to check coercivity of bilinear forms on subspaces of Hilbert-spaces and Banach-spaces. The presented criterion allows to derive many standard and non-standard variants of Poincaré- and Friedrichs-type inequalities with very little effort.
We introduce a simple criterion to check coercivity of bilinear forms on subspaces of Hilbert-spaces and Banach-spaces. The presented criterion allows to derive many standard and non-standard variants of Poincaré- and Friedrichs-type inequalities with very little effort.
△ Less
Submitted 18 January, 2023; v1 submitted 9 December, 2015;
originally announced December 2015.
-
A Variational Approach to Particles in Lipid Membranes
Authors:
Charles M. Elliott,
Carsten Gräser,
Graham Hobbs,
Ralf Kornhuber,
Maren-Wanda Wolf
Abstract:
A variety of models for the membrane-mediated interaction of particles in lipid membranes, mostly well-established in theoretical physics, is reviewed from a mathematical perspective. We provide mathematically consistent formulations in a variational framework, relate apparently different modelling approaches in terms of successive approximation, and investigate existence and uniqueness. Numerical…
▽ More
A variety of models for the membrane-mediated interaction of particles in lipid membranes, mostly well-established in theoretical physics, is reviewed from a mathematical perspective. We provide mathematically consistent formulations in a variational framework, relate apparently different modelling approaches in terms of successive approximation, and investigate existence and uniqueness. Numerical computations illustrate that the new variational formulations are directly accessible to effective numerical methods.
△ Less
Submitted 3 May, 2016; v1 submitted 8 December, 2015;
originally announced December 2015.