Hi,
1) It depends on what you have set as your rhoS. If you exactly know the rhoS you want to use for injection, you can simply add it to the well as .rhoS and the simulator will use that value. For most users, flashing at the surface to get injection densities
will be a reasonable default and this is in line with what most commercial simulators do.
2) Specifying rhoS from flash will give the same result. With the old models you are limited to a single density for all wells, though, so not all cases are possible to simulate with the old model.
3) It is irrelevant if you do not have an aqueous phase just as the initial saturations in the reservoir are irrelevant for the same setup. If a third phase is present, the sum of the two phases governed by the EoS matter.
4) If a well is producing then the surface density in the fluid is used (if no separator is specified). Components are then associated with a phase stream based on their LV-partition as a pure mixture at surface conditions, which is a simplification to avoid
re-flashing (but a separator can be used to get a more realistic behavior).
5) The generic models do not output components, but instead output the more general total masses per second in the wellSol which can be used to reconstruct the molar fractions. W.components is still used to set the injection composition.