ZEMAX has the very powerful, and often overlooked, ability to accurately calculate Frustrated Total Internal Reflection (FTIR) in Non-Sequential Mode.  Through the use of coatings, FTIR across relatively complex interface shapes can be calculated.