Solution Procedures

The constraints given above complete the mathematical formulation for both cases considered. The constraints given above contain a number of nonlinearities, which complicates the solution of resulting models. This was dealt with through the solution procedure discussed below. 

One would notice that there are a number of nonlinearities in the above constraints, more specifically in the contaminant mass balances around a unit and the central storage vessel. The nonlinearities arise due to the fact that the outlet concentration of each contaminant may not necessarily be at its respective maximum. Unlike the single contaminant case where one could replace the outlet concentration with the maximum outlet concentration, in the multiple contaminant case the outlet concentration of each contaminant remains a variable. Furthermore, the concentration within the central storage vessel is always variable, since the contaminant mass and mass of water within the vessel changes each time a stream enters or exits the vessel. To deal with this situation the following procedure is considered. 

In this section the application of the multiple contaminant methodology is demonstrated though a number of illustrative examples. The first example is solved for both the case where there is no central storage vessel and the case where there is a central storage vessel. The second example is an adapted literature example. Due to the nature of the example it was only solved considering a central storage vessel. 

The axial temperature profiles can be developed beginning from either end of the kiln by means of Equations (8.30) and (8.33). However, because direct-fired kilns operate in a countercurrent flow mode, we can employ a shooting method (Tscheng and Watkinson, 1979) to 

At this point, the finite difference method described for convective-diffusion equations (Patanker, 1980) might be used for the solution of the governing bed/wall equations derived for the kiln cross section 

Chapter 8 Heat Transfer Processes in the Rotary Kiln Bed 

The heat diffusion flux, D, the convection flux, F, and the associated Peclet number are 

Within a numerical model, the function A P is represented by a power law scheme as (Patanker, 1980) 

---

Fig. 1. The Split Integration Symplectic Method (SISM) solution procedure.

In principle, Chen, given the flux relations there is no difficulty in constructing differencial equations to describe the behavior of a catalyst pellet in steady or unsteady states. In practice, however, this simple procedure is obstructed by the implicit nature of the flux relations, since an explicit solution of usefully compact form is obtainable only for binary mixtures- In steady states this impasse is avoided by using certain, relations between Che flux vectors which are associated with the stoichiometry of Che chemical reaction or reactions taking place in the pellet, and the major part of Chapter 11 is concerned with the derivation, application and limitations of these stoichiometric relations. Fortunately they permit practicable solution procedures to be constructed regardless of the number of substances in the reaction mixture, provided there are only one or two stoichiomeCrically independent chemical reactions. 

Though the solution procedure sounds straightforward, if tedious, practice difficulty is encountered immediately because of the implicit nature of the available flux models. As we saw in Chapter 5 even the si lest of these, the dusty gas model, has solutions which are too cumbersc to be written down for more than three components, while the ternary sol tion itself is already very complicated. It is only for binary mixtures therefore, that the explicit formulation and solution of equations (11. Is practicable. In systems with more than two components, we rely on... 

All numerical computations inevitably involve round-off errors. This error increases as the number of calculations in the solution procedure is increased. Therefore, in practice, successive mesh refinements that increase the number of finite element calculations do not necessarily lead to more accurate solutions. However, one may assume a theoretical situation where the rounding error is eliminated. In this case successive reduction in size of elements in the mesh should improve the accuracy of the finite element solution. Therefore, using a P C" element with sufficient orders of interpolation and continuity, at the limit (i.e. when element dimensions tend to zero), an exact solution should be obtaiiied. This has been shown to be true for linear elliptic problems (Strang and Fix, 1973) where an optimal convergence is achieved if the following conditions are satisfied ... 

Nonlinear Programming The most general case for optimization occurs when both the objective function and constraints are nonlinear, a case referred to as nonlinear programming. While the idea behind the search methods used for unconstrained multivariable problems are applicable, the presence of constraints complicates the solution procedure. 

The variational energy principles of classical elasticity theory are used in Section 3.3.2 to determine upper and lower bounds on lamina moduli. However, that approach generally leads to bounds that might not be sufficiently close for practical use. In Section 3.3.3, all the principles of elasticity theory are invoked to determine the lamina moduli. Because of the resulting complexity of the problem, many advanced analytical techniques and numerical solution procedures are necessary to obtain solutions. However, the assumptions made in such analyses regarding the interaction between the fibers and the matrix are not entirely realistic. An interesting approach to more realistic fiber-matrix interaction, the contiguity approach, is examined in Section 3.3.4. The widely used Halpin-Tsai equations are displayed and discussed in Section 3.3.5. 

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. 

Solution of the equations in (2.8-2.18) proceeds with an adaptive nonlinear boundary value method. The solution procedure has been discussed in detail elsewhere (10) and we outline only the essential features here. Our goal is to obtain a discrete solution of the governing equations on the mesh Af... 

In this section we apply the adaptive boundary value solution procedure and the pseudo-arclength continuation method to a set of strained premixed hydrogen-air flames. Our goal is to predict accurately and efficiently the extinction behavior of these flames as a function of the strain rate and the equivalence ratio. Detailed transport and complex chemical kinetics are included in all of the calculations. The reaction mechanism for the hydrogen-air system is listed in Table... 

Franks has suggested that the solution procedure should be based on the concentration value, Yi orX], whichever having the greatest magnitude. Thus for a linear equilibrium in which Yj = m X], the equilibrium equation can be differentiated to give... 

As discussed by Franks (1972), in order to solve this system of equations, a value of temperature T must be found to satisfy the condition that the difference term 6 = P - Zpj is very small, i.e., that the equilibrium condition is satisfied. This is known as a bubble point calculation. The above system of defining equations, however represent, an implicit algebraic loop and the trial and error solution procedure can be very time consuming, especially when incorporated into a dynamic simulation program. 

Solution of the required column height is achieved by integrating the two component balance equations and the heat balance equation, down the column from the known conditions Xi , yout and TLin, until the condition that either Y is greater than or X is greater than Xqui is achieved. In this solution approach, variations in the overall mass transfer capacity coefficient both with respect to temperature and to concentration, if known, can also be included in the model as required. The solution procedure is illustrated by the simulation example AMMON AB. 

The reaction is endothermic, and heat transfer to the reactor is required in order to accomplish the decomposition of the acetylated oil, to liberate acetic acid vapour. The example has been considered previously by Perona (1972), Smith (1972), Cooper and Jeffreys (1971) and Froment and Bischoff (1990), although their solution procedures differ from that here. Data values are based on those used by Froment and Bischoff. 

Especially for the electrons, the fluid model has the advantage of a lower computational effort than the PIC/MC method. Their low mass (high values of the transport coefficients) and consequent high velocities give rise to small time steps in the numerical simulation (uAf < Aa) if a so-called explicit method is used. This restriction is easily eliminated within the fluid model by use of an implicit method. Also, the electron density is strongly coupled with the electric field, which results in numerical Instabilities. This requires a simultaneous implicit solution of the Poisson equation for the electric field and the transport equation for the electron density. This solution can be deployed within the fluid model and gives a considerable reduction of computational effort as compared to a nonsi-multaneous solution procedure [179]. Within the PIC method, only fully explicit methods can be applied. 

Principles and Characteristics Fractional solution procedures usually consist of consecutive extractions with solvents of increasing solvent power. These labour intensive methods benefit from a larger surface area to mass ratio. Other methods for fractionation by solubility rely on fractional precipitation through addition of a nonsolvent, lowering the temperature or solvent volatilisation (Section 3.7). 

MULTICOMPONENT SYSTEMS RIGOROUS SOLUTION PROCEDURES (COMPUTER METHODS)... 

The basic steps in any rigorous solution procedure will be ... 

Without a solution, formulated mathematical systems (models) are of little value. Four solution procedures are mainly followed the analytical, the numerical (e.g., finite different, finite element), the statistical, and the iterative. Numerical techniques have been standard practice in soil quality modeling. Analytical techniques are usually employed for simplified and idealized situations. Statistical techniques have academic respect, and iterative solutions are developed for specialized cases. Both the simulation and the analytic models can employ numerical solution procedures for their equations. Although the above terminology is not standard in the literature, it has been used here as a means of outlining some of the concepts of modeling. 

The most representative characteristics are given, since for example statistical formulations can be subject to statistical analytic or numerical solution procedures. 

The solution procedure for the problem stated above involves two optimisation stages in which freshwater and reusable water storage capacity are minimised sequentially, as explained in Section 5.4.1 below. 

---