I've read a bunch of articles claiming that wall functions, commonly coupled with 1st order turbulence models, are not suitable for free convection in thermo aeraulic. They are most likely suitable for forced convection situations.
I wonder if it would be a good approach to couple integral laws (semi empiric correlations adapted for free convection) to a CFD code in order to estimate hcv near the walls?