Traditional precise orbit determination (POD) for low Earth orbit (LEO) satellites relies on observations from ground stations and onboard receivers. Although the accuracy can reach centimeter level, there are still problems such as insufficient autonomous operation capability. The inter-satellite link (ISL) is a link used for communication between satellites and has a function of dual-way ranging. Numerous studies have shown that observational data using ISLs can be adopted for POD of navigation satellites. In this contribution, we mainly focus on LEO satellites POD with ISLs. First, we design LEO constellations with different numbers of satellites and ISL measurements, based on which the constellations are simulated. Then rough tests of POD using different link topologies are carried out. The results show that in the 60-LEO constellation the average 3-dimensional (3D) orbital errors are 0.112 m using “4-connected” link topology with constant 4 links per satellite and 0.069 m using “all-connected” link topology with theoretically maximum numbers of links. After that, we carry out refined POD experiments with several sets of satellite numbers and different observation accuracy. The results show the higher link ranging accuracy and the more numbers of links bring higher POD precision. POD with ISLs gets bad performance in the case of center of gravity reference when link ranging accuracy is poor and numbers of links are small. When the link accuracy is 40 cm, average 3D orbital errors of 60-LEO constellation are 0.358 m, which can only meet the demand of autonomous navigation. With the constraint of the right ascension of the ascending node (RAAN), POD using ISLs reaches an extremely high precision when adopting a spatial reference provided by navigation satellites. For 120-LEO constellation, the average 3D orbital errors are 0.010 m; for 192-LEO constellation, the errors are 0.006 m.