Differencing

A very simple procedure for time evolving the wavepacket is the second order differencing method. Here we illustrate how this method is used in conjunction with a fast Fourier transfonn method for evaluating the spatial coordinate derivatives in the Hamiltonian. 

S.A. Silling, Stability and Accuracy of Differencing Schemes for Viscoplastic Models in Wavecodes, SAND91-0141, Sandia National Laboratories, Albuquerque, NM, 1991. 

Emery, A.F., An Evaluation of Several Differencing Methods for Inviscid Fluid-Flow Problems, Sandia Corporation, Livermore Laboratory Monograph No. SCL-DC-66-78, Livermore, CA, 57 pp., March 1967. 

Errors in advection may completely overshadow diffusion. The amplification of random errors with each succeeding step causes numerical instability (or distortion). Higher-order differencing techniques are used to avoid this instability, but they may result in sharp gradients, which may cause negative concentrations to appear in the computations. Many of the numerical instability (distortion) problems can be overcome with a second-moment scheme (9) which advects the moments of the distributions instead of the pollutants alone. Six numerical techniques were investigated (10), including the second-moment scheme three were found that limited numerical distortion the second-moment, the cubic spline, and the chapeau function. 

In order to quantitatively trace behavior as a function of A, it is clear that we need to look at statistical measures that distinguish between ordered and random behavior. To this end, consider the spreading rates of differenc e patterns and entropy. 

Example 8.12 Use the backward differencing method to solve the heat transfer problem of Example 8.3. Select A-t = 0.25 and A = 0.0625. 

There are three adjustable parameters in a Box-Jenkins analysis, one each for autoregression, differencing, and moving average terms. Corrections for cyclical behavior may be added as three optional terms. The approach is flexible, and provides much information. 

The ARIMA analysis evaluates the autocorrelation functions to determine the order of the appropriate moving average and the need for differencing. An appropriate model is chosen and the fit to the data is constructed followed by a careful analysis of the residuals. The parameters are adjusted and the fit is checked again. The process is applied iteratively until the errors are minimized or the model fails to converge. 

Examine the autocorrelation function. The high autocorrelations will indicate the order of the autoregressive part if any. The rate of decay of the autocorrelations will indicate a need for differencing. 

Starting with Eq. (10), the CMVC equation, we assuaied that first-order differencing (d-1) applies and that terms higher than quadratic In B can be safely Ignored. These approximations, generally applicable for most chemical processes, give the equation ... 

The features due to adsorbed water and carbonates observed on the boehmite and y-alumina deserve further attention as they differ from results published by previous investigators. Figure 4 shows a series of difference spectra for adsorption on y-alumina. Spectra were taken after drying the y-alumina at 350 C, cooling to room temperature and carrying out room temperature adsorption. The spectra are the difference of the sample before and after adsorption. Spectrum 4e is the spectrum for the as received alumina differenced with the dried alumina. The positive band at 3400 cm" is due to adsorbed water, and the small negative feature at 3740 cm" is due to isolated hydroxyls on the dried surface. Besides the three... 

In order to increase the accuracy of the approximation to the convective term, not only the nearest-neighbor nodes, but also more distant nodes can be included in the sum appearing in Eq. (37). An example of such a higher order differencing scheme is the QUICK scheme, which was introduced by Leonard [82]. Within the QUICK scheme, an interpolation parabola is fitted through two downstream and one upstream nodes in order to determine O on the control volume face. The un-... 

The QUICK scheme has a truncation error of order h. However, similarly as in the case of the central differencing scheme, at high flow velocities some of the coupling coefficients of Eq. (37) become negative. 

In order to minimize numerical diffusion, Boris and Book [131] formulated the idea of blending a low-order stable differencing scheme with a higher order, potentially unstable, scheme in such a way that steep concentration gradients are maintained as well as possible. The algorithm they proposed consists of the following steps ... 

In order indicate the contribution of each reaction to the overall model more clearly, the matrix of reaction extents is differenced, giving the reaction extents in each interval ... 

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

---