-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathparameters.edp
102 lines (100 loc) · 3.01 KB
/
parameters.edp
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
// Parameters:
//
// real gamma: specific heats ratio.
//
// real alpha: the angle of attack, in radiants.
//
// real cinf: the speed of sound in the medium.
//
// real Ma: the infinite flow Mach number.
//
//
// Inlet conditions:
//
// real uinf: far-field velocity
//
// real rhoinf: far-field density
//
// real pinf: far-field pressure
//
// real rhoEinf: far-field total energy
//
real gamma = 1.4;
real gammam1 = gamma - 1.0;
real Ma = 1.8;
real alphadeg = 0.0;
real alpha = (alphadeg*2*pi)/360.0;
func rhoinf = 1.2;
func pinf = 101325.0;
func cinf = sqrt(gamma*pinf/rhoinf);
func uinf = Ma*cinf;
func uxinf = uinf*cos(alpha);
func uyinf = uinf*sin(alpha);
func rhouxinf = rhoinf*uxinf;
func rhouyinf = rhoinf*uyinf;
func rhoEinf = pinf/gammam1 + 0.5*rhoinf*(uxinf^2+uyinf^2);
real[int] colormapHSVjet = [
0.6667, 1.0000, 0.5625,
0.6667, 1.0000, 0.6250,
0.6667, 1.0000, 0.6875,
0.6667, 1.0000, 0.7500,
0.6667, 1.0000, 0.8125,
0.6667, 1.0000, 0.8750,
0.6667, 1.0000, 0.9375,
0.6667, 1.0000, 1.0000,
0.6563, 1.0000, 1.0000,
0.6458, 1.0000, 1.0000,
0.6354, 1.0000, 1.0000,
0.6250, 1.0000, 1.0000,
0.6146, 1.0000, 1.0000,
0.6042, 1.0000, 1.0000,
0.5938, 1.0000, 1.0000,
0.5833, 1.0000, 1.0000,
0.5729, 1.0000, 1.0000,
0.5625, 1.0000, 1.0000,
0.5521, 1.0000, 1.0000,
0.5417, 1.0000, 1.0000,
0.5313, 1.0000, 1.0000,
0.5208, 1.0000, 1.0000,
0.5104, 1.0000, 1.0000,
0.5000, 1.0000, 1.0000,
0.4889, 0.9375, 1.0000,
0.4762, 0.8750, 1.0000,
0.4615, 0.8125, 1.0000,
0.4444, 0.7500, 1.0000,
0.4242, 0.6875, 1.0000,
0.4000, 0.6250, 1.0000,
0.3704, 0.5625, 1.0000,
0.3333, 0.5000, 1.0000,
0.2963, 0.5625, 1.0000,
0.2667, 0.6250, 1.0000,
0.2424, 0.6875, 1.0000,
0.2222, 0.7500, 1.0000,
0.2051, 0.8125, 1.0000,
0.1905, 0.8750, 1.0000,
0.1778, 0.9375, 1.0000,
0.1667, 1.0000, 1.0000,
0.1563, 1.0000, 1.0000,
0.1458, 1.0000, 1.0000,
0.1354, 1.0000, 1.0000,
0.1250, 1.0000, 1.0000,
0.1146, 1.0000, 1.0000,
0.1042, 1.0000, 1.0000,
0.0938, 1.0000, 1.0000,
0.0833, 1.0000, 1.0000,
0.0729, 1.0000, 1.0000,
0.0625, 1.0000, 1.0000,
0.0521, 1.0000, 1.0000,
0.0417, 1.0000, 1.0000,
0.0313, 1.0000, 1.0000,
0.0208, 1.0000, 1.0000,
0.0104, 1.0000, 1.0000,
0, 1.0000, 1.0000,
0, 1.0000, 0.9375,
0, 1.0000, 0.8750,
0, 1.0000, 0.8125,
0, 1.0000, 0.7500,
0, 1.0000, 0.6875,
0, 1.0000, 0.6250,
0, 1.0000, 0.5625,
0, 1.0000, 0.5000];