In the present study, the problems of the gravitational settling of a single porous particle and a pair of porous particles in a vertical channel have been numerically simulated using the multiple relaxation time lattice Boltzmann method (MRT-LBM). The reason behind using the MRT-LBM is twofold. First, with the standard MRT-LBM, there is no need for employing the moving and body-fitted grids, which are computationally demanding. Second, the momentum exchange method specific to the LBM employed here calculates the hydrodynamic forces using the mesoscopic information without any need to deal with the stress integration process. As such, four studies were designed here, the first two of which are related to the settling of a single porous particle, released from on-axis and off-axis positions, respectively. The remaining two investigations pertain to the concurrent settling of a pair of permeable particles. These particles are released either from positions aligned with the central axis but separated vertically, or from off-axis positions with both vertical and horizontal offsets. The outcomes of the initial study reveal a noteworthy enhancement in the terminal settling velocity as porosity increases. This escalation in porosity corresponds to a rise in permeability, as dictated by the Carman-Kozeny equation. For the second study, the findings show the presence of distinct settling modes for porous particles across various porosity ranges. These modes comprise “steady equilibrium with transient overshoot”, “steady equilibrium with monotonic approach”, and a damped oscillatory mode. The third and fourth studies focus on the Drafting-Kissing-Tumbling (DKT) behavior exhibited by two porous particles during concurrent settling. In the third study, it was noted that the particles consistently remain aligned with the symmetry axis. For porosity values up to 0.9, the particles exhibit Drafting-Kissing-Separating (DKS) behavior, analogous to DKT behavior minus the tumbling ste