Modelling conditional probability distributions for periodic variables

Most of the common techniques for estimating conditional probability densities are inappropriate for applications involving periodic variables. In this paper we apply two novel techniques to the problem of extracting the distribution of wind vector directions from radar catterometer data gathered by a remote-sensing satellite.

Bayesian invariant measurements of generalisation for discrete distributions

Neural network learning rules can be viewed as statistical estimators. They should be studied in Bayesian framework even if they are not Bayesian estimators. Generalisation should be measured by the divergence between the true distribution and the estimated distribution. Information divergences are invariant measurements of the divergence between two distributions. The posterior average information divergence is used to measure the generalisation ability of a network. The optimal estimators for multinomial distributions with Dirichlet priors are studied in detail. This confirms that the definition is compatible with intuition. The results also show that many commonly used methods can be put under this unified framework, by assume special priors and special divergences.

CFD modelling of the fast pyrolysis of an in-flight cellulosic particle subjected to convective heat transfer

The pyrolysis of a freely moving cellulosic particle inside a 41.7mgs -1 source continuously fed fluid bed reactor subjected to convective heat transfer is modelled. The Lagrangian approach is adopted for the particle tracking inside the reactor, while the flow of the inert gas is treated with the standard Eulerian method for gases. The model incorporates the thermal degradation of cellulose to char with simultaneous evolution of gases and vapours from discrete cellulosic particles. The reaction kinetics is represented according to the Broido–Shafizadeh scheme. The convective heat transfer to the surface of the particle is solved by two means, namely the Ranz–Marshall correlation and the limit case of infinitely fast external heat transfer rates. The results from both approaches are compared and discussed. The effect of the different heat transfer rates on the discrete phase trajectory is also considered.

Computational modelling of the impact of particle size to the heat transfer coefficient between biomass particles and a fluidised bed

The fluid–particle interaction and the impact of different heat transfer conditions on pyrolysis of biomass inside a 150 g/h fluidised bed reactor are modelled. Two different size biomass particles (350 µm and 550 µm in diameter) are injected into the fluidised bed. The different biomass particle sizes result in different heat transfer conditions. This is due to the fact that the 350 µm diameter particle is smaller than the sand particles of the reactor (440 µm), while the 550 µm one is larger. The bed-to-particle heat transfer for both cases is calculated according to the literature. Conductive heat transfer is assumed for the larger biomass particle (550 µm) inside the bed, while biomass–sand contacts for the smaller biomass particle (350 µm) were considered unimportant. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Biomass reaction kinetics is modelled according to the literature using a two-stage, semi-global model which takes into account secondary reactions. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of User Defined Function (UDF).

Sequence distributions in free-radical polymers

This thesis is concerned with demonstrating how the visual representation of the sequence distribution of individual monomer units, of a polymer, that would be observed upon polymerisation, may be utilised in designing and synthesizing polymers with relatively low cell adhesion characteristics, The initial part of this thesis is concerned with demonstrating the use of a computer simulation technique, in illustrating the sequence distribution that would be observed upon the polymerisation of a set of monomers. The power of the computer simulation technique has been demonstrated through the simulation of the sequence distributions of some generic contact lens materials. These generic contact lens materials were chosen simply because in the field of biomaterials their compositions are amongst the most systematically regulated and they present a wide range of compositions. The validity of the computer simulation technique has been assessed through the synthesis and analysis of linear free-radical polymers at different conversions. Two main parameters were examined, that of composition and the number-average sequence lengths of individual monomer units, at various conversions. The polymers were synthesized through the solution polymerisation process. The monomer composition was determined by elemental analysis and 13C nuclear magnetic analysis (NMR). Number-average sequence lengths were determined exclusively through 13C NMR. Although the computer simulation technique provides a visual representation of the monomer sequence distribution up to 100% conversion, these assessments were made on linear polymers at a reasonably high conversion (above 50%) but below 100% conversion of ease for analysis. The analyses proved that the computer simulation technique was reasonably accurate in predicting the sequence distribution of monomer units, upon polymerisation, in the polymer.An approach has been presented which allows one to manipulate the use of monomers, with their reactivity ratios, thereby enabling us to design polymers with controlled sequence distributions.Hydrogel membranes, with relatively controlled sequence distributions and polymerised to 100% conversion, were synthesized to represent prospective biomaterials. Cell adhesion studies were used as a biological probe to investigate the susceptibility of the surface of these membranes to cell adhesion. This was necessary in order to assess the surface biocompatibility or biotolerance of these prospective biomaterials.

Hydrodynamic modeling of mineral wool fiber suspensions in a two-dimensional channel flow

A consequence of a loss of coolant accident is the damage of adjacent insulation materials (IM). IM may then be transported to the containment sump strainers where water is drawn into the ECCS (emergency core cooling system). Blockage of the strainers by IM lead to an increased pressure drop acting on the operating ECCS pumps. IM can also penetrate the strainers, enter the reactor coolant system and then accumulate in the reactor pressure vessel. An experimental and theoretical study that concentrates on mineral wool fiber transport in the containment sump and the ECCS is being performed. The study entails fiber generation and the assessment of fiber transport in single and multi-effect experiments. The experiments include measurement of the terminal settling velocity, the strainer pressure drop, fiber sedimentation and resuspension in a channel flow and jet flow in a rectangular tank. An integrated test facility is also operated to assess the compounded effects. Each experimental facility is used to provide data for the validation of equivalent computational fluid dynamic models. The channel flow facility allows the determination of the steady state distribution of the fibers at different flow velocities. The fibers are modeled in the Eulerian-Eulerian reference frame as spherical wetted agglomerates. The fiber agglomerate size, density, the relative viscosity of the fluid-fiber mixture and the turbulent dispersion of the fibers all affect the steady state accumulation of fibers at the channel base. In the current simulations, two fiber phases are separately considered. The particle size is kept constant while the density is modified, which affects both the terminal velocity and volume fraction. The relative viscosity is only significant at higher concentrations. The numerical model finds that the fibers accumulate at the channel base even at high velocities; therefore, modifications to the drag and turbulent dispersion forces can be made to reduce fiber accumulation.

The separation of secondary oil-water dispersions in particulate beds

The literature relating to haze formation, methods of separation, coalescence mechanisms, and models by which droplets <100 μm are collected, coalesced and transferred, have been reviewed with particular reference to particulate bed coalescers. The separation of secondary oil-water dispersions was studied experimentally using packed beds of monosized glass ballotini particles. The variables investigated were superficial velocity, bed depth, particle size, and the phase ratio and drop size distribution of inlet secondary dispersion. A modified pump loop was used to generate secondary dispersions of toluene or Clairsol 350 in water with phase ratios between 0.5-6.0 v/v%.Inlet drop size distributions were determined using a Malvern Particle Size Analyser;effluent, coalesced droplets were sized by photography. Single phase flow pressure drop data were correlated by means of a Carman-Kozeny type equation. Correlations were obtained relating single and two phase pressure drops, as (ΔP2/μc)/ΔP1/μd) = kp Ua Lb dcc dpd Cine A flow equation was derived to correlate the two phase pressure drop data as, ΔP2/(ρcU2) = 8.64*107 [dc/D]-0.27 [L/D]0.71 [dp/D]-0.17 [NRe]1.5 [e1]-0.14 [Cin]0.26  In a comparison between functions to characterise the inlet drop size distributions a modification of the Weibull function provided the best fit of experimental data. The general mean drop diameter was correlated by: q_p q_p p_q /β      Γ ((q-3/β) +1) d qp = d fr  .α        Γ ((P-3/β +1 The measured and predicted mean inlet drop diameters agreed within ±15%. Secondary dispersion separation depends largely upon drop capture within a bed. A theoretical analysis of drop capture mechanisms in this work indicated that indirect interception and London-van der Waal's mechanisms predominate. Mathematical models of dispersed phase concentration m the bed were developed by considering drop motion to be analogous to molecular diffusion.The number of possible channels in a bed was predicted from a model in which the pores comprised randomly-interconnected passage-ways between adjacent packing elements and axial flow occured in cylinders on an equilateral triangular pitch. An expression was derived for length of service channels in a queuing system leading to the prediction of filter coefficients. The insight provided into the mechanisms of drop collection and travel, and the correlations of operating parameters, should assist design of industrial particulate bed coalescers.

The effect of catchment characteristics on sewage settling velocity grading

A procedure has been developed which measures the settling velocity distribution of particles within a complete sewage sample. The development of the test method included observations of particle and liquid interaction using both synthetic media and sewage. Comparison studies with two other currently used settling velocity test procedures was undertaken. The method is suitable for use with either DWF or storm sewage. Information relating to the catchment characteristics of 35 No. wastewater treatment works was collected from the privatised water companies in England and Wales. 29 No. of these sites were used in an experimental programme to determine the settling velocity grading of 33 No. sewage samples. The collected data were analysed in an attempt to relate the settling velocity distribution to the characteristics of the contributing catchment. Statistical analysis of the catchment data and the measured settling velocity distributions was undertaken. A curve fitting exercise using an S-shaped curve which had the same physical characteristics as the settling velocity distributions was performed. None of these analyses found evidence that the settling velocity distribution of sewage had a significant relationship with the chosen catchment characteristics. The regression equations produced from the statistical analysis cannot be used to assist in the design of separation devices. However, a grading curve envelope was produced, the limits of which were clearly defined for the measured data set. There was no evidence of a relationship between settling velocity grading and the characteristics of the contributing catchment, particularly the catchment area. The present empirical approach to settling tank design cannot be improved upon at present by considering the variation in catchment parameters. This study has provided a basis for future research into the settling velocity measurement and should be of benefit to future workers within this field.

A study of the micro-mechanics of granular media

The development of more realistic constitutive models for granular media, such as sand, requires ingredients which take into account the internal micro-mechanical response to deformation. Unfortunately, at present, very little is known about these mechanisms and therefore it is instructive to find out more about the internal nature of granular samples by conducting suitable tests. In contrast to physical testing the method of investigation used in this study employs the Distinct Element Method. This is a computer based, iterative, time-dependent technique that allows the deformation of granular assemblies to be numerically simulated. By making assumptions regarding contact stiffnesses each individual contact force can be measured and by resolution particle centroid forces can be calculated. Then by dividing particle forces by their respective mass, particle centroid velocities and displacements are obtained by numerical integration. The Distinct Element Method is incorporated into a computer program 'Ball'. This program is effectively a numerical apparatus which forms a logical housing for this method and allows data input and output, and also provides testing control. By using this numerical apparatus tests have been carried out on disc assemblies and many new interesting observations regarding the micromechanical behaviour are revealed. In order to relate the observed microscopic mechanisms of deformation to the flow of the granular system two separate approaches have been used. Firstly a constitutive model has been developed which describes the yield function, flow rule and translation rule for regular assemblies of spheres and discs when subjected to coaxial deformation. Secondly statistical analyses have been carried out using data which was extracted from the simulation tests. These analyses define and quantify granular structure and then show how the force and velocity distributions use the structure to produce the corresponding stress and strain-rate tensors.

HTHP syngas cleaning concept of two stage biomass gasification for FT synthesis

The EU intends to increase the fraction of fuels from biogenic energy sources from 2% in 2005 to 8% in 2020. This means a minimum of 30 million TOE/a of fuels from biomass. This makes technical-scale generation of syngas from high-grade biomass, e.g. straw, hay, bark, or paper/cardboard waste, and the production of synthetic fuels by Fischer-Tropsch (FT) synthesis highly attractive. The BTL concept (Biomass to Liquids) of the Karlsruhe Research Center, labeled bioliq, focuses on this challenge by locally concentrating the biomass energy content by fast pyrolysis in a coke/oil slurry followed by slurry conversion to syngas in a central entrained flow gasifier at 1200C and pressures above 4MPa. FT synthesis generates intermediate products for synthetic fuels. To prevent the sensitive catalysts from being poisoned the syngas must be free of tar and particulates. Trace concentrations of H2S, COS, CS2, HCl, NH3, and HCN must be on the order of a few ppb. Moreover, maximum conversion efficiency will be achieved by cleaning the gas above the synthesis conditions. (T>350C, P>4MPa). The concept of an innovative dry HTHP syngas cleaning process is presented. Based on HT particle filtration and suitable sorption and catalysis processes for the relevant contaminants, an overall concept will be derived, which leads to a syngas quality required for FT synthesis in only two combined stages. Results of filtration experiments on a pilot scale are presented. The influence of temperature on the separation and conversion, respectively, of particulates and gaseous contaminants is discussed on the basis of experimental results obtained on a laboratory and pilot scale. Extensive studies of this concept are performed in a scientific network comprising the Karlsruhe Research Center and five universities; funding is provided by the Helmholtz Association of National Research Centers in Germany.

Preparation of particle-stabilized emulsions using membrane emulsification

There is currently significant interest in particle-stabilized emulsions for a variety of applications and as precursors to other materials such as micro-capsules or colloidosomes. A prerequisite for many applications is the ability to produce stable droplets with a well-controlled size. The preparation of oil-in-water (o/w) emulsions stabilized by silica colloids has been demonstrated here using membrane ulsification techniques. Emulsions were produced using both a cross-flow membrane device and a rotating membrane reactor. Under the correct conditions, highly stable emulsions with very narrow droplet size distributions can be produced. Investigations into the effects of changing the cross-flow shear rate at a fixed droplet production rate illustrate the fine control over mean droplet size that is possible with these emulsification techniques. Evidence for the importance of particle adsorption kinetics onto growing droplets prior to detachment from the membrane surface was obtained by varying the droplet production rate under fixed shear conditions. The presence of a critical surface coverage by the stabilizing particles to prevent droplet coalescence was clearly seen. Comparison with samples produced using conventional high-shear homogenization highlights the improved control over size distribution available from these membrane techniques.

A study of the properties and applications of electrostatic charged particle oscillators

This investigation originated from work by Dr. A.H. McIlraith of the National Physical Laboratory who, in 1966, described a new type of charged particle oscillator. This makes use of two equal cylindrical electrodes to constrain the particles in such a way that they follow extremely long oscillatory paths between the electrodes under the influence of an electrostatic field alone. The object of this work has been to study the principle of the oscillator in detail and to investigate its properties and applications. Any device which is capable of creating long electron trajectories has potential application in the field of ultra high vacuum technology. It was therefore considered that a critical review of the problems associated with the production and measurement of ultra high vacuum was relevant in the initial stages of the work. The oscillator has been applied with a considerable degree of success as a high energy electrostatic ion source. This offers several advantages over existing ion sources. It can be operated at much lower pressures without the need of a magnetic field. The oscillator principle has also been applied as a thermionic ionization gauge and has been compared with other ionization gauges to pressures as low as 5 x 10- 11 torr.. This new gauge exhibited a number of advantages over most of the existing gauges. Finally the oscillator has been used in an evaporation ion pump and has exhibited fairly high pumping speeds for argon gas relative to those for nitrogen. This investigation supports the original work of Dr. A.H. McIlraith and shows that his proposed oscillator has considerable potential in the fields of vacuum technology and electron physics.

Monte Carlo and molecular dynamics simulations of methane in potassium montmorillonite clay hydrates at elevated pressures and temperatures

The structure and dynamics of methane in hydrated potassium montmorillonite clay have been studied under conditions encountered in sedimentary basin and compared to those of hydrated sodium montmorillonite clay using computer simulation techniques. The simulated systems contain two molecular layers of water and followed gradients of 150 barkm-1 and 30 Kkm-1 up to a maximum burial depth of 6 km. Methane particle is coordinated to about 19 oxygen atoms, with 6 of these coming from the clay surface oxygen. Potassium ions tend to move away from the center towards the clay surface, in contrast to the behavior observed with the hydrated sodium form. The clay surface affinity for methane was found to be higher in the hydrated K-form. Methane diffusion in the two-layer hydrated K-montmorillonite increases from 0.39×10-9 m2s-1 at 280 K to 3.27×10-9 m2s-1 at 460 K compared to 0.36×10-9 m2s-1 at 280 K to 4.26×10-9 m2s-1 at 460 K in Na-montmorillonite hydrate. The distributions of the potassium ions were found to vary in the hydrates when compared to those of sodium form. Water molecules were also found to be very mobile in the potassium clay hydrates compared to sodium clay hydrates. © 2004 Elsevier Inc. All All rights reserved.

Tensile properties of as fabricated, aged, and stretched SiC whisker and particle reinforced 2009 aluminium alloy

Tensile tests were carried out using specimens of 2009 aluminium alloy reinforced by either SiC whiskers or particles. The size distributions of the whiskers and particles in the matrix were obtained by image analysis. It was found that failure was a result of uniform void nucleation and coalescence in the as fabricated composites, or a result of fast crack propagation initiated by a flaw developed at clusters of SiC in the aged or stretched and aged composites. The strengths of the as fabricated composites were estimated based on the results of image analysis using continuum mechanics and dislocation theories. The estimation indicated that the tensile strengths are largely contributed to by composite strengthening, supplemented by residual dislocation strengthening and work hardening. Owing to the flaw controlled failure, the tensile strengths of the aged or stretched and aged composites were independent of aging time, aging temperature, and the amount of stretching. The elastic moduli of the composites were estimated using the Halpin-Tsai model and a good correlation was found between the measured and estimated moduli. © 1996 The Institute of Materials.

Fatigue crack growth behaviour in molten-metal processed SiC particle-reinforced aluminium alloys

Fatigue crack initiation and subsequent short crack growth behaviour of 2014-5wt%SiC aluminium alloy composites has been examined in 4-point bend loading using smooth bar specimens. The growth rates of long fatigue cracks have also been measured at different stress ratios using pre-cracked specimens. The distributions of SiC particles and of coarse constituent particles in the matrix (which arise as a result of the molten-metal processing and relatively slow cooling rate) have been investigated. Preferential crack initiation sites were found to be SiC-matrix interfaces, SiC particles associated with constituent particles and the coarse constituent particles themselves. For microstructurally short cracks the dispersed SiC particles also act as temporary crack arresters. In the long crack growth tests, higher fatigue crack growth rates were obtained than for monolithic alloys. This effect is attributed to the contribution of void formation, due to the decohesion of SiC particles, to the fatigue crack growth process in the composite. Above crack depths of about 200 μm 'short' crack growth rates were in good agreement with the long crack data, showing a Pris exponent, m = 4 in both cases. For the long crack and short crack growth tests little effect of specimen orientation and grain size was observed on fatigue crack growth rates, but, specimen orientation affected the toughness. No effect of stress ratio in the range R = 0.2-0.5 was seen for long crack data in the Paris region.