The main objective of this thesis is to do the study of the computational fluid dynamics of the flapping wings of 20 cm wingspan by estimating three dimensional(3D) aerodynamic values along with the study of the flow field using the software COMSOL Multiphysics. The COMSOL Multiphysics is a comprehensive simulation software environment for wide array of applications designed to provide the most accurate results which gives the user access to choose most of the conditions by lowering the assumptions its user must make. Most basically the modification of the computational fluid dynamics simulation of a single flapping wing is done by altering the condition of the upper boundary from free to enclosed just like the real wind tunnel boundary condition. The study done is broadly done in three different cases for the single wing in which the lower velocity value of 1m/s and 2 m/s is set up in the laminar model setup and the higher velocity of 3 m/s is done in the turbulent model. Most importantly the comparison between the modified and the previous data is done. And later on, the study of single wing is extended to the wing to wing case which is modeled in a way that the two wings are set between the two-supporting plane in the wing tunnel and the study of the flow field along with the aerodynamic value is done. The simulation of the formation flight study will be done near in future after the successful completion of the study of the wing to wing simulation case.