Complex column simulations

More complex situations where ideal behaviour can no longer be assumed require the incorporation of activity coefficient terms in the calculation of the equilibrium vapour compositions. Assuming ideal behaviour in the gas phase, the equilibrium relation for component i is 

Multicomponent equilibria combined with distillation heat effects are discussed in more detail in Section 3.3.4 below. 

---

Schneider R, Kenig EY, Gorak A. Complex reactive absorption processes model optimization and dynamic column simulation. Proceedings of ESCAPE 11 European Symposium on Computer-Aided Process Engineering, Kolding, Denmark, 2001. 

Several of the commercial simulation programs offer preconfigured complex column rigorous models for petroleum fractionation. These models include charge heaters, several side strippers, and one or two pump-around loops. These fractionation column models can be used to model refinery distillation operations such as crude oil distillation, vacuum distillation of atmospheric residue oil, fluidized catalytic cracking (FCC) process main columns, and hydrocracker or coker main columns. Aspen Plus also has a shortcut fractionation model, SCFrac, which can be used to configure fractionation columns in the same way that shortcut distillation models are used to initialize multicomponent rigorous distillation models. 

It is easier to formulate and converge a computer simulation for a sequence of simple columns than a complex column. 

To redistribute the stages in the remaining sections, a shortcut simulation is used to find out the required number of trays, the feed tray location and the minimum reflux ratio for each column in the sequence. To make use of the existing column with the same number of trays (24 trays) iterations are required to adjust the sum of the rectifying sections in each column equal to 24 (number of trays in the main column). Finally, the sequence of the simple columns is merged into a complex column. The main column is not changed, but the side strippers and pump arounds need to be relocated or adjusted. 

The dynamic model proposed proved to represent well the separation in a DWC of a ternary hydrocarbon mixture. The values of internal flows and temperature distributions along the trays reached at steady state were in good agreement with the simulations obtained in the frame of commercial simulators. The use as control variables the reflux ratio or the side-stream flowrate proved to enable a reduction of the startup time with about 70 % compared with classical startup procedures. The complex technique developed can be a useful tool in studying dynamic behavior and startup optimization for complex columns and can be easily extended to various mixtures. 

The equilibrium-stage model can be made more realistic by including tray efficiencies (see Sect. 5.4.1.4) in the E-equations (5.2-68). Another possibility for improving column simulation is the use of the concept of transfer units (see Sect. 5.2.1.1) instead of the concept of equihbrium stages. This concept is well suited for packed columns. Such models, called rate-based models, simultaneously solve the relevant thermodynamic and mass transfer equations describing the complex mechanisms in a column (e.g., Taylor et al. 1994 Kloecker et al. 2005). Thus, rate-based models are much more complex since, for instance, the mass transfer coefficients and the interfacial area (see Sect. 5.4.1.4) must a priori be known. 

In the example distillation system considered in Chapters 3 and 4, we studied the binary propane/isobutane separation in a single distillation column. This is a fairly ideal system from the standpoint of vapor-liquid equilibrium (VLE), and it has only two components, a single feed and two product streams. In this chapter, we will show that the steady-state simulation methods can be extended to multicomponent nonideal systems and to more complex column configurations. 

Now if the reflux is increased back to 0.09 kmol/s, the column does not converge to the same steady state that it had previously at this flow rate. The flow rate must be increased to about 0.11 kmol/s to reestablish the desired low water content in the bottoms. This multiple steady-state phenomenon is one of the severe complexities that simulations of distillation columns experience when highly nonideal VLLE relationships are involved. 

The MMSP programs were developed for the simulation, design, and control of Multicomponent Multistage Reparation Processes. These include the operations listed in Figure 7. These processes may be carried out in conventional columns (one feed and two products), in complex columns (multiple feeds, multiple products, and multiple interstage coolers or heaters), or in a series of conventional or complex columns. 

General Properties of Computerized Physical Property System. Flow-sheeting calculations tend to have voracious appetites for physical property estimations. To model a distillation column one may request estimates for chemical potential (or fugacity) and for enthalpies 10,000 or more times. Depending on the complexity of the property methods used, these calculations could represent 80% or more of the computer time requited to do a simulation. The design of the physical property estimation system must therefore be done with extreme care. 

A special mention is in order of high-resolution electron microscopy (HREM), a variant that permits columns of atoms normal to the specimen surface to be imaged the resolution is better than an atomic diameter, but the nature of the image is not safely interpretable without the use of computer simulation of images to check whether the assumed interpretation matches what is actually seen. Solid-state chemists studying complex, non-stoichiometric oxides found this image simulation approach essential for their work. The technique has proved immensely powerful, especially with respect to the many types of defect that are found in microstructures. 

The principle of the perfectly-mixed stirred tank has been discussed previously in Sec. 1.2.2, and this provides essential building block for modelling applications. In this section, the concept is applied to tank type reactor systems and stagewise mass transfer applications, such that the resulting model equations often appear in the form of linked sets of first-order difference differential equations. Solution by digital simulation works well for small problems, in which the number of equations are relatively small and where the problem is not compounded by stiffness or by the need for iterative procedures. For these reasons, the dynamic modelling of the continuous distillation columns in this section is intended only as a demonstration of method, rather than as a realistic attempt at solution. For the solution of complex distillation problems, the reader is referred to commercial dynamic simulation packages. 

Wiggins et al. [456] used neutrons from the thermal column of a 10 kW pool-type research reactor and from a 120 pg Cf source to study the prompt photon emission resulting from neutron capture in magnesium nodules (ter-romanganese oxides) from the ocean floor. Spectra were recorded with a Ce(Ii) detector and a 1024-channel analyser. Complex spectra were obtained by irradiation of seawater, but it was possible to detect and estimate manganese in nodules in a simulated marine environment by means of the peaks at 7.00, 6.55, 6.22, and 6.04 pV. 

In system 1, the 3-D dynamic bubbling phenomena in a gas liquid bubble column and a gas liquid solid fluidized bed are simulated using the level-set method coupled with an SGS model for liquid turbulence. The computational scheme in this study captures the complex topological changes related to the bubble deformation, coalescence, and breakup in bubbling flows. In system 2, the hydrodynamics and heat-transfer phenomena of liquid droplets impacting upon a hot flat surface and particle are analyzed based on 3-D level-set method and IBM with consideration of the film-boiling behavior. The heat transfers in... 

---