Approximation techniques finite differences

The concentration gradient may have to be approximated in finite difference terms (finite differencing techniques are described in more detail in Secs. 4.2 to 4.4). Calculating the mass diffusion rate requires a knowledge of the area, through which the diffusive transfer occurs, since... 

It was soon realised that at least unequal intervals, crowded closely around the UMDE edge, might help with accuracy, and Heinze was the first to use these in 1986 [300], as well as Bard and coworkers [71] in the same year. Taylor followed in 1990 [545]. Real Crank-Nicolson was used in 1996 [138], in a brute force manner, meaning that the linear system was simply solved by LU decomposition, ignoring the sparse nature of the system. More on this below. The ultimate unequal intervals technique is adaptive FEM, and this too has been tried, beginning with Nann [407] and Nann and Heinze [408,409], and followed more recently by a series of papers by Harriman et al. [287,288,289, 290,291,292,293], some of which studies concern microband electrodes and recessed UMDEs. One might think that FEM would make possible the use of very few sample points in the simulation space however, as an example, Harriman et al. [292] used up to about 2000 nodes in their work. This is similar to the number of points one needs to use with conformal mapping and multi-point approximations in finite difference methods, for similar accuracy. 

These coupled second-order partial differential equations do not have a closed-form solution. Accordingly, the approximate numerical technique of finite differences is employed. First, however, the boundary conditions must be prescribed in order to complete the formulation of the problem. Symmetry of the laminate about several planes permits reduction of the region of consideration to a quarter of the laminate cross section in the y-z plane at any value of x as shown in Figure 4-52. There, along the stress-free upper surface. 

This then provides a physical derivation of the finite-difference technique and shows how the solution to the differential equations can be propagated forward in time from a knowledge of the concentration profile at a series of mesh points. Algebraic derivations of the finite-difference equations can be found in most textbooks on numerical analysis. There are a variety of finite-difference approximations ranging from the fully explicit method (illustrated above) via Crank-Nicolson and other weighted implicit forward. schemes to the fully implicit backward method, which can be u.sed to solve the equations. The methods tend to increase in stability and accuracy in the order given. The difference scheme for the cylindrical geometry appropriate for a root is... 

The simulation of a continuous, evaporative, crystallizer is described. Four methods to solve the nonlinear partial differential equation which describes the population dynamics, are compared with respect to their applicability, accuracy, efficiency and robustness. The method of lines transforms the partial differential equation into a set of ordinary differential equations. The Lax-Wendroff technique uses a finite difference approximation, to estimate both the derivative with respect to time and size. The remaining two are based on the method of characteristics. It can be concluded that the method of characteristics with a fixed time grid, the Lax-Wendroff technique and the transformation method, give satisfactory results in most of the applications. However, each of the methods has its o%m particular draw-back. The relevance of the major problems encountered are dicussed and it is concluded that the best method to be used depends very much on the application. 

Use of Computer Simulation to Solve Differential Equations Pertaining to Diffusion Problems. As shown earlier (Section 4.2.11), differential equations used in the solutions of Fick s second law can often be solved analytically by the use of Laplace transform techniques. However, there are some cases in which the equations can be solved more quickly by using an approximate technique known as the finite-difference method (Feldberg, 1968). 

The variational principle has not been widely used in diffusion kinetic problems. Nevertheless, it is such a powerful technique that it is suitable for discussing the many-body problems which have still to be tackled. Wherever approximate methods are necessary, the variational principle should be considered. The trial function(s) should be chosen with care, based on a good idea of the nature of the trial function from its behaviour in certain asymptotic limits. The only application known to the author of the variation principle to a numerical study of a diffusion kinetic problem on a molecular system is that of Delair et al. [377]. They used the variational principle to generate an implicit finite difference scheme for solving the Debye—Smoluchowski equation. Interesting comments have been made by Brykalski and Krason more in the context of heat diffusion [510]. 

It is interesting to compare the digital waveguide simulation technique to the recursion produced by the finite difference approximation (FDA) applied to the wave equation. Recall from (10.10) that the time update recursion for the ideal string digitized via the... 

To analyze this phenomenon further, 2D numerical simulations of (49) and (50) were performed using a central finite difference approximation of the spatial derivatives and a fourth order Runge-Kutta integration of the resulting ordinary differential equations in time. Details of the simulation technique can be found in [114, 119]. The material parameters of the polymer blend PDMS/PEMS were used and the spatial scale = (K/ b )ll2 and time scale r = 2/D were established from the experimental measurements of the structure factor evolution under a homogeneous temperature quench. 

The value of the parameter n in Eqs. (9.13)—(9.14) is different from unity, and these differential equations are nonlinear and cannot be solved analytically. Therefore, Eqs. (9.13)-(9.16) subject to the conditions of Eqs. (9.17)-(9.21) were solved using numerical analysis techniques. Selim et al. (1976a) used the explicit-implicit finite-difference approximations as the method of solution. This was successfully used by Selim et al. (1975) for steady water flow conditions and by Selim et al. (1976a) for transient... 

The most common methodology when solving transient problems using the finite element method, is to perform the usual Garlerkin weighted residual formulation on the spatial derivatives, body forces and time derivative terms, and then using a finite difference scheme to approximate the time derivative. The development, techniques and limitations that we introduced in Chapter 8 will apply here. The time discretization, explicit and implicit methods, stability, numerical diffusion etc., have all been discussed in detail in that chapter. For a general partial differential equation, we can write... 

A common method for solving partial differential equations (PDEs) is known as the method of lines. Here, finite difference approximations for spatial derivatives are used to convert a PDE model to a large set of ordinary differential equations, which are then solved using any of the ODE integration techniques discussed earlier. 

Wajge et al. (1997) attempted to develop rigorous PDAE model for packed batch distillation with and without chemical reaction and used finite difference and orthogonal collocation techniques to solve such model. The main purpose of the study was to investigate the efficiencies of the numerical methods employed. The authors observed that the collocation techniques are computationally more efficient compared to the finite difference method, however the order of approximating polynomial needs to be carefully chosen to achieve a right balance between accuracy and efficiency. See the original reference for further details. 

There has been little recent work on methods for differentiable functions which avoid explicit evaluation of derivatives. Powell s conjugate direction method 36 is still used, but the generally accepted approach is now to use standard quasi-Newton methods with finite-difference approximations to the derivatives. On the other hand there has been considerable interest in methods for nondifferentiable functions, as shown by the collection of papers edited by Balinski and Wolfe 37, in which the technique described by Lemarechal is of particular interest. Other contributions in this difficult field are due to Shor 38, ... 

In the following sections some background information on stiff ordinary differential equations will be given and the general finite difference approximations for particle temperatures will be derived. Later, the technique will be applied to coal pyrolysis in a transport reactor where the difference equations for reaction kinetics will be discussed and the calculation results will be compared with those obtained by the previously established techniques. 

As an alternative to the simultaneous solution of stiff differential equations through an implicit technique a method is described here which approximates the solution by successive computations of the corresponding finite difference equations. The successive nature of this method essentially decouples the K(N + 1)... 

The Finite Difference Approximations and the Iterative Technique. The reaction kinetics are dependent on the rest of the equations through particle temperatures. If we have an initial guess on the coal particle temperature history during its total residence in the reactor then the relations given by Eq. (15) are decoupled from the rest and can be solved separately. Consequently, the volatilization estimates, v j(t) s, can be used to update particle temperatures with Eqs. (11) through (14). 

This procedure had converged in 4 or 5 iterations to four significant figures for all cases tried in this study. The accuracy of the calculations depends on the time increment At because the finite difference approximations become more accurate as At gets smaller. A summary of some iteration results and a comparison between this technique and the numerical integration with Gear s method will be presented after the following discussion on the stability of the temperature equation. 

---