Runge Kutta

It is therefore concluded that this method, using explicit RK as described in Chap. 4, is not worthwhile mainly because of the A limitation. 

There are, however, implicit variants of RK, and these may have promise. There are several classes of these, see a thorough text on the subject [284,286]. One of these classes, the Rosenbrock method, has been recently examined [100,113, and see references therein] and found very efficient. This is described in its own Sect. 9.4, below. 

---

Fig. 2. The time evolution of the total energy of four water molecules (potential-energy details are given in [48]) as propagated by the symplectic Verlet method (solid) and the nonsymplectic fourth-order Runge-Kutta method (dashed pattern) for Newtonian dynamics at two timestep values.

The IE and IM methods described above turn out to be quite special in that IE s damping is extreme and IM s resonance patterns are quite severe relative to related symplectic methods. However, success was not much greater with a symplectic implicit Runge-Kutta integrator examined by Janezic and coworkers [40],... 

D. Janezic and B. Orel. Implicit Runge-Kutta method for molecular dynamics integration. J. Chem. Info. Comp. Sd., 33 252-257, 1993. 

Janezic, D., Orel, B. Implicit Runge-Kutta Method for Molecular Dynamics Integration. J. Chem. Inf. Comput. Sci. 33 (1993) 252-257 Janezic, D., Orel, B. Improvement of Methods for Molecular Dynamics Integration. Int. J. Quant. Chem. 51 (1994) 407-415... 

Thus we find that the choice of quaternion variables introduces barriers to efficient symplectic-reversible discretization, typically forcing us to use some off-the-shelf explicit numerical integrator for general systems such as a Runge-Kutta or predictor-corrector method. 

The reaction rate equations give differential equations that can be solved with methods such as the Runge-Kutta [14] integration or the Gear algorithm [15]. 

Runge-Kutta methods are explicit methods that use several function evaluations for each time step. Runge-Kutta methods are traditionally written for/(f, y). The first-order Runge-Kutta method is Euler s method. A second-order Runge-Kutta method is... 

A popular fourth-order Runge-Kutta method is the Runge-Kutta-Feldberg formulas (Ref. Ill), which have the property that the method is fourth-order but achieves fifth-order accuracy. The popular integration package RKF45 is based on this method. 

Runge-Kutta, 2nd order, A At < 2 Runge-Kutta-Feldberg, A At < 3.0... 

Solution of Batch-Mill Equations In general, the grinding equation can be solved by numerical methods—for example, the Luler technique (Austin and Gardner, l.st Furopean Symposium on Size Reduction, 1962) or the Runge-Kutta technique. The matrix method is a particiilarly convenient fornmlation of the Euler technique. 

To check the effect of integration, the following algorithms were tried Euler, explicit Runge-Kutta, semi-implicit and implicit Runge-Kutta with stepwise adjustment. All gave essentially identical results. In most cases, equations do not get stiff before the onset of temperature runaway. Above that, results are not interesting since tubular reactors should not be... 

Therefore, the slope of the linear plot Cg versus gives the ratio kj/kj. Knowing kj -i- kj and kj/kj, the values of kj and kj ean be determined as shown in Figure 3-10. Coneentration profiles of eom-ponents A, B, and C in a bateh system using the differential Equations 3-95, 3-96, 3-97 and the Runge-Kutta fourth order numerieal method for the ease when Cgg =Cco = 0 nd kj > kj are reviewed in Chapter 5. 

Coneentration profiles ean he developed witli time using tlie differential Equations 3-180, 3-181, and 3-182, respeetively, with the Runge-Kutta fouith order mediod at known values of kj and kj for a hateh system. 

The following details establish reactor performance, considers the overall fractional yield, and predicts the concentration profiles with time of complex reactions in batch systems using the Runge-Kutta numerical method of analysis. 

Cbo = Cco = 0- For known values of kj and kj, simulate die eoneen-tradons of A, B, and C for 10 minutes at a dme interval of t = 0.5 min. A eomputer program has been developed using die Runge-Kutta fourdi order mediod to determine die eoneentrations of A, B, and C. The differential Equations 5-64, 5-65, and 5-66 are expressed, respeetively, in die form of X-arrays and funetions in die eomputer program as... 

Equations 5-81, 5-82 and 5-83 are first order differential equations that ean be solved simultaneously using the Runge-Kutta fourth order method. Consider two eases ... 

Equations 5-88, 5-89, and 5-90 are first order differential equations and the Runge-Kutta fourth order method with the boundary eonditions is used to determine the eoneentrations versus time of the eomponents. 

The following examples review some complex reactions and determine the concentrations history for a specified period using the Runge-Kutta fourth order numerical method. 

Equations 5-110, 5-112, 5-113, and 5-114 are first order differential equations and the Runge-Kutta fourth order numerical method is used to determine the concentrations of A, B, C, and D, with time, with a time increment h = At = 0.5 min for a period of 10 minutes. The computer program BATCH57 determines the concentration profiles at an interval of 0.5 min for 10 minutes. Table 5-6 gives the results of the computer program and Figure 5-16 shows the concentration profiles of A, B, C, and D from the start of the batch reaction to the final time of 10 minutes. 

Equations 5-118, 5-120, 5-121, and 5-122 are first order differential equations. A simulation exereise on the above equations using the Runge-Kutta fourth order method, ean determine the numher of moles with time inerement h = At = 0.2 hr for 2 hours. Computer program BATCH58 evaluates the numher of moles of eaeh eomponent as a funetion of time. Table 5-7 gives the results of the simulation, and Eigure 5-17 shows the plots of the eoneentrations versus time. 

The differential Equations 5-127 and 5-133 are solved using the Runge-Kutta fourth order while eomponents C and D are ealeulated from the mass balanee of Equations 5-130 and 5-132, respeetively, at a time inerement of h = At = 0.5 min. The eomputer program BATCH59 ealeulates the eoneentrations of A, B, C, and D as funetion... 

Equations 5-146, 5-149, and 5-152 are first order differential equations. The eoneentration profiles of A, B, C, and the volume V of the bateh using Equation 5-137 is simulated with respeet to time using the Runge-Kutta fourth order numerieal method. 

---