Dear OpenFOAM users. I am u
 

Dear OpenFOAM users.

I am using OpenFOam-1.5 I wanted to solve the transport equation for temperature variance for my research purpose. I have added Temperature equation and Trms equation in the simplFoam code for the same purpose.( I have replaced buoyantSimpleRadiation code with SimpleFOam code sincemy original simpleFoam doce is running with some other modifications in energy equation)

My TEqn.H file reads as

aplhaEff = aplha + turbulence->nut()/Prt;

tmp<fvscalarmatrix>TEqn
(
fvm::div(phi, T)
- fvm::laplacian(aphaEff, T)
);
TEqn().relax();
eqnResidual() = solve(TEqn()).intialResidual();
maxResidual = max(eqnResidual, maxResidual);


My TrmsEqn.H file reads as

alphavariance = alpha + (0.13*(turbulence->k()*turbulence->k()/turbulence->epsilon()));

tmp<fvscalarmatrix>
(
fvm::div(phi, Trms)
-fvm::laplacian(alphavariance, Trms)
+fvm::Sp(fvc::div(aplhaEff, T), Trms)
+0.134*(turbulence->epsilon()/turbulence->k())*Trms
);
TrmsEqn().relax();
eqnResidual = solve(TrmsEqn()).initialResidual();
maxResidual = max(eqnResidual, maxResidual);

On compilation i get following error

Ameya@Ameya:~/OpenFOAM/OpenFOAM-1.5/applications/solvers/heatTransfer/buoyantSimpleRadiatio nFoam> wclean
Ameya@Ameya:~/OpenFOAM/OpenFOAM-1.5/applications/solvers/heatTransfer/buoyantSimpleRadiatio nFoam> wmake
Making dependency list for source file buoyantSimpleRadiationFoam.C
SOURCE=buoyantSimpleRadiationFoam.C ; g++ -m32 -Dlinux -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-40 -I/home/Ameya/OpenFOAM/OpenFOAM-1.5/src/finiteVolume/lnInclude -I/home/Ameya/OpenFOAM/OpenFOAM-1.5/src/turbulenceModels/RAS -I/home/Ameya/OpenFOAM/OpenFOAM-1.5/src/transportModels -IlnInclude -I. -I/home/Ameya/OpenFOAM/OpenFOAM-1.5/src/OpenFOAM/lnInclude -I/home/Ameya/OpenFOAM/OpenFOAM-1.5/src/OSspecific/Unix/lnInclude -fPIC -pthread -c $SOURCE -o Make/linuxGccDPOpt/buoyantSimpleRadiationFoam.o
In file included from buoyantSimpleRadiationFoam.C:72:
TrmsEqn.H: In function 'int main(int, char**)':
TrmsEqn.H:8: error: no matching function for call to 'div(Foam::volScalarField&, Foam::volScalarField&)'
make: *** [Make/linuxGccDPOpt/buoyantSimpleRadiationFoam.o] Error 1



Any ideas on what the error actually means????
I know there is some mistake in the way in which the source term of Trms equation is written but i am not able to find out what and exactly how that term should be written. Can anyone please help me?????????

Hello Ameya, I have a feeling
 

Hello Ameya,
I have a feeling that your nomenclature is a little mixed up: one time you write "alpha" and the other "aplha", the same for alphaEff.

Do the error-messages disappear when correcting that? Otherwise after a first overview I see no other problems in your code.

Hey, TrmsEqn.H:8: error: no
 

Hey,

TrmsEqn.H:8: error: no matching function for call to 'div(Foam::volScalarField&, Foam::volScalarField&)'

it refers to:
+fvm::Sp(fvc::div(aplhaEff, T), Trms)

what is the mathematical expression you want to express, the divergence of a vector I presume?
also fvc::div generally needs a surfaceScalarField to operate.

see for example:
fvScalarMatrix bEqn
(
fvm::ddt(rho, b)
+ mvConvection->fvmDiv(phi, b)
+ fvm::div(phiSt, b, "div(phiSt,b)")
- fvm::Sp(fvc::div(phiSt), b)
- fvm::laplacian(turbulence->muEff(), b)
);

hope this helps

Pierre

I forgot, try to interpolate a
 

I forgot, try to interpolate aplhaEff to the faces of your cell.

Pierre