Open sashahafner opened 1 year ago
@afeilberg see demo 12. I had to reduce the reaction and mass transfer rates to see a clear effect. The liquid profile plot looks right, I think:
The blue line has recirc = True
, red does not (others have no reaction, with and without recirculation).
The blue line starts (at inlet) a bit above the outlet concentration.
This did not require much new code.
In the rates()
function the inlet liquid concentration clin
is set to the concentration in the last layer around line 93:
# Get outlet liquid phase concentration to use at inlet if recirc = True
if recirc:
if counter:
oi = 0
else:
oi = nc - 1
clin = ccl[oi]
It is meant to work for countercurrent flow also, but I have not tested it.
Comment from @afeilberg : coudl be nice to have a delay to represent retention time in recirculation system. Or buffer tank.
Buffer tank could actually be easy(ish). Could include as new state variables.
Anders says we could add reactions in buffer tank.
Try to work on this over the next week.