<div dir="ltr"><div>Hello everyone,<br><br></div><div>I've encountered a strange issue with my electrostatic computation. I'm using the V formulation with the electric potential v saved as a post-operation output (see pro file attached).<br></div><div><br></div><div>The geometry is desperately simple: two conductors  in a square surrounding medium, each with its insulation. When the conductors are round, everything's fine and the electric potential is computed correctly. But when their shape is changed to rectangular, for some reason I get nan values everywhere in the domain.<br></div><div><br></div><div>My msh file is generated using the Gmsh Python API. Attached is the corresponding Python script, where one can switch the shape from circles to rectangles. For debugging purposes, the script is written so that the minimal amount of code that is not shared by the two cases is isolated behind a conditional statement.</div><div><br></div><div>Apparently, the problem arises when and only when one or several of the left or right insulation outer curves are straight lines. I have checked that there is nothing wrong with physical regions.<br></div><div><br></div><div>In both cases I'm using the same pro file, and the command to run GetDP is the following:</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>getdp -pro <a href="http://electrostatics.pro">electrostatics.pro</a> -msh test_geom.msh -solve EleSta_vf -pos get_v<br></div></blockquote><div><br></div><div> Note that I was able to carry out the very same computation successfully  in a different, much more complicated geometry with rectangular conductors. I sincerely have no idea as to the root of this seemingly intractable problem, and therefore any help would be greatly appreciated. If needed, I can also provide the msh files I get when I'm running my Python script.</div><div><br></div><div>Other potentially relevant information below<br></div><div>- Python API from Gmsh 4.4.1 Linux 64 bits SDK <br></div><div>- GetDP 3.2.0 Linux 64 bits</div><div>- Python 3.7</div><div>- Ubuntu 18.04.3<br></div><div><br></div><div></div><div>Best regards,</div><div><br></div><div>Théo Mathurin</div><div><br></div>Research engineer<br><div>Laboratoire Systèmes Electrotechniques et Environnement (LSEE) - Université d'Artois</div><div><br></div></div>