<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style type="text/css" id="owaParaStyle"></style>

</head>

<body fpstyle="1" ocsi="0">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">

<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">

<div>Hello everyone,</div>

<div><br>

</div>

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

<div><br>

<div>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 ?</div>

<div><br>

</div>

<div>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) ?</div>

<div><br>

</div>

<div>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 ?</div>

</div>

<div><br>

</div>

<div>I would sincerely appreciate your help on those questions. Faithfully yours,</div>

<div><br>

</div>

<div>Guillaume Dilasser</div>

</div>

</div>

</body>

</html>