Advances in time domain numerical modeling for electromagnetic wave propagation in bi-isotropic media