Measuring the flow properties of small powder samples using an avalanche tester

The feasibility of using a small-scale avalanche tester to measure the flow properties of pharmaceutical lactose powders was examined. The modes of behaviour observed in larger systems were displayed and showed a clear distinction between angular, free-flowing particles and more spherical particles of similar flow characteristics. Angular Lactohale LH100 particles showed slumping behaviour at a rotational frequency of 0.33Hz which disappeared at higher frequencies. Spherical lactose powder with a similar flow function to LH100 only showed rolling behaviour under the same conditions, as did more cohesive powders LH200 and LH300. Further investigation of the LH100 data using Fast Fourier analysis showed that the slumping frequency was 1/10th of the rotational frequency.

The impact of powder flowing properties on dosator filling system operation

Purpose: To study the impact of powder flow properties on dosator filling systems, with particular focus on improvements in dose weight accuracy and repeatability. Method: This study evaluates a range of critical powder flow properties such as: flow function, cohesion, wall friction, adhesion to wall surfaces, density/compressibility data, stress ratio “K” and gas permeability. The characterisations of the powders considered in this study were undertaken using an annular shear cell using a sample size of 0.5 litres. This tester also incorporated the facility to measure bed expansion during shear in addition to contraction under consolidation forces. A modified Jenike type linear wall friction tester was used to develop the failure loci for the powder sample in conjunction with multiple wall samples (representing a variety of material types and surface finishes). Measurements of the ratio of applied normal stress versus lateral stress were determined using a piece of test equipment specifically designed for the purpose. Results: The correct characterisation of powders and the incorporation of this data into the design of process equipment are recognised as critical for reliable and accurate operation. An example of one aspect of this work is the stress ratio “K”. This characteristic is not well understood or correctly interpreted in many cases – despite its importance. Fig 1 [Omitted] (illustrates a sample of test data. The slope of the line gives the stress ratio in a uniaxial compaction system – indicating the behaviour of the material under compaction during dosing processes. Conclusions: A correct assessment of the bulk powder properties for a given formulation can allow prediction of: cavity filling behaviour (and hence dosage), efficiency of release from dosator, and strength and stability of extruded dose en route to capsule filling Influences over the effectiveness of dosator systems have been shown to be impacted upon by: bed pre-compaction history, gas permeability in the bed (with respect to local density effects), and friction effects for materials of construction for dosators

Modelling of dosator operation for improved manufacturing performance

Purpose: To develop an improved mathematical model for the prediction of dose accuracy of Dosators - based upon the geometry of the machine in conjunction with measured flow properties of the powder. Methods: A mathematical model has been created, based on a analytical method of differential slices - incorporating measured flow properties. The key flow properties of interest in this investigation were: flow function, effective angle of wall friction, wall adhesion, bulk density, stress ratio K and permeability. To simulate the real process and (very importantly) validate the model, a Dosator test-rig has been used to measure the forces acting on the Dosator during the filling stage, the force required to eject the dose and the dose weight. Results: Preliminary results were obtained from the Dosator test-rig. Figure 1 [Omitted] shows the dose weight for different depths to the bottom of the powder bed at the end of the stroke and different levels of pre-compaction of the powder bed. A strong influence over dose weight arising from the proximity between the Dosator and the bottom of the powder bed at the end of the stroke and the conditions of the powder bed has been established. Conclusions: The model will provide a useful tool to predict dosing accuracy and, thus, optimise the future design of Dosator based equipment technology – based on measured bulk properties of the powder to be handled. Another important factor (with a significant influence) on Dosator processes, is the condition of the powder bed and the clearance between the Dosator and the bottom of the powder bed.

Modelling Cohesive-Frictional Particulate Solids for Bulk Handling Applications

Many powders and particulate solids are cohesive in nature and the strength often exhibits dependence on the consolidation stress. As a result, the stress history in the material leading up to a handling scenario needs to be considered when evaluating its handleability. This paper outlines the development of a DEM contact model accounting for plasticity and adhesion force, which is shown to be suitable for modelling the stress history dependent cohesive strength. The model was used to simulate the confined consolidation and the subsequent unconfined loading of iron ore fines with particle sizes up to 1.18mm. The predicted flow function was found to be comparable to the experimental results. 

Micromechanical analysis of cohesive granular materials using the discrete element method with an adhesive elasto-plastic contact model

An adhesive elasto-plastic contact model for the discrete element method with three dimensional non-spherical particles is proposed and investigated to achieve quantitative prediction of cohesive powder flowability. Simulations have been performed for uniaxial consolidation followed by unconfined compression to failure using this model. The model has been shown to be capable of predicting the experimental flow function (unconfined compressive strength vs. the prior consolidation stress) for a limestone powder which has been selected as a reference solid in the Europe wide PARDEM research network. Contact plasticity in the model is shown to affect the flowability significantly and is thus essential for producing satisfactory computations of the behaviour of a cohesive granular material. The model predicts a linear relationship between a normalized unconfined compressive strength and the product of coordination number and solid fraction. This linear relationship is in line with the Rumpf model for the tensile strength of particulate agglomerate. Even when the contact adhesion is forced to remain constant, the increasing unconfined strength arising from stress consolidation is still predicted, which has its origin in the contact plasticity leading to microstructural evolution of the coordination number. The filled porosity is predicted to increase as the contact adhesion increases. Under confined compression, the porosity reduces more gradually for the load-dependent adhesion compared to constant adhesion. It was found that the contribution of adhesive force to the limiting friction has a significant effect on the bulk unconfined strength. The results provide new insights and propose a micromechanical based measure for characterising the strength and flowability of cohesive granular materials. 

A DEM study of silo discharge of a cohesive solid

Bulk handling of powders and granular solids is common in many industries and often gives rise to handling difficulties especially when the material exhibits complex cohesive behaviour. For example, high storage stresses in a silo can lead to high cohesive strength of the stored solid, which may in turn cause blockages such as ratholing or arching near the outlet during discharge. This paper presents a Discrete Element Method study of discharge of a granular solid with varying levels of cohesion from a flat-bottomed silo. The DEM simulations were conducted using the commercial EDEM code with a recently developed DEM contact model for cohesive solids implemented through an API. The contact model is based on an elasto-plastic contact with adhesion and uses hysteretic non-linear loading and unloading paths to model the elastic-plastic contact deformation. The adhesion parameter is a function of the maximum contact overlap. The model has been shown to be able to predict the stress history dependent behaviour depicted by a flow function of the material. The effects of cohesion on the discharge rate and flow pattern in the silo are investigated. The predicted discharge rates are compared for the varying levels of cohesion and the effect of adhesion is evaluated. The ability of the contact model to qualitatively predict the phenomena that are present in the discharge of a silo has been shown with the salient feature of mixed flow from a flat bottomed hopper identified in the simulation.

Experimental and discrete element modelling of cohesive iron ore fines

In this study, the behaviour of iron ore fines with varying levels of adhesion was investigated using a confined compression test and a uniaxial test. The uniaxial test was conducted using the semi-automated uniaxial EPT tester in which the cohesive strength of a bulk solid is evaluated from an unconfined compression test following a period of consolidation to a pre-defined vertical stress. The iron ore fines were also tested by measuring both the vertical and circumferential strains on the cylindrical container walls under vertical loading in a separate confined compression tester - the K0 tester, to determine the lateral pressure ratio. Discrete Element Method simulations of both experiments were carried out and the predictions were compared with the experimental observations. A recently developed DEM contact model for cohesive solids, an Elasto-Plastic Adhesive model, was used. This particle contact model uses hysteretic non-linear loading and unloading paths and an adhesion parameter which is a function of the maximum contact overlap. The model parameters for the simulations are phenomenologically based to reproduce the key bulk characteristics exhibited by the solid. The simulation results show a good agreement in capturing the stress history dependent behaviour depicted by the flow function of the cohesive iron ore fines while also providing a reasonably good match for the lateral pressure ratio observed during the confined compression K0 tests. This demonstrates the potential for the DEM model to be used in the simulation of bulk handling applications.

Option pricing in market models driven by telegraph processes with jumps

Esta tesis está dividida en dos partes: en la primera parte se presentan y estudian los procesos telegráficos, los procesos de Poisson con compensador telegráfico y los procesos telegráficos con saltos. El estudio presentado en esta primera parte incluye el cálculo de las distribuciones de cada proceso, las medias y varianzas, así como las funciones generadoras de momentos entre otras propiedades. Utilizando estas propiedades en la segunda parte se estudian los modelos de valoración de opciones basados en procesos telegráficos con saltos. En esta parte se da una descripción de cómo calcular las medidas neutrales al riesgo, se encuentra la condición de no arbitraje en este tipo de modelos y por último se calcula el precio de las opciones Europeas de compra y venta.

Fast decay of the velocity autocorrelation function in dense shear flow of inelastic hard spheres

We find in complementary experiments and event-driven simulations of sheared inelastic hard spheres that the velocity autocorrelation function psi(t) decays much faster than t(-3/2) obtained for a fluid of elastic spheres at equilibrium. Particle displacements are measured in experiments inside a gravity-driven flow sheared by a rough wall. The average packing fraction obtained in the experiments is 0.59, and the packing fraction in the simulations is varied between 0.5 and 0.59. The motion is observed to be diffusive over long times except in experiments where there is layering of particles parallel to boundaries, and diffusion is inhibited between layers. Regardless, a rapid decay of psi(t) is observed, indicating that this is a feature of the sheared dissipative fluid, and is independent of the details of the relative particle arrangements. An important implication of our study is that the non-analytic contribution to the shear stress may not be present in a sheared inelastic fluid, leading to a wider range of applicability of kinetic theory approaches to dense granular matter.

Velocity distribution function for a dilute granular material in shear flow

The velocity distribution function for the steady shear flow of disks (in two dimensions) and spheres (in three dimensions) in a channel is determined in the limit where the frequency of particle-wall collisions is large compared to particle-particle collisions. An asymptotic analysis is used in the small parameter epsilon, which is naL in two dimensions and na(2)L in three dimensions, where; n is the number density of particles (per unit area in two dimensions and per unit volume in three dimensions), L is the separation of the walls of the channel and a is the particle diameter. The particle-wall collisions are inelastic, and are described by simple relations which involve coefficients of restitution e(t) and e(n) in the tangential and normal directions, and both elastic and inelastic binary collisions between particles are considered. In the absence of binary collisions between particles, it is found that the particle velocities converge to two constant values (u(x), u(y)) = (+/-V, O) after repeated collisions with the wall, where u(x) and u(y) are the velocities tangential and normal to the wall, V = (1 - e(t))V-w/(1 + e(t)), and V-w and -V-w, are the tangential velocities of the walls of the channel. The effect of binary collisions is included using a self-consistent calculation, and the distribution function is determined using the condition that the net collisional flux of particles at any point in velocity space is zero at steady state. Certain approximations are made regarding the velocities of particles undergoing binary collisions :in order to obtain analytical results for the distribution function, and these approximations are justified analytically by showing that the error incurred decreases proportional to epsilon(1/2) in the limit epsilon --> 0. A numerical calculation of the mean square of the difference between the exact flux and the approximate flux confirms that the error decreases proportional to epsilon(1/2) in the limit epsilon --> 0. The moments of the velocity distribution function are evaluated, and it is found that [u(x)(2)] --> V-2, [u(y)(2)] similar to V-2 epsilon and -[u(x)u(y)] similar to V-2 epsilon log(epsilon(-1)) in the limit epsilon --> 0. It is found that the distribution function and the scaling laws for the velocity moments are similar for both two- and three-dimensional systems.

Flow-mediated dilatation of the brachial artery is a poorly reproducible indicator of microvascular function in Type I diabetes mellitus

Aim: Two Type I diabetes and control group comparator studies were conducted to assess the reproducibility of FMD and to analyse blood flow data normally discarded during FMD measurement.

Design: The studies were sequential and differed only with regard to operator and ultrasound machine. Seventy-two subjects with diabetes and 71 controls were studied in total.

Methods: Subjects had FMD measured conventionally. Blood velocity waveforms were averaged over 10 pulses post forearm ischaemia and their component frequencies analysed using the wavelet transform, a mathematical tool for waveform analysis. The component frequencies were grouped into 11 bands to facilitate analysis.

Results: Subjects were well-matched between studies. In Study 1, FMD was significantly impaired in subjects with Type I diabetes vs. controls (median 4.35%, interquartile range 3.10-4.80 vs. 6.50, 4.79-9.42, P < 0.001). No differences were detected between groups in Study 2, however. However, analysis of blood velocity waveforms yielded significant differences between groups in two frequency bands in each study.

Conclusions: This report highlights concerns over the reproducibility of FMD measures. Further work is required to fully elucidate the role of analysing velocity waveforms after forearm ischaemia.

Susceptibility of streptozotocin-induced diabetic rat retinal function and ocular blood flow to acute intraocular pressure challenge.

PURPOSE: To consider whether STZ-induced hyperglycemia renders rat retinal function and ocular blood flow more susceptible to acute intraocular pressure (IOP) challenge.

METHODS: Retinal function (electroretinogram, ERG) was measured during acute IOP challenge (10-100 mmHg, 5 mmHg increments, 3 min/step, vitreal cannulation) in adult Long-Evans rats (6-week old, citrate: n=6, STZ: n=10) 4 weeks after citrate buffer or streptozotocin (STZ, 65 mg/kg, blood glucose > 15 mmol/l) injection. At each IOP, dim and bright flash (-4.56, -1.72 log cd.s.m^-2) ERG responses were recorded to measure inner retinal and ON-bipolar cell function, respectively. Ocular blood flow (laser Doppler flowmetry, citrate; n=6, STZ; n=10) was also measured during acute IOP challenge. Retinae were isolated for qPCR analysis of nitric oxide synthase mRNA expression endothelial, eNos; inducible, iNos; neuronal, nNos).

RESULTS: STZ-induced diabetes increased the susceptibility of inner retinal (IOP at 50% response, 60.1, CI: 57.0-62.0 mmHg vs. citrate: 67.5, CI: 62.1-72.4 mmHg) and ON-bipolar cell function (STZ: 60.3, CI: 58.0-62.8 mmHg vs. citrate: 65.1, CI: 58.0-62.78 mmHg) and ocular blood flow (43.9, CI: 40.8-46.8 vs. citrate: 53.4, CI: 50.7-56.1 mmHg) to IOP challenge. Citrate eyes showed elevated eNos mRNA (+49.7%) after IOP stress, an effect not found in STZ-diabetic eyes (-5.7%, P<0.03). No difference was observed for iNos or nNos (P>0.05) following IOP elevation.

CONCLUSIONS: STZ-induced diabetes increased functional susceptibility during acute IOP challenge. This functional vulnerability is associated with a reduced capacity for diabetic eyes to upregulate eNOS expression and to autoregulate blood flow in response to stress.

Effects of Cerebral Blood Flow and PETCO2 on Cognitive Function During Passive Heat Stress

This thesis tested whether cognitive performance during passive heat stress may be affected by changes in cerebrovascular variables as opposed to strictly thermally-induced changes. A pharmacological reduction in cerebral blood flow (CBF) using indomethacin along with a hypocapnia-induced CBF reduction during passive heat stress (Tre ~1.5°C above baseline) were used to investigate any cerebrovascular-mediated changes in cognitive performance. Repeated measures analysis of variance indicated that One-Touch Stockings of Cambridge (OTS) performance was not affected by a significant reduction in CBF during passive heat stress. More specifically, OTS accuracy measures did not change as a result of either a reduction in CBF or increasing passive heat stress. However, it was found that OTS response time indices improved with increasing passive heat stress independent of CBF changes. In conclusion, a significant reduction in CBF does not cause additional changes in performance of an executive functioning task during severe passive heat stress.

Optimal reactive power flow via the modified barrier Lagrangian function approach

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)