GalerkinToolkit.jl
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
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
.
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.