Given a conic metric along a simple normal crossing divisor, we consider the conical Kahler-Ricci flow which preserves the conic structure. We use approximation method to construct the solutions on the largest time interval. Then we establish high order estimates for such flow solution. Finally we briefly introduce some convergence results.