The aim of this thesis is three-fold: i) to investigate the performance of both the Eulerian-Lagrangian model and the Eulerian-Eulerian model to simulate the turbulent gas-particle flow; ii) to investigate the indoor airflows and contaminant particle flows using the Eulerian-Lagrangian model; iii) to develop and validate particle-wall collision models and a wall roughness model for the Eulerian-Lagrangian model and to utilize these models to investigate the effects of wall roughness on the particle flows. Firstly, the Eulerian-Lagrangian model in the software package FLUENT (FLUENT Inc.) and the Eulerian-Eulerian model in an in-house research code were employed to simulate the gas-particle flows. The validation against the measurement for two-phase flow over backward facing step and in a 90-degree bend revealed that both CFD approaches provide reasonably good prediction for both the gas and particle phases. Then, the Eulerian-Lagrangian model was employed to investigate the indoor airflows and contaminant particle concentration in two geometrically different rooms. For the first room configuration, the performances of three turbulence models for simulating indoor airflow were evaluated and validated against the measured air phase velocity data. All the three turbulence models provided good prediction of the air phase velocity, while the Large Eddy Simulation (LES) model base on the Renormalization Group theory (RNG) provided the best agreement with the measurements. As well, the RNG LES model is able to provide the instantaneous air velocity and turbulence that are required for the evaluation and design of the ventilation system. In the other two-zone ventilated room configuration, contaminant particle concentration decay within the room was simulated and validated against the experimental data using the RNG LES model together with the Lagrangian model. The numerical results revealed that the particle-wall coll ision model has a considerable effect on the particle concentration prediction in the room. This research culminates with the development and implementation of particle-wall collision models and a stochastic wall roughness model in the Eulerian-Lagrangian model. This Eulerian-Lagrangian model was therefore used to simulate the gas-particle flow over an in-line tube bank. The numerical predictions showed that the wall roughness has a considerable effect by altering the rebounding behaviours of the large particles and consequently affecting the particles motion downstream along the in-line tube bank and particle impact frequency on the tubes. Also, the results demonstrated that for the large particles the particle phase velocity fluctuations are not influenced by the gas-phase fluctuations, but are predominantly determined by the particle-wall collision. For small particles, the influence of particle-wall collisions on the particle fluctuations can be neglected. Then, the effects of wall roughness on the gas-particle flow in a two-dimensional 90-degree bend were investigated. It was found that the wa ll roughness considerably altered the rebounding behaviours of particles by significantly reducing the 'particle free zone' and smoothing the particle number density profiles. The particle mean velocities were reduced and the particle fluctuating velocities were increased when taking into consideration the wall roughness, since the wall roughness produced greater randomness in the particle rebound velocities and trajectories.