[chnrdu]

I am interested in the viscoelastic flow with free surface.

Thanks Favero, who gives me a way to solve it.

But I have tried the numerical algorithm, and have not get the right solution.

I combined the interFoam and viscoelasticFluidFoam, and using the Giesekus constitutive model for viscoelastic flow.

The attachment is my combined source. I think the critical part is the constitutive correction of tauP ( in the source, it is tau_), and divTau(U)

the tauP correction is
fvSymmTensorMatrix tauEqn
(
etaP()*lambda()*fvm::ddt(tau_)
+ etaP()*lambda()*fvm::div(phi(), tau_)
==
sqr(etaP())*twoD
+ etaP()*lambda()*twoSymm(C)
- lambda()*(alpha())*(tau_ & tau_)
- etaP()*fvm::Sp(scalar(1), tau_)
)

where the etaP, lambda and alpha is respectively the combined two phase value with respectively gamma.

the divTau(U) is

fvc::div(tau_)
- fvm::laplacian(etaPf(), U)
+ fvm::laplacian( (etaPf()+etaSf()), U)

Please help me to review the source.

Thank you very much!

Martin