Reference

Abstract supertype for advection/integration algorithms.

Uniform 2D Cartesian grid descriptor.

Build a 2D Cartesian grid from domain limits and cell counts.

Uniform 3D Cartesian grid descriptor.

Build a 3D Cartesian grid from domain limits and cell counts.

ELVIRA-style normal selected by local volume-fraction stencil scoring. Currently implemented for 2D split grids.

Height-function normal. Currently implemented for 2D split grids.

Abstract supertype for interface representations.

First-order Lie directional splitting advection scheme.

Abstract supertype for normal-vector reconstruction algorithms.

Piecewise-linear interface representation.

Composable reconstruction selector.

Abstract supertype for interface-reconstruction algorithms.

Youngs gradient direction reduced to an axis-aligned SLIC interface.

2D split-VOF state container.

3D split-VOF state container.

Simulation parameters for split geometric VOF on a Cartesian grid.

Key fields:

• nx, ny, nz: number of cells
• xlim, ylim, zlim: domain bounds
• cfl: target CFL for compute_dt/integrate!
• tol: interface tolerance for empty/full detection
• reconstruction, advection: algorithm selectors

Directional Strang splitting advection scheme.

Youngs volume-fraction gradient normal.

Youngs gradient-based normal reconstruction with PLIC volume enforcement.

Advance one split-advection step in 2D using first-order Lie splitting.

Advance one split-advection step in 3D using first-order Lie splitting.

Periodic index wrap helper: maps indices outside 1:n back into range.

Return (xL, xR, yB, yT) bounds of 2D cell (i, j).

Return (xL, xR, yB, yT, zD, zU) bounds of 3D cell (i, j, k).

Return center coordinates of a 2D cell as SVector{2}.

Return center coordinates of a 3D cell as SVector{3}.

Build a circular implicit function for 2D initialization.

Alias for compute_dt.

Compute a stable time step from face-velocity CFL constraints.

Supports velocity call signatures (x,t), (x,y,t)/(x,y,z,t), (x), or coordinate-only variants.

Return spatial dimension (2 or 3) of a split-VOF state.

Compute ELVIRA-style interface normals in 2D and store in nxint, nyint.

Return the volume-fraction array stored in the split-VOF state.

Compute height-function interface normals in 2D and store in nxint, nyint.

Initialize 2D volume fractions from an implicit level-set function.

f can accept either (x, y) or SVector{2} and is interpreted as:

• f > 0: inside liquid
• f < 0: outside liquid

Initialize 3D volume fractions from an implicit level-set function.

f can accept either (x, y, z) or SVector{3} and is interpreted as:

• f > 0: inside liquid
• f < 0: outside liquid

Integrate to tf using algorithm choice in params.

Integrate from current time vg.t to final time tf.

Return true if cell is empty (α <= tol).

Return true if cell is full (α >= 1 - tol).

Return true if cell is mixed (tol < α < 1 - tol).

L1 error of the volume-fraction field against a reference field.

L2 error of the volume-fraction field against a reference field.

L∞ error of the volume-fraction field against a reference field.

Reconstruct 3D PLIC planes (cplic) from the current volume fractions.

Reconstruct 2D PLIC planes (cplic) from the current volume fractions.

Dispatch reconstruction based on params.reconstruction.

Compute bulk shape descriptors from the volume-fraction field.

Returns a named tuple:

• volume
• centroid
• second_moments (central second-moment tensor)

Build a spherical implicit function for 3D initialization.

Allocate and initialize a split-VOF state from SplitVOFParams.

Returns SplitVOFGrid when nz <= 1, otherwise SplitVOFGrid3D.

Allocate a 2D split-VOF state from an existing Cartesian grid.

Allocate a 3D split-VOF state from an existing Cartesian grid.

Advance one split-advection step in 2D using Strang splitting.

Advance one split-advection step in 3D using symmetric split sweeps.

Advance one advection step using algorithm choice in params.

Alias for step!.

Alias for splitgrid(params).

Total material volume (area in 2D) represented by the volume-fraction field.

Compute Youngs interface normals in 3D and store in nxint, nyint, nzint.

Compute Youngs interface normals in 2D and store in nxint, nyint.

Build the standard Zalesak slotted-disk implicit function.