Field propagator

TwoFluidQuasiNeutralFull

Advances u, ue and phi.

What it solves

FEEC discretization of the following equations:

find $u \in H (div)$ , $u_{e} \in H (div)$ and $ϕ \in L^{2}$ such that

\int_{Ω} \partial_{t} u \cdot v d x 00 = \int_{Ω} ϕ \nabla \cdot v d x + \int_{Ω} u \times B_{0} \cdot v d x + ν \int_{Ω} \nabla u : \nabla v d x + \int_{Ω} f v d x \forall v \in H (div) . = - \int_{Ω} ϕ \nabla \cdot v_{e} d x - \int_{Ω} u_{e} \times B_{0} \cdot v_{e} d x + ν_{e} \int_{Ω} \nabla u_{e} : \nabla v_{e} d x + \int_{Ω} f_{e} v_{e} d x \forall v_{e} \in H (div) . = \int_{Ω} ψ \nabla \cdot (u - u_{e}) d x \forall ψ \in L^{2} .

time_discret: fully implicit.

The variables this step updates over a time step dt, and the discrete spaces each one accepts.

Set these on the propagator in a parameter file, for example TwoFluidQuasiNeutralFull.Options(nu=...).

Option	Default	Choices	Description
`nu`	`1`	float	Ion viscosity coefficient.
`nu_e`	`1`	float	Electron viscosity coefficient.
`eps_norm`	`0.001`	float	Normalization/scaling parameter in Lorentz coupling terms.
`boundary_data_u`	`None`	dict[tuple[int, int], typing.Callable] \| None	Inhomogeneous Dirichlet data for ion velocity faces.
`boundary_data_ue`	`None`	dict[tuple[int, int], typing.Callable] \| None	Inhomogeneous Dirichlet data for electron velocity faces.
`source_u`	`None`	Optional	Source term for ion momentum equation.
`source_ue`	`None`	Optional	Source term for electron momentum equation.
`stab_sigma`	`None`	float \| None	Optional stabilization coefficient for electron block.
`solver`	`gmres`	`pbicgstabbicgstabgmres`	Linear/saddle-point solver used for the global system.
`solver_params`	`None`	struphy.linear_algebra.solver.SolverParameters \| None	Solver controls.

Each model splits its time step into a sequence of propagators. This one appears in 1 model, at the position shown.