velocity_bounceback_incomp Subroutine

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

Inlet Velocity Bounce Back boundary condition with qvalues

This is taken from the paper: [1] S. Izquierdo and N. Fueyo, "Characteristic non-reflecting boundary conditions for open boundaries in lattice Boltzmann methods," Physical Review E, vol. 78, no. 46707, 2008 The incoming densities $\bar{f_{alpha}}$ are reconstructed from the already reflected densities, enforcing the given velocity me\%ubb\%velocity. [2] M. Junk and Z. Yang. One-point boundary condition for the lattice Boltzmann method, Physical Review E, vol 72, issue 8, year 2005.

This boundary condition has error of 1st order pressure and 2nd order velocity only if wall is located exactly at q=1/2 else the accuracy of both pressure and velocity reduce by order 1.

Usage

boundary_condition = {
  { label = 'inlet',
    kind = 'velocity_bounceback',
    velocity = 'inlet_vel',
  }
}
variable = {
  name = 'inlet_vel',
  ncomponents = 3,
  vartype = 'st_fun',
  st_fun = {0.06, 0.0, 0.0}
}

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