Nucleation regime map for liquid bound granules

Nucleation is the first step in granulation where the powder and liquid first contact. Two types of nucleation in wet granulation processes are proposed. Drop controlled nucleation, where one drop forms one nucleus, occurs when drops hitting the powder surface do not overlap (low spray flux Psi(a)) and the drop must wet quickly into the bed (short drop penetration time t(p)). If either criterion is not met, powder mixing characteristics will dominate (mechanical dispersion regime). Granulation experiments were performed with lactose powder, water, PEG200, and 7% HPC solution in a 6 L and a 25 L mixer granulator. Size distributions were measured as the drop penetration time and spray flux were varied. At short penetration times, decreasing Psi(a) caused the nuclei distribution to become narrower. When drop penetration time was high, the nuclei size distribution was broad independent of changes in dimensionless spray flux. Nucleation regime maps were plotted for each set of experiments in each mixer as a function of the dimensionless distribution width delta. The nucleation regime map demonstrates the interaction between drop penetration time and spray flux in nucleation. The narrowest distribution consistently occurred at low spray flux and low penetration time, proving the existence of the drop controlled regime. The nucleation regime map provides a rational basis for design and scale-up of nucleation and wetting in wet granulation.

Flow patterns in granulating systems

Particle flow patterns were investigated for wet granulation and dry powder mixing in ploughshare mixers using Positron Emission Particle Tracking (PEPT). In a 4-1 mixer, calcium carbonate with mean size 45 mum was granulated using a 50 wt.% solution of glycerol and water as binding fluid, and particle movement was followed using a 600-mum calcium hydroxy-phosphate tracer particle. In a 20-1 mixer, dry powder flow was studied using a 600-mum resin bead tracer particle to simulate the bulk polypropylene powder with mean size 600 mum. Important differences were seen between particle flow patterns for wet and dry systems. Particle speed relative to blade speed was lower in the wet system than in the dry system, with the ratios of average particle speed to blade tip speed for all experiments in the range 0.01-015. In the axial plane, the same particle motion was observed around each blade; this provides a significant advance for modelling flow in ploughshare mixers. For the future, a detailed understanding of the local velocity, acceleration and density variations around a plough blade will reveal the effects of flow patterns in granulating systems on the resultant distribution of granular product attributes such as size, density and strength. (C) 2002 Elsevier Science B.V All rights reserved.

Shapes and Dynamics of Blood Cells in Poiseuille and Shear Flows

The dynamics, shape, deformation, and orientation of red blood cells in microcirculation affect the rheology, flow resistance and transport properties of whole blood. This leads to important correlations of cellular and continuum scales. Furthermore, the dynamics of RBCs subject to different flow conditions and vessel geometries is relevant for both fundamental research and biomedical applications (e.g drug delivery). In this thesis, the behaviour of RBCs is investigated for different flow conditions via computer simulations. We use a combination of two mesoscopic particle-based simulation techniques, dissipative particle dynamics and smoothed dissipative particle dynamics. We focus on the microcapillary scale of several μm. At this scale, blood cannot be considered at the continuum but has to be studied at the cellular level. The connection between cellular motion and overall blood rheology will be investigated. Red blood cells are modelled as viscoelastic objects interacting hydrodynamically with a viscous fluid environment. The properties of the membrane, such as resistance against bending or shearing, are set to correspond to experimental values. Furthermore, thermal fluctuations are considered via random forces. Analyses corresponding to light scattering measurements are performed in order to compare to experiments and suggest for which situations this method is suitable. Static light scattering by red blood cells characterises their shape and allows comparison to objects such as spheres or cylinders, whose scattering signals have analytical solutions, in contrast to those of red blood cells. Dynamic light scattering by red blood cells is studied concerning its suitability to detect and analyse motion, deformation and membrane fluctuations. Dynamic light scattering analysis is performed for both diffusing and flowing cells. We find that scattering signals depend on various cell properties, thus allowing to distinguish different cells. The scattering of diffusing cells allows to draw conclusions on their bending rigidity via the effective diffusion coefficient. The scattering of flowing cells allows to draw conclusions on the shear rate via the scattering amplitude correlation. In flow, a RBC shows different shapes and dynamic states, depending on conditions such as confinement, physiological/pathological state and cell age. Here, two essential flow conditions are studied: simple shear flow and tube flow. Simple shear flow as a basic flow condition is part of any more complex flow. The velocity profile is linear and shear stress is homogeneous. In simple shear flow, we find a sequence of different cell shapes by increasing the shear rate. With increasing shear rate, we find rolling cells with cup shapes, trilobe shapes and quadrulobe shapes. This agrees with recent experiments. Furthermore, the impact of the initial orientation on the dynamics is studied. To study crowding and collective effects, systems with higher haematocrit are set up. Tube flow is an idealised model for the flow through cylindric microvessels. Without cell, a parabolic flow profile prevails. A single red blood cell is placed into the tube and subject to a Poiseuille profile. In tube flow, we find different cell shapes and dynamics depending on confinement, shear rate and cell properties. For strong confinements and high shear rates, we find parachute-like shapes. Although not perfectly symmetric, they are adjusted to the flow profile and maintain a stationary shape and orientation. For weak confinements and low shear rates, we find tumbling slippers that rotate and moderately change their shape. For weak confinements and high shear rates, we find tank-treading slippers that oscillate in a limited range of inclination angles and strongly change their shape. For the lowest shear rates, we find cells performing a snaking motion. Due to cell properties and resultant deformations, all shapes differ from hitherto descriptions, such as steady tank-treading or symmetric parachutes. We introduce phase diagrams to identify flow regimes for the different shapes and dynamics. Changing cell properties, the regime borders in the phase diagrams change. In both flow types, both the viscosity contrast and the choice of stress-free shape are important. For in vitro experiments, the solvent viscosity has often been higher than the cytosol viscosity, leading to a different pattern of dynamics, such as steady tank-treading. The stress-free state of a RBC, which is the state at zero shear stress, is still controversial, and computer simulations enable direct comparisons of possible candidates in equivalent flow conditions.

Numerical simulation of flow in a photocatalytic reactor containing baffles

In this study a new immobilized flat plate photocatalytic reactor for wastewater treatment has been investigated using computational fluid dynamics (CFD). The reactor consists of a reactor inlet, a reactive section where the catalyst is coated, and outlet parts. For simulation, the reactive section of the reactor was modelled with an array of baffles. In order to optimize the fluid mixing and reactor design, this study attempts to investigate the influence of baffles with differing heights on the flow field of the flat plate reactor. The results obtained from the simulation of a baffled flat plate reactor hydrodynamics for differing baffle heights for certain positions are presented. Under the conditions simulated, the qualitative flow features, such as the distribution of local stream lines, velocity contours, and high shear region, boundary layers separation, vortex formation, and the underlying mechanism are examined. At low and high Re numbers, the influence of baffle heights on the distribution of species mass fraction of a model pollutant are also highlighted. The simulation of qualitative and quantitative properties of fluid dynamics in a baffled reactor provides valuable insight to fully understand the effect of baffles and their role on the flow pattern, behaviour, and features of wastewater treatment using a photocatalytic reactor.

CFD Simulation of turbulence promoters in a water treatment reactor

A new immobilized flat plate photocatalytic reactor for wastewater treatment has been proposed in this study to avoid subsequent catalyst removal from the treated water. The reactor consists of an inlet, reactive section where catalyst is coated and an outlet parts. In order to optimize the fluid mixing and reactor design, this study aims to investigate the influence of baffles and its arrangement on the flat plate reactor hydrodynamics using computational fluid dynamics (CFD) simulation. For simulation, an array of baffles acting as turbulence promoters is inserted in the reactive zone of the reactor. In this regard, results obtained from the simulation of a baffled- flat plate photoreactor hydrodynamics for different baffle positions, heights and intervals are presented utilizing RNG k-ε turbulence model. Under the conditions simulated, the qualitative flow features, such as the development and separation of boundary layers, vortex formation, the presence of high shear regions and recirculation zones, and the underlying mechanism are examined. The influence of various baffle sizes on the distribution of pollutant concentration is also highlighted. The results presented here indicate that the spanning of recirculation increases the degree of interfacial distortion with a larger interfacial area between fluids which results in substantial enhancement in fluid mixing. The simulation results suggest that the qualitative and quantitative properties of fluid dynamics in a baffled reactor can be obtained which provides valuable insight to fully understand the effect of baffles and its arrangements on the flow pattern, behaviour, and feature.

Power-law viscous materials for analogue experiments: New data on the rheology of highly-filled silicone polymers

The selection of appropriate analogue materials is a central consideration in the design of realistic physical models. We investigate the rheology of highly-filled silicone polymers in order to find materials with a power-law strain-rate softening rheology suitable for modelling rock deformation by dislocation creep and report the rheological properties of the materials as functions of the filler content. The mixtures exhibit strain-rate softening behaviour but with increasing amounts of filler become strain-dependent. For the strain-independent viscous materials, flow laws are presented while for strain-dependent materials the relative importance of strain and strain rate softening/hardening is reported. If the stress or strain rate is above a threshold value some highly-filled silicone polymers may be considered linear visco-elastic (strain independent) and power-law strain-rate softening. The power-law exponent can be raised from 1 to ~3 by using mixtures of high-viscosity silicone and plasticine. However, the need for high shear strain rates to obtain the power-law rheology imposes some restrictions on the usage of such materials for geodynamic modelling. Two simple shear experiments are presented that use Newtonian and power-law strain-rate softening materials. The results demonstrate how materials with power-law rheology result in better strain localization in analogue experiments.

Simulation of the interaction between blood flow and atherosclerotic plaque

It has been well accepted that over 50% of cerebral ischemic events are the result of rupture of vulnerable carotid atheroma and subsequent thrombosis. Such strokes are potentially preventable by carotid interventions. Selection of patients for intervention is currently based on the severity of carotid luminal stenosis. It has been, however, widely accepted that luminal stenosis alone may not be an adequate predictor of risk. To evaluate the effects of degree of luminal stenosis and plaque morphology on plaque stability, we used a coupled nonlinear time-dependent model with flow-plaque interaction simulation to perform flow and stress/strain analysis for stenotic artery with a plaque. The Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations for the fluid. The Ogden strain energy function was used for both the fibrous cap and the lipid pool. The plaque Principal stresses and flow conditions were calculated for every case when varying the fibrous cap thickness from 0.1 to 2mm and the degree of luminal stenosis from 10% to 90%. Severe stenosis led to high flow velocities and high shear stresses, but a low or even negative pressure at the throat of the stenosis. Higher degree of stenosis and thinner fibrous cap led to larger plaque stresses, and a 50% decrease of fibrous cap thickness resulted in a 200% increase of maximum stress. This model suggests that fibrous cap thickness is critically related to plaque vulnerability and that, even within presence of moderate stenosis, may play an important role in the future risk stratification of those patients when identified in vivo using high resolution MR imaging.

Viscoelastic Properties of Nanocrystalline Films of Semiconducting Chalcogenides at Liquid/Liquid Interface

The interfacial shear rheological properties of a continuous single-crystalline film of CuS and a 3D particulate gel of CdS nanoparticles (3−5 nm in diameter) formed at toluene−water interfaces have been studied. The ultrathin films (50 nm in thickness) are formed in situ in the shear cell through a reaction at the toluene−water interface between a metal−organic compound in the organic layer and an appropriate reagent for sulfidation in the aqueous layer. Linear viscoelastic spectra of the nanofilms reveal solid-like rheological behavior with the storage modulus higher than the loss modulus over the range of angular frequencies probed. Large strain amplitude sweep measurements on the CdS nanofilms formed at different reactant concentrations suggest that they form a weakly flocculated gel. Under steady shear, the films exhibit a yield stress, followed by a steady shear thinning at high shear rates. The viscoelastic and flow behavior of these films that are in common with those of many 3D “soft” materials like gels, foams, and concentrated colloidal suspensions can be described by the “soft” glassy rheology model.

Individual variation in in vitro efficacy of antiplatelet medication - aspirin and clopidogrel

Antiplatelet medication is known to decrease adverse effects in patients with atherothrombotic disease. However, despite ongoing antiplatelet medication considerable number of patients suffer from atherothrombotic events. The aims of the study were 1) to evaluate the individual variability in platelet functions and compare the usability of different methods in detecting it, 2) to assess variability in efficacy of antiplatelet medication with aspirin (acetylsalicylic acid) or the combination of aspirin and clopidogrel and 3) to investigate the main genetic and clinical variables as well as potential underlying mechanisms of variability in efficacy of antiplatelet medication. In comparisons of different platelet function tests in 19 healthy individuals PFA-100® correlated with traditional methods of measuring platelet function and was thus considered appropriate for testing individual variability in platelet activity. Efficacy of ongoing 100mg aspirin daily was studied in 101 patients with coronary artery disease (CAD). Aspirin response was measured with arachidonic acid (AA)-induced platelet aggregation, which reflects cyclo-oxygenase (COX)-1 dependent thromboxane (Tx) A2 formation, and PFA-100®, which evaluates platelet activation under high shear stress in the presence of collagen and epinephrine. Five percent of patients failed to show inhibition of AA-aggregation and 21% of patients had normal PFA-100® results despite aspirin and were thus considered non-responders to aspirin. Interestingly, the two methods of assessing aspirin efficacy, platelet aggregation and PFA-100®, detected different populations as being aspirin non-responders. It could be postulated that PFA-100® actually measures enhanced platelet function, which is not directly associated with TxA2 inhibition exerted by aspirin. Clopidogrel efficacy was assessed in 50 patients who received a 300mg loading dose of clopidogrel 2.5 h prior to percutaneous coronary intervention (PCI) and in 51 patients who were given a loading dose of 300mg combined with a five day treatment of 75mg clopidogrel daily mimicking ongoing treatment. Clopidogrel response was assessed with ADP-induced aggregations, due to its mechanism of action as an inhibitor of ADP-induced activation. When patients received only a loading dose of clopidogrel prior to PCI, 40% did not gain measurable inhibition of their ADP-induced platelet activity (inhibition of 10% or less). Prolongation of treatment so that all patients had reached a plateau of inhibition exerted by clopidogrel, decreased the incidence of non-responders to 20%. Polymorphisms of COX-1 and GP VI, as well as diabetes and female gender, were associated with decreased in vitro aspirin efficacy. Diabetes also impaired the in vitro efficacy of short-term clopidogrel. Decreased response to clopidogrel was associated with limited inhibition by ARMX, an antagonist of P2Y12-receptor, suggesting the reason for clopidogrel resistance to be receptor-dependent. Conclusions: Considerable numbers of CAD patients were non-responders either to aspirin, clopidogrel or both. In the future, platelet function tests may be helpful to individually select effective and safe antiplatelet medication for these patients.

A study on boundary-layer transition induced by free-stream turbulence

Boundary-layer transition at different free-stream turbulence levels has been investigated using the particle-image velocimetry technique. The measurements show organized positive and negative fluctuations of the streamwise fluctuating velocity component, which resemble the forward and backward jet-like structures reported in the direct numerical simulation of bypass transition. These fluctuations are associated with unsteady streaky structures. Large inclined high shear-layer regions are also observed and the organized negative fluctuations are found to appear consistently with these inclined shear layers, along with highly inflectional instantaneous streamwise velocity profiles. These inflectional velocity profiles are similar to those in the ribbon-induced boundary-layer transition. An oscillating-inclined shear layer appears to be the turbulent spot-precursor. The measurements also enabled to compare the actual turbulent spot in bypass transition with the simulated one. A proper orthogonal decomposition analysis of the fluctuating velocity field is carried out. The dominant flow structures of the organized positive and negative fluctuations are captured by the first few eigenfunction modes carrying most of the fluctuating energy. The similarity in the dominant eigenfunctions at different Reynolds numbers suggests that the flow prevails its structural identity even in intermittent flows. This analysis also indicates the possibility of the existence of a spatio-temporal symmetry associated with a travelling wave in the flow.

Total internal reflection Raman spectroscopy of poly(alpha-olefin) oils in a lubricated contact

The rheology of a poly(alpha-olefin) base oil (PAO) in a sliding point contact has been investigated by total internal reflection (TIR) Raman spectroscopy. TIR Raman has the sensitivity to analyse nanometer-thick lubricant films in a tribological contact. The Raman signal generated from the sliding contact was used to determine the lubricant film thickness. The experimentally obtained film thicknesses were compared with theoretical calculations and a transition from Newtonian to non-Newtonian behaviour was observed at high shear rates. The Raman spectra showed no significant changes in the conformation of the PAO chains under the applied conditions of pressure and shear, but the polarisation dependence of the spectra revealed a preferred orientation of the hydrocarbon side chains in the shear-thinned region. Monolayers formed by a boundary lubricant, arachidic acid, dissolved in the PAO could be detected on the surfaces in the elastohydrodynamic regime.

Secondary instability as a possible mechanism for clear-air turbulence: a case study

Many theories and mechanisms have been proposed to explain the phenomenon of clear-air turbulence (CAT), and some of them have been successful in predicting light, moderate and, in some cases, severe turbulence. It is only recently that skill in the forecasting of the severe form of CAT, which could lead to injuries to passengers and damage to aircraft, has improved. Recent observations and simulations suggest that some severe to extreme turbulence could be caused by horizontal vortex tubes resulting from secondary instabilities of regions of high shear in the atmosphere. We have conducted direct numerical simulations to understand the scale relationship between primary structures (larger-scale structures related to one of the causes mentioned above) and secondary structures (smaller-sized, shear structures of the size of aircraft). From shear layer simulations, we find that the ratio of sizes of primary and secondary vortices is of the right order to generate aircraft-scale vortex tubes from typical atmospheric shear layers. We have also conducted simulations with a mesoscale atmospheric model, to understand possible causes of turbulence experienced by a flight off the west coast of India. Our simulations show the occurrence of primary flow structures related to synoptic conditions around the time of the incident. The evidence presented for this mechanism also has implications for possible methods of detection and avoidance of severe CAT.

Large carbon cluster thin film gauges for measuring aerodynamic heat transfer rates in hypersonic shock tunnels

Different types of Large Carbon Cluster (LCC) layers are synthesized by a single-step pyrolysis technique at various ratios of precursor mixture. The aim is to develop a fast responsive and stable thermal gauge based on a LCC layer which has relatively good electrical conduction in order to use it in the hypersonic flow field. The thermoelectric property of the LCC layer has been studied. It is found that these carbon clusters are sensitive to temperature changes. Therefore suitable thermal gauges were developed for blunt cone bodies and were tested in hypersonic shock tunnels at a flow Mach number of 6.8 to measure aerodynamic heating. The LCC layer of this thermal gauge encounters high shear forces and a hostile environment for test duration in the range of a millisecond. The results are favorable to use large carbon clusters as a better sensor than a conventional platinum thin film gauge in view of fast responsiveness and stability.

Effect of a mesh on boundary layer transitions induced by free-stream turbulence and an isolated roughness element

Streamwise streaks, their lift-up and streak instability are integral to the bypass transition process. An experimental study has been carried out to find the effect of a mesh placed normal to the flow and at different wall-normal locations in the late stages of two transitional flows induced by free-stream turbulence (FST) and an isolated roughness element. The mesh causes an approximately 30% reduction in the free-stream velocity, and mild acceleration, irrespective of its wall-normal location. Interestingly, when located near the wall, the mesh suppresses several transitional events leading to transition delay over a large downstream distance. The transition delay is found to be mainly caused by suppression of the lift-up of the high-shear layer and its distortion, along with modification of the spanwise streaky structure to an orderly one. However, with the mesh well away from the wall, the lifted-up shear layer remains largely unaffected, and the downstream boundary layer velocity profile develops an overshoot which is found to follow a plane mixing layer type profile up to the free stream. Reynolds stresses, and the size and strength of vortices increase in this mixing layer region. This high-intensity disturbance can possibly enhance transition of the accelerated flow far downstream, although a reduction in streamwise turbulence intensity occurs over a short distance downstream of the mesh. However, the shape of the large-scale streamwise structure in the wall-normal plane is found to be more or less the same as that without the mesh.

3D Finite Volume Scheme for Czochralski Single Crystal Growth

Czochralski (Cz) technique, which is used for growing single crystals, has dominated the production of single crystals for electronic applications. The Cz growth process involves multiple phases, moving interface and three-dimensional behavior. Much has been done to study these phenomena by means of numerical methods as well as experimental observations. A three-dimensional curvilinear finite volume based algorithm has been developed to model the Cz process. A body-fitted transformation based approach is adopted in conjunction with a multizone adaptive grid generation (MAGG) technique to accurately handle the three-dimensional problems of phase-change in irregular geometries with free and moving surfaces. The multizone adaptive model is used to perform a three-dimensional simulation of the Cz growth of silicon single crystals.Since the phase change interface are irregular in shape and they move in response to the solution, accurate treatment of these interfaces is important from numerical accuracy point of view. The multizone adaptive grid generation (MAGG) is the appropriate scheme for this purpose. Another challenge encountered is the moving and periodic boundary conditions, which is essential to the numerical solution of the governing equations. Special treatments are implemented to impose the periodic boundary condition in a particular direction and to determine the internal boundary position and shape varying with the combination of ambient physicochemical transport process and interfacial dynamics. As indicated above that the applications and processes characterized by multi-phase, moving interfaces and irregular shape render the associated physical phenomena three-dimensional and unsteady. Therefore a generalized 3D model rather than a 2D simulation, in which the governing equations are solved in a general non-orthogonal coordinate system, is constructed to describe and capture the features of the growth process. All this has been implemented and validated by using it to model the low pressure Cz growth of silicon. Accuracy of this scheme is demonstrated by agreement of simulation data with available experimental data. Using the quasi-steady state approximation, it is shown that the flow and temperature fields in the melt under certain operating conditions become asymmetric and unsteady even in the absence of extrinsic sources of asymmetry. Asymmetry in the flow and temperature fields, caused by high shear initiated phenomena, affects the interface shape in the azimuthal direction thus results in the thermal stress distribution in the vicinity, which has serious implications from crystal quality point of view.