Reactors nonlinearity

Chicone, C. and Retzloff, D. G., 1981, Dynamics of the CR equations modelling a constant flow stirred tank reactor, Nonlinear Anal. 6, 983-1000. 

Research Areas Modeling, Simulation and Optimization of Chemical and Biological Processes, Clean Fuels (Hydrogen, Biodiesel and Ethanol), Fixed and Fluidized Bed Catalytic Reactors, Nonlinear Dynamics, Bifurcation and Chaos,... 

Theorem 6 has a certain appeal because the monotonic character of the reactor nonlinearity, exp(x) — 1, is taken into account in the proof. The Popov-Kalman result is of course valid for a much larger class of nonlinearities. 

First, we do think that future work on Eqs. (1) should take the particular features of the reactor nonlinearity into consideration more than has so far generally been done. Results worked out for more general differential or integral equations should not be blindly applied. Instead one should make more efficient use of the monotonicity and boundedness from below of exp(x) — 1. This may sound trivial but there are far too many papers in reactor dynamics which only try to apply some more general criteria to Eqs. (1) and which do not take the special features of these equations into account. Usually this approach produces overly restrictive results. 

Another important reaction supporting nonlinear behaviour is the so-called FIS system, which involves a modification of the iodate-sulfite (Landolt) system by addition of ferrocyanide ion. The Landolt system alone supports bistability in a CSTR the addition of an extra feedback chaimel leads to an oscillatory system in a flow reactor. (This is a general and powerfiil technique, exploiting a feature known as the cross-shaped diagram , that has led to the design of the majority of known solution-phase oscillatory systems in flow... 

Consequently, when D /Dj exceeds the critical value, close to the bifurcation one expects to see the appearance of chemical patterns with characteristic lengtli i= In / k. Beyond the bifurcation point a band of wave numbers is unstable and the nature of the pattern selected (spots, stripes, etc.) depends on the nonlinearity and requires a more detailed analysis. Chemical Turing patterns were observed in the chlorite-iodide-malonic acid (CIMA) system in a gel reactor [M, 59 and 60]. Figure C3.6.12(a) shows an experimental CIMA Turing spot pattern [59]. 

Computer simulation of the reactor kinetic hydrodynamic and transport characteristics reduces dependence on phenomenological representations and idealized models and provides visual representations of reactor performance. Modem quantitative representations of laminar and turbulent flows are combined with finite difference algorithms and other advanced mathematical methods to solve coupled nonlinear differential equations. The speed and reduced cost of computation, and the increased cost of laboratory experimentation, make the former increasingly usehil. 

The second classification is the physical model. Examples are the rigorous modiiles found in chemical-process simulators. In sequential modular simulators, distillation and kinetic reactors are two important examples. Compared to relational models, physical models purport to represent the ac tual material, energy, equilibrium, and rate processes present in the unit. They rarely, however, include any equipment constraints as part of the model. Despite their complexity, adjustable parameters oearing some relation to theoiy (e.g., tray efficiency) are required such that the output is properly related to the input and specifications. These modds provide more accurate predictions of output based on input and specifications. However, the interactions between the model parameters and database parameters compromise the relationships between input and output. The nonlinearities of equipment performance are not included and, consequently, significant extrapolations result in large errors. Despite their greater complexity, they should be considered to be approximate as well. 

This matrix will contain information regarding loading characteristics such as flooding hmits, exchanger areas, pump curves, reactor volumes, and the like. While this matrix may be adjusted during the course of model development, it is a boundary on any possible interpretation of the measurements. For example, distillation-column performance markedly deteriorates as flood is approached. Flooding represents a boundary. These boundaries and nonlinearities in equipment performance must be accounted for. 

The sum of squares as defined by Equation 7.8 is the general form for the objective function in nonlinear regression. Measurements are made. Models are postulated. Optimization techniques are used to adjust the model parameters so that the sum-of-squares is minimized. There is no requirement that the model represent a simple reactor such as a CSTR or isothermal PER. If necessary, the model could represent a nonisothermal PFR with variable physical properties. It could be one of the distributed parameter models in Chapters 8 or 9. The model... 

Almost all flows in chemical reactors are turbulent and traditionally turbulence is seen as random fluctuations in velocity. A better view is to recognize the structure of turbulence. The large turbulent eddies are about the size of the width of the impeller blades in a stirred tank reactor and about 1/10 of the pipe diameter in pipe flows. These large turbulent eddies have a lifetime of some tens of milliseconds. Use of averaged turbulent properties is only valid for linear processes while all nonlinear phenomena are sensitive to the details in the process. Mixing coupled with fast chemical reactions, coalescence and breakup of bubbles and drops, and nucleation in crystallization is a phenomenon that is affected by the turbulent structure. Either a resolution of the turbulent fluctuations or some measure of the distribution of the turbulent properties is required in order to obtain accurate predictions. 

Parameter estimation. Integral reactor behavior was used for the interpretation of the experimental data, using N2O conversion levels up to 70%. The temperature dependency of the rate parameters was expressed in the Arrhenius form. The apparent rate parameters have been estimated by nonlinear least-squares methods, minimizing the sum of squares of the residual N2O conversion. Transport limitations could be neglected. 

The methods concerned with differential equation parameter estimation are, of course, the ones of most concern in this book. Generally reactor models are non-linear in their parameters and therefore we are concerned mostly with nonlinear systems. 

Rippin, D.W.T., L.M. Rose and C. Schifferli. "Nonlinear Experimental Design with Approximate Models in Reactor Studies for Process Development", Chem. Eng. 5c/., 35, 356 (1980). 

In the operation of BWRs, especially when operating near the threshold of instability, the stability margin of the stable system and the amplitude of the limit cycle under unstable condition become of importance. A number of nonlinear dynamic studies of BWRs have been reported, notably in an International Workshop on Boiling Water Reactor Stability (1990). The following references are mentioned for further study. 

Abstract Acoustic cavitation is the formation and collapse of bubbles in liquid irradiated by intense ultrasound. The speed of the bubble collapse sometimes reaches the sound velocity in the liquid. Accordingly, the bubble collapse becomes a quasi-adiabatic process. The temperature and pressure inside a bubble increase to thousands of Kelvin and thousands of bars, respectively. As a result, water vapor and oxygen, if present, are dissociated inside a bubble and oxidants such as OH, O, and H2O2 are produced, which is called sonochemical reactions. The pulsation of active bubbles is intrinsically nonlinear. In the present review, fundamentals of acoustic cavitation, sonochemistry, and acoustic fields in sonochemical reactors have been discussed. 

A control algorithm has been derived that has improved the dynamic decoupling of the two outputs MW and S while maintaining a minimum "cost of operation" at the steady state. This algorithm combines precompensation on the flow rate to the reactor with state variable feedback to improve the overall speed of response. Although based on the linearized model, the algorithm has been demonstrated to work well for the nonlinear reactor model. 

---