Convergence methods systems

Another potential advancement is permitted in the ASPEN system. Tear streams can be designated as desired, so that a user might define blocks or series of blocks and simulate these sets as quasi-linear blocks. The convergence method could utilize this information and solve the material (and energy) balances explicitly. In this way, a simultaneous modular architecture could be utilized. Implementation of these programs will be for later enhancements of ASPEN, not the initial version. 

The tandem aldol condensation-radical cyclization sequence for the elaboration of functionalized bicyclo[3.3.0]octane systems has been developed [48]. Conjugate addition of Me.AlSePh to dimethylcyclopentenone followed by trapping of the resultant enolate with aldehyde afforded the trans.erythro aldol predominantly which then underwent radical cyclization with Bu SnH and catalytic AIBN yielding the bicyclic ketol stereospecifically. This approach represents a highly convergent method for the annulation of carbocycles leading to the polyquinane sesquiterpenes. 

Exact Solutions Given by the Method of Convergence for Systems of Distillation Columns... 

The 0 method of convergence for systems is readily extended to include other types of units such as mixers and proportional dividers which are commonly found in systems of distillation columns. To demonstrate the application of the method to systems containing units such as these, consider the system shown in Fig. 3-13. Suppose the specifications for the system are taken to be D2 and B3. In addition to these, specifications such as the total-flow rate, the thermal condition and composition of the feed F, the reflux rate, the column pressure, the type of condenser, and the plate configuration for unit 2 are made. The remaining flow rates, Du D3, and B2 are computed from the set of three overall material balance equations. For any component i, the component-material balances are as follows... 

A Convergence Method for Distillation Systems, Ph.D. Dissertation, Texas A M... 

The computations shown in the left graph in Fig. 2 all started from a saturated bed. For nonlinear systems it can be expected that the performance of a method depends greatly on the the initial conditions. For System C, however, it turned out that when the computations were started from an almost empty bed, results similar to Fig. 2, left graph, were obtained. Thus the convergence of System C does not depend much on the initial conditions. For System D this will be no longer the case. 

For choices of 0 < < 1, we have under-relaxation methods, which are successful for some systems that are not convergent for Gauss-Seidel. Those methods associated with > 1 are called over-relaxation methods and are useful in accelerating the convergence for systems that are already convergent by Gauss-Seidel. These over-relaxation methods are also named successive over-relaxation (SOR), and find application in the numerical solution of certain partial differential equations. 

As a process simulator, we used Aspen HYSYS (see Figure 8.15). In all the scenarios, as initial values, we use the stream values obtained when the simulation is done in open loop. Moreover, instead of using the simnlator tools for converging the system (the recycle unit operation in HYSYS), we connect the simulator with external modules developed in MATLAB . In that way, we have a complete control over the numerical methods used for converging the system. In all the cases, a termination tolerance is eqnal to 10 , using a norm 1. 

The convergence method used in this ternary case with inerts is outlined in the following list and is quite different that those used in the other systems. 

These systems are solved by a step-limited Newton-Raphson iteration, which, because of its second-order convergence characteristic, avoids the problem of "creeping" often encountered with first-order methods (Law and Bailey, 1967) ... 

The Hemian-Kluk method has been developed further [153-155], and used in a number of applications [156-159]. Despite the formal accuracy of the approach, it has difficulties, especially if chaotic regions of phase space are present. It also needs many trajectories to converge, and the initial integration is time consuming for large systems. Despite these problems, the frozen Gaussian approximation is the basis of the spawning method that has been applied to... 

---