The current research proposed the theoretical model for ship twin-propeller jet based on the axial momentum theory and Gaussian normal distribution. The twin-propeller jet model is compared to the more matured single propeller jet model with good agreement. Computational Fluid Dynamics (CFD) method is used to acquire the velocity distribution within the twin-propeller jet for understanding of flow characteristics and validation purposes. Efflux velocity is the maximum velocity within the entire jet with strong influences by the geometrical profiles of the blades. Twin-propeller jet model showed four-peaked profile at the initial plane downstream to the propeller compared to the two-peaked profile from a single-propeller. The four-peaked profile merges to be two-peaked velocity profile and then one-peaked profile due to the fluid mixing. Entrainment occurs between the ambient still water outside and the rotating flow within jet due to the high velocity gradient. The research proposes a twin-propeller jet theory with a serial of equations enabling the predictions of velocity magnitude within the jet.