The development, verification, and validation of a n ew Eulerian-Lagrangian DNS solver MGLET-LaParT for solving particle-laden flows is presented in this work. The proposed solver is applicable to perform one-way coupled simulations of dilute suspensions of small, inertial particles in laminar/turbulent carrier flow. We also discuss various parallelization techniques implemented using the Message Passing Interface (MPI) in MGLET-LaParT and assess their computational efficiency for performing resource-intensive Eulerian-Lagrangian simulations of particle-laden flows. MGLET-LaParT is employed to analyze dilute particulate suspension flows to understand the dispersion of suspended particles in various flow configurations. The present work investigates planar shear inflow and proximity interference effects in bluff body wake flows on the dispersion of tiny inertial particles over a range of responses to the carrier fluid. Three different classes of particles are tracked in each case, where each class is characterized by its Stokes num ber (St), namely, St= 0.1, 1, and 10. We examine the dispersion behavior of particles using instantaneous (e.g., Voronoi diagram) and mean statistical quantities (e.g., mean local concentration, velocity statistics, etc.). The effect of planar shear inflow on the dispersion of inertial particles in the unsteady laminar wake characterized by Karman vortex street behind a normal flat plate is investigated. To explore the effect of shear, we perform a comparative assessment between two different cases, distinguished based on their inflow profiles, namely, uniform and linear shear flow (with shear parameter K=0.1). We find that shear in the inflow, albeit characterized by a moderate shear parameter (K=0.1), significantly changes the distribution and dispersion of particles, especially at streamwise positions away from the plate. At such positions, the impact of shear on the concentration of St = 1 particles is the most significant among all classes, where the mean particle concentrat ion near the wake centerline is significantly increased. Furthermore, local peaks in particle concentration profiles and minima in particle velocity profiles shift laterally in their positions towards the low-stream region due to the shear. We also investigate the dispersion of tiny inertial particles in the two-dimensional laminar wake from a pair of square cylinders placed side-by-side to each other. The flow Reynolds number (defined based on the cylinder sh and uniform inflow velocity Uc) in our study is fixed at 75. We examine the effect of varying the cross-stream spacing between the cylinders on the body impaction and dispersion of particles in the wake flow over a range of their responses. Three different spacing ratios, s/h = 0.3, 2, and 4 are considered. We find that the impaction efficiency of high inertia particles over the cylinder increases with cross-stream spacing. The non-uniform clustering and local concentration of St = 1 and 10 particles are significantly infl uenced at streamwise positions away from cylinders by cross-stream spacing. Change in the cross-stream spacing between the cylinders modifies magnitudes and cross-stream positions of minimum and peak local particle concentrations.