[hiuluom]

Hi everybody,

When I tried wallHeatFlux the buoyantSimpleFoam OF 2.1.1 in hotroom tutorial it was no problem with thermotype:

thermoType      hPsiThermo<pureMixture<constTransport<specieThermo<hConstThermo<perfectGas>>>>>;

pRef            100000;

mixture
{
    specie
    {
        nMoles          1;
        molWeight       28.9;
    }
    thermodynamics
    {
        Cp              1000;
        Hf              0;
    }
    transport
    {
        mu              1.8e-05;
        Pr              0.7;
    }
}

But in sonicFoam with prism tutorial has difference thermo:

thermoType      ePsiThermo<pureMixture<constTransport<specieThermo<eConstThermo<perfectGas>>>>>;

mixture
{
    specie
    {
        nMoles          1;
        molWeight       28.9;
    }
    thermodynamics
    {
        Cv              717.5;
        Hf              0;
    }
    transport
    {
        mu              1.8e-05;
        Pr              0.7;
    }
}

I used the wallHeatFluxRho and wallHeatflux (default) which was provide in forum, OF 2.1.1 and had been a problem:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : wallHeatFlux -latestTime
Date   : Jul 25 2015
Time   : 06:53:06
Host   : "CompEng"
PID    : 2979
Case   : /home/huynh/OpenFOAM/huynh-2.1.1/run/prism
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0.0003
Time = 0.0003
Selecting thermodynamics package ePsiThermo<pureMixture<constTransport<specieThermo<eConstThermo<perfectGas>>>>>

--> FOAM FATAL ERROR: 
Not implemented
    From function basicThermo::h()
    in file basicThermo/basicThermo.C at line 260.
FOAM aborting
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Foam::basicThermo::h() in "/opt/openfoam211/platforms/linuxGccDPOpt/lib/libbasicThermophysicalModels.so"
#3  
 in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/wallHeatFlux"
#4  __libc_start_main in "/lib/tls/i686/cmov/libc.so.6"
#5  
 in "/opt/openfoam211/platforms/linuxGccDPOpt/bin/wallHeatFlux"
Aborted

therefore, I copy the thermotype from buoyantSimpleFoam to sonicFoam calculation wallHeatFlux again. It was run but the result equal to zero.

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : wallHeatFluxRho -latestTime
Date   : Jul 25 2015
Time   : 22:05:09
Host   : "CompEng"
PID    : 4590
Case   : /home/huynh/OpenFOAM/huynh-2.1.1/run/prism
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0.0003

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

Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
    Prt             1;
}


Wall heat fluxes [W]
prismWall 0

End

I saw the source code wallHeatFlux which was compute by gradient of enthalpy h corresponding hPsiThermo. Otherwise, the ePsiThermo with internal energy e was not described in equation wallHeatFlux. I replace the enthalpy h by internal energy e and compile it unsuccessful.

fvc::interpolate(RASModel->alphaEff())*fvc::snGrad(h)

I would like to calculate wallHeatFlux for sonicFoam but it cannot. Could anybody can show me the solve for this, please?

Thank you so much