From electric field of a point charge: $$ \vec{E} = \frac{k Q \vec{r}}{\gamma^{2}r^3(1-\beta^2sin^2\theta)^{\frac{3}{2}}}, \vec{B} = \frac{\vec{u} \times \vec{E}}{c^2} $$ taking curl of B gives $$ \nabla \times B = {\mu_0J}$$ Taking divergence of above gives $$ \nabla . J = \frac{\nabla .(\nabla \times B)}{\mu_0} =0 = -\frac{\partial q}{\partial t} $$ which implies stationary charge distribution despite the point charge being in motion. Expressions for E and B derived using special relativity where the charge Q moves across space with constant velocity to inertial observer. So why exactly does the expressions for E and B for moving field fail?
I have found answers on this post claim it is due to not considering speed of light limit but assuming we let enough time for the charge in its rest frame to produce electric field in all of its space, it doesn't seem to be a factor since we only need to lorentz transform four force experienced in charge's rest frame to get above same equation for E and B, where the charge will be in motion. Other answer on the same post attribute the error to lack of consideration of charge conservation which if we let q in its rest frame be constant (avoids need to re-consider speed of light limit) would also be constant in inertial frame where the charge is moving, by invariance of charge (already assumed for the derivation of above equations).