This paper proposes a new synchronous machine model based on the application of Park's transformation to the discretized equations of the phase-domain model. The proposed approach maintains the precision of the phase-domain model and eliminates its computational inefficiencies through a constant admittance matrix. The model includes magnetic saturation. Precision analysis is performed within the accuracy constraints of the surrounding network.