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: 63
Rep Power: 3
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: 24
Rep Power: 3
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: 63
Rep Power: 3
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: 5
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: 63
Rep Power: 3
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: 5
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, 25 views)
Pekka is offline   Reply With Quote

Old   February 25, 2010, 03:45
Default
  #7
Member
 
Join Date: Apr 2009
Location: Karlsruhe, Germany
Posts: 83
Rep Power: 3
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: 5
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

Old   May 28, 2010, 14:28
Default
  #9
New Member
 
Join Date: Feb 2010
Posts: 5
Rep Power: 2
Pekka is on a distinguished road
Hello Foamers,


The problem with wall flux calculation seems to be solved. On incompressible simulation modified wallHeatFlux utility gives “floating point exception” error when kOmegaSST model is used. I change pressure values from zero to 1e-12 on whole calculation field on first time step and then the wallHeatFluxRho utility works.


Can any OF guru tell me, is it possible to solve wall heat flux without any wall functions? What is the reason for changing the wall function to compressible wall function during run the wallHeatFlux utility?


Have I understood correctly that the RAS-model sets the wall functions on the wallHeatFlux utility? For example, if I use the kOmegaSST model the wallHeatFlux utility changes the wall functions to compressible:: omegaWallFunction; and low-Re models zeroGradient; boundary conditions be kept up, is it OK?


Thanks in advance
Pekka is offline   Reply With Quote

Old   June 23, 2010, 09:43
Default
  #10
New Member
 
M.H.Sedaghat
Join Date: Jun 2010
Posts: 1
Rep Power: 0
mh.sedagaht@gmail.com is on a distinguished road
hello
my new username is the next one

Last edited by mh.sedagaht@gmail.com; June 23, 2010 at 16:45.
mh.sedagaht@gmail.com is offline   Reply With Quote

Old   June 23, 2010, 10:37
Default
  #11
New Member
 
Mohammad Hadi Sedaghat
Join Date: Jun 2010
Location: Iran
Posts: 1
Rep Power: 0
M.H.Sedaghat is on a distinguished road
hello
I change the wallHeatFlux files for laminar force and free convection heat transfer solvers(e.g. icoFoam with temperature for laminar force convection and boussinesqBuoyantFoam for laminar natural convection).
It works properly. I rename it wallHeatFluxLaminar. you can download it .

please visit our website:
http://sarreshtehdari.net/
Attached Files
File Type: zip wallHeatFluxLaminar.zip (2.3 KB, 24 views)
M.H.Sedaghat is offline   Reply With Quote

Old   August 12, 2010, 08:55
Default heat flux on multiRegion cases?
  #12
Senior Member
 
maddalena's Avatar
 
maddalena
Join Date: Mar 2009
Posts: 185
Rep Power: 3
maddalena is on a distinguished road
... and what about solvers that require multiple regions, as for example chtMultiRegionFoam? It can be interesting to evaluate the convective heat transfer at the interfaces...
Does someone has an application for this kind of problem as well?
regards,

mad
maddalena 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
problem with sampling Utility in openFOAM 1.6 carmir OpenFOAM Post-Processing 2 August 31, 2010 13:52
wallHeatFlux BC not constant after restart eelcovv OpenFOAM Running / Solving / CFD 20 August 16, 2010 15:00
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


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