CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Transport equation for Temperatue variace

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   December 5, 2008, 02:06
Default Dear OpenFOAM users. I am u
  #1
New Member
 
Ameya Durve
Join Date: Mar 2009
Location: Mumbai, Maharashtra, India
Posts: 20
Rep Power: 8
ameya is on a distinguished road
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?????????
ameya is offline   Reply With Quote

Old   December 15, 2008, 09:57
Default Hello Ameya, I have a feeling
  #2
New Member
 
Tammo Wenterodt
Join Date: Mar 2009
Posts: 24
Rep Power: 8
wenterodt is on a distinguished road
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.
wenterodt is offline   Reply With Quote

Old   December 15, 2008, 12:55
Default Hey, TrmsEqn.H:8: error: no
  #3
Member
 
Pierre Le Fur
Join Date: Mar 2009
Location: UK
Posts: 60
Rep Power: 8
pierre is on a distinguished road
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
pierre is offline   Reply With Quote

Old   December 15, 2008, 12:58
Default I forgot, try to interpolate a
  #4
Member
 
Pierre Le Fur
Join Date: Mar 2009
Location: UK
Posts: 60
Rep Power: 8
pierre is on a distinguished road
I forgot, try to interpolate aplhaEff to the faces of your cell.

Pierre
pierre is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
UDF for concentration transport equation heri FLUENT 1 August 15, 2011 11:23
Poisson equation combined with transport equation Se-Hee CFX 0 December 27, 2007 02:00
transport equation Arman Main CFD Forum 1 August 28, 2007 12:48
Transport equation using UDF wassim FLUENT 5 February 14, 2007 23:56
transport equation CMB CD-adapco 2 July 9, 2004 05:19


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