1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
#!/bin/sh
#------------------------------------------------------------------------------
# ========= |
# \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
# \\ / O peration |
# \\ / A nd | Copyright (C) 2016 OpenCFD Ltd.
# \\/ M anipulation |
#-------------------------------------------------------------------------------
# License
# This file is part of OpenFOAM.
#
# OpenFOAM is free software: you can redistribute it and/or modify it
# under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
# FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
# for more details.
#
# You should have received a copy of the GNU General Public License
# along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
#
# Script
# createGraphs
#
# Description
# Creates .eps graphs of OpenFOAM results vs experiment for the buoyant
# cavity case
#
#------------------------------------------------------------------------------
createEpsT()
{
OF=$1
OFPROBE=$2
EXPT=$3
EXPPROBE=$4
OUTPUT=$5
gnuplot<<EOF
set terminal postscript default
set output "($OUTPUT) OF_vs_EXPT($EXPPROBE).eps"
set xlabel "time [sec]"
set ylabel "$OUTPUT"
set grid
plot \
'$EXPT' using 1:$4 title "Exp (probe $EXPPROBE)" , \
'$OF' using 1:(\$$OFPROBE - 273) title "OpenFOAM"
EOF
}
# Test if gnuplot exists on the system
command -v gnuplot >/dev/null 2>&1 || {
echo "gnuplot not found - skipping graph creation" 1>&2
exit 1
}
SETSDIRT="../postProcessing/thermoCouple"
[ -d $SETSDIRT ] || echo "No probe data available to create graphs" && exit 1
arr=(`ls $SETSDIRT`)
OFDATAROOTT="$SETSDIRT/${arr[0]}"
# paths to data
EXPTDATAROOT=./Probes/Temperature
DATAFILET="$EXPTDATAROOT/temperature.dat"
# generate temperature profiles
OF="$OFDATAROOTT/T"
EXPT="$DATAFILET"
echo $OF
echo $EXPT
#
#createEpsT $OF 2 $EXPT 12 Temperature
#createEpsT $OF 3 $EXPT 10 Temperature
#createEpsT $OF 5 $EXPT 3 Temperature
createEpsT $OF 6 $EXPT 18 Temperature
createEpsT $OF 7 $EXPT 19 Temperature
createEpsT $OF 8 $EXPT 20 Temperature
#createEpsT $OF 9 $EXPT 13 Temperature
#createEpsT $OF 10 $EXPT 2 Temperature
#createEpsT $OF 11 $EXPT 17 Temperature
#createEpsT $OF 12 $EXPT 16 Temperature
#createEpsT $OF 13 $EXPT 15 Temperature
#createEpsT $OF 14 $EXPT 14 Temperature
SETSDIRT="../postProcessing/probes_O2"
arr=(`ls $SETSDIRT`)
OFDATAROOTT="$SETSDIRT/${arr[0]}/O2"
# paths to data
EXPTDATAROOT=./Probes/Oxygen/oxygen.dat
#createEpsT $OFDATAROOTT 3 $EXPTDATAROOT 2 O2
#createEpsT $OFDATAROOTT 2 $EXPTDATAROOT 3 O2
#SETSDIRT="../postProcessing/patchInlet_phi"
#arr=(`ls $SETSDIRT`)
#OFDATAROOTT="$SETSDIRT/${arr[0]}/faceSource.dat"
# paths to data
#EXPTDATAROOT=./Probes/phi/massInlet.dat
#createEpsT $OFDATAROOTT 2 $EXPTDATAROOT 3 MassInlet
echo Done
#------------------------------------------------------------------------------
|
