FCOFFS/venturi_testing.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
from FCOFFS.utilities.units import *
#from FCOFFS.utilities.units import * # .. \
from FCOFFS.components import *
from FCOFFS.interfaces.interface import *
from FCOFFS.systems.steady import *


PS = SteadySolver(ref_p=UnitValue("IMPERIAL", "PRESSURE", "psi", 15))

#can be dynamically initialized now
interface1 = Interface("INTER1")
interface2 = Interface("INTER2")
interface3 = Interface("INTER3")
interface4 = Interface("INTER4")

fluid = "C2H6O"

inlet = pressure_inlet.PressureInlet(PS, UnitValue.create_unit("in", 0.25), fluid, UnitValue.create_unit("psi", 800), UnitValue.create_unit("c",25), UnitValue.create_unit("m/s", 0.5), "pressure_inlet")
p = pipe.Pipe(PS, UnitValue.create_unit("in", 0.25), fluid, UnitValue.create_unit("ft", 2))
venturi = cavitating_venturi.CavitatingVenturi(PS, UnitValue.create_unit("in", 0.25), UnitValue.create_unit("in", 0.25), UnitValue.create_unit("in", 0.01), fluid, 0.95)
#tube = pipe.Pipe(PS, UnitValue.create_unit("in", 0.25), fluid, UnitValue.create_unit("in", 10))
#injector_manifold = injector.Injector(PS,)
p2 = pipe.Pipe(PS, UnitValue.create_unit("in", 0.25), fluid, UnitValue.create_unit("m", 0.5))
outlet = pressure_outlet.PressureOutlet(PS, UnitValue.create_unit("in",0.25) , fluid, UnitValue.create_unit("psi", 300), "pressure_outlet")

inlet.set_connection(downstream=interface1)
p.set_connection(interface1, interface2)
venturi.set_connection(interface2, interface3)
p2.set_connection(interface3, interface4)
#tube.set_connection(interface2, interface3)
outlet.set_connection(upstream=interface4)

'''
 inter_face
'''

PS.initialize([inlet, p, venturi, p2, outlet])

#PS.Output.show_tree()
PS.Output.toggle_steady_state_output()
PS.Output.set_ouput_unit("psi")
#PS.Output.toggle_convergence_output()
PS.solve()