So I think i have found a solution. First i set up the model with incompressible ideal gas and have all boundary conditions as pressure outlets. Then I initialised with ambient air temperature so make fluent calculate the average fluid density. I then input this average density back into the operating density and run the simulation using this. However, I am not sure why this gives me the correct solution.