Pre-processing of GroupOfEdgesOf Cohomology Basis for Complex Geometry
Hi Christophe,
We are trying to run a transient simulation of two (or more) insulated spiral pancake coils connected in series using an \vec{H}-\phi
formulation, see below for the geometry. Unfortunately, we could not reproduce the problem with easier geometries...
This model setup then has two cohomology basis functions. We follow Section 5 of the Gmsh cohomology paper to create a cohomology basis of 2 elements with evident physical interpretation:
We then set up a very simple .pro file that only calls InitSolution
in the Resolution
and enforce a current of 1A in the first cut and 0A in the second (effectively "disabling" it). This should yield the following initial solution (log scale) with currents only in the conducting parts:
The above picture is for a spiral of 12 turns. When increasing the number of turns to 13, we have non-zero currents in the air domain:
A few more symptoms:
- Removing the second cut region, i.e.
DOM_airHoleCut = Region[{ }];
, there are no currents in non-conducting domains. - Switching to the cohomology basis originally created from the
HomologyComputation
(Regions 8 and 9 as commented in the.pro
) without the rearranging to make homology and cohomology compatible, again no currents in the non-conducting domain. - Below 13 turns (exclusive), no air currents. Above 13 turns, air currents.
Do you have an idea what's going on? Thanks for the help!
Please find below the simple .pro and two meshes for the two cases. We tested it with the official GetDP 3.5.0
release.
Thank you!
Erik