A Numerical Framework for Diffusive Transport in Rotational Symmetric Systems with Discontinuous Interlayer Conditions

A numerical method based on a finite element approach for modelling interface conditions in the context of diffusive transport is presented. Such interface conditions may result in concentration jumps due to drug partitioning and/or mass transfer effects. The paper extends previous work in a one-dimensional multilayer setting to the two-dimensional one of cylindrical coordinates with rotational symmetry, applicable e.g. to drug transport in and out of optical lenses with possible surface barriers. Simulations of such a case are presented. The method derived is also applicable to a two-dimensional Cartesian setting.