[giovannibonino]

Hello everyone,

I'm a new to OF and I have some issues imposing convergence conditions in my transient runs. Please mind that I am using a customized multiphysics solver that is not part of the OF package and I have OF 4.0.

In short, I have this field named "flux" (i.e. neutron flux) and I want the run to stop when the relative variation of the flux between two subsequent time steps is less than 10e-5 everywhere in my domain.
Mathematically, this condition could be expressed as: max | flux(i+1,j)-flux(i,j) | / | flux(i,j) | < 10e-5 where "i" is the time step index and "j" is the spatial index.
Is it possible to do that by changing something in the controlDict file or would I have to modify the solver itself?

Thanks.