We present a mathematical framework for simulation of optical fields in complex gravitational-wave interferometers. The simulation framework uses the two-photon formalism for optical fields and includes radiation pressure effects, an important addition required for simulating signal and noise fields in next-generation interferometers with high circulating power. We present a comparison of results from the simulation with analytical calculation and show that accurate agreement is achieved
