[superkelle]

Hi,
I had some trouble to implement a scalarCodedSource in fvOptions for rhoPimpleFoam. But I got it working and I want to share it with you, maybe it is helpful for someone, since I read some other posts that had trouble with it. Here is my example for the sqareBendLiq tutorial ( /rhoPimpleFoamTutorials/RAS/squareBendLiq/ ) , you can just create an fvOptions file in system directory and run rhoPimpleFoam after blockMesh.

It will create an sinodial moving heat source upward the stream. It has no physical meaning and was just for testing. Have fun, I will attach a gif.

https://media.giphy.com/media/lrc8Xa...G8aE/giphy.gif

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      fvOptions;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
energySource
{
    type            scalarCodedSource;    //scalarSemiImplicitSource
    active          true; //for activation either don't use any "active" or use value "true", for deactivation "false"
    name            sourceTime;

    scalarCodedSourceCoeffs
    {
        fields          (e);

        fieldNames        (e);

        name            sourceTime;
       
        selectionMode   all;

        codeInclude
        #{

        #};

        codeCorrect
        #{
        #};
       
        codeAddSup
        #{
            const Time& time = mesh().time();
           
           
           
            const scalarField& V = mesh_.V();
           
            const vectorField& C = mesh_.C();
                       
            const scalar endTime = 0.5;
                       
            const scalar xt = -0.45 + time.value()*2;

            const scalar yt = -0.05+0.01*sin(time.value()*100.0);
           
            const scalar rt = 0.00005;
           
           
            //const scalar  zt = -0.05;
                       
            scalarField& heSource = eqn.source();

            forAll (C,i)
            {   
                const scalar & Cx =C[i].component(0);
                const scalar & Cy =C[i].component(1);
                const scalar & Cz =C[i].component(2);
               
                if ( (sqr(Cx-xt) + sqr(Cy-yt) ) <= rt && time.value() < endTime)
                {   
                    heSource[i] = - 2e10 * V[i];
                }
               
                /*
                if (Cx <= xt && Cy >= 0.8*yt && Cy <= 1.2*yt && time.value() < endTime)
                {   
                    heSource[i] = - 5e9 * V[i];
                }
                */
               
           
            }
        #};
       
        codeSetValue
        #{
        #};
       
        code
        #{
            $codeInclude
            $codeCorrect
            $codeAddSup
            $codeSetValue
        #};
   
    }

    sourceTimeCoeffs
    {
       selectionMode   all;
    }

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

[jdw135]

Hi,

I am trying to get a scalarCodedSource working as well, but I am getting a warning during runtime that says "Source ICM defined for field qq but never used". Below is my fvOptions code. I would greatly appreciate any help!

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      fvOptions;
}

energySource
    {
        type            scalarCodedSource;
        name            sourceTime;

        active          yes;

        scalarCodedSourceCoeffs
        {
            selectionMode   all;
            fields          (qq);

            codeInclude
            #{
            #};

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

            codeAddSup
            #{
                const Time& time = mesh().time();
                const volScalarField& rhoS = mesh_.lookupObject<volScalarField>("rho");
                const volVectorField& uS = mesh_.lookupObject<volVectorField>("U");
                const volScalarField& pS = mesh_.lookupObject<volScalarField>("p");
                const vectorField& C = mesh_.C();
                const scalarField& V = mesh_.V();
                scalarField& SoS = sqrt(1.4*(pS/rhoS));
                scalarField& mach = uS/SoS;
                scalarField& rhoQ = (pS*(1-(mach*mach))1005)/(1.4*(mach*mach)*287.058);
                scalarField& qq = eqn.source();
                qq -= rhoQ*V;
            #};

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

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