[maruthamuthu_venkatraman]

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

[ulli]

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

[maruthamuthu_venkatraman]

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.

[Pekka]

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

[maruthamuthu_venkatraman]

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
-----------------------------------------------------------------------

[Pekka]

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?

[Thomas Baumann]

Hi,

have a look at the discussion at
HeatFlux in buoyantBoussinesq

Regards Thomas

[Pekka]

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]

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

[mh.sedagaht@gmail.com]

hello
my new username is the next one

[M.H.Sedaghat]

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/

[maddalena]

... 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