Commit a4c9e32b authored by Ned Coltman's avatar Ned Coltman
Browse files

[output] add scripts for plotting and reference comparison

parent 6e2462b2
rm *.png
gnuplot scripts/coefficientoffriction.gp
gnuplot scripts/velocitydistribution_m4.gp
gnuplot scripts/velocitydistribution_1.gp
gnuplot scripts/velocitydistribution_4.gp
gnuplot scripts/velocitydistribution_6.gp
gnuplot scripts/velocitydistribution_10.gp
eog *.png
# Exp data from Driver and Seegmiller, bottom wall
# x means x/H
#variables="x","cf","error"
-3.956, 2.88e-3, 0.202e-3
-1.804, 2.85e-3, 0.171e-3
-0.804, 3.11e-3, 0.311e-3
0.4840, 0.09e-3, 0.014e-3
1.8040, -0.11e-3, 0.019e-3
2.8040, -0.52e-3, 0.109e-3
3.8040, -1.02e-3, 0.163e-3
4.8040, -0.82e-3, 0.123e-3
5.8820, -0.22e-3, 0.062e-3
7.0900, 0.45e-3, 0.167e-3
7.9780, 1.02e-3, 0.122e-3
9.1960, 1.18e-3, 0.106e-3
11.196, 1.54e-3, 0.123e-3
13.196, 1.73e-3, 0.121e-3
16.196, 1.91e-3, 0.134e-3
20.196, 2.02e-3, 0.141e-3
24.136, 2.04e-3, 0.143e-3
28.130, 2.05e-3, 0.205e-3
32.196, 2.12e-3, 0.148e-3
35.994, 2.30e-3, 0.184e-3
# Exp data drom Driver and Seegmiller,
# y is y/H, u is u/Uref, v is v/Uref,
#variables="XID","y","u", "XID","y","u", "XID","y","u", "XID","y","u", "XID","y","u"
-4.0, 1.10, 0.657, 1.0, 0.10, -0.008, 4.0, 0.10, -0.171, 6.0, 0.10, 0.147, 10.0, 0.10, 0.418,
-4.0, 1.15, 0.696, 1.0, 0.15, -0.017, 4.0, 0.15, -0.154, 6.0, 0.15, 0.169, 10.0, 0.15, 0.433,
-4.0, 1.20, 0.719, 1.0, 0.20, -0.027, 4.0, 0.20, -0.134, 6.0, 0.20, 0.188, 10.0, 0.20, 0.447,
-4.0, 1.30, 0.760, 1.0, 0.30, -0.035, 4.0, 0.30, -0.057, 6.0, 0.30, 0.248, 10.0, 0.30, 0.471,
-4.0, 1.40, 0.790, 1.0, 0.40, -0.035, 4.0, 0.40, 0.052, 6.0, 0.40, 0.306, 10.0, 0.40, 0.496,
-4.0, 1.50, 0.818, 1.0, 0.50, -0.035, 4.0, 0.50, 0.169, 6.0, 0.50, 0.369, 10.0, 0.50, 0.526,
-4.0, 1.70, 0.870, 1.0, 0.60, -0.027, 4.0, 0.60, 0.279, 6.0, 0.60, 0.435, 10.0, 0.60, 0.560,
-4.0, 2.00, 0.926, 1.0, 0.70, -0.005, 4.0, 0.70, 0.381, 6.0, 0.70, 0.501, 10.0, 0.70, 0.596,
-4.0, 2.40, 0.982, 1.0, 0.80, 0.081, 4.0, 0.80, 0.473, 6.0, 0.80, 0.569, 10.0, 0.80, 0.638,
-4.0, 2.80, 1.003, 1.0, 0.90, 0.292, 4.0, 0.90, 0.568, 6.0, 0.90, 0.632, 10.0, 0.90, 0.676,
-4.0, 3.20, 1.003, 1.0, 1.00, 0.519, 4.0, 1.00, 0.656, 6.0, 1.00, 0.686, 10.0, 1.00, 0.712,
-4.0, 3.60, 1.000, 1.0, 1.10, 0.704, 4.0, 1.10, 0.730, 6.0, 1.10, 0.730, 10.0, 1.10, 0.745,
-4.0, 4.00, 1.001, 1.0, 1.15, 0.740, 4.0, 1.15, 0.754, 6.0, 1.15, 0.747, 10.0, 1.15, 0.758,
-4.0, 5.00, 1.001, 1.0, 1.20, 0.762, 4.0, 1.20, 0.778, 6.0, 1.20, 0.763, 10.0, 1.20, 0.772,
-4.0, 6.00, 1.000, 1.0, 1.30, 0.793, 4.0, 1.30, 0.811, 6.0, 1.30, 0.794, 10.0, 1.30, 0.798,
-4.0, 7.50, 1.001, 1.0, 1.40, 0.827, 4.0, 1.40, 0.835, 6.0, 1.40, 0.815, 10.0, 1.40, 0.819,
-4.0, 8.20, 0.943, 1.0, 1.50, 0.851, 4.0, 1.50, 0.856, 6.0, 1.50, 0.834, 10.0, 1.50, 0.839,
-4.0, 8.20, 0.943, 1.0, 1.70, 0.891, 4.0, 1.70, 0.895, 6.0, 1.70, 0.869, 10.0, 1.70, 0.871,
-4.0, 8.20, 0.943, 1.0, 2.00, 0.942, 4.0, 2.00, 0.940, 6.0, 2.00, 0.915, 10.0, 2.00, 0.908,
-4.0, 8.20, 0.943, 1.0, 2.40, 0.990, 4.0, 2.40, 0.978, 6.0, 2.40, 0.949, 10.0, 2.40, 0.929,
-4.0, 8.20, 0.943, 1.0, 2.80, 1.009, 4.0, 2.80, 0.986, 6.0, 2.80, 0.954, 10.0, 2.80, 0.933,
-4.0, 8.20, 0.943, 1.0, 3.20, 1.003, 4.0, 3.20, 0.984, 6.0, 3.20, 0.956, 10.0, 3.20, 0.937,
-4.0, 8.20, 0.943, 1.0, 3.60, 1.002, 4.0, 3.60, 0.985, 6.0, 3.60, 0.959, 10.0, 3.60, 0.939,
-4.0, 8.20, 0.943, 1.0, 4.00, 1.000, 4.0, 4.00, 0.985, 6.0, 4.00, 0.961, 10.0, 4.00, 0.940,
-4.0, 8.20, 0.943, 1.0, 5.00, 0.997, 4.0, 5.00, 0.983, 6.0, 5.00, 0.963, 10.0, 5.00, 0.944,
-4.0, 8.20, 0.943, 1.0, 6.00, 0.994, 4.0, 6.00, 0.980, 6.0, 6.00, 0.965, 10.0, 6.00, 0.947,
-4.0, 8.20, 0.943, 1.0, 7.50, 0.993, 4.0, 7.50, 0.980, 6.0, 7.50, 0.969, 10.0, 7.50, 0.949,
-4.0, 8.20, 0.943, 1.0, 8.20, 0.919, 4.0, 8.20, 0.892, 6.0, 8.20, 0.870, 10.0, 8.20, 0.827,
reset
set term pngcairo size 800,600 solid
set output "./FrictionCoefficient.png"
set datafile separator ','
DATA='./'
set xlabel "distance from the step [-]"
set ylabel "Coefficient of Friction at Base Wall"
set yrange [-0.15e-2:0.25e-2]
set xrange [0:40]
set key right bottom samplen 1
plot './references/bfs_experimentaldata_friction.csv' u 1:2 w p t 'DriverSeegmiller Experimental Data', \
'./references/bfs_referencemodel_friction.csv' u 1:2 w l lc rgb "black" t 'CFL3D Comparison Model',\
DATA.'backwardfacingstep_friction_kepsilon.csv' u 24:(($15+$16)/($5*0.5*44.2*44.2)) w l lc rgb "red" t 'KEpsilon with u_{tau}',\
DATA.'backwardfacingstep_friction_lowrekepsilon.csv' u 24:(($15+$16)/($5*0.5*44.2*44.2)) w l lc rgb "yellow" t 'Low Re KEpsilon',\
DATA.'backwardfacingstep_friction_komega.csv' u 24:(($15+$16)/($5*0.5*44.2*44.2)) w l lc rgb "green" t 'KOmega',\
DATA.'backwardfacingstep_friction_oneeq.csv' u 24:(($15+$16)/($5*0.5*44.2*44.2)) w l lc rgb "cyan" t 'OneEq',\
DATA.'backwardfacingstep_friction_zeroeq.csv' u 24:(($15+$16)/($5*0.5*44.2*44.2)) w l lc rgb "purple" t 'ZeroEq'
DATA="."
SCRIPT="../../../../../bin/postprocessing/extractlinedata.py"
runSim () {
./$1 $INPUT -Problem.Name $2 | tee -a logfile.out
input=`ls -ltr bfs*vtu | tail -n 1 | awk '{print $9}'`
echo $input" -> "$2 | tee -a logfile.out
pvpython $SCRIPT -f $input -o $DATA -of $2"_velocity_m4" -p1 $VProfile1_P1 -p2 $VProfile1_P2 -v 2 -r 10000 | tee -a logfile.out
pvpython $SCRIPT -f $input -o $DATA -of $2"_velocity_1" -p1 $VProfile1_P1 -p2 $VProfile1_P2 -v 2 -r 10000 | tee -a logfile.out
pvpython $SCRIPT -f $input -o $DATA -of $2"_velocity_4" -p1 $VProfile4_P1 -p2 $VProfile4_P2 -v 2 -r 10000 | tee -a logfile.out
pvpython $SCRIPT -f $input -o $DATA -of $2"_velocity_6" -p1 $VProfile6_P2 -p2 $VProfile6_P2 -v 2 -r 10000 | tee -a logfile.out
pvpython $SCRIPT -f $input -o $DATA -of $2"_velocity_10" -p1 $VProfile10_P2 -p2 $VProfile10_P2 -v 2 -r 10000 | tee -a logfile.out
pvpython $SCRIPT -f $input -o $DATA -of $2"_BaseProfile" -p1 $BaseProfile_P1 -p2 $BaseProfile_P2 -v 2 -r 10000 | tee -a logfile.out
}
### backwardsfacingstep
INPUT=scripts/params_verbose.input
VProfilem4_P1="-4.0 1.0 0.0"
VProfilem4_P2="-4.0 9.0 0.0"
VProfile1_P1="1.0 0.0 0.0"
VProfile1_P2="1.0 9.0 0.0"
VProfile4_P1="4.0 0.0 0.0"
VProfile4_P2="4.0 9.0 0.0"
VProfile6_P1="6.0 0.0 0.0"
VProfile6_P2="6.0 9.0 0.0"
VProfile10_P1="10.0 0.0 0.0"
VProfile10_P2="10.0 9.0 0.0"
BaseProfile_P1="0.0 0.0 0.0"
BaseProfile_P2="40.0 0.0 0.0"
runSim test_ff_rans_bfs_kepsilon bfs_kepsilon
runSim test_ff_rans_bfs_komega bfs_komega
runSim test_ff_rans_bfs_lowrekepsilon bfs_lowrekepsilon
runSim test_ff_rans_bfs_oneeq bfs_oneeq
runSim test_ff_rans_bfs_zeroeq bfs_zeroeq
gnuplot scripts/coefficientoffriction.gp
gnuplot scripts/velocitydistribution_m4.gp
gnuplot scripts/velocitydistribution_1.gp
gnuplot scripts/velocitydistribution_4.gp
gnuplot scripts/velocitydistribution_6.gp
gnuplot scripts/velocitydistribution_10.gp
eog *.png
reset
set term pngcairo size 800,600 solid
set output "./VelocityDistribution_x1.png"
set datafile separator ','
DATA='./'
set xlabel "v_x/v_{x,max} [-]"
set ylabel "y [-]"
set yrange [0:5]
set xrange [-0.2:1.01]
set key left center samplen 1
plot './references/bfs_experimentaldata_velocity.csv' u 6:5 w p t 'DriverSeegmiller Experimental Data' ,\
'./references/bfs_referencemodel_velocity.csv' u 6:5 w l lc rgb "black" t 'CFL3D Comparison Model'
reset
set term pngcairo size 800,600 solid
set output "./VelocityDistribution_x10.png"
set datafile separator ','
DATA='./'
set xlabel "v_x/v_{x,max} [-]"
set ylabel "y [-]"
set yrange [0:5]
set xrange [-0.2:1.01]
set key left center samplen 1
plot './references/bfs_experimentaldata_velocity.csv' u 15:14 w p t 'DriverSeegmiller Experimental Data' ,\
'./references/bfs_referencemodel_velocity.csv' u 15:14 w l lc rgb "black" t 'CFL3D Comparison Model'
reset
set term pngcairo size 800,600 solid
set output "./VelocityDistribution_x4.png"
set datafile separator ','
DATA='./'
set xlabel "v_x/v_{x,max} [-]"
set ylabel "y [-]"
set yrange [0:5]
set xrange [-0.2:1.01]
set key left center samplen 1
plot './references/bfs_experimentaldata_velocity.csv' u 9:8 w p t 'DriverSeegmiller Experimental Data' ,\
'./references/bfs_referencemodel_velocity.csv' u 9:8 w l lc rgb "black" t 'CFL3D Comparison Model'
reset
set term pngcairo size 800,600 solid
set output "./VelocityDistribution_x6.png"
set datafile separator ','
DATA='./'
set xlabel "v_x/v_{x,max} [-]"
set ylabel "y [-]"
set yrange [0:5]
set xrange [-0.2:1.01]
set key left center samplen 1
plot './references/bfs_experimentaldata_velocity.csv' u 12:11 w p t 'DriverSeegmiller Experimental Data' ,\
'./references/bfs_referencemodel_velocity.csv' u 12:11 w l lc rgb "black" t 'CFL3D Comparison Model'
reset
set term pngcairo size 800,600 solid
set output "./VelocityDistribution_x-4.png"
set datafile separator ','
DATA='./'
set xlabel "v_x/v_{x,max} [-]"
set ylabel "y [-]"
set yrange [1:5]
set xrange [0:1.025]
set key left center samplen 1
plot './references/bfs_experimentaldata_velocity.csv' u 3:2 w p t 'DriverSeegmiller Experimental Data' ,\
'./references/bfs_referencemodel_velocity.csv' u 3:2 w l lc rgb "black" t 'CFL3D Comparison Model'
add_input_file_links()
dune_symlink_to_source_files(FILES references)
dune_symlink_to_source_files(FILES references scripts)
add_executable(test_ff_rans_lauferpipe_kepsilon EXCLUDE_FROM_ALL main.cc)
target_compile_definitions(test_ff_rans_lauferpipe_kepsilon PUBLIC "KEPSILON=1")
......
DATA="."
SCRIPT="../../../../bin/postprocessing/extractlinedata.py"
SCRIPT="../../../../../bin/postprocessing/extractlinedata.py"
runSim () {
./$1 $INPUT -Problem.Name $2 | tee -a logfile.out
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment