The kOmegaSSTLM transition model cannot calculate derivatives during field inversion

[6Mrzhang]

Dear professor He,
The version is dafoam 3.1.2, when I used the kOmegaSSTLM transition model for field inversion, the initial solution converged very well, but when using the runAdjoint task, there was a situation where the derivatives were all 0.

Then I modified the production term Pk (G) in the turbulence kinetic energy equation of the kOmegaSSTLM model

to

Then the full derivative appeared, and it was an intuitive full derivative. Below is a comparison of the derivatives of the two,

This is my runScript. py,

runScript_v2.txt
The only difference between the two cases is the production item, nothing else has changed, I don't know why, and using gammaInt instead of gammaIntEff is not a truly effective transition model. I want to know how to use gammaIntEff to calculate the correct full derivative

[friedenhe]

Is this channel flow fully laminar? If yes, the deriv wrt beta should be zero because it will not impact the flow.