adambarfi |
May 10, 2013 01:58 |
problem in adding dissipation term to T equation!!!
Hi everybody,
I want to add the dissipation term to my T equation. I did it as follows:
Code:
dimensionedScalar kappa=(etaS+etaP)/(rho*Pr);
fvScalarMatrix TEqn
(
fvm::div(phi, T)
- fvm::Sp(fvc::div(phi), T)
- fvm::laplacian(kappa, T)
- (visco.tau() && (fvc::curl(U))) //Dissipation Term
);
where visco is my viscoelastic model and visco.tau is a symmetric tensor.
when I wmake my solver the following errros appeared:
Code:
In file included from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/Field.C:735:0,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/Field.H:360,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/labelField.H:39,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/primitiveFields.H:37,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/pointField.H:36,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/edge.H:40,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/edgeList.H:32,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/primitiveMesh.H:57,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/polyMesh.H:44,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/finiteVolume/lnInclude/fvMesh.H:50,
from /home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/finiteVolume/lnInclude/fvCFD.H:7,
from bbvffPerfectTdepDiss2.C:36:
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/FieldFunctions.C: In function ‘void Foam::dotdot(Foam::Field<typename Foam::scalarProduct<Type1, Type2>::type>&, const Foam::UList<T>&, const Foam::UList<Key>&) [with Type1 = Foam::SymmTensor<double>, Type2 = Foam::Vector<double>, typename Foam::scalarProduct<Type1, Type2>::type = double]’:
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/GeometricFieldFunctions.C:956:1: instantiated from ‘void Foam::dotdot(Foam::GeometricField<typename Foam::scalarProduct<Type1, Type2>::type, PatchField, GeoMesh>&, const Foam::GeometricField<TypeR, PatchField, GeoMesh>&, const Foam::GeometricField<Type1, PatchField, GeoMesh>&) [with Type1 = Foam::SymmTensor<double>, Type2 = Foam::Vector<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh, typename Foam::scalarProduct<Type1, Type2>::type = double]’
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/GeometricFieldFunctions.C:956:1: instantiated from ‘Foam::tmp<Foam::GeometricField<typename Foam::scalarProduct<Type1, Type2>::type, PatchField, GeoMesh> > Foam::operator&&(const Foam::tmp<Foam::GeometricField<TypeR, PatchField, GeoMesh> >&, const Foam::tmp<Foam::GeometricField<Type1, PatchField, GeoMesh> >&) [with Type1 = Foam::SymmTensor<double>, Type2 = Foam::Vector<double>, PatchField = Foam::fvPatchField, GeoMesh = Foam::volMesh, typename Foam::scalarProduct<Type1, Type2>::type = double]’
TEqn.H:22:38: instantiated from here
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/FieldFunctions.C:705:1: error: no match for ‘operator&&’ in ‘*(f2P ++) && *(f3P ++)’
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/FieldFunctions.C:705:1: note: candidates are: operator&&(bool, bool) <built-in>
/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/src/OpenFOAM/lnInclude/dimensionSet.H:293:29: note: Foam::dimensionSet Foam::operator&&(const Foam::dimensionSet&, const Foam::dimensionSet&)
make: *** [Make/linuxGccDPOpt/bbvffPerfectTdepDiss2.o] Error 1
anybody knows where is the problem? I'm completely confused.
thanks
Mostafa
|