pub fn cell_problem_1d(cell: &UnitCell, max_iter: usize, tol: f64) -> Vec<f64>Expand description
Solve the scalar periodic cell problem via finite differences to obtain the effective conductivity in the x-direction.
The cell problem: -∇·(k ∇χ) = ∂k/∂x with periodic BCs. Here we use a simple 1-D iterative solve across a column of voxels.
§Arguments
cell– Unit cell with nx × ny voxels (nz=1 assumed for 2-D).max_iter– Maximum Gauss-Seidel iterations.tol– Convergence tolerance on relative residual.