[Getdp] Various questions about superconductors, cohomology & procedural constraints in GetDP

[DILASSER Guillaume]

Hello everyone,

I am writing to (hopefully) have some response to the questions I have been asking myself about getDP and the simulations of superconducting materials.

A) In every simulation of superconducting behavior I make, I am more or less obliged to add some resistive (ohmic) elements in order to get results. Otherwise, the computation wil sooner or later end with an external PETSc error "singular matrix". The model I am currently experimenting on is derived from Christophe's helix.pro (the non linear material law is the same), in which there is also always a resistive matrix that you can't remove. Could anyone with numerical know-how explain to me why this error is happening ? Would there be another way to consider only the superconductor ?

B) In the enclosed project file I tried to automate the implementation of current constraints yet I am not sure I made it right. The system I want to work on is a multi-layered superconducting coil made of REBCO tape, since I want to be able to isolate what is happening in each separate turn, I need to proceed to a lot of cuts (one between each layer) and generate the corresponding cohomology constraints. Could anyone tell me if what I implemented is ok(-ish) ?

C) Obviously, my current implementation is not working. With more than one turn I get NaN non-linear residue, but even with just one turn I get wrong results (even is I made the same computation in another project file earlier and the results were fine). It seem that (see enclosed picture) the cohomology cut hitting the superconducting surface is the cause. One can see that I have a weird (and surely wrong) behavior at the intersection. It is almost as if (despite the current constraint I try to enforce) the total current flowing through the conductor was set to zero... Can anyone help me to understand what is happening there ?

I would sincerely appreciate your help on those questions. Faithfully yours,

Guillaume Dilasser
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20151207/ea35daa3/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Bobinette_supra_2.geo
Type: application/octet-stream
Size: 7851 bytes
Desc: Bobinette_supra_2.geo
URL: <http://www.geuz.org/pipermail/getdp/attachments/20151207/ea35daa3/attachment-0001.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Bobinette_supra_2.par
Type: application/octet-stream
Size: 3042 bytes
Desc: Bobinette_supra_2.par
URL: <http://www.geuz.org/pipermail/getdp/attachments/20151207/ea35daa3/attachment-0001.par>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Bobinette_supra_2.pro
Type: application/octet-stream
Size: 15623 bytes
Desc: Bobinette_supra_2.pro
URL: <http://www.geuz.org/pipermail/getdp/attachments/20151207/ea35daa3/attachment-0001.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: GDilasserMesh&Cut.jpg
Type: image/jpeg
Size: 107459 bytes
Desc: GDilasserMesh&Cut.jpg
URL: <http://www.geuz.org/pipermail/getdp/attachments/20151207/ea35daa3/attachment-0003.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: GDilasserStrangeBfield.jpg
Type: image/jpeg
Size: 124284 bytes
Desc: GDilasserStrangeBfield.jpg
URL: <http://www.geuz.org/pipermail/getdp/attachments/20151207/ea35daa3/attachment-0004.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: GDilasserJdensity.jpg
Type: image/jpeg
Size: 80793 bytes
Desc: GDilasserJdensity.jpg
URL: <http://www.geuz.org/pipermail/getdp/attachments/20151207/ea35daa3/attachment-0005.jpg>