Finite differences numerical issues

The solution of the boundary value problems for the different electrode mechanisms has been carried out by semi-numerical method and by a variety of finite difference numerical techniques. The particular cases of the modified electrodes represent a very complex issue. They require approaches that vary not only as a function of the nature of the modifier, but also of the specific modification, since the mathematics should accounts for the thickness and structure of it. Hence, the problem to solve depends on the exact nature and peculiar experimental conditions under which the modification has been created. 

The numerical methods for partial differential equations can be classified according to the type of equation (see Partial Differential Equations ) parabolic, elliptic, and hyperbolic. This section uses the finite difference method to illustrate the ideas, and these results can be programmed for simple problems. For more complicated problems, though, it is common to rely on computer packages. Thus, some discussion is given to the issues that arise when using computer packages. 

In this approach, the finite volume methods discussed in the previous chapter can be applied to simulate the continuous fluid (in a Eulerian framework). Various algorithms for treating pressure-velocity coupling, and the discussion on other numerical issues like discretization schemes are applicable. The usual interpolation practices (discussed in the previous chapter) can be used. When solving equations of motion for a continuous fluid in the presence of the dispersed phase, the major differences will be (1) consideration of phase volume fraction in calculation of convective and diffusive terms, and (2) calculation of additional source terms due to the presence of dispersed phase particles. For the calculation of phase volume fraction and additional source terms due to dispersed phase particles, it is necessary to calculate trajectories of the dispersed phase particles, in addition to solving the equations of motion of the continuous phase. 

There is an abimdant literature on the comparison between experimental and calculated band profiles for binary mixtures. The most popular methods used have been the forward-backward finite difference scheme and the OCFE method. The former lends itself readily to numerical calculations in many cases representative of the present preoccupations in preparative chromatography. We present first a comparison between the band profiles obtained with the ideal and the equilibrium-dispersive model to illustrate the dispersive influence of the column efficiency. Related to the comparison between these two models is the issue of the use of the hodograph transform of experimental results discussed in Section 11.2.2. Computer experiments are easy to carry out and most instructive because it is possible to show e effects of the change of a single parameter at a time. Some... 

The alternative interpretation is that modeling and numerical issues should deliberately be combined. The governing equations are then solved in physical space often using a second order accurate finite difference or finite volume algorithm. 

Some other versions of the DFT method like the Beijing Density Functional method (BDF) (see the chapter of C. van Wuellen in this issue) were also used for small compounds of the heaviest elements like 111 and 114 [115-117]. There, four-component numerical atomic spinors obtained by finite-difference atomic calculations are used for cores, while basis sets for valence spinors are a combination of numerical atomic spinors and kinetically balanced Slater-type functions. The non-relativistic GGA for F is used there. 

Finally, computational aspects and numerical issues are considered. This includes a short description of the finite difference, the finite volume, and the finite element discretization methods. The chapter ends with some general comments. 

Hydro-mechanical coupling in fractured rock mass is an important issue for many rock mechanics and hydrogeology apphcations (Rutqvist Stephansson 2003). Various numerical methods, i.e. distinct element method (DEM), finite element method (FEM), finite difference method (FDM), etc, are widely used to simulate and analyze the rock hydro-mechanical coupling behaviors. UDEC is DEM software and can be used for modeling the hydro-mechanical coupling behavior of fractured rock masses. 

Unlike FEP, there is no finite difference term in the quantity whose ensemble average we must determine, so differences between the potential surfaces of neighboring X states are not explicitly an issue. Instead, for TI the main concern with respect to the X pathway is that we select enough X points so that the numerical integration over these points is reasonably accurate. For smoothly and slowly varying AG versus X curves, a modest number of points will usually suffice. Fortunately, such a curve is characteristic of most free energy calculations. 

---