FCOFFS/two_phase_transient_test.py at main · UTATRocketry/FCOFFS · 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
44
45
46
47
48
49
50
51
52

import os

from FCOFFS.utilities.units import *
from FCOFFS.components import *
from FCOFFS.interfaces.interface import *
from FCOFFS.systems.transient import *

TS = TransientSolver(ref_p=UnitValue("IMPERIAL", "PRESSURE", "psi", 850))

interface1 = Interface("INTER1")
interface2 = Interface("INTER2")
interface3 = Interface("INTER3")
interface4 = Interface("INTER4")
interface5 = Interface("INTER5")
interface6 = Interface("INTER6")
interface7 = Interface("INTER7")

inlet = pressure_tank_inlet.PressurantTank(TS.quasi_steady_solver, UnitValue.create_unit("inch", 0.5), "N2", UnitValue.create_unit("psi", 1000), UnitValue.create_unit("C", 20), UnitValue.create_unit("m^3", 1), UnitValue.create_unit("m/s", 10))
# reg = pressure_regulator.PressureRegulator(TS.quasi_steady_solver, UnitValue.create_unit("inch", 0.25), "N2", os.path.join(os.getcwd(), "DOCS", "Components", "Pressure Regulator", "Regulator Curves", "KPF.csv"), UnitValue.create_unit("psi", 800))
p = pipe.Pipe(TS.quasi_steady_solver, UnitValue.create_unit("in", 0.5), "N2", UnitValue.create_unit("m", 0.6))
#inlet = mass_flow_inlet.MassFlowInlet(TS.quasi_steady_solver, UnitValue.create_unit("inch", 0.25), "N2", UnitValue.create_unit("kg/s", 0), UnitValue.create_unit("C", 20),  UnitValue.create_unit("m/s", 10))
tp_tank = two_phase_tank.TwoPhaseTank(TS.quasi_steady_solver, UnitValue.create_unit("inch", 0.5), UnitValue.create_unit("inch", 0.25), "N2", "C2H6O", UnitValue.create_unit("kg", 1.2), UnitValue.create_unit("C", 20), UnitValue.create_unit("psi", 950), UnitValue.create_unit("m", 0.05), UnitValue.create_unit("m", 0.24))
p2 = pipe.Pipe(TS.quasi_steady_solver, UnitValue.create_unit("in", 0.25), "C2H6O", UnitValue.create_unit("m", 0.4))
#venturri = cavitating_venturi.CavitatingVenturi(TS.quasi_steady_solver,  UnitValue.create_unit("in", 0.25), UnitValue.create_unit("in", 0.25), UnitValue.create_unit("in", 0.076), "C2H6O", 0.95)
#p3 = pipe.Pipe(TS.quasi_steady_solver, UnitValue.create_unit("in", 0.25), "C2H6O", UnitValue.create_unit("m", 0.4))
out = mass_flow_outlet.MassFlowOutlet(TS.quasi_steady_solver, UnitValue.create_unit("in", 0.25), "C2H6O", UnitValue.create_unit("kg/s", 0.7))
# out = pressure_outlet.PressureOutlet(TS.quasi_steady_solver, UnitValue.create_unit("in", 0.25), "C2H6O", UnitValue.create_unit("psi", 500))

# reg = pressure_regulator.PressureRegulator(TS.quasi_steady_solver, UnitValue.create_unit("in", 0.25), "N2", os.path.join(os.getcwd(), "DOCS", "Components", "Pressure Regulator", "Regulator Curves", "KPF.csv"), UnitValue.create_unit("psi", 800))


inlet.set_connection(downstream=interface1)
#reg.set_connection(interface1, interface2)
p.set_connection(interface1, interface2)
tp_tank.set_connection(interface2, interface3)
p2.set_connection(interface3, interface4)
#venturri.set_connection(interface4, interface5)
#p3.set_connection(interface5, interface6)
out.set_connection(upstream=interface4)

TS.initialize([inlet, p, tp_tank, p2, out])

#1.8kg

#SS.Output.show_tree()
TS.Output.add_probes([(tp_tank, "liquid_height"), (interface2, "mdot"), (tp_tank, "liquid_mass"), (inlet, "p"), (tp_tank, "tank_pressure"), (interface3, "mdot")])
TS.Output.toggle_steady_state_output()
TS.Output.toggle_convergence_output()
TS.Output.set_ouput_unit("psi")
#TS.Output.show_config()
TS.solve(1, 0.01)