I have a hopefully simple question, given a nonlinear DAE that takes itself as a parameter: j(j, stuff) = 0 is there a way to represent this with FEniCS?
If it’s important, the definition is:
eta = phis - phie - Uocp - F * Rfilm * j
j = i0 * (fem.exp((1 - alpha) * F * eta / (R * Tref)) - fem.exp(-alpha * F * eta / (R * Tref)))
Where phis and phie are trial functions and i0 and Uocp are functions of other trial functions.