wallFunctions submodule updated (#217)

* wallFunctions submodule updated

* update name wallfunction lib

Author: CyrilleBonamy

TurbulenceModels/wallFunctions

@@ -54,7 +54,7 @@ maxAlphaCo      5;
 maxDeltaT       1e0;
 
 libs (
-     "libWilcoxOmegaWallFunction.so"
+     "libroughWallFunctions.so"
      );
 
 // ************************************************************************* //

test-ci/1DBoundaryLayer/system/controlDict

@@ -5,8 +5,9 @@
 import fluidfoam
 
 
-def rms(x):
-    return np.sqrt(x.dot(x) / x.size)
+def rms(y, x):
+    u = np.trapz(y**2, x)
+    return np.sqrt(u / np.abs(x[-1]-x[0]))
 
 
 #
@@ -95,11 +96,11 @@ def rms(x):
 # zpprodDNS    = npzfiles['arr_9']
 
 u_interp = np.interp(zuw, z[:] / H, U[0, :] / Umax)
-rms_u = rms(u_interp - uw)
+rms_u = rms(u_interp[:-1] - uw[:-1], zuw[:-1])
 print("rms_u=", rms_u)
 assert rms_u <= 0.05
 
 k_interp = np.interp(zkw, z[:] / H, k[:] / utau**2)
-rms_k = rms(k_interp - kw)
+rms_k = rms(k_interp[:-1] - kw[:-1], zkw[:-1])
 print("rms_k=", rms_k)
 assert rms_k <= 0.16

test-ci/1DBoundaryLayer/test_1DBoundaryLayer.py

@@ -54,7 +54,7 @@ maxAlphaCo      5;
 maxDeltaT       1e0;
 
 libs (
-     "libWilcoxOmegaWallFunction.so"
+     "libroughWallFunctions.so"
      );
 
 // ************************************************************************* //

tutorials/RAS/1DBoundaryLayer/system/controlDict

@@ -56,7 +56,7 @@ maxDeltaT       5e-5;
 libs (
      "libgroovyBC.so"
      "libgroovyStandardBCs.so"
-     "libWilcoxOmegaWallFunction.so"
+     "libroughWallFunctions.so"
      );
 
 // ************************************************************************* //

tutorials/RAS/1DOscill/system/controlDict

@@ -56,7 +56,7 @@ maxDeltaT       1e-3;
 libs (
       "libgroovyBC.so"
       "libgroovyStandardBCs.so"
-      "libWilcoxOmegaWallFunction.so"
+      "libroughWallFunctions.so"
      );
 
 // ************************************************************************* //