Hello,
I am trying to model a heterogeneous reaction by specifying a Flux boundary condition that depends on the surface concentration, gamma, which changes with time: gamma(t) = gamma(t=0) - gamma_reacted(t). gamma_reacted(t) is just the integral from 0 to t of the surface flux. I see that the function timeint can only be evaluated after the simulation is finished, so I wanted to use a time discrete study and define gamma_reacted as a variable by adding timestep*flux to it at each step. Unfortunately I cannot define gamma_reacted = gamma_reacted + timestep*flux, but that is what I would like to do.
I am trying to model a heterogeneous reaction by specifying a Flux boundary condition that depends on the surface concentration, gamma, which changes with time: gamma(t) = gamma(t=0) - gamma_reacted(t). gamma_reacted(t) is just the integral from 0 to t of the surface flux. I see that the function timeint can only be evaluated after the simulation is finished, so I wanted to use a time discrete study and define gamma_reacted as a variable by adding timestep*flux to it at each step. Unfortunately I cannot define gamma_reacted = gamma_reacted + timestep*flux, but that is what I would like to do.