Help with Setting Up Transport Equation for Advected Scalar

[JhonCordova]

Hello,

I need assistance with setting up the transport equation for an advected quantity in PeleC. I have included the directive PELEC_NUM_ADV = 1 in the GNUmakefile. Subsequently, in the prob.H file, I initialized the values of this advected quantity within a sphere.

My objective is to maintain this quantity constant over time. To achieve this, I intend to introduce a source term that counteracts the advection term. I am using a specific case as a reference for adding the source term.

Could anyone provide guidance or suggestions on how to implement this effectively?

Thanks!

[baperry2]

I'm not sure if I understand what you're trying to do. Why set up an advected scalar, but then force it not to advect? You can add whatever you want in the custom source term function, so you could add a forcing that pushes the new state towards the old state, but I'm not sure if that's exactly your goal.