In order to establish the next-generation reactor physics calculation method based on the numerical nuclear reactor technology and realize high-fidelity modeling and calculation, a new numerical nuclear reactor neutronics code SHARK is developed. The code is based on the direct transport method with construct solid geometry (CSG) method, advanced subgroup resonance method, direct transport MOC method in rectangle and hexagonal geometry, large-scale parallel, and CMFD acceleration method. The C5G7, macro BEAVRS and VERA benchmarks are verified to show the accuracy of the code and method. Numerical results show good accuracy and calculation performance of SHARK, and the direct transport method can be adopted on numerical nuclear reactor calculation.