trublent flow with solidification

Hello all..

I am solving a turblent flow coupled with heat transfer and solidification.

In the model library to model the same problem for continuous casting,, laminar flow is used.
I am including the turblence equations..

while solving CC in model libarary , to dampen the velocity in ths solid region source term F is used.

now in turblence equation to dampen the turblent KE generation and KE dissipation in the solid region..
i want to use similar kind of source term in the governing equations of k and e.

rho*u*div k = ...........+.......+....+....+Fk
rho*u*div e = ...........+.......+....+....+Fe

where Fk and Fe are the terms to add in governing equations of k and e..
how can i d o it..i am trying to use weak contribution but not succesful..

Fk and Fe are similar to Fr an dFz in CC prsent in model library
Fr and Fz are added in navier stokes ..
i need Fk and Fe to add in k and e equations..

can anyone suggest how to do that

i know that i am solving a highly nonliner problem and my system is crying to get the solution..

please suggest the better solver steps to get the solution faster with desired accuracy

thakns in advance