Open jguterl opened 5 months ago
Brendan wants just Faraday induction, not the full MHD model
We have enough information: Problem
eq1 = ∂ₜ(B, Δt) - ∇ × (η × (∇ × B)) #, :sf => Δt/B0
J = (∇ × B)
vars = variables(B)
eqs = equations(eq1)
bcs = boundary_conditions(grid.n ⋅ B - 0.0, grid.n ⋅ J - 0.0, grid.e.z ⋅ B - 0.0)