Superior Center of Technological Education 1,2 , Unicamp, Limeira/SP leobravo@ceset.unicamp.br and patroc@ceset.unicamp.br, P.O. 13484-370, Brazil, Federal Center of Technological Teaching 3 , São Luis/ MA, bzuza@dee.cefet-ma.br, P.O. 6506-500, Brazil, Fluminense Federal University 4 , Niteroi/RJ pablo@telecom.uff.br, , P.O. 24210-220, , Brazil, University of Campinas (Unicamp) 1,2,5 , Campinas/SP, P.O. 13083-970, hugo@dmo.fee.unicamp.br, Brazil. Abstract -A new and efficient spatio-temporal formulation based on the finite element method 2D is applied to simulate optical devices with cylindrical symmetry. One spatial coordinate r, radial component, and one temporal coordinate , transversal component, constitute the discretization domain. The propagation occurs along zdirection, being incremented by Crank-Nicholson scheme. The finite-element in conjunction with Galerkn's method it's been applied to discretizate the transversal domain and perfect matched layers (PML's) are placed to truncate the transversal boundaries. Index Terms -Beam-propagation method (BPM), Finiteelement method, optical divices, spatio-temporal, perfect matched layers.