Finished test case, added reference data for a few typical schemes and plotting thereof

Author: UnnamedMoose

OFtutorial15_discretisation/OFtutorial15.C

@@ -26,6 +26,35 @@ License
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
+
+/*
+scalarTransportFoam solver solves this:
+
+    while (simple.loop(runTime))
+    {
+        Info<< "Time = " << runTime.timeName() << nl << endl;
+
+        while (simple.correctNonOrthogonal())
+        {
+            fvScalarMatrix TEqn
+            (
+                fvm::ddt(T)
+              + fvm::div(phi, T)
+              - fvm::laplacian(DT, T)
+             ==
+                fvOptions(T)
+            );
+
+            TEqn.relax();
+            fvOptions.constrain(TEqn);
+            TEqn.solve();
+            fvOptions.correct(T);
+        }
+
+        runTime.write();
+    }
+*/
+
 int main(int argc, char *argv[])
 {
 	// Set up the case, parse command line options and create the grid

OFtutorial15_discretisation/testCase/Allclean

@@ -1,2 +1,2 @@
 #!/bin/bash
-rm -rf constant/polyMesh/* log.* 0.* 1* 2*
+rm -rf constant/polyMesh/* log.* 0.* 1* 2* postProcessing

OFtutorial15_discretisation/testCase/plotData.py

@@ -0,0 +1,35 @@
+#!/usr/bin/env python3
+# -*- coding: utf-8 -*-
+"""
+Created on Sun Jun  6 11:04:35 2021
+
+@author: artur
+"""
+
+import pandas
+import numpy as np
+import matplotlib
+import matplotlib.pyplot as plt
+import os
+
+font = {'family' : 'serif',
+        'weight' : 'normal',
+        'size'   : 16}
+
+matplotlib.rc('font', **font)
+
+# %% Read the data.
+refData = {f.split("_")[1].replace(".csv", ""): pandas.read_csv(os.path.join("refData", f))
+           for f in os.listdir("./refData")}
+
+# %% Plot solutions obtained using each scheme.
+fig, ax = plt.subplots(1)
+ax.set_xlabel("x [m]")
+ax.set_ylabel(r"$\phi$ [-]")
+for s in refData:
+    ax.plot(refData[s]["x"], refData[s]["T"], "o--", lw=1, ms=4, mew=2, label=s)
+xlim = ax.get_xlim()
+ax.hlines([0, 1], xlim[0], xlim[1], color="k", ls="--", alpha=0.5, zorder=-10)
+ax.set_xlim(xlim)
+ax.legend(loc="lower left", framealpha=1)
+plt.tight_layout()

OFtutorial15_discretisation/testCase/refData/centrelineData_linear.csv

@@ -0,0 +1,21 @@
+x,y,z,T
+0.025,0.05,0.05,1.00021
+0.075,0.05,0.05,0.999031
+0.125,0.05,0.05,0.99892
+0.175,0.05,0.05,1.00479
+0.225,0.05,0.05,1.01202
+0.275,0.05,0.05,1.00212
+0.325,0.05,0.05,0.953314
+0.375,0.05,0.05,0.855654
+0.425,0.05,0.05,0.717507
+0.475,0.05,0.05,0.560317
+0.525,0.05,0.05,0.407872
+0.575,0.05,0.05,0.277604
+0.625,0.05,0.05,0.177362
+0.675,0.05,0.05,0.106822
+0.725,0.05,0.05,0.060899
+0.775,0.05,0.05,0.0329906
+0.825,0.05,0.05,0.0170445
+0.875,0.05,0.05,0.00841831
+0.925,0.05,0.05,0.00402652
+0.975,0.05,0.05,0.00165676

OFtutorial15_discretisation/testCase/refData/centrelineData_linearUpwind.csv

@@ -0,0 +1,21 @@
+x,y,z,T
+0.025,0.05,0.05,0.999996
+0.075,0.05,0.05,1.00001
+0.125,0.05,0.05,1.00022
+0.175,0.05,0.05,1.0004
+0.225,0.05,0.05,0.997676
+0.275,0.05,0.05,0.983597
+0.325,0.05,0.05,0.94426
+0.375,0.05,0.05,0.866203
+0.425,0.05,0.05,0.745505
+0.475,0.05,0.05,0.593092
+0.525,0.05,0.05,0.43145
+0.575,0.05,0.05,0.284822
+0.625,0.05,0.05,0.169655
+0.675,0.05,0.05,0.0907059
+0.725,0.05,0.05,0.0432548
+0.775,0.05,0.05,0.0182241
+0.825,0.05,0.05,0.00667237
+0.875,0.05,0.05,0.00205313
+0.925,0.05,0.05,0.000487457
+0.975,0.05,0.05,4.4752e-05

OFtutorial15_discretisation/testCase/refData/centrelineData_upwind.csv

@@ -0,0 +1,21 @@
+x,y,z,T
+0.025,0.05,0.05,0.999308
+0.075,0.05,0.05,0.996427
+0.125,0.05,0.05,0.989001
+0.175,0.05,0.05,0.973398
+0.225,0.05,0.05,0.945365
+0.275,0.05,0.05,0.901116
+0.325,0.05,0.05,0.838635
+0.375,0.05,0.05,0.758637
+0.425,0.05,0.05,0.664766
+0.475,0.05,0.05,0.562933
+0.525,0.05,0.05,0.460049
+0.575,0.05,0.05,0.362635
+0.625,0.05,0.05,0.275731
+0.675,0.05,0.05,0.202337
+0.725,0.05,0.05,0.143416
+0.775,0.05,0.05,0.0982881
+0.825,0.05,0.05,0.0652061
+0.875,0.05,0.05,0.041932
+0.925,0.05,0.05,0.0262306
+0.975,0.05,0.05,0.0165926

OFtutorial15_discretisation/testCase/system/controlDict

@@ -23,13 +23,13 @@ startTime       0;
 
 stopAt          endTime;
 
-endTime         1.5;
+endTime         0.99;
 
 deltaT          0.01;
 
 writeControl    timeStep;
 
-writeInterval   10;
+writeInterval   50;
 
 purgeWrite      0;
 
@@ -45,6 +45,31 @@ timePrecision   6;
 
 runTimeModifiable true;
 
+functions
+{
+    samplingLine
+    {
+        type            sets;
+        libs            ("libsampling.so");
+        writeControl    writeTime;
+
+        fields          (T);
+        interpolationScheme cellPoint;
+        setFormat   csv;
+
+        sets
+        (
+            centrelineData
+            {
+                type    lineCell;
+                axis    xyz;
+                nPoints 100;
+                start       (0 0.05 0.05);
+                end         (1 0.05 0.05);
+            }
+        );
+    }
+}
 
 // ************************************************************************* //
 

OFtutorial15_discretisation/testCase/system/fvSchemes

@@ -28,7 +28,9 @@ gradSchemes
 divSchemes
 {
     default         none;
-    div(phi,T)      Gauss linearUpwind grad(T);
+//    div(phi,T)      Gauss linearUpwind grad(T);
+//    div(phi,T)      Gauss upwind;
+    div(phi,T)      Gauss linear;
 }
 
 laplacianSchemes