cls_2010 (CLS with mapped normal, 2010)

← Components

intsharp/sharpening.py

CLS with mapped normal (2010)

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

RHS for \( \alpha \) (non-conservative form):

\[ \text{RHS} = \mathbf{n} \cdot \nabla\bigl[\epsilon|\nabla\alpha| - \alpha(1-\alpha)\bigr] \]

Normal from a mapped field: \( \phi(\alpha) = \alpha^a / (\alpha^a + (1-\alpha)^a) \), \( \mathbf{n} = \nabla\phi/|\nabla\phi| \).

Update:

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

Parameters:

  • eps_target: Target interface thickness \( \epsilon \)
  • strength: Sharpening strength \( \Gamma \)
  • mapping_alpha: Mapping exponent \( a \) (default 2); set under sharpening.method_params in the YAML config