Unconfined multiple interacting confluent round jets are interesting from a purely scientific point of view, as interaction between neighboring jets brings additional complexity to the flow field. Unconfined confluent round jets also exist in various engineering applications, such as ventilation supply devices, sewage disposal systems, combustion burners, chemical mixing or chimney stacks. Even so, little scientific attention has been paid to unconfined confluent round jets, and detailed investigations on the complete three-dimensional turbulent flow field of confluent jet configurations issuing into an unconfined stagnant environment are scarce.The present work uses both advanced measurement techniques and computational models to provide deeper understanding of the turbulent flow field development of unconfined confluent round jets. Both Laser Doppler Anemometry (LDA) and Particle Image Velocimetry (PIV) have been used to measure mean velocity and turbulence properties within two setups, consisting of a single row of 1×6 jets and a square in-line array of 6×6 confluent jets. LDA was used for measurements very close to the nozzle inlets and to verify the main findings made with PIV. A comparison between LDA and PIV results showed that the two measurement techniques produced results that are in good agreement with each other.Simulations using Computational Fluid Dynamics (CFD) of the 6×6 setup were conducted using three different Reynolds Averaged Navier-Stokes (RANS) turbulence models: the standard k-ε, the RNG k-ε and the Reynolds Stress Model (RSM). The results from the CFD simulations were compared with the experimental PIV data. The employed RANS turbulence models were all capable of accurately predicting mean velocities and turbulent properties in the investigated confluent jet array. In general the RSM and the standard k-ε model provided smaller deviations between numerical and experimental results than the RNG k-ε model. In terms of mean velocity the second-order closure model (RSM) was not found to be superior to the less complex standard k-ε model.The validated CFD results provide data at all locations in the solution domain and the computational results were used to further analyze the investigated confluent jet array. The validated CFD model was also employed in a parametrical investigation, including five independent variables: inlet velocity, nozzle diameter, nozzle edge-toedge spacing, nozzle height and the number of jets in the array. The parametrical investigations made use of statistical methods in the form of Response Surface Methodology and Box-Behnken Design of Experiments. The derived response surface models provided information on the principal influence and relative importance of the investigated parameters within the design limits.ii The positions of the jets within the array strongly influence both mean velocity and turbulence. In all investigated setups the jets experience merging and combining. Square arrays also include considerable jet convergence, which was not present in th...