[divyesh]

Hello Foamers,

I want to add a source Term in laplacianFoam. exp(time)*(Sqr(valueOfMesh(X))-2)

Anyone have idea how to do that?

Can I use funkySetFields?
But still, I have to define something in solver...

[simrego]

Hi!


I think you are looking for this:
https://www.openfoam.com/documentati...dedSource.html


EDIT:
I did not compiled it but I hope it'll work. I wrote it for you with loop so you can see some general way to set the source. You can do it without loop but maybe you will need it in the future so you can do it like this. I hope it helps.

 codeAddSup
#{
              const Time& time = mesh().time();
             const vectorField& C = mesh().C();

              const scalarField& V = mesh_.V();
              scalarField& Source = eqn.source();
             forAll(C, i)
             {
                            Source[i] -= exp(time)*(sqr(C[i].x())-2.0)*V[i];
             }

  #};

[divyesh]

source
{
    type    scalarCodedSource;
    
    scalarCodedSourceCoeffs
    {
        selectionMode    all;
        fields    (T);
        //name    internalField;
    codeInclude
    #{
    #};

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


 codeAddSup
#{
    const Time& time = mesh().time();
    const vectorField& C = mesh().C();
    const scalarField& V = mesh_.V();
    scalarField& source = eqn.source();
             forAll(C, i)
             {
                            source[i] -= exp(time)*(sqr(C[i].x())-2.0)*V[i];
             }

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

I put this code in fvOptions...It is showing error...

Creating finite volume options from "system/fvOptions"

Selecting finite volume options model type scalarCodedSource
    Source: source
    - selecting all cells
    - selected 50 cell(s) with volume 0.01


--> FOAM FATAL IO ERROR: 
keyword name is undefined in dictionary "/system/fvOptions.source"

[HPE]

the error message is quite clear? comment out 'name' entry?

[simrego]

Hi!


Sorry for the late answer. HPE is correct, it should be clear from the error. But there is an another problem which is the missing '}' after the scalarCodedSourceCoeffs.

[divyesh]

thanks. problem solved.
another thing i add is time.value() instead of just time.