Volumetric scattering prevents imaging modalities in biomedical optics from imaging deep inside tissue. The optimization of the incident wavefront has the potential to improve these imaging modalities. To investigate the optimization and light propagation of such beams inside scattering media rigorously, full-vectorial simulations based on solutions of Maxwell’s equations are necessary. In this publication, we present a versatile two-step beam synthesis method to efficiently simulate the scanning and phase optimization of a focused beam inside a static scattering medium. We present four different approaches to the phase optimization of the energy density and the absolute value of the Poynting vector. We find that these quantities have two regions with different, almost exponential decays over depth for a non-optimized beam. Optimization by conjugating the phase of the projected electric field in various directions at the focus shows an improvement below a certain penetration depth. Seeking global solutions to the optimization problems reveals an even better enhancement in the energy density and the absolute value of the Poynting vector in the focus. For Poynting vector optimization, the differences between the presented optimization approaches are more significant than for the energy density. With the presented method, it is possible to efficiently simulate different imaging methods improved by wavefront shaping to investigate their possible penetration depths.