Solving the NEQ Model Equations

Solving the NEQ Model Equations In general, a nonequilibrium model of a column has many more equations than does an equivalent equilibrium-stage model. Nevertheless, we use may essentially the same computational approaches to solve the nonequilibrium model equations simultaneous convergence (Krishnamurthy and Taylor, op. cit.) and continuation methods [Powers et al., Comput Chem. Engng., 12, 1229 (1988)]. Convergence of a nonequilibrium model is likely to be slower than that of the equilibrium model because of the greater... 

The NEQ model requires thermodynamic properties, not only for calculation of phase equilibrium but also for calculation of driving forces for mass transfer. In addition, physical properties such as surface tension, diffusion coefficients, and viscosities, for calculation of mass (and heat) transfer coefficients and interfacial areas are required. The steady-state model equations most often are solved using Newton s method or by homotopy-continuation. A review of early applications of NEQ models is available [5]. 

Lee and Dudukovic [18] described an NEQ model for homogeneous RD in tray columns. The Maxwell-Stefan equations are used to describe interphase transport, with the AIChE correlations used for the binary (Maxwell-Stefan) mass-transfer coefficients. Newton s method and homotopy continuation are used to solve the model equations. Close agreement between the predictions of EQ and NEQ models were found only when the tray efficiency could correctly be predicted for the EQ model. In a subsequent paper Lee and Dudukovic [19] presented a dynamic NEQ model of RD in tray columns. The DAE equations were solved by use of an implicit Euler method combined with homotopy continuation. Murphree efficiencies calculated from the results of an NEQ simulation of the production of ethyl acetate were not constant with time. 

---