Wall functions are a set of semi-empirical formulas used in computational fluid dynamics (CFD) to approximate the behavior of turbulent flows near solid boundaries, where the flow is significantly affected by viscosity and turbulent shear stresses. They bridge the gap between the fully turbulent flow and the viscous sublayer, allowing for more efficient and less computationally intensive simulations by avoiding the need to resolve the entire boundary layer in detail.