Well, it is true that turbulence modelling is needed to solve the closure problem when applying Reynolds averaging (space-time filtering in general) to the NSE, but Reynolds Averaging is not a condition for the numerical solution of the NSE. These two are separate concepts. You can run any simulation using PISO/SIMPLE without a turbulence model. In FOAM, that means the turbulence fluxes and viscosity are identically zero.