WIP: ENH: renovate v2-f turbulence closure model

corrects strain rate and tke production term to enable consistent compressible formulation

corrects Ts and Ls function expressions according to the original paper

corrects the transport equation terms to enable consistent compressible formulation

adds protections against floating-point exceptions

Usage-related:

• updates the header documentation
• allows the hard-coded constants to be tuned by users
• adds tags to easily track the implementations' expressions in the original papers