Implementation of wall function with SA model
 

Hi, I'm trying to implement the wall function in a 3d RANS compressible c.c.f.v solver. Without adiabatic density correction the implementation is working fine. However the wall function for the temperature T+ implementation is a bit confusing. Can anyone suggest me the theory and the methodology. I followed the NT Frink's paper "Assessment of unstructured grid ..".