[jdw135]

Originally Posted by HPE

Let's see what we can do during the next week. (At least something for you to progress further rather than a complete solution.)

I've now done quite a bit more reading and I think I can do what I need with a scalarCodedSource. My only problem is that I'm not sure how to implement it, and if my code is doing what I want it to do.

First, do I make a file called "fvOptions" in the /system folder, or do I add this code to the controlDict file in the functions section? Looking at the OpenFOAM API here

https://www.openfoam.com/documentati...e.html#details

it says to put it in the controlDict file. When I do this, the case runs as normal, but there is no difference in the output. I am wondering if there is a problem with my code that is causing an issue. Here is my code:

energySource
    {
        type            scalarCodedSource;

        active          yes;

        scalarCodedSourceCoeffs
        {
            selectionMode   all;
            fieldNames      (h);
            name            sourceTime;

            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& qq = -(pS*(1-(mach*mach))1005)/(1.4*(mach*mach)*287.058);
                scalarField& heSource = eqn.source();
                heSource += fvc::ddt(qq);
            #};

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

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

I tried making an "fvOptions" file, but I kept getting an error about "name" not being found. Can someone please help with this? There isn't a lot of documentation and I'm pulling my hair out here.