Internal API

solvepressure(electrolyte)

Check if pressure is to be solved for using the momentum balance equation, or if the pressure is obtained from the (Navier)-Stokes solver.

edgevelocity(electrolyte, iedge)

Obtain the velocity projection onto a simplex edge of index i (normals to Voronoi cell boundaries).

value_or_getindex(v,i)

Return number or i-th component of Union{Number, Vector}.

chargedensity!(q, sol,electrolyte)

Calculate charge density from solution (on the whole grid), writing the result into q and returning this vector.

wnorm(u,w,p)

Weighted norm with respect to columns

pnpreaction!(f, u, node, electrolyte)

Finite volume reaction term

pnpflux(f, u, edge, electrolyte)

Finite volume flux. It calls either μex_flux!, cent_flux! or act_flux!.

pnpstorage!(f, u, node, electrolyte)

Finite volume storage term

pnpbstorage!(f, u, node, electrolyte)

Finite volume boundary storage term

 dμex(γk, γl, electrolyte)

Calculate differences of excess chemical potentials from activity coefficients

l pbreaction!(f, u, node, electrolyte)

Reaction expression for Poisson-Boltzmann

pbflux!(f, u, edge, electrolyte)

Flux expression for Poisson-Boltzmann

splitz(range::AbstractRange)

If range contains zero, split it into two parts, one with values <=0 and one with values >=0. Otherwise, return the range or its reverse, such that first value always is the one with the smallest absolute value.

splitz(range::Vector)

Version of splitz(range::AbstractRange) for vectors.

showstruct(io::IO, this)

Print struct with field names and field values.

myround(x; kwargs...)

Rounding for use in showstruct.