Magnetic-fusion energy (MFE) devices are large, complex systems to develop a major new consumer energy source. Many MFE plasma models incorporate large-scale nonlinear partial differential equations (PDEs) describing the confinement of hot plasma in electromagnetic fields. This presentation will introduce the use of parallel preconditioned Newton-Krylov methods in PETSc for the solution of nonlinear PDEs arising in edge plasma modeling. This parallel edge application is a key part of a coupled core-edge-wall modeling capability being developed by the Framework Application for Core-Edge Transport Simulation (FACETS) SciDAC project.