CFD Online Logo CFD Online URL
www.cfd-online.com
Home > Forums > OpenFOAM

wallHeatFlux utility in OpenFoam1.6

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 9, 2010, 09:39
Default wallHeatFlux utility in OpenFoam1.6
  #1
Member
 
Maruthamuthu Venkatraman
Join Date: Mar 2009
Location: Norway
Posts: 50
Rep Power: 2
maruthamuthu_venkatraman is on a distinguished road
Hello Foamers,
In order to check the heat flux balance i tried to use the utility wallHeatFlux. It seems it only works for Combustion solvers, when i run this utility for all those tutorial problems in heatTransfer it says as below:
-----------------------------------------------------------------------

Create time
Create mesh for time = 0
Time = 0
Selecting thermodynamics package hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>

Unknown hCombustionThermo type hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Valid hCombustionThermo types are:
9
(
hPsiMixtureThermo<homogeneousMixture<sutherlandTra nsport<specieThermo<janafThermo<perfectGas>>>>>
hPsiMixtureThermo<reactingMixture<gasThermoPhysics >>
hPsiMixtureThermo<inhomogeneousMixture<constTransp ort<specieThermo<hConstThermo<perfectGas>>>>>
hPsiMixtureThermo<multiComponentMixture<gasThermoP hysics>>
hPsiMixtureThermo<veryInhomogeneousMixture<constTr ansport<specieThermo<hConstThermo<perfectGas>>>>>
hPsiMixtureThermo<homogeneousMixture<constTranspor t<specieThermo<hConstThermo<perfectGas>>>>>
hPsiMixtureThermo<inhomogeneousMixture<sutherlandT ransport<specieThermo<janafThermo<perfectGas>>>>>
hPsiMixtureThermo<veryInhomogeneousMixture<sutherl andTransport<specieThermo<janafThermo<perfectGas>> >>>
hPsiMixtureThermo<dieselMixture<sutherlandTranspor t<specieThermo<janafThermo<perfectGas>>>>>
)

From function hCombustionThermo::New(const fvMesh&)
in file combustionThermo/hCombustionThermo/newhCombustionThermo.C at line 66.
FOAM exiting
-----------------------------------------------------------------------

Then i tried to look at the wallHeatFlux utility and found "#include basicRhoThermo.H " is missing. I added to it and compiled the same. Again this doesnt solve the problem.

Any hints are welcome.

Thanks
maruthamuthu_venkatraman is offline   Reply With Quote

Old   February 10, 2010, 09:17
Default
  #2
New Member
 
Ulrich Heck
Join Date: Mar 2009
Location: Krefeld, Germany
Posts: 19
Rep Power: 2
ulli is on a distinguished road
Hello Maruthamuthu,
I did this like this
1. I copied the wallHeatFlux-dir into a new dir e.g. wallHeatFluxRho
2. I renamed the *c-File in dir wallHeatFluxRho to wallHeatFluxRho.c
3. I modified the name in "files" (Make –dir)
4. I replaced (with an editor) hCombustionThermo by basicPsiThermo (which is the base class, see Doxygen) in createFields.h and wallHeatFluxRho.c. I think this class has be renamed in 1.6
5. Than I complied with wclean/wmake
This should work.
Hope it helps
Ulrich
ulli is offline   Reply With Quote

Old   February 16, 2010, 13:11
Default
  #3
Member
 
Maruthamuthu Venkatraman
Join Date: Mar 2009
Location: Norway
Posts: 50
Rep Power: 2
maruthamuthu_venkatraman is on a distinguished road
Thanks for your explanations.

Yes, it worked. I can able to run this utility for simpleRadiationFoam. I compiled BuoyantPisoRadiationFoam and tried to use wallHeatFluxRho utility. It doesnt work for that case. I will check it out in detail later.
maruthamuthu_venkatraman is offline   Reply With Quote

Old   February 20, 2010, 14:16
Default Problem with incompressible wall heat flux calculation
  #4
New Member
 
Join Date: Feb 2010
Posts: 4
Rep Power: 2
Pekka is on a distinguished road
Hello Sir/Madam

I carried out benchmark calculations by using OF 1.6.x for cylinder on the cross flow and one interesting point is on the wall heat convection. Based on previous tips I do “wallHeatFluxRho” function so that I change “hCombustionThermo” to “basicRhoThermo” on wallHeatFluxRho.C and createFields.H files and complied with wclean/wmake. WallHeatFluxRho works on compressible e.g. BuoyantPisoFoam many different RASmodels but it works only realizableKE and LaunderSharmaKE models on the incompressible flow simulation (buoyantBoussinesqPisoFoam) . I use on the thermo physical models

“thermoType hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>;”.

How can I calculate wall heat flux on the incompressible case? WallHeatFluxRho function does not work with for example kOmegaSST model.

Thanks in advance
Pekka is offline   Reply With Quote

Old   February 23, 2010, 05:33
Default
  #5
Member
 
Maruthamuthu Venkatraman
Join Date: Mar 2009
Location: Norway
Posts: 50
Rep Power: 2
maruthamuthu_venkatraman is on a distinguished road
Hi Pekka,
I couldn't able to run thios utility wallHeatFluxRho for buoyantPisoFoam. The following error is showing up. Since you run this utility for buoyantPisoFoam , could you upload the sources. I will recompile and see my mistake.

Thanks.

-------------------------------Error Message---------------------------

Create time
Create mesh for time = 0
Time = 0
Selecting thermodynamics package hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>

Unknown basicPsiThermo type hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Valid basicPsiThermo types are:
25
(
ePsiThermo<pureMixture<sutherlandTransport<specieT hermo<janafThermo<perfectGas>>>>>
hhuMixtureThermo<veryInhomogeneousMixture<sutherla ndTransport<specieThermo<janafThermo<perfectGas>>> >>
hhuMixtureThermo<egrMixture<sutherlandTransport<sp ecieThermo<janafThermo<perfectGas>>>>>
hhuMixtureThermo<homogeneousMixture<constTransport <specieThermo<hConstThermo<perfectGas>>>>>
hPsiMixtureThermo<homogeneousMixture<sutherlandTra nsport<specieThermo<janafThermo<perfectGas>>>>>
ePsiThermo<pureMixture<sutherlandTransport<specieT hermo<eConstThermo<perfectGas>>>>>
hhuMixtureThermo<veryInhomogeneousMixture<constTra nsport<specieThermo<hConstThermo<perfectGas>>>>>
hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
hPsiMixtureThermo<reactingMixture<gasThermoPhysics >>
hhuMixtureThermo<inhomogeneousMixture<sutherlandTr ansport<specieThermo<janafThermo<perfectGas>>>>>
hPsiMixtureThermo<inhomogeneousMixture<constTransp ort<specieThermo<hConstThermo<perfectGas>>>>>
hPsiMixtureThermo<multiComponentMixture<gasThermoP hysics>>
hPsiMixtureThermo<veryInhomogeneousMixture<constTr ansport<specieThermo<hConstThermo<perfectGas>>>>>
hhuMixtureThermo<egrMixture<constTransport<specieT hermo<hConstThermo<perfectGas>>>>>
hPsiMixtureThermo<homogeneousMixture<constTranspor t<specieThermo<hConstThermo<perfectGas>>>>>
hPsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>
hPsiMixtureThermo<inhomogeneousMixture<sutherlandT ransport<specieThermo<janafThermo<perfectGas>>>>>
hPsiMixtureThermo<veryInhomogeneousMixture<sutherl andTransport<specieThermo<janafThermo<perfectGas>> >>>
hPsiMixtureThermo<dieselMixture<sutherlandTranspor t<specieThermo<janafThermo<perfectGas>>>>>
ePsiThermo<pureMixture<constTransport<specieThermo <eConstThermo<perfectGas>>>>>
hhuMixtureThermo<homogeneousMixture<sutherlandTran sport<specieThermo<janafThermo<perfectGas>>>>>
ePsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>
ePsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
hPsiThermo<pureMixture<sutherlandTransport<specieT hermo<janafThermo<perfectGas>>>>>
hhuMixtureThermo<inhomogeneousMixture<constTranspo rt<specieThermo<hConstThermo<perfectGas>>>>>
)

From function basicPsiThermo::New(const fvMesh&)
in file psiThermo/basicPsiThermo/newBasicPsiThermo.C at line 64.
FOAM exiting
-----------------------------------------------------------------------
maruthamuthu_venkatraman is offline   Reply With Quote

Old   February 23, 2010, 11:53
Default
  #6
New Member
 
Join Date: Feb 2010
Posts: 4
Rep Power: 2
Pekka is on a distinguished road
Hi Maruthamuthu,


take a copy of wallHeatFlux folder and rename it, replace wallHeatFlux.C, createFiles.H and files to the new one and then compile. It should work with the BuoyantPisoFoam solver.


Has anybody any ideas how wall heat flux calculations are done for incompressible flow?
Attached Files
File Type: gz wallHeatFluxRho.tar.gz (1.6 KB, 4 views)
Pekka is offline   Reply With Quote

Old   February 25, 2010, 03:45
Default
  #7
Member
 
Join Date: Apr 2009
Location: Karlsruhe, Germany
Posts: 59
Rep Power: 2
Thomas Baumann is on a distinguished road
Hi,

have a look at the discussion at
HeatFlux in buoyantBoussinesq

Regards Thomas
Thomas Baumann is offline   Reply With Quote

Old   February 27, 2010, 14:48
Default
  #8
New Member
 
Join Date: Feb 2010
Posts: 4
Rep Power: 2
Pekka is on a distinguished road
Hi,


thanks for the link and it clarified the problematics, but I have still the same problem with incompressible flow calculation. I change all kappaEff to alphaEff in the Tegn.H and recompiled it, but without success.

When I try post processing heat flux on the walls with wallHeatFluxRho function I get the following error:

----------Error message-------

Time = 0
Selecting thermodynamics package hRhoThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>
Reading/calculating face flux field phi

Selecting RAS turbulence model kOmegaSST
--> Upgrading k to employ run-time selectable wall functions
Backup original k to k.old
Writing updated k
--> Upgrading omega to employ run-time selectable wall functions
Backup original omega to omega.old
Writing updated omega
--> Creating mut to employ run-time selectable wall functions
Writing new mut
--> Creating alphat to employ run-time selectable wall functions
Writing new alphat
#0 Foam::error::rintStack(Foam::Ostream&) in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:erator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#6 Foam::compressible::RASModels::kOmegaSST::F2() const in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#7 Foam::compressible::RASModels::kOmegaSST::kOmegaSS T(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#8 Foam::compressible::RASModel::adddictionaryConstru ctorToTable<Foam::compressible::RASModels::kOmegaS ST>::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#9 Foam::compressible::RASModel::New(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::basicThermo const&) in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#10 main in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/wallHeatFluxRho"
#11 __libc_start_main in "/lib64/libc.so.6"
#12 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/Pekka/OpenFOAM/OpenFOAM-1.6.x/applications/bin/linux64GccDPOpt/wallHeatFluxRho"
Floating point exception


Any ideas?


BR Pekka
Pekka 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
Sample Utility not working in OpenFoam 1.6 titio OpenFOAM 4 March 5, 2010 10:55
Sample Utility not working in OpenFoam 1.6 titio OpenFOAM Post-Processing 0 February 5, 2010 12:12
How to compile a new utility rudy OpenFOAM 3 December 14, 2009 08:43
problem with sampling Utility in openFOAM 1.6 carmir OpenFOAM Post-Processing 1 September 3, 2009 06:07
wallHeatFlux BC not constant after restart eelcovv OpenFOAM Running / Solving / CFD 6 July 31, 2009 09:04


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