Computational mixing models

The algorithms of the mixed classical-quantum model used in HyperChem are different for semi-empirical and ab mi/io methods. The semi-empirical methods in HyperChem treat boundary atoms (atoms that are used to terminate a subset quantum mechanical region inside a single molecule) as specially parameterized pseudofluorine atoms. However, HyperChem will not carry on mixed model calculations, using ab initio quantum mechanical methods, if there are any boundary atoms in the molecular system. Thus, if you would like to compute a wavefunction for only a portion of a molecular system using ab initio methods, you must select single or multiple isolated molecules as your selected quantum mechanical region, without any boundary atoms. 

ZINDO/S is different from ZINDO/I because they use different algorithms in computing the Coulomb integrals. Hence the two equations used in the mixed model in ZINDO/1 are also employed... 

The mixed solvent models, where the first solvation sphere is accounted for by including a number of solvent molecules, implicitly include the solute-solvent cavity/ dispersion terms, although the corresponding tenns between the solvent molecules and the continuum are usually neglected. Once discrete solvent molecules are included, however, the problem of configuration sampling arises. Nevertheless, in many cases the first solvation shell is by far the most important, and mixed models may yield substantially better results than pure continuum models, at the price of an increase in computational cost. 

Airlift loop reactor (ALR), basically a specially structured bubble column, has been widely used in chemical industry, biotechnology and environmental protection, due to its high efficiency in mixing, mass transfer, heat transfer etc [1]. In these processes, multiple reactions are commonly involved, in addition to their complicated aspects of mixing, mass transfer, and heat transfer. The interaction of all these obviously affects selectivity of the desired products [2]. It is, therefore, essential to develop efficient computational flow models to reveal more about such a complicated process and to facilitate design and scale up tasks of the reactor. However, in the past decades, most involved studies were usually carried out in air-water system and the assumed reactor constructions were oversimplified which kept itself far away from the real industrial conditions [3] [4]. 

The classical CRE model for a perfectly macromixed reactor is the continuous stirred tank reactor (CSTR). Thus, to fix our ideas, let us consider a stirred tank with two inlet streams and one outlet stream. The CFD model for this system would compute the flow field inside of the stirred tank given the inlet flow velocities and concentrations, the geometry of the reactor (including baffles and impellers), and the angular velocity of the stirrer. For liquid-phase flow with uniform density, the CFD model for the flow field can be developed independently from the mixing model. For simplicity, we will consider this case. Nevertheless, the SGS models are easily extendable to flows with variable density. 

Before discussing in detail specific molecular mixing models, it is useful to first state a few important constraints that can be derived by computing expected values. The first constraint follows from45... 

Complete Mix Reactor - The complete mix reactor is also labeled a completely stirred tank reactor. It is a container that has an inhnite diffusion coefficient, such that any chemical that enters the reactor is immediately mixed in with the solvent. In Example 2.8, we used the complete mix reactor assumption to estimate the concentration of three atmospheric pollutants that resulted from an oil spill. We will use a complete mix reactor (in this chapter) to simulate the development of high salt content in dead-end lakes. A series of complete mix reactors may be placed in series to simulate the overall mixing of a one-dimensional system, such as a river. In fact, most computational transport models are a series of complete mix reactors. 

Wolfinger RD, Gibson G, Wollinger ED, Bennett L, Hamadeh H, Bushel P, Afshari C, Paules RS. Assessing gene significance from cDNA microarray expression data via mixed models. J Comput Biol 2001 8 625-637. 

In between the implicit and explicit solvent models, there are mixed models, such as the solvation shell approximation.67-69 This model describes explicitly only the first solvation shell molecules and treats as implicit the solvent region beyond the first solvation shell. Such treatment both provides the information about the solvent structure near the solute and allows for faster computation. 

Rawlings and co-workers proposed to carry out parameter estimation using Newton s method, where the gradient can be cast in terms of the sensitivity of the mean (Haseltine, 2005). Estimation of one parameter in kinetic, well-mixed models showed that convergence was attained within a few iterations. As expected, the parameter values fluctuate around some average values once convergence has been reached. Finally, since control problems can also be formulated as minimization of a cost function over a control horizon, it was also suggested to use Newton s method with relatively smooth sensitivities to accomplish this task. The proposed method results in short computational times, and if local optimization is desired, it could be very useful. 

For the computation of E(t) and T(p), in the case of a perfect mixing model, we use the representation and notation given in Fig. 3.24. Including the mass balance of the species in the signal, we derive the following differential equation ... 

Let us consider a shallow fluidized bed combustor with multiple coal feeders which are used to reduce the lateral concentration gradient of coal (11). For simplicity, let us assume that the bed can be divided into N similar cylinders of radius R, each with a single feed point in the center. The assumption allows us to use the symmetrical properties of a cylindrical coordinate system and thus greatly reduce the difficulty of computation. The model proposed is based on the two phase theory of fluidization. Both diffusion and reaction resistances in combustion are considered, and the particle size distribution of coal is taken into account also. The assumptions of the model are (a) The bed consists of two phases, namely, the bubble and emulsion phases. The voidage of emulsion phase remains constant and is equal to that at incipient fluidization, and the flow of gas through the bed in excess of minimum fluidization passes through the bed in the form of bubbles (12). (b) The emulsion phase is well mixed in the axial... 

The most expensive part of a simulation of a system with explicit solvent is the computation of the long-range interactions because this scales as Consequently, a model that represents the solvent properties implicitly will considerably reduce the number of degrees of freedom of the system and thus also the computational cost. A variety of implicit water models has been developed for molecular simulations [56-60]. Explicit solvent can be replaced by a dipole-lattice model representation [60] or a continuum Poisson-Boltzmann approach [61], or less accurately, by a generalised Bom (GB) method [62] or semi-empirical model based on solvent accessible surface area [59]. Thermodynamic properties can often be well represented by such models, but dynamic properties suffer from the implicit representation. The molecular nature of the first hydration shell is important for some systems, and consequently, mixed models have been proposed, in which the solute is immersed in an explicit solvent sphere or shell surrounded by an implicit solvent continuum. A boundary potential is added that takes into account the influence of the van der Waals and the electrostatic interactions [63-67]. 

These models require information about mean velocity and the turbulence field within the stirred vessels. Computational flow models can be developed to provide such fluid dynamic information required by the reactor models. Although in principle, it is possible to solve the population balance model equations within the CFM framework, a simplified compartment-mixing model may be adequate to simulate an industrial reactor. In this approach, a CFD model is developed to establish the relationship between reactor hardware and the resulting fluid dynamics. This information is used by a relatively simple, compartment-mixing model coupled with a population balance model (Vivaldo-Lima et al., 1998). The approach is shown schematically in Fig. 9.2. Detailed polymerization kinetics can be included. Vivaldo-Lima et a/. (1998) have successfully used such an approach to predict particle size distribution (PSD) of the product polymer. Their two-compartment model was able to capture the bi-modal behavior observed in the experimental PSD data. After adequate validation, such a computational model can be used to optimize reactor configuration and operation to enhance reactor performance. 

Translating batch data for continuous reactors In most cases, laboratory-and bench-scale experiments required to validate the reactor concept are carried out in batch mode. It is then necessary to translate (or to use) the data obtained in these experiments to design continuous reactors. The location of feed pipes, outlets and their influence on mixing and performance needs to be understood. Computational flow models can be of great help in this regard. 

Tubular reactors are also used to carry out some multiphase reactions. Wamecke et al. (1999) reported use of a computational flow model to simulate an industrial tubular reactor carrying out a gas-liquid reaction (propylene oxide manufacturing process). In this process, liquid is a dispersed phase and gas is a continuous phase. The two-fluid model discussed earlier may be used to carry out simulations of gas-liquid flow through a tubular reactor. Warnecke et al. (1999) applied such a model to evaluate the influence of bends etc. on flow distribution and reactor performance. The model may be used to evolve better reactor configurations. In many tubular reactors, static mixers are employed to enhance mixing and other transport processes. Computational flow models can also make significant contributions to understanding the role of static mixers and for their optimization. Visser et al. (1999) reported CFD... 

---