The presence of air bubbles in ship wakes and their dynamic interaction with the turbulent vortical flows create security problems by modifying the hydro-acoustic properties of ship navigation. This study concerns understanding and controlling the mechanisms, which may lead to selective concentration of bubbles to form clusters or clouds and to predict their size distribution and motion. Numerical simulations were conducted using the Large Eddy Simulation (LES) technique in conjunction with a Lagrangian particle tracking (LPT) technique appropriate for dispersed two-phase turbulent flows. In order to cut down on the execution time of LES, the simulation for the flow around the ship model are not considered, instead the simulations started just after the ship where the inlet conditions are prescribed with the help of a newly developed Random Flow Generation (RFG) procedure, and to compute bubble distributions. Moreover, to improve turnaround time of the computations and allow LES of the developing wake in a larger domain, parallel simulation tools are developed and adopted in this study. The turbulence characteristics of ship wakes on a straight tract and a circular tract are investigated using the above mentioned techniques, and processes (e.g. free surface effects, anisotropy, etc.) contributing to turbulence generation are identified, and appropriate sub-grid scale (SGS) models are developed and applied. For the first time, a relatively long developing near wake of three ship lengths was simulated using LES on parallel computers with more than six million nodes.