1
Chapter 5 Bistability and oscillations in flow reactors Closed system: batch From non-equilibrium to equilibrium Open system: to maintain non-equilibrium state reactants flow in and products flow out If the reactor is well stirred, we call it CSTR(continuous-flow stirred tank reactor) Molecules can spend different time at the CSTR because of flowing. The average time spent in reactor is called the mean residence time T res = volume/ flow rate. The system can display steady state, bistable state, oscillations, chaos in a CSTR. http://www.grc.org/programs.aspx?year=2 008&program=oscillat reactants stirrer

2
5.1 Steady state and bistability Example: Iodate - arsenite system IO 3 - +5I - +3H3AsO 3 → 6 I - +3 H3AsO 4 IO 3 - +5I - +3H3AsO 3 → 6 I - +3 H3AsO 4 At high flow rate, the concentrations of reactants in a CSTR have not bigger different with the inflow concentrations, this branch is known as flow branch( 流动分枝 ). At high flow rate, the concentrations of reactants in a CSTR have not bigger different with the inflow concentrations, this branch is known as flow branch( 流动分枝 ). At low flow rate, the system in a CSTR approach the thermodynamic equilibrium, the curve is called thermodynamic branch( 热力学分枝 ) At low flow rate, the system in a CSTR approach the thermodynamic equilibrium, the curve is called thermodynamic branch( 热力学分枝 ) As the flowrate is decreased, the system changes from flow branch to thermodynamic branch Attention: when flowrate change, the system evolve to new state for short(low flowrate) or long (high flow rate)time. Attention: when flowrate change, the system evolve to new state for short(low flowrate) or long (high flow rate)time. reacta nts stirrer

3
Situations: Situations: 1. When changing flowrate up or down, the plot of steady concentrations have only one curve. But during the some region of flowrate, there is a quick change of concentration from one branch to another branch. flow diagram

4
2. Thermodynamic branch and flow branch overlap as the direction of changing flowrate is opposite. This is phenomena is named bistabilty. Within the region of bistability, the actual state selected depends on not only the parameter but also the operational history. When the flowrate is up or down, the state jumps or fall off in delay. This phenomena is called hysteresis When the flowrate is up or down, the state jumps or fall off in delay. This phenomena is called hysteresis This abrupt changes in compositions is discontinuous responded to the continuous changes in the operating conditions. This is bifurcation. This abrupt changes in compositions is discontinuous responded to the continuous changes in the operating conditions. This is bifurcation.

8
5.4 Turning points and tangancies Line2,4 are tangential to R ， The system is bistable between line2 and line 4 Line2,4 are tangential to R ， The system is bistable between line2 and line 4 The condition for tangancy F=R dF/d α =dR/d α F=R dF/d α =dR/d α The stable concentration of α α ss ±=0.25{3±(1-8β0)0.5} To makeα ss real, β0<1/8 The flowrate of points of tangancies( Turning points) Phase diagram Phase diagram

9
5.5 Nodes and saddles: from bistable states to saddle-node bifurcation In bistable region, there are three steady states, (or three branchs), two are stable, the middle is unstable, How to understand it? For their stability, We use the potential rate=-dα/dt= α(1+β0-α) 2 -К0(1-α) rate=-dα/dt= α(1+β0-α) 2 -К0(1-α) V=1/2 m Rate 2 V=1/2 m Rate 2 dV/drate=Rate dV/drate=Rate V=∫Rate d(Rate) V=∫Rate d(Rate) V=0.5(1+β0)2α2-2/3(1+β0)α3+1/4α4- К 0 α(1-1/2α)+V0 At one parameter, Vα1, α3 has the minimum as node, and Vα2 has the maxmum as saddle. when the parameter changes to another, the maxmum as saddle. when the parameter changes to another, Vα3 merges with Vα2 from node to saddle. This process is called saddle-node bifurcation.

13
C To understand dynamics from nullclines β nullcline γ nullcline intersection situation have four possibilities. β nullcline γ nullcline intersection situation have four possibilities. Bistable high β low β oscillations Bistable high β low β oscillations Cross-shapes diagrams Parameter κ0 γ0 ‘ Cross-shapes diagrams Parameter κ0 γ0 ‘ a bistable d oscillations b high β c lowβ a bistable d oscillations b high β c lowβ near cusp perturbation oscllations or go to stable state near cusp perturbation oscllations or go to stable state If b c display oscillations, then a is the field of birhythmicity and d is the field of complex If b c display oscillations, then a is the field of birhythmicity and d is the field of complex oscillations or chaos oscillations or chaos