A new methodology for finding hybrid RANS/LES damping functions is proposed. The Flow Simulation Methodology framework is modified from Speziale's original formulation [AIAA 36,173 1998] to allow for fewer grid points away from solid boundaries. The proposed functional form is regressed using an evolutionary algorithm. The solution is trained by considering DNS data of a separated flow. This sophisticated approach to building the damping function allows for the development of a hybrid methodology adaptable to any level of required modelling. The regressed model is compared to highly resolved reference LES on the classical two dimensional periodic hills case, for which the agreement is very good.