def test_two_recycle_loops_with_complete_overlap():
f.set_flowsheet('two_recycle_loops_with_complete_overlap')
settings.set_thermo(['Water'], cache=True)
feedstock = Stream('feedstock', Water=1000)
water = Stream('water', Water=10)
recycle = Stream('recycle')
inner_recycle = Stream('inner_recycle')
product = Stream('product')
inner_water = Stream('inner_water', Water=10)
P1 = Pump('P1', feedstock)
P2 = Pump('P2', water)
P3 = Pump('P3', inner_water)
M1 = Mixer('M1', [P1 - 0, P2 - 0, recycle])
M2 = Mixer('M2', [M1 - 0, P3 - 0, inner_recycle])
S2 = Splitter('S2', M2 - 0, ['', inner_recycle], split=0.5)
S1 = Splitter('S1', S2 - 0, [product, recycle], split=0.5)
recycle_loop_sys = f.create_system('recycle_loop_sys')
network = recycle_loop_sys.to_network()
actual_network = Network([
P1, P2, P3,
Network([M1, Network([M2, S2], recycle=inner_recycle), S1],
recycle=recycle)
])
assert network == actual_network
recycle_loop_sys.simulate()
x_nested_solution = np.vstack([recycle.mol, inner_recycle.mol])
recycle_loop_sys.flatten()
assert recycle_loop_sys.path == (P1, P2, P3, M1, M2, S2, S1)
recycle_loop_sys.empty_recycles()
recycle_loop_sys.simulate()
x_flat_solution = np.vstack([recycle.mol, inner_recycle.mol])
assert_allclose(x_nested_solution, x_flat_solution, rtol=1e-2)
def test_feed_forward_recycle_loop():
f.set_flowsheet('feed_forward_recycle_loop')
settings.set_thermo(['Water'], cache=True)
feedstock = Stream('feedstock', Water=1000)
water = Stream('water', Water=10)
recycle = Stream('recycle')
inner_recycle = Stream('inner_recycle')
product = Stream('product')
P1 = Pump('P1', feedstock)
P2 = Pump('P2', water)
M1 = Mixer('M1', [P1 - 0, recycle, inner_recycle])
S1 = Splitter('S1', M1 - 0, ['', recycle], split=0.5)
M2 = Mixer('M2', [P2 - 0, S1 - 0])
S2 = Splitter('S2', M2 - 0, [inner_recycle, product], split=0.5)
P3 = Pump('P3', product)
recycle_loop_sys = f.create_system('recycle_loop_sys')
network = recycle_loop_sys.to_network()
actual_network = Network(
[P1, P2, Network([M1, S1, M2, S2], recycle=S2 - 0), P3])
assert network == actual_network
recycle_loop_sys.simulate()
x_nested_solution = np.vstack([recycle.mol, inner_recycle.mol])
recycle_loop_sys.flatten()
assert recycle_loop_sys.path == (P1, P2, M1, S1, M2, S2, P3)
recycle_loop_sys.empty_recycles()
recycle_loop_sys.simulate()
x_flat_solution = np.vstack([recycle.mol, inner_recycle.mol])
assert_allclose(x_nested_solution, x_flat_solution, rtol=1e-2)
def test_sugarcane_ethanol_biorefinery_network():
from biorefineries.sugarcane import flowsheet as f
sugarcane_sys = f.create_system('sugarcane_sys')
globals().update(f.unit.__dict__)
network = sugarcane_sys.to_network()
actual_network = Network([
U101, U102, U103, T201,
Network([U201, S201, M201],
recycle=M201 - 0), T202, H201, T203, P201, T204, T205, P202,
Network([M202, H202, T206, C201, C202, P203], recycle=P203 - 0), S202,
S301, F301, P306, M301, H301, T305, R301, T301, C301, M302, P301,
Network([H302, D302, P302], recycle=P302 - 0),
Network([M303, D303, H303, U301], recycle=U301 - 0), H304, T302, P304,
T303, P305, M304, T304, D301, P303, M305, U202
])
assert network == actual_network
sugarcane_sys.empty_recycles()
sugarcane_sys.simulate()
bst.process_tools.default()
def test_recycle_loop():
f.set_flowsheet('simple_recycle_loop')
settings.set_thermo(['Water'], cache=True)
feedstock = Stream('feedstock', Water=1000)
water = Stream('water', Water=10)
recycle = Stream('recycle')
product = Stream('product')
P1 = Pump('P1', feedstock)
P2 = Pump('P2', water)
M1 = Mixer('M1', [P1 - 0, P2 - 0, recycle])
S1 = Splitter('S1', M1 - 0, [product, recycle], split=0.5)
recycle_loop_sys = f.create_system('recycle_loop_sys')
network = recycle_loop_sys.to_network()
actual_network = Network([P1, P2, Network([M1, S1], recycle=recycle)])
assert network == actual_network
recycle_loop_sys.simulate()
x_nested_solution = recycle.mol.copy()
recycle_loop_sys.flatten()
assert recycle_loop_sys.path == (P1, P2, M1, S1)
recycle_loop_sys.empty_recycles()
recycle_loop_sys.simulate()
x_flat_solution = recycle.mol.copy()
assert_allclose(x_nested_solution, x_flat_solution, rtol=1e-2)
def test_separate_recycle_loops():
f.set_flowsheet('separate_recycle_loops')
settings.set_thermo(['Water'], cache=True)
feedstock_a = Stream('feedstock_a', Water=1000)
water_a = Stream('water_a', Water=10)
recycle_a = Stream('recycle_a')
product_a = Stream('product_a')
P1_a = Pump('P1_a', feedstock_a)
P2_a = Pump('P2_a', water_a)
M1_a = Mixer('M1_a', [P1_a - 0, P2_a - 0, recycle_a])
S1_a = Splitter('S1_a', M1_a - 0, [product_a, recycle_a], split=0.5)
feedstock_b = Stream('feedstock_b', Water=1000)
water_b = Stream('water_b', Water=10)
recycle_b = Stream('recycle_b')
product_b = Stream('product_b')
P1_b = Pump('P1_b', feedstock_b)
P2_b = Pump('P2_b', water_b)
M1_b = Mixer('M1_b', [P1_b - 0, P2_b - 0, recycle_b])
S1_b = Splitter('S1_b', M1_b - 0, [product_b, recycle_b], split=0.5)
recycles = [recycle_a, recycle_b]
recycle_loop_sys = f.create_system('recycle_loop_sys')
recycle_loop_sys.simulate()
network = recycle_loop_sys.to_network()
actual_network = Network([
P1_a, P2_a,
Network([M1_a, S1_a], recycle=S1_a - 1), P1_b, P2_b,
Network([M1_b, S1_b], recycle=S1_b - 1)
])
assert network == actual_network
x_nested_solution = np.vstack([i.mol for i in recycles])
recycle_loop_sys.flatten()
recycle_loop_sys.empty_recycles()
recycle_loop_sys.simulate()
x_flat_solution = np.vstack([i.mol for i in recycles])
assert_allclose(x_nested_solution, x_flat_solution, rtol=1e-2)
assert recycle_loop_sys.path == (P1_a, P2_a, M1_a, S1_a, P1_b, P2_b, M1_b,
S1_b)
def test_nested_recycle_loops():
f.set_flowsheet('feed_forward_recycle_loop')
settings.set_thermo(['Water'], cache=True)
feedstock = Stream('feedstock', Water=1000)
recycle_1, recycle_2, recycle_3, recycle_4, recycle_5 = recycles = [
Stream('recycle_1'),
Stream('recycle_2'),
Stream('recycle_3'),
Stream('recycle_4'),
Stream('recycle_5'),
]
feed_1 = Stream('feed_1', Water=10)
feed_2 = Stream('feed_2', Water=10)
feed_3 = Stream('feed_3', Water=10)
feed_4 = Stream('feed_4', Water=10)
inner_recycle = Stream('inner_recycle')
product = Stream('product')
P1 = Pump('P1', feedstock)
M1 = Mixer('M1', [P1 - 0, recycle_1])
P2 = Pump('P2', M1 - 0)
H1 = HXprocess('H1', [P2 - 0, recycle_2])
P3 = Pump('P3', feed_1)
M3 = Mixer('M3', [H1 - 0, P3 - 0], recycle_2)
P4 = Pump('P4', H1 - 1)
M4 = Mixer('M4', [P4 - 0, recycle_3])
P5 = Pump('P5', feed_2)
P6 = Pump('P6', feed_3)
M5 = Mixer('M5', [M4 - 0, P5 - 0, P6 - 0])
H2 = HXprocess('H2', [M5 - 0, recycle_4])
P7 = Pump('P7', feed_4)
M6 = Mixer('M6', [H2 - 0, P7 - 0])
S1 = Splitter('S1', M6 - 0, [recycle_4, ''], split=0.5)
S2 = Splitter('S2', H2 - 1, split=0.5)
P8 = Pump('P8', S2 - 1)
S3 = Splitter('S3', P8 - 0, split=0.5)
H3 = HXprocess('H3', [S2 - 0, recycle_5], ['', recycle_3])
M7 = Mixer('M7', [H3 - 0, S3 - 0])
S4 = Splitter('S4', M7 - 0, ['', recycle_5], split=0.5)
S5 = Splitter('S5', S4 - 0, ['', recycle_1], split=0.5)
M8 = Mixer('M8', [S3 - 1, S1 - 1, S5 - 0], product)
recycle_loop_sys = f.create_system('recycle_loop_sys')
recycle_loop_sys.simulate()
network = recycle_loop_sys.to_network()
actual_network = Network([
P1, P3, P5, P6, P7,
Network([
M1, P2,
Network([H1, M3], recycle=M3 - 0), P4,
Network([
M4, M5,
Network([H2, M6, S1], recycle=S1 - 0), S2,
Network([M7, S4, H3], recycle=H3 - 0), P8, S3
],
recycle=H3 - 1), S5
],
recycle=S5 - 1), M8
])
assert network == actual_network
x_nested_solution = np.vstack([i.mol for i in recycles])
recycle_loop_sys.flatten()
recycle_loop_sys.empty_recycles()
recycle_loop_sys.simulate()
x_flat_solution = np.vstack([i.mol for i in recycles])
assert_allclose(x_nested_solution, x_flat_solution, rtol=1e-2)