Abstract: The derivation of Gao-Yong turbulence equations was based on a partial statistical average scheme proposed in the past by the first author in an attempt to capture the drift flow, a first-order statistical moment of turbulent fluctuations. Continuity and momentum equations of turbulent fluctuations were obtained as a result of this valuable first statistical information. Orthotropic turbulence and momentum transfer chain were formulated in the modeling of correlation terms, eventually, led to a complete set of equations of turbulent flow with no empirical coefficients and wall functions. Numerical simulations in the past preliminarily verified the adaptability of Gao-Yong equations to a wide range of complex benchmark turbulent flows. In order to further application of these equations to engineering flows, the Gao-Yong turbulence equations were implemented into OpenFOAM and applied to the calculation of wing-body junction flow. The wing-body junction flow involves very complex feathers, such as a 3-D separation and a horseshoe vortex. The results of numerical simulation of this flow compared to experimental data show that all these main characteristics have been captured, both qualitatively and quantitatively.