Related Differential Equations

The application of the reverse Euler method of solution to a system of coupled differential equations yields a system of coupled algebraic equations that can be solved by the method of Gaussian elimination and back substitution. In this chapter I demonstrated the solution of simultaneous algebraic equations by means of this method and showed how the solution of algebraic equations can be used to solve the related differential equations. In the process, I presented subroutine GAUSS, the computational engine of all of the programs discussed in the chapters that follow. 

Equations (5.95), (5.96) and (5.97) are suitable for constant critical melting porosity. In a one dimensional steady state melting column as a result of decompression melting, the porosity may increase from the bottom to the top of the column. If melting porosities change as a function of the spatial position, the related differential equations need to be solved numerically. More details of various melt transport models by porous flow have been given by Spiegelman and Elliott (1993), Iwamori, (1994), and Lundstrom (2000). 

Elementary creep and relaxation tests as a means to experimentally characterize polymers were discussed in Chapter 3. Further, elementary mechanical models and the related differential equations were discussed as... 

Comparison of various viscoelastic models and related differential equations, rrstress parameter, e strain parameter, E elastic parameter, i viscous parameter, X. relaxation time... 

It turns out that there is another branch of mathematics, closely related to tire calculus of variations, although historically the two fields grew up somewhat separately, known as optimal control theory (OCT). Although the boundary between these two fields is somewhat blurred, in practice one may view optimal control theory as the application of the calculus of variations to problems with differential equation constraints. OCT is used in chemical, electrical, and aeronautical engineering where the differential equation constraints may be chemical kinetic equations, electrical circuit equations, the Navier-Stokes equations for air flow, or Newton s equations. In our case, the differential equation constraint is the TDSE in the presence of the control, which is the electric field interacting with the dipole (pemianent or transition dipole moment) of the molecule [53, 54, 55 and 56]. From the point of view of control theory, this application presents many new features relative to conventional applications perhaps most interesting mathematically is the admission of a complex state variable and a complex control conceptually, the application of control teclmiques to steer the microscopic equations of motion is both a novel and potentially very important new direction. 

The fimdamental kinetic master equations for collisional energy redistribution follow the rules of the kinetic equations for all elementary reactions. Indeed an energy transfer process by inelastic collision, equation (A3.13.5). can be considered as a somewhat special reaction . The kinetic differential equations for these processes have been discussed in the general context of chapter A3.4 on gas kmetics. We discuss here some special aspects related to collisional energy transfer in reactive systems. The general master equation for relaxation and reaction is of the type [H, 12 and 13, 15, 25, 40, 4T ] ... 

Hence, in order to contract extended BO approximated equations for an N-state coupled BO system that takes into account the non-adiabatic coupling terms, we have to solve N uncoupled differential equations, all related to the electronic ground state but with different eigenvalues of the non-adiabatic coupling matrix. These uncoupled equations can yield meaningful physical... 

This part of our chapter has shown that the use of the two variables, moduli and phases, leads in a direct way to the derivation of the continuity and Hamilton-Jacobi equations for both scalar and spinor wave functions. For the latter case, we show that the differential equations for each spinor component are (in the nearly nomelativistic limit) approximately decoupled. Because of this decoupling (mutual independence) it appears that the reciprocal relations between phases and moduli derived in Section III hold to a good approximation for each spinor component separately, too. For velocities and electromagnetic field strengths that ate nomrally below the relativistic scale, the Berry phase obtained from the Schrddinger equation (for scalar fields) will not be altered by consideration of the Dirac equation. 

Another difference is related to the mathematical formulation. Equation (1) is deterministic and does not include explicit stochasticity. In contrast, the equations of motion for a Brownian particle include noise. Nevertheless, similar algorithms are adopted to solve the two differential equations as outlined below. The most common approach is to numerically integrate the above differential equations using small time steps and preset initial values. 

For example, the SHAKE algorithm [17] freezes out particular motions, such as bond stretching, using holonomic constraints. One of the differences between SHAKE and the present approach is that in SHAKE we have to know in advance the identity of the fast modes. No such restriction is imposed in the present investigation. Another related algorithm is the Backward Euler approach [18], in which a Langevin equation is solved and the slow modes are constantly cooled down. However, the Backward Euler scheme employs an initial value solver of the differential equation and therefore the increase in step size is limited. 

In the non-linear differential equation Eq. (43), k is related to the inverse Debye-Hiickel length. The method briefly outlined above is implemented, e.g., in the pro-... 

In unsteady states the situation is less satisfactory, since stoichiometric constraints need no longer be satisfied by the flux vectors. Consequently differential equations representing material balances can be constructed only for binary mixtures, where the flux relations can be solved explicitly for the flux vectors. This severely limits the scope of work on the dynamical equations and their principal field of applicacion--Che theory of stability of steady states. The formulation of unsteady material and enthalpy balances is discussed in Chapter 12, which also includes a brief digression on stability problems. 

It is only for smooth field models, in this sense, that partial differential equations relating species concentrations to position in space can be written down. However, a pore geometry which is consistent with the smooth... 

To reduce the material balance conditions (11,1) to differential equations for the composition and pressure, flux relations must be used to relate the vectors to the gradients of the composition and pressure... 

They represent the generalization of our earlier equations (11.12) and include an additional equation for the pressure, since this can no longer be assumed to take the constant value p throughout the pellet. The differential equation relating X to f is derived exactly as before and has the same form, except that p cannot be replaced by p In the present case. Thus we have... 

Of course the algebraic stoichiometric relations (11.3) must be distinguished from the corresponding differential equations (11.4) which are always satisfied in steady states, since they are direct consequences of... 

The bead and spring model is clearly based on mechanical elements just as the Maxwell and Voigt models were. There is a difference, however. The latter merely describe a mechanical system which behaves the same as a polymer sample, while the former relates these elements to actual polymer chains. As a mechanical system, the differential equations represented by Eq. (3.89) have been thoroughly investigated. The results are somewhat complicated, so we shall not go into the method of solution, except for the following observations ... 

The vertical displacements w are described by the fourth order differential equation according to the equilibrium and the constitutive laws. The following relations for w,... 

The analysis of steady-state and transient reactor behavior requires the calculation of reaction rates of neutrons with various materials. If the number density of neutrons at a point is n and their characteristic speed is v, a flux effective area of a nucleus as a cross section O, and a target atom number density N, a macroscopic cross section E = Na can be defined, and the reaction rate per unit volume is R = 0S. This relation may be appHed to the processes of neutron scattering, absorption, and fission in balance equations lea ding to predictions of or to the determination of flux distribution. The consumption of nuclear fuels is governed by time-dependent differential equations analogous to those of Bateman for radioactive decay chains. The rate of change in number of atoms N owing to absorption is as follows ... 

There are special numerical analysis techniques for solving such differential equations. New issues related to the stabiUty and convergence of a set of differential equations must be addressed. The differential equation models of unsteady-state process dynamics and a number of computer programs model such unsteady-state operations. They are of paramount importance in the design and analysis of process control systems (see Process control). 

The natural laws in any scientific or technological field are not regarded as precise and definitive until they have been expressed in mathematical form. Such a form, often an equation, is a relation between the quantity of interest, say, product yield, and independent variables such as time and temperature upon which yield depends. When it happens that this equation involves, besides the function itself, one or more of its derivatives it is called a differential equation. 

A relation between the variables, involving no derivatives, is called a solution of the differential equation if this relation, when substituted in the equation, satisfies the equation. A solution of an ordinaiy differential equation which includes the maximum possible number of arbitrary constants is called the general solution. The maximum number of arbitrai y constants is exactly equal to the order of the dif-... 

A solution of a difference equation is a relation between the variables which satisfies the equation. If the difference equation is of order n, the general solution involves n arbitraty constants. The techniques for solving difference equations resemble techniques used for differential equations. 

The number of independent rate equations is the same as the number of independent stoichiometric relations. In the present example. Reactions (1) and (2) are reversible reactions and are not independent. Accordingly, C,. and C, for example, can be eliminated from the equations for and which then become an integrable system. Usually only systems of linear differential equations with constant coefficients are solvable analytically. 

Energy Laws Several laws have been proposed to relate size reduction to a single variable, the energy input to the mill. These laws are encompassed in a general differential equation (Walker, Lewis, McAdams, and Gilliland, Principles of Chemical Engineering, 3d ed., McGraw-HiU, New York, 1937) ... 

The unsteady material balances of tracer tests are represented by linear differential equations with constant coefficients that relate an input function Cj t) to a response function of the form... 

The solution of this partial differential equation is recorded in the literature (Otake and Kunigata, Kngaku Kogaku, 22, 144 [1958]). The plots of E(t ) against t are bell-shaped, resembling the corresponding Erlang plots. A relation is cited later between the Peclet number,... 

Find the differential equation relating the displaeements X[ t) and Xo t) for the spring-mass-damper system shown in Figure 2.5. What would be the effeet of negleeting the mass ... 

A flywheel of moment of inertia / sits in bearings that produee a frietional moment of C times the angular veloeity uj t) of the shaft as shown in Figure 2.7. Find the differential equation relating the applied torque T t) and the angular veloeity uj t). 

Figure 2.8 shows a reduetion gearbox being driven by a motor that develops a torque T tn(t). It has a gear reduetion ratio of and the moments of inertia on the motor and output shafts are and /q, and the respeetive damping eoeffieients Cm and Cq. Find the differential equation relating the motor torque CmfO and the output angular position 6a t). 

Find the differential equation relating v t) and viit) for the RC network shown in Figure 2.11. 

Heat flows from a heat source at temperature 6 t) through a wall having ideal thermal resistance Ri to a heat sink at temperature 62(1) having ideal thermal capacitance Ct as shown in Figure 2.14. Find the differential equation relating 6 t) and 02(0-... 

---