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

Sequential two level mesh for 3D periodic puzzle piece is incorrect, and solver fails due to parallel boundary face labeling

Open jfdev001 opened this issue 10 months ago • 12 comments

As can be seen below, and based on our discussion on 2024-04-02, the sequential two level mesh for the periodic puzzle piece appears to have issues at the interface. The visualize_mesh helper function, and # TODO: comments in example.jl commit 1fa5b349c814dbcb02d4d6feccac3600cde9d17e indicate areas of the code that produce the preliminary diagnostic information we discussed today. The top of example.jl describes the naming conventions, but for reading the below files it's worth noting final_mesh vtk files are named as final_mesh_unit_cell_mesh_<info about unit cell mesh>_coarse_cell_mesh_<info about coarse cell mesh>.

The image below is from the output/final_mesh_unit_cell_mesh_3D_periodic_gmsh_puzzlepiece_geometry_triangular_refcell_coarse_cell_mesh_3D_periodic_glk_box_geometry_quad_2x2x2_refcell.vtu.

image

It is also worth inspecting output/coarse_cell_id_final_mesh_unit_cell_mesh_3D_periodic_gmsh_box_geometry_triangular_refcell_coarse_cell_mesh_3D_periodic_glk_box_geometry_quad_4x4x4_refcell.vtu since the coarse cell ids appear to be labeled correctly.

jfdev001 avatar Apr 02 '24 14:04 jfdev001