AbstractNumerical simulations of star formation faces challenges including the huge spatial dynamic range and the presence of multiply coupled highly non-linear physics, such as magnetic field, supersonic turbulence, gravitation, radiation and protostellar outflow feedback. We present in this symposium our latest high resolution adaptive mesh refinement numerical simulations of the formation of filamentary dark molecular clouds from 4.55 pc size scale down to the formation of a protostellar cluster with maximum resolution at 28 AU. The physical properties of the long braided filamentary dark cloud formed in the simulation, the magnetic field properties of the cloud clumps, and the protostellar mass function in the simulations match well with the latest observations.