Kutta method

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.

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

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. 

Starting with an initial value of and knowing c t), Eq. (8-4) can be solved for c t + At). Once c t + At) is known, the solution process can be repeated to calciilate c 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. To achieve accurate solutions with an Eiiler approach, one often needs to take small steps in time. At. A number of 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 Bunge Kutta method, which involves the following calculations ... 

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

The simulated concentration-time values are displayed graphically in Fig. 5-1. These values could not have been obtained from analytical expressions. Most schemes that comprised a few steps can be handled by the Runge-Kutta methods, provided the time scales of the various steps are comparable. 

There are four equations in four dependent variables, a, d, T, and T. They can be integrated using the Runge-Kutta method as outlined in Appendix 2. Note that they are integrated in the reverse direction e.g., a = at) — similarly for 2 and in Equations (2.47). 

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

The radial equations was then solved using the Runge-Kutta method (7). 

Eqn. (5.5-101) solved numerically by Runge-Kutta method, least squares 0.0207 0.01... 

Figure 4 shows the output of this program, which consists of concentrations of o-methylol, p-methylol and methylene ether groups at various reaction times. Although many integration routines can be used in CSMP calculations, the variable interval Runge-Kutta method was used in this case since that is the option selected when no other method is specified. 

More complicated numerical methods, such as the Runge-Kutta method, yield more accurate solutions, and for precisely formulated problems requiring accurate solutions these methods are helpful. Examples of such problems are the evolution of planetary orbits or the propagation of seismic waves. But the more accurate numerical methods are much harder to understand and to implement than is the reverse Euler method. In the following chapters, therefore, I shall show the wide range of interesting environmental simulations that are possible with simple numerical methods. 

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. 

---