Pneumatic conveying systems are employed across different industries with its routing flexibility being one of the reasons it is favoured over other systems. In pneumatic conveying, material fed into the system at the inlet goes through a series of closed piping in a very clean operation. However, the use of pipe bends for routing can cause a phenomenon called particle roping which may decrease the overall efficiency of the system. Due to inertial forces acting on the conveyed Greek alphabets Volume fraction Θ Collisional dissipation of energy Θ Granular temperature Bulk viscosity / Second viscosity Dynamic viscosity Eddy viscosity Density Turbulent Prandtl number ̿ Stress tensor Subscripts Gas phase Arbitrary phase Solid phase [4] on the pickup velocity of particles and de Moraes et al. [5] on head loss coefficient for pipe fittings, both of which deal with properties fundamental to pneumatic conveying processes. While experimental studies form a significant portion of pneumatic conveying investigations, the number of numerical studies are increasing as well. There are a few numerical approaches to study pneumatic conveying flows that build on top of established methods in calculating single-Constitutive relations are usually derived for certain flow conditions and may not be adequate beyond the range that it was originally derived for [10]. Therefore, this thesis first aims to demonstrate the ability of an Eulerian scheme in predicting the flow characteristics of a granular phase in a pneumatic pipe configuration Particle ropes, a region of the pipe cross-sectional area where there is an increased concentration of particles, are caused by pipe bends [14]. When particles enter a pipe bend, inertial effects will tend to cluster the particles into a region of higher particle concentration. This area is usually located on the outer wall of the pipe bend. Particle ropes are sometimes seen as an unwanted occurrence in operations,