[ChrisBa]

Must I define something in the controlDict?

[ChrisBa]

I changed the fvOptions a little bit and now have this error:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.1
Exec   : buoyantBoussinesqSimpleFoam
Date   : Feb 14 2017
Time   : 13:51:56
Host   : "cfd-MEDION"
PID    : 13396
Case   : /home/cfd/OpenFOAM/cfd-4.1/run/elektronen
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


SIMPLE: convergence criteria
    field p_rgh     tolerance 0.01
    field U     tolerance 0.0001
    field T     tolerance 0.01
    field G     tolerance 1e-06
    field "(k|epsilon|omega)"     tolerance 0.001

Reading thermophysical properties

Reading field T

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

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

Reading field alphat


Reading g

Reading hRef
Calculating field g.h

No MRF models present

Radiation model not active: radiationProperties not found
Selecting radiationModel none
Creating finite volume options from "constant/fvOptions"

Selecting finite volume options model type scalarCodedSource
    Source: heatSource
    - selecting all cells
    - selected 17600 cell(s) with volume 1.8e-10


--> FOAM FATAL IO ERROR: 
keyword fields is undefined in dictionary "/home/cfd/OpenFOAM/cfd-4.1/run/elektronen/constant/fvOptions.heatSource.scalarCodedSourceCoeffs"

file: /home/cfd/OpenFOAM/cfd-4.1/run/elektronen/constant/fvOptions.heatSource.scalarCodedSourceCoeffs from line 27 to line 88.

    From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const
    in file db/dictionary/dictionary.C at line 441.

FOAM exiting

The fvOption:

heatSource
{
    type            scalarCodedSource;

    active          true;
   

    scalarCodedSourceCoeffs
    {
        selectionMode   all;
        fieldNames      (T);
        name            sourceTime;
         
        
        codeAddSup
        {
            scalarField& TSource = eqn.source();


           
            const scalar power = 1000;
            const scalar Radius = 0.5;
            const scalar rho = 1000;
            const scalar Cp = 4.19;
            const scalar Uvol = 400;
            const scalar Iel = 1000;
            const scalar electrons = 10000;
            const scalar xCenter = 3;
            const scalar yCenter = 0.01;
            const scalar zCenter = 2;             
           
           

          
            const List<point>& cf = p.Cf();
            const scalar xCF = cf[c][0];
            const scalar XF = xCenter - xCF;
            const scalar check+ = xCenter + Radius; 
            const scalar check- = xCenter - Radius;

            if ((xCF < check+) && (xCF > check-))
            {
                const scalar r2 = XF * XF; 
                const scalar re2 = Radius * Radius; 
                const scalar factor = r2/re2; 
                const scalar pre = exp(factor); 
                const scalar qfa_ [c] = power * exp;
                const scalar TzuF_ [c] = (1/(rho*Cp))*qfa   
            }
            else 
            {
                const scalar qfa_ [c] = 0;
            };
            
           
            
            
            
            
            
                       
         

            TSource -= Tzuf; 
            Pout << "***codeAddSup***" << endl;

        };
 
        codeSetValue
        {
            Pout<< "**codeSetValue**" << endl;
        };

        // Dummy entry. Make dependent on above to trigger recompilation
        code
        {
           $codeInclude
           $codeCorrect
           $codeAddSup
           $codeSetValue
        };
    }

    sourceTimeCoeffs
    {
        // Dummy entry
    }
}

I don't add something in the controlDict.
I don't use setFields or something like that, because I want to use all cells.

Maybe somebody can help me.

[tomf]

Hi,

This part explains what OpenFOAM is looking for:

--> FOAM FATAL IO ERROR:
keyword fields is undefined in dictionary "/home/cfd/OpenFOAM/cfd-4.1/run/elektronen/constant/fvOptions.heatSource.scalarCodedSourceCoeffs"

It looks for a keyword, that is called fields. I only see a fieldNames keyword in your fvOptions file. It should probably be changed to fields and than that error will not show up anymore.

Regards,
Tom

[ChrisBa]

Thanks for your reply tom.
I changed it in fields.

heatSource
{
    name            sourceTime; //Nur der Name 
    type            scalarCodedSource; // ginge auch Vector
    active          yes;

    scalarCodedSourceCoeffs 
    {
        selectionMode   all; //Alle Zellen werden ausgewählt

        fields          (T); // In welches Feld er etwas hinzufügen soll 
/*
Valid functionEntries are for codeCorrect:

8

inputMode  
codeStream
include
calc
remove
includeIfPresent
includeEtc
includeFunc 
*/
        codeCorrect
        {
        
        };
        
 
        codeAddSup
        {
            scalarField& TSource = eqn.source();


           
            const scalar power = 1000;
            const scalar Radius = 0.5;
            const scalar rho = 1000;
            const scalar Cp = 4.19;
            const scalar Uvol = 400;
            const scalar Iel = 1000;
            const scalar electrons = 10000;
            const scalar xCenter = 3;
            const scalar yCenter = 0.01;
            const scalar zCenter = 2;             
           
           

          
            const List<point>& cf = p.Cf();
            const scalar xCF = cf[c][0];
            const scalar XF = xCenter - xCF;
            const scalar check+ = xCenter + Radius; 
            const scalar check- = xCenter - Radius;

            if ((xCF < check+) && (xCF > check-))
            {
                const scalar r2 = XF * XF; 
                const scalar re2 = Radius * Radius; 
                const scalar factor = r2/re2; 
                const scalar pre = exp(factor); 
                const scalar qfa_ [c] = power * exp;
                const scalar TzuF_ [c] = (1/(rho*Cp))*qfa   
            }
            else 
            {
                const scalar qfa_ [c] = 0;
            };
            
           
            
            
            
            
            
                       
         

            TSource -= Tzuf; 
            Pout << "***codeAddSup***" << endl;

        };
 
        codeSetValue
        {
            Pout<< "**codeSetValue**" << endl;
        };

        // Dummy entry. Make dependent on above to trigger recompilation
        code
        {
           $codeInclude
           $codeCorrect
           $codeAddSup
           $codeSetValue
        };
    }

    sourceTimeCoeffs
    {
        // Dummy entry
    }
}

But I have a new problem.

Selecting finite volume options model type scalarCodedSource
    Source: heatSource
    - selecting all cells
    - selected 17600 cell(s) with volume 1.8e-10


--> FOAM FATAL IO ERROR: 
Attempt to return dictionary entry as a primitive

file: /home/cfd/OpenFOAM/cfd-4.1/run/elektronen/constant/fvOptions.heatSource.scalarCodedSourceCoeffs.codeCorrect

    From function virtual Foam::ITstream& Foam::dictionaryEntry::stream() const
    in file db/dictionary/dictionaryEntry/dictionaryEntry.C at line 83.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::IOerror::abort() at ??:?
#2  Foam::dictionaryEntry::stream() const at ??:?
#3  Foam::fv::CodedSource<double>::read(Foam::dictionary const&) at ??:?
#4  Foam::fv::option::adddictionaryConstructorToTable<Foam::fv::CodedSource<double> >::New(Foam::word const&, Foam::word const&, Foam::dictionary const&, Foam::fvMesh const&) at ??:?
#5  Foam::fv::option::New(Foam::word const&, Foam::dictionary const&, Foam::fvMesh const&) at ??:?
#6  Foam::fv::optionList::reset(Foam::dictionary const&) at ??:?
#7  Foam::fv::optionList::optionList(Foam::fvMesh const&, Foam::dictionary const&) at ??:?
#8  Foam::fv::options::options(Foam::fvMesh const&) at ??:?
#9  Foam::fv::options::New(Foam::fvMesh const&) at ??:?
#10  ? at ??:?
#11  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12  ? at ??:?
Abgebrochen (Speicherabzug geschrieben)

If i write in fvOptions:

 
codeCorrect
        {
            # 
        };

I got this message:

Radiation model not active: radiationProperties not found
Selecting radiationModel none
Creating finite volume options from "constant/fvOptions"



--> FOAM FATAL ERROR: 
Unknown functionEntry '' in "/home/cfd/OpenFOAM/cfd-4.1/run/elektronen/constant/fvOptions" near line 46

Valid functionEntries are :

8
(
inputMode
codeStream
include
calc
remove
includeIfPresent
includeEtc
includeFunc
)


    From function static bool Foam::functionEntry::execute(const Foam::word&, Foam::dictionary&, Foam::Istream&)
    in file db/dictionary/functionEntries/functionEntry/functionEntry.C at line 79.

FOAM exiting

Can someone explain me how I use codeCorrect?

Regards
Chris

[ChrisBa]

Hello all,

I hope somebody can help me with the scalarCodedSource.
I tried many things, but it doesn't work and I can't figure out why.
my fvOptions file:

heatSource
{
    name            sourceTime; 
    type            scalarCodedSource; 
    active          yes;

    scalarCodedSourceCoeffs 
    {
        selectionMode   all;

        fields          (T);

        codeCorrect
        #{
            Pout << "***codeCorrect***" << endl;
 
        #};      
 
        codeAddSup
        #{
            scalarField& TSource = eqn.source();


           
            const scalar power = 1000;
            const scalar Radius = 0.5;
            const scalar rho = 1000;
            const scalar Cp = 4.19;
            const scalar Uvol = 400;
            const scalar Iel = 1000;
            const scalar electrons = 10000;
            const scalar xCenter = 3;
            const scalar yCenter = 0.01;
            const scalar zCenter = 2;             
           
           

          
            const List<point>& cf = p.Cf();
            const scalar xCF = cf[c][0];
            const scalar XF = xCenter - xCF;
            const scalar check+ = xCenter + Radius; 
            const scalar check- = xCenter - Radius;

            if ((xCF < check+) && (xCF > check-))
            {
                const scalar r2 = XF * XF; 
                const scalar re2 = Radius * Radius; 
                const scalar factor = r2/re2; 
                const scalar pre = exp(factor); 
                const scalar qfa_ [c] = power * exp;
                const scalar TzuF_ [c] = (1/(rho*Cp))*qfa   
            }
            else 
            {
                const scalar qfa_ [c] = 0;
            };
            
           
            
            
            
            
            
                       
         

            TSource -= Tzuf; 
            Pout << "***codeAddSup***" << endl;

        #};
 
        codeSetValue
        #{
            Pout<< "**codeSetValue**" << endl;
        #};

        // Dummy entry. Make dependent on above to trigger recompilation
        code
        {
           $codeInclude
           $codeCorrect
           $codeAddSup
           $codeSetValue
        };
    }

    sourceTimeCoeffs
    {
        // Dummy entry
    }
}

// ************************************************************************* //

The Error is this one:

 
Create time

Create mesh for time = 0


SIMPLE: convergence criteria
    field p_rgh     tolerance 0.01
    field U     tolerance 0.0001
    field T     tolerance 0.01
    field G     tolerance 1e-06
    field "(k|epsilon|omega)"     tolerance 0.001

Reading thermophysical properties

Reading field T

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Creating turbulence model

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

Reading field alphat


Reading g

Reading hRef
Calculating field g.h

No MRF models present

Radiation model not active: radiationProperties not found
Selecting radiationModel none
Creating finite volume options from "constant/fvOptions"



--> FOAM FATAL ERROR: 
Attempt to return primitive entry ITstream : /home/cfd/OpenFOAM/cfd-4.1/run/elektronen/constant/fvOptions.heatSource.scalarCodedSourceCoeffs.codeCorrect, line 31, IOstream: Version 2.0, format ASCII, line 0, OPENED, GOOD
    primitiveEntry 'codeCorrect' comprises 
        on line 31 the verbatim string "\
            Pout << \"***codeCorrect***\" << endl;\
 \
        "
 as a sub-dictionary

    From function virtual const Foam::dictionary& Foam::primitiveEntry::dict() const
    in file db/dictionary/primitiveEntry/primitiveEntry.C at line 189.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::primitiveEntry::dict() const at primitiveEntry.C:?
#3  Foam::dictionary::substituteScopedKeyword(Foam::word const&) at ??:?
#4  Foam::entry::New(Foam::dictionary&, Foam::Istream&) at ??:?
#5  Foam::dictionary::read(Foam::Istream&, bool) at ??:?
#6  Foam::dictionary::dictionary(Foam::fileName const&, Foam::dictionary const&, Foam::Istream&) at ??:?
#7  Foam::dictionaryEntry::dictionaryEntry(Foam::keyType const&, Foam::dictionary const&, Foam::Istream&) at ??:?
#8  Foam::entry::New(Foam::dictionary&, Foam::Istream&) at ??:?
#9  Foam::dictionary::read(Foam::Istream&, bool) at ??:?
#10  Foam::dictionary::dictionary(Foam::fileName const&, Foam::dictionary const&, Foam::Istream&) at ??:?
#11  Foam::dictionaryEntry::dictionaryEntry(Foam::keyType const&, Foam::dictionary const&, Foam::Istream&) at ??:?
#12  Foam::entry::New(Foam::dictionary&, Foam::Istream&) at ??:?
#13  Foam::dictionary::read(Foam::Istream&, bool) at ??:?
#14  Foam::dictionary::dictionary(Foam::fileName const&, Foam::dictionary const&, Foam::Istream&) at ??:?
#15  Foam::dictionaryEntry::dictionaryEntry(Foam::keyType const&, Foam::dictionary const&, Foam::Istream&) at ??:?
#16  Foam::entry::New(Foam::dictionary&, Foam::Istream&) at ??:?
#17  Foam::dictionary::read(Foam::Istream&, bool) at ??:?
#18  Foam::operator>>(Foam::Istream&, Foam::dictionary&) at ??:?
#19  Foam::IOdictionary::readFile(bool) at ??:?
#20  Foam::IOdictionary::IOdictionary(Foam::IOobject const&) at ??:?
#21  Foam::fv::options::options(Foam::fvMesh const&) at ??:?
#22  Foam::fv::options::New(Foam::fvMesh const&) at ??:?
#23  ? at ??:?
#24  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#25  ? at ??:?
Abgebrochen (Speicherabzug geschrieben)

I can't understand this because I used the example of the headerfile.
Thanks for xour help.

Regards,
Chris

[ChrisBa]

Morning all,

i solved the problem with the "primitive entry". Maybe somebody can help me with the code.
I need 3 zones in z. The first one should be one face of my mesh and for the other two there should be different values for each cell.

Regards
Chris

[ChrisBa]

Hy all,

sorry to bother you, but I have a problem.
I only want a zone for a heatSource. So I used topoSetDict and cellSet, but in ParaView the whole mesh is used.

Here is my fvOptions:

harmonic
{   name            sourceTime;
    type            scalarCodedSource;
    active          true;
    selectionMode   cellZone;
    cellZone        heatSource;
    scalarCodedSourceCoeffs
    { 

        fields      (T);
        selectionMode   cellZone;
        cellZone        heatSource;
        codeInclude
        #{

        #};

        codeCorrect
        #{
            Pout<< "**codeCorrect**" << endl;
        #};

        codeAddSup
        #{
            const Time& time = mesh().time();
            const scalarField& V = mesh_.V();
            const vectorField& C = mesh_.C();
            scalarField& TSource = eqn.source();
            forAll(C, i)
            {
                TSource += 1*V[i];
            }
            
            Pout << "***codeAddSup***" << endl;
        #};

        codeSetValue
        #{
            Pout<< "**codeSetValue**" << endl;
        #};

        // Dummy entry. Make dependent on above to trigger recompilation
        code
        #{
            $codeInclude
            $codeCorrect
            $codeAddSup
            $codeSetValue
        #};
    }

    sourceTimeCoeffs
    {
        // Dummy entry
    }
}


// ************************************************************************* //

Here is the topoSetDict:

actions
(

    {
       name    heatSource;
       type cellSet;
       action  new;
       source boxToCell;
       sourceInfo
       {
            box   (0 0 0) (0.2 0.2 0.2);
       }

    }

    {   
        name    heatSource;
        type    cellZoneSet;
        action  new;
        source  setToCellZone;
        sourceInfo
        {
            set heatSource;
        }
    }

);

I'm using buoyantBoussinesqSimpleFoam and Openfoam 4.1

Can somebody help me or explain me, why this is the result?

regards

[ChrisBa]

Hy All

Can I select the surfaces of my mesh?
Because i want to define a source only on the surface of the mesh.

regards,
Chris

[student666]

Instead of a volumetric heat source, I would rather say that you need a heatflux boundary condition on a surface.

Maybe these two bounday condition on the following link can be of help:

buoyantSimpleFoam: fixed wall heat flux

If you use the banana trick on fvOptions for scalarSemiImplctSource

energySource
{
    type            scalarSemiImplicitSource;
    active on; //to switch or not heat source 
    //timeStart 0;
    //duration  20;
    selectionMode cellZone;
    cellZone c0;

    scalarSemiImplicitSourceCoeffs{
        selectionMode   banana_all;//cellZone;
        cellZone        c0;
    
        volumeMode      absolute;//specific;//
        injectionRateSuSp
        {
            h           (0 0.1); //   q in [W]; or in [W/m³] if you use specific mode
        }
    }
}

openfoam gives you other choices, and you may see that you can select points

Selecting finite volume options model type scalarSemiImplicitSource
    Source: energySource


--> FOAM FATAL IO ERROR: 
banana_all is not in enumeration: 
4
(
all
cellSet
cellZone
points
)

but It is a volumetric heat source, and I think that "points" means internal points.
I'm not able to figure it out, how manage it on a boundary surface, maybe my thinking is not correct but:
try to think you set only a point, with a volumetric heat source on a surface that is a boundary of your computational domain too; then others points what would be, zeroGradient?
ok, but how the heat will be transfered?
Because other points has zeroGradient so no heat flux along the surface as well as outside of your domain.
So you have a heat flux from your heat point on the surface (that belongs to a volumetric cell) toward the center point of the volumetric cell that it belongs to, but as above you're setting a heat flux, not a heat source. Maybe it works for your purpose.

You may even think to use thermaBaffles:

https://openfoamwiki.net/index.php/F...hermal_baffles
https://openfoam.org/release/2-3-0/thermal/

Hope it helps.

Regards.

[sitajeje]

Dear Foamers,

I read through the thread and a few other webpages, and found the following to options in injectionRateSuSp to create a heat source by fvOption with scalarSemiImplicitSource:

injectionRateSuSp
{
T (0.1 0); //DT = PHI''' / (rho * cp)
}

or

injectionRateSuSp
{
h (0.1 0); //W absolute, W/m3 specific
}

I wonder which one to choose for chtMultiRegionFoam?

I chose the second one get the computation result as shown in the photo. The case with OpenFOAM 4.0 is also attached. It is a long thin cuboid with heat source of 10W and a flat large cuboid lying above. There is no thermal contact resistance between each other. The solid is located in the middle of a closed (noSlip for U boundary condition) and insulated box (zeroGradient for T). Only a quarter is simulated, thus two faces of the computation domain are symmetry, and the internal temperature and velocity field can be observed.

I simulated for two seconds, so that 20J heat was generated. However the increase of the total internal energy of the two solid cuboids increased only 5J (taking the temperature of the solids is uniform at the highest temperature in domain, and the thermal properties are given in the case). I wonder where are the other 75% of energy? Is it due to the symmetry boundary condition that only 25% of the heat was generated?

Thank you very much for any help in advance!

Sitajeje

[piu58]

I don't think your problem is due to the b.c. Symmetric simply generates a slip boundary and does not expand the region in any way.

Sent from my HTC One SV using CFD Online Forum mobile app

[sitajeje]

Hallo Uwe,

thank you very much for your quick reply! I figured out by trial and error that it is the second option is correct:

injectionRateSuSp
{
h (0.1 0); //W absolute, W/m3 specific
}

It is true that it is nothing to do with the boundary condition. I made the symmetryPlane as Wall, so that it is a complete insulated box. But the temperature is the same, i.e. only 5J emerges in the solid domain, and 15J was missing.

Sorry everyone,

It was a wrong estimation of the internal energy. I checked with a simpler geometry, it was correct.

Sitajeje

[zfaraday]

Hi Sitajeje,

First of all, both implementations should give you the same results since the only difference between them is the comment (the characters behind //), so it has no impact on the calculations.

Can you, please, upload your fvSolution specification? Or at least tell us the tolerancies you used for enthalpy for all your regions. Convergence is usually very slow in chtMultiregion* solvers and you would need to set a very low tolerance and let the solver run during several iterations. You can use non ortogonal correctors in order to force the solver to take more iterations.

Best regards,

Alex

[sitajeje]

Alex,

thank you very much for your reply! I tried
injectionRateSuSp
{
T (0.1 0); //DT = PHI''' / (rho * cp)
}
but there was no heat generated. I checked the solveSolid.H in applications/solvers/heatTransfer/chtMultiRegionFoam/solid, and it is written "fvScalarMatrix hEqn" there. I suppose that is the reason. What is your experience?

I used the following:
h
tolerance 1e-7;
relTol 0.1

hFinal
tolerance 1e-7;
relTol 0;

Can you download my file? I uploaded it in my last input.

The velocity and temperature didn't need many iterations. But the pressure frequently needed 1000 iterations. So I changed the "p_rghFinal relTol" from 0 to 0.01. The results were the same.

What is your experience? I didn't understand your comments on tolerance and convergence.

Thank you very much in advance!

Best regards,

Sitajeje

[ChrisBa]

Hy Foamers,

I work with scalarCodedSource.
My solver is buoyantBoussinesqSimpleFoam and I use the version OpenFoam 4.1.
I would like to create a volumentric source that uses a constant power.
I have deduced from the energy equation that I need the specific power first.
This value I then divide by the density rho [kg / m^3] and the heat capacity Cp [J / (kgK)].
By integrating the equation, I would have to multiply this value by the cell volume.

My questions are the following:

1. How does the solver handle the fvOptions value?
2. Does the fvOptions(T) function return the calculated value?
3. Is my derivation correct?
4. Is it correct if I use eqn.source() as my variable, if my performance is constant?
5. Can I use this fvOptions also for the solver buoyantBoussinesqPimpleFoam?

Thanks in advance for your help.

Regards,
ChrisBa

[gu1]

Hello, I hope someone can help me ...

If I want to set a fixed pressure gradient, can I? to a cyclic domain (pipe).

Could anyone teach me how to configure these options in my fvOptions?

Thanks.

[CalumR]

Hi there


I am trying to implement a heat sink through scalarCodedSource using a table of calculated Heat fluxes as my input. Can anyone give me any advice on the implementation of this? I am familiar with using tables to specify boundary conditions but am unsure how to implement them in the source.


Any help would be appreciated!

[Vishsel]

Hi all,

My case was porous medium with heat transfer. But I only need to generate heat source on wall of the solid body (ie. on surface). Not from fluid body.

These are my BC for wall surface which i have to create heat source (at this BC's)
wall-mp-1
wall-mp-2
wall-dp-1
wall-dp-2

And I have used cyclicAMI type for interfacing the fluid regions and for measure the pressure values.

So I succeed with cyclicAMI, now i have to create a heat source at wall(ie @surface). I need to only check heat distribution from those BC's.

Solver i am using is porousSimpleFoam..
Running mpirun -np no.ofprocessor porousSimpleFoam -parallel

Now my case is with heat transfer,incompressible fluid, steady state

I have an input like heat flux (W/m2) and heat generation (W) value.

I am getting this error while running. I know i have to include energy equation.. I don't have an idea about header file creation..

--> FOAM Warning :
From function void option::checkApplied() const
in file fvOptions/fvOption.c at line 368
Source wall-mp-1 defined for field h but never used

please view this comment for fvOptions,fvSchemes & fvSolution files
energy equation for porousSimpleFoam

Please correct me what i have to modify in this EEqn.H file

{
    volScalarField& he = thermo.he();

    fvScalarMatrix EEqn
    (
        fvm::div(phi, he)
      + fvc::div(phi, K)
      + (
            he.name() == "e"
          ? fvc::div
            (
                fvc::absolute(phi/fvc::interpolate(rho), U),
                p,
                "div(phiv,p)"
            )
          : -dpdt
        )
      - fvm::laplacian(turbulence->alphaEff(), he)
     ==
      + fvOptions(rho, he)
    );

    EEqn.relax();

    fvOptions.constrain(EEqn);

    EEqn.solve();

    fvOptions.correct(he);

    thermo.correct();
    radiation->correct();
}

Thankyou in advance
Vishsel

[saidc.]

Hi,

I'm trying to import a source term which depends on height of the rode.

z: height (0 to 0.3 meter)
profile: 0.00000444*z^4 - 0.00025507*z^3 + 0.00144848*z^2 + 0.07409711*z + 1.64833117

codedSource
{
    type            scalarCodedSource;
    selectionMode   cellZone;
    cellZone        MYREGION;

    fields          (h);
    name            sourceTime;

    codeInclude
    #{
    #};
    codeCorrect
    #{
    Pout<< "**codeCorrect**" << endl;
    #};

    codeAddSup
    #{
    const scalarField& Zcoord = mesh.C().component(Vector::Z);
    scalarField& heSource = eqn.source();
    forAll(Zcoord, i)
    {
        heSource += (0.00000444*pow(Zcoord[i],4)) - (0.00025507*pow(Zcoord[i],3)) + (0.00144848*pow(Zcoord[i],2)) + (0.07409711*Zcoord[i]) + 1.64833117;    
    }
    #};

    codeConstrain
    #{
    Pout<< "**codeCorrect**" << endl;
    #};
    code
    #{
            $codeInclude
            $codeCorrect
            $codeAddSup
            $codeSetValue
    #};
    sourceTimeCoeffs
    {
    // Dummy entry
    }

It gives this error:

Selecting finite volume options type scalarCodedSource
    Source: codedSource
    - selecting cells using cellZone MYREGION
    - selected 80339 cell(s) with volume 0.074947679



...



 Solving for solid region heatPipesRegion
DICPCG:  Solving for h, Initial residual = 1, Final residual = 8.5586193e-07, No Iterations 50
Min/max T:293.14761 305.44315
Using dynamicCode for fvOption::sourceTime at line 24 in "/mnt/c/Users/saidc/Desktop/codedTest/constant/MYREGION/fvOptions.codedSource"
Could not load "/mCould not load "/mnt/c/Users/saidCould not load "/mnt/c/Users/saidc/DeskCould not load "/mnt/c/Users/saidc/Desktop/codedTnt/c/Users/saidc/Desktop/codedTest/dynamicCc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPIntop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/liode/platforms/linux64GccDPInt32Opt/lib/est/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b0libsourceTimeb/libsourceTime_t32Opt/lib/libsourceT5a445f885e6ec8482cb3441c896bime_ecc90b05a445f885e6ec8482cb3441c896bd_ecc90b05a445f885e6ec8482cb3441c896bd535d.so"
/mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so: cannot open shared object file: No such file or directory
ecc90b05a445f885e6ec8482cb3441c896bd535d.so"
/mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so: cannot open shared object file: No such file or directory
d535d.so"
/mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so: cannot open shared object file: No such file or directory
535d.so"
/mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so: cannot open shared object file: No such file or directory
Creating new library in "dynamicCode/sourceTime/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so"
Invoking wmake libso /mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/sourceTime
wmake libso /mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/sourceTime
    ln: ./lnInclude
    dep: codedFvOptionTemplate.C
    Ctoo: codedFvOptionTemplate.C
[DESKTOP-5Q0IUMJ:04407] 3 more processes have sent help message help-btl-vader.txt / cma-permission-denied
[DESKTOP-5Q0IUMJ:04407] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
/mnt/c/Users/saidc/Desktop/codedTest/constant/MYREGION/fvOptions.codedSource: In member function ‘virtual void Foam::fv::sourceTimeFvOptionscalarSource::addSup(Foam::fvMatrix<double>&, Foam::label)’:
/mnt/c/Users/saidc/Desktop/codedTest/constant/MYREGION/fvOptions.codedSource:41:29: error: invalid use of member function ‘const Foam::fvMesh& Foam::fv::option::mesh() const’ (did you forget the ‘()’ ?)
make: *** [/usr/lib/openfoam/openfoam2112/wmake/rules/General/transform:35: Make/linux64GccDPInt32Opt/codedFvOptionTemplate.o] Error 1
[0] 
[0] 
[0] --> FOAM FATAL IO ERROR: (openfoam-2112 patch=220610)
[0] Failed wmake "dynamicCode/sourceTime/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so"
[0] 
[0] 
[0] file: constant/MYREGION/fvOptions.codedSource at line 24.
[0] 
[0]     From void Foam::codedBase::createLibrary(Foam::dynamicCode&, const Foam::dynamicCodeContext&) const
[0]     in file db/dynamicLibrary/codedBase/codedBase.C at line 240.
[0] 
FOAM parallel run exiting

The code didn't compiled. Any help will be appreciated.

Kind regards,

[massive_turbulence]

Originally Posted by saidc.

Hi,

I'm trying to import a source term which depends on height of the rode.

z: height (0 to 0.3 meter)
profile: 0.00000444*z^4 - 0.00025507*z^3 + 0.00144848*z^2 + 0.07409711*z + 1.64833117

codedSource
{
    type            scalarCodedSource;
    selectionMode   cellZone;
    cellZone        MYREGION;

    fields          (h);
    name            sourceTime;

    codeInclude
    #{
    #};
    codeCorrect
    #{
    Pout<< "**codeCorrect**" << endl;
    #};

    codeAddSup
    #{
    const scalarField& Zcoord = mesh.C().component(Vector::Z);
    scalarField& heSource = eqn.source();
    forAll(Zcoord, i)
    {
        heSource += (0.00000444*pow(Zcoord[i],4)) - (0.00025507*pow(Zcoord[i],3)) + (0.00144848*pow(Zcoord[i],2)) + (0.07409711*Zcoord[i]) + 1.64833117;    
    }
    #};

    codeConstrain
    #{
    Pout<< "**codeCorrect**" << endl;
    #};
    code
    #{
            $codeInclude
            $codeCorrect
            $codeAddSup
            $codeSetValue
    #};
    sourceTimeCoeffs
    {
    // Dummy entry
    }

It gives this error:

Selecting finite volume options type scalarCodedSource
    Source: codedSource
    - selecting cells using cellZone MYREGION
    - selected 80339 cell(s) with volume 0.074947679



...



 Solving for solid region heatPipesRegion
DICPCG:  Solving for h, Initial residual = 1, Final residual = 8.5586193e-07, No Iterations 50
Min/max T:293.14761 305.44315
Using dynamicCode for fvOption::sourceTime at line 24 in "/mnt/c/Users/saidc/Desktop/codedTest/constant/MYREGION/fvOptions.codedSource"
Could not load "/mCould not load "/mnt/c/Users/saidCould not load "/mnt/c/Users/saidc/DeskCould not load "/mnt/c/Users/saidc/Desktop/codedTnt/c/Users/saidc/Desktop/codedTest/dynamicCc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPIntop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/liode/platforms/linux64GccDPInt32Opt/lib/est/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b0libsourceTimeb/libsourceTime_t32Opt/lib/libsourceT5a445f885e6ec8482cb3441c896bime_ecc90b05a445f885e6ec8482cb3441c896bd_ecc90b05a445f885e6ec8482cb3441c896bd535d.so"
/mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so: cannot open shared object file: No such file or directory
ecc90b05a445f885e6ec8482cb3441c896bd535d.so"
/mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so: cannot open shared object file: No such file or directory
d535d.so"
/mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so: cannot open shared object file: No such file or directory
535d.so"
/mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so: cannot open shared object file: No such file or directory
Creating new library in "dynamicCode/sourceTime/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so"
Invoking wmake libso /mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/sourceTime
wmake libso /mnt/c/Users/saidc/Desktop/codedTest/dynamicCode/sourceTime
    ln: ./lnInclude
    dep: codedFvOptionTemplate.C
    Ctoo: codedFvOptionTemplate.C
[DESKTOP-5Q0IUMJ:04407] 3 more processes have sent help message help-btl-vader.txt / cma-permission-denied
[DESKTOP-5Q0IUMJ:04407] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
/mnt/c/Users/saidc/Desktop/codedTest/constant/MYREGION/fvOptions.codedSource: In member function ‘virtual void Foam::fv::sourceTimeFvOptionscalarSource::addSup(Foam::fvMatrix<double>&, Foam::label)’:
/mnt/c/Users/saidc/Desktop/codedTest/constant/MYREGION/fvOptions.codedSource:41:29: error: invalid use of member function ‘const Foam::fvMesh& Foam::fv::option::mesh() const’ (did you forget the ‘()’ ?)
make: *** [/usr/lib/openfoam/openfoam2112/wmake/rules/General/transform:35: Make/linux64GccDPInt32Opt/codedFvOptionTemplate.o] Error 1
[0] 
[0] 
[0] --> FOAM FATAL IO ERROR: (openfoam-2112 patch=220610)
[0] Failed wmake "dynamicCode/sourceTime/platforms/linux64GccDPInt32Opt/lib/libsourceTime_ecc90b05a445f885e6ec8482cb3441c896bd535d.so"
[0] 
[0] 
[0] file: constant/MYREGION/fvOptions.codedSource at line 24.
[0] 
[0]     From void Foam::codedBase::createLibrary(Foam::dynamicCode&, const Foam::dynamicCodeContext&) const
[0]     in file db/dynamicLibrary/codedBase/codedBase.C at line 240.
[0] 
FOAM parallel run exiting

The code didn't compiled. Any help will be appreciated.

Kind regards,

If you could post the snippets of the lines below that might help more. Can you run this on a single CPU or thread? Start small and work up from it.

file: constant/MYREGION/fvOptions.codedSource at line 24.
[0]
[0] From void Foam::codedBase::createLibrary(Foam::dynamicCode&, const Foam::dynamicCodeContext&) const
[0] in file db/dynamicLibrary/codedBase/codedBase.C at line 240.