Crank-Nicholson finite-difference implicit method

When faster reactions are dealt with, it may be profitable to remove the At/Ay2 < 0.5 condition and use an implicit method such as the Crank-Nicholson method.15 17 The finite difference approximation is then applied at the value of t corresponding to the middle of the j to j + 1 interval, leading to... 

In the last section we considered explicit expressions which predict the concentrations in elements at (t + At) from information at time t. An error is introduced due to asymmetry in relation to the simulation time. For this reason implicit methods, which predict what will be the next value and use this in the calculation, were developed. The version most used is the Crank-Nicholson method. Orthogonal collocation, which involves the resolution of a set of simultaneous differential equations, has also been employed. Accuracy is better, but computation time is greater, and the necessity of specifying the conditions can be difficult for a complex electrode mechanism. In this case the finite difference method is preferable7. 

Methods applying reverse differences in time are called implicit. Generally these implicit methods, as e.g. the Crank-Nicholson method, show high numerical stability. On the other side, there are explicit methods, and the methods of iterative solution algorithms. Besides the strong attenuation (numeric dispersion) there is another problem with the finite differences method, and that is the oscillation. 

The numerical solution to the advection-dispersion equation and associated adsorption equations can be performed using finite difference schemes, either in their implicit and/or explicit form. In the one-dimensional MRTM model (Selim et al., 1990), the Crank-Nicholson algorithm was applied to solve the governing equations of the chemical transport and retention in soils. The web-based simulation system for the one-dimensional MRTM model is detailed in Zeng et al. (2002). The alternating direction-implicit (ADI) method is used here to solve the three-dimensional models. 

The governing equations are discretized by using the finite difference method. The Reynolds equation solution leads to solving a tri-diagonal system of linear equations. Using the semi-implicit scheme of Crank-Nicholson solves the energy equations in the film and in the rings. 

---