We develop and validate a three-dimensional particle method based on an Enhanced Weakly-Compressible MPS approach for modeling immersed dense granular flows. For this purpose, we adopt a generalized rheological model, using a regularized visco-inertial rheology, for all regimes of multiphase granular flow. Moreover, we propose a new consistent formulation to estimate the effective pressure of the solid skeleton based on the continuity equation of the pore-water. To improve the accuracy of the multiphase particle methods, especially near boundaries and interfaces, we introduce a modified high-order diffusive term by employing the convergent form of the Laplacian operator. The effectiveness of the new diffusive term is particularly demonstrated by modeling the hydrostatic pressure of two fluid phases. Further, coupling the generalized rheology model with the flow equations, we investigate the gravity-driven granular flows in the immersed granular collapse and slide in three dimensions. As a part of this study, we represent the experiment on the immersed granular collapse to validate the model. The evolution and runout length of the granular bulk are compared with those from experiments confirming good compatibility. Overall, the qualitative and quantitative results justify the proposed developments shown to be essential for predicting different states of the immersed granular flows.