Simulation of Dynamic Models

Simulation of Dynamic Models Linear dynamic models are particularly useful for analyzing control system behavior. The insight gained through linear analysis is invaluable. However, accurate dynamic process models can involve large sets of nonlinear equations. Analytical solution of these models is not possible. Thus, in these cases, one must turn to simulation approaches to study process dynamics and the effect of process control. Equation (8-3) will be used to illustrate the simulation of nonlinear processes. If dcjdt on the left-hand side of Eq. (8-3) is replaced with its finite difference approximation, one gets 

Starting with an initial value of cA and given c,(t), Eq. (8-4) can be solved for cA(t + At). Once cA(t + At) is known, the solution process can be repeated to calculate cA(t + 2At), and so on. This approach is called the Euler integration method while it is simple, it is not necessarily the best approach to numerically integrating nonlinear differential equations. As discussed in Sec. 3, more sophisticated approaches are available that allow much larger step sizes to be taken but require additional calculations. One widely used approach is the fourth-order Runge Kutta method, which involves the following calculations  

To illustrate how Laplace transforms work, consider the problem of solving Eq. (8-2), subject to the initial condition that cA = c, = 0 at t = 0. If cA were not initially zero, one would define a deviation variable between cA and its initial value cA0. Then the transfer function would be developed by using this deviation variable. If c, changes from zero to r, taking the Laplace transform of both sides of Eq. (8-2) gives 

By using the entries in Table 8-1, Eq. (8-13) can be inverted to give the transient response of cA as 

---

These apparent restrictions in size and length of simulation time of the fully quantum-mechanical methods or molecular-dynamics methods with continuous degrees of freedom in real space are the basic reason why the direct simulation of lattice models of the Ising type or of solid-on-solid type is still the most popular technique to simulate crystal growth processes. Consequently, a substantial part of this article will deal with scientific problems on those time and length scales which are simultaneously accessible by the experimental STM methods on one hand and by Monte Carlo lattice simulations on the other hand. Even these methods, however, are too microscopic to incorporate the boundary conditions from the laboratory set-up into the models in a reahstic way. Therefore one uses phenomenological models of the phase-field or sharp-interface type, and finally even finite-element methods, to treat the diffusion transport and hydrodynamic convections which control a reahstic crystal growth process from the melt on an industrial scale. 

In a few instances, quantum mechanical calculations on the stability and reactivity of adsorbates have been combined with Monte Carlo simulations of dynamic or kinetic processes. In one example, both the ordering of NO on Rh(lll) during adsorption and its TPD under UHV conditions were reproduced using a dynamic Monte Carlo model involving lateral interactions derived from DFT calculations and different adsorption... 

The experiments and the simulation of CSTR models have revealed a complex dynamic behavior that can be predicted by the classical Andronov-Poincare-Hopf theory, including limit cycles, multiple limit cycles, quasi-periodic oscillations, transitions to chaotic dynamic and chaotic behavior. Examples of self-oscillation for reacting systems can be found in [4], [17], [18], [22], [23], [29], [30], [32], [33], [36]. The paper of Mankin and Hudson [17] where a CSTR with a simple reaction A B takes place, shows that it is possible to drive the reactor to chaos by perturbing the cooling temperature. In the paper by Perez, Font and Montava [22], it has been shown that a CSTR can be driven to chaos by perturbing the coolant flow rate. It has been also deduced, by means of numerical simulation, that periodic, quasi-periodic and chaotic behaviors can appear. 

Only deterministic models for cellular rhythms have been discussed so far. Do such models remain valid when the numbers of molecules involved are small, as may occur in cellular conditions Barkai and Leibler [127] stressed that in the presence of small amounts of mRNA or protein molecules, the effect of molecular noise on circadian rhythms may become significant and may compromise the emergence of coherent periodic oscillations. The way to assess the influence of molecular noise on circadian rhythms is to resort to stochastic simulations [127-129]. Stochastic simulations of the models schematized in Fig. 3A,B show that the dynamic behavior predicted by the corresponding deterministic equations remains valid as long as the maximum numbers of mRNA and protein molecules involved in the circadian clock mechanism are of the order of a few tens and hundreds, respectively [128]. In the presence of molecular noise, the trajectory in the phase space transforms into a cloud of points surrounding the deterministic limit cycle. 

Lassiter, R. R., and D. W. Hayne. A finite difference model for simulation of dynamic processes in ecosystems, pp. 367-440. In B. C. Patton, Ed. Systems Analysis and Simulation In Ecology. Vol. 1. New York Academic Press. 1971. 

Dynetica California Institute of Technology Simulator of dynamic networks written in Java that does model building for systems expressed as reaction networks (http // www.duke.edu/ you/Dynetica page.htm)... 

A problem with this interpretation relates to electrostriction, a process in which the density of the solvent changes about a solute. Shim et al. [243] noted evidence of electrostriction in molecular dynamics simulations of a model chromophore in an IL, and the degree of electrostriction was sensitive to the charge distribution of the solute. This observation does not necessarily contradict the framework above, as some local disruption of solvent structure due to dispersive interactions is inevitable. However, it is desirable to obtain a clearer understanding of the competition between these local interactions and the need to maintain a uniform charge distribution in the liquid. 

Modeling and Simulation of Dynamic and Steady-State Characteristics of Shallow Fluidized Bed Combustors... 

Wistrom, A. and J. Earrell (1998). Simulation and system identification of dynamic models for flocculation control. Water Science Technol. Proc. 7th Int. Workshop on Instrumentation, Control and Automation of Water and Wastewater Treatment and Transport Syst., July 6-9, 1997, Brighton, England, 37, 12, 181-192. Elsevier Science Ltd., Exeter, England. 

The regions la Ila where both types of oscillations are possible were of main interest. The investigation of the character of periodic oscillations was done only by dynamic simulation of the model equations. Three modes of osci-... 

Our main concern here is the dynamics of these molecule numbers A, of the species i in relationship with the condition of the recursive growth of the (proto)cell. In our model there are four basic parameters the total number of molecules N, the total number of molecular species k, the mutation rate p, and the reaction path rate p. By carrying out simulations of this model (choosing a variety of parameter values N, k, p, and p) and also by taking various random networks, we have found that the behaviors are classified into the following three phases [32,33] ... 

P were fixed, the dynamics are reduced entirely to a many-body relaxation problem. None of the theories mentioned in the NY Times article have considered many-body relaxation either at all or directly. Molecular dynamics simulations starting from some interaction potential as well as Monte Carlo simulations of toy models necessarily have captured the effects of many-body relaxation, but these are computer experiments and not theoretical solution of the problem. 

---