potential_neumann Subroutine

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

Linkwise neumann potential non-equilibrium boundary condition for curved and straight walls (zero gradient). For straight wall, values are extrapolated along boundary normal instead of along the link. The accuracy of straight wall depends on the qVal defined in config file and default is set to 0.5

The pdf is decomposed into equilibrium (eq) and non-equilibrium (neq) part: - 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 - boundary potential (pot_b) has to be extrapolated using two fluid neighb: pot_b = (4*pot_f - pot_ff)/3 (typo in paper) - Linkwise: as much as possible is outsourced to the subroutine "mus_set_nonEqExpol"

NOTE: - possibility to extend the equation to extraplotate pot_b by a gradient if a non-zero one is desired (see: Huang H, Lee T S, and Shu C. “Thermal curved boundary treatment for the thermal lattice Boltzmann equation,” Int. J. Mod. Phys. C, vol. 17(05), pp. 631-643, 2006) - more accurate schemes to determine pot_b can be found in: Chen Q, Zhang X, and Zhang J. “Improved treatments for general boundary conditions in the lattice Boltzmann method for convection- diffusion and heat transfer processes,” Phys Rev E, vol. 88(3), 033304, 2013

Usage

boundary_condition = {
  { label = 'outer',
    kind = 'potential_neumann_link',
    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