Partial differential equations, numerical computational methods

C. Johnson. Numerical solution of Partial Differential Equations by the Finite Element Method. Cambridge University Press, Cambridge, 1990. An excellent introductory book on the finite element method. The text assumes mathematical background of first year graduate students in applied mathematics and computer science. An excellent introduction to the theory of adaptive methods. 

Bertoluzza S. A wavelet collocation method for the numerical solution of partial differential equations. Appl Comput Harmonic Anal 1996 3 1-9. 

At this point we can summarize the steps required to implement the Gauss-Newton method for PDE models. At each iteration, given the current estimate of the parameters, ky we obtain w(t,z) and G(t,z) by solving numerically the state and sensitivity partial differential equations. Using these values we compute the model output, y(t k(i)), and the output sensitivity matrix, (5yr/5k)T for each data point i=l,...,N. Subsequently, these are used to set up matrix A and vector b. Solution of the linear equation yields Ak(jH) and hence k°M) is obtained. The bisection rule to yield an acceptable step-size at each iteration of the Gauss-Newton method should also be used. 

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. 

The partial differential equations representing material and energy balances of a reaction in a packed bed are rarely solvable by analytical means, except perhaps when the reaction is of zero or first order. Two examples of derivation of the equations and their analytical solutions are P8.0.1.01 and P8.01.02. In more complex cases analytical, approximations can be made (by "Collocation" or "Perturbation", for instance), but these usually are quite sophisticated to apply. Numerical solutions, on the other hand, are simple in concept and are readily implemented on a computer. Two such methods that are suited to nonlinear kinetics problems will be described. 

CONSTANTTNIDEE, Applied Numerical Methods with Personal Computers, McGraw-Hill, 1987. Nonlinear regression, partial differential equations, matrix manipulations, and a more flexible program for simultaneous ODEs. 

While the boundary layer equations that were derived in the preceding sections are much simpler than the general equations from which they were derived, they still form a complex set of simultaneous partial differential equations. Analytical solutions to this set of equations have been obtained in a few important cases. For the majority of flows, however, a numerical solution procedure must be adopted. Such solutions are readily obtained today using modest modem computing facilities. This was, however, not always so. For this reason, approximate solutions to the boundary layer equations have in the past been quite widely used. While such methods of solution are less important today, they are still used to some extent. One such approach will, therefore, be considered in the present text. 

P. Concus, G. H. Golub, and D. P. O Leary, in Sparse Matrix Computations, J. R. Bunch and D. J. Rose, Eds., Academic Press, New York, 1976, pp. 309—332. A Generalized Conjugate Gradient Method for the Numerical Solution of Elliptic Partial Differential Equations. 

The method of moments possibly new numerical methods for the solution of the partial differential equations for g(z) and computer simulation of the kinetic process are the possible ways of solving the problem. 

An orthogonal collocation method for elliptic partial differential equations is presented and used to solve the equations resulting from a two-phase two-dimensional description of a packed bed. Comparisons are made between the computational results and experimental results obtained from earlier work. Some qualitative discrimination between rival correlations for the two-phase model parameters is possible on the basis of these comparisons. The validity of the numerical method is shown by applying it to a one-phase packed-bed model for which an analytical solution is available problems arising from a discontinuity in the wall boundary condition and from the semi-infinite domain of the differential operator are discussed. 

Equations (8) form an infinite system of coupled non-linear partial differential equations for the fransformed potentials,, . For computation purposes, system (8) is also truncated at the Ntii row and colimm, with N sufficiently large for the required convergence. A few automatic numerical integrators for tiiis class of one-dimensional partial differential systems are now readily available, such as those based on tiie Method of Lines [41, 52]. Once the transformed potentials have been computed from numerical solution of system (8), tiie inversion formula Eq.(7.b) is recalled to reconstruct the original potentials, in explicit form along thejc v -iables. 

---