CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   problem in adding dissipation term to T equation!!! (http://www.cfd-online.com/Forums/openfoam-programming-development/117527-problem-adding-dissipation-term-t-equation.html)

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

pashanxiao June 25, 2015 02:14

add the dissipation term to TEqn
 
Hi Mostafa,
I have the same problem. I add the dissipation term[ - (visco.tau() && (fvc::grad(U)))] to TEqn, when I wmake the solver , the following errros appeared just as your problem. Did your solution solve the problem?

I am a newbie in this field, could you give me tips how could I implement it.

Thanks in advance.

pashanxiao

adambarfi June 25, 2015 03:05

Quote:

Originally Posted by pashanxiao (Post 552055)
Hi Mostafa,
I have the same problem. I add the dissipation term[ - (visco.tau() && (fvc::grad(U)))] to TEqn, when I wmake the solver , the following errros appeared just as your problem. Did your solution solve the problem?

I am a newbie in this field, could you give me tips how could I implement it.

Thanks in advance.

pashanxiao

Greetings Qingwu,

It's for 2 years ago and I can't remember what had I done:confused::confused::confused:

Perhaps I've expanded the dissipation term and added it term-by-term.

Regards,
Mostafa


All times are GMT -4. The time now is 09:09.