mp-operational-envelope/main.py at master · wdoppenberg/mp-operational-envelope · GitHub

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
from analysis import constants as c, OperationalEnvelope, Experiment
import matplotlib.pyplot as plt

input_vars1 = dict(
    l_tube=0.4,
    d_tube=1.57e-3,
    propellant=c.Water,
    p_0=1.1e5,
    V_fraction=0.12,
    C_d=0.66,
    xi_s=0.3
)

input_vars2 = dict(
    l_tube=0.3,
    d_tube=1.57e-3,
    propellant=c.Water,
    p_0=2e5,
    V_fraction=0.12,
    C_d=0.66,
    xi_s=0.3
)

input_vars3 = dict(
    l_tube=0.3,
    d_tube=1.57e-3,
    propellant=c.Water,
    p_0=1.5e5,
    V_fraction=0.3,
    C_d=0.66,
    xi_s=0.6
)

oe1 = OperationalEnvelope(input_vars1, name='oe1')
oe2 = OperationalEnvelope(input_vars2, name='oe2')
oe3 = OperationalEnvelope(input_vars3, name='oe3')

dt, t_end = 0.1, 800
ex = Experiment(oe1, oe2, oe3)

ex.verification(dt, t_end)
ex.comparison(dt, t_end)
plt.show()