Although lattice gases composed by kNN particles, forbidding up to their kth nearest neighbors of being occupied, have been widely investigated in literature, the location and the universality class of the fluid-columnar transition in the 2NN model on the square lattice is still a topic of debate. Here, we present grand-canonical solutions of this model on Husimi lattices built with diagonal square lattices, with 2L(L + 1) sites, for L 7. The systematic sequence of mean-field solutions confirms the existence of a continuous transition in this system and extrapolations of the critical chemical potential µ2,c(L) and particle density ρ2,c(L) to L → ∞ yield estimates of these quantities in close agreement with previous results for the 2NN model on the square lattice. To confirm the reliability of this approach we employ it also for the 1NN model, where very accurate estimates for the critical parameters µ1,c and ρ1,c -for the fluid-solid transition in this model on the square lattice -are found from extrapolations of data for L 6. The non-classical critical exponents for these transitions are investigated through the coherent anomaly method (CAM), which in the 1NN case yields β and ν differing by at most 6% from the expected Ising exponents. For the 2NN model, the CAM analysis is somewhat inconclusive, because the exponents sensibly depend on the value of µ2,c used to calculate them. Notwithstanding, our results suggest that β and ν are considerably larger than the Ashkin-Teller exponents reported in numerical studies of the 2NN system.