[Asghari_M]

Hello everyone!
There is a transient heat source in my problem according to the following relation:

S(heat)=N × sqrt(t) ;

That t is time and N is an integer multiplier.
I can use the following code in the heat source term definition indicated at scalarCodedSource in fvOptions only for one multiplier (for example N=1):

N=1;
heSource -= N* sqrt(time.value())*V;

However, I want to do 30 runs for a geometry case. Every run is corresponding to one n.
Really, multiplier (N) Varies from one to 30.After the simulation has ended, the data of the last time-step is saved to a separate directory. Similar task has been described for inlet velocity at following link by PyFoam (and in 2.2 parameter variation):
https://openfoamwiki.net/index.php/Contrib/PyFoam ------->2.2parameter variation
In the other words, I want to vary source term value without direct contact with OpenFoam2.3.1 in a loop structure. How can I do this operation?
Thanks in advance for any help.