Hello,

I am looking into solving the 2 following coupled PDEs:

(i) Div[Grad(V)] = -q/cste

(ii) Div[q.Grad(V)] = 0

With q and V being functions: q(x;y) & V(x;y)

With boundary conditions on V all around the 2D domain

- (i) is a basic Poisson problem and (ii) is a variable coefficient Poisson problem

Independant resolution of (i) and (ii) works fine with FEniCS, thanks to relevant documentation.

For both cases: Assuming q is known as a starting point, I can get V as a result on my mesh/domain - I tried to couple (i) & (ii) as follows

STEP1: Assume a first q, set BCs on V and solve (ii) => V1 as an output

STEP2: Using V1, calculate q with (i) => q2 as an output

STEP3: Using q2, set (same) BCs on V and solve (ii) => V3 as an output

Etc…

Unfortunately, result doesn’t converge - I also tried to work with some coupling FEniCS tutorials, but they are more complicated than my need and I couldn’t adapt them

Does anyone know a simple way to couple (i) and (ii) in the same variational statement, and get q and v efficiently ?

In advance Thank you for your help !