CFD modeling of the Wurster bed coater

In the Wurster bed coater, the wetting, drying, and circulation of particles are combined to produce a high quality coating. The drying and wetting conditions in a laboratory scale Wurster bed coater are modeled and compared with experimental data. A model combining multiphase fluid dynamics with heat and mass transfer is developed to model the particle and gas motion and the transport of thermal energy and moisture. A wetting region is defined, where a specified moisture content is set in the particle phase, above the jet inlet, to describe the injection of coating liquid. The simulation shows the characteristic circulation of particles in the equipment, as well as the behavior of the moisture in the system and agrees well with measurements. The simulation indicates how different process conditions influence the drying regions. The results show that most of the drying, under typical operating conditions, takes place in the Wurster tube. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

Measurement of the particle movement in the fountain region of a Wurster type bed

A non-intrusive technique to follow single particles in the fountain region of a spouted bed with draft tube, the Wurster coater, is developed and tested. Compared to other techniques found in literature, the measurement technique presented here is easy to handle, inexpensive and suitable for particles with a diameter down to 500 μm. The measurements are performed in a laboratory scale Wurster bed constructed of Plexiglas to allow optical access to the flow. A small quantity of the particles in the bed are marked with fluorescence. A UV lamp is used to excite the fluorescence and the marked particles are followed with a high-speed video camera equipped with an optical filter. A single marked particle in the fountain can be followed in 5-30 images in sequence at a frame rate of 125 fps. The particle position and velocity are calculated. The distribution of the particle trajectories maximum height is studied for different particle loading, jet air velocity and position of the Wurster tube. The technique can be used for the characterization of differences in process dynamics due to different operating conditions.     

CFD study of droplet atomisation using a binary nozzle in fluidised bed coating

An Eulerian–Lagrangian computational fluid dynamics (CFD) model was built to describe two-fluid atomisation in a tapered fluidised bed coater using the air-blast/air-assisted atomiser model. Atomisation was modelled both with and without the inclusion of the solid phase (i.e. gas–liquid and gas–solid–liquid multiphase modelling). In addition, a multi-fluid flow model (Eulerian–Eulerian framework) combined with a population balance model was used as an alternative approach for modelling the spray produced by a two-fluid nozzle. In this approach, the CFD solver couples the population balance equation along with the Navier–Stokes equations for predicting the droplet diameter and mass fraction distribution. Comparison between simulated spray pattern (gas–liquid model) and that experimentally visualised by means of UV illumination was made and a good agreement was obtained. Parametric studies were done in order to investigate the effects of operating conditions on spray cone and liquid mass fraction inside the reactor. Furthermore, comparison of time-averaged fluidised bed behaviour with the inclusion of sprays obtained by both gas–solid–liquid multiphase modelling methods is presented.     

固定床反应器内气体预分布器研究

研究了直径1 000 mm,高3 000 mm的固定床冷模装置中气体预分布器对反应器内气流分布的影响。结果表明:气体分布器可改变床层内气流流形并使径向气流的速度分布趋于均匀;随着表观气速的增加,反应器内气流的不均匀程度增加;分布器的环隙高度在一定的范围可使反应器内气流的不均匀程度相对较好。应用计算流体力学软件CFX对固定床反应器内的流场进行模拟计算,并与大型冷模试验测试结果进行比较,模型计算值和冷模试验测量值吻合较好。     

Powder coating efficiency of small particles and their agglomeration in circulating fluidized bed

The coating efficiency of fluidizing small particles and their agglomeration were investigated to evaluate the possibility of powder coating by the use of a circulating fluidized bed. Glass beads, whose mean diameter was 43 Μm, and silica powder of 1 Μm were used as a core and a coating material. Polyvinyl alcohol was used as a binder and its solution was supplied together with silica powder from a spray nozzle equipped in the circulating fluidized bed. Glass beads of 43 Μm, which had been impossible to coat in a conventional fluidized bed coater, were successfully coated with silica powder in a circulating fluidized bed, and agglomeration among core particles was prevented. From this result, it was confirmed that a circulating fluidized bed performs excellently as a coater, especially for fine core particles, so a circulating fluidized bed coater has bright prospects for particle coating.     

Local analysis of the tablet coating process: Impact of operation conditions on film quality

Spray coating is frequently used in the pharmaceutical industry to control the release of the active pharmaceutical ingredient of a tablet or to mask its taste. The uniformity of the coating is of significant importance, as the coating usually has critical functional properties. However, coating uniformity is difficult to predict without significant experimental work, and even advanced particle simulations need to be augmented by CFD models to fully describe the coating uniformity on a single tablet.In this study we analyze the coating process by using detailed computational fluid dynamics (CFD) multiphase spray simulations. The impact and the deposition of droplets on tablets with different shape, as well as the production and evolution of the liquid film on the surface of the tablets are numerically modeled. Spray droplets are simulated with a Discrete Droplets Method (DDM) Euler–Lagrange approach. Models for multi-component evaporation and particle/wall interaction are taken into account. The wall film is treated with a two-dimensional model incorporating submodels for interfacial shear force, film evaporation and heat transfer between film, solid wall and gas phase. Our simulations show how different physical parameters of the coating spray affect the coating process on a single tablet. For example, we analyze for the first time the deposition behavior of the droplets on the tablet. The outcome of our work provides a deeper understanding of the local interaction between the spray and the tablet bed, allowing a step forward in the design, scale-up, optimization and operation of industrial coating devices. Furthermore, it may serve as a basis for the combination with state-of-the-art DEM particle simulation tools.     

Effect of solid circulation rate on coating efficiency and agglomeration in circulating fluidized bed type coater

Circulating fluidized bed was proposed to be used as a coater, and coating experiments of glass beads with silica powder were performed in a circulating fluidized bed. Glass beads and silica powder were chosen as model particles, because their shape was almost spherical. The respective effects of gas flow rates supplied from a distributor and from an air nozzle for solid circulation, feed rate of powder suspension and particle content in the bed on coating efficiency and agglomeration are mainly discussed. Coating efficiency in circulating fluidized bed coater was correlated well with solid circulation time rather than with gas flow rates or solid circulation rate, while the agglomeration among core particles was mainly governed by solid circulation rate.     

Statistical modelling of the spouted bed coating process using positron emission particle tracking (PEPT) data

Coating of particles larger than about 1 mm can be achieved in a spouted bed, a particle mobilization device in which a strong particle circulation occurs, rapidly upwards in a lean central “spout” region and downwards in a slowly moving annular settled bed. In a spouted bed coater, a spray nozzle is placed at the base of the spout, spraying upwards into a distinct coating zone. The coating formation in a spouted bed is inter alia a function of (i) the particle motion, that is, how often and where particles enter and traverse the coating zone and (ii) the extent of droplet collection by individual particles passing through the coating zone. The coating model proposed here is based on the statistical history of individual particles, whose projected area governs the collection of spray droplets in the coating zone. Positron emission particle tracking (PEPT) has been used to determine the particle trajectories, the distribution of cycle times and the size and voidage of the spout. Whilst the model is not capable of delivering absolute values of coating mass a priori, it can predict deviations from a mean, which can itself be determined from an overall mass balance. To validate the model, a spouted bed coating process was studied in which coarse PVC spheres were coated with the hot‐melt coating material polyethylene glycol (PEG) 1500. Coating mass distributions, derived from the weight data of individual particles before and after manual coating removal, compared (for the studied conditions) very well with the predictions of the model.     

Hydrodynamic studies on three-phase fluidized bed using CFD analysis

Three-phase fluidization refers to fluidization of solid particles by co-current, upward flow of gas and liquid-phases for the purpose of bringing three-phases in contact in a single operation. Due to complications in understanding hydrodynamics of three-phase fluidized bed, CFD analysis is used to predict the hydrodynamics of it. In this study, liquid-phase is water which flows continuously, where as the gas phase is air which is distributed discretely throughout the bed. Ceramic particle of 1 mm diameter, density of 2650 kg/m³ is used as a solid phase. Excellent mixing, heat and mass transfer rates are the unique features of three-phase fluidized bed. The selection of distributor plays an important role in the quality of fluidization [1]. CFD model is created as the realistic representation of actual fluidized bed. The liquid and solid flow is represented by the mixture model. The air is injected from the bottom of the fluidized by means of discrete phase method (DPM). Simulation results are obtained by using porous jump and porous zone model to represent the distributor. It is found that porous zone model is best applicable in the industries, since stability of operating conditions is achieved even with non-uniform air, water flowrates and with different bed heights(100 mm, 200 mm, 300 mm, 400 mm and 500 mm).Simulated Pressure drop values of the fluidized bed have good agreement with the experimental findings. As the gas flowrate increases, the pressure drop in the column is decreases, provided the initial bed height, diameter of the column, and liquid flowrate are constant. This is due to decrease in density of the fluid medium in the bed by means of more gas hold up. The approach of the simulated values to the experimental values can be reduced with better understanding the nature of the fluidized bed.     

用VOF法模拟导向立体传质塔板罩内两相流

利用计算流体力学中的VOF方法对导向立体传质塔板（CTST-8）罩内气液两相流场进行了数值模拟,用Fluent的前处理软件gambit建立物理模型，在速度梯度大的地方采用局部加密网格，根据计算结果进行网格的自适应。用Fluent 6.1 对CTST-8在板孔气速为8.4 m/s，清液层高度为25 mm的工况下进行数值模拟。模拟结果显示了罩内流场的三维特性，并反映出罩内流场的相含率分布，速度分布以及压强分布等。将罩内压强的测量结果与模拟结果进行对比，两者吻合较好。说明本文数值模型具有较好的精度，可以用于CTST-8罩内两相流场的预测。     

DEM simulation of continuous tablet coating: Effects of tablet shape and fill level on inter-tablet coating variability

Tablet coating is a common pharmaceutical technique of applying a thin polymer-based film to a tablet or a granule containing active pharmaceutical ingredients (APIs). Inter- and intra-tablet variability of film coating is a critical issue in the production of solid oral dosage forms. In fact, inhomogeneity in the coating thickness can lead to significant variations in the delivery rate of active pharmaceutical ingredients and compromise the functional attributes of the tablet film. Although attempts have been made to use numerical approaches to analyze this complex problem, at present the uniformity of coating thickness is difficult to predict without expensive experimental work.The aim of this work is to analyze and understand the effects of tablet form and fill volume on the intra-tablet coating variability in a semi-continuous coating device. To this end, the Discrete Element Method was used to numerically reproduce the tablet motion inside a chamber of the coating pan. First, the material attributes of a sample placebo tablet were experimentally quantified in detail. Thereafter, three different tablet shapes, namely bi-convex, oval, and round, were modeled by means of the “glued spheres” method. The effect of three different fill volumes was then analyzed in terms of RT of the tablets under the coating spray, leading to a quantification of the intra-tablet coating variability for each particle shape. A detailed analysis of the tablets' velocities, both translational and rotational, on top of the tablet bed is presented. These results help to understand the dynamical behavior of the tablets under a spray gun that is essential for a satisfactory intra-tablet coating homogeneity. Finally, the various behaviors observed during the numerical simulations were addressed through a detailed analysis of the tablets' flow on the bed in terms of mean velocities and granular temperatures. The aim of this work is to demonstrate how a numerical simulation may be used for the development and design of continuous pharmaceutical tablet coating processes.     

大颗粒三相环隙气升式环流反应器流体力学行为

研究了大颗粒体系气升式环流反应器的流体力学行为,考察了表观气速和颗粒质量分数对床层膨胀高度、循环液速和固含率分布的影响。实验结果表明,按颗粒的运动状态不同可以将反应器内的流动分为3个区域,即固定床区域、膨胀床区域和循环床区域,各流动区域内的流动行为存在显著差异。随着颗粒质量浓度的增大,起始流化气速和最小循环气速均显著增大。基于三相流化床的流化模型和环流反应器的特点建立了相应的数学模型,对大颗粒三相气升式环流反应器的起始流化气速和最小循环气速进行了预测,模型预测值与实验测量值吻合良好。     

Measurement of particle concentration in a Wurster fluidized bed by electrical capacitance tomography sensors

It is essential to measure and monitor the particle flow characteristics in a Wurster fluidized bed to understand and optimize the coating processes. In this article, two electrical capacitance tomography (ECT) sensors are used to measure the particle concentration in different regions in a Wurster fluidized bed for the “cold” particle flows. One ECT sensor has a 12‐4 internal‐external electrodes and another has eight electrodes. The 12‐4‐electrode ECT sensor is used to measure the particle concentration in the annular fluidization region (outside of the Wurster tube) and the eight‐electrode ECT sensor is used to measure the particle flow in the central region (inside the Wurster tube). The effect of particle type, particle moisture, fluidization velocity, and geometrical parameters on the Wurster fluidization process is studied based on the two ECT measurements. The radial particle concentration profiles in the annular fluidization and central flow regions with different operation parameters are given. Fast Fourier Transform analysis of the particle concentration in the Wurster tube is performed with different superficial air velocities. The optimum operating ranges of the Wurster fluidization process for different particles are given. In the end of the article, computational fluids dynamics simulation results are given and used to compare with the measurement results by ECT for a typical Wurster fluidized bed. © 2014 American Institute of Chemical Engineers AIChE J 60: 4051–4064, 2014     

A two-scale model for discrete cell gravure roll coating

A two-dimensional model for predicting the fluid pickout and coated film thickness characteristics of a discrete cell direct gravure roll coater operating in reverse mode is derived. A novel multiscale approach is adopted for this purpose and the resulting equations solved numerically for inertia-less flow conditions. A system of stiff ordinary differential equations is found to be sufficient to capture the major gross flow features, while at the cell level the analysis is based on a finite element solution of the momentum and continuity equations. It represents the first such predictive model of its kind, with particular interest placed on the nature of both the pressure distribution and web-to-roll gap profile spanning the coating bead. The effect of key operating parameters, web-to-roll speed ratio, web-tension, wrap-angle, capillary number and cell-geometry, on the degree of fluid pickout from gravure cells and the coated film thickness is explored. Although an idealised model, the trends observed show qualitative agreement with existing experimental data collected on a small-scale gravure coating rig and point the way forward to the eventual formulation of a full three-dimensional predictive model of the process.     

Batch top-spray fluid bed coating: Scale-up insight using dynamic heat- and mass-transfer modelling

A mathematical model was developed for batch top-spray fluid bed coating processes based on Ronsse et al. [2007a, b. Combined population balance and thermodynamic modelling of the batch top-spray fluidised bed coating process. Part I—model development and validation. Journal of Food Engineering 78, 296-307; Combined population balance and thermodynamic modelling of the batch top-spray fluidised bed coating process. Part II—model and process analysis. Journal of Food Engineering 78, 308-322]. The model is based on one-dimensional discretisation of the fluid bed into a number of well-mixed control volumes. In each control volume, dynamic heat and mass balances were set up allowing the simulation of the contents of water vapour, water on core particles and deposited coating mass as well as fluidisation gas, particle and chamber wall temperature. The model was used to test different scale-up principles by comparing simulation results with experimental temperature and humidity data obtained from inorganic salt coating of placebo cores in three pilot fluid bed scales being a 0.5 kg small-scale (GEA Aeromatic-Fielder Strea-1), 4 kg medium-scale (GEA Niro MP-1) and 24 kg large-scale (GEA MP-2/3). Results show good agreement between simulated and experimental outlet fluidisation air temperature and humidity as well as bed temperature. Simulations reveal that vertical temperature and humidity gradients increase significantly with increasing scale and that in fluid beds as the simulated 900 kg (RICA-TEC Anhydro) production-scale, the gradients become too large to use the simple combined drying force/relative droplet size scale-up approach without also increasing the inlet fluidisation air temperature significantly. Instead, scale-up in terms of combinations of the viscous Stokes theory with simulated particle liquid layer profiles (obtained with the model) is suggested. In this way, the given fluid bed scale may be optimised in terms of low agglomeration tendency for a given process intensity across scale.     

颗粒形状对包衣设备内药片颗粒运动特性影响的数值模拟

为研究颗粒形状对包衣设备内药片颗粒运动特性的影响,基于离散单元法及自行编写的喷雾区颗粒检测算法,采用数值模拟的方法对五种不同形状(棒状、长椭球、扁椭球、双凸形和球形)的药片颗粒在包衣设备内的运动行为规律进行了研究。分析了颗粒形状对颗粒系统的能量、床面颗粒平动速度、颗粒温度及颗粒流在喷雾区域停留时间的分布及其相对标准差的影响。结果表明,颗粒形状对颗粒的平均动能、颗粒床面速度、颗粒温度、喷雾区域停留时间分布及颗粒间包衣均匀性有重要影响。除双凸形颗粒系统外,对于其他四种形状的颗粒系统,随着颗粒球形度增大,颗粒系统具有的动能、床面速度和颗粒温度均呈减小趋势。除棒状颗粒系统外,对于其余四种形状的颗粒系统,随着颗粒球形度增大,颗粒系统在包衣喷雾区域内平均停留时间减小,平均停留时间的相对标准差增大,包衣均匀性变差。与球形颗粒系统相比,非球形颗粒系统的包衣均匀性更好;药片颗粒形状对包衣设备内颗粒运动特性及颗粒之间的包衣均匀性有重要影响。     

A simplified mathematical modeling for thermo-catalytic decomposition of methane over carbon black catalyst in a fluidized bed reactor

A mathematical model for thermo-catalytic decomposition of methane over carbon black catalysts in a fluidized bed was proposed. The simplified isothermal, uniform flow model was considered and implemented into a computer code to predict the reactor performance. The experiment of methane decomposition into hydrogen and carbon was carried out in a fluidized bed of I.D of 0.055 m and height of 1.0 m. The range of reaction temperature was 850–900 °C, gas velocity was 1.0–3.0 U _mf , and catalyst loading was 50–200 g. The reaction parameters for model equation were determined from the curve fittings and the comparison of experimental data with simulation results showed good agreement for fluidized bed reactor system. From the simulation results, the fluidized bed performance with different operating conditions were obtained, and this simple model can be used to predict the performance of a larger scale fluidized bed reactor and also in determining the optimum operating conditions.     

Transient response of a solid-liquid model biological fluidised bed to a step change in fluid superficial velocity

The evolution of a solid-liquid model biological fluidised bed under a step change in fluid superficial velocity is described. During a transient step change, the fluidised bed divides into a top zone which remains at the initial porosity and a bottom zone which settles at the final porosity. The interface of discontinuity in porosity moves progressively upwards through the fluidised bed. The velocity at which the top of the fluidised bed expands or contracts and the upward velocity of the porosity transition interface depend only upon the initial and final states of the bed porosity and the fluid superficial velocity. This results in a linear evolution with time of the total bed height and the height of porosity transition interface. The proposed model is well suited to describe the transient response of low-density particles in a fluidised bed, such as encountered in biological systems, to a sudden change of liquid superficial velocity. The model was validated experimentally.