mus_PBLinear_advRel_generic Subroutine

public subroutine mus_PBLinear_advRel_generic(fieldProp, inState, outState, auxField, neigh, nElems, nSolve, level, layout, params, varSys, derVarPos)

Advection relaxation routine for the linear poisson boltzmann equation with an explicit calculation of all equilibrium quantities. Slow and simple. $$\nabla^2 \phi = k^2 \phi$$ Where k^2 is inverse of debye length and in this kernel refered as RHS_coeff $$k^2 = \sum_i \frac{ c_{\infty}z_i^2 e^2}{\epsilon_r \epsilon_0 k_b T}$$

This subroutine interface must match the abstract interface definition kernel in scheme/mus_scheme_type_module.f90 in order to be callable via compute function pointer.

\todo KM 20170821: change omega for multilevel Calculate equilibrium distribution functions fEq

Calculate source term

Generic fetching step: Streaming for pull Local copy for push

Relaxation

Arguments

Type	Intent	Optional	Attributes		Name
type(mus_field_prop_type),	intent(in)			::	fieldProp(:)	Array of field properties (fluid or species)
real(kind=rk),	intent(in)			::	inState(nElems*varSys%nScalars)	input pdf vector
real(kind=rk),	intent(out)			::	outState(nElems*varSys%nScalars)	output pdf vector
real(kind=rk),	intent(inout)			::	auxField(nElems*varSys%nAuxScalars)	Auxiliary field computed from pre-collision state Is updated with correct velocity field for multicomponent models
integer,	intent(in)			::	neigh(nElems*layout%fStencil%QQ)	connectivity vector
integer,	intent(in)			::	nElems	number of elements in state Array
integer,	intent(in)			::	nSolve	number of elements solved in kernel
integer,	intent(in)			::	level	current level
type(mus_scheme_layout_type),	intent(in)			::	layout	current layout
type(mus_param_type),	intent(in)			::	params	global parameters
type(tem_varSys_type),	intent(in)			::	varSys	variable system definition
type(mus_derVarPos_type),	intent(in)			::	derVarPos(:)	position of derived quantities in varsys for all fields

Contents

Variables

iElem iDir QQ nScalars pdfTmp potential rhs_coeff eqState source omega pot_diff

Variables

Type	Visibility	Attributes		Name		Initial
integer,	private		::	iElem
integer,	private		::	iDir
integer,	private		::	QQ
integer,	private		::	nScalars
real(kind=rk),	private		::	pdfTmp
real(kind=rk),	private		::	potential
real(kind=rk),	private		::	rhs_coeff
real(kind=rk),	private		::	eqState
real(kind=rk),	private		::	source
real(kind=rk),	private		::	omega
real(kind=rk),	private		::	pot_diff