Numerical computational methods considerations

Methods for numerical analyses such as tlris can be obtained from commercial software, and the advent of the computer has considerably eased the work required to obtain numerical values for heat distribution and profiles in a short time, or even continuously if a monitor supplies the boundary values of heat content or temperature during an operation. 

Comparison of Figure 4.1.10a and b demonsfrafes fhaf despite the quantitative differences in fhe deduced values, bofh fhe extraction methods yield a similar trend in the range of equivalence ratios investigated. The overall activation energy is observed to peak close to the stoichiometric condition and decrease on both the lean and rich sides. In addition, the overall activation energy values for n-heptane/air mixtures are observed to be lower when compared with iso-octane/air mixtures for all equivalence ratios under consideration. This similarity of trend and the differences in absolute values using two different extraction methods are also observed in the numerical computations with the available detailed... 

Thus even approximate analytical solutions are often more instructive than the more accurate numerical solutions. However considerable caution must be used in this approach, since some of the approximations, employed to make the equations tractable, can lead to erroneous answers. A number of approximate solution for the hot spot system (Eq 1) are reviewed by Merzhanov and their shortcomings are pointed out (Ref 14). More recently, Friedman (Ref 15) has developed approximate analytical solutions for a planar (semi-infinite slab) hot spot. These were discussed in Sec 4 of Heat Effects on p H39-R of this Vol. To compare Friedman s approximate solutions with the exact numerical solution of Merzhanov we computed r, the hot spot halfwidth, of a planar hot spot by both methods using the same thermal kinetic parameters in both calculations. Over a wide range of input variables, the numerical solution gives values of r which are 33 to 43% greater than the r s of the approximate solution. Thus it appears that the approximate solution, from which the effect of the process variables are much easier to discern than from the numerical solution, gives answers that differ from the exact numerical solution by a nearly constant factor... 

Kolb, M. and Thiel, W., Beyond the MNDO model methodical considerations and numerical results, J. Computational Chem., 14, 775-789, 1993. 

The principal considerations in% choosing a finite-difference method for (7) are accuracy, stability, computation time, and computer storage requirements. Accuracy of a method refers to the degree to which the numerically computed temporal and spatial derivatives approximate the true derivatives. Stability considerations place restrictions on the maxi-... 

The last three years have seen considerable interest in the development of semiclassical methods for treating complex molecular collisions, i.e. those which involve inelastic or reactive processes. One of the reasons for this activity is that the recent work, primarily that of Miller and that of Marcus,2 has shown how numerically computed classical trajectories can be used as input to the semiclassical theory, so that it is not necessary to make any dynamical approximations when applying these semiclassical approaches... 

Numerical methods used to solve a system of ODEs are widely available in computational libraries and through texts such as Numerical Recipes. Certain considerations arise in the use of these standard techniques for nonlinear systems, particularly in models of chemical systems, which often consist of systems of stiff equations that require special care. Stiff equations are characterized by the presence of widely differing time scales, which leads to eigenvalues of the Jacobian matrix differing by many orders of magnitude. 

As discussed in more detail in Sect. 1.1.5, this volume of the Encyclopedia is divided into three broad sections. The first section, of which this chapter is an element, is concerned with introducing some of the basic concepts of electroanalytical chemistry, instrumentation - particularly electronic circuits for control and measurements with electrochemical cells - and an overview of numerical methods. Computational techniques are of considerable importance in treating electrochemical systems quantitatively, so that experimental data can be analyzed appropriately under realistic conditions [8]. Although analytical solutions are available for many common electrochemical techniques and processes, extensions to more complex chemical systems and experimental configurations requires the availability of computational methods to treat coupled reaction-mass transport problems. 

Several algorithms do exist in the literature for numerical computation of fluid flow problems on the basis of primitive variables, in a finite volume framework. One of the most commonly used algorithms of this kind is the SIMPLE (semi-implicit method for pressure-linked equations) algorithm [2]. With reference to a generic staggered control volume for solution of the momentum equation for u (see Fig. 3) and with similar considerations for Ihe other velocity components, major steps of the SIMPLE algorithm can be summarized as follows ... 

The atomic case is included to illustrate how the numerical difficulties escalate as the number of nuclei increases. An atom has just one nucleus. The one-electron problem is separable and exactly soluble in the familiar spherical polar coordinates r, 9, general form R r)Y 9,spherical harmonic). This reduces equation (3) from a 3D PDE to a ID ordinary differential equation (ODE) for the function R r). This is solved numerically. Atomic HF calculations were being performed by Hartree himself and his co-workers in the 1930s, well before the advent of digital computers Further consideration of atomic HF methods is beyond the scope of this article, but a thorough discussion of one FDA implementation is available. ... 

The finite element method (FEM) has become the dominant computational method in structural engineering. In general, the input parameters in the standard FEM assume deterministic values. In earthquake engineering, at least the excitation is often random. However, considerable uncertainties might be involved not only in the excitation of a structure but also in its material and geometric properties. A rational treatment of these uncertainties needs a mathematical concept similar to that underlying the standard FEM. Thus, FEM as a numerical method for solving boundary value problems has to be extended to stochastic boundary value problems. The extension of the FEM to stochastic boundary value problems is called stochastic finite element method (SEEM). 

Figure 9 In the naive multipole method, (a) the molecule is divided into boxes of equal sizes (for simplicity, only a few selected boxes are shown). For calculating interactions between remote boxes, e.g., A and Q—C4, one can use much larger boxes (b) and still achieve high numerical accuracy. This reduces the computational complexity considerably and is an important stepstone on the way to linear scaling.

---