Update the Spherical Kernel Toy model

Verified that this gives consistent results on hyperion for 16 / 32 / 48 / 96 cores and, just in case, 53. Not sure if the variations should be considered significant, to be honest, but there is nothing especially suspicious in the solutions. Grid has 480000 tetrahedral elements, P2 for velocity, P1 for pressure

16: Pressure (C1): (83577, 4.501686375938372e-05, -0.7372135020091626, 0.7297416090999129, 3.762374422418013, 40.86046689261976,  np.float64(0.14133832355990378))

32: Pressure (C1): (83577, 0.00011456943015901214, -0.7382004092043332, 0.7033917668241267, 9.575369264399757, 40.878594801302334, np.float64(0.1414010288694044))

96: Pressure (C1): (83577, 0.0003194554969466519, -0.7394435695222795, 0.7056720965038505, 26.699132068310327, 41.07756585313935, np.float64(0.1420892793237494))

53: Pressure (C1): (83577, 8.779219727421255e-05, -0.7371671638813321, 0.6970817607318897, 7.337408471586862, 40.92046307139055, np.float64(0.14154585323277188))

Author: lmoresi

Notebooks/Examples-Stokes_Kernels/Ex_Stokes_Spherical_Kernels.py

@@ -227,7 +227,7 @@
 # Optional - read in a different solution if you want to just use this for
 # visualisation - set the `expt_name` here
 
-expt_name="Spherical_Kernel_np_7"
+expt_name="Spherical_Kernel_np_96"
 
 import os
 
@@ -258,7 +258,8 @@
     pvmesh.point_data["V"] = vis.vector_fn_to_pv_points(pvmesh, v_soln_ckpt.sym)
     pvmesh.point_data["R"] = vis.scalar_fn_to_pv_points(pvmesh, radius_fn)
 
-    pvmesh_c = pvmesh.clip_scalar(scalars="R", value=0.95).clip()
+    # pvmesh_c = pvmesh.clip_scalar(scalars="R", value=0.95).clip()
+    pvmesh_c = pvmesh.clip(normal="z")
 
     skip = 10
     points = np.zeros((meshball._centroids[::skip].shape[0], 3))
@@ -279,14 +280,14 @@
     pl = pv.Plotter(window_size=(750, 750))
 
     pl.add_mesh(
-        pvmesh,
+        pvmesh_c,
         cmap="coolwarm",
-        edge_color="Grey",
-        edge_opacity=0.33,
+        edge_color="Black",
+        # edge_opacity=0.33,
         scalars="P",
         show_edges=True,
         use_transparency=False,
-        opacity=0.9,
+        opacity=1.0,
         show_scalar_bar=False
     )
 
@@ -296,16 +297,16 @@
     vsol_rms = np.sqrt(velocity_points.point_data["V"][:, 0] ** 2 + velocity_points.point_data["V"][:, 1] ** 2).mean()
     # print(vsol_rms)
 
-    pl.export_html("stokes_sphere_plot.html")
+    pl.export_html("stokes_sphere_plot_96.html")
     # pl.show(cpos="xy", jupyter_backend="trame")
 
 
 # +
 from IPython.display import IFrame
-IFrame(src="./stokes_sphere_plot.html", width=750, height=750)
+IFrame(src="./stokes_sphere_plot_96.html", width=750, height=750)
 
 
 # -
 
 
-
+p_soln.stats()