Shear effects on the properties and separation characteristics of whey protein precipitates

Selective isoelectric whey protein precipitation and aggregation is carried out at laboratory scale in a standard configuration batch agitation vessel. Geometric scale-up of this operation is implemented on the basis of constant impeller power input per unit volume and subsequent clarification is achieved by high speed disc-stack centrifugation. Particle size and fractal geometry are important in achieving efficient separation while aggregates need to be strong enough to resist the more extreme levels of shear that are encountered during processing, for example through pumps, valves and at the centrifuge inlet zone. This study investigates how impeller agitation intensity and ageing time affect aggregate size, strength, fractal dimension and hindered settling rate at laboratory scale in order to determine conditions conducive for improved separation. Particle strength is measured by observing the effects of subjecting aggregates to moderate and high levels of process shear in a capillary rig and through a partially open ball-valve respectively. The protein precipitate yield is also investigated with respect to ageing time and impeller agitation intensity. A pilot scale study is undertaken to investigate scale-up and how agitation vessel shear affects centrifugal separation efficiency. Laboratory scale studies show that precipitates subject to higher impeller shear-rates during the addition of the precipitation agent are smaller but more compact than those subject to lower impeller agitation and are better able to resist turbulent breakage. They are thus more likely to provide a better feed for more efficient centrifugal separation. Protein precipitation yield improves significantly with ageing, and 50 minutes of ageing is required to obtain a 70 - 80% yield of α-lactalbumin. Geometric scale-up of the agitation vessel at constant power per unit volume results in aggregates of broadly similar size exhibiting similar trends but with some differences due to the absence of dynamic similarity due to longer circulation time and higher tip speed in the larger vessel. Disc stack centrifuge clarification efficiency curves show aggregates formed at higher shear-rates separate more efficiently, in accordance with laboratory scale projections. Exposure of aggregates to highly turbulent conditions, even for short exposure times, can lead to a large reduction in particle size. Thus, improving separation efficiencies can be achieved by the identification of high shear zones in a centrifugal process and the subsequent elimination or amelioration of such.

The formulation, testing and stability of 16% fat cream liqueurs

Cream liqueurs manufactured by a one-step process, where alcohol was added before homogenisation, were more stable than those processed by a two -step process which involved addition of alcohol after homogenisation. Using the one-step process, it was possible to produce creaming-stable liqueurs by using one pass through a homogeniser (27.6 MPa) equipped with "liquid whirl" valves. Test procedures to characterise cream liqueurs and to predict shelf life were studied in detail. A turbidity test proved simple, rapid and sensitive for characterising particle size and homogenisation efficiency. Prediction of age thickening/gelation in cream liqueurs during incubation at 45 °C depended on the age of the sample when incubated. Samples that gelled at 45 °C may not do so at ambient temperature. Commercial cream liqueurs were similar in gross chemical composition, and unlike experimentally produced liqueurs, these did not exhibit either age-gelation at ambient or elevated temperatures. Solutions of commercial sodium caseinates from different sources varied in their calcium sensitivity. When incorporated into cream liqueurs, caseinates influenced the rate of viscosity increase, coalescence and, possibly, gelation during incubated storage. Mild heat and alcohol treatment modified the properties of caseinate used to stabilise non-alcoholic emulsions, while the presence of alcohol in emulsions was important in preventing clustering of globules. The response to added trisodium citrate varied. In many cases, addition of the recommended level (0.18%) did not prevent gelation. Addition of small amounts of NaOH with 0.18 % trisodium citrate before homogenisation was beneficial. The stage at which citrate was added during processing was critical to the degree of viscosity increase (as opposed to gelation) in the product during 45 °C incubation. The component responsible for age-gelation was present in the milk-solids non fat portion of the cream and variations in the creams used were important in the age-gelation phenomenon Results indicated that, in addition to possibly Ca++, the micellar casein portion of serum may play a role in gelation. The role of the low molecular weight surfactants, sodium stearoyl lactylate and monodiglycerides in preventing gelation, was influenced by the presence of trisodium citrate. Clustering of fat globules and age-gelation were inhibited when 0.18 % citrate was included. Inclusion of sodium stearoyl lactylate, but not monodiglycerides, reduced the extent of viscosity increase at 45 °C in citrate containing liqueurs.

Preparation and characterisation of ceria particles

Cerium dioxide (ceria) nanoparticles have been the subject of intense academic and industrial interest. Ceria has a host of applications but academic interest largely stems from their use in the modern automotive catalyst but it is also of interest because of many other application areas notably as the abrasive in chemical-mechanical planarisation of silicon substrates. Recently, ceria has been the focus of research investigating health effects of nanoparticles. Importantly, the role of non-stoichiometry in ceria nanoparticles is implicated in their biochemistry. Ceria has well understood non-stoichiometry based around the ease of formation of anion vacancies and these can form ordered superstructures based around the fluorite lattice structure exhibited by ceria. The anion vacancies are associated with localised or small polaron states formed by the electrons that remain after oxygen desorption. In simple terms these electrons combine with Ce4+ states to form Ce3+ states whose larger ionic radii is associated with a lattice expansion compared to stoichiometric CeO2. This is a very simplistic explanation and greater defect chemistry complexity is suggested by more recent work. Various authors have shown that vacancies are mobile and may result in vacancy clustering. Ceria nanoparticles are of particular interest because of the high activity and surface area of small particulates. The sensitivity of the cerium electronic band structure to environment would suggest that changes in the properties of ceria particles at nanoscale dimensions might be expected. Notably many authors report a lattice expansion with reducing particle size (largely confined to sub-10 nm particles). Most authors assign increased lattice dimensions to the presence of a surface stable Ce2O3 type layer at low nanoparticle dimensions. However, our understanding of oxide nanoparticles is limited and their full and quantitative characterisation offers serious challenges. In a series of chemical preparations by ourselves we see little evidence of a consistent model emerging to explain lattice parameter changes with nanoparticle size. Based on these results and a review of the literature it is worthwhile asking if a model of surface enhanced defect concentration is consistent with known cerium/cerium oxide chemistries, whether this is applicable to a range of different synthesis methods and if a more consistent description is possible. In Chapter one the science of cerium oxide is outlined including the crystal structure, defect chemistry and different oxidation states available. The uses and applications of cerium oxide are also discussed as well as modelling of the lattice parameter and the doping of the ceria lattice. Chapter two describes both the synthesis techniques and the analytical methods employed to execute this research. Chapter three focuses on high surface area ceria nano-particles and how these have been prepared using a citrate sol-gel precipitation method. Changes to the particle size have been made by calcining the ceria powders at different temperatures. X-ray diffraction methods were used to determine their lattice parameters. The particles sizes were also assessed using transmission electron microscopy (TEM), scanning electron microscopy (SEM), and BET, and, the lattice parameter was found to decrease with decreasing particle size. The results are discussed in light of the role played by surface tension effects. Chapter four describes the morphological and structural characterization of crystalline CeO2 nanoparticles prepared by forward and reverse precipitation techniques and compares these by powder x-ray diffraction (PXRD), nitrogen adsorption (BET) and high resolution transmission electron microscopy (HRTEM) analysis. The two routes give quite different materials although in both cases the products are essentially highly crystalline, dense particulates. It was found that the reverse precipitation technique gave the smallest crystallites with the narrowest size dispersion. This route also gave as-synthesised materials with higher surface areas. HRTEM confirmed the observations made from PXRD data and showed that the two methods resulted in quite different morphologies and surface chemistries. The forward route gives products with significantly greater densities of Ce3+ species compared to the reverse route. Data are explained using known precipitation chemistry and kinetic effects. Chapter five centres on the addition of terbia to ceria and has been investigated using XRD, XRF, XPS and TEM. Good solid solutions were formed across the entire composition range and there was no evidence for the formation of mixed phases or surface segregation over either the composition or temperature range investigated. Both Tb3+ and Tb4+ ions exist within the solution and the ratios of these cations are consistent with the addition of Tb8O15 to the fluorite ceria structure across a wide range of compositions. Local regions of anion vacancy ordering may be visible for small crystallites. There is no evidence of significant Ce3+ ion concentrations formed at the surface or in the bulk by the addition of terbia. The lattice parameter of these materials was seen to decrease with decreasing crystallite size. This is consistent with increased surface tension effects at small dimension. Chapter six reviews size related lattice parameter changes and surface defects in ceria nanocrystals. Ceria (CeO2) has many important applications, notably in catalysis. Many of its uses rely on generating nanodimensioned particles. Ceria has important redox chemistry where Ce4+ cations can be reversibly reduced to Ce3+ cations and associated anion vacancies. The significantly larger size of Ce3+ (compared with Ce4+) has been shown to result in lattice expansion. Many authors have observed lattice expansion in nanodimensioned crystals (nanocrystals), and these have been attributed to the presence of stabilized Ce3+ -anion vacancy combinations in these systems. Experimental results presented here show (i) that significant, but complex changes in the lattice parameter with size can occur in 2-500 nm crystallites, (ii) that there is a definitive relationship between defect chemistry and the lattice parameter in ceria nanocrystals, and (iii) that the stabilizing mechanism for the Ce3+ -anion vacancy defects at the surface of ceria nanocrystals is determined by the size, the surface status, and the analysis conditions. In this work, both lattice expansion and a more unusual lattice contraction in ultrafine nanocrystals are observed. The lattice deformations seen can be defined as a function of both the anion vacancy (hydroxyl) concentration in the nanocrystal and the intensity of the additional pressure imposed by the surface tension on the crystal. The expansion of lattice parameters in ceria nanocrystals is attributed to a number of factors, most notably, the presence of any hydroxyl moieties in the materials. Thus, a very careful understanding of the synthesis combined with characterization is required to understand the surface chemistry of ceria nanocrystals.

Synthesis and characterisation of novel superficially porous silica based stationary phases for use in liquid chromatography

The concept of pellicular particles was suggested by Horváth and Lipsky over fifty years ago. The reasoning behind the idea of these particles was to improve column efficiency by shortening the pathways analyte molecules can travel, therefore reducing the effect of the A and C terms. Several types of shell particles were successfully marketed around this time, however with the introduction of high quality fully porous silica under 10 μm, shell particles faded into the background. In recent years a new generation of core shell particles have become popular within the separation science community. These particles allow fast and efficient separations that can be carried out on conventional HPLC systems. Chapter 1 of this thesis introduces the chemistry of chromatographic stationary phases, with an emphasis on silica bonded phases, particularly focusing on the current state of technology in this area. The main focus is on superficially porous silica particles as a support material for liquid chromatography. A summary of the history and development of these particles over the past few decades is explored, along with current methods of synthesis of shell particles. While commercial shell particles have a rough outer surface, Chapter 2 focuses on the novel approach to growth of smooth surface superficially porous particles in a step-by-step manner. From the Stöber methodology to the seeded growth technique, and finally to the layer-bylayer growth of the porous shell. The superficially porous particles generated in this work have an overall diameter of 2.6 μm with a 350 nm porous shell; these silica particles were characterised using SEM, TEM and BET analysis. The uniform spherical nature of the particles along with their surface area, pore size and particle size distribution are examined in this chapter. I discovered that these smooth surface shell particles can be synthesised to give comparable surface area and pore size in comparison to commercial brands. Chapter 3 deals with the bonding of the particles prepared in Chapter 2 with C18 functionality; one with a narrow and one with a wide particle size distribution. This chapter examines the chromatographic and kinetic performance of these silica stationary phases, and compares them to a commercial superficially porous silica phase with a rough outer surface. I found that the particle size distribution does not seem to be the major contributor to the improvement in efficiency. The surface morphology of the particles appears to play an important role in the packing process of these particles and influences the Van Deemter effects. Chapter 4 focuses on the functionalisation of 2.6 μm smooth surface superficially porous particles with a variety of fluorinated and phenyl silanes. The same processes were carried out on 3.0 μm fully porous silica particles to provide a comparison. All phases were accessed using elemental analysis, thermogravimetric analysis, nitrogen sorption analysis and chromatographically evaluated using the Neue test. I observed comparable results for the 2.6 μm shell pentaflurophenyl propyl silica when compared to 3.0 μm fully porous silica. Chapter 5 moves towards nano-particles, with the synthesis of sub-1 μm superficially porous particles, their characterisation and use in chromatography. The particles prepared are 750 nm in total with a 100 nm shell. All reactions and testing carried out on these 750 nm core shell particles are also carried out on 1.5 μm fully porous particles in order to give a comparative result. The 750 nm core shell particles can be synthesised quickly and are very uniform. The main drawback in their use for HPLC is the system itself due to the backpressure experienced using sub – 1 μm particles. The synthesis of modified Stöber particles is also examined in this chapter with a range of non-porous silica and shell silica from 70 nm – 750 nm being tested for use on a Langmuir – Blodgett system. These smooth surface shell particles have only been in existence since 2009. The results displayed in this thesis demonstrate how much potential smooth surface shell particles have provided more in-depth optimisation is carried out. The results on packing studies reported in this thesis aims to be a starting point for a more sophisticated methodology, which in turn can lead to greater chromatographic improvements.

Novel structured emulsions for delivering of food flavours

Flavour release from food is determined by the binding of flavours to other food ingredients and the partition of flavour molecules among different phases. Food emulsions are used as delivery systems for food flavours, and tailored structuring in emulsions provides novel means to better control flavour release. The current study investigated four structured oil-in-water emulsions with structuring in the oil phase, oil-water interface, and water phase. Oil phase structuring was achieved by the formation of monoglyceride (MG) liquid crystals in the oil droplets (MG structured emulsions). Structured interface was created by the adsorption of a whey protein isolate (WPI)-pectin double layer at the interface (multilayer emulsion). Water phase structured emulsions referred to emulsion filled protein gels (EFP gels), where emulsion droplets were embedded in WPI gel network, and emulsions with maltodextrins (MDs) of different dextrose-equivalent (DE) values. Flavour compounds with different physicochemical properties were added into the emulsions, and flavour release (release rate, headspace concentration and air-emulsion partition coefficient) was described by GC headspace analysis. Emulsion structures, including crystalline structure, particle size, emulsion stability, rheology, texture, and microstructures, were characterized using differential scanning calorimetry and X-ray diffraction, light scattering, multisample analytical centrifuge, rheometry, texture analysis, and confocal laser scanning microscopy, respectively. In MG structured emulsions, MG self-assembled into liquid crystalline structures and stable β-form crystals were formed after 3 days of storage at 25 °C. The inclusion of MG crystals allowed tween 20 stabilized emulsions to present viscoelastic properties, and it made WPI stabilized emulsions more sensitive to the change of pH and NaCl concentrations. Flavour compounds in MG structured emulsions had lower initial headspace concentration and air-emulsion partition coefficients than those in unstructured emulsions. Flavour release can be modulated by changing MG content, oil content and oil type. WPI-pectin multilayer emulsions were stable at pH 5.0, 4.0, and 3.0, but they presented extensive creaming when subjected to salt solutions with NaCl ≥ 150 mM and mixed with artificial salivas. Increase of pH from 5.0 to 7.0 resulted in higher headspace concentration but unchanged release rate, and increase of NaCl concentration led to increased headspace concentration and release rate. The study also showed that salivas could trigger higher release of hydrophobic flavours and lower release of hydrophilic flavours. In EFP gels, increases in protein content and oil content contributed to gels with higher storage modulus and force at breaking. Flavour compounds had significantly reduced release rates and air-emulsion partition coefficients in the gels than the corresponding ungelled emulsions, and the reduction was in line with the increase of protein content. Gels with stronger gel network but lower oil content were prepared, and lower or unaffected release rates of the flavours were observed. In emulsions containing maltodextrins, water was frozen at a much lower temperature, and emulsion stability was greatly improved when subjected to freeze-thawing. Among different MDs, MD DE 6 offered the emulsion the highest stability. Flavours had lower air-emulsion partition coefficients in the emulsions with MDs than those in the emulsion without MD. Moreover, the involvement of MDs in the emulsions allowed most flavours had similar release profiles before and after freeze-thaw treatment. The present study provided information about different structured emulsions as delivery systems for flavour compounds, and on how food structure can be designed to modulate flavour release, which could be helpful in the development of functional foods with improved flavour profile.

Dry mixing of spice powders - investigation of effect of powder properties on mixture quality of binary powder mixtures

Dry mixing of binary food powders was conducted in a 2L lab-scale paddle mixer. Different types of food powders such as paprika, oregano, black pepper, onion powder and salt were used for the studies. A novel method based on a digital colour imaging system (DCI) was developed to measure the mixture quality (MQ) of binary food powder mixtures. The salt conductivity method was also used as an alternative method to measure the MQ. In the first part of the study the DCI method was developed and it showed potential for assessing MQ of binary powder mixes provided there was huge colour difference between the powders. In the second and third part of the study the effect of composition, water content, particle size and bulk density on MQ was studied. Flowability of powders at various moisture contents was also investigated. The mixing behaviour was assessed using coefficient of variation. Results showed that water content and composition influence the mixing behavior of powders. Good mixing was observed up to size ratios of 4.45 and at higher ratios MQ disimproved. The bulk density had a larger influence on the MQ. In the final study the MQ evaluation of binary and ternary powder mixtures was compared by using two methods – salt conductivity method and DCI method. Two binary food and two quaternary food powder mixtures with different coloured ingredients were studied. Overall results showed that DCI method has a potential for use by industries and it can analyse powder mixtures with components that have differences in colour and that are not segregating in nature.

Biomethane production from grass silage: laboratory assessment to maximise yields

On-farm biogas production is typically associated with forage maize as the biomass source. Digesters are designed and operated with the focus of optimising the conditions for this feedstock. Thus, such systems may not be ideally suited to the digestion of grass. Ireland has ca. 3.85 million ha of grassland. Annual excess grass, surplus to livestock requirements, could potentially fuel an anaerobic digestion industry. Biomethane associated with biomass from 1.1 % of grassland in Ireland, could potentially generate over 10 % renewable energy supply in transport. This study aims to identify and optimise technologies for the production of biomethane from grass silage. Mono-digestion of grass silage and co-digestion with slurry, as would occur on Irish farms, is investigated in laboratory trials. Grass silage was shown to have 7 times greater methane potential than dairy slurry on a fresh weight basis (107 m3 t-1 v 16 m3 t-1). However, comprehensive trace element profiles indicated that cobalt, iron and nickel are deficient in mono-digestion of grass silage at a high organic loading rate (OLR) of 4.0 kg VS m-3 d-1. The addition of a slurry co-substrate was beneficial due to its wealth of essential trace elements. To stimulate hydrolysis of high lignocellulose grass silage, particle size reduction (physical) and rumen fluid addition (biological) were investigated. In a continuous trial, digestion of grass silage of <1 cm particle size achieved a specific methane yield of 371 L CH4 kg-1 VS when coupled with rumen fluid addition. The concept of demand driven biogas was also examined in a two-phase digestion system (leaching with UASB). When demand for electricity is low it is recommended to disconnect the UASB from the system and recirculate rumen fluid to increase volatile fatty acid (VFA) and soluble chemical oxygen demand (SCOD) production whilst minimising volatile solids (VS) destruction. At times of high demand for electricity, connection of the UASB increases the destruction of volatiles and associated biogas production. The above experiments are intended to assess a range of biogas production options from grass silage with a specific focus on maximising methane yields and provide a guideline for feasible design and operation of on-farm digesters in Ireland.

Heterogeneities in fullerene nanoparticle aggregates affecting reactivity, bioactivity, and transport.

Properties of nanomaterial suspensions are typically summarized by average values for the purposes of characterizing these materials and interpreting experimental results. We show in this work that the heterogeneity in aqueous suspensions of fullerene C(60) aggregates (nC(60)) must be taken into account for the purposes of predicting nanomaterial transport, exposure, and biological activity. The production of reactive oxygen species (ROS), microbial inactivation, and the mobility of the aggregates of the nC(60) in a silicate porous medium all increased as suspensions were fractionated to enrich with smaller aggregates by progressive membrane filtration. These size-dependent differences are attributed to an increasing degree of hydroxylation of nC(60) aggregates with decreasing size. As the quantity and influence of these more reactive fractions may increase with time, experiments evaluating fullerene transport and toxicity end points must take into account the evolution and heterogeneity of fullerene suspensions.

The synchronization of superparamagnetic beads driven by a micro-magnetic ratchet.

We present theoretical, numerical, and experimental analyses on the non-linear dynamic behavior of superparamagnetic beads exposed to a periodic array of micro-magnets and an external rotating field. The agreement between theoretical and experimental results revealed that non-linear magnetic forcing dynamics are responsible for transitions between phase-locked orbits, sub-harmonic orbits, and closed orbits, representing different mobility regimes of colloidal beads. These results suggest that the non-linear behavior can be exploited to construct a novel colloidal separation device that can achieve effectively infinite separation resolution for different types of beads, by exploiting minor differences in their bead's properties. We also identify a unique set of initial conditions, which we denote the "devil's gate" which can be used to expeditiously identify the full range of mobility for a given bead type.

A method of assessing air toxics concentrations in urban areas using mobile platform measurements.

The objective of this paper is to demonstrate an approach to characterize the spatial variability in ambient air concentrations using mobile platform measurements. This approach may be useful for air toxics assessments in Environmental Justice applications, epidemiological studies, and environmental health risk assessments. In this study, we developed and applied a method to characterize air toxics concentrations in urban areas using results of the recently conducted field study in Wilmington, DE. Mobile measurements were collected over a 4- x 4-km area of downtown Wilmington for three components: formaldehyde (representative of volatile organic compounds and also photochemically reactive pollutants), aerosol size distribution (representing fine particulate matter), and water-soluble hexavalent chromium (representative of toxic metals). These measurements were,used to construct spatial and temporal distributions of air toxics in the area that show a very strong temporal variability, both diurnally and seasonally. An analysis of spatial variability indicates that all pollutants varied significantly by location, which suggests potential impact of local sources. From the comparison with measurements at the central monitoring site, we conclude that formaldehyde and fine particulates show a positive correlation with temperature, which could also be the reason that photochemically generated formaldehyde and fine particulates over the study area correlate well with the fine particulate matter measured at the central site.

Extraction of polydispersity information from photon correlation data

Photon correlation spectroscopy (PCS) is a light-scattering technique for particle size diagnosis. It has been used mainly in the investigation of hydrosol particles since it is based on the measurement of the correlation function of the light scattered from the Brownian motion of suspended particles. Recently this technique also proved useful for studying soot particles in flames and similar aerosol systems. In the case of a polydispersed system the problem of recovering the particle size distribution can be reduced to the problem of inverting the Laplace transform. In this paper we review several methods introduced by the authors for the solution of this problem. We present some numerical results and we discuss the resolution limits characterizing the reconstruction of the size distributions. © 1989.

Inversion of light scattering data by singular system analysis

We consider the problem of inverting experimental data obtained in light scattering experiments described by linear theories. We discuss applications to particle sizing and we describe fast and easy-to-implement algorithms which permit the extraction, from noisy measurements, of reliable information about the particle size distribution. © 1987, SPIE.

A comprehensive copper stockpile leach model: background and model formulation

The first phase in the sign, development and implementation of a comprehensive computational model of a copper stockpile leach process is presented. The model accounts for transport phenomena through the stockpile, reaction kinetics for the important mineral species, oxgen and bacterial effects on the leach reactions, plus heat, energy and acid balances for the overall leach process. The paper describes the formulation of the leach process model and its implementation in PHYSICA+, a computational fluid dynamic (CFD) software environment. The model draws on a number of phenomena to represent the competing physical and chemical features active in the process model. The phenomena are essentially represented by a three-phased (solid liquid gas) multi-component transport system; novel algorithms and procedures are required to solve the model equations, including a methodology for dealing with multiple chemical species with different reaction rates in ore represented by multiple particle size fractions. Some initial validation results and application simulations are shown to illustrate the potential of the model.

A comprehensive copper stockpile leach model: background and parameter sensitivity

The design and development of a comprehensive computational model of a copper stockpile leach process is summarized. The computational fluid dynamic software framework PHYSICA+ and various phenomena were used to model transport phenomena, mineral reaction kinetics, bacterial effects, and heat, energy and acid balances for the overall leach process. In this paper, the performance of the model is investigated, in particular its sensitvity to particle size and ore permeability. A combination of literature and laboratory sources was used to parameterize the model. The simulation results from the leach model are compared with closely controlled column pilot scale tests. The main performance characteristics (e.g. copper recovery rate) predicted by the model compare reasonably well with the experimental data and clearly reflect the qualitiative behavior of the process in many respects. The model is used to provide a measure of the sensitivity of ore permeability on leach behavior, and simulation results are examined for several different particle size distributions.

Analysis and modeling of heaping behavior of granular mixtures within a computational mechanics framework

This paper presents a continuum model of the flow of granular material during filling of a silo, using a viscoplastic constitutive relation based on the Drucker-Prager plasticity yield function. The performed simulations demonstrate the ability of the model to realistically represent complex features of granular flows during filling processes, such as heap formation and non-zero inclination angle of the bulk material-air interface. In addition, micro-mechanical parametrizations which account for particle size segregation are incorporated into the model. It is found that numerical predictions of segregation phenomena during filling of a binary granular mixture agree well with experimental results. Further numerical tests indicate the capability of the model to cope successfully with complex operations involving granular mixtures.