Differential equations, special forms

Mathematical and Computational Implementation. Solution of the complex systems of partial differential equations governing both the evolution of pollutant concentrations and meteorological variables, eg, winds, requires specialized mathematical techniques. Comparing the two sets of equations governing pollutant dynamics (eq. 5) and meteorology (eqs. 12—14) shows that in both cases they can be put in the form ... 

Cauchy Momentum and Navier-Stokes Equations The differential equations for conservation of momentum are called the Cauchy momentum equations. These may be found in general form in most fliiid mechanics texts (e.g., Slatteiy [ibid.] Denu Whitaker and Schlichting). For the important special case of an incompressible Newtonian fluid with constant viscosity, substitution of Eqs. (6-22) and (6-24) lead to the Navier-Stokes equations, whose three Cartesian components are... 

These equations form a fourth-order system of differential equations which cannot be solved analytically in almost all interesting nonseparable cases. Further, according to these equations, the particle slides from the hump of the upside-down potential — V(Q) (see fig. 24), and, unless the initial conditions are specially chosen, it exercises an infinite aperiodic motion. In other words, the instanton trajectory with the required periodic boundary conditions,... 

An alternative procedure is the dynamic programming method of Bellman (1957) which is based on the principle of optimality and the imbedding approach. The principle of optimality yields the Hamilton-Jacobi partial differential equation, whose solution results in an optimal control policy. Euler-Lagrange and Pontrya-gin s equations are applicable to systems with non-linear, time-varying state equations and non-quadratic, time varying performance criteria. The Hamilton-Jacobi equation is usually solved for the important and special case of the linear time-invariant plant with quadratic performance criterion (called the performance index), which takes the form of the matrix Riccati (1724) equation. This produces an optimal control law as a linear function of the state vector components which is always stable, providing the system is controllable. 

The boundary conditions were used to obtain special forms of these equations at the boundary nodes. The complete pelletizer model contained a total of 207 differential and algebraic equations which were solved simultaneously. The differential/algebraic program, DASSL, developed at Sandia National Laboratories 2., .) was used. The solution procedure is outlined in Figure 5. 

Let us consider the special class of problems where all state variables are measured and the parameters enter in a linear fashion into the governing differential equations. As usual, we assume that x is the n-dimemional vector of state variables and k is the p-dimemional vector of unknown parameters. The structure of the ODE mode is of the form... 

Several years ago Baer proposed the use of a matrix A, that transforms the adiabatic electronic set to a diabatic one [72], (For a special twofold set this was discussed in [286,287].) Computations performed with the diabatic set are much simpler than those with the adiabatic set. Subject to certain conditions, A is the solution of a set of first order partial differential equations. A is unitary and has the form of a path-ordered phase factor, in which the phase can be formally written as... 

So the principal idea of the method of ansatzes, as well as of the direct method of reduction of partial differential equations is a special choice of the class of functions to which the solution to be found should belong. Within the framework of the preceding methods, a solution of system (46) is sought in the form... 

J is the vector of the transformed inhomogeneous terms. With the special form of A, Eq. 4 decays into n uncoupled differential equations ... 

Exercise. A special but not uncommon form of At(t) is (t)B with constant matrix B and scalar random function (t). If, moreover (2.8) holds it is again possible to write (3.3) without time ordering. Give the resulting differential equation for u(t) to all orders in arc. 

This differential equation is generally nonlinear [depending upon the form of D(c), and solutions therefore can be obtained analytically only in certain special cases which are not discussed here [4]. 

It is surprising how often some of the special forms of easily soluble differential equations turn up. For single equations in which the derivative is a function only of the state, dxldt = f(x), x(0) = X, the solution by quadratures is immediate... 

Equation (11) is written in the form of Newton s second law and states that the mass times acceleration of a fluid particle is equal to the sum of the forces causing that acceleration. In flow problems that are accelerationless (Dx/Dt = 0) it is sometimes possible to solve Eq. (11) for the stress distribution independently of any knowledge of the velocity field in the system. One special case where this useful feature of these equations occurs is the case of rectilinear pipe flow. In this special case the solution of complex fluid flow problems is greatly simplified because the stress distribution can be discovered before the constitutive relation must be introduced. This means that only a first-order differential equation must be solved rather than a second-order (and often nonlinear) one. The following are the components of Eq. (11) in rectangular Cartesian, cylindrical polar, and spherical polar coordinates ... 

In the presence of multiple states, the right-hand-side term consists of sums, products, and nesting of elementary functions such as logy, expy, and trigonometric functions, called the S -system formalism [602]. Using it as a canonical form, special numerical methods were developed to integrate such systems [603]. The simple example of the diffusion-limited or dimensionally restricted homodimeric reaction was presented in Section 2.5.3. Equation 2.23 is the traditional rate law with concentration squared and time-varying time constant k (t), whereas (2.22) is the power law (c7 (t)) in the state differential equation with constant rate. 

The balance equations (7.12 and 7.13) form a set of coupled ordinary differential equations, which has to be solved numerically. Analytical integration is possible for special cases only. 

Hinshelwood model developed in general form by Dagonnier et al. (295) as described in Section IV,A has also been applied to the special case of CO oxidation by the same authors (324). This model is essentially composed of the same differential equations as in the general mechanism. An oscillation... 

Consider here the special case of linear second-order differential equations with constant coefficients, written in general form as... 

In the application of the heat conduction equation in its general form (2.8) a series of simplifying assumptions are made, through which a number of special differential equations, tailor made for certain problems, are obtained. A significant simplification is the assumption of constant material properties A and c. The linear partial differential equations which emerge in this case are discussed in the next section. Further simple cases are... 

Model representations in Laplace transform form are mainly used in control theory. This approach is limited to linear differential equation systems or their linearized approximations and is achieved by a combination of first-order lag function and time delays. This limitation together with additional complications of modelling procedures are the main reasons for not using this method here. Specialized books in control theory as mentioned above use this approach and are available to the interested reader. 

Recall from chap. 2 that often in the solution of differential equations, useful strategies are constructed on the basis of the weak form of the governing equation of interest in which a differential equation is replaced by an integral statement of the same governing principle. In the previous chapter, we described the finite element method, with special reference to the theory of linear elasticity, and we showed how a weak statement of the equilibrium equations could be constructed. In the present section, we wish to exploit such thinking within the context of the Schrodinger equation, with special reference to the problem of the particle in a box considered above and its two-dimensional generalization to the problem of a quantum corral. 

The Mass-Action representation is clearly a special case of the GMA representation in which all exponents are positive integers. The Michaelis-Menten representation is, in turn, a special case of the traditional Mass-Action representation in which two important restrictions have been imposed (Savageau, 1992). First, it is assumed that the mechanism is in quasi-steady state. The derivatives of the dependent state variables in the Mass-Action Formalism can then be set to zero, thereby reducing the description from differential equations to algebraic equations. Second, it is assumed that complexes do not occur between different forms of an enzyme or between different enzymes. The algebraic equations will then be linear in the concentrations of the various enzyme forms, and one can derive the rational function that is the representation of the rate law within the Michaelis-Menten Formalism. 

---