Lie symmetry analysis, new group invariant for the (3 + 1)-dimensional and variable coefficients for liquids with gas bubbles models