Two-layer turbulent model
 

Hi,
I am trying to add a two-layer turbulent model into OpenFOAM.
The near wall region use k-l model and central region use k-epsilon model.
At first, I was trying to match epsilon, so that I only need to solve k-eqn on time. But I am not sure how to manipulate the wall function of epsilon.
So, now in my code, it solves 3 equations: k_wall-eqn, epsilon-eqn, and k-eqn. And k_wall-eqn and k-eqn is the same.
Could any one give me some advice?
Thanks.