Simulation with isothermal conditions

The tube simulation can be ran under either isothermal conditions (by excluding Equation (6.24) from the numerical integration and using a pre-defined field of temperatures) or adiabatic conditions (in this case, Equation (6.24) is used to evaluate the temperature profile). In the isothermal case, it is possible to simulate the electrochemical performance (in terms of I-V curves) at different imposed operational temperatures the option of running the simulation with isothermal sohd temperature distribution is often apphed, since this is the condition under which some experimental data are obtained at RRFCS (see Figure 6.11). All the other results reported in this section have been obtained under adiabatic conditions (i.e. perfect insulation of the vessel where the tube is contained), since this is a realistic practical operating condition for the tube when included into the plant. 

As shown in Figure 5.23, the concentration profile of the precursor gas is related to the value of Dah The concentration gradient of the precursor gas becomes steeper with increasing Dah Due to complex coupling effects between the pressure and temperature of F-CVI it is very difficult to model this phenomenon. A large body of research work has been undertaken under isothermal conditions to simplify the simulation. Under isothermal conditions the concentration profile represents the deposition gradient. In such a case densification always occurs more rapidly at the precursor gas entrance region of the preform. 

Equilibrium Compositions for Single Reactions. We turn now to the problem of calculating the equilibrium composition for a single, homogeneous reaction. The most direct way of estimating equilibrium compositions is by simulating the reaction. Set the desired initial conditions and simulate an isothermal, constant-pressure, batch reaction. If the simulation is accurate, a real reaction could follow the same trajectory of composition versus time to approach equilibrium, but an accurate simulation is unnecessary. The solution can use the method of false transients. The rate equation must have a functional form consistent with the functional form of K,i,ermo> e.g., Equation (7.38). The time scale is unimportant and even the functional forms for the forward and reverse reactions have some latitude, as will be illustrated in the following example. 

The reactions are assumed to take place under isothermal conditions at 130°C at 10 bar. The liquid feed of BA is 0.0133 kmol-s 1 and the gaseous feed of chlorine is 0.1 kmol-s-1. The objective is to maximize the fractional yield of a-monochlorobutanoic acid with respect to butanoic acid. Specialized software is required to perform the calculations, in this case, using simulated annealing. 

Sucessful simulations have been performed with computerized fluid dynamics programs (CFD), based on the fundamental Navier-Stokes equations, with appropriate volume element grids and/or a finite elemet approach, and in some cases were backed up by isothermal physical modeling (hydraulic modeling) experiments [455]. Examples for the CFD software used are FLUENT (Fluent Inc.) [450] and CFDS-FLOW3D [458] and others, usually modified by the contractor or licensor [448] to adapt them to the conditions in a secondary reformer. Discussion of simulation with CFD can be found in [444], [448], [449]. 

The initial state of the simulations consisted of RDX perfect crystals using simulation cells containing 8 molecules (one unit cell, 168 atoms) and 3D periodic conditions. After relaxing the atomic positions at each density with low temperature MD, we studied the time evolution of the system at the desired temperature with isothermal isochoric (NVT ensemble) MD simulations (using a Berendsen thermostat the relaxation time-scale associated with the coupling between the thermostat and the atomistic system was 200 femtoseconds). 

In order to simulate the concentration-time data and evaluate the kinetic parameters in equations 6 - 8, a batch reactor model was used. Under isothermal conditions, the variation of the liquid phase concentration of crotonaldehyde and butyraldehyde with time can be described by the following set of mass balance equations. 

The rubber used was the EPDM compound described in Table 3.2 with 2% peroxide. The Monsanto 2000 E MDR is used under isothermal conditions, while the simulation is done by calculation. The thickness of the sample is taken as equal to 0.2 cm. A calorimeter C 80 was used in scanning mode with a heating rate of 0.2°C/min and a sample of around 6 g. 

It is possible to assume that kinetic and hydrodynamic methods of reactors analysis and design are well advanced at present. Methods of computer simulation and modelling are widely used. So, we can say that if we know processes kinetic and hydrodynamic parameters and fundamental particularities of reactor functioning we can calculate all process characteristics and its stmcture, we also can predict effectiveness of apparatus operation and consumer properties of chemical production. Meanwhile criterion function development for calculation and organization of novel processes and optimization of present productions require overcoming of a big number of problems in production practice. This principle is satisfactory enough for processes with low or medium reactions rates, when creation of isothermal conditions in apparatus is easy. In this case it is easy to calculate and reproduce in working conditions all characteristics of chemical process and to control the last ones. 

Reports are also available on CO2 selective membrane reactors for WGS reaction. Zou et al. [40] first time synthesized polymeric C02-selective membrane by incorporating fixed and mobile carriers in cross-linked poly vinyl alcohol. Micro-porous Teflon was used as support. They used Cu0/Zn0/Al203 catalyst for low temperature WGS reaction. They investigated the effect of water content on the CO2 selectivity and CO2/H2 selectivity. As the water concentration in the sweep gas increased, both CO2 permeability and CO2/H2 selectivity increased significantly. Figure 6.18 shows the influence of temperature on CO2 permeability and CO2/H2 selectivity. Both CO2 permeability and CO2/ H2 selectivity decrease with increasing reactimi temperature. After the catalyst activation, the synthesis gas feed containing 1% CO, 17% CO2, 45% H2 and 37% N2 was pumped into the membrane reactor. They are able to achieve almost 100% CO conversion. They also developed a one-dimensional non-isothermal model to simulate the simultaneous reaction and transport process and verified the model experimentally under an isothermal condition. 

Kinetic and hydrodynamic analyses, and methods for the calculation of the parameters of industrial reactors are sufficiently developed today [2-6]. Computer simulation is also popular because if we know the kinetic and hydrodynamic parameters of processes and the principles of reactor behaviour, it is not a problem to calculate process characteristics and final product performance. This principle is an adequate tool for the description of low and medium rate chemical transformations with uniform concentration fields and isothermic conditions which are easy to achieve. In this case, it is easy to calculate and control all the characteristics of a chemical process under real conditions. 

---