[YannGuyot]

Dear all,
I am writing this message since I am struggling with an issue.
I want to add steam in my buoyantsimpleFoam simulations. For this I added a scalartransport term in the controldict. I want to have the source (a specific volume in the mesh) dependent on the computed temperature T, so I used the scalarCodedSource f in the fvOptions.
No error is coming out, the problem is that the coded source term is not taken into account.
Here is the codes I used, I know the toposet defining my source volume is working fine since it is also a source term for heat and this part works fine.


output :

OAM Warning :
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source steam defined for field H2O but never used

controlDict :

functions
{
    steam
    {
          type                 scalarTransport;
          functionObjectLibs   ("libsolverFunctionObjects.so");
          field                H2O;
          writeControl         timeStep;
          schemesField         k;
          writeInterval        1;
          resetOnStartUp       false;
          autoSchemes          true;
      }
}

fvOptions:

steam
{
    type            scalarCodedSource;
    active          yes;
    name            sourceTime;

    scalarCodedSourceCoeffs
    {
        selectionMode   cellSet;
        cellSet         Porosity;
        fields          (H2O);
        codeInclude
        #{
        #};
        codeCorrect
        #{
        #};
        codeAddSup
        #{
              const volScalarField& Tm = mesh_.lookupObject<volScalarField>("T");
              const volScalarField Tcc = Tm.mesh().C() & vector(1,1,1);
              const scalarField& V = mesh_.V();
              scalarField& H2OSource = eqn.source();
              const labelList& cellIDs = cells();
              forAll(cellIDs, i)
              {
                  label cellI = cellIDs[i];
                  H2OSource[cellI] = -0.00000245*(HEATBROILER-(0.61*HEATBROILER*(1000+20*(20-(Tm[cellI]-273.15)))/1000-(0.228/1000)*(Tm[cellI]-273.15)*(Tm[cellI]-273.15)))*DENSITY/0.3*V[cellI];
              }
        #};
        codeSetValue
        #{
        #};
    }
    sourceTimeCoeffs
    {
        $scalarCodedSourceCoeffs;
    }
}

Thank you very much for any help
Yann

[Manouu]

Hi Yann, I have the same problem did you solve it please???

[SHUBHAM9595]

Originally Posted by YannGuyot

Dear all,
I am writing this message since I am struggling with an issue.
I want to add steam in my buoyantsimpleFoam simulations. For this I added a scalartransport term in the controldict. I want to have the source (a specific volume in the mesh) dependent on the computed temperature T, so I used the scalarCodedSource f in the fvOptions.
No error is coming out, the problem is that the coded source term is not taken into account.
Here is the codes I used, I know the toposet defining my source volume is working fine since it is also a source term for heat and this part works fine.


output :

OAM Warning :
    From function virtual void Foam::fv::option::checkApplied() const
    in file cfdTools/general/fvOptions/fvOption.C at line 125
    Source steam defined for field H2O but never used

controlDict :

functions
{
    steam
    {
          type                 scalarTransport;
          functionObjectLibs   ("libsolverFunctionObjects.so");
          field                H2O;
          writeControl         timeStep;
          schemesField         k;
          writeInterval        1;
          resetOnStartUp       false;
          autoSchemes          true;
      }
}

fvOptions:

steam
{
    type            scalarCodedSource;
    active          yes;
    name            sourceTime;

    scalarCodedSourceCoeffs
    {
        selectionMode   cellSet;
        cellSet         Porosity;
        fields          (H2O);
        codeInclude
        #{
        #};
        codeCorrect
        #{
        #};
        codeAddSup
        #{
              const volScalarField& Tm = mesh_.lookupObject<volScalarField>("T");
              const volScalarField Tcc = Tm.mesh().C() & vector(1,1,1);
              const scalarField& V = mesh_.V();
              scalarField& H2OSource = eqn.source();
              const labelList& cellIDs = cells();
              forAll(cellIDs, i)
              {
                  label cellI = cellIDs[i];
                  H2OSource[cellI] = -0.00000245*(HEATBROILER-(0.61*HEATBROILER*(1000+20*(20-(Tm[cellI]-273.15)))/1000-(0.228/1000)*(Tm[cellI]-273.15)*(Tm[cellI]-273.15)))*DENSITY/0.3*V[cellI];
              }
        #};
        codeSetValue
        #{
        #};
    }
    sourceTimeCoeffs
    {
        $scalarCodedSourceCoeffs;
    }
}

Thank you very much for any help
Yann

Hi Yann
Have you tried

active          true;

instead of

active          yes;

Also, if you look at the syntax implementation can u try replacing

 H2OSource[cellI] = -0.00000245*(HEATBROILER-(0.61*HEATBROILER*(1000+20*(20-(Tm[cellI]-273.15)))/1000-(0.228/1000)*(Tm[cellI]-273.15)*(Tm[cellI]-273.15)))*DENSITY/0.3*V[cellI];

with

 H2OSource[cellI] += -0.00000245*(HEATBROILER-(0.61*HEATBROILER*(1000+20*(20-(Tm[cellI]-273.15)))/1000-(0.228/1000)*(Tm[cellI]-273.15)*(Tm[cellI]-273.15)))*DENSITY/0.3*V[cellI];