cl (Chiu-Lin)

← Components

intsharp/sharpening.py

Chiu-Lin sharpening

Supports: 1D and 2D (selected automatically based on domain)

Flux formulation:

\[ \mathbf{q} = \epsilon \nabla\alpha - \alpha(1-\alpha) \hat{\mathbf{n}} \]

where \( \hat{\mathbf{n}} = \nabla\alpha / |\nabla\alpha| \) is the interface unit normal.

RHS = \( \nabla \cdot \mathbf{q} \) (divergence of the flux).

Update:

\[ \alpha^{n+1} = \alpha^n + \Delta t \cdot \Gamma \cdot \text{RHS} \]

Parameters:

  • eps_target: Target interface thickness \( \epsilon \)
  • strength: Sharpening strength \( \Gamma \)

Source

A conservative phase field method for solving incompressible two-phase flows (PDF)

BibTeX:

@article{chiu2011conservative,
  title   = {A conservative phase field method for solving incompressible two-phase flows},
  author  = {Chiu, Pao-Hsiung and Lin, Yan-Ting},
  journal = {Journal of Computational Physics},
  volume  = {230},
  number  = {1},
  pages   = {185--204},
  year    = {2011},
  publisher = {Elsevier}
}