Numerical modeling of the SNS H− ion source