Update boundary conditions and logging in HeatConduction2DFinGmsh.html; adjust mesh density in rect.geo; add rect_quad2.msh file - Still in debugging phase

nikoscham · nikoscham · commit 789ced2c2065 · 2025-05-08T15:51:21.000+03:00
diff --git a/tutorials/HeatConduction2DFinGmsh.html b/tutorials/HeatConduction2DFinGmsh.html
@@ -149,17 +149,11 @@ <h2 id="results"><a name="results"></a>Results</h2>
         // Parse the mesh file first
         const result = await importGmshQuadTri(meshFile);
 
+        console.log("GMSH Boundary Conditions:", result.boundaryConditions);
         console.log("GMSH Boundary Elements:", result.boundaryElements);
-
-        // Fix nodal numbering format for the solver
-        if (result.nodalNumbering.quadElements && result.nodalNumbering.quadElements.length > 0) {
-          result.nodalNumbering = result.nodalNumbering.quadElements;
-        } else if (
-          result.nodalNumbering.triangleElements &&
-          result.nodalNumbering.triangleElements.length > 0
-        ) {
-          result.nodalNumbering = result.nodalNumbering.triangleElements;
-        }
+        console.log("GMSH Nodal Numbering:", result.nodalNumbering);
+        console.log("x Coordinates:", result.nodesXCoordinates);
+        console.log("y Coordinates:", result.nodesYCoordinates);
 
         // Define mesh configuration with the parsed result
         gmshModel.setMeshConfig({
@@ -169,10 +163,10 @@ <h2 id="results"><a name="results"></a>Results</h2>
         });
 
         // Define boundary conditions using Gmsh physical group tags
-        gmshModel.addBoundaryCondition("1", ["constantTemp", 200]); // bottom boundary
-        gmshModel.addBoundaryCondition("4", ["symmetry"]); // left boundary
-        gmshModel.addBoundaryCondition("3", ["convection", 1, 20]); // top boundary
-        gmshModel.addBoundaryCondition("2", ["constantTemp", 200]); // right boundary
+        gmshModel.addBoundaryCondition("0", ["constantTemp", 200]); // bottom boundary
+        gmshModel.addBoundaryCondition("1", ["constantTemp", 200]); // right boundary
+        gmshModel.addBoundaryCondition("2", ["constantTemp", 100]); // top boundary
+        gmshModel.addBoundaryCondition("3", ["symmetry"]); // left boundary
 
         console.log("GMSH Boundary Conditions:", gmshModel.boundaryConditions);
 
diff --git a/tutorials/rect.geo b/tutorials/rect.geo
@@ -2,7 +2,7 @@
 // With Physical Lines for boundary labeling
 // Use the command "gmsh rect.geo -2" to generate the mesh
 
-lc = 0.8; // Characteristic length (mesh density)
+lc = 0.3; // Characteristic length (mesh density)
 
 // Points (x, y, z, mesh size)
 Point(1) = {0, 0, 0, lc};     // Bottom left
diff --git a/tutorials/rect_quad2.msh b/tutorials/rect_quad2.msh
