We assume both are correct.
div D = ρ applies to emw propagation in infinite boundless space.
And Div D = ρ + α dV/dt applies to emw scattering under Diriclet boundary conditions on the voltage V in bounded space.
But the question arises, can it be proven theoretically or experimentally?