pressure_eq Subroutine

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

Outlet Pressure Equilibrium type boundary conditions

The incoming densities are set to the equilibrium distribution with macroscopic values used from: * Velocity: Linear extrapolation of velocity of precollision of current time step * Pressure as defined in the configuration file

Usage:\n

  boundary_condition = {
    { label = 'outlet',
       kind = 'pressure_eq',
       pressure = 'pressure_out',
    }
  }
 variable = {
   name = 'pressure_out',
   ncomponents = 1,
   vartype = 'st_fun',
   st_fun = 1.0
 }

This is described in the paper:\n S. Izquierdo, P. Martinez-Lera, and N. Fueyo,\n Analysis of open boundary effects in unsteady lattice Boltzmann simulations Comput. Math. Applications, vol. 58, no. 5, pp. 914–921, Sep. 2009.\n

The velocity is extrapolated (1st order) from the flow domain See equation (6)

DISCLAIMER: This BC requires pre-collision PDF from 1st and 2nd neighbor which are not available if those neighbors are boundary halo on different process. So, if you encounter a problem with this BC, please check the domain distribution to verify whether 2 neighbors of boundary elements are in one process. The domain distribution can be changed by shifting the boundary by one element or by changing the origin of bounding cube or by using different number of MPI processes.

\todo 20200410, KM: Check if this BC works if neighbufferPre is replaced by neighbufferPost

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

Arguments