[Ali Blues]

Hi,
I'm trying to add an explicit source term (Sp) to the scalarTransportFoam, say for the field T. The source Term involves the product of T with another volScalarField object, lets call it E, which actually remains constant through the simulation. I'm trying to use the type scalarCodedSource in fvOptions as in https://github.com/OpenFOAM/OpenFOAM.../CodedSource.H


Since I'm using an explicit source term, T.oldTime() is used as an estimate.
For E, since it remains constant during the simulations, I thought of putting it in the constant directory and read it from there, and that's where I'm having getting difficulties. Here is the code:

energySource 
     { 
         type            scalarCodedSource; 
         active          yes; 
         scalarCodedSourceCoeffs 
         { 
             selectionMode   all; 
             fieldNames      (T); 
             redirectType    sourceTime; 
  
             codeInclude 
             #{ 
             #}; 
  
             codeCorrect 
             #{ 
             #}; 
  
             codeAddSup 
             #{ 
                 const Time& time = mesh().time();             
                 const scalarField& V = mesh_.V(); 
  
                const volScalarField& Temp = mesh_.lookupObject<volScalarField>("T"); 

                volScalarField E  
             (  
             IOobject  
             (  
                 "E",  
                 runTime.constant(),  
                 mesh_,  
                 IOobject::MUST_READ,  
                 IOobject::NO_WRITE  
             ),  
             mesh_  
             ); 
                 
             scalarField& heSource = eqn.source(); 
             heSource -=E*Temp.oldTime()*V;           
  
             #}; 
  
             codeSetValue 
             #{                
             #}; 
              
             code 
             #{ 
                 $codeInclude 
                 $codeCorrect 
                 $codeAddSup 
                 $codeSetValue 
             #}; 

         } 
         sourceTimeCoeffs 
         { 
             
         } 
     }

Q1) I get the following error: ‘runTime’ was not declared in this scope. So I appreciate any suggestions on how to fix this?
Q2) Following the link for the CodedSource.H, I was wondering why can I interchangeably use mesh().time() and mesh_.time() ?

Thanks,
Ali

[JackW]

add

const Time& time = mesh().time();

then replace

runTime.timeName()

with

time.timeName()

As for your second question, mesh_ is declared in the compilation. mesh() is a function that returns the same object.

[Ali Blues]

Thanks Jack for the reply.
As I noted, field E remains constant during the simulations, and it is computed prior to running the scalarTransportFoam solver. So I prefer to read it from a different path say the constant/ directory.
Giving time.timeName() as path would only work for the first time step if I include the field E in the 0/ directory. but for the later time steps in would give error since field E is non-existent.

Best
Ali

[Ali Blues]

just realized time is actually of object Time. So accessing the constant directory is possible via time.constant(). So reading field E is as follows:

volScalarField E                
(
          IOobject                
          ( 
            "E", 
            time.constant(),                     
            mesh_,                    
            IOobject::MUST_READ,                    
            IOobject::NO_WRITE                
           ),                
          mesh_                
);

[jpeter3]

You need to get mesh and time objects, they do not exist in every class under this names.

Maybe try to Access them through velocity field, something like:

U_.mesh() instead of mesh

and

U_.time() instead of runTime