[Tobi]

hmm ... don't know.

but i think that's a problem, couse you can not decline some thermophysicalProperties for a solid. Just for fluids you can declare it.

good question.
I 'm trying on a case today afternoon (chtMulti) with one fluid and solid region. Hope i can calculate the wallFlux with the post tool.

i ll give an replay about it.
see you later.

Tobi

[Tobi]

Question:

how did you successful start wallHeatFlux on the fluid region?

you put every file into another folder or? And generate a "virtual" case for you post tool?


-> i did a simulation yesterday with one fluid and solid region. Solved it an put every fluid file in another folder. After using wallHeatFluxLaminar i got those message:

--> FOAM FATAL IO ERROR: 
Unknown patchField type compressible::turbulentTemperatureCoupledBaffleMixed for patch type directMappedWallValid patchField types are :
42
(
advective
buoyantPressure
calculated
cyclic
directMapped
directionMixed

then i changed the patchField type to directMapped with the tool changeDictionary.

After that i got the message

--> FOAM FATAL ERROR: 
compound has already been transfered from token
    on line 20 the empty compound of type List<scalar>

    From function token::transferCompoundToken()
    in file lnInclude/token.C at line 95.

FOAM aborting

How did you get your wallHeatFlux in you fluid region?
Did you use another wallHeatFlux - post Tool?

regards Tobi

[NicolasB]

I've tried the wallHeatFluxRho utility from a previous post in this thread, and I've copied all files from fluid zone to a new directory (0, 5000, constant, system) in which I use the utility.

Hope this helps you.

[Tobi]

sure, ... i 've tried it with that tool, too.

What 's your T - File ?
can you post it ?

Thx Tobi

[NicolasB]

Hi Tobi,
Since I'm not at the office, I haven't got the needed file. However, I've posted my case here (if you have any hint on this subject, I'd be glad to have your advice )

Regards

[Tobi]

Hey Nicolas,

short question. Did you get the wallHeatFlux calculated at the surfaces:

" fluid_air_to_solid_cable " and " fluid_air_to_solid_beton ".


Tobi

[NicolasB]

Here is what I get in the console:

caelinux@caelinux-desktop:~/OpenFOAM/caelinux-1.7.0/Freyssinet3/air$ wallHeatFluxRho 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.7.0-113391ee57bd
Exec   : wallHeatFluxRho
Date   : Apr 27 2011
Time   : 13:11:10
Host   : caelinux-desktop
PID    : 25175
Case   : /home/caelinux/OpenFOAM/caelinux-1.7.0/Freyssinet3/air
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

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

Selecting RAS turbulence model laminar

Wall heat fluxes [W]
fluid_air_to_solid_cable 0
fluid_air_to_solid_beton 0

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

Selecting RAS turbulence model laminar

Wall heat fluxes [W]
fluid_air_to_solid_cable -981.77797
fluid_air_to_solid_beton -1298.0191

End

NB: I work with OF 1.7.0

Nicolas

[ishaninair]

Hi Ulrich Heck,
I did the same thing as told by you. But I am getting following error -

..........................
Create mesh for time = 0

Time = 0
Selecting thermodynamics package ePsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>


--> FOAM FATAL ERROR:
Not implemented

From function basicThermo::h()
in file basicThermo/basicThermo.C at line 259.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::basicThermo::h() in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/libbasicThermophysicalModels.so"
#3
in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/wallHeatFluxRho"
#4 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#5
in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/wallHeatFluxRho"
Aborted

..............

Please Help

Regards,
Ishan

[ishaninair]

Hi,

I tried solution given above for calculating heat flux in rhoCentralFoam. But it gives following error-

....
Create time

Create mesh for time = 0

Time = 0
Selecting thermodynamics package ePsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>


--> FOAM FATAL ERROR:
Not implemented

From function basicThermo::h()
in file basicThermo/basicThermo.C at line 259.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::basicThermo::h() in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/lib/linuxGccDPOpt/libbasicThermophysicalModels.so"
#3
in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/wallHeatFluxRho"
#4 __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#5
in "/home/ishan/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linuxGccDPOpt/wallHeatFluxRho"
Aborted

...................

Please Help!!

Regards,
Ishan

[maddalena]

Hello everyone,
I would like to have the convective heat transfer on a wall patch. Following from the definition: alpha = wallHeatFlux / (T_wall-T_ref). This can be obtained with the use of ParaFoam, see here. What about if I modify wallHeatFluxLaminar in order to calculate alpha at the end of my simulation? This is what I would do:

forAll(wallAlphaCoeff.boundaryField(), patchi)
        {
                wallAlphaCoeff.boundaryField()[patchi] =
                    wallHeatFlux.boundaryField()[patchi]
                        /(T.boundaryField()[patchi] - Tref);
        }

where:

dimensionedScalar Tref 
( 
transportProperties.lookup("Tref") 
);

However, this has two main problems:

1. OpenFOAM complains about the T.boundaryField()[patchi] - Tref implementation. It does not like the - operator, but how should I write it?
2. If I will have a division by zero. Thus I should add something like:

forAll(wallAlphaCoeff.boundaryField(), patchi)
if T.boundaryField()[patchi] == Tref
      wallAlphaCoeff.boundaryField()[patchi] = 0.0;
else
      {
                wallAlphaCoeff.boundaryField()[patchi] =
                    wallHeatFlux.boundaryField()[patchi]
                        /(T.boundaryField()[patchi] - Tref);
        }

Of course, this is not written in a correct way. Is there anyone that can help me to understand what I am doing wrong?

mad