Algebraic equations, numerical solution

One can always solve these differential or algebraic equations numerically for any values of the kjS, but analytical solutions assuming aU ks equal still provide useful first approximations. 

Note that fugacities are calculated on the basis of thermodynamic equations with use of one of the state equations. The system of equations (5.169)-(5.172) represents nonlinear system of algebraic equations. Its solution can be found by using appropriate numerical methods. 

Families of finite elements and their corresponding shape functions, schemes for derivation of the elemental stiffness equations (i.e. the working equations) and updating of non-linear physical parameters in polymer processing flow simulations have been discussed in previous chapters. However, except for a brief explanation in the worked examples in Chapter 2, any detailed discussion of the numerical solution of the global set of algebraic equations has, so far, been avoided. We now turn our attention to this important topic. 

NUMERICAL SOLUTION OF THE GLOBAL SYSTEMS OF ALGEBRAIC EQUATIONS... 

As mentioned in Chapter 2, the numerical solution of the systems of algebraic equations is based on the general categories of direct or iterative procedures. In the finite element modelling of polymer processing problems the most frequently used methods are the direet methods. 

As mentioned earlier, overall accuracy of finite element computations is directly detennined by the accuracy of the method employed to obtain the numerical solution of the global system of algebraic equations. In practical simulations, therefore, computational errors which are liable to affect the solution of global stiffness equations should be carefully analysed. 

Opposing reactions. Calculate half-times for equilibration in the triphenyl methyl system starting with all A or with the equilibrated mixture, for the conditions given in Table 3-3. Use algebraic equations, not the tabulated numerical values. Compare the latter with the t fi from the approximate solution given in Eq. (3-39). Compare the values of 4AT-15o and a, to assess whether Eq. (3-39) provides an adequate representation. 

This equation is coupled to the component balances in Equation (3.9) and with an equation for the pressure e.g., one of Equations (3.14), (3.15), (3.17). There are A +2 equations and some auxiliary algebraic equations to be solved simultaneously. Numerical solution techniques are similar to those used in Section 3.1 for variable-density PFRs. The dependent variables are the component fluxes , the enthalpy H, and the pressure P. A necessary auxiliary equation is the thermodynamic relationship that gives enthalpy as a function of temperature, pressure, and composition. Equation (5.16) with Tref=0 is the simplest example of this relationship and is usually adequate for preliminary calculations. 

Following these procedures, we are led to a system of algebraic equations, thereby reducing numerical solution of an initial (linear) differential equation to solving an algebraic system. 

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. 

Within esqjlicit schemes the computational effort to obtain the solution at the new time step is very small the main effort lies in a multiplication of the old solution vector with the coeflicient matrix. In contrast, implicit schemes require the solution of an algebraic system of equations to obtain the new solution vector. However, the major disadvantage of explicit schemes is their instability [84]. The term stability is defined via the behavior of the numerical solution for t —> . A numerical method is regarded as stable if the approximate solution remains bounded for t —> oo, given that the exact solution is also bounded. Explicit time-step schemes tend to become unstable when the time step size exceeds a certain value (an example of a stability limit for PDE solvers is the von-Neumann criterion [85]). In contrast, implicit methods are usually stable. 

Algebraic equations Steady state of CSTR with first-order kinetics. Algebraic solution and optimisation (least squares. Draper and Smith, 1981). Steady state of CSTR with complex kinetics. Numerical solution and optimisation (least squares or likelihood function). 

Thus Y1 is obtained not as the result of the numerical integration of a differential equation, but as the solution of an algebraic equation, which now requires an iterative procedure to determine the equilibrium value, Xj. The solution of algebraic balance equations in combination with an equilibrium relation has again resulted in an implicit algebraic loop. Simplification of such problems, however, is always possible, when Xj is simply related to Yi, as for example... 

In this chapter we concentrate on dynamic, distributed systems described by partial differential equations. Under certain conditions, some of these systems, particularly those described by linear PDEs, have analytical solutions. If such a solution does exist and the unknown parameters appear in the solution expression, the estimation problem can often be reduced to that for systems described by algebraic equations. However, most of the time, an analytical solution cannot be found and the PDEs have to be solved numerically. This case is of interest here. Our general approach is to convert the partial differential equations (PDEs) to a set of ordinary differential equations (ODEs) and then employ the techniques presented in Chapter 6 taking into consideration the high dimensionality of the problem. 

Leung 1993). PBPK models for a particular substance require estimates of the chemical substance-specific physicochemical parameters, and species-specific physiological and biological parameters. The numerical estimates of these model parameters are incorporated within a set of differential and algebraic equations that describe the pharmacokinetic processes. Solving these differential and algebraic equations provides the predictions of tissue dose. Computers then provide process simulations based on these solutions. 

---