We present a detailed computer simulation study of the formation and evolution of topological defects in thin biaxial nematic films. We have used the Straley generalized Hamiltonian and a wide range of biaxial parameters to generate and explore a variety of mesophases. The differences obtained for the various cases are analyzed for thin biaxial nematic films with random planar (Schlieren) surface alignment.