Problem with boundary conditions in Gridap and GridapGmsh

General Usage
1

Dear all,
I am trying to solve some 2D electrostatic problem (Poisson Equation) using Gripap.jl and gmsh for mesh generation. I solve for the potential in space bounded by conductors, which have surface potentials prescribed as Dirichlet boundary conditions U=1 on the outer boundary and U=0 on the inner electrode.

fig1
fig1800×600 115 KB

When the inserted electrode is “small” (only one mesh cell per boundary segment), the solver seems to ignore the potential assigned to the small electrode (U=0) as shown in the figure below.

fig2
fig2800×600 59.8 KB

Apparently, the mesh cell number on each side is important as also illustrated by the figure below, where only the sides with two cells per side are correctly represented.

fig3
fig3800×600 107 KB

My source code is available in the following gist:

Is there some solver parameter to fix this issue?

I am using the gmsh boundary markers to assign the boundary conditions. Note that in Fenicsx ( fenicsproject.org ) I can solve on the same mesh without issues.

My software versions are
Julia 1.8.2
Gridap 0.17.14
GridapGmsh 0.6.1

Thanks in advance for any help.

2

Dear @rouckas,

Physical Curve("C0", 10) = {bnd1()};
Physical Curve("C1", 20) = {bnd2()};

Have you tried assigning the tag C0 also to the physical points in the closure of bnd1() (see this small example on a disk)?

Sorry for the late answer. The community is more active on gitter.