CFD Online Discussion Forums - Howto use scalarCodedSource in fvOptions

Page 3 of 4

Show 40 post(s) from this thread on one page

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - Howto use scalarCodedSource in fvOptions (https://www.cfd-online.com/Forums/openfoam-solving/123670-howto-use-scalarcodedsource-fvoptions.html)

ChrisBa

February 14, 2017 05:55

Must I define something in the controlDict?

ChrisBa

February 14, 2017 07:58

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

Code:

/*---------------------------------------------------------------------------*\

| =========                 |                                                 |

| \\      /  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:

Code:

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	February 15, 2017 04:23

Hi,

This part explains what OpenFOAM is looking for:

Quote:

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

February 15, 2017 04:33

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

Code:

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.

Code:

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:

Code:

 

codeCorrect

        {

            # 

        };

I got this message:

Code:

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

February 15, 2017 11:03

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:

Code:

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:

Code:

 

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

February 16, 2017 01:51

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

February 18, 2017 16:17

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:

Code:

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:

Code:

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

February 24, 2017 06:27

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

February 24, 2017 20:45

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:

https://www.cfd-online.com/Forums/op...tml#post583906

If you use the banana trick on fvOptions for scalarSemiImplctSource

Code:

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

Code:

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

March 11, 2017 17:11

How to set "injectionRateSuSp" in "scalarSemiImplicitSource"?

2 Attachment(s)

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

March 12, 2017 01:50

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

March 12, 2017 10:51

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

March 13, 2017 09:32

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

March 13, 2017 11:27

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

March 31, 2017 02:02

Question about fvOptions

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	November 1, 2018 13:20

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

November 30, 2018 04:37

Coded Source table input

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

September 5, 2019 08:06

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

Code:

--> 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
https://www.cfd-online.com/Forums/op...tml#post743697

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

Code:

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

November 12, 2022 15:38

Hi,

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

Code:

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

Code:

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:

Code:

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

November 14, 2022 10:53

Quote:

Originally Posted by saidc. (Post 839255)

Hi,

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

Code:

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

Code:

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:

Code:

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.

All times are GMT -4. The time now is 10:48.

Page 3 of 4

Show 40 post(s) from this thread on one page