Implicit like time discretization for the one-phase Hele-Shaw problem with surface tension

Ido Lavi,Nicolas Meunier,Olivier Pantz

2023-05-10

Abstract:In this work, we propose and compare three numerical methods to handle the one-phase Hele-Shaw problem with surface tension in dimension two by using three variational approaches in the spirit of the seminal works \cite{Otto, Gia_Otto}.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is about the numerical solution method for the one - phase Hele - Shaw problem with surface tension. Specifically, the author proposes and compares three numerical methods to deal with the one - phase Hele - Shaw problem with surface tension in two - dimensional space, and these methods are based on the idea of the variational method.

### Problem Background

The Hele - Shaw problem is a classical problem that describes the movement of fluids in narrow gaps (such as between two parallel glass plates). Due to the special geometric structure, the movement of the fluid is strongly damped, and it is assumed to be governed by Darcy's law. The effect of surface tension makes the droplet have a circular cross - section at rest.

### Mathematical Model

Mathematically, the one - phase Hele - Shaw problem with surface tension can be represented as the following continuous - time moving boundary - value problem:

\[
\begin{cases}
u+\nabla p = 0 & \text{in }\Omega(t)\\
\Delta p = 0 & \text{in }\Omega(t)\\
p=\sigma\kappa & \text{on }\Gamma(t)
\end{cases}
\]

where:
- \(u(t, x)\) is the velocity field,
- \(p(t, x)\) is the pressure field,
- \(\Omega(t)\) is the region that describes the droplet,
- \(\Gamma(t)=\partial\Omega(t)\) is the boundary of the droplet,
- \(\kappa\) is the mean curvature of the free boundary (positive for a circle),
- \(V_{n}\) and \(n\) are the normal velocity field on the boundary and the outer unit normal vector respectively,
- \(\sigma > 0\) is a given surface tension constant.

In addition, the moving boundary condition is also satisfied:

\[
V_{n}=V_{\Gamma(t)}\cdot n = u\cdot n\quad \text{on }\Gamma(t)
\]

### Research Objectives

The main objective of the paper is to numerically discretize the above problem through three different variational methods and compare the effects of these methods. Specifically, the author proposes the following methods:
1. **Explicit Time Discretization**: Directly start from the known geometric configuration and explicitly handle the surface tension term.
2. **Boundary Variation Problem**: Solve the problem by variational minimization of the boundary.
3. **Minimization Problem with Curl Penalty**: Introduce a curl penalty term to handle the incompressible constraint.
4. **Fully Nonlinear Problem**: Consider the fully nonlinear incompressible constraint.

### Significance

The solution of these problems helps to more accurately simulate and understand complex fluid - dynamic behaviors, especially in cases involving interface dynamics and surface tension. This is of great significance for many practical applications (such as microfluidic technology, biomedical engineering, etc.).