In the present study, blended models, which combine hybrid RANS/LES models and a correlation-based transition model, have been developed for simulating the flow containing massive flow separation and laminar-turbulence transition inside the boundary layer. To ensure that the base models work well in their own flow regimes of interest, and to achieve smooth transition between the flow regimes, a new blending function was developed based on an empirical analogy from Menter's k-ω SST turbulence model. The new blended models were implemented in a three-dimensional incompressible flow solver based on unstructured meshes. Validations of the blending function and the blended models were conducted for a circular cylinder at a supercritical Reynolds number, an airfoil at high angle of attack, and a dynamic stall problem of an airfoil. The results of the blended models were compared with those of the base models and the experiments. It was found that the new blended models works well by simultaneously capturing the boundary layer transition and the massive flow separation, compared to the base turbulence models. Nomenclature c = chord length of airfoil CD SST = cross diffusion term in the SST model D = diameter of circular cylinder D i = destruction term (i=k, ω, γ, θ t ) d 1 = distance to nearest wall k = turbulent kinetic energy L = turbulent length scale L vK = von Karman length scale P i = production term (i=k, ω, γ, θ t ) R T = eddy viscosity ratio (ρk/( ω μ l )) S = magnitude of strain rate t = unsteady time step U = local velocity U ∞ = freestream velocity x j = Cartesian coordinate in j-direction ρ = density ω = specific turbulent dissipation rate μ l = molecular viscosity μ t = eddy viscosity τ ij = turbulent shear stress γ = local turbulent intermittency γ eff = effective turbulent intermittency Re t = local transition onset momentum thickness Reynolds number Ψ = azimuthal angle measured from windward stagnation point(degree)