Runge-Kutta method order

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.

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. 

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. 

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... 

This improved procedure is an example of the Runge-Kutta method of numerical integration. Because the derivative was evaluated at two points in the interval, this is called a second-order Runge-Kutta process. We chose to evaluate the mean derivative at points Pq and Pi, but because there is an infinite number of points in the interval, an infinite number of choices for the two points could have been made. In calculating the average for such choices appropriate weights must be assigned. 

More than two points can be used in the Runge-Kutta method, and the fourth-order Runge-Kutta integration is commonly employed. Obviously computers are... 

The Runge-Kutta method takes the weighted average of the slope at the left end point of the interval and at some intermediate point. This method can be extended to a fourth-order procedure with error 0 (Ax) and is given by... 

Pag), where y o mole fraction of A in bulk gas phase can be determined iteratively, yAi = mole fiaction of A in gas inlet. Equations (1) to (6) were solved using fourth order Runge-Kutta method [1, 8]. The value of enhancement factor, E, was predicted using equation of Van Krevelen and Hoftijzer [2]. 

Solve Eq. (3) to obtain B+1 using the second-order TVD-Runge-Kutta method presented as follows ... 

Solve the convection equation of high order (3rd order) essentially non-oscillatory (ENO) upwind scheme (Sussman et al., 1994) is used to calculate the convective term V V

velocity field P". The time advancement is accomplished using the second-order total variation diminishing (TVD) Runge-Kutta method (Chen and Fan, 2004). 

Both POLYMATH and CONSTANTINIDES use this method and also a fourth order Runge-Kutta method. Other methods are available in other software, but these two are adequate for the present book. 

The higher order ODEs are reduced to systems of first-order equations and solved by the Runge-Kutta method. The missing condition at the initial point is estimated until the condition at the other end is satisfied. After two trials, linear interpolation is applied after three or more, Lagrange interpolation is applied. 

Fixed step, 2nd-order, Runge-Kutta method (RK2). 

Figure 3.6 Numerical solution of ordinary differential equations sketch of the four steps of the Runge-Kutta method to the order four giving the n+1 th estimate y(n+1) from the nth estimate y ". ...

The situation is different for those readers who do not have access to Matlab and rely completely on Excel. In the following, we explain how a fourth order Runge-Kutta method can be incorporated into a spreadsheet and used to solve non-stiff ODE s. 

Fourth Order Runge-Kutta Method in Excel... 

---