results.jl

using PyPlot;


souTriP = [
7143	0
7143	0.01
7143	0.02
7143	0.03
7143	0.04
7143	0.05
7143	0.06
7143	0.07
7143	0.08
7143	0.09
7143	0.1
7143	0.11
7143	0.12
7143	0.13
7143	0.14
7143	0.15
7143	0.16
7143	0.17
7143	0.18
7143	0.19
7143	0.2
7143	0.21
7143	0.22
7143	0.23
7143	0.24
7143	0.25
7143	0.26
7143	0.27
7143	0.28
7143	0.29
7143	0.3
7143	0.31
7143	0.32
7143	0.33
7143	0.34
7143	0.35
7143	0.36
7143	0.37
7143	0.38
7143	0.39
7143	0.4
7143	0.41
7143	0.42
7143	0.43
7143	0.44
7143	0.45
7143.1	0.46
7143.1	0.47
7143.1	0.48
7143.1	0.49
7143.2	0.5
7143.3	0.51
7143.5	0.52
7143.7	0.53
7143.8	0.54
7144.2	0.55
7145	0.56
7145.9	0.57
7146.7	0.58
7147.5	0.59
7148.7	0.6
7152.4	0.61
7155.7	0.62
7156.8	0.63
7158	0.64
7160.9	0.65
7164.1	0.66
7167.7	0.67
7175.4	0.68
7197.3	0.69
7220.7	0.7
7250.4	0.71
7280.2	0.72
7310	0.73
7354.5	0.74
7490.8	0.75
7582.7	0.76
7674.5	0.77
7766.4	0.78
7858.3	0.79
8040.1	0.8
8292	0.81
8544	0.82
8796	0.83
9048	0.84
9324.6	0.85
9700.1	0.86
10076	0.87
10451	0.88
10827	0.89
11186	0.9
11552	0.91
11930	0.92
12308	0.93
12686	0.94
13047	0.95
13282	0.96
13540	0.97
13799	0.98
14057	0.99
14311	1
14423	1.01
14553	1.02
14684	1.03
14814	1.04
14944	1.05
14990	1.06
15044	1.07
15097	1.08
15151	1.09
15204	1.1
15218	1.11
15233	1.12
15248	1.13
15263	1.14
15278	1.15
15286	1.16
15295	1.17
15304	1.18
15314	1.19
15323	1.2
15327	1.21
15330	1.22
15334	1.23
15338	1.24
15342	1.25
15341	1.26
15340	1.27
15339	1.28
15338	1.29
15337	1.3
15336	1.31
15337	1.32
15337	1.33
15338	1.34
15338	1.35
15340	1.36
15342	1.37
15344	1.38
15346	1.39
15348	1.4
15348	1.41
15349	1.42
15349	1.43
15349	1.44
15349	1.45
15349	1.46
15348	1.47
15347	1.48
15347	1.49
15346	1.5
15346	1.51
15347	1.52
15347	1.53
15348	1.54
15348	1.55
15349	1.56
15350	1.57
15350	1.58
15351	1.59
15352	1.6
15351	1.61
15351	1.62
15351	1.63
15350	1.64
15350	1.65
15349	1.66
15349	1.67
15349	1.68
15349	1.69
15348	1.7
15349	1.71
15349	1.72
15349	1.73
15349	1.74
15350	1.75
15349	1.76
15349	1.77
15349	1.78
15349	1.79
15349	1.8
15349	1.81
15348	1.82
15348	1.83
15348	1.84
15347	1.85
15347	1.86
15347	1.87
15347	1.88
15346	1.89
15346	1.9
15346	1.91
15346	1.92
15346	1.93
15346	1.94
15346	1.95
15345	1.96
15345	1.97
15345	1.98
15345	1.99
15344	2
15344	2.01
15344	2.02
15343	2.03
15343	2.04
15343	2.05
15343	2.06
15343	2.07
15342	2.08
15342	2.09
15342	2.1
15342	2.11
15341	2.12
15341	2.13
15341	2.14
15341	2.15
15341	2.16
15341	2.17
15340	2.18
15340	2.19
15340	2.2
15340	2.21
15340	2.22
15340	2.23
15340	2.24
15340	2.25
15339	2.26
15339	2.27
15339	2.28
15339	2.29
15339	2.3
15339	2.31
15339	2.32
15339	2.33
15339	2.34
15339	2.35
15339	2.36
15339	2.37
15338	2.38
15338	2.39
15338	2.4
15338	2.41
15339	2.42
15339	2.43
15340	2.44
15340	2.45
15340	2.46
15340	2.47
15341	2.48
15341	2.49
15342	2.5
15342	2.51
15343	2.52
15345	2.53
15347	2.54
15348	2.55
15348	2.56
15348	2.57
15349	2.58
15349	2.59
15350	2.6
15352	2.61
15356	2.62
15359	2.63
15363	2.64
15366	2.65
15378	2.66
15402	2.67
15426	2.68
15449	2.69
15469	2.7
15485	2.71
15514	2.72
15542	2.73
15570	2.74
15599	2.75
15616	2.76
15705	2.77
15869	2.78
16033	2.79
16207	2.8
16433	2.81
16660	2.82
16886	2.83
17126	2.84
17346	2.85
17601	2.86
17856	2.87
18111	2.88
18366	2.89
18764	2.9
19368	2.91
19972	2.92
20579	2.93
21193	2.94
21794	2.95
22395	2.96
22996	2.97
23597	2.98
24199	2.99
24725	3
25198	3.01
25672	3.02
26146	3.03
26619	3.04
27002	3.05
27318	3.06
27635	3.07
27951	3.08
28268	3.09
28504	3.1
28658	3.11
28812	3.12
28967	3.13
29121	3.14
29232	3.15
29294	3.16
29357	3.17
29419	3.18
29482	3.19
29526	3.2
29544	3.21
29562	3.22
29581	3.23
29599	3.24
29612	3.25
29615	3.26
29619	3.27
29622	3.28
29626	3.29
29628	3.3
29627	3.31
29627	3.32
29626	3.33
29626	3.34
29625	3.35
29624	3.36
29623	3.37
29622	3.38
29620	3.39
29619	3.4
29618	3.41
29617	3.42
29616	3.43
29615	3.44
29614	3.45
29614	3.46
29613	3.47
29613	3.48
29612	3.49
29611	3.5
29611	3.51
29610	3.52
29610	3.53
29609	3.54
29608	3.55
29607	3.56
29606	3.57
29605	3.58
29604	3.59
29603	3.6
29601	3.61
29598	3.62
29595	3.63
29592	3.64
29589	3.65
29586	3.66
29584	3.67
29582	3.68
29580	3.69
29578	3.7
29576	3.71
29575	3.72
29574	3.73
29574	3.74
29573	3.75
29573	3.76
29573	3.77
29572	3.78
29572	3.79
29572	3.8
29572	3.81
29573	3.82
29573	3.83
29573	3.84
29573	3.85
29572	3.86
29572	3.87
29572	3.88
29571	3.89
29570	3.9
29570	3.91
29569	3.92
29569	3.93
29568	3.94
29567	3.95
29566	3.96
29565	3.97
29564	3.98
29564	3.99
29563	4
29563	4.01
29563	4.02
29563	4.03
29563	4.04
29564	4.05
29565	4.06
29565	4.07
29565	4.08
29565	4.09
29564	4.1


];

fouTriP = [
7143	0
7143	0.01
7143	0.02
7143	0.03
7143	0.04
7143	0.05
7143	0.06
7143	0.07
7143	0.08
7143	0.09
7143	0.1
7143	0.11
7143	0.12
7143	0.13
7143	0.14
7143	0.15
7143	0.16
7143	0.17
7143	0.18
7143	0.19
7143	0.2
7143	0.21
7143	0.22
7143	0.23
7143	0.24
7143	0.25
7143	0.26
7143	0.27
7143	0.28
7143	0.29
7143	0.3
7143	0.31
7143	0.32
7143	0.33
7143.1	0.34
7143.1	0.35
7143.2	0.36
7143.2	0.37
7143.2	0.38
7143.3	0.39
7143.9	0.4
7144.3	0.41
7144.6	0.42
7144.9	0.43
7145.2	0.44
7145.8	0.45
7147.2	0.46
7148.6	0.47
7149.9	0.48
7151.3	0.49
7152.6	0.5
7157.6	0.51
7162.6	0.52
7167.6	0.53
7172.5	0.54
7180.2	0.55
7192.1	0.56
7204.3	0.57
7216.5	0.58
7228.6	0.59
7243.2	0.6
7279	0.61
7311.8	0.62
7329.2	0.63
7346.7	0.64
7379.2	0.65
7413	0.66
7455.2	0.67
7512.5	0.68
7607	0.69
7700	0.7
7786.9	0.71
7873.8	0.72
7960.7	0.73
8061.8	0.74
8221.6	0.75
8353.6	0.76
8485.5	0.77
8617.5	0.78
8749.4	0.79
8925	0.8
9149.9	0.81
9374.8	0.82
9599.7	0.83
9824.6	0.84
10053	0.85
10306	0.86
10559	0.87
10812	0.88
11065	0.89
11311	0.9
11548	0.91
11790	0.92
12032	0.93
12273	0.94
12508	0.95
12698	0.96
12897	0.97
13097	0.98
13296	0.99
13493	1
13627	1.01
13772	1.02
13918	1.03
14063	1.04
14209	1.05
14295	1.06
14392	1.07
14488	1.08
14584	1.09
14681	1.1
14734	1.11
14794	1.12
14853	1.13
14912	1.14
14972	1.15
15003	1.16
15038	1.17
15073	1.18
15107	1.19
15142	1.2
15159	1.21
15179	1.22
15198	1.23
15218	1.24
15237	1.25
15247	1.26
15257	1.27
15268	1.28
15278	1.29
15289	1.3
15294	1.31
15299	1.32
15305	1.33
15310	1.34
15316	1.35
15319	1.36
15321	1.37
15324	1.38
15327	1.39
15330	1.4
15331	1.41
15332	1.42
15334	1.43
15335	1.44
15337	1.45
15337	1.46
15338	1.47
15339	1.48
15339	1.49
15340	1.5
15340	1.51
15340	1.52
15341	1.53
15341	1.54
15341	1.55
15341	1.56
15341	1.57
15341	1.58
15342	1.59
15342	1.6
15342	1.61
15342	1.62
15342	1.63
15342	1.64
15342	1.65
15342	1.66
15342	1.67
15342	1.68
15342	1.69
15342	1.7
15341	1.71
15341	1.72
15341	1.73
15341	1.74
15341	1.75
15341	1.76
15341	1.77
15341	1.78
15341	1.79
15341	1.8
15341	1.81
15341	1.82
15341	1.83
15341	1.84
15341	1.85
15341	1.86
15341	1.87
15341	1.88
15340	1.89
15340	1.9
15340	1.91
15340	1.92
15340	1.93
15340	1.94
15340	1.95
15340	1.96
15340	1.97
15340	1.98
15340	1.99
15340	2
15340	2.01
15340	2.02
15340	2.03
15340	2.04
15340	2.05
15341	2.06
15341	2.07
15341	2.08
15341	2.09
15341	2.1
15341	2.11
15342	2.12
15342	2.13
15342	2.14
15343	2.15
15343	2.16
15344	2.17
15344	2.18
15345	2.19
15346	2.2
15347	2.21
15348	2.22
15349	2.23
15350	2.24
15351	2.25
15352	2.26
15355	2.27
15357	2.28
15360	2.29
15362	2.3
15365	2.31
15370	2.32
15375	2.33
15381	2.34
15386	2.35
15392	2.36
15394	2.37
15396	2.38
15397	2.39
15399	2.4
15400	2.41
15415	2.42
15438	2.43
15459	2.44
15475	2.45
15485	2.46
15490	2.47
15511	2.48
15532	2.49
15553	2.5
15573	2.51
15604	2.52
15644	2.53
15684	2.54
15709	2.55
15731	2.56
15748	2.57
15765	2.58
15781	2.59
15811	2.6
15855	2.61
15921	2.62
15987	2.63
16053	2.64
16104	2.65
16196	2.66
16354	2.67
16512	2.68
16670	2.69
16805	2.7
16898	2.71
17031	2.72
17164	2.73
17298	2.74
17431	2.75
17532	2.76
17728	2.77
17972	2.78
18217	2.79
18464	2.8
18724	2.81
18985	2.82
19246	2.83
19516	2.84
19727	2.85
19943	2.86
20158	2.87
20374	2.88
20589	2.89
20874	2.9
21253	2.91
21631	2.92
22010	2.93
22392	2.94
22757	2.95
23103	2.96
23449	2.97
23794	2.98
24140	2.99
24462	3
24772	3.01
25082	3.02
25392	3.03
25702	3.04
25978	3.05
26240	3.06
26502	3.07
26764	3.08
27026	3.09
27247	3.1
27450	3.11
27652	3.12
27855	3.13
28057	3.14
28221	3.15
28363	3.16
28504	3.17
28646	3.18
28788	3.19
28899	3.2
28986	3.21
29072	3.22
29159	3.23
29246	3.24
29312	3.25
29357	3.26
29402	3.27
29447	3.28
29492	3.29
29526	3.3
29544	3.31
29562	3.32
29579	3.33
29597	3.34
29611	3.35
29614	3.36
29617	3.37
29620	3.38
29624	3.39
29626	3.4
29623	3.41
29621	3.42
29618	3.43
29615	3.44
29613	3.45
29609	3.46
29605	3.47
29601	3.48
29596	3.49
29592	3.5
29589	3.51
29586	3.52
29582	3.53
29579	3.54
29575	3.55
29572	3.56
29570	3.57
29567	3.58
29564	3.59
29561	3.6
29559	3.61
29557	3.62
29555	3.63
29553	3.64
29551	3.65
29549	3.66
29548	3.67
29546	3.68
29545	3.69
29543	3.7
29542	3.71
29538	3.72
29535	3.73
29531	3.74
29527	3.75
29524	3.76
29522	3.77
29520	3.78
29519	3.79
29516	3.8
29512	3.81
29508	3.82
29505	3.83
29501	3.84
29498	3.85
29495	3.86
29491	3.87
29488	3.88
29485	3.89
29482	3.9
29480	3.91
29477	3.92
29473	3.93
29469	3.94
29464	3.95
29460	3.96
29455	3.97
29452	3.98
29448	3.99
29444	4
29440	4.01
29436	4.02
29429	4.03
29423	4.04
29417	4.05
29412	4.06
29411	4.07
29410	4.08
29410	4.09
29409	4.1

];


theta = 61.0;
thetaRad = 61.0*pi/180.0;
halfSL = 0.5*tan(thetaRad);

ksiRef = 23.279
ksiRefRad = ksiRef*pi/180.0;
halfSLRef = tan(thetaRad) + 0.5/tan(ksiRefRad);

Minf = 2.9;
Pinf = 7143.0;
gamma = 1.4;

P1_exact  = 7143;
P2_exact  = 15282.0;
P3_exact =  29340.0;

Cp1_exact = 2.0/gamma/Minf/Minf*(P1_exact/Pinf - 1.0);
Cp2_exact = 2.0/gamma/Minf/Minf*(P2_exact/Pinf - 1.0);
Cp3_exact = 2.0/gamma/Minf/Minf*(P3_exact/Pinf - 1.0);


Cp_exact =[

0.0    Cp1_exact
halfSL Cp1_exact 
halfSL Cp2_exact 
halfSLRef  Cp2_exact 
halfSLRef  Cp3_exact 
4.1 	Cp3_exact 
];


Cp_fouTri = 2.0/gamma/Minf/Minf*(fouTriP[:,1]./Pinf .- 1.0);
Cp_souTri = 2.0/gamma/Minf/Minf*(souTriP[:,1]./Pinf .- 1.0);


figure(2)
clf()

plot(Cp_exact[:,1],Cp_exact[:,2], "-k", linewidth = 1.5, label = "exact" );

plot(fouTriP[:,2],Cp_fouTri, "--r", linewidth = 2.0, label = "fouTri" );
plot(souTriP[:,2],Cp_souTri, "-.g", linewidth = 2.0, label = "souTri" );

legend();
xlabel("x")
ylabel("Cp")
xlim(0.0,4.1)