The evaporation flux distribution of sessile drops is investigated by molecular dynamic simulations. Three evaporating modes are classified, including the diffusion dominant mode, the substrate heating mode, and the environment heating mode. Both hydrophilic and hydrophobic drop-substrate interactions are considered. To count the evaporation flux distribution, which is position dependent, we proposed an azimuthal-angle-based division method under the assumption of spherical crown shape of drops. The modeling results show that the, edge evaporation, i.e., near the contact line, is enhanced for hydrophilic drops in all the three modes. The surface diffusion of liquid molecular absorbed on solid substrate for hydrophilic cases plays an important role as well as the space diffusion on the enhanced evaporation rate at the edge. For hydrophobic drops, the edge evaporation flux is higher for the substrate heating mode, but lower than elsewhere of the drop for the diffusion dominant mode; however, a nearly uniform distribution is found for the environment heating mode. The evidence shows that the temperature distribution inside drops plays a key role in the position-dependent evaporation flux.