This subroutine calculates the flux of the Nernst-Planck equation on the reference cubic face.
| 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 where and 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 where and 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(3) |
The flux between left and right cell. The order of this vector is the same as the input arguments. |