[tyhero]

Hello, everyone, I compiled a heat source with fvmodels. I want to know if I input the heat flux density in this way, that is, the quantity in W. Here is my code.

energySource
{
    type            coded;
    selectionMode   all;

    field          h;

    codeInclude
    #{
        #include "fvCFD.H"
    #}

    codeAddSup
    #{
    #};

    codeAddRhoSup
    #{
        scalar q0 = 3e6;
        scalar l = 0.04;
        scalar au = 0.04;

        scalarField& heSource = eqn.source();

        const scalarField& V = mesh().V();
        const scalarField& celly = mesh().C().component(1);
        heSource -= q0*(exp((-1)*(celly+au+0.028)/l))*V;
        
    #};
}