  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.
