A CFD simulation solver that uses a Meshfree method and eliminates complicated grid generation processes was developed. This method uses a moving least squares interpolation in order to improve the analysis accuracy of the virtual boundary method. It was applied to the numerical simulation of high Reynolds number flow around a two-dimensional circular cylinder and a three-dimensional automobile to verify the accuracy of the flow analysis. The simulation results show good agreement with results from wind tunnel experiments.