nerplanck_numflux_concentration Subroutine

private subroutine nerplanck_numflux_concentration(left, right, diffusivity, flux)

of the Nernst-Planck equations.

This subroutine calculates the flux of the Nernst-Planck equation on the reference cubic face.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in)			::	left(4)	Left state vector (as conservative variables). The order of this vector has to be $(u, \rho_x, \rho_y, \rho_z)$ where $u$ and $\rho$ denoted concentration (scalar) and diffusive fluxes (vector).
real(kind=rk),	intent(in)			::	right(4)	Right state vector (as conservative variables). The order of this vector has to be $(u, \rho_x, \rho_y, \rho_z)$ where $u$ and $\rho$ denoted concentration (scalar) and diffusive fluxes (vector).
real(kind=rk),	intent(in)			::	diffusivity	The diffusivity in the left and right cell. Since this is assumed to be equal in both cells, this flux function cannot be used to calculate the flux for cells with different material properties
real(kind=rk),	intent(out)			::	flux	The flux between left and right cell. The order of this vector is the same as the input arguments.

Called by

Contents

Variables

diffusivitySqrt

Variables

Type	Visibility	Attributes		Name		Initial
real(kind=rk),	private		::	diffusivitySqrt