I'm simulating a liquid metal flow through a rod bundle. I first run a simulation with the kOmegaSST model and all went good. Then, since I'm interested in resolving secondary flows and streamline curvature, I set up a run with the SSG model. I calculated R and epsilon at the end of the kOmega run using the "turbulenceFields" functionObject (I'm using OF v9), and copied them into the 0 folder of the SSG run. Here come the problems: in my 1st attempt I changed the BCs at the wall for both R and epsilon into "fixedValue" (with very small numbers, since R must approach zero at the wall) and "epsilonWallFunction", respectively. The simulation converged well but the wall shear stress was overestimated (about the double with respect to the kOmega run and to empirical correlations). In a 2nd attempt, I did not change the wall boundary condition for epsilon, leaving the one that was present in the original epsilon file calculated at the end of the kOmega run (I attach a screenshot for clarity). In this case, the simulation did not converge, velocities kept oscillating, but the wall shear stress magnitude is closer the expected one (although itself oscillating, of course). I also tried to simulate a pseudo transient to help convergence, but I did not reach it. Am I missing something in the procedure to run the SSG model, particularly in the BCs? Thank you in advance.