<p>In this presentation, we present a recently developed particle simulation too PIFE-PIC, which is a novel three-dimensional (3-D) Parallel Immersed-Finite-Element (IFE) Particle-in-Cell (PIC) simulation model for particle simulations of plasma-material interactions. This framework is based on the ecently developed non-homogeneous electrostatic IFE-PIC method, which is designed to handle complex plasma-material interface conditions associated with irregular geometries using a Cartesian-mesh-based PIC. Three-dimensional domain decomposition is utilized for both the electrostatic field solver with IFE and the particle operations in PIC to distribute the computation among multiple processors. A simulation of the orbital-motion-limited (OML) sheath of a dielectric sphere immersed in a stationary plasma is carried out to validate PIFE-PIC and profile the parallel performance of the code package. Parallel efficiency up to approximately 110 superlinear speedup was achieved for strong scaling test. Furthermore, a large-scale simulation of plasma charging at a lunar crater containing 2 million PIC cells (10 million FE/IFE cells) and about 1 billion particles, running for 20,000 PIC steps in about 154 wall-clock hours, is presented to demonstrate the high-performance computing capability of</p>