#r "BoSSSpad.dll"
using System;
using System.Collections.Generic;
using System.Linq;
using System.IO;
using System.Data;
using System.Globalization;
using System.Threading;
using ilPSP;
using ilPSP.Utils;
using BoSSS.Platform;
using BoSSS.Foundation;
using BoSSS.Foundation.Grid;
using BoSSS.Foundation.Grid.Classic;
using BoSSS.Foundation.IO;
using BoSSS.Solution;
using BoSSS.Solution.Control;
using BoSSS.Solution.GridImport;
using BoSSS.Solution.Statistic;
using BoSSS.Solution.Utils;
using BoSSS.Solution.Gnuplot;
using BoSSS.Application.BoSSSpad;
using BoSSS.Application.XNSE_Solver;
using static BoSSS.Application.BoSSSpad.BoSSSshell;
using BoSSS.Foundation.Grid.RefElements;
using BoSSS.Platform.LinAlg;
using BoSSS.Application.XNSE_Solver.PhysicalBasedTestcases.PrintingNip;
Init();
Only if the output directories were previously created via Part0_PrintingNip_Setup output is stored to the subdirectories.
Otherwise, datatables are stored directly to the working directory and plots displayed in the notebook.
The plots are not adjusted to look "nice" in the notebook.
This is not part of the ValidationTestRunner.
static bool GenerateOutput = Directory.Exists(@"./PrintingNip") && Directory.Exists(@"./PrintingNip/Figures") && Directory.Exists(@"./PrintingNip/Files") && Directory.Exists(@"./PrintingNip/Output");
string ProjectName = "PrintingNip_Part2";
BoSSSshell.WorkflowMgm.Init(ProjectName);
Project name is set to 'PrintingNip_Part2'. Default Execution queue is chosen for the database. Opening existing database '\\fdygitrunner\ValidationTests\databases\PrintingNip_Part2'.
static var myDb = BoSSSshell.WorkflowMgm.DefaultDatabase;
var sessions = myDb.Sessions.Where(s => s.ProjectName == ProjectName).ToArray();
sessions.Count()
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Clear();
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("PressureRange", Postprocessing.PressureRange);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("VelocityXMax", Postprocessing.VelocityXMax);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("Massflux", Postprocessing.Massflux);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("NipShearRate", Postprocessing.NipShearRate);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("NipShearRateCylinder", Postprocessing.NipShearRateCylinder);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("ShearStressCylinder", Postprocessing.ShearStressCylinder);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("NipShearRateSubstrate", Postprocessing.NipShearRateSubstrate);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("ShearStressSubstrate", Postprocessing.ShearStressSubstrate);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("ViscousDissipation", Postprocessing.ViscousDissipation);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("PositionOfSynchronousFlow", Postprocessing.PositionOfSynchronousFlow);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("PositionOfStagnatingFlow", Postprocessing.PositionOfStagnatingFlow);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("dPdXatSynchronousPoint", Postprocessing.dPdXatSynchronousPoint);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("dPdXatStagnationPoint", Postprocessing.dPdXatStagnationPoint);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("dPdXatNip", Postprocessing.dPdXatNip);
BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Add("dPdXatConstantX", Postprocessing.dPdXatConstantX);
DataTable tab;
if(GenerateOutput){
if(!File.Exists("./PrintingNip/Files/PrintingNip-Part2Export.csv")){
tab = sessions.GetSessionTable(BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Select(kv => new Tuple<string, Func<ISessionInfo, object>>(kv.Key, kv.Value)).ToArray());
tab.ToCSVFile("./PrintingNip/Files/PrintingNip-Part2Export.csv", ColSep: '%');
} else{
tab = TableExtensions.FromCSVFile("./PrintingNip/Files/PrintingNip-Part2Export.csv", ColSep: '%');
}
} else {
// in this case always recreate the datatable
tab = sessions.GetSessionTable(BoSSSshell.WorkflowMgm.AdditionalSessionTableColums.Select(kv => new Tuple<string, Func<ISessionInfo, object>>(kv.Key, kv.Value)).ToArray());
tab.ToCSVFile("./PrintingNip-Part2Export.csv", ColSep: '%');
}
Loaded Fields for J2000_delta1E-06_V10_xm0.01 Loaded Fields for J2000_delta1E-06_V5_xm0.01 Loaded Fields for J2000_delta1E-06_V1_xm0.01 Loaded Fields for J2000_delta1E-06_V0.5_xm0.01 Loaded Fields for J2000_delta1E-06_V0.1_xm0.01 Loaded Fields for J2000_delta1E-05_V10_xm0.01 Loaded Fields for J2000_delta1E-05_V5_xm0.01 Loaded Fields for J2000_delta1E-05_V1_xm0.01 Loaded Fields for J2000_delta1E-05_V0.5_xm0.01 Loaded Fields for J2000_delta1E-05_V0.1_xm0.01 Loaded Fields for J2000_delta0.0001_V10_xm0.01 Loaded Fields for J2000_delta0.0001_V5_xm0.01 Loaded Fields for J2000_delta0.0001_V1_xm0.01 Loaded Fields for J2000_delta0.0001_V0.5_xm0.01 Loaded Fields for J2000_delta0.0001_V0.1_xm0.01 Loaded Fields for J2000_delta0.001_V10_xm0.01 Loaded Fields for J2000_delta0.001_V5_xm0.01 Loaded Fields for J2000_delta0.001_V1_xm0.01 Loaded Fields for J2000_delta0.001_V0.5_xm0.01 Loaded Fields for J2000_delta0.001_V0.1_xm0.01
string[] Columns = new string[] {"PressureRange", "VelocityXMax", "Massflux", "NipShearRate", "NipShearRateCylinder", "ShearStressCylinder", "NipShearRateSubstrate", "ShearStressSubstrate", "ViscousDissipation", "PositionOfSynchronousFlow", "PositionOfStagnatingFlow", "dPdXatSynchronousPoint", "dPdXatStagnationPoint", "dPdXatNip", "dPdXatConstantX"};
string[] Labels = new string[] {"Pressure", "Velocity", "Massflux", "Shearrate (Nip)", "Shearrate (Cylinder)", "Shearforce (Cylinder)", "Shearrate (Substrate)", "Shearforce (Substrate)", "Dissipation", "$x_{sync}$", "$x_{stag}$", @"$\\\frac{\\\partial p}{\\\partial x}_{sync}$", @"$\\\frac{\\\partial p}{\\\partial x}_{stag}$", @"$\\\frac{\\\partial p}{\\\partial x}_{nip}$", @"$\\\frac{\\\partial p}{\\\partial x}$"};
string[] Units = new string[] {@"$\\\left[\\\frac{N}{m^2}\\right]$", @"$\\\left[\\\frac{m}{s}\\right]$", @"$\\\left[\\\frac{kg}{s}\\right]$", @"$\\\left[\\\frac{1}{s}\\right]$", @"$\\\left[\\\frac{1}{s}\\right]$", @"$\\\left[N\\right]$", @"$\\\left[\\\frac{1}{s}\\right]$", @"$\\\left[N\\right]$", @"$\\\left[W\\right]$", @"$\\\left[m\\right]$", @"$\\\left[m\\right]$", @"$\\\left[\\\frac{N}{m^3}\\right]$", @"$\\\left[\\\frac{N}{m^3}\\right]$", @"$\\\left[\\\frac{N}{m^3}\\right]$", @"$\\\left[\\\frac{N}{m^3}\\right]$"}; // some escape character nonsense
var tab2export = tab.ExtractColumns(Columns.ToList().Prepend("id:V_Wall").Prepend("id:P_Diff").Prepend("id:Radius").Prepend("id:delta").Prepend("SessionName").ToArray());
if(GenerateOutput){
tab2export.ToHTMLFile("PrintingNip-Part2.html", Path.GetFullPath("./PrintingNip/Files"));
tab2export.ToCSVFile("./PrintingNip/Files/PrintingNip-Part2.csv", ColSep: ';');
} else {
tab2export.ToHTMLFile("PrintingNip-Part2.html", Path.GetFullPath("./"));
tab2export.ToCSVFile("./PrintingNip-Part2.csv", ColSep: ';');
}
Plot all Properties, we looked at before as functions of nip width. Remember, this time around the pressure b.c. is such that the stagnation point should remain constant
var tab2plt = tab.ExtractRows((i,row) => Convert.ToDouble(row["id:P_Diff"]) > 0.0);
List<Plot2Ddata> plts = new List<Plot2Ddata>();
foreach(string col in Columns){
int k = Columns.IndexOf(col);
var plt = tab2plt.ToPlot("id:delta", col, ColName_GroupSelection: new string[] { "id:V_Wall" });
plt.LogX = true;
plt.LogY = true;
plt.Xlabel = @"Nip Width $\\\left[m\\right]$";
plt.Ylabel = Labels[k] + " " + Units[k];
plt.LabelTitleFont = 32;
plt.LabelFont = 24;
plt.LegendFont = 18;
plt.lmargin = 12;
plt.bmargin = 4;
plt.ShowLegend = true;
plt.LegendBox = true;
if(plt.Regression().Average(kv => kv.Value) > 0){
plt.LegendAlignment = new string[]{"i", "l", "t"};
plt.LegendSwap = true;
} else{
plt.LegendAlignment = new string[]{"i", "r", "t"};
plt.LegendSwap = false;
}
for(int n= 0; n < plt.dataGroups.Count(); n++){
plt.dataGroups[n].Format.LineWidth = 2;
plt.dataGroups[n].Format.PointSize = 1;
plt.dataGroups[n].Name = "$V_W = " + Convert.ToDouble(plt.dataGroups[n].Name.Split("id:V_Wall").Last()).ToString("N1") + @" \\\frac{m}{s}$";
}
plt.ModFormat();
plts.Add(plt);
}
int k = 0;
foreach(var plt in plts){
var gp = new Gnuplot();
gp.PlotLogSlope(plt, format: new PlotFormat(lineColor: (LineColors)7));
gp.Cmd("set key box width -5");
plt.ToGnuplot(gp);
if(GenerateOutput){
gp.PlotCairolatex().SaveTo("./PrintingNip/Figures/Part2/"+Columns[k]+".tex");
} else {
gp.PlotNow().Display();
}
k++;
}
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!
Using gnuplot: C:\Program Files (x86)\FDY\BoSSS\bin\native\win\gnuplot-gp510-20160418-win32-mingw\gnuplot\bin\gnuplot.exe Note: In a Jupyter Worksheet, you must NOT have a trailing semicolon in order to see the plot on screen; otherwise, the output migth be surpressed.!