Implicit solvers

Equation (7-54) allows calculation of the residence time required to achieve a given conversion or effluent composition. In the case of a network of reactions, knowing the reaction rates as a function of volumetric concentrations allows solution of the set of often nonlinear algebraic material balance equations using an implicit solver such as the multi variable Newton-Raphson method to determine the CSTR effluent concentration as a function of the residence time. As for batch reactors, for a single reaction all compositions can be expressed in terms of a component conversion or volumetric concentration, and Eq. (7-54) then becomes a single nonlinear algebraic equation solved by the Newton-Raphson method (for more details on this method see the relevant section this handbook). 

A more efficient simulation tool would be to use an implicit solver, but that solver would require improved implicit contact algorithms. A scheme for extending the capabilities of the simulation from using a pin-jointed analogy to consider yarn crossovers to allowing for yam sliding would be another positive step in the simulation moving closer to the observed mechanical behavior of fabrics. 

Operations with the Jacobian matrix are known to consume high computational time in most simulations involving implicit solvers. Based on the sensitivity analysis of the Jacobian matrix, around two-third of reactions and one-third of chemical species can be eliminated from the complete mechanism without significant loss in the quality of results. 

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. 

ODE solver. Relative to non-stiff ODE solvers, stiff ODE solvers typically use implicit methods, which require the numerical inversion of an Ns x Ns Jacobian matrix, and thus are considerably more expensive. In a transported PDF simulation lasting T time units, the composition variables must be updated /Vsm, = T/At 106 times for each notional particle. Since the number of notional particles will be of the order of A p 106, the total number of times that (6.245) must be solved during a transported PDF simulation can be as high as A p x A sim 1012. Thus, the computational cost associated with treating the chemical source term becomes the critical issue when dealing with detailed chemistry. 

For these time periods, the ODEs and active algebraic constraints influence the state and control variables. For these active sets, we therefore need to be able to analyze and implicitly solve the DAE system. To represent the control profiles at the same level of approximation as for the state profiles, approximation and stability properties for DAE (rather than ODE) solvers must be considered. Moreover, the variational conditions for problem (16), with different active constraint sets over time, lead to a multizone set of DAE systems. Consequently, the analogous Kuhn-Tucker conditions from (27) must have stability and approximation properties capable of handling ail of these DAE systems. 

M is a singular Matrix. Zero entries on the main diagonal of this matrix identify the algebraic equations, and all other entries which have the value 1 represent the differential equation. The vector x describes the state of the system. As numeric tools for the solution of the DAE system, MATLAB with the solver odel5s was used. In this solver, a Runge Kutta procedure is coupled with a BDF procedure (Backward Difference Formula). An implicit numeric scheme is used by the solver. 

Formally, the above process is equivalent to (6.4), extended for any n and solving that system. The u-v device is a more efficient way of solving it than any linear equation solver that might otherwise have been used, as n becomes larger. The u-v device will be extensively used in this book, even with implicit methods for coupled equation systems, where we must solve for a number of concentration profiles (see below). There are practitioners who believe that n = 2, that is the two-point G-approximation, is good enough. This is justified in cases where H is very small, as it often is, at least near the electrode, when unequal intervals are used (see Chap. 9). In that case, one can simply use (6.5). 

The Colebrook equation is implicit inf, and thus the determination of the friction factor requires some iteration unless an equation solver such as EES is used. All approximate explicit relation for/was given by S. E. Haaland in 1983 as... 

Equation 1 provides a comprehensive description of the equilibrium on liquid-liquid interfaces, but it cannot be solved explicitly for A cp without some assumption when more than two ions are involved. For more than two ions the implicit form can be solved numerically. We used a TK-Plus Solver (Universal Technical Systems, Rockford, IL) software for a personal computer. 

A few examples will be demonstrated in the following section. Each feature was incorporated in the software because it has been found necessary or useful in some modelling project. All the modelling tasks that previously required tailor made solutions in each project can now be solved in a unified manner in the ModEst environment. ModEst is able to deal with explicit algebraic, implicit algebraic (systems of nonlinear equations) and ordinary differential equations. As the standard way to handle PDE systems, the Numerical Method of Lines, which transforms a PDE system to a number of ODE components is used. In addition, any model with a solver provided may be dealt with as an algebraic system. The basic numerical tools are contained in the well tested public domain software (Bias, Linpack, Eispack, LSODE). 

---