The numerical simulation of flow around a three dimensional moving body faces different problems in several methods, such as disruption of the structure of the grid, the need for deletion and insertion of nodes, interpolation, and data transfer between different parts of grid. In order to tackle the above-mentioned problems, a new configuration has been developed for meshing domain, which besides providing the body with the capability of rotational and oscillatory motions in large displacements, saves the grid's primitive quality. In the introduced method, the grid connections are manipulated with the motion of the body, but the general form of the grid is not changed or disrupted. This needs a special form for nodes of the grid, which is explained in this paper. The three dimensional unsteady form of the Euler equations is solved and the properties over each cell faces are evaluated using an averaging method. For time integration of the equations an implicit dual time method is used. It can prove that the volume of all elements is constant in the introduced grid. Therefore, there is no need to calculate elements volume in every time step. Several test cases are solved and the results are compared with experimental or other numerical data.