[qi.yang@polimi.it]

Originally Posted by romant

However, I had a harder time understanding the implementation of the momentum term that in Versteeg and Malalasekera is just given as tau_w * A_cell = S and should be somehow introduced in the momentum equation.

Hi Roman,

Now I need to implement the wall shear stress for solid phase (twophaseeulerfoam) through a source term in the momentum equation. The wall shear stress can be obtained via the following formular:

Sp=Tauw=alpha1*rho1*st*(U1.component(2).boundaryFi eld ()[patchi].patchInternalField())
*(U1.boundaryField()[patchi].patchInternalField() - U1.boundaryField()[patchi])

where "st" is the friction factor which should be calculated ieratively through
for (int i=0; i<10; i++)
{
st = sqr(kappa)/Log(E*Res*sqrt(st));
}
where "Res" is the Reynolds number of solid which calculated by Res=rho1*U1.component(2).boundaryField()[patchi].patchInternalField()*y[facei]/mus

However, now I wrote the codes in the Uequation but cannot be compiled. Anyone can give me some suggestions?