[kiamey]

Hi
I have developed a code which uses adi method to simulate shallow water equations. For turbulence modelling I have used a basic mixing length concept and it works.
This turbulence model is a zero-equation model, so no additional equations shall be solved and I have no problem in solving the tridiagonal systems of equations.
But k-epsilon model adds two other unknowns.
I have no idea how to add these equations to my model. How can I descritize K-epsilon equations? How can I solve the new system of equations which includes 5 unknowns (water depth, velocities in two x and y directions,k and epsilon).
Does someone have any idea or can somebody introduce a reference on the subject.
Thank you