[Nayeem12]

I am investigating a subchannel of Molten Salt Reactor. The design is essentially a graphite pipe, through which, heat generating molten salt is flowing.
For the heat generation, in ANSYS, I have simply added following equation in the flowing molten salt cell zone:
1.57 * 213.5e6 * sin((y*PI)/(3.96 [m])) * 1 [W m^-3]

Now for validation, I have extracted the exact mesh and imported to OpenFOAM, added exactly same boundary conditions, and for adding the heat source to the fluid, I have added an fvOptions to the fluid region with following lines:

energySource
{
    type            scalarCodedSource;
    active          true;
    name            sourceTime; 

    scalarCodedSourceCoeffs
    {
        selectionMode   cellZone;
        cellZone        fuel_region; 
        volumeMode      absolute; 
        fields          (h);          
        
        fieldNames      (h);
        name            sourceTime;
        
        codeInclude
        #{

        #};

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

        codeAddSupRho
        #{
            const scalarField& V = mesh_.V();
            const scalarField& celly = mesh_.C().component(1);
            scalarField& hSource = eqn.source();

            // Initialize min and max with the first element of celly
            scalar minCelly = celly[0];
            scalar maxCelly = celly[0];

            forAll(celly, i)
            {
                hSource[i] = -(1.57*213.5e6*sin((3.14159*celly[i])/3.96))*V[i];

                // Update min and max values
                if (celly[i] < minCelly)
                {
                    minCelly = celly[i];
                }
                if (celly[i] > maxCelly)
                {
                    maxCelly = celly[i];
                }
            }

            // Print out the min and max celly values
            Pout << "Minimum celly value: " << minCelly << endl;
            Pout << "Maximum celly value: " << maxCelly << endl;

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

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

limitTemperature
{
    type    limitTemperature;
    active  true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        min    839;
        max    2400;
    }
}

But the result is far from what ANSYS predicted. Even the temperature profile is completely different.

To check if I have added the boundary conditions properly, I have changed the scalarCodedSource to scalarSemiImplicitSource and added a constant specific heat generation term to the fluid region with following fvOptions:

scalarSource
{
    active  true;
    type    scalarSemiImplicitSource;
    scalarSemiImplicitSourceCoeffs
    {
        volumeMode  specific;
        selectionMode   all;
        injectionRateSuSp
        {
            h   (213.5e6 0);
        }
    }
}


limitTemperature
{
    type    limitTemperature;
    active  true;
    limitTemperatureCoeffs
    {
        selectionMode all;
        min    839;
        max    2400;
    }
}

And did the same for ANSYS, and the temperature profile came out as perfectly same.

This led me to believe I am possibly doing something wrong with the sinusoidal volumetric heat source I am using in fvOptions?