This paper investigates distributed beam focusing for coordinated satellite constellations with phased arrays, motivated by future non-terrestrial network (NTN) systems. A geometric and channel model is developed by incorporating satellite positions, array orientations, antenna directivity, and polarization effects. Under ideal synchronization, the achievable coherent combining gain is analyzed for different constellation geometries, showing that maximum ratio transmission (MRT) enables quadratic scaling of the received power with the number of satellites. The impact of phase errors caused by residual synchronization, timing, mobility, and localization mismatches is then investigated. Closed-form expressions for the average coherent gain are derived for uniformly distributed timing offsets, demonstrating the transition from coherent to non-coherent combining. The results show that synchronization and timing mismatches reduce the coherent combining gain, while geometry dependent effects govern the resulting spatial focusing behavior. Numerical results further show that linear and circular constellations provide different focusing characteristics and spatial separation capabilities. However, MRT-based focusing results in strong sidelobes and limited spatial division capability, motivating the need for joint analog beamforming and digital precoding optimization to improve spatial selectivity and robustness against mobility and localization errors.