getEquilibrium_forElemfromState Function

private pure function getEquilibrium_forElemfromState(state, elem, layout, varPos, nScalars, neigh) result(equil)

Calculate the equilibrium distribution function in all directions

The equilibrim distribution function is:\n $$f^{eq}_i = w_i \rho ( 1 + \frac{\vec c_i \cdot \vec u}{c^2_s} + \frac{ {(\vec c_i \cdot \vec u)}^2}{2c^4_s} - \frac{\vec u \cdot \vec u}{2c^2_s})$$ \n where $w_i$ is the weight in each direction,\n $\rho$ is the macroscopic value of density,\n $c_s$ is the speed of sound,\n $\vec c_i$ is the lattice unit velocity in each direction,\n $\vec u$ is the macroscopic value of velocity.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=rk),	intent(in)			::	state(:)
integer,	intent(in)			::	elem
type(mus_scheme_layout_type),	intent(in)			::	layout
integer,	intent(in)			::	varPos(:)
integer,	intent(in)			::	nScalars
integer,	intent(in)			::	neigh(:)

Return Value real(kind=rk)(layout%fStencil%QQ)

Contents

Variables

rho vel ucx usq iDir

Variables

Type	Visibility	Attributes		Name		Initial
real(kind=rk),	private		::	rho
real(kind=rk),	private		::	vel(3)
real(kind=rk),	private		::	ucx
real(kind=rk),	private		::	usq
integer,	private		::	iDir