General Systems of Differential Equations

Schneider et al. [66] differentiated Equation (15.8) with respect to time and obtained a system of differential equations, called the rate equations, enabling the creation of auxiliary functions 0i(t) interrelated in the following way  

The functions and describe the total volume and surface of all the spherulites,per unit of volume, neglecting, however, spherulite impingement and truncation while taking into account phantom spherulites. With the same assumptions, 02/8 rand 8 r represent the sum of radii and the number of spherulites in unit volume. 

More recently, Haudin and Chenot [67] also differentiated Equation (15.8), with additional considerations on nucleation based on Avrami s work [58, 59] active nuclei originate from potential nuclei that are activated at the frequency (t). These potential nuclei may disappear either by activation or by absorption by a growing entity. Conversely, new potential nuclei may be generated during cooling. In the 3D case, the authors arrive to a nonlinear system of seven differential equations with seven unknown functions  

The main variables are N, a, Na, Na and three auxiliary functions T, P and Q are added to get a first-order ordinary differential system. N and Na are the number of potential and activated nuclei per unit volume at time respectively Na is the extended number of activated nuclei per unit volume. The model predicts crystallization using three physical parameters the initial density of potential nuclei the frequency of activation q of these nuclei, and the growth rate G. 

A major interest of the differential forms is that they are more suitable for numerical simulation. Both Schneider et al. [66] and Haudin and Chenot [67] approaches have been extended to integrate flow effects [68, 69]. Flow introduces additional nuclei. Two types of morphologies can be considered spherulites (point-like nuclei) or shish-kebabs (thread-like nuclei).The number density Nf of flow-induced nuclei is given by the following type of equation  

---

The use of the new approximation for the Uehling potential permits to decrease the computation errors for this term down to 0.5- 1%. Besides, using such a simple analytical expression for approximating the Uehling potential allows its easy inclusion into the general system of differential equations. This system includes also the Dirac equation and the equations for the matrix elements. 

General first-order kinetics also play an important role for the so-called local eigenvalue analysis of more complicated reaction mechanisms, which are usually described by nonlinear systems of differential equations. Linearization leads to effective general first-order kinetics whose analysis reveals infomiation on the time scales of chemical reactions, species in steady states (quasi-stationarity), or partial equilibria (quasi-equilibrium) [M, and ]. 

The error in Runge-Kutta calculations depends on h, the step size. In systems of differential equations that are said to be stiff, the value of h must be quite small to attain acceptable accuracy. This slows the calculation intolerably. Stiffness in a set of differential equations arises in general when the time constants vary widely in magnitude for different steps. The complications of stiffness for problems in chemical kinetics were first recognized by Curtiss and Hirschfelder.27... 

Crosslinking of many polymers occurs through a complex combination of consecutive and parallel reactions. For those cases in which the chemistry is well understood it is possible to define the general reaction scheme and thus derive the appropriate differential equations describing the cure kinetics. Analytical solutions have been found for some of these systems of differential equations permitting accurate experimental determination of the individual rate constants. 

Solutions are presented in the form of equations, tables, and graphs—most often the last. Serious numerical results generally have to be obtained with computers or powerful calculators. The introductory chapter describes the numerical procedures that are required. Inexpensive software has been used here for integration, differentiation, nonlinear equations, simultaneous equations, systems of differential equations, data regression, curve fitting, and graphing. 

The equivalent to the law of mass action, as encountered in the previous chapter (e.g. in equation (3.22)), are systems of differential equations, defined by the chemical model or the reaction mechanism and the corresponding rate constants. We start with a general chemical reaction, just to practise the notation — it is not a realistic example ... 

We have discussed stresses and strain rates. A critical objective is to relate the two, leading to equations of motion governing how fluid packets are accelerated by the forces acting on them. Generally, we are working toward a differential-equation description of a momentum balance, F = ma. The approach is to represent both the forces and the accelerations as functions of the velocity field. The result will be a system of differential equations in which velocities are the dependent variables and the spatial coordinates and time are the independent variables (i.e., the Navier-Stokes equations). 

Industrial problems are usually more complicated than the earlier problems in this book. But their solutions generally require the same steps, tools and procedures. Therefore, an engineer needs to learn how to handle these problems in both a direct and an integrated way. A typical industrial problem might involve solving one system of differential equations and then solving an algebraic equation (or another DE) at each point of the solution profile. 

The T mapping is not given analytically from the beginning. It is determined with the help of the solution of a system of differential equations. In this sense we do not have it, but we can specify and research general properties of T. 

The model-based controller-observer scheme requires to solve online the system of differential equations of the observer. The phenol-formaldehyde reaction model is characterized by 15 differential equations, and it is, thus, unsuitable for online computations. To overcome this problem, one of the reduced models developed in Sect. 3.8.1 can be adopted. In order to be consistent with the general form of nonchain reactions (2.27) adopted to develop the controller-observer scheme, the reduced model (3.57) with first-order kinetics has been used to design the observer. The mass balances of the reduced model are given by... 

Scheme (3.28) may be taken for the basic conjugation mechanism of chemical reactions, because it generally reproduces die specificity of chemical induction. Before we determine the determinant for scheme (3.28), a system of differential equations taking into account consumption of the components should be composed ... 

The decomposition (2.42) of one exponential into a product of exponentials is a generalization of the formula due to Weyl and Glauber.187 One finds the expressions for /, g, h, ip by differentiating the two members of (2.42) with respect to r and using (A.26) to eliminate the exponentials. One finds a system of differential equations in the variable r, whose solution satisfying the required limit conditions is (2.43). 

A general numerical algorithm of the boundary layer type for stiff systems of differential equations has been proposed by Miranker [173] and applied to a few kinetic problems by Aiken and Lapidus [174,175]. The principle of the method will be briefly described in the case of the following system of differential equations, involving stiff variable x and non-stiff variable y. 

Metabolic control analysis and related theories are based on examining how a system responds to infinitesimally small perturbations. Thus in general one can make use of a linearized set of governing kinetic equations. To be specific, an A-dimensional system governed by a non-linear system of differential equations... 

The system of differential equations is too complex to be solved analytically. Assumptions of a linear adsorption isotherm can be used to obtain analytical solutions, but this approach is generally not applicable to describe affinity chromatography experiments. Several numerical techniques arc used to solve the system of partial differential equations. The other method is to use an analytical solution with simplifying approaches [32] that describe the adsorption process with a single step and a lumped mass transfer coefficient [27],... 

Briefly, the aim of Lie transformations in Hamiltonian theory is to generate a symplectic (that is, canonical) change of variables depending on a small parameter as the general solution of a Hamiltonian system of differential equations. The method was first proposed by Deprit [75] (we follow the presentation in Ref. 76) and can be stated as follows. 

A more sophisticated mathematical approach is required for the results of this chapter because the equations considered cannot be reduced to planar systems. The theorem described in this section and used throughout this chapter is one of great power. When applied to a particular system of differential equations, it generally provides an elegant, simple. 

In order to visualize the molecular selection process in the more general context of optimization of replication rates, we consider the simple case of replication with ultimate accuracy first. In this case we have = di, the value matrix W is diagonal (= An — and the corresponding system of differential equations is weakly coupled by the (t) term only ... 

To "solve" this system of simultaneous equations, we want to be able to calculate the value of [A], [B] and [C] for any value of t. For all but the simplest of these systems of equations, obtaining an exact or analytical expression is difficult or sometimes impossible. Such problems can always be solved by numerical methods, however. Numerical methods are completely general. They can be applied to systems of differential equations of any complexity, and they can be applied to any set of initial conditions. Numerical methods require extensive calculations but this is easily accomplished by spreadsheet methods. 

The influence of activity changes on the dynamic behavior of nonisothermal pseudohomogeneoiis CSTR and axial dispersion tubular reactor (ADTR) with first order catalytic reaction and reversible deactivation due to adsorption and desorption of a poison or inert compound is considered. The mathematical models of these systems are described by systems of differential equations with a small time parameter. Thereforej the singular perturbation methods is used to study several features of their behavior. Its limitations are discussed and other, more general methods are developed. 

Examination and discrimination of hypotheses Each of the proposed hypotheses can pose of detecting those features that mig t be expected for the dif ferent classes of mechanisms. Kinetic equations are derived at this stage for the stationary reaction route rates and (in case of linear mechanisms) equations are obtained for the rate of substance formation or depletion or, in the most general case, systems of differential equations are written. Examination of the kinetic models allows us to devise a plan for model (hypothesis) discrimination using chemical, physicochemical, and kinetic methods. It is worthwhile discussing in more detail these methods f< r discriminating reaction mechanism hypotheses. 

Because of the relative simplicity of this system, we will consider its simulation by CSMP in detail, to acquaint the reader with the general features of CSMP programming. The following system of differential equations describes the above reaction sequence. 

The resulting numerical scheme for solving this system of differential equations can be written in the following general form ... 

---