Skip to content
Snippets Groups Projects
Commit 7db337b4 authored by Thomas Fetzer's avatar Thomas Fetzer
Browse files

[test][rans] Remove old PDELab-implementation based comparison files

parent d0755697
No related branches found
No related tags found
1 merge request!1110Freeflow/fix eddy thermal conductivity
add_input_file_links()
dune_symlink_to_source_files(FILES laufer_re50000.csv laufer_re50000_u+y+.csv
pdelab-kepsilon_twolayers.csv pdelab-kepsilon_wallfunction.csv
pdelab-komega.csv pdelab-lowrekepsilon.csv pdelab-zeroeq.csv pdelab-oneeq.csv)
dune_symlink_to_source_files(FILES laufer_re50000.csv laufer_re50000_u+y+.csv)
add_executable(test_pipe_laufer_kepsilon EXCLUDE_FROM_ALL test_pipe_laufer.cc)
target_compile_definitions(test_pipe_laufer_kepsilon PUBLIC "KEPSILON=1")
......
source diff could not be displayed: it is too large. Options to address this: view the blob.
source diff could not be displayed: it is too large. Options to address this: view the blob.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
......@@ -241,22 +241,16 @@ int main(int argc, char** argv) try
gnuplot_lawOfTheWall.setOption("set xrange [1:3000]");
gnuplot_lawOfTheWall.addFileToPlot("laufer_re50000_u+y+.csv", "u 1:2 w p t 'Laufer 1954, Re=50000'");
#if LOWREKEPSILON
gnuplot_lawOfTheWall.addFileToPlot("pdelab-lowrekepsilon.csv", "u 23:22 w l lw 2 t 'PDELab Low-Re k-epsilon'");
gnuplot_lawOfTheWall.addFileToPlot(std::string(fileName) + ".csv", "u 12:13 w l lc 7");
#elif KEPSILON
gnuplot_lawOfTheWall.addFileToPlot("pdelab-kepsilon_twolayers.csv", "u 34:33 w l lw 2 t 'PDELab k-epsilon (two layers)'");
gnuplot_lawOfTheWall.addFileToPlot("pdelab-kepsilon_wallfunction.csv", "u 36:35 w l lw 2 t 'PDELab k-epsilon (wall function)'");
gnuplot_lawOfTheWall.addFileToPlot(std::string(fileName) + ".csv", "u 12:13 w l lc 7 t 'with u_{tau}'");
gnuplot_lawOfTheWall.addFileToPlot(std::string(fileName) + ".csv", "u 16:17 w l lc 8 t 'with u_{tau,nom}'");
gnuplot_lawOfTheWall.addFileToPlot(std::string(fileName) + ".csv", "u 12:13 w l lc 7");
#elif KOMEGA
gnuplot_lawOfTheWall.addFileToPlot("pdelab-komega.csv", "u 24:23 w l lw 2 t 'PDELab k-omega'");
gnuplot_lawOfTheWall.addFileToPlot(std::string(fileName) + ".csv", "u 12:13 w l lc 7");
#elif ONEEQ
gnuplot_lawOfTheWall.addFileToPlot("pdelab-oneeq.csv", "u 22:21 w l lw 2 t 'PDELab 1-Eq.'");
gnuplot_lawOfTheWall.addFileToPlot(std::string(fileName) + ".csv", "u 12:13 w l lc 7");
#else
gnuplot_lawOfTheWall.addFileToPlot("pdelab-zeroeq.csv", "u 22:21 w l lw 2 t 'PDELab 0-Eq.'");
gnuplot_lawOfTheWall.addFileToPlot(std::string(fileName) + ".csv", "u 12:13 w l lc 7");
#endif
gnuplot_lawOfTheWall.plot(std::string(gnuplotFileName));
......@@ -270,20 +264,14 @@ int main(int argc, char** argv) try
gnuplot_velocityProfile.setYlabel("y [-]");
gnuplot_velocityProfile.addFileToPlot("laufer_re50000.csv", "u 2:1 w p t 'Laufer 1954, Re=50000'");
#if LOWREKEPSILON
gnuplot_velocityProfile.addFileToPlot("pdelab-lowrekepsilon.csv", "u 5:($27/0.2456) w l lw 2 t 'PDELab Low-Re k-epsilon'");
gnuplot_velocityProfile.addFileToPlot(std::string(fileName) + ".csv", "u 7:($26/0.2456) w l lc 7");
#elif KEPSILON
gnuplot_velocityProfile.addFileToPlot("pdelab-kepsilon_twolayers.csv", "u ($5):($40/0.2469) w l lw 2 t 'PDELab k-epsilon (two layers)'");
gnuplot_velocityProfile.addFileToPlot("pdelab-kepsilon_wallfunction.csv", "u ($5):($40/0.2469) w l lw 2 t 'PDELab k-epsilon (wall function)'");
gnuplot_velocityProfile.addFileToPlot(std::string(fileName) + ".csv", "u 7:($28/0.2456) w l lc 7");
#elif KOMEGA
gnuplot_velocityProfile.addFileToPlot("pdelab-komega.csv", "u 5:($29/0.2456) w l lw 2 t 'PDELab k-omega'");
gnuplot_velocityProfile.addFileToPlot(std::string(fileName) + ".csv", "u 7:($26/0.2456) w l lc 7");
#elif ONEEQ
gnuplot_velocityProfile.addFileToPlot("pdelab-oneeq.csv", "u 5:($26/0.2456) w l lw 2 t 'PDELab 1-Eq.'");
gnuplot_velocityProfile.addFileToPlot(std::string(fileName) + ".csv", "u 7:($25/0.2456) w l lc 7");
#else
gnuplot_velocityProfile.addFileToPlot("pdelab-zeroeq.csv", "u 5:($26/0.2456) w l lw 2 t 'PDELab 0-Eq.'");
gnuplot_velocityProfile.addFileToPlot(std::string(fileName) + ".csv", "u 7:($24/0.2456) w l lc 7");
#endif
gnuplot_velocityProfile.plot(std::string(gnuplotFileName));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment