Critical reactors, computing

We also described how concrete equations for critical DSR and CSTRs may be computed. These expressions are complicated to compute analytically, which are derived from geometric controllability arguments developed by Feinberg (2000a, 2000b). These conditions are intricate, and thus it is often not possible to compute analytic solutions to the equations that describe critical reactors. For three-dimensional systems, a shortcut method involving the vDelR condition may be used to find critical a policies. Irrespective of the method used, the conditions for critical reactors are well defined, irrespective of the legitimacy of the kinetics studied, and thus these conditions must be enforced if we wish to attain points on the true AR boundary. 

The AR is composed of mixing lines and manifolds of PFR trajectories. The final approach to the extreme points of the AR boundary is achieved using PFR solution trajectories—if a desired operating point resides on the AR boundary, a PFR must be incorporated into the reactor structure in order to reach it, and thus PFRs are often the best terminating reactor to use in practice (for any kinetics and feed point). Only combinations of PFRs, CSTRs, and DSRs are required to form the AR. This result is true for all dimensions. Distinct expressions may be derived to compute critical a policies for the DSR profile and critical CSTR residence times. These expressions are intricate and complex in nature, which are ultimately based on the lack of controllability in a critical reactor. This idea is important in understanding the nature of the AR and how to achieve points on the true AR boundary. 

The flux distributions for the case ife = 0.2, tpR = 0, and B — 0.8 are shown in Fig. 8.19. It is of interest to note the discontinuity in the flux due to the application of the Serber-Wilson condition. Also, note that the critical mass computed by the continuity-of-flux condition exceeds the Serber-Wilson estimate by as much as 40 per cent for small reactors. 

Alternative procedures for determining equivalent systems is given by F. G. Prohammer, A Comparison of One-dimensional Critical Mass Computations for Completely Reflected Reactors, ORNL-2007, Mar. 1, 1956. 

Previous computations (189) show that the critical value of Rat for non-Boussinesq conditions is approximately the same as that for a Boussinesq fluid in a box heated from below, at least when H2 is the carrier gas. Thus, results from the stability analysis of the classical Rayleigh-Benard problem of a two-dimensional fluid layer heated from below (see reference 190 for a review) may be used to indicate the type of behavior to be expected in a horizontal reactor with insulated side walls. As anticipated from this analysis, an increase in the reactor height from 2 to 4 cm raises the value of Rat to 4768, which is beyond the stability limit, Rat critical = 2056, for a box of aspect ratio 2 (188). The trajectories show the development of buoyancy-driven axial rolls that are symmetric about the midplane and rotating inward. For larger values of Rat (>6000), transitions to three-dimensional or time-de-... 

This is the most common mode of addition. For safety or selectivity critical reactions, it is important to guarantee the feed rate by a control system. Here instruments such as orifice, volumetric pumps, control valves, and more sophisticated systems based on weight (of the reactor and/or of the feed tank) are commonly used. The feed rate is an essential parameter in the design of a semi-batch reactor. It may affect the chemical selectivity, and certainly affects the temperature control, the safety, and of course the economy of the process. The effect of feed rate on heat release rate and accumulation is shown in the example of an irreversible second-order reaction in Figure 7.8. The measurements made in a reaction calorimeter show the effect of three different feed rates on the heat release rate and on the accumulation of non-converted reactant computed on the basis of the thermal conversion. For such a case, the feed rate may be adapted to both safety constraints the maximum heat release rate must be lower than the cooling capacity of the industrial reactor and the maximum accumulation should remain below the maximum allowed accumulation with respect to MTSR. Thus, reaction calorimetry is a powerful tool for optimizing the feed rate for scale-up purposes [3, 11]. 

In theory, by feeding the MWD and experimental rate data into a mathematical model containing a variety of polymerization mechanisms, it should be possible to find the mechanism which explains all the experimental phenomena and to evaluate any unknown rate constants. As pointed out by Zeman (58), as long as there are more independent experimental observations than rate parameters, the solution should, in principle, be unique. This approach involves critical problems in choice of experiments and in experimental as well as computational techniques. We are not aware of its having yet been successfully employed. The converse— namely, predicting MWD from different reactor types on the basis of mathematical models and kinetic data—has been successfully demonstrated, however, as discussed above. The recent series of interesting papers by Hamielec et al. is a case in point. 

At the other extreme, it may be argued that the traditional low-dimensional models of reactors (such as the CSTR, PFR, etc.) should be abandoned in favor of the detailed models of these systems and numerical solution of the full convection-diffusion reaction (CDR) equations using computational fluid dynamics (CFD). While this approach is certainly feasible (at least for singlephase systems) due to the recent availability of computational power and tools, it may be computationally prohibitive, especially for multi-phase systems with complex chemistry. It is also not practical when design, control and optimization of the reactor or the process is of main interest. The two main drawbacks/criticisms of this approach are (i) It leads to discrete models of very high dimension that are difficult to incorporate into design and control schemes. 

The modelling of super critical water oxidation (SCWO), up to now not been used in large scale industrial applications, is important for design of pilot plants and, later, industrial plants. The applied programme to model the continuous flow in a reactor is called CAST (Computer Aided Simulation of Turbulent Flows [8]) and is based on the method of the finite volume. That means that the balance equations were integrated over the surfaces of each control volume. 

Despite these limitations computer simulations have become a key method in dealing with the various critical issues related to plasma edge phenomena on the path towards economical fusion power reactors, by quantifying at least the known parts of edge science, in a most detailed and complete way possible today. 

In critical cases it may well be worthwhile to make a complete analysis of stability. In many cases, however, enough can be learned by studying what Bilous and Amundson (B7) called parametric sensitivity. These authors derived formulas for calculating the amplification or attenuation of disturbances imposed on an unpacked tubular reactor originally in a steady state, with the idea that if the disturbances grow unduly the performance of the reactor is too sensitive to the conditions imposed on it, that is, to the parameters of the system. The effect of feedback from a control system was not considered. As pointed out by the authors, it would be a much more complicated task to apply their procedure to a packed reactor, but it still would entail far less computation than a study of the transient response. 

Fluidized bed reactors have received increased interest in recent years owing to their application in coal gasification. The section on fluidized beds discusses critical areas in fluid bed reactor modeling. Computer simulation of both solid-catalyzed gas phase reactions as well as gas-solid reactions are included. 

It was suggested that computer-based data analysis techniques (often involving multivariate statistical methods) can aid in this classification or simplification, as has been so profitable in other thermochemical conversion endeavours, for example, as applied to coal and petroleum. Again, it was emphasized that there is a need for a critical synthesis of the wealth of experimental data into regimes of behaviour, and simpler predictive equations or simulations, that are useful to the technologists in industry who are designing industrial scale reactors. 

Nowadays, improved computing facilities and, more importantly, the availability of the Chemkin package (Kee and Rupley, 1990) and similar kinetic compilers and processors have made these complex kinetic schemes more user-friendly and allows the study of process alternatives as well as the design and optimization of pyrolysis coils and furnaces. In spite of their rigorous, theoretical approach, these kinetic models of pyrolysis have always been designed and used for practical applications, such as process simulation, feedstock evaluations, process alternative analysis, reactor design and optimization, process control and so on. Despite criticisms raised recently by Miller et al. (2005), these detailed chemical kinetic models constitute an excellent tool for the analysis and understanding of the chemistry of such systems. 

Rieck [Rl] has used the computer codes LEOPARD IBl] and SIMULATE [FI] to predict the power distribution in the fuel and poison arrangement shown in Fig. 3.19 for the first fuel cycle for this reactor, and the amount of thermal energy produced by each assembly up to the time when the reactor ceases to be critical with all soluble boron removed from the cooling water. Figure 3.20 is a horizontal cross section of one-quarter of the core of this reactor. Each square represents one fuel assembly. The core arrangement has 90° rotational symmetry, about the central assembly 1AA at the upper left of the figure. 

Let us compute the critical cluster size and the number of NH4CI molecules in a specific case. Countess and Heicklen (1973) studied the growth of NH4CI particles from the reaction of NH3 and HCl in a flow reactor. The experimental procedure involved mixing 60 ppm of NH 3 and 60 ppm of HCl in 1 atm of nitrogen at 293 K. The critical cluster size can be computed from (10.128). First, we need to know Kp at 293 K. We may use... 

The requirements of ANSI/ANS 8.1 specify that calculational methods for away-from-reactor criticality safety analyses be validated against experimental measurements. If credit is to be taken for the reduced reactivity of burned or spent fuel relative to its original fresh composition, it is necessary to benchmark computational methods used in determining such reactivity worth against spent fuel reactivity measurements. This report summarizes a portion of the ongoing effort to benchmark away-from-reactor criticality analysis methods using critical configurations from commercial pressurized- water reactors (PWR). 

---