Skip to content

Custom membrane flux

Description

In the current state, species with the same name in different compartments get to have a Dirichlet-Dirichlet interface condition to make them continuous. This is not optimal because many interesting things may happen at the intersection of the domains. Dirichlet-Dirichlet transmission conditions may also result in bugs where continuity is not required but equal fluxes (#33 (closed) related?).

Proposal

Other possible settings are possible. Using the math parser, we could give the user the power to decide the transmission conditions for species between two domains.

How to test the implementation?

  • Compare a manual Dirichlet-Dirichlet BC with past results.
  • Neumann-Neumann BC should have similar results to Dirichlet-Dirichlet BC on simple cases.
  • Add a test with a reaction happening on the interface.

Related issues

See #33 (closed)