  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.
