Prediction of slurry transport in SAG mills using SPH fluid flow in a dynamic DEM based porous media

DEM modelling of the motion of coarse fractions of the charge inside SAG mills has now been well established for more than a decade. In these models the effect of slurry has broadly been ignored due to its complexity. Smoothed particle hydrodynamics (SPH) provides a particle based method for modelling complex free surface fluid flows and is well suited to modelling fluid flow in mills. Previous modelling has demonstrated the powerful ability of SPH to capture dynamic fluid flow effects such as lifters crashing into slurry pools, fluid draining from lifters, flow through grates and pulp lifter discharge. However, all these examples were limited by the ability to model only the slurry in the mill without the charge. In this paper, we represent the charge as a dynamic porous media through which the SPH fluid is then able to flow. The porous media properties (specifically the spatial distribution of porosity and velocity) are predicted by time averaging the mill charge predicted using a large scale DEM model. This allows prediction of transient and steady state slurry distributions in the mill and allows its variation with operating parameters, slurry viscosity and slurry volume, to be explored. (C) 2006 Published by Elsevier Ltd.

Eliminating non-renewable CO2 emissions from sewage treatment: An anaerobic migrating bed reactor pilot plant study

The aim of this work was to demonstrate at pilot scale a high level of energy recovery from sewage utilising a primary Anaerobic Migrating Bed Reactor (AMBR) operating at ambient temperature to convert COD to methane. The focus is the reduction in non-renewable CO2 emissions resulting from reduced energy requirements for sewage treatment. A pilot AMBR was operated on screened sewage over the period June 2003 to September 2004. The study was divided into two experimental phases. In Phase 1 the process operated at a feed rate of 10 L/h (HRT 50 h), SRT 63 days, average temperature 28 degrees C and mixing time fraction 0.05. In Phase 2 the operating parameters were 20 L/h, 26 days, 16 degrees C and 0.025. Methane production was 66% of total sewage COD in Phase 1 and 23% in Phase 2. Gas mixing of the reactor provided micro-aeration which suppressed sulphide production. Intermittent gas mixing at a useful power input of 6 W/m(3) provided satisfactory process performance in both phases. Energy consumption for mixing was about 1.5% of the energy conversion to methane in both operating phases. Comparative analysis with previously published data confirmed that methane supersaturation resulted in significant losses of methane in the effluent of anaerobic treatment systems. No cases have been reported where methane was considered to be supersaturated in the effluent. We have shown that methane supersaturation is likely to be significant and that methane losses in the effluent are likely to have been greater than previously predicted. Dissolved methane concentrations were measured at up to 2.2 times the saturation concentration relative to the mixing gas composition. However, this study has also demonstrated that despite methane supersaturation occurring, microaeration can result in significantly lower losses of methane in the effluent (< 11% in this study), and has demonstrated that anaerobic sewage treatment can genuinely provide energy recovery. The goal of demonstrating a high level of energy recovery in an ambient anaerobic bioreactor was achieved. An AMBR operating at ambient temperature can achieve up to 70% conversion of sewage COD to methane, depending on SRT and temperature. (c) 2006 Wiley Periodicals, Inc.

Eigenvalue sensitivity analysis for probabilistic small signal stability assessment

This paper presents a method to analyze the first order eigenvalue sensitivity with respect to the operating parameters of a power system. The method is based on explicitly expressing the system state matrix into sub-matrices. The eigenvalue sensitivity is calculated based on the explicitly formed system state matrix. The 4th order generator model and 4th order exciter system model are used to form the system state matrix. A case study using New England 10-machine 39-bus system is provided to demonstrate the effectiveness of the proposed method. This method can be applied into large scale power system eigenvalue sensitivity with respect to operating parameters.

Scaleup of wet granulation processes: science not art

Significant advances have been made in the last decade to quantify the process of wet granulation. The attributes of product granules from the granulation process are controlled by a combination of three groups of processes occurring in the granulator: (1) wetting and nucleation, (2) growth and consolidation and (3) breakage and attrition. For the first two of these processes, the key controlling dimensionless groups are defined and regime maps are presented and validated with data from tumbling and mixer granulators. Granulation is an example of particle design. For quantitative analysis, both careful characterisation of the feed formulation and knowledge of operating parameters are required. A key thesis of this paper is that the design, scaleup and operation of granulation processes can now be considered as quantitative engineering rather than a black art. Résumé

Coalescence of secondary dispersions in fibrous beds

SUMMARY A study has been made of the coalescence of secondary dispersions in a fibrous bed. The literature pertaining to the formation, hydrodynamic behaviour and methods of separation of droplets less than one hundred micrometres in diameter has been reviewed with particular reference to fibrous bed coalescers. The main operating parameters were identified as inlet drop size distribution, phase ratio, superficial velocity and the thickness and voidage of the bed . A recirculatory rig with interchangeable fibrous bed pads was designed and operated with toluene-water dispersions generated by a combination of centrifugal pumps . Inlet drop sizes were analysed using a Coulter Counter and outlet drops were sized photographically. A novel technique, involving conductivity measur ements at different planes in the bed, was developed to measure hold up distribution. Single phase flow and two phase flow pressure drops were correlated by a Blake-Kozeny type equation. Exit drop size was independent of inlet drop size distribution and phase ratio but a function of superficialvelocity and packing thickness. Average bed hold up was independent of inlet drop size distribution and phase ratio, but decreased with increase in superficial velocity. Hold up was not evenly distributed in the bed, the highest value occurred at the inlet followed by a sharp -2 drop at approximately 1.2 x 10 m. Hold up remained constant throughout the rest of the bed until the exit plane, where it increased. From the results, a mechanism is postulated involving: (a) Capture of the inlet drops followed by interdrop coalescence until an equilibrium value is reached. (b) Equilibrium size droplets flowing as rivulets through the intermediate portion of the bed, and (c) Each rivulet forms droplets at the exit face, which detach by a 'drip point' mechanism.

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.

Investigation of chaotic mixing in laminar fluid systems by the use of computational fluid dynamics

This work presents significant development into chaotic mixing induced through periodic boundaries and twisting flows. Three-dimensional closed and throughput domains are shown to exhibit chaotic motion under both time periodic and time independent boundary motions, A property is developed originating from a signature of chaos, sensitive dependence to initial conditions, which successfully quantifies the degree of disorder withjn the mixing systems presented and enables comparisons of the disorder throughout ranges of operating parameters, This work omits physical experimental results but presents significant computational investigation into chaotic systems using commercial computational fluid dynamics techniques. Physical experiments with chaotic mixing systems are, by their very nature, difficult to extract information beyond the recognition that disorder does, does not of partially occurs. The initial aim of this work is to observe whether it is possible to accurately simulate previously published physical experimental results through using commercial CFD techniques. This is shown to be possible for simple two-dimensional systems with time periodic wall movements. From this, and subsequent macro and microscopic observations of flow regimes, a simple explanation is developed for how boundary operating parameters affect the system disorder. Consider the classic two-dimensional rectangular cavity with time periodic velocity of the upper and lower walls, causing two opposing streamline motions. The degree of disorder within the system is related to the magnitude of displacement of individual particles within these opposing streamlines. The rationale is then employed in this work to develop and investigate more complex three-dimensional mixing systems that exhibit throughputs and time independence and are therefore more realistic and a significant advance towards designing chaotic mixers for process industries. Domains inducing chaotic motion through twisting flows are also briefly considered. This work concludes by offering possible advancements to the property developed to quantify disorder and suggestions of domains and associated boundary conditions that are expected to produce chaotic mixing.

Combined bioreaction and separation in a simulated counter-current chromatographic bioreactor-separator system

The objective of this work has been to investigate the principle of combined bioreaction and separation in a simulated counter-current chromatographic bioreactor-separator system (SCCR-S). The SCCR-S system consisted of twelve 5.4cm i.d x 75cm long columns packed with calcium charged cross-linked polystyrene resin. Three bioreactions, namely the saccharification of modified starch to maltose and dextrin using the enzyme maltogenase, the hydrolysis of lactose to galactose and glucose in the presence of the enzyme lactase and the biosynthesis of dextran from sucrose using the enzyme dextransucrase. Combined bioreaction and separation has been successfully carried out in the SCCR-S system for the saccharification of modified starch to maltose and dextrin. The effects of the operating parameters (switch time, eluent flowrate, feed concentration and enzyme activity) on the performance of the SCCR-S system were investigated. By using an eluent of dilute enzyme solution, starch conversions of up to 60% were achieved using lower amounts of enzyme than the theoretical amount required by a conventional bioreactor to produce the same amount of maltose over the same time period. Comparing the SCCR-S system to a continuous annular chromatograph (CRAC) for the saccharification of modified starch showed that the SCCR-S system required only 34.6-47.3% of the amount of enzyme required by the CRAC. The SCCR-S system was operated in the batch and continuous modes as a bioreactor-separator for the hydrolysis of lactose to galactose and glucose. By operating the system in the continuous mode, the operating parameters were further investigated. During these experiments the eluent was deionised water and the enzyme was introduced into the system through the same port as the feed. The galactose produced was retarded and moved with the stationary phase to be purge as the galactose rich product (GalRP) while the glucose moved with the mobile phase and was collected as the glucose rich product (GRP). By operating at up to 30%w/v lactose feed concentrations, complete conversions were achieved using only 48% of the theoretical amount of enzyme required by a conventional bioreactor to hydrolyse the same amount of glucose over the same time period. The main operating parameters affecting the performance of the SCCR-S system operating in the batch mode were investigated and the results compared to those of the continuous operation of the SCCR-S system. . During the biosynthesis of dextran in the SCCR-S system, a method of on-line regeneration of the resin was required to operate the system continuously. Complete conversion was achieved at sucrose feed concentrations of 5%w/v with fructose rich. products (FRP) of up to 100% obtained. The dextran rich products were contaninated by small amounts of glucose and levan formed during the bioreaction. Mathematical modelling and computer simulation of the SCCR-S. system operating in the continuous mode for the hydrolysis of lactose has been carried out. .

The atomisation of detergent slurries

A study was made to determine the conditions under which the optimum droplet size distribution (ie., narrowest size range with a minimum of fines and over-sized agglomerates), is generated in sprays from centrifugal pressure nozzles. A range of non-Newtonian detergent slurries were tested but the results are of wider application and parallel work was undertaken with water, ionic solutions and chalk slurries. Six centrifugal pressure nozzles were used and the drop-size distributions correlated as a function of fluid properties, pressure, fiowrate, feed temperature, and nozzle characteristics. Measurements were made using a Malvern Particle Size Anayser slung across a specially-designed transparent tower section of approximately 1.2m diameter in order to reduce obscuration caused by the spray and improve existing droplet sizing techniques. The results obtained were based upon the Rosin-Rammler distribution model and the Size Analyser provided a direct print-out of size distribution and the parameters characterising it. A Spraying System nozzle, AAASSTC8-8, produced the optimum spray distribution with the detergent slurry at a temperature of 60°C whilst operating at 1200 psi. With other fluids the Delevan 2.2SJ nozzle produced the optimum spray distribution operating at 1200 psi but with the Spraying Systems nozzles there was no clear-cut optimum set of conditions, ie. the nozzle and pressure varied depending upon the fluid being sprayed. The mechanisms of liquid sheet break-up and droplet dispersion were investigated in specially-constructed, scaled-up, transparent nozzles. A mathematical model of centrifugal pressure nozzle atomisation was developed based upon fundamental operating parameters rather than resorting to empirical correlations. This enabled theoretical predictions to be made over a wide range of operating conditions and nozzle types. The model predictions for volumetric fiowrate, liquid sheet length and air core diameter showed good agreement with the experimentally determined results. However, the model predicted smaller droplet sizes than were produced experimentally due to inaccuracies identified in the initial assumptions.

The separation of secondary oil-water dispersions in particulate beds

The literature relating to haze formation, methods of separation, coalescence mechanisms, and models by which droplets <100 μm are collected, coalesced and transferred, have been reviewed with particular reference to particulate bed coalescers. The separation of secondary oil-water dispersions was studied experimentally using packed beds of monosized glass ballotini particles. The variables investigated were superficial velocity, bed depth, particle size, and the phase ratio and drop size distribution of inlet secondary dispersion. A modified pump loop was used to generate secondary dispersions of toluene or Clairsol 350 in water with phase ratios between 0.5-6.0 v/v%.Inlet drop size distributions were determined using a Malvern Particle Size Analyser;effluent, coalesced droplets were sized by photography. Single phase flow pressure drop data were correlated by means of a Carman-Kozeny type equation. Correlations were obtained relating single and two phase pressure drops, as (ΔP2/μc)/ΔP1/μd) = kp Ua Lb dcc dpd Cine A flow equation was derived to correlate the two phase pressure drop data as, ΔP2/(ρcU2) = 8.64*107 [dc/D]-0.27 [L/D]0.71 [dp/D]-0.17 [NRe]1.5 [e1]-0.14 [Cin]0.26  In a comparison between functions to characterise the inlet drop size distributions a modification of the Weibull function provided the best fit of experimental data. The general mean drop diameter was correlated by: q_p q_p p_q /β      Γ ((q-3/β) +1) d qp = d fr  .α        Γ ((P-3/β +1 The measured and predicted mean inlet drop diameters agreed within ±15%. Secondary dispersion separation depends largely upon drop capture within a bed. A theoretical analysis of drop capture mechanisms in this work indicated that indirect interception and London-van der Waal's mechanisms predominate. Mathematical models of dispersed phase concentration m the bed were developed by considering drop motion to be analogous to molecular diffusion.The number of possible channels in a bed was predicted from a model in which the pores comprised randomly-interconnected passage-ways between adjacent packing elements and axial flow occured in cylinders on an equilateral triangular pitch. An expression was derived for length of service channels in a queuing system leading to the prediction of filter coefficients. The insight provided into the mechanisms of drop collection and travel, and the correlations of operating parameters, should assist design of industrial particulate bed coalescers.

Improvements to the arc-spraying process

After a brief review of the various forms of thermal spraying equipment and processes, descriptions of the basic principles involved and the general functions for which thermally sprayed coatings are used are given. The background of the collaborating company, Metallisation, is described and their position in the overall market discussed, providing a backdrop against which the appropriateness of various project options might be judged. Current arc-spraying equipment is then examined, firstly in terms of the workings of their constituent parts and subsequently by examining the effects of changes in design and in operating parameters both upon equipment operation and the coatings produced. Published literature relating to these matters is reviewed. Literature relating to the production, comminution and propulsion of the particles which form the spray is discussed as are the mechanisms involved at impact with the substrate. Literature on the use of rockets for thermal spraying and induction heating as a process for feedstock melting are also reviewed. Three distinct options for further study are derived and preliminary tests and costings made to allow one option alone, the use of rocket acceleration, to go forward to the experimental phase. A suitable rocket burner was developed, tested and incorporated into an arc-spray system so that the sprayability of the whole could be assessed. Coatings were made using various parameters and these are compared with coatings produced by a standard system. Coatings were examined for macro and micro hardness, cohesive strength, porosity and by microstructural examination. The results indicate a high degree of similarity between the coatings produced by the standard system and the high velocity system. This was surprising in view of the very different atomising media and velocities. Possible causes for this similarity and the general behaviour of this new system and the standard system are discussed before the study reaches its conclusions in not proving the hypothesis that an increase in particle velocity would improve the mechanical properties of arc-sprayed steel coatings. KEY WORDS: Sprayed metal coatings, Electric arc spraying, High velocity flame spraying, Sprayed coating properties

Effects of engine cooling water temperature on performance and exhaust emission characteristics of a multi-cylinder CI engine operated with Jatropha-Butanol blend

Renewable alternatives such as biofuels and optimisation of the engine operating parameters can enhance engine performance and reduce emissions. The temperature of the engine coolant is known to have significant influence on engine performance and emissions. Whereas much existing literature describes the effects of coolant temperature in engines using fossil derived fuels, very few studies have investigated these effects when biofuel is used as an alternative fuel. Jatropha oil is a non-edible biofuel which can substitute fossil diesel for compression ignition (CI) engine use. However, due to the high viscosity of Jatropha oil, technique such as transesterification, preheating the oil, mixing with other fuel is recommended for improved combustion and reduced emissions. In this study, Jatropha oil was blended separately with ethanol and butanol, at ratios of 80:20 and 70:30. The fuel properties of all four blends were measured and compared with diesel and jatropha oil. It was found that the 80% jatropha oil + 20% butanol blend was the most suitable alternative, as its properties were closest to that of diesel. A 2 cylinder Yanmar engine was used; the cooling water temperature was varied between 50°C and 95°C. In general, it was found that when the temperature of the cooling water was increased, the combustion process enhanced for both diesel and Jatropha-Butanol blend. The CO2 emissions for both diesel and biofuel blend were observed to increase with temperature. As a result CO, O2 and lambda values were observed to decrease when cooling water temperature increased. When the engine was operated using diesel, NOX emissions correlated in an opposite manner to smoke opacity; however, when the biofuel blend was used, NOX emissions and smoke opacity correlated in an identical manner. The brake thermal efficiencies were found to increase slightly as the temperature was increased. In contrast, for all fuels, the volumetric efficiency was observed to decrease as the coolant temperature was increased. Brake specific fuel consumption was observed to decrease as the temperature was increased and was higher on average when the biofuel was used, in comparison to diesel. The study concludes that the effects of engine coolant temperature on engine performance and emission characteristics differ between biofuel blend and fossil diesel operation. The coolant temperature needs to be optimised depending on the type of biofuel for optimum engine performance and reduced emissions.

High rate growth of nanocrystalline diamond films using high microwave power and pure nitrogen/methane/hydrogen plasma

In this work, we investigate the impact of minute amounts of pure nitrogen addition into conventional methane/hydrogen mixtures on the growth characteristics of nanocrystalline diamond (NCD) films by microwave plasma assisted chemical vapour deposition (MPCVD), under high power conditions. The NCD films were produced from a gas mixture of 4% CH4/H2 with two different concentrations of N2 additive and microwave power ranging from 3.0 kW to 4.0 kW, while keeping all the other operating parameters constant. The morphology, grain size, microstructure and texture of the resulting NCD films were characterized by using scanning electron microscope (SEM), micro-Raman spectroscopy and X-ray diffraction (XRD) techniques. N2 addition was found to be the main parameter responsible for the formation and for the key change in the growth characteristics of NCD films under the employed conditions. Growth rates ranging from 5.4 μm/h up to 9.6 μm/h were achieved for the NCD films, much higher than those usually reported in the literature. The enhancing factor of nitrogen addition on NCD growth rate was obtained by comparing with the growth rate of large-grained microcrystalline diamond films grown without nitrogen and discussed by comparing with that of single crystal diamond through theoretical work in the literature. This achievement on NCD growth rate makes the technology interesting for industrial applications where fast coating of large substrates is highly desirable.

Rearrangement of α-pinene oxide using a surface catalysed spinning disc reactor (SDR)

Isomerisation of α-pinene oxide to campholenic aldehyde was performed by immobilising zinc triflate based catalysts on the surface of a spinning disc reactor (SDR). Two types of catalyst have been studied and the influence of operating parameters such as rotational speed, feed flow rate and reaction temperature on conversion and selectivity towards campholenic aldehyde has been investigated in considerable detail. The findings of the study suggest that immobilising the catalyst on the reactor surface and performing the reaction in continuous mode has potential for achieving benefits of Green Chemical Technology (GCT).

Concentration polarization model of spiral-wound membrane modules with application to batch-mode RO desalination of brackish water

Batch-mode reverse osmosis (batch-RO) operation is considered a promising desalination method due to its low energy requirement compared to other RO system arrangements. To improve and predict batch-RO performance, studies on concentration polarization (CP) are carried out. The Kimura-Sourirajan mass-transfer model is applied and validated by experimentation with two different spiral-wound RO elements. Explicit analytical Sherwood correlations are derived based on experimental results. For batch-RO operation, a new genetic algorithm method is developed to estimate the Sherwood correlation parameters, taking into account the effects of variation in operating parameters. Analytical procedures are presented, then the mass transfer coefficient models are developed for different operation processes, i.e., batch-RO and continuous RO. The CP related energy loss in batch-RO operation is quantified based on the resulting relationship between feed flow rates and mass transfer coefficients. It is found that CP increases energy consumption in batch-RO by about 25% compared to the ideal case in which CP is absent. For continuous RO process, the derived Sherwood correlation predicted CP accurately. In addition, we determined the optimum feed flow rate of our batch-RO system.