{"references": ["W. M. Stacey, Nuclear Reactor Physics, Wiley-VCH, Weinheim,\nGermany (2004)", "H. Finnemann, F. Benneritz and M. R. Wagner, \"Interface current\ntechniques for multidimensional reactor calculations\" Atomkernenergie\n30, 123 (1977)", "F. C. Silva, A. C. M. Alvim and A. S. Martinez, \"An Alternative Solver\nfor the Nodal Expansion Method Equations\" PHYSOR 2010 - Advances\nin Reactor Physics to Power the Nuclear Renaissance, Pittsburgh,\nPennsylvania, USA (2010)", "A. S. Martinez, V. Pereira and F. C. da Silva, \"A system for the\nprediction and determination of subcritical multiplication condition\"\nKerntechinik 64, 230 (1999)", "M. R. Wagner, K. Koebke and H. J. Winter, \"A nonlinear extension of\nthe nodal expansion method\" Proceedings of ANS/ENS International\nTopical Meeting, Munich, Germany, 2, 43 (1981)", "H. D. Fischer and H. Finnemann, \"The Nodal Integration Method - A\nDiverse Solver for Neutron Diffusion problems\" Atomkernenergie-\nKerntechnik, 30, 229 (1981)", "H. K. Joo, C. H. Kim, J. M. Noh and S. Kim, \"A Semi-Analytic\nMultigroup Nodal Method\" Annals of Nuclear Energy, 26, 699 (1999)", "B. Christensen, \"Three-dimensional static and dynamic reactor\ncalculations by the nodal expansion method\" RISO National Laboratory\nDK-4000, Roskild, Denmark (1985)", "A. C. M. Alvim, F. C. Silva and A. S. Martinez, \"Depletion Calculation\nfor a Nodal Reactor Physics Code\" ICONE 18 - 18th International\nConference on Nuclear Engineering, Xi'an, China (2010)"]}

For a quick and accurate calculation of spatial neutron distribution in nuclear power reactors 3D nodal codes are usually used aiming at solving the neutron diffusion equation for a given reactor core geometry and material composition. These codes use a second order polynomial to represent the transverse leakage term. In this work, a nodal method based on the well known nodal expansion method (NEM), developed at COPPE, making use of this polynomial expansion was modified to treat the transverse leakage term for the external surfaces of peripheral reflector nodes. The proposed method was implemented into a computational system which, besides solving the diffusion equation, also solves the burnup equations governing the gradual changes in material compositions of the core due to fuel depletion. Results confirm the effectiveness of this modified treatment of peripheral nodes for practical purposes in PWR reactors.