Postprocessing in multiple subdomains
5 weeks ago by
I'm currently working on magnetostatics with a convenient vector potential formulation in a multi-material domain. When we recover the magnetic induction vector by projecting the curl of the vector potential into a CG1 space, obviously all components of the magnetic induction are continuous across the material interface. But I want to allow a jump of the vector components, which are not orientated normally to the interface.
I'm not keen on projecting the magnetic induction vector into a DG0 space, because I want the vector components to be continuous within the material subdomains.
Is there a simple way of recovering the magnetic induction vector which fulfills these requirements (continuous within material subdomain, jump allowed at interface)?
Community: FEniCS Project
Please login to add an answer/comment or follow this question.