The transportation of particles in circular pipes holds great significance across industries like mining, oil, gas, chemicals, and food, emphasizing the importance of comprehending two-phase flows. This study aimed to assess the impact of fluid-particle interaction on axial and radial velocity profiles within a laminar boundary layer in a pipe, utilizing COMSOL Multiphysics software with a bidirectional coupling analysis. Initial steps defined parameters influencing velocity profile development. A computational model was constructed, and mesh quality was validated by comparing velocity profiles with existing literature. Radial velocity profiles were examined at 2% and 10% of the theoretical entry length (Le), while axial velocity profiles were studied along the pipe's centerline. Key particle variables, including volume fraction, Stokes number, relative size, and release position, were scrutinized within a range of values. Results highlighted the considerable influence of particle volumetric flow variation on velocity profiles. Particle relative size demonstrated an inversely proportional effect, with larger particles reducing particle concentration and influencing velocity profiles. Increasing the Stokes number rendered particles more independent from the fluid, impacting velocity profiles by the increased particle inertia. Notably, particles' release at the entire pipe inlet caused the most pronounced interferences. In conclusion, this research delved into fluid-particle interaction effects on laminar pipe flow. The study's findings contribute to understanding multiphase systems, impacting various industries and aiding in optimizing particle transport processes.