Modeling Methods and Pulse Interactions in Epsilon-Near-Zero Multilayer