Method Rung-Kutta

The numerical methods that are available in Stella are Euler s method, Runge-Kutta second order, or Runge-Kutta fourth order. One of the menu items allows the user to specify the length of simulation time, as well as the time increment dt, and the type of numerical method that is to be employed. 

Sometimes the ODEs that arise in studies in nonlinear dynamics can be solved using explicit methods (such as the forward Euler) which require less computations per step and are thus cheaper and ter to implement. The Runge-Kutta femily of algorithms are a popular implementation of the explicit methods. Runge—Kutta methods begin with a Taylor series expansion the order of the particular Runge-Kutta method used is simply the highest order term retained in the Taylor series. 

To accelerate the restart for multistep methods Runge-Kutta methods can be used to generate the necessary starting values, [SB95]. 

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. 

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. 

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. 

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. 

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. 

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. 

Equations 6-94 and 6-97 are first order differential equations, and it is possible to solve for both the eonversion and temperature of hydrogenation of nitrobenzene relative to the reaetor length of 25 em. A eomputer program PLUG61 has been developed employing the Runge-Kutta fourth order method to determine the temperature and eonversion using a eatalyst bed step size of 0.5 em. Table 6-6 shows... 

---