Add tutorial for 2D Stokes flow in a lid-driven cavity

Author: nikoscham

assets/lid-driven-cavity-2d-stokes-u-results.png

@@ -171,6 +171,14 @@ <h1 class="top">
         Heat conduction through a wall
         <a href="https://feascript.com/tutorials/heat-conduction-1d-wall.html">tutorial</a>
       </p>
+      <img
+        src="./assets/lid-driven-cavity-2d-stokes-u-results.png"
+        alt="Results of lid-driven cavity 2D Stokes flow tutorial"
+      />
+      <p class="image-caption">
+        Lid-driven cavity (2D Stokes flow)
+        <a href="https://feascript.com/tutorials/lid-driven-cavity-2d-stokes.html">tutorial</a>
+      </p>
       <img
         src="./assets/solidification-front-2d-results.png"
         alt="Results of solidification front propagation in a 2D domain tutorial"
@@ -412,6 +420,7 @@ <h2 id="features"><a name="Features"></a>Features</h2>
         <ul>
           <li>Front propagation</li>
           <li>Heat conduction</li>
+          <li>Incompressible Stokes flow</li>
         </ul>
       </li>
       <li>
@@ -666,6 +675,12 @@ <h2 id="tutorials"><a name="Tutorials"></a>Tutorials</h2>
             style="width: 16px; height: 16px; vertical-align: middle; margin-left: 2px; margin-bottom: 4px"
         /></a>
       </li>
+      <li>
+        <span class="tutorial-name">Lid-Driven Cavity (2D Stokes Flow)</span>
+        <a href="https://feascript.com/tutorials/lid-driven-cavity-2d-stokes.html" class="badge badge-basic"
+          >Basic HTML</a
+        >
+      </li>
       <li>
         <span class="tutorial-name">Solidification Front Propagation in a 2D Domain</span>
         <a href="https://feascript.com/tutorials/solidification-front-2d.html" class="badge badge-basic"

assets/stokes-2d-lid-driven.png

@@ -219,8 +219,8 @@ <h2 id="solvingwithfeascript"><a name="Solving with FEAScript"></a>Solving with
       });
 
       // Define boundary conditions
-      model.addBoundaryCondition("0", ["constantValue", 1]); // Left boundary, u(0) = 1
-      model.addBoundaryCondition("1", "zeroGradient"); // Right boundary, zero gradient (du/dx = 0)
+      model.addBoundaryCondition("0", ["constantValue", 1]); // Left boundary
+      model.addBoundaryCondition("1", "zeroGradient"); // Right boundary
 
       // Solve
       model.setSolverMethod("lusolve");

assets/stokes-2d-lid-driven.svg

@@ -207,8 +207,8 @@ <h2 id="solvingwithfeascript"><a name="Solving with FEAScript"></a>Solving with
       });
 
       // Apply boundary conditions
-      model.addBoundaryCondition("0", ["convection", 1, 25]);
-      model.addBoundaryCondition("1", ["constantTemp", 5]);
+      model.addBoundaryCondition("0", ["convection", 1, 25]); // Left boundary
+      model.addBoundaryCondition("1", ["constantTemp", 5]); // Right boundary
 
       // Solve
       model.setSolverMethod("lusolve");

index.html

@@ -194,10 +194,10 @@ <h2 id="gmshfileimport"><a name="gmshfileimport"></a>Gmsh File Import</h2>
       });
 
       // Apply boundary conditions using Gmsh physical group tags
-      model.addBoundaryCondition("0", ["constantTemp", 200]); // bottom boundary
-      model.addBoundaryCondition("1", ["constantTemp", 200]); // right boundary
-      model.addBoundaryCondition("2", ["convection", 1, 20]); // top boundary
-      model.addBoundaryCondition("3", ["symmetry"]); // left boundary
+      model.addBoundaryCondition("0", ["constantTemp", 200]); // Bottom boundary
+      model.addBoundaryCondition("1", ["constantTemp", 200]); // Right boundary
+      model.addBoundaryCondition("2", ["convection", 1, 20]); // Top boundary
+      model.addBoundaryCondition("3", ["symmetry"]); // Left boundary
 
       // Solve
       model.setSolverMethod("lusolve");

tutorials/advection-diffusion-1d.html

@@ -175,10 +175,10 @@ <h2 id="webgpuimplementation">WebGPU Implementation</h2>
     });
 
     // Apply boundary conditions
-    model.addBoundaryCondition("0", ["constantTemp", 200]);
-    model.addBoundaryCondition("1", ["symmetry"]);
-    model.addBoundaryCondition("2", ["convection", 1, 20]);
-    model.addBoundaryCondition("3", ["constantTemp", 200]);
+    model.addBoundaryCondition("0", ["constantTemp", 200]); // Bottom boundary
+    model.addBoundaryCondition("1", ["symmetry"]); // Left boundary
+    model.addBoundaryCondition("2", ["convection", 1, 20]); // Top boundary
+    model.addBoundaryCondition("3", ["constantTemp", 200]); // Right boundary
 
     // Set solver method to use the GPU-accelerated Jacobi solver
     model.setSolverMethod("jacobi-gpu");

tutorials/heat-conduction-1d-wall.html

@@ -187,10 +187,10 @@ <h2 id="webworkerimplementation"><a name="Web Worker Implementation"></a>Web Wor
       });
 
       // Apply boundary conditions
-      await model.addBoundaryCondition("0", ["constantTemp", 200]);
-      await model.addBoundaryCondition("1", ["symmetry"]);
-      await model.addBoundaryCondition("2", ["convection", 1, 20]);
-      await model.addBoundaryCondition("3", ["constantTemp", 200]);
+      await model.addBoundaryCondition("0", ["constantTemp", 200]); // Bottom boundary
+      await model.addBoundaryCondition("1", ["symmetry"]); // Left boundary
+      await model.addBoundaryCondition("2", ["convection", 1, 20]); // Top boundary
+      await model.addBoundaryCondition("3", ["constantTemp", 200]); // Right boundary
 
       // Solve
       await model.setSolverMethod("lusolve");