In transonic interfacial crack propagating fracture problem, the generation-phase simulations were done using the moving finite element method based on Delaunay automatic mesh generation. And the contact function based on the penalty method was newly developed to consider the crack face contact near the propagating interfacial crack tip. It was succeeded to visualize in 3-dimensions the Mach shock wave emanated from the propagating crack tip. And it was tried for the transonically propagating crack problem that solving the energy flows through the contact zone or along the Mach shock wave line emitted from the crack tip. The energy flow patterns into the crack tip were also visualized. Furthermore, from the values of the separated dynamic J integrals, it was found that the dynamic J integral is non-zero even for transonic fracture region and the most of the energy release rate is provided from the more compliant material epoxy.