def time_differential_pressure_meter_solver_D2_Miller_orifice(self):
differential_pressure_meter_solver(D=0.07366,
P1=200000.0,
P2=183000.0,
m=7.702338035732167,
rho=999.1,
mu=0.0011,
k=1.33,
meter_type='Miller orifice',
taps='corner')
def time_differential_pressure_meter_solver_m_Miller_orifice(self):
differential_pressure_meter_solver(D=0.07366,
D2=0.05,
P1=200000.0,
P2=183000.0,
rho=999.1,
mu=0.0011,
k=1.33,
meter_type='Miller orifice',
taps='corner')
def time_differential_pressure_meter_solver_D2_Hollingshead(self):
differential_pressure_meter_solver(D=0.07366,
P1=200000.0,
P2=183000.0,
m=7.702338035732167,
rho=999.1,
mu=0.0011,
k=1.33,
meter_type='Hollingshead orifice',
taps='D')
def time_differential_pressure_meter_solver_m_Hollingshead(self):
differential_pressure_meter_solver(D=0.07366,
D2=0.05,
P1=200000.0,
P2=183000.0,
rho=999.1,
mu=0.0011,
k=1.33,
meter_type='Hollingshead orifice',
taps='D')
def time_differential_pressure_meter_solver_P1(self):
differential_pressure_meter_solver(D=0.07366,
D2=0.05,
P2=183000.0,
m=7.702338035732167,
rho=999.1,
mu=0.0011,
k=1.33,
meter_type='ISO 5167 orifice',
taps='D')
def run():
D = 0.25 * u.inch
# rho = 2.11 * u.lb / u.ft ** 3 # density of helium at 3000 psi
rho = 7.3935 * u.lb / u.ft**3
Q = 0.336 * u.ft**3 / u.s
Do = (fl.differential_pressure_meter_solver(
D=D.to('m').magnitude,
rho=rho.to('kg / m ** 3').magnitude,
mu=N2_DynamicViscosity_DEFAULT.magnitude,
k=1.4,
P1=(3000 * u.psi).to('Pa').magnitude,
P2=(2900 * u.psi).to('Pa').magnitude,
m=(Q * rho).to('kg / s').magnitude,
meter_type='ISO 5167 orifice',
taps='flange') * u.m).to('in')
# DP = 80 * u.psi
# Y = 0.99
# guessarr = (0.5, 1.0, 0.2)
# data = (D, DP, rho, Q, Y)
# sol = sp.fsolve(equations, guessarr, args=data)
# d = (1 / ((1 / D ** 4) + DP / rho * (0.065607 * C * Y / Q) ** 2)) ** 0.25
# print(str(sol))
print(str(Do))
print('Done')