[idrees khan]

Dear foamer's


I'm trying to add the Temperature equation(screenshort of the eqaution from the paper is attached) to viscoelastiFluidFoam solver
when i add and compile it with "wmake' I'm getting errors(screenshort is attached).


please note that here solvent viscosity (etaS) polymeric viscosity (etaP) and relaxation time Lambda depinding upon temperature.


humble request kindly guide me.


IWhatsApp Image 2020-12-07 at 5.00.12 PM.jpg

Screenshot from 2020-12-07 10-00-11.jpg

Screenshot from 2020-11-11 12-27-33.jpg

regards
Idrees

[rtyme06]

It's complaining about a mismatched type and coming from viscoelasticFluidFoam.C line 35. Can you share what that line is?

[idrees khan]

Originally Posted by rtyme06

It's complaining about a mismatched type and coming from viscoelasticFluidFoam.C line 35. Can you share what that line is?

Dear raihan

Actually i add temp equation to a file name Teqn.H and add it to viscoelasticFluidFoam.C on line 35

so my equation of temp is





fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi,T)
- fvm::laplacian(DT,T)
- 1.0/rho/Cv * (gamma*visco.tau() && Foam::symm(fvc::grad(U))+ Foam::exp( -C11*(T - To) / (C21 + T - To) )*etaS* Foam::fvc::symm(grad(U))&& Foam::fvc::grad(U))
);

TEqn.relax();
TEqn.solve();


this is what given in the paper (screenshort is attached in previous msg).
there is also some other term which i've to add later.but I'm adding terms step by step.
regards
idrees

[rtyme06]

Based on your error log and the equation, here is my advice for you. Simplify the temperature equation. Create new variables for the double dot product terms and use those variables in the equation. This way you know exactly which term is causing the problem.

[idrees khan]

Originally Posted by rtyme06

Based on your error log and the equation, here is my advice for you. Simplify the temperature equation. Create new variables for the double dot product terms and use those variables in the equation. This way you know exactly which term is causing the problem.

Dear raihan.


fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi,T)
- fvm::laplacian(DT,T)
- 1.0/rho/Cv * (gamma*visco.tau() && Foam::symm(fvc::grad(U))+ (Foam::exp( -C11*(T - To) / (C21 + T - To) )*etaS*Foam::symm(fvc::grad(U)))&& Foam:fvc::grad(U)))
);

TEqn.relax();
TEqn.solve();


when i add and compile with "wmake" getting the same error's


after this i removed the term "&& Foam:fvc::grad(U)))" and compile with wmake it work.
i don't know why
Foam::symm(fvc::grad(U)) is scalar it should be tensor which is basically causing error.
regards
idrees

[rtyme06]

I checked the symm() function. Looks like it returns volSymmTensorField rather than volTensorField. So that may be your problem.
//- Return the symmetric part of a tensor
template<class Cmpt>
inline SymmTensor<Cmpt> symm(const Tensor<Cmpt>& t)
{
return SymmTensor<Cmpt>
(
t.xx(), 0.5*(t.xy() + t.yx()), 0.5*(t.xz() + t.zx()),
t.yy(), 0.5*(t.yz() + t.zy()),
t.zz()
);
}

I guess the possible solution would be to simplify the steps, use additional variables, as I said before.