WALL FUNCTION: Implement process

boling2002

I am currently working on implementing the wall function to Wilcox wall function.Now I am trying the case for flat plate. My code is using JST scheme. The main steps I am using are: 1) Compute u_{\tau} by using y^{+}=1/k ln( u^{+} + B with given velocity from the first interior cell. 2) Compute \tau_{w} using u_{\tau} and \rho. 3) Using the integral x-component of \tau_{w} multiply the maginitude of surface vector. 4) Set the engergy item to zero since I am working for steady adibatic case.

Is there anything wrong in my implementation, I could not get the resonable results. Could anyone give me your hand ?

Thanks