A new optimization techniques based on the adjoint lattice Boltzmann method is derived for airfoil inverse design in laminar compressible flows. In this study, the developed adjoint lattice Boltzmann scheme based on the circular function (CF) is extended for airfoil inverse design problems in laminar incompressible and compressible flows. New mathematical derivation based on compressible lattice Boltzmann equations (LBE) is developed which can find target shape of an airfoil with available desired pressure distribution. The adjoint lattice Boltzmann method is extended for both the incompressible and compressible flows by selecting the circular function idea for calculating the equilibrium distribution functions. So, the adjoint equation is also expanded based on CF idea for calculation of objective function gradient vector. The steepest decent technique is utilized as gradient optimizer. Also, a novel solution is presented to remove singularity problem of the adjoint boundary condition. In order to validate the developed optimization algorithm, results are presented for both incompressible and compressible inverse problem in steady and unsteady flow and accurate results are obtained.