A numerical approach to ionized nonequilibrium boundary layers