Fracture prediction in Orthotropic Functionally Graded Materials (FGM) poses a formidable challenge due to the intricate interplay of spatially varying material properties. Orthotropic FGMs, characterized by distinct material behaviors along different axes, introduce heightened complexity in modeling, particularly when addressing fracture phenomena. The need for accurate and computationally efficient fracture models becomes crucial for understanding and predicting the structural response of such materials. In this context, phase field modeling emerges as a highly feasible approach, offering an effective means to simulate crack initiation, propagation, and their intricate interactions in a numerically stable manner. This paper details implementing a phase field model for fracture prediction in Orthotropic Functionally Graded Materials (FGM) within the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source finite element framework. Leveraging MOOSE's flexibility, the study incorporates a hybrid split of strain energy proposed by Ambati and introduces an adaptive mesh refinement scheme to enhance computational efficiency. The primary emphasis is on the meticulous implementation aspects, showcasing how MOOSE accommodates the complexities of orthotropic materials. Through reproducing crack paths from literature using in-house codes, the implemented model's reliability is demonstrated, validating its accuracy. The innovative use of adaptive mesh refinement not only reduces computational costs significantly but also highlights the efficiency gains achieved in simulating fracture in orthotropic FGMs. The paper underscores the advantages of utilizing MOOSE as an open-source finite element framework for transparency, collaboration, and accessibility within the research community, showcasing its role in advancing accurate and efficient fracture predictions.