GridapMHD.jl icon indicating copy to clipboard operation
GridapMHD.jl copied to clipboard

Using Gmsh transfinite mehses with curved geometries

Open FRUrgorri opened this issue 2 years ago • 7 comments

I am trying to launch the function "main" in debug mode using transfine meshes designed with Gmsh.

I have prepared two simple examples: A square-section channel and a pipe. In the first case everything seems to work ok. However, in the second case I obtain a bound error when main is defining the test space V_p (line 83 of src/Main.jl V_p = TestFESpace(Ω,reffe_p))

Perhaps it is the way in which I define the pipe mesh but I do not see the difference between the two cases other than the curvature (which naturally introduce non orthogonal cells).

You can reproduce the error with the test.jl file (you need to add GridapGmsh to the project first). I include the .msh files. I also include the .geo files I used to generete them in case they help
Test.zip

Thanks for the help

FRUrgorri avatar Mar 10 '22 09:03 FRUrgorri