[PL/CT] Extend residual output for more than one component and to Picard solver

Extensions of residual output for ComponentTransport process

• Output of hydraulic residuum
• Output for more than one component mass/molar balance
• Output, if Picard solver is selected

Tests:

• A test for the combination staggered&Newton solver was already present and is updated.
• A test for the combination staggered & Picard solver has been given a reference by using the old version.
• An existing test for multiple components and monolithic/picard was extendend for the flow-rate tests.
• I have not analytically verified the results. If someone knows well one of the tests for ComponentTransport it would be nice, if you could verify analytically the results and add the corresponding variable to the testing.

1. Feature description was added to the changelog
2. Tests covering your feature were added?
3. Any new feature or behavior change was documented?