A novel two-equations model for computing the flow properties of a spatially-developing, incompressible, zero-pressure-gradient, turbulent boundary layer over a smooth, flat wall is developed. The mean streamwise velocity component inside the boundary layer is described by the Reynolds-averaged Navier–Stokes equation where the Reynolds shear stress is given by an extended mixing-length model. The nondimensional form of the mixing length is described by a polynomial function in terms of the nondimensional wall normal coordinate. Moreover, a stream function approach is applied with a leading-order term described by a similarity function. Two ordinary differential equations are derived for the solution of the similarity function along the wall normal coordinate and for its streamwise location. A numerical integration scheme of the model equations is developed and enables the solution of flow properties. The coefficients of the mixing-length polynomial function are modified at each streamwise location as part of solution iterations to satisfy the wall and far-field boundary conditions and adjust the local boundary layer thickness, , to a location where streamwise speed is 99.4% of the far-field streamwise velocity. The elegance of the present approach is established through the successful solution of the various flow properties across the boundary layer (i.e., mean streamwise velocity, viscous stress, Reynolds shear stress, skin friction coefficient, and growth rate of boundary layer among others) from the laminar regime all the way to the fully turbulent regime. It is found that results agree with much available experimental data and direct numerical simulations for a wide range of based on the momentum thickness () from up to , except for the transition region from laminar to turbulent flow. Furthermore, results shed light on the von Kármán constant as a function of , the possible four-layer nature of the mean streamwise velocity profile, the universal profiles of the streamwise velocity and the Reynolds shear stress at high , and the scaling laws at the outer region.