Mathematical methods differential equations

An understanding of chemical oscillations and wave patterns in the Belousov-Zhabotinskii reaction requires some familiarity with the language and methods of chemical kinetics on one hand and some facility with the mathematics of differential equations on the other. Since not every reader can be expected to know both fields to the extent which we will need later,... 

B. GarciarArchilla, J. M. Sanz-Serna, and R. Skeel. Long-time-step methods for oscillatory differential equations. Applied Mathematics and Computation Reports 1996/7, Universidad de Valladolid, 1996. 

The weighted residual method provides a flexible mathematical framework for the construction of a variety of numerical solution schemes for the differential equations arising in engineering problems. In particular, as is shown in the followmg section, its application in conjunction with the finite element discretizations yields powerful solution algorithms for field problems. To outline this technique we consider a steady-state boundary value problem represented by the following mathematical model... 

In applying quantum mechanics to real chemical problems, one is usually faced with a Schrodinger differential equation for which, to date, no one has found an analytical solution. This is equally true for electronic and nuclear-motion problems. It has therefore proven essential to develop and efficiently implement mathematical methods which can provide approximate solutions to such eigenvalue equations. Two methods are widely used in this context- the variational method and perturbation theory. These tools, whose use permeates virtually all areas of theoretical chemistry, are briefly outlined here, and the details of perturbation theory are amplified in Appendix D. 

Computer simulation of the reactor kinetic hydrodynamic and transport characteristics reduces dependence on phenomenological representations and idealized models and provides visual representations of reactor performance. Modem quantitative representations of laminar and turbulent flows are combined with finite difference algorithms and other advanced mathematical methods to solve coupled nonlinear differential equations. The speed and reduced cost of computation, and the increased cost of laboratory experimentation, make the former increasingly usehil. 

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. 

Traditional control systems are in general based on mathematical models that describe the control system using one or more differential equations that define the system response to its inputs. In many cases, the mathematical model of the control process may not exist or may be too expensive in terms of computer processing power and memory. In these cases a system based on empirical rules may be more effective. In many cases, fuzzy control can be used to improve existing controller systems by adding an extra layer of intelligence to the current control method. 

Calculational problems with the Runge-Kutta technique also surface if the reaction scheme consists of a large number of steps. The number of terms in the rate expression then grows enormously, and for such systems an exact solution appears to be mathematically impossible. One approach is to obtain a solution by an approximation such as the steady-state method. If the investigator can establish that such simplifications are valid, then the problem has been made tractable because the concentrations of certain intermediates can be expressed as the solution of algebraic equations, rather than differential equations. On the other hand, the fact that an approximate solution is simple does not mean that it is correct.28,29... 

Usually the finite difference method or the grid method is aimed at numerical solution of various problems in mathematical physics. Under such an approach the solution of partial differential equations amounts to solving systems of algebraic equations. 

Also, we consider the total approximation method as a constructive method for creating economical difference schemes for the multidimensional equations of mathematical physics. The notion of additive scheme is introduced as a system of operator difference equations that approximates the original differential equation in the total sense. Two quite general heuristic methods (proposed earlier by the author) for obtaining additive economical schemes are discussed in full details. The additive schemes require a new technique for investigating convergence and a new type of a priori estimates that take into account the definition of the property of approximation. 

The radial functions Sni p) and R i(r) may be expressed in terms of the associated Laguerre polynomials L p), whose definition and mathematical properties are discussed in Appendix F. One method for establishing the relationship between Sniip) and L p) is to relate Sni p) in equation (6.50) to the polynomial L p) in equation (F.15). That process, however, is long and tedious. Instead, we show that both quantities are solutions of the same differential equation. 

This closure property is also inherent to a set of differential equations for arbitrary sequences Uk in macromolecules of linear copolymers as well as for analogous fragments in branched polymers. Hence, in principle, the kinetic method enables the determination of statistical characteristics of the chemical structure of noncyclic polymers, provided the Flory principle holds for all the chemical reactions involved in their synthesis. It is essential here that the Flory principle is meant not in its original version but in the extended one [2]. Hence under mathematical modeling the employment of the kinetic models of macro-molecular reactions where the violation of ideality is connected only with the short-range effects will not create new fundamental problems as compared with ideal models. 

The limitation of transfer function representation becomes plain obvious as we tackle more complex problems. For complex systems with multiple inputs and outputs, transfer function matrices can become very clumsy. In the so-called modem control, the method of choice is state space or state variables in time domain—essentially a matrix representation of the model equations. The formulation allows us to make use of theories in linear algebra and differential equations. It is always a mistake to tackle modem control without a firm background in these mathematical topics. For this reason, we will not overreach by doing both the mathematical background and the control together. Without a formal mathematical framework, we will put the explanation in examples as much as possible. The actual state space control has to be delayed until after tackling classical transfer function feedback systems. 

That chemistry and physics are brought together by mathematics is the raison d etre" of tbe present volume. The first three chapters are essentially a review of elementary calculus. After that there are three chapters devoted to differential equations and vector analysis. The remainder of die book is at a somewhat higher level. It is a presentation of group theory and some applications, approximation methods in quantum chemistry, integral transforms and numerical methods. 

References Brown, J. W., and R. V. Churchill, Fourier Series and Boundary Value Problems, 6th ed., McGraw-Hill, New York (2000) Churchill, R. V, Operational Mathematics, 3d ed., McGraw-Hill, New York (1972) Davies, B., Integral Transforms and Their Applications, 3d ed., Springer (2002) Duffy, D. G., Transform Methods for Solving Partial Differential Equations, Chapman Hall/CRC, New York (2004) Varma, A., and M. Morbidelli, Mathematical Methods in Chemical Engineering, Oxford, New York (1997). 

In the following section, we only consider the integration of the equation of linear motion Eq. (20) the procedure for the equation of rotational motion, Eq. (21), will be completely analogous. Mathematically, Eq. (20) represents an initial-value ordinary differential equation. The evolution of particle positions and velocities can be traced by using any kind of method for ordinary differential equations. The simplest method is the first-order integrating scheme, which calculates the values at a time t + 5t from the initial values at time t (which are indicated by the superscript 0 ) via ... 

Ideal reactors can be classified in various ways, but for our purposes the most convenient method uses the mathematical description of the reactor, as listed in Table 14.1. Each of the reactor types in Table 14.1 can be expressed in terms of integral equations, differential equations, or difference equations. Not all real reactors can fit neatly into the classification in Table 14.1, however. The accuracy and precision of the mathematical description rest not only on the character of the mixing and the heat and mass transfer coefficients in the reactor, but also on the validity and analysis of the experimental data used to model the chemical reactions involved. 

As previously mentioned, the QSSA is a common method for eliminating intermediates from the kinetic models of complex catalytic reactions and corresponding transformation of these models. Mathematically, it is a zero-order approximation of the original (singularly perturbed) system of differential equations, which describes kinetics of the complex reaction. We simply replace... 

Although the mathematical methods in this book include algebra, calculus, differential equations, matrix, statistics, and numerical analyses, students with background in algebra and calculus alone are able to understand most of the contents. In addition, since simple models are presented before more complex models and additional parameters are added gradually, students should not worry about the difficulties in mathematics. 

---