Mathematical modelling and control of agitated extraction columns

Liquid-liquid extraction has long been known as a unit operation that plays an important role in industry. This process is well known for its complexity and sensitivity to operation conditions. This thesis presents an attempt to explore the dynamics and control of this process using a systematic approach and state of the art control system design techniques. The process was studied first experimentally under carefully selected. operation conditions, which resembles the ranges employed practically under stable and efficient conditions. Data were collected at steady state conditions using adequate sampling techniques for the dispersed and continuous phases as well as during the transients of the column with the aid of a computer-based online data logging system and online concentration analysis. A stagewise single stage backflow model was improved to mimic the dynamic operation of the column. The developed model accounts for the variation in hydrodynamics, mass transfer, and physical properties throughout the length of the column. End effects were treated by addition of stages at the column entrances. Two parameters were incorporated in the model namely; mass transfer weight factor to correct for the assumption of no mass transfer in the. settling zones at each stage and the backmixing coefficients to handle the axial dispersion phenomena encountered in the course of column operation. The parameters were estimated by minimizing the differences between the experimental and the model predicted concentration profiles at steady state conditions using non-linear optimisation technique. The estimated values were then correlated as functions of operating parameters and were incorporated in·the model equations. The model equations comprise a stiff differential~algebraic system. This system was solved using the GEAR ODE solver. The calculated concentration profiles were compared to those experimentally measured. A very good agreement of the two profiles was achieved within a percent relative error of ±2.S%. The developed rigorous dynamic model of the extraction column was used to derive linear time-invariant reduced-order models that relate the input variables (agitator speed, solvent feed flowrate and concentration, feed concentration and flowrate) to the output variables (raffinate concentration and extract concentration) using the asymptotic method of system identification. The reduced-order models were shown to be accurate in capturing the dynamic behaviour of the process with a maximum modelling prediction error of I %. The simplicity and accuracy of the derived reduced-order models allow for control system design and analysis of such complicated processes. The extraction column is a typical multivariable process with agitator speed and solvent feed flowrate considered as manipulative variables; raffinate concentration and extract concentration as controlled variables and the feeds concentration and feed flowrate as disturbance variables. The control system design of the extraction process was tackled as multi-loop decentralised SISO (Single Input Single Output) as well as centralised MIMO (Multi-Input Multi-Output) system using both conventional and model-based control techniques such as IMC (Internal Model Control) and MPC (Model Predictive Control). Control performance of each control scheme was. studied in terms of stability, speed of response, sensitivity to modelling errors (robustness), setpoint tracking capabilities and load rejection. For decentralised control, multiple loops were assigned to pair.each manipulated variable with each controlled variable according to the interaction analysis and other pairing criteria such as relative gain array (RGA), singular value analysis (SVD). Loops namely Rotor speed-Raffinate concentration and Solvent flowrate Extract concentration showed weak interaction. Multivariable MPC has shown more effective performance compared to other conventional techniques since it accounts for loops interaction, time delays, and input-output variables constraints.

A study of the dynamic behaviour of a pulsed column solvent extraction plant

A study has been made of the dynamic behaviour of a nuclear fuel reprocessing plant utilising pulsed solvent extraction columns. A flowsheet is presented and the choice of an extraction device is discussed. The plant is described by a series of modules each module representing an item of equipment. Each module consists of a series of differential equations describing the dynamic behaviour of the equipment. The model is written in PMSP, a language developed for dynamic simulation models. The differential equations are solved to predict plant behaviour with time. The dynamic response of the plant to a range of disturbances has been assessed. The interactions between pulsed columns have been demonstrated and illustrated. The importance of auxillary items of equipment to plant performance is demonstrated. Control of the reprocessing plant is considered and the effect of control parameters on performance assessed.

A study of liquid-liquid extraction in a sieve plate column

The literature relating to sieve plate liquid extraction columns and relevant hydrodynamic phenomena have been surveyed. Mass transfer characteristics during drop formation, rise and coalescence, and related models were also reviewed. Important design parameters i.e. flooding, dispersed phase hold-up, drop size distribution, mean drop size, coalescence/flocculation zone height beneath a plate and jetting phenomena were investigated under non-mass transfer and mass transfer conditions in a 0.45m diameter, 2.3m high sieve plate column. This column had provision for four different plate designs, and variable plate spacing and downcomer heights, and the system used was Clairsol `350' (dispersed) - acetone - deionised water (continuous) with either direction of mass transfer. Drop size distributions were best described by the functions proposed by Gal-or, and then Mugele-Evans. Using data from this study and the literature, correlations were developed for dispersed phase hold-up, mean drop size in the preferred jetting regime and in the non-jetting regime, and coalescence zone height. A method to calculate the theoretical overall mass transfer coefficient allowing for the range of drop sizes encountered in the column gave the best fit to experimental data. This applied the drop size distribution diagram to estimate the volume percentage of stagnant, circulating and oscillating drops in the drop population. The overall coefficient Kcal was then calculated as the fractional sum of the predicted individual single drop coefficients and their proportion in the drop population. In a comparison between the experimental and calculated overall mass transfer coefficients for cases in which all the drops were in the oscillating regime (i.e. 6.35mm hole size plate), and for transfer from the dispersed(d) to continuous(c) phase, the film coefficient kd predicted from the Rose-Kintner correlation together with kc from that of Garner-Tayeban gave the best representation. Droplets from the 3.175mm hole size plate, were of a size to be mainly circulating and oscillating; a combination of kd from the Kronig-Brink (circulating) and Rose-Kintner (oscillating) correlations with the respective kc gave the best agreement. The optimum operating conditions for the SPC were identified and a procedure proposed for design from basic single drop data.

The hydrodynamic behaviour and mass transfer characteristics of single droplets in a pulsed sieve plate column

The literature relating to the performance of pulsed sieve plate liquid-liquid extraction columns and the relevant hydrodynamic phenomenon have been surveyed. Hydrodynamic behaviour and mass transfer characteristics of droplets in turbulent and non-turbulent conditions have also been reviewed. Hydrodynamic behaviour, i.e. terminal and characteristic velocity of droplets, droplet size and droplet breakup processes, and mass transfer characteristics of single droplets (d≤0.6 cm) were investigated under pulsed (mixer-settler & transitional regimes) and non-pulsed conditions in a 5.0 cm diameter, 100 cm high, pulsed sieve plate column with three different sieve plate types and variable plate spacing. The system used was toluene (displaced) - acetone - distilled water. Existing photographic techniques for following and recording the droplet behaviour, and for observing the parameters of the pulse and the pulse shape were further developed and improved. A unique illumination technique was developed by which a moving droplet could be photographed using cine or video photography with good contrast without using any dye. Droplet size from a given nozzle and droplet velocity for a given droplet diameter are reduced under pulsing condition, and it was noted that this effect is enhanced in the presence of sieve plate. The droplet breakup processes are well explained by reference to an impact-breakup mechanism. New correlations to predict droplet diameter based on this mechanism are given below.vskip 1.0cm or in dimensionless groups as follows:- (We)crit= 3.12 - 1.79 (Eo)crit A correlation based on the isotropic turbulence theory was developed to calculate droplet diameter in the emulsion regime.vskip 1.0cm Experimental results show that in the mixer-settler and transitional regimes, pulsing parameters had little effect on the overall dispersed phase mass transfer coefficient during the droplet formation and unhindered travel periods.

BTX extraction by sulfolane in a mixer-settler cascade

The literature relating to the extraction of the aromatics, benzene, toluene and xylene (BTX) using different commercial solvents, and to mixer-settler design and performance, has been reviewed. Liquid-liquid equilibria of the ternary systems: hexane-benzene-sulfolane, n-heptane-toluene-sulfolane, and octane-xylene-sulfolane were determined experimentally at temperatures of 30oC, 35oC, and 40oC. The work was then extended to a multicomponent system. The data were correlated by using Hand's method and were found to be in a good agreement with theoretical predictions using the UNIFAC method. A study was made of the performance of a 10-stage laboratory mixer-settler cascade for the extraction of BTX from a synthetic reformate utilizing sulfolane as a solvent. Murphree stage efficiency decreased with stage number but 99% extraction was achievable within 4 stages. The effects of temperature, phase ratio, and agitator speed were investigated. The efficiency increased with agitator speed but > 1050 rpm resulted in secondary haze formation. An optimum temperature of 30oC was selected from the phase equilibria; the optimum solvent: feed ratio was 3:1 for 4 stages. The experimental overall mass transfer coefficients were compared with those predicted from single drop correlations and were in all cases greater, by a factor of 1.5 to 3, due to the surface renewal associated with drop break-up and coalescence promoted by agitation. A similar investigation was performed using real reformate from the Kuwait Oil Company. The phase ratios were in the range 0.5 to 1 to 3.25 to 1, the agitator speed 1050 rpm, and the operating temperature 30oC. A maximum recovery of 99% aromatics was achieved in 4 stages at a phase ratio of 3.25 to 1. A backflow model was extended to simulate conditions in the mixer-settler cascade with this multicomponent system. Overall mass transfer coefficients were estimated by obtaining the best fit between experimental and predicted concentration profiles. They were up to 10% greater than those with the synthetic feed but close agreement was not possible because the distribution coefficient and phase ratio varied with stage number. Sulfolane was demonstrated to be an excellent solvent for BTX recovery and a mixer-settler cascade was concluded to be a technically viable alternative to agitated columns for this process.