Numerical determination of wall functions for the turbulent natural convection boundary layer