A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma-wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remains an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.