Optimisation of granola breakfast cereal manufacturing process by wet granulation and pneumatic conveying

This study has considered the optimisation of granola breakfast cereal manufacturing processes by wet granulation and pneumatic conveying. Granola is an aggregated food product used as a breakfast cereal and in cereal bars. Processing of granola involves mixing the dry ingredients (typically oats, nuts, etc.) followed by the addition of a binder which can contain honey, water and/or oil. In this work, the design and operation of two parallel wet granulation processes to produce aggregate granola products were incorporated: a) a high shear mixing granulation process followed by drying/toasting in an oven. b) a continuous fluidised bed followed by drying/toasting in an oven. In high shear granulation the influence of process parameters on key granule aggregate quality attributes such as granule size distribution and textural properties of granola were investigated. The experimental results show that the impeller rotational speed is the single most important process parameter which influences granola physical and textural properties. After that binder addition rate and wet massing time also show significant impacts on granule properties. Increasing the impeller speed and wet massing time increases the median granule size while also presenting a positive correlation with density. The combination of high impeller speed and low binder addition rate resulted in granules with the highest levels of hardness and crispness. In the fluidised bed granulation process the effect of nozzle air pressure and binder spray rate on key aggregate quality attributes were studied. The experimental results show that a decrease in nozzle air pressure leads to larger in mean granule size. The combination of lowest nozzle air pressure and lowest binder spray rate results in granules with the highest levels of hardness and crispness. Overall, the high shear granulation process led to larger, denser, less porous and stronger (less likely to break) aggregates than the fluidised bed process. The study also examined the particle breakage of granola during pneumatic conveying produced by both the high shear granulation and the fluidised bed granulation process. Products were pneumatically conveyed in a purpose built conveying rig designed to mimic product conveying and packaging. Three different conveying rig configurations were employed; a straight pipe, a rig consisting two 45° bends and one with 90° bend. Particle breakage increases with applied pressure drop, and a 90° bend pipe results in more attrition for all conveying velocities relative to other pipe geometry. Additionally for the granules produced in the high shear granulator; those produced at the highest impeller speed, while being the largest also have the lowest levels of proportional breakage while smaller granules produced at the lowest impeller speed have the highest levels of breakage. This effect clearly shows the importance of shear history (during granule production) on breakage during subsequent processing. In terms of the fluidised bed granulation, there was no single operating parameter that was deemed to have a significant effect on breakage during subsequent conveying. Finally, a simple power law breakage model based on process input parameters was developed for both manufacturing processes. It was found suitable for predicting the breakage of granola breakfast cereal at various applied air velocities using a number of pipe configurations, taking into account shear histories.

Climate and energy scenarios for Ireland to 2050 using the Irish TIMES energy systems model

Due to growing concerns regarding the anthropogenic interference with the climate system, countries across the world are being challenged to develop effective strategies to mitigate climate change by reducing or preventing greenhouse gas (GHG) emissions. The European Union (EU) is committed to contribute to this challenge by setting a number of climate and energy targets for the years 2020, 2030 and 2050 and then agreeing effort sharing amongst Member States. This thesis focus on one Member State, Ireland, which faces specific challenges and is not on track to meet the targets agreed to date. Before this work commenced, there were no projections of energy demand or supply for Ireland beyond 2020. This thesis uses techno-economic energy modelling instruments to address this knowledge gap. It builds and compares robust, comprehensive policy scenarios, providing a means of assessing the implications of different future energy and emissions pathways for the Irish economy, Ireland’s energy mix and the environment. A central focus of this thesis is to explore the dynamics of the energy system moving towards a low carbon economy. This thesis develops an energy systems model (the Irish TIMES model) to assess the implications of a range of energy and climate policy targets and target years. The thesis also compares the results generated from the least cost scenarios with official projections and target pathways and provides useful metrics and indications to identify key drivers and to support both policy makers and stakeholder in identifying cost optimal strategies. The thesis also extends the functionality of energy system modelling by developing and applying new methodologies to provide additional insights with a focus on particular issues that emerge from the scenario analysis carried out. Firstly, the thesis develops a methodology for soft-linking an energy systems model (Irish TIMES) with a power systems model (PLEXOS) to improve the interpretation of the electricity sector results in the energy system model. The soft-linking enables higher temporal resolution and improved characterisation of power plants and power system operation Secondly, the thesis develops a methodology for the integration of agriculture and energy systems modelling to enable coherent economy wide climate mitigation scenario analysis. This provides a very useful starting point for considering the trade-offs between the energy system and agriculture in the context of a low carbon economy and for enabling analysis of land-use competition. Three specific time scale perspectives are examined in this thesis (2020, 2030, 2050), aligning with key policy target time horizons. The results indicate that Ireland’s short term mandatory emissions reduction target will not be achieved without a significant reassessment of renewable energy policy and that the current dominant policy focus on wind-generated electricity is misplaced. In the medium to long term, the results suggest that energy efficiency is the first cost effective measure to deliver emissions reduction; biomass and biofuels are likely to be the most significant fuel source for Ireland in the context of a low carbon future prompting the need for a detailed assessment of possible implications for sustainability and competition with the agri-food sectors; significant changes are required in infrastructure to deliver deep emissions reductions (to enable the electrification of heat and transport, to accommodate carbon capture and storage facilities (CCS) and for biofuels); competition between energy and agriculture for land-use will become a key issue. The purpose of this thesis is to increase the evidence-based underpinning energy and climate policy decisions in Ireland. The methodology is replicable in other Member States.

Analysis of soil carbon transit times and age distributions using network theories

The long-term soil carbon dynamics may be approximated by networks of linear compartments, permitting theoretical analysis of transit time (i.e., the total time spent by a molecule in the system) and age (the time elapsed since the molecule entered the system) distributions. We compute and compare these distributions for different network. configurations, ranging from the simple individual compartment, to series and parallel linear compartments, feedback systems, and models assuming a continuous distribution of decay constants. We also derive the transit time and age distributions of some complex, widely used soil carbon models (the compartmental models CENTURY and Rothamsted, and the continuous-quality Q-Model), and discuss them in the context of long-term carbon sequestration in soils. We show how complex models including feedback loops and slow compartments have distributions with heavier tails than simpler models. Power law tails emerge when using continuous-quality models, indicating long retention times for an important fraction of soil carbon. The responsiveness of the soil system to changes in decay constants due to altered climatic conditions or plant species composition is found to be stronger when all compartments respond equally to the environmental change, and when the slower compartments are more sensitive than the faster ones or lose more carbon through microbial respiration. Copyright 2009 by the American Geophysical Union.

Scale-wise evolution of rainfall probability density functions fingerprints the rainfall generation mechanism

© 2010 by the American Geophysical Union.The cross-scale probabilistic structure of rainfall intensity records collected over time scales ranging from hours to decades at sites dominated by both convective and frontal systems is investigated. Across these sites, intermittency build-up from slow to fast time-scales is analyzed in terms of heavy tailed and asymmetric signatures in the scale-wise evolution of rainfall probability density functions (pdfs). The analysis demonstrates that rainfall records dominated by convective storms develop heavier-Tailed power law pdfs toward finer scales when compared with their frontal systems counterpart. Also, a concomitant marked asymmetry build-up emerges at such finer time scales. A scale-dependent probabilistic description of such fat tails and asymmetry appearance is proposed based on a modified q-Gaussian model, able to describe the cross-scale rainfall pdfs in terms of the nonextensivity parameter q, a lacunarity (intermittency) correction and a tail asymmetry coefficient, linked to the rainfall generation mechanism.

Universal Quake Statistics: From Compressed Nanocrystals to Earthquakes.

Slowly-compressed single crystals, bulk metallic glasses (BMGs), rocks, granular materials, and the earth all deform via intermittent slips or "quakes". We find that although these systems span 12 decades in length scale, they all show the same scaling behavior for their slip size distributions and other statistical properties. Remarkably, the size distributions follow the same power law multiplied with the same exponential cutoff. The cutoff grows with applied force for materials spanning length scales from nanometers to kilometers. The tuneability of the cutoff with stress reflects "tuned critical" behavior, rather than self-organized criticality (SOC), which would imply stress-independence. A simple mean field model for avalanches of slipping weak spots explains the agreement across scales. It predicts the observed slip-size distributions and the observed stress-dependent cutoff function. The results enable extrapolations from one scale to another, and from one force to another, across different materials and structures, from nanocrystals to earthquakes.

Scaling laws of ambush predator 'waiting' behaviour are tuned to a common ecology

The decisions animals make about how long to wait between activities can determine the success of diverse behaviours such as foraging, group formation or risk avoidance. Remarkably, for diverse animal species, including humans, spontaneous patterns of waiting times show random ‘burstiness’ that appears scale-invariant across a broad set of scales. However, a general theory linking this phenomenon across the animal kingdom currently lacks an ecological basis. Here, we demonstrate from tracking the activities of 15 sympatric predator species (cephalopods, sharks, skates and teleosts) under natural and controlled conditions that bursty waiting times are an intrinsic spontaneous behaviour well approximated by heavy-tailed (power-law) models over data ranges up to four orders of magnitude. Scaling exponents quantifying ratios of frequent short to rare very long waits are species-specific, being determined by traits such as foraging mode (active versus ambush predation), body size and prey preference. A stochastic–deterministic decision model reproduced the empirical waiting time scaling and species-specific exponents, indicating that apparently complex scaling can emerge from simple decisions. Results indicate temporal power-law scaling is a behavioural ‘rule of thumb’ that is tuned to species’ ecological traits, implying a common pattern may have naturally evolved that optimizes move–wait decisions in less predictable natural environments.

Examination of the flow rheological and textural properties of polymer gels composed of poly(methylvinylether-co-maleic anhydride) and poly(vinylpyrrolidone): rheological and mathematical interpretation of textural parameters.

The purpose of this study was to mathematically characterize the effects of defined experimental parameters (probe speed and the ratio of the probe diameter to the diameter of sample container) on the textural/mechanical properties of model gel systems. In addition, this study examined the applicability of dimensional analysis for the rheological interpretation of textural data in terms of shear stress and rate of shear. Aqueous gels (pH 7) were prepared containing 15% w/w poly(methylvinylether-co-maleic anhydride) and poly(vinylpyrrolidone) (PVP) (0, 3, 6, or 9% w/w). Texture profile analysis (TPA) was performed using a Stable Micro Systems texture analyzer (model TA-XT 2; Surrey, UK) in which an analytical probe was twice compressed into each formulation to a defined depth (15 mm) and at defined rates (1, 3, 5, 8, and 10 mm s-1), allowing a delay period (15 s) between the end of the first and beginning of the second compressions. Flow rheograms were performed using a Carri-Med CSL2-100 rheometer (TA Instruments, Surrey, UK) with parallel plate geometry under controlled shearing stresses at 20.0°?±?0.1°C. All formulations exhibited pseudoplastic flow with no thixotropy. Increasing concentrations of PVP significantly increased formulation hardness, compressibility, adhesiveness, and consistency. Increased hardness, compressibility, and consistency were ascribed to enhanced polymeric entanglements, thereby increasing the resistance to deformation. Increasing probe speed increased formulation hardness in a linear manner, because of the effects of probe speed on probe displacement and surface area. The relationship between formulation hardness and probe displacement was linear and was dependent on probe speed. Furthermore, the proportionality constant (gel strength) increased as a function of PVP concentration. The relationship between formulation hardness and diameter ratio was biphasic and was statistically defined by two linear relationships relating to diameter ratios from 0 to 0.4 and from 0.4 to 0.563. The dramatically increased hardness, associated with diameter ratios in excess of 0.4, was accredited to boundary effects, that is, the effect of the container wall on product flow. Using dimensional analysis, the hardness and probe displacement in TPA were mathematically transformed into corresponding rheological parameters, namely shearing stress and rate of shear, thereby allowing the application of the power law (??=?k?n) to textural data. Importantly, the consistencies (k) of the formulations, calculated using transformed textural data, were statistically similar to those obtained using flow rheometry. In conclusion, this study has, firstly, characterized the relationships between textural data and two key instrumental parameters in TPA and, secondly, described a method by which rheological information may be derived using this technique. This will enable a greater application of TPA for the rheological characterization of pharmaceutical gels and, in addition, will enable efficient interpretation of textural data under different experimental parameters.

Single-ionization of helium at Ti:sapphire wavelengths: rates and scaling laws

We present a numerical and theoretical study of intense-field single-electron ionization of helium at 390 nm and 780 nm. Accurate ionization rates (over an intensity range of (0.175-34) X10^14 W/ cm^2 at 390 nm, and (0.275 - 14.4) X 10^14 W /cm^2 at 780 nm) are obtained from full-dimensionality integrations of the time-dependent helium-laser Schroedinger equation. We show that the power law of lowest order perturbation theory, modified with a ponderomotive-shifted ionization potential, is capable of modelling the ionization rates over an intensity range that extends up to two orders of magnitude higher than that applicable to perturbation theory alone. Writing the modified perturbation theory in terms of scaled wavelength and intensity variables, we obtain to first approximation a single ionization law for both the 390 nm and 780 nm cases. To model the data in the high intensity limit as well as in the low, a new function is introduced for the rate. This function has, in part, a resemblance to that derived from tunnelling theory but, importantly, retains the correct frequency-dependence and scaling behaviour derived from the perturbative-like models at lower intensities. Comparison with the predictions of classical ADK tunnelling theory confirms that ADK performs poorly in the frequency and intensity domain treated here.

The 2003 November 14 occultation by Titan of TYC 1343-1865-1

We observed a stellar occultation by Titan on 2003 November 14 from La Palma Observatory using ULTRACAM with three Sloan filters: u, g, and i (358, 487, and 758 nm, respectively). The occultation probed latitudes 2°?S and 1°?N during immersion and emersion, respectively. A prominent central flash was present in only the i filter, indicating wavelength-dependent atmospheric extinction. We inverted the light curves to obtain six lower-limit temperature profiles between 335 and 485 km (0.04 and 0.003 mb) altitude. The i profiles agreed with the temperature measured by the Huygens Atmospheric Structure Instrument [Fulchignoni, M., and 43 colleagues, 2005. Nature 438, 785 791] above 415 km (0.01 mb). The profiles obtained from different wavelength filters systematically diverge as altitude decreases, which implies significant extinction in the light curves. Applying an extinction model [Elliot, J.L., Young, L.A., 1992. Astron. J. 103, 991 1015] gave the altitudes of line of sight optical depth equal to unity: 396±7 and 401±20 km (u immersion and emersion); 354±7 and 387±7 km (g immersion and emersion); and 336±5 and 318±4 km (i immersion and emersion). Further analysis showed that the optical depth follows a power law in wavelength with index 1.3±0.2. We present a new method for determining temperature from scintillation spikes in the occulting body's atmosphere. Temperatures derived with this method are equal to or warmer than those measured by the Huygens Atmospheric Structure Instrument. Using the highly structured, three-peaked central flash, we confirmed the shape of Titan's middle atmosphere using a model originally derived for a previous Titan occultation [Hubbard, W.B., and 45 colleagues, 1993. Astron. Astrophys. 269, 541 563].

NUMERICAL SIMULATIONS OF BURST PROCESSES IN FIBER-BUNDLES

We discuss a very effective numerical method for simulating fibre-bundle models with equal load-sharing and with local load-sharing. Particular attention is paid to the case of the local load-sharing model, in which the critical load x(c) is defined as the average load per fibre that causes the final complete failure. It is shown that x(c) --> 0 when the size of the system N --> infinity. We also show analytically that the power law of the burst size distribution, D(Delta) alpha Delta(-xi), is approximately correct. The exponent xi in the local load-sharing case is not universal, since it depends on the strength distribution as well on as the size of the system.

CRITICAL PROCESSES, LANGEVIN EQUATION AND UNIVERSALITY

In nature there are ubiquitous systems that can naturally approach critical states, The Langevin equation in the discrete version can be used to describe a class of critical processes, which are characterized by power-law behaviors and scaling relations. As an example, we present a simple model for a clinical thermometer, whose reading cannot fall even when its temperature decreases. The fibers bundle model and the spring-block model are also shown to belong to such a class.

BURST-SIZE DISTRIBUTION IN FIBER-BUNDLES WITH LOCAL LOAD-SHARING

A critical load x(c) is introduced into the fiber-bundle model with local load-sharing. The critical load is defined as the average load per fiber that causes the final complete failure. It is shown that x(c) --> 0 when the size of the system N --> infinity. A power law for the burst-size distribution, D(DELTA) is-proportional-to DELTA(-xi) is approximately correct. The exponent xi is not universal, since it depends on the strength distribution as well as the size of the system.

Beyond the ideal free distribution: More general models of predator distribution

The ideal free distribution model which relates the spatial distribution of mobile consumers to that of their resource is shown to be a limiting case of a more general model which we develop using simple concepts of diffusion. We show how the ideal free distribution model can be derived from a more general model and extended by incorporating simple models of social influences on predator spacing. First, a free distribution model based on patch switching rules, with a power-law interference term, which represents instantaneous biased diffusion is derived. A social bias term is then introduced to represent the effect of predator aggregation on predator fitness, separate from any effects which act through intake rate. The social bias term is expanded to express an optimum spacing for predators and example solutions of the resulting biased diffusion models are shown. The model demonstrates how an empirical interference coefficient, derived from measurements of predator and prey densities, may include factors expressing the impact of social spacing behaviour on fitness. We conclude that empirical values of log predator/log prey ratio may contain information about more than the relationship between consumer and resource densities. Unlike many previous models, the model shown here applies to conditions without continual input. (C) 1997 Academic Press Limited.</p>

Are power laws that estimate fractal dimension a good descriptor of soil structure and its link to soil biological properties?

The study of interrelationships between soil structure and its functional properties is complicated by the fact that the quantitative description of soil structure is challenging. Soil scientists have tackled this challenge by taking advantage of approaches such as fractal geometry, which describes soil architectural complexity through a scaling exponent (D) relating mass and numbers of particles/aggregates to particle/aggregate size. Typically, soil biologists use empirical indices such as mean weight diameters (MWD) and percent of water stable aggregates (WSA), or the entire size distribution, and they have successfully related these indices to key soil features such as C and N dynamics and biological promoters of soil structure. Here, we focused on D, WSA and MWD and we tested whether: D estimated by the exponent of the power law of number-size distributions is a good and consistent correlate of MWD and WSA; D carries information that differs from MWD and WSA; the fraction of variation in D that is uncorrelated with MWD and WSA is related to soil chemical and biological properties that are thought to establish interdependence with soil structure (e.g., organic C, N, arbuscular mycorrhizal fungi). We analysed observational data from a broad scale field study and results from a greenhouse experiment where arbuscular mycorrhizal fungi (AMF) and collembola altered soil structure. We were able to develop empirical models that account for a highly significant and large portion of the correlation observed between WSA and MWD but we did not uncover the mechanisms that underlie this correlation. We conclude that most of the covariance between D and soil biotic (AMF, plant roots) and abiotic (C. N) properties can be accounted for by WSA and MWD. This result implies that the ecological effects of the fragmentation properties described by D and generally discussed under the framework of fractal models can be interpreted under the intuitive perspective of simpler indices and we suggest that the biotic components mostly impacted the largest size fractions, which dominate MWD, WSA and the scaling exponent ruling number-size distributions. (C) 2010 Elsevier Ltd. All rights reserved.

GIANT LOBES OF CENTAURUS A RADIO GALAXY OBSERVED WITH THE SUZAKU X-RAY SATELLITE

We report on Suzaku observations of selected regions within the southern giant lobe of the radio galaxy Centaurus A. In our analysis we focus on distinct X-ray features detected with the X-ray Imaging Spectrometer within the range 0.5-10 keV, some of which are likely associated with fine structure of the lobe revealed by recent high-quality radio intensity and polarization maps. With the available photon statistics, we find that the spectral properties of the detected X-ray features are equally consistent with thermal emission from hot gas with temperatures kT > 1 keV, or with a power-law radiation continuum characterized by photon indices Gamma similar to 2.0 +/- 0.5. However, the plasma parameters implied by these different models favor a synchrotron origin for the analyzed X-ray spots, indicating that a very efficient acceleration of electrons up to greater than or similar to 10 TeV energies is taking place within the giant structure of Centaurus A, albeit only in isolated and compact regions associated with extended and highly polarized radio filaments. We also present a detailed analysis of the diffuse X-ray emission filling the whole field of view of the instrument, resulting in a tentative detection of a soft excess component best fitted by a thermal model with a temperature of kT similar to 0.5 keV. The exact origin of the observed excess remains uncertain, although energetic considerations point to thermal gas filling the bulk of the volume of the lobe and mixed with the non-thermal plasma, rather than to the alternative scenario involving a condensation of the hot intergalactic medium around the edges of the expanding radio structure. If correct, this would be the first detection of the thermal content of the extended lobes of a radio galaxy in X-rays. The corresponding number density of the thermal gas in such a case is n(g) similar to 10(-4) cm(-3), while its pressure appears to be in almost exact equipartition with the volume-averaged non-thermal pressure provided by the radio-emitting electrons and the lobes' magnetic field. A prominent large-scale fluctuation of the Galactic foreground emission, resulting in excess foreground X-ray emission aligned with the lobe, cannot be ruled out. Although tentative, our findings potentially imply that the structure of the extended lobes in active galaxies is likely to be highly inhomogeneous and non-uniform, with magnetic reconnection and turbulent acceleration processes continuously converting magnetic energy to internal energy of the plasma particles, leading to possibly significant spatial and temporal variations in the plasma beta parameter around the volume-averaged equilibrium condition beta similar to 1.