potential_nonEqExpol_curved Subroutine

public subroutine potential_nonEqExpol_curved(me, state, bcBuffer, globBC, levelDesc, tree, nSize, iLevel, sim_time, neigh, layout, fieldProp, varPos, nScalars, varSys, derVarPos, physics, iField, mixture)

Linkwise Dirichlet potential non-equilibrium boundary condition for curved wall

The pdf is decomposed into equilibrium (eq) and non-equilibrium (neq) part: f = f_eq + f_neq - f_eq is calculated by weighting a fictitious potential, which is obtained by an extrapolation using the fluid neighbor(s) - f_neq is approximated by second-order extrapolation using the fluid neighbor(s) - for qVal < 0.75 even the second neighbor is used for the extrapolations - Dirichlet: potential on wall is directly known - Linkwise: as much as possible is outsourced to the subroutine "mus_set_nonEqExpol"

Usage

boundary_condition = {
  { label = 'inner',
    kind = 'potential_noneq_expol',
    potential = pot_inner,
    curved = true
  }
}

This is described in the paper: Luo, K.; Wu, J.; Yi HL. & Tan HP. (2016). A lattice Boltzmann method for electric field-space charge coupled problems. Proceedings of the 2016 Electrostatics Joint Conference (June 2016).

This subroutine's interface must match the abstract interface definition boundaryRoutine in bc/mus_bc_header_module.f90 in order to be callable via fnct function pointer.

Arguments