Struct LevelSetEvolution 

pub struct LevelSetEvolution {
    pub dt: f64,
    pub n_steps: usize,
}

Level-set evolution by advection and mean curvature flow.

Implements the upwind finite difference scheme for advection and the explicit mean-curvature-flow equation for smoothing.

Fields

dt: f64

Time step for explicit integration.

n_steps: usize

Number of time steps per evolution call.

Implementations

impl LevelSetEvolution

pub fn new(dt: f64, n_steps: usize) -> Self

Create a new LevelSetEvolution solver.

Arguments

• dt – explicit time step (should satisfy CFL: dt ≤ dx/|v|_max)
• n_steps – number of time steps per advect or mean_curvature_flow call

pub fn advect(&self, field: &mut LevelSetField, velocity: [f64; 3])

Advect the level-set field by a uniform velocity field using the upwind scheme.

Solves ∂φ/∂t + v · ∇φ = 0 with first-order upwind spatial discretisation.

Arguments

• field – level-set field (modified in place)
• velocity – uniform velocity vector [vx, vy, vz]

pub fn mean_curvature_flow(&self, field: &mut LevelSetField)

Evolve the level-set by mean curvature flow.

Solves ∂φ/∂t = κ |∇φ| where κ is the mean curvature. Uses explicit Euler integration.

Arguments

• field – level-set field (modified in place)

pub fn normal_velocity_extension( &self, field: &mut LevelSetField, normal_velocity: f64, )

Normal velocity extension: propagate the interface normal velocity outward.

For each node computes the update φ_new = φ - dt * F |∇φ| where F is a prescribed scalar normal velocity.

Arguments

• field – level-set field (modified in place)
• normal_velocity – scalar normal speed (positive = expansion)