Finite difference simulations of neutral gas-MHD interactions in partially ionized plasmas