Time trajectories

Abstract. A stochastic path integral is used to obtain approximate long time trajectories with an almost arbitrary time step. A detailed description of the formalism is provided and an extension that enables the calculations of transition rates is discussed. 

In accordance with the one-dimensional periodic orbit theory, any orbit contributing to g E) is supposedly constructed from closed classical orbits in the well and subbarrier imaginary-time trajectories. These two classes of trajectories are bordering on the turning points. For the present model the classical motion in the well is separable, and the harmonic approximation for classical motion is quite reasonable for more realistic potentials, if only relatively low energy levels are involved. 

In many simulations that use a time trajectory to sample membrane properties, it is the equilibrium situation that one is interested in, rather than the dynamics themselves. The dynamics are then just a by-product that is only used to judge the degree of equilibration. 

It was interesting to find that the branching ratio depends on the elapsed time trajectories were more selective when they reached product region faster. Thus, the fate of the trajectories is primarily determined by dynamics effect, or momentum, and the branching ratio depends on the shape of the surface. 

If, for example, the reactor temperature is disturbed from 2/3, the high temperature steady state, to y4 that is not a steady-state temperature, then the static diagram helps us to determine the direction of temperature change towards a steady state of the system as follows. At y4 the heat generation exceeds the heat removal since G(y4) > R(y4) from the graph. Therefore, the temperature will increase. A quantitative analysis of the temperature-time trajectory can of course only be determined from a dynamic model of the system. 

The classical batch reactor is a perfectly mixed vessel in which reactants are converted to products during the course of a batch cycle. All variables change dynamically with time. The reactants are charged into the vessel. Heat and/or catalyst is added to initiate reaction. Reactant concentrations decrease and product concentrations increase with time. Temperature or pressure is controlled according to some desired time trajectory. Batch time is also a design and operating variable, which has a strong impact on productivity. 

Putting these important issues aside, the production of ethanol by batch fermentation is an important example of a batch reactor. The basic regulatory control of a batch ethanol fermentor is not a difficult problem because the heat removal requirements are modest and there is no need for very intense mixing. In this section we develop a very simple dynamic model and present the predicted time trajectories of the important variables such as the concentrations of the cells, ethanol, and glucose. The expert advice of Bjom Tyreus of DuPont is gratefully acknowledged. Sources of models and parameter values are taken from three publications.1 3... 

However, high temperatures mean large reaction rates and high heat removal requirements that can exceed the heat transfer capacity of the condenser. Thus there is an optimum pressure-time trajectory that maximizes yield while not violating heat removal constraints. The normal ramp rate of the pressure controller setpoint is about0.il atm/min. 

The variability of the sizes of the solid sodium-lead alloy particles introduces another factor. The smaller the particle size, the more surface area is exposed to reaction, which increases the effective reaction rate. Big particles react slowly small particles react quickly. So, for the same pressure-time trajectory, significantly different heat load-time requirements can occur. If the alloy charge contains very small particles, a runaway reaction can occur if the pressure controller setpoint is ramped up too quickly. 

Dynamic surface tension is the time trajectory of surface tension before equilibrium is reached. Dynamic surface tension tracks the changes during surface formation when surfactants are added. The bubble pressure method is the one most commonly used for the determination of dynamic surface tension. The details of this method are described in ASTM D3825-90 (2000) [ 19]. In this method a capillary tube is immersed in a sample liquid and a constant flow of gas is maintained through the tube forming bubbles in the sample liquids. The surface tension of the sample is calculated from the pressure difference inside and outside the bubble and the radius of the bubble. 

Fig. 29. Origin of systematic errors in spite of potentially error-free analysis. On-line sampling setups (top) and time trajectories of limiting substrate concentration during sample preparation in the two paradigmatic setups depending on the actual culture density (bottom). Either a filter in bypass loop is used for the preparation of cell-free supernatant (upper part in top insert) or an aliquot of the entire culture is removed using an automatic sampler valve and a sample bus for further inactivation and transport of the samples taken (lower part). Both methods require some finite time for sample transportation from the reactor outlet (at z = 0) to the location where separation of cells from supernatant or inactivation by adding appropriate inactivators (at z = L) takes place. During transport from z = 0 to z = L, the cells do not stop consuming substrate. A low substrate concentration in the reactor (namely s KS) and a maximal specific substrate consumption rate of 3 g g h 1 were assumed in the simulation example to reflect the situation of either a fed-batch or a continuous culture of an industrially relevant organism such as yeast. The actual culture density (in g 1 1) marks some trajectories in the mesh plot. Note that the time scale is in seconds...

Fig. 30. Schematic design of a simple but very useful and efficient data reduction algorithm. Data representing the time trajectory of an individual variable are only kept (= recorded, stored) when the value leaves a permissive window which is centered around the last stored value. If this happens, the new value is appended to the data matrix and the window is re-centered around this value. This creates a two-column matrix for each individual variable with the typical time stamps in the first column and the measured (or calculated) values in the second column. In addition, the window width must be stored since it is typical for an individual variable. This algorithm assures that no storage space is wasted whenever the variable behaves as a parameter (i.e. does not change significantly with time, is almost constant) but also assures that any rapid and/or singular dynamic behavior is fully documented. No important information is then lost...

Fig. 7.31 Position-time trajectories for the compression and release waves.

Figure E7-2.1 shows a comparison of the concentration-time trajectory for ethane oaicnlated from the PSSH (CPI) with the ethane trajectory (Cl) calculated from solving the mole balance Equations (E7-2.14) through (E7-2.20). Figure E7-2.2 shows a similar comparison for ethylene (CPS) and (C5). One notes that the curves are identical, indicating the validity of the PSSH under these conditions. Figure E7-2.3 shows a comparison the concentration-time trajectories for methane (C3) and butane (C8). Problem P7-2(a) explores the temperature for which the PSSH is valid for the cracking of ethane.

Figure E7-2.3 Comparison of concentration-time trajectories for aumhane (C3)...

Figure E7-2.I Comparison of concentraiion-time trajectories for ethane.

Figure E7-2.2 Comparison for temperature-time trajectory for ethylene.

Adapted from the problem by Ronald Willey, Seminar on a Nitroanaline Reactor Rupture. Prepared for SAGHE, Center for Chemical Process Safety, American Institute of Chemical Engineers, New York (1994). The values of and UA were estimated in the plant data of the temperature-time trajectory in the article by G. C. Vincent, Loss Prevention, Vol. 5. p. 46-52, AIChE, New York NY. 

Plot the temperature-time trajectory up to a period of 120 min after the reactants were mixed and brought up to 175°C. Showthat the following three conditions had to have been present for the explosion to occur (1) increased ONCB charge, (2) reactor stopped for 10 min, and (3) relief system failure. 

Return of the cooling occurs at 55 min. The values at the end of the period of adiabatic operation (T = 468 K, X = 0.04.23) become the initial conditions for the period of operation with heat exchange. The cooling is turned on at its maximum capacity, Q = UA(29S — T),at 55 min. Table E9-2.1 gives the POLYMATH program to determine the temperature-time trajectory. 

The complete temperature-time trajectory is shown in Figure E9-2.2. One notes the long plateau after the cooling is turned back on. Using the values of Qg and 2, at 55 min and substituting into Equation (E9-2.S), we find that... 

Figure E9-3.2 Concentration-time trajectories in a semibatch reactor.

Slow decay - Temperature-Time Trajectories (10.7.2) Moderate decay -Moving-Bed Reactors (10.7.3)... 

---