Refactor code structure for improved readability

Author: nikoscham

assets/gmsh_rect_quad_struct.png

@@ -188,11 +188,6 @@ <h2 id="features"><a name="Features"></a>Features</h2>
                   margin-bottom: 4px;
                 " /></a
             >) files
-            <i
-              ><a href="https://github.com/FEAScript/FEAScript-core/issues/5" target="_blank"
-                >(under development)</a
-              ></i
-            >
           </li>
         </ul>
       </li>
@@ -246,6 +241,10 @@ <h2 id="tutorials"><a name="Tutorials"></a>Tutorials</h2>
           >API (web worker)</a
         >
         |
+        <a href="https://feascript.com/tutorials/HeatConduction2DFinGmsh.html" target="_blank"
+          >API (Gmsh mesh)</a
+        >
+        |
         <a href="https://feascript.com/tutorials/HeatConduction2DFinPlatform.html" target="_blank"
           >XML (FEAScript platform)</a
         >

assets/gmsh_rect_quad_unstruct.png

@@ -20,4 +20,8 @@
     <loc>https://feascript.com/tutorials/HeatConduction2DFinPlatform.html</loc>
     <priority>0.8</priority>
   </url>
+  <url>
+    <loc>https://feascript.com/tutorials/HeatConduction2DFinGmsh.html</loc>
+    <priority>0.8</priority>
+  </url>
 </urlset>

index.html

@@ -145,7 +145,7 @@ <h2 id="gmshfileimport"><a name="gmshfileimport"></a>Gmsh File Import</h2>
     model.setMeshConfig({
       parsedMesh: result,
       meshDimension: "2D",
-      elementOrder: "linear",
+      elementOrder: "quadratic",
     });
 
     // Define boundary conditions using Gmsh physical group tags
@@ -213,7 +213,7 @@ <h3 id="gmshgeofile"><a name="gmshgeofile"></a>Example Gmsh .geo File</h3>
 // With Physical Lines for boundary labeling
 // Use the command "gmsh rect.geo -2" to generate the mesh
 
-lc = 0.2; // Characteristic length (mesh density)
+lc = 0.7; // Characteristic length (mesh density)
 
 // Points (x, y, z, mesh size)
 Point(1) = {0, 0, 0, lc};     // Bottom left
@@ -297,7 +297,11 @@ <h2 id="generatedmesh"><a name="generatedmesh"></a>Generated Mesh</h2>
     <h2 id="results"><a name="results"></a>Results</h2>
     <p>
       Below is the 2D contour plot of the computed temperature distribution. This plot is generated in real
-      time using FEAScript.
+      time using FEAScript. Please note that solutions computed on unstructured meshes like this may exhibit
+      small numerical differences compared to solutions on structured orthogonal meshes. This occurs because
+      derivative calculations in non-orthogonal elements are inherently less precise due to the Jacobian
+      transformation process. These small differences are expected and acceptable for most engineering
+      applications, but may be noticeable in regions with steep temperature gradients.
     </p>
 
     <div id="solutionPlot"></div>
@@ -327,9 +331,13 @@ <h2 id="results"><a name="results"></a>Results</h2>
       //  importGmshQuadTri,
       //  plotSolution,
       //  printVersion,
+      //  logSystem,
       //} from "../../FEAScript-core/src/index.js";
 
       window.addEventListener("DOMContentLoaded", async () => {
+        // Enable debug logging mode
+        //logSystem("debug");
+
         // Print FEAScript version in the console
         printVersion();
 
@@ -352,15 +360,11 @@ <h2 id="results"><a name="results"></a>Results</h2>
         // Parse the mesh file first
         const result = await importGmshQuadTri(meshFile);
 
-        //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({
           parsedMesh: result,
           meshDimension: "2D",
-          elementOrder: "linear",
+          elementOrder: "quadratic",
         });
 
         // Define boundary conditions using Gmsh physical group tags
@@ -369,8 +373,6 @@ <h2 id="results"><a name="results"></a>Results</h2>
         gmshModel.addBoundaryCondition("2", ["convection", 1, 20]); // top boundary
         gmshModel.addBoundaryCondition("3", ["symmetry"]); // left boundary
 
-        //console.log("GMSH Boundary Conditions:", gmshModel.boundaryConditions);
-
         // Set solver method
         //model.setSolverMethod("jacobi");
         gmshModel.setSolverMethod("lusolve");
@@ -379,8 +381,6 @@ <h2 id="results"><a name="results"></a>Results</h2>
         const { solutionVector: gmshSolutionVector, nodesCoordinates: gmshNodesCoordinates } =
           gmshModel.solve();
 
-        //console.log("GMSH Solution Vector:", gmshSolutionVector);
-
         // Plot the GMSH solution
         plotSolution(
           gmshSolutionVector,

sitemap.xml

@@ -2,7 +2,7 @@
 // With Physical Lines for boundary labeling
 // Use the command "gmsh rect.geo -2" to generate the mesh
 
-lc = 0.2; // Characteristic length (mesh density)
+lc = 0.7; // Characteristic length (mesh density)
 // Note: In the structured case, the lc = 0.2 no longer directly determines the mesh size,
 // as the Transfinite Curve commands explicitly set the number of elements along each boundary