volterra-solver

Beris-Edwards nematohydrodynamics solver for active nematics on 3D Cartesian grids.

Part of the volterra workspace.

Overview

volterra-solver implements the equations of motion for active liquid crystals. The Q-tensor evolves under the Beris-Edwards equation with co-rotation, strain coupling, and Landau-de Gennes molecular field. An incompressible velocity field is obtained via an FFT-based Stokes pressure solve. Concentration dynamics follow Cahn-Hilliard with exponential time differencing (ETD).

Topological disclination lines (charge +1/2 and -1/2) are detected by computing the holonomy of the nematic frame field around mesh plaquettes, using cartan-geo.

Modules

Example

use volterra_solver::run_dry_active_nematic;
use volterra_core::ActiveNematicParams;

let params = ActiveNematicParams { ..Default::default() };
let snapshots = run_dry_active_nematic(&params, 10_000)?;

License

MIT