[letoppina]

Hello,

I am using OpenFOAM version 6 and chtmuliregionfoam. I am using the following code to impose a sinusoidal heat source on my solid region throughout my fvOptions file:

energySource
    {
        type            scalarCodedSource;
        active          yes;
        timeStart       0;
        duration        1e6;
        name            sourceTime; 
    selectionMode   all;

    scalarCodedSourceCoeffs
        {
        selectionMode   all; 
                volumeMode      specific;           
        fields          (h);
            codeInclude
            #{
            #};

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

            codeAddSup
            #{
                    const Time& time = mesh().time();
                    const scalarField& V = mesh_.V();
                    const vectorField& C = mesh_.C();
                    scalarField& heSource = eqn.source();

            heSource += 2.729e6*sin(2*3.14*0.1*time.value())*V;  

                Pout << "***codeAddSup***" << endl;
            #};

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

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

Is it possible to impose a square function instead of a sinusoidal one? Or impose a table of values?

Thank you in advance!!