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

interface1 = Interface("INTER1")
interface2 = Interface("INTER2")
interface3 = Interface("INTER3")
interface4 = Interface("INTER4")

SS = SteadySolver(ref_p=UnitValue.create_unit("psi", 700), ref_T=UnitValue.create_unit("C", 20))

pipe_diameter = 0.18
orifice_diameter = 0.001

inlet = pressure_inlet.PressureInlet(SS, UnitValue.create_unit("in", pipe_diameter), "N2", UnitValue.create_unit("psi", 800), UnitValue.create_unit("C", 20), UnitValue.create_unit("m/s", 5))
p = pipe.Pipe(SS, UnitValue.create_unit("in", pipe_diameter), "N2", UnitValue.create_unit("in", 10), roughness=0.015)
orifice = critical_orifice.CriticalOrifice(SS, UnitValue.create_unit("in", pipe_diameter), UnitValue.create_unit("in", pipe_diameter), UnitValue.create_unit("in", orifice_diameter), "N2")
#p2 = pipe.Pipe(SS, UnitValue.create_unit("in", pipe_diameter), "N2", UnitValue.create_unit("in", 10), roughness=0.015)
#outlet = mass_flow_outlet.MassFlowOutlet(SS, UnitValue.create_unit("in", pipe_diameter), "N2", )
outlet = pressure_outlet.PressureOutlet(SS, UnitValue.create_unit("in", pipe_diameter), "N2", UnitValue.create_unit("psig", 0))
#outlet = flowrate_outlet.FlowRateOutlet(SS, UnitValue.create_unit("in", 0.18), "N2", UnitValue.create_unit("SCFM", 106.5))

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

SS.initialize([inlet, p, orifice, outlet])
SS.Output.set_ouput_unit("psi")
SS.Output.toggle_steady_state_output()
SS.solve()

Re = outlet.interface_in.state.rho*outlet.interface_in.state.u*outlet.diameter/Fluid.dynamic_viscosity("N2", outlet.interface_in.state.rho, outlet.interface_in.state.T, outlet.interface_in.state.p)
flow = outlet.interface_in.state.mdot/outlet.interface_in.state.rho
flow2 = orifice.interface_in.state.mdot/orifice.interface_in.state.rho
print("Reynolds Number:", Re)
print(flow2, flow)
print("Outlet Flow Rate:", flow2.to("SCFM", temperature=orifice.interface_in.state.T, pressure=orifice.interface_in.state.p), flow.to("SCFM", temperature=outlet.interface_in.state.T, pressure=outlet.interface_in.state.p))