Abstract. During the last decades considerable efforts have been exerted for the development of micro air vehicles as well as microelectromechanical systems in general, for a wide range of applications. However, such systems involve microscale rarefied gas flows, which appear to be significantly different comparing to flows at macroscale and continuum regime; it is this the reason the Navier-Stokes equations fail to simulate such phenomena without further modification. To this end, the enhancement of the in-house academic Computational Fluid Dynamics solver Galatea to encounter such simulations is reported in this study. In case of rarefied gas flows and particularly for fluids in slip flow regime (Knudsen number greater than 0.01) the no-slip condition on solid wall surfaces is no longer valid; hence, velocity slip conditions as well as temperature jump ones have to be included instead. Furthermore, to increase accuracy at the same region the second-order accurate spatial slip model of Beskok and Karniadakis has been incorporated, which avoids the numerical difficulties, entailed by the evaluation of the second derivative of slip velocity when complex geometries along with unstructured hybrid grids are encountered. Due to oscillations that might appear, especially during the initial steps of the iterative procedure, a normalization scheme is additionally employed, to allow for the gradual increase of the corresponding slip/jump values. Galatea has been validated against a benchmark test case concerning rarefied laminar flow (inside the slip flow regime) over a wing with a NACA0012 airfoil in different angles of attack. The obtained results were compared with those of a reference solver, and with those obtained with the paralleld open-source kernel SPARTA, based on the Direct Simulation Monte-Carlo method. According to this last approach, the flow domain is divided into a finite number of computational cells, while the required sample macroscopic flow properties are retrieved assuming intermolecular collisions of the simulated particles inside such cells. An excellent agreement was achieved between the results obtained by Galatea and SPARTA as well.