Solve diffusion equation and Laplace equation together

[th509]

Thanks for your contribution!

I want to solve diffusion equation with electro-migration term.

Current code with custom pde class is written without electromigration term

def evolution_rate(self, state, t=0):
   return self.diffusivity * state.laplace(bc=["periodic", [{"value": 0}, {"value": 1}]])

The electric potential (\phi) should be solved simultaneously following equation,

where c1 = 0.1 and c2 = 0.9, epsilon2 = epsilon1 * 1000

How can I solve the time-dependent diffusion equation and steady state Laplace equation together?

[david-zwicker]

I'm afraid that this is not easily possible in the current version of the package since it focuses more on solving dynamical equations (whereas in your case \phi is not a dynamical field anymore).