[harryk]

Dear cfd-online board,
I am quite new to the OpenFOAM-world and I have already used ANSYS Fluent. I want to model a heat conduction problem in a solid body, where I prescribe a temperature-dependent source term using fvOptions. The solver "laplacianFoam" has so far been sufficient. In the next step, I want to limit the total amount of energy that can be generated by the source term in each individual cell (modelling a pyrolysis process where the amount of available solid fuel is used up after some time). For that purpose, I would like to have a variable (e.g. Qtot), that is initialized with = 0 at every cell at t=0 and where the current source term can be added at the end of every timestep. In ANSYS, I would solve this using 2 UDMs (user defined memory): one to store the latest source term for each cell (stored whenever the source term is calculated, i.e. for each iteration) and one for Qtot, with an additional routine that is called at the end of each timestep to add the source term to Qtot.
How can I do something similar in OpenFOAM? (store latest source term using the routine in fvOptions, add source term to Qtot at the end of the timestep)
Afterwards, I would like to be able to plot Qtot using paraFoam.
Thank you for your help in advance - I hope I chose the right forum for this question
Best regards,
harryk