Thermodynamically consistent and positivity-preserving discretization of the thin-film equation with thermal noise

by   Benjamin Gess, et al.

In micro-fluidics not only does capillarity dominate but also thermal fluctuations become important. On the level of the lubrication approximation, this leads to a quasi-linear fourth-order parabolic equation for the film height h driven by space-time white noise. The gradient flow structure of its deterministic counterpart, the thin-film equation, which encodes the balance between driving capillary and limiting viscous forces, provides the guidance for the thermodynamically consistent introduction of fluctuations. We follow this route on the level of a spatial discretization of the gradient flow structure. Starting from an energetically conformal finite-element (FE) discretization, we point out that the numerical mobility function introduced by Grün and Rumpf can be interpreted as a discretization of the metric tensor in the sense of a mixed FE method with lumping. While this discretization was devised in order to preserve the so-called entropy estimate, we use this to show that the resulting high-dimensional stochastic differential equation (SDE) preserves pathwise and pointwise strict positivity, at least in case of the physically relevant mobility function arising from the no-slip boundary condition. As a consequence, this discretization gives rise to a consistent invariant measure, namely a discretization of the Brownian excursion (up to the volume constraint), and thus features an entropic repulsion. The price to pay over more naive discretizations is that when writing the SDE in Itô's form, which is the basis for the Euler-Mayurama time discretization, a correction term appears. To conclude, we perform various numerical experiments to compare the behavior of our discretization to that of the more naive finite difference discretization of the equation.



There are no comments yet.


page 1

page 2

page 3

page 4


A linear implicit Euler method for the finite element discretization of a controlled stochastic heat equation

We consider a numerical approximation of a linear quadratic control prob...

On a deterministic particle-FEM discretization to micro-macro models of dilute polymeric fluids

In this paper, we propose a deterministic particle-FEM discretization to...

An adaptive time-stepping full discretization for stochastic Allen–Cahn equation

It is known in [1] that a regular explicit Euler-type scheme with a unif...

Free energy diminishing discretization of Darcy-Forchheimer flow in poroelastic media

In this paper, we develop a discretization for the non-linear coupled mo...

Structure preserving discretization of time-reparametrized Hamiltonian systems with application to nonholonomic mechanics

We propose a discretization of vector fields that are Hamiltonian up to ...

An Introduction to Multi-Point Flux (MPFA) and Stress (MPSA) Finite Volume Methods for Thermo-Poroelasticity

We give a unified introduction to the MPFA- and MPSA-type finite volume ...

Q-tensor gradient flow with quasi-entropy and discretizations preserving physical constraints

We propose and analyze numerical schemes for the gradient flow of Q-tens...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.