Accuracy and optimal sampling in Monte Carlo solution of population balance equations

Implementation of a Monte Carlo simulation for the solution of population balance equations (PBEs) requires choice of initial sample number (N0), number of replicates (M), and number of bins for probability distribution reconstruction (n). It is found that Squared Hellinger Distance, H2, is a useful measurement of the accuracy of Monte Carlo (MC) simulation, and can be related directly to N0, M, and n. Asymptotic approximations of H2 are deduced and tested for both one-dimensional (1-D) and 2-D PBEs with coalescence. The central processing unit (CPU) cost, C, is found in a power-law relationship, C= aMNb0, with the CPU cost index, b, indicating the weighting of N0 in the total CPU cost. n must be chosen to balance accuracy and resolution. For fixed n, M × N0 determines the accuracy of MC prediction; if b > 1, then the optimal solution strategy uses multiple replications and small sample size. Conversely, if 0 < b < 1, one replicate and a large initial sample size is preferred. © 2015 American Institute of Chemical Engineers AIChE J, 61: 2394–2402, 2015

Modeling and simulation of population balances for particulate processes

This work focuses on the modeling and numerical approximations of population balance equations (PBEs) for the simulation of different phenomena occurring in process engineering. The population balance equation (PBE) is considered to be a statement of continuity. It tracks the change in particle size distribution as particles are born, die, grow or leave a given control volume. In the population balance models the one independent variable represents the time, the other(s) are property coordinate(s), e.g., the particle volume (size) in the present case. They typically describe the temporal evolution of the number density functions and have been used to model various processes such as granulation, crystallization, polymerization, emulsion and cell dynamics. The semi-discrete high resolution schemes are proposed for solving PBEs modeling one and two-dimensional batch crystallization models. The schemes are discrete in property coordinates but continuous in time. The resulting ordinary differential equations can be solved by any standard ODE solver. To improve the numerical accuracy of the schemes a moving mesh technique is introduced in both one and two-dimensional cases ...

Implementation of population balances in simulation of gas-liquid reactors using a CFD-PBM coupled approach

In bubbly flow simulations, bubble size distribution is an important factor in determination of hydrodynamics. Beside hydrodynamics, it is crucial in the prediction of interfacial area available for mass transfer and in the prediction of reaction rate in gas-liquid reactors such as bubble columns. Solution of population balance equations is a method which can help to model the size distribution by considering continuous bubble coalescence and breakage. Therefore, in Computational Fluid Dynamic simulations it is necessary to couple CFD and Population Balance Model (CFD-PBM) to get reliable distribution. In the current work a CFD-PBM coupled model is implemented as FORTRAN subroutines in ANSYS CFX 10 and it has been tested for bubbly flow. This model uses the idea of Multi Phase Multi Size Group approach which was previously presented by Sha et al. (2006) [18]. The current CFD-PBM coupled method considers inhomogeneous flow field for different bubble size groups in the Eulerian multi-dispersed phase systems. Considering different velocity field for bubbles can give the advantageof more accurate solution of hydrodynamics. It is also an improved method for prediction of bubble size distribution in multiphase flow compared to available commercial packages.

Balance del modelo agroexportador en América Latina al comenzar el siglo XXI

El propósito de este escrito es ofrecer un panorama del momento por el que atraviesa hoy la agricultura regional, revisando algunos de sus saldos agregados. Nos interesa debatir el tono optimista y acrítico con que se divulgan las cifras y poner en evidencia aspectos menos difundidos, pero de suma importancia para el presente y futuro de la región. Para comenzar nos referimos al contexto general en que se suscitan las transformaciones en el modelo de crecimiento de la región y sus pretensiones. Enseguida nos ocupamos de la agricultura, enfocados en el contenido común de los proyectos nacionales, para revisar luego algunas de sus secuelas. Concluimos con algunas reflexiones sobre el significado de estos cambios para la expansión del capital en tanto sistema y para la población rural involucrada.

Balance del modelo agroexportador en América Latina al comenzar el siglo XXI

El propósito de este escrito es ofrecer un panorama del momento por el que atraviesa hoy la agricultura regional, revisando algunos de sus saldos agregados. Nos interesa debatir el tono optimista y acrítico con que se divulgan las cifras y poner en evidencia aspectos menos difundidos, pero de suma importancia para el presente y futuro de la región. Para comenzar nos referimos al contexto general en que se suscitan las transformaciones en el modelo de crecimiento de la región y sus pretensiones. Enseguida nos ocupamos de la agricultura, enfocados en el contenido común de los proyectos nacionales, para revisar luego algunas de sus secuelas. Concluimos con algunas reflexiones sobre el significado de estos cambios para la expansión del capital en tanto sistema y para la población rural involucrada.

Balance del modelo agroexportador en América Latina al comenzar el siglo XXI

El propósito de este escrito es ofrecer un panorama del momento por el que atraviesa hoy la agricultura regional, revisando algunos de sus saldos agregados. Nos interesa debatir el tono optimista y acrítico con que se divulgan las cifras y poner en evidencia aspectos menos difundidos, pero de suma importancia para el presente y futuro de la región. Para comenzar nos referimos al contexto general en que se suscitan las transformaciones en el modelo de crecimiento de la región y sus pretensiones. Enseguida nos ocupamos de la agricultura, enfocados en el contenido común de los proyectos nacionales, para revisar luego algunas de sus secuelas. Concluimos con algunas reflexiones sobre el significado de estos cambios para la expansión del capital en tanto sistema y para la población rural involucrada.

Determination of coalescence kernels for high-shear granulation using DEM simulations

Discrete element method (DEM) modeling is used in parallel with a model for coalescence of deformable surface wet granules. This produces a method capable of predicting both collision rates and coalescence efficiencies for use in derivation of an overall coalescence kernel. These coalescence kernels can then be used in computationally efficient meso-scale models such as population balance equation (PBE) models. A soft-sphere DEM model using periodic boundary conditions and a unique boxing scheme was utilized to simulate particle flow inside a high-shear mixer. Analysis of the simulation results provided collision frequency, aggregation frequency, kinetic energy, coalescence efficiency and compaction rates for the granulation process. This information can be used to bridge the gap in multi-scale modeling of granulation processes between the micro-scale DEM/coalescence modeling approach and a meso-scale PBE modeling approach.

Physical function interfering with pain and symptoms in fibromyalgia patients

Objectives. The aim of this study was to assess the relationship between variables of physical assessment - muscular strength, flexibility and dynamic balance - with pain, pain threshold, and fibromyalgia symptoms (FM). Methods. Our sample consists of 55 women, with age ranging from 30 to 55 years (mean of 46.5, (standard deviation, SD=6.6)), mean body mass index (BMI) of 28.7(3.8) and diagnosed for FM according to the American College of Rheumatology criteria. Pain intensity was measured using a visual analogue scale (VAS) and pain threshold (PT) using Fisher`s dolorimeter. FM symptoms were assessed by the Fibromyalgia Impact Questionnaire (FIQ); flexibility by the third finger to floor test (3FF); the muscular strength index (MSI) by the maximum volunteer isometric contraction at flexion and extension of right knee and elbow using a force transducer, dynamic balance by the time to get up and go (TUG) test and the functional reach test (FRT). Data were analysed using Pearson`s correlation, as well as simple and multivariate regression tests, with significance level of 5%. Results. PT and FIQ were weakly but significantly correlated with the TUG, MSI and 3FF as well as VAS with the TUG and MSI (p<0.05). VAS, PT and FIQ was not correlated with FRT. Simple regression suggests that, alone, TUG, FR, MSI and 3FF are low predictors of VAS, PT and FIQ. For the VAS, the best predictive model includes TUG and MSI, explaining 12.6% of pain. variability. For TP and total symptoms, as obtained by the FIQ, most predictive model includes 3FF and MSI, which respectively respond by 30% and 21% of the variability. Conclusion. Muscular strength, flexibility and balance are associated with pain, pain threshold, and symptoms in FM patients.

Polymer product engineering utilising oscillatory baffled reactors

The Oscillatory baffled reactor (OBR) can be used to produce particles with controlled size and morphology, in batch or continuous flow. This is due to the effect of the superimposed oscillations that radially mixes fluid but still allows plug-flow (or close to plug flow) behaviour in a continuous system. This mixing, combined with a close to a constant level of turbulence intensity in the reactor, leads to tight droplet and subsequent product particle size distributions. By applying population balance equations together with experimental droplet size distributions, breakage rates of droplets can be determined and this is a useful tool for understanding the product engineering in OBRs. (C) 2002 Elsevier Science B.V All rights reserved.

Review and future directions in the modelling and control of continuous drum granulation

Many granulation plants operate well below design capacity, suffering from high recycle rates and even periodic instabilities. This behaviour cannot be fully predicted using the present models. The main objective of the paper is to provide an overview of the current status of model development for granulation processes and suggest future directions for research and development. The end-use of the models is focused on the optimal design and control of granulation plants using the improved predictions of process dynamics. The development of novel models involving mechanistically based structural switching methods is proposed in the paper. A number of guidelines are proposed for the selection of control relevant model structures. (C) 2002 Published by Elsevier Science B.V.

Control and Simulation of Batch Crystallization

Concerning process control of batch cooling crystallization the present work focused on the cooling profile and seeding technique. Secondly, the influence of additives on batch-wise precipitation process was investigated. Moreover, a Computational Fluid Dynamics (CFD) model for simulation of controlled batch cooling crystallization was developed. A novel cooling model to control supersaturation level during batch-wise cooling crystallization was introduced. The crystallization kinetics together with operating conditions, i.e. seed loading, cooling rate and batch time, were taken into account in the model. Especially, the supersaturation- and suspension density- dependent secondary nucleation was included in the model. The interaction between the operating conditions and their influence on the control target, i.e. the constant level of supersaturation, were studied with the aid of a numerical solution for the cooling model. Further, the batch cooling crystallization was simulated with the ideal mixing model and CFD model. The moment transformation of the population balance, together with the mass and heat balances, were solved numerically in the simulation. In order to clarify a relationship betweenthe operating conditions and product sizes, a system chart was developed for anideal mixing condition. The utilization of the system chart to determine the appropriate operating condition to meet a required product size was introduced. With CFD simulation, batch crystallization, operated following a specified coolingmode, was studied in the crystallizers having different geometries and scales. The introduced cooling model and simulation results were verified experimentallyfor potassium dihydrogen phosphate (KDP) and the novelties of the proposed control policies were demonstrated using potassium sulfate by comparing with the published results in the literature. The study on the batch-wise precipitation showed that immiscible additives could promote the agglomeration of a derivative of benzoic acid, which facilitated the filterability of the crystal product.

Gas-liquid mass transfer in bubbly flow: Estimation of mass transfer, bubble size and reactor performance in various applications

Gas-liquid mass transfer is an important issue in the design and operation of many chemical unit operations. Despite its importance, the evaluation of gas-liquid mass transfer is not straightforward due to the complex nature of the phenomena involved. In this thesis gas-liquid mass transfer was evaluated in three different gas-liquid reactors in a traditional way by measuring the volumetric mass transfer coefficient (kLa). The studied reactors were a bubble column with a T-junction two-phase nozzle for gas dispersion, an industrial scale bubble column reactor for the oxidation of tetrahydroanthrahydroquinone and a concurrent downflow structured bed.The main drawback of this approach is that the obtained correlations give only the average volumetric mass transfer coefficient, which is dependent on average conditions. Moreover, the obtained correlations are valid only for the studied geometry and for the chemical system used in the measurements. In principle, a more fundamental approach is to estimate the interfacial area available for mass transfer from bubble size distributions obtained by solution of population balance equations. This approach has been used in this thesis by developing a population balance model for a bubble column together with phenomenological models for bubble breakage and coalescence. The parameters of the bubble breakage rate and coalescence rate models were estimated by comparing the measured and calculated bubble sizes. The coalescence models always have at least one experimental parameter. This is because the bubble coalescence depends on liquid composition in a way which is difficult to evaluate using known physical properties. The coalescence properties of some model solutions were evaluated by measuring the time that a bubble rests at the free liquid-gas interface before coalescing (the so-calledpersistence time or rest time). The measured persistence times range from 10 msup to 15 s depending on the solution. The coalescence was never found to be instantaneous. The bubble oscillates up and down at the interface at least a coupleof times before coalescence takes place. The measured persistence times were compared to coalescence times obtained by parameter fitting using measured bubble size distributions in a bubble column and a bubble column population balance model. For short persistence times, the persistence and coalescence times are in good agreement. For longer persistence times, however, the persistence times are at least an order of magnitude longer than the corresponding coalescence times from parameter fitting. This discrepancy may be attributed to the uncertainties concerning the estimation of energy dissipation rates, collision rates and mechanisms and contact times of the bubbles.

Ultrafine grinding of FGD and phosphogypsum with an attrition bead mill and a jet mill : optimisation and modelling of grinding and mill comparison

Fine powders of minerals are used commonly in the paper and paint industry, and for ceramics. Research for utilizing of different waste materials in these applications is environmentally important. In this work, the ultrafine grinding of two waste gypsum materials, namely FGD (Flue Gas Desulphurisation) gypsum and phosphogypsum from a phosphoric acid plant, with the attrition bead mill and with the jet mill has been studied. The ' objective of this research was to test the suitability of the attrition bead mill and of the jet mill to produce gypsum powders with a particle size of a few microns. The grinding conditions were optimised by studying the influences of different operational grinding parameters on the grinding rate and on the energy consumption of the process in order to achieve a product fineness such as that required in the paper industry with as low energy consumption as possible. Based on experimental results, the most influential parameters in the attrition grinding were found to be the bead size, the stirrer type, and the stirring speed. The best conditions, based on the product fineness and specific energy consumption of grinding, for the attrition grinding process is to grind the material with small grinding beads and a high rotational speed of the stirrer. Also, by using some suitable grinding additive, a finer product is achieved with a lower energy consumption. In jet mill grinding the most influential parameters were the feed rate, the volumetric flow rate of the grinding air, and the height of the internal classification tube. The optimised condition for the jet is to grind with a small feed rate and with a large rate of volumetric flow rate of grinding air when the inside tube is low. The finer product with a larger rate of production was achieved with the attrition bead mill than with the jet mill, thus the attrition grinding is better for the ultrafine grinding of gypsum than the jet grinding. Finally the suitability of the population balance model for simulation of grinding processes has been studied with different S , B , and C functions. A new S function for the modelling of an attrition mill and a new C function for the modelling of a jet mill were developed. The suitability of the selected models with the developed grinding functions was tested by curve fitting the particle size distributions of the grinding products and then comparing the fitted size distributions to the measured particle sizes. According to the simulation results, the models are suitable for the estimation and simulation of the studied grinding processes.

Does the Ingestion of Isoxaflutole Herbicide Affect the Midgut and Salivary Glands of Pentatomidae Predators?

Podisus nigrispinus is a generalist predator naturally occurring in agricultural and forestry systems that effectively contributes to the population balance of phytophagous insects, especially defoliating caterpillars. Histological changes were evaluated in the salivary glands and midgut of P. nigrispinus caused by ingestion of systemic herbicide isoxaflutole. These predator females were fed with leaves of eucalyptus plants, Tenebrio molitor pupae or water, contaminated or not by herbicide. Salivary glands and midguts were dissected, processed and analyzed under a light microscope. Activity level and cell morphology of the salivary glands and midgut showed differences among insects fed on plants, contaminated water or pupae. The epithelia of the salivary gland and midgut of individuals which had no contact with the herbicide showed homogeneous cytoplasm, nucleus with predominance of decondensed chromatin and evident nucleoli, intense cell activity features. As for the insects in contact with contaminated food, they presented undeveloped nucleus and condensed chromatin. The luminal contents of the salivary glands in the contaminated insects had become more acidophilus than in insects without poisoning, as well as having heterogeneous and granular secretion, being more evident in the bioassay in which the insects fed on contaminated water. There was a marked morphological change in the midgut cells in contaminated insects. High degree of apoptosis, disorganization and secretory vacuoles in the epithelial cytoplasm were observed. The apical portion of the midgut cells proved undeveloped, irregular and partially destroyed. It is concluded that isoxaflutole causes morphological changes in the digestive system of the predator P. nigrispinus.