Intergranular or intragranular anisotropic pores can be easily observed in the FCC structure of nuclear reactor core structural materials, such as austenitic stainless steel or nickel-based alloys. Austenitic stainless steel contains a certain amount of nickel(Ni), Ni undergo transmutation reaction under neutron irradiation to produce helium. Helium combines with vacancy and continuously absorbs more helium and vacancy, evolving into underpressure pores filled with a small amount of helium. The morphology of pores is influenced by both the surface anisotropy of the crystal and grain boundary characteristic because pores nucleation predominantly occurs at grain boundaries. The swelling effect caused by pores and the embrittlement effect of high temperature helium are related to the morphology, size and distribution of pores. The phase field method can couple multiple physical fields and accurately describe the effects of material microscopic defects on pores. In this study, we use establish a free energy functional coupling between crystal plane anisotropy and pore-grain boundary interactions, using phase field methods to simulate pores’ evolution and morphology. Our results demonstrate that helium gas induces pore nucleation, with higher concentrations leading to shorter incubation period, faster nucleation rate, and greater growth rate. Grain boundaries act as heterogeneous nucleation sites for helium pores, leading to the formation of pores along these boundaries and high-density diffusion pores within grains. The intragranular pores exhibit anisotropic characteristics regulated by interfacial energy's anisotropic modulus, the strength of the anisotropy, and crystal orientation. The high-density intergranular pores interact significantly and are influenced by grain boundaries, while the anisotropic morphology is negligible. Additionally, it has been observed that the pores located at the middle of grain boundaries tend to exhibit an elliptical shape. The stress inside the pores that contain a small amount of helium gas is negative, which is lower than the value in the matrix. The findings presented herein align well with experimental results and have implications for life prediction models for service components as well as core material design.