Merge pull request #19 from bcube-project/dev_stokes

Addition of a Stokes flow example (steady and unsteady)
bcube-project · Nov 13, 2024 · f9ad529 · f9ad529
2 parents e5ef3e2 + 1b9d21b
commit f9ad529
Show file tree

Hide file tree

Showing 8 changed files with 3,296 additions and 0 deletions.
diff --git a/docs/make.jl b/docs/make.jl
@@ -89,6 +89,7 @@ for name in (
     "transport_supg",
     "heat_equation_sphere",
     "transport_hypersurface",
+    "stokes_flow",
 )
     gen_markdown_with_literate(joinpath(example_src, name), "$(name).jl", example_dir)
 end
@@ -119,6 +120,7 @@ makedocs(;
             "example/heat_equation_sphere.md",
             "example/heat_equation_two_layers.md",
             "example/transport_hypersurface.md",
+            "example/stokes_flow.md",
         ],
     ],
     # remotes = nothing, # tmp fix for bmxam windows

diff --git a/docs/src/assets/Stokes_flow_Moffat_vortices.png b/docs/src/assets/Stokes_flow_Moffat_vortices.png
diff --git a/docs/src/assets/Stokes_flow_unsteady.gif b/docs/src/assets/Stokes_flow_unsteady.gif
diff --git a/docs/src/assets/Stokes_steady_Moffat_vortices_1.png b/docs/src/assets/Stokes_steady_Moffat_vortices_1.png
diff --git a/input/mesh/domainTriangle_tri.geo b/input/mesh/domainTriangle_tri.geo
@@ -0,0 +1,34 @@
+//+
+//L = 10.0;
+//a = 3.25;
+//b = 13.5;
+//x1 = L*a/b;
+//y1 = Sqrt(L*L-x1*x1);
+y1 = 1.0;
+x1 = Tan(Pi*14./180.);
+Point(1) = {0, 0, 0, 1.0};
+//+
+Point(2) = {x1, y1, 0, 1.0};
+//+
+Point(3) = {-x1, y1, 0, 1.0};
+//+
+Line(1) = {1, 2};
+//+
+Line(2) = {2, 3};
+//+
+Line(3) = {3, 1};
+//+
+Curve Loop(1) = {1, 2, 3};
+//+
+Plane Surface(1) = {1};
+//+
+Physical Curve("edge_right") = {1};
+//+
+Physical Curve("edge_top") = {2};
+//+
+Physical Curve("edge_left") = {3};
+//+
+Physical Surface("Domain") = {1};
+//+
+Transfinite Curve {1, 3} = 51 Using Progression 1;
+Transfinite Curve {2} = 21 Using Progression 1;