Partial differential equations numerical analysis

See also Numerical Analysis and Approximate Methods and General References References for General and Specific Topics—Advanced Engineering Mathematics for additional references on topics in ordinary and partial differential equations. 

The partial differential equations describing the catalyst particle are discretized with central finite difference formulae with respect to the spatial coordinate [50]. Typically, around 10-20 discretization points are enough for the particle. The ordinary differential equations (ODEs) created are solved with respect to time together with the ODEs of the bulk phase. Since the system is stiff, the computer code of Hindmarsh [51] is used as the ODE solver. In general, the simulations progressed without numerical problems. The final values of the rate constants, along with their temperature dependencies, can be obtained with nonlinear regression analysis. The differential equations were solved in situ with the backward... 

The modeling of steady-state problems in combustion and heat and mass transfer can often be reduced to the solution of a system of ordinary or partial differential equations. In many of these systems the governing equations are highly nonlinear and one must employ numerical methods to obtain approximate solutions. The solutions of these problems can also depend upon one or more physical/chemical parameters. For example, the parameters may include the strain rate or the equivalence ratio in a counterflow premixed laminar flame (1-2). In some cases the combustion scientist is interested in knowing how the system mil behave if one or more of these parameters is varied. This information can be obtained by applying a first-order sensitivity analysis to the physical system (3). In other cases, the researcher may want to know how the system actually behaves as the parameters are adjusted. As an example, in the counterflow premixed laminar flame problem, a solution could be obtained for a specified value of the strain... 

Stone, H.L., "Iterative Solution of Implicit Approximations of Multi-Dimensional Partial Differential Equations", SIAM J. Numerical Analysis, 5, 530-558 (1968). 

Burns and Curtiss (1972) and Burns et al. (1984) have used the Facsimile program developed at AERE, Harwell to obtain a numerical solution of simultaneous partial differential equations of diffusion kinetics (see Eq. 7.1). In this procedure, the changes in the number of reactant species in concentric shells (spherical or cylindrical) by diffusion and reaction are calculated by a march of steps method. A very similar procedure has been adopted by Pimblott and La Verne (1990 La Verne and Pimblott, 1991). Later, Pimblott et al. (1996) analyzed carefully the relationship between the electron scavenging yield and the time dependence of eh yield through the Laplace transform, an idea first suggested by Balkas et al. (1970). These authors corrected for the artifactual effects of the experiments on eh decay and took into account the more recent data of Chernovitz and Jonah (1988). Their analysis raises the yield of eh at 100 ps to 4.8, in conformity with the value of Sumiyoshi et al. (1985). They also conclude that the time dependence of the eh yield and the yield of electron scavenging conform to each other through Laplace transform, but that neither is predicted correctly by the diffusion-kinetic model of water radiolysis. 

This approach of subdividing space into an increasing number of discrete pieces provides the basis for many numerical computer models (e.g., the so-called finite difference models) an example will be discussed in Chapter 23. Although these models are extremely powerful and convenient for the analysis of field data, they often conceal the basic principles which are responsible for a given result. Therefore, in the next chapter we will discuss models which are not only continuous in time, but also continuous along one or several space axes. In this context continuous in space means that the concentrations are given not only as steadily varying functions in time [QY)], but also as functions in space [C,(r,x) or C,(t,x,y,z)]. Such models lead to partial differential equations. A prominent example is Fick s second law (Eq. 18-14). 

In many cases ordinary differential equations (ODEs) provide adequate models of chemical reactors. When partial differential equations become necessary, their discretization will again lead to large systems of ODEs. Numerical methods for the location, continuation and stability analysis of periodic and quasi-periodic trajectories of systems of coupled nonlinear ODEs (both autonomous and nonautonomous) are extensively used in this work. We are not concerned with the numerical description of deterministic chaotic trajectories where they occur, we have merely inferred them from bifurcation sequences known to lead to deterministic chaos. Extensive literature, as well as a wide choice of algorithms, is available for the numerical analysis of periodic trajectories (Keller, 1976,1977 Curry, 1979 Doedel, 1981 Seydel, 1981 Schwartz, 1983 Kubicek and Hlavacek, 1983 Aluko and Chang, 1984). 

For many applications, interpolations of functions of two or three variables defined in two-and three-dimensional domains must be considered. For example, global interpolations in two- and three-dimensional systems are analogous to polynomial interpolation in onedimensional systems however, global interpolants do not exist in 2- and 3D. This is a big drawback in numerical analysis because a basic tool available for one variable is not available for multivariable approximation [21], The best developed aspect of this theory is that of piecewise polynomial approximation, associated with finite element and finite volume approximations for partial differential equations, which will be examined in detail in Chapters 9 and 10. 

From a mathematical point of view the PCM models can be unified according to the approach they use to solve the linear partial differential equations determining the electrostatic interactions between solute and solvent. This analysis is presented by Cances who reviews both the mathematical and the numerical aspects of such an integral equation approach when applied to PCM models. 

In order to have theoretical relationships with which experimental data can be compared for analysis it is necessary to obtain solutions to the partial differential equations describing the diffusion-kinetic behaviour of the electrode process. Only a very brief account f the theoretical methods is given here and this is done merely to provide a basis for an appreciation of the problems involved and to point out where detailed treatments can be found. A very lucid introduction to the theoretical methods of dealing with transient electrochemical response has appeared (MacDonald, 1977) which is highly recommended in addition to the classic detailed treatment (Delahay, 1954). Analytical solutions of the partial differential equations are possible only in the most simple cases. In more complex cases either numerical methods are used to solve the equations or they are transformed into finite difference forms and solved by digital simulation. 

Unlike the radiant loss from an optically thin flame, conductive or convective losses never can be consistent exactly with the plane-flame assumption that has been employed in our development. Loss analyses must consider non-one-dimensional heat transfer and should also take flame shapes into account if high accuracy is to be achieved. This is difficult to accomplish by methods other than numerical integration of partial differential equations. Therefore, extinction formulas that in principle can be used with an accuracy as great as that of equation (21) for radiant loss are unavailable for convective or conductive loss. The most convenient approach in accounting for convective or conductive losses appears to be to employ equation (24) with L(7 ) estimated from an approximate analysis. The accuracy of the extinction prediction then depends mainly on the accuracy of the heat-loss estimate. Rough heat-loss estimates are readily obtained from overall balances. 

In order to build up and implement efficient numerical schemes for partial differential equations, it is necessary to have informations on the mathematical properties of the system of equations—this has been done in the previous sections—as well as on the stability and the convergence properties of the schemes this is the purpose of numerical analysis. 

Integral Equation Solutions. As a consequence of the quasi-steady approximation for gas-phase transport processes, a rigorous simultaneous solution of the governing differential equations is not necessary. This mathematical simplification permits independent analytical solution of each of the ordinary and partial differential equations for selected boundary conditions. Matching of the remaining boundary condition can be accomplished by an iterative numerical analysis of the solutions to the governing differential equations. 

Finite-Difference Methods. The numerical analysis literature abounds with finite difference methods for the numerical solution of partial differential equations. While these methods have been successfully applied in the solution of two-dimensional problems in fluid mechanics and diffusion (24, 25), there is little reported experience in the solution of three-dimensional, time-dependent, nonlinear problems. Application of these techniques, then, must proceed by extending methods successfully applied in two-dimensional formulations to the more complex problem of solving (7). The various types of finite-difference methods applicable in the solution of partial differential equations and their advantages and disadvantages are discussed by von Rosenberg (26), Forsythe and Wasow (27), and Ames (2S). 

Analytical integration and differentiation, linear algebra, statistics, optimization, numerical integration, Fourier analysis, filtering, ordinary differential equations, partial differential equations, and matrix manipulations... 

Forsythe, G.E. and Rosenbloom, P.C. 1958. Numerical Analysis and Partial Differential Equations. NY John Wiley and Sons. 

In this study a simulation of such processes for various arrangements of capillaries was performed to attain a better understanding of oxygen supply conditions in the brain cortex. Oxygen tension values were calculated for steady-state conditions by numerical solution of partial differential equations for different capillary-tissue systems. The analysis helps to explain the experimental results obtained from the microstructure of rat... 

The analysis of batch crystallizers normally requires the consideration of the time-dependent, batch conservation equations (e.g., population, mass, and energy balances), together with appropriate nucleation and growth kinetic equations. The solution of these nonlinear partial differential equations is relatively difficult. Under certain conditions, these batch conservation equations can be solved numerically by a moment technique. Several simple and useful techniques to study crystallization kinetics and CSDs are discussed. These include the thermal response technique, the desupersaturation curve technique, the cumulative CSD method, and the characterization of CSD maximum. 

Computational techniques are centrally important at every stage of investigation of nonlinear dynamical systems. We have reviewed the main theoretical and computational tools used in studying these problems among these are bifurcation and stability analysis, numerical techniques for the solution of ordinary differential equations and partial differential equations, continuation methods, coupled lattice and cellular automata methods for the simulation of spatiotemporal phenomena, geometric representations of phase space attractors, and the numerical analysis of experimental data through the reconstruction of phase portraits, including the calculation of correlation dimensions and Lyapunov exponents from the data. 

This transfer function can now be studied in the frequency domain. It should be noted that these are linear partial differential equations and that the process of frequency domain analysis is appropriate. The range of values of e = 0.01 to 0.2, M = 5 to 20, and R = 0.75 have been established [Grant and Cotton, 1991] in a numerical finite difference solution of the governing equations. Having established these values the frequency response can be completed. 

---