Simulations of agglomerate sedimentation and suspension

Mineral wool insulation material applied to the primary cooling circuit of a nuclear reactor maybe damaged in the course of a loss of coolant accident (LOCA). The insulation material released by the leak may compromise the operation of the emergency core cooling system (ECCS), as it maybe transported together with the coolant in the form of mineral wool fiber agglomerates (MWFA) suspensions to the containment sump strainers, which are mounted at the inlet of the ECCS to keep any debris away from the emergency cooling pumps. In the further course of the LOCA, the MWFA may block or penetrate the strainers. In addition to the impact of MWFA on the pressure drop across the strainers, corrosion products formed over time may also accumulate in the fiber cakes on the strainers, which can lead to a significant increase in the strainer pressure drop and result in cavitation in the ECCS. Therefore, it is essential to understand the transport characteristics of the insulation materials in order to determine the long-term operability of nuclear reactors, which undergo LOCA. An experimental and theoretical study performed by the Helmholtz-Zentrum Dresden-Rossendorf and the Hochschule Zittau/Görlitz is investigating the phenomena that maybe observed in the containment vessel during a primary circuit coolant leak. The study entails the generation of fiber agglomerates, the determination of their transport properties in single and multi-effect experiments and the long-term effects that particles formed due to corrosion of metallic containment internals by the coolant medium have on the strainer pressure drop. The focus of this presentation is on the numerical models that are used to predict the transport of MWFA by CFD simulations. A number of pseudo-continuous dispersed phases of spherical wetted agglomerates can represent the MWFA. The size, density, the relative viscosity of the fluid-fiber agglomerate mixture and the turbulent dispersion all affect how the fiber agglomerates are transported. In the cases described here, the size is kept constant while the density is modified. This definition affects both the terminal velocity and volume fraction of the dispersed phases. Application of such a model to sedimentation in a quiescent column and a horizontal flow are examined. The scenario also presents the suspension and horizontal transport of a single fiber agglomerate phase in a racetrack type channel.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors. Part B: Heat, momentum and mass transport in bubbling fluidised beds

The fluid–particle interaction inside a 150 g/h fluidised bed reactor is modelled. The biomass particle is injected into the fluidised bed and the heat, momentum and mass transport from the fluidising gas and fluidised sand is modelled. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Heat transfer from the bubbling bed to the discrete biomass particle, as well as biomass reaction kinetics are modelled according to the literature. 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). The study completes the fast pyrolysis modelling in bubbling fluidised bed reactors.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors, Part A: Eulerian computation of momentum transport in bubbling fluidised beds

The fluid–particle interaction inside a 150 g/h fluidised bed reactor is modelled. The biomass particle is injected into the fluidised bed and the momentum transport from the fluidising gas and fluidised sand is modelled. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase, according to the literature. FLUENT 6.2 has been used as the modelling framework of the simulations with a completely revised drag model, in the form of user defined function (UDF), to calculate the forces exerted on the particle as well as its velocity components. 2-D and 3-D simulations are tested and compared. The study is the first part of a complete pyrolysis model in fluidised bed reactors.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors: modelling the impact of biomass shrinkage

The fluid–particle interaction and the impact of shrinkage on pyrolysis of biomass inside a 150 g/h fluidised bed reactor is modelled. Two 500 View the MathML sourcem in diameter biomass particles are injected into the fluidised bed with different shrinkage conditions. The two different conditions consist of (1) shrinkage equal to the volume left by the solid devolatilization, and (2) shrinkage parameters equal to approximately half of particle volume. The effect of shrinkage is analysed in terms of heat and momentum transfer as well as product yields, pyrolysis time and particle size considering spherical geometries. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Heat transfer from the bubbling bed to the discrete biomass particle, as well as biomass reaction kinetics are modelled according to the literature. 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).

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).

Use of computational fluid dynamics in the design of bubble column reactors

Investigations into the modelling techniques that depict the transport of discrete phases (gas bubbles or solid particles) and model biochemical reactions in a bubble column reactor are discussed here. The mixture model was used to calculate gas-liquid, solid-liquid and gasliquid-solid interactions. Multiphase flow is a difficult phenomenon to capture, particularly in bubble columns where the major driving force is caused by the injection of gas bubbles. The gas bubbles cause a large density difference to occur that results in transient multi-dimensional fluid motion. Standard design procedures do not account for the transient motion, due to the simplifying assumptions of steady plug flow. Computational fluid dynamics (CFD) can assist in expanding the understanding of complex flows in bubble columns by characterising the flow phenomena for many geometrical configurations. Therefore, CFD has a role in the education of chemical and biochemical engineers, providing the examples of flow phenomena that many engineers may not experience, even through experimentation. The performance of the mixture model was investigated for three domains (plane, rectangular and cylindrical) and three flow models (laminar, k-e turbulence and the Reynolds stresses). mThis investigation raised many questions about how gas-liquid interactions are captured numerically. To answer some of these questions the analogy between thermal convection in a cavity and gas-liquid flow in bubble columns was invoked. This involved modelling the buoyant motion of air in a narrow cavity for a number of turbulence schemes. The difference in density was caused by a temperature gradient that acted across the width of the cavity. Multiple vortices were obtained when the Reynolds stresses were utilised with the addition of a basic flow profile after each time step. To implement the three-phase models an alternative mixture model was developed and compared against a commercially available mixture model for three turbulence schemes. The scheme where just the Reynolds stresses model was employed, predicted the transient motion of the fluids quite well for both mixture models. Solid-liquid and then alternative formulations of gas-liquid-solid model were compared against one another. The alternative form of the mixture model was found to perform particularly well for both gas and solid phase transport when calculating two and three-phase flow. The improvement in the solutions obtained was a result of the inclusion of the Reynolds stresses model and differences in the mixture models employed. The differences between the alternative mixture models were found in the volume fraction equation (flux and deviatoric stress tensor terms) and the viscosity formulation for the mixture phase.

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 influence of non-metallic inclusions upon the properties of linepipe steels

The principal aim of this work was to determine the role of non-metallic inclusions in the process of hydrogen stepwise cracking (SWC). Additionally, the influence of inclusions upon the notch ductility of hydrogen charged (HC) and uncharged (UN) tensile specimens was examined. To obtain a basis for experiment a series of low carbon-manganese steels were prepared by induction melting. In order to produce variations in the composition, morphology, volume fraction, size and distribution of the inclusions the steel chemistry was adjusted prior to casting by additions of deoxidiser and Ca-Si injection. Sections of each ingot were hot rolled. Metallography, image analysis, mechanical tests and hydrogen SWC tests were then carried out. The volume fraction, morphology, and shape of inclusions influenced the tensile ductility of the steels. Marked anisotropy was found in the steels containing type II MnS inclusions at all rolling temperatures, whereas the fully Ca treated steel was isotropic. It was found that several inclusion parameters (projected length PL, mean free distance MFD, nearest-neighbour distance NND) correlated with fracture strain. An increase in inclusion volume fraction and/or the dimension of inclusions on a plane parallel to the plane of fracture led to a decrease in fracture strain. The inclusion parameters did not correlate with the fracture strains for the HC tensile specimens. However, large or clusters of inclusions acted as the principal sites for crack initiation. `Fisheyes' or areas of `flat' fracture were often found on these fracture surfaces. The criteria for SWC initiation was found to be either large inclusions or clusters of inclusions. As the PL of inclusions increased the probability of large SWCs occurring increased. SWC initiation at inclusions was believed to occur at a critical concentration of hydrogen. Factors which assisted the concentration of hydrogen at inclusions were discussed. None of the proposed mechanisms of hydrogen embrittlement could be identified as the single cause of SWC.

The mechanical properties and formability of dual-phase steel

In recent years dual phase steels comprising of 5-20% martensite in a ferrite matrix have come into the limelight of high strength cold formable steels because of their potential for vehicle weight saving. They show the following features: no yield point; relatively low initial flow stress; high initial workhardening rate; well sustained work hardening. As a consequence of these characteristics, dual phase steels exhibit a better combination of strength and elongation than other HSLA steels. In this thesis, a broad view of the factors which influence their properties is presented. Mechanical properties and forming ability of a commercially available dual phase steel and an AL-Si killed steel processed to dual phase form are investigated to ascertain the effect of their microstructure on their properties. It is found that the yield phenomena are masked by the transformation induced stresses present during processing and so yield point could be recovered under suitable ageing treatment; that apart from giving the above properties dual phasing gives rise to very low strain-rate sensitivity and a low R value ~ 1; that the mechanical response under rolling conditions is not different from those under tension; that there is a danger of damage to tooling during forming operations of these steels if fracture should precede instability as a result of grain size dependent strength found for these steels. It is also found that very little deformation of the martensite islands took place during deformation except at high strains. The work-hardening and the strength levels can be controlled by either decreasing the grain size or increasing the martensite volume fraction, but it is found that increasing martensite has a detrimental effect on ductility and the ductility and fracture strength can be controlled better by refining the grain size. A remarkable effect found in the dual phase steel tested is that the compressive strength is higher than the tensile strength. The reason for this observation is not yet clear but it is suggested that it might be due to the introduction of emissary type dislocations into the ferrite lattice as a result of twins formed in the martensite during transformation from austenite. The twins are envisaged to be {111} <112> in character.

Verification and validation of numerical models of the transport of insulation debris

Damage to insulation materials located near to a primary circuit coolant leak may compromise the operation of the emergency core cooling system (ECCS). Insulation material in the form of mineral wool fiber agglomerates (MWFA) maybe transported to the containment sump strainers, where they may block or penetrate the strainers. Though the impact of MWFA on the pressure drop across the strainers is minimal, corrosion products formed over time may also accumulate in the fiber cakes on the strainers, which can lead to a significant increase in the strainer pressure drop and result in cavitation in the ECCS. An experimental and theoretical study performed by the Helmholtz-Zentrum Dresden-Rossendorf and the Hochschule Zittau/Görlitz is investigating the phenomena that maybe observed in the containment vessel during a primary circuit coolant leak. The study entails the generation of fiber agglomerates, the determination of their transport properties in single and multi-effect experiments and the long-term effect that corrosion and erosion of the containment internals by the coolant has on the strainer pressure drop. The focus of this paper is on the verification and validation of numerical models that can predict the transport of MWFA. A number of pseudo-continuous dispersed phases of spherical wetted agglomerates represent the MWFA. The size, density, the relative viscosity of the fluid-fiber agglomerate mixture and the turbulent dispersion all affect how the fiber agglomerates are transported. In the cases described here, the size is kept constant while the density is modified. This definition affects both the terminal velocity and volume fraction of the dispersed phases. Note that the relative viscosity is only significant at high concentrations. Three single effect experiments were used to provide validation data on the transport of the fiber agglomerates under conditions of sedimentation in quiescent fluid, sedimentation in a horizontal flow and suspension in a horizontal flow. The experiments were performed in a rectangular column for the quiescent fluid and a racetrack type channel that provided a near uniform horizontal flow. The numerical models of sedimentation in the column and the racetrack channel found that the sedimentation characteristics are consistent with the experiments. For channel suspension, the heavier fibers tend to accumulate at the channel base even at high velocities, while lighter phases are more likely to be transported around the channel.

Numerical models used for the modelling of the transport of fibrous insulation debris

Mineral wool insulation material applied to the primary cooling circuit of a nuclear reactor maybe damaged in the course of a loss of coolant accident (LOCA). The insulation material released by the leak may compromise the operation of the emergency core cooling system (ECCS), as it maybe transported together with the coolant in the form of mineral wool fiber agglomerates (MWFA) suspensions to the containment sump strainers, which are mounted at the inlet of the ECCS to keep any debris away from the emergency cooling pumps. In the further course of the LOCA, the MWFA may block or penetrate the strainers. In addition to the impact of MWFA on the pressure drop across the strainers, corrosion products formed over time may also accumulate in the fiber cakes on the strainers, which can lead to a significant increase in the strainer pressure drop and result in cavitation in the ECCS. Therefore, it is essential to understand the transport characteristics of the insulation materials in order to determine the long-term operability of nuclear reactors, which undergo LOCA. An experimental and theoretical study performed by the Helmholtz-Zentrum Dresden-Rossendorf and the Hochschule Zittau/Görlitz1 is investigating the phenomena that maybe observed in the containment vessel during a primary circuit coolant leak. The study entails the generation of fiber agglomerates, the determination of their transport properties in single and multi-effect experiments and the long-term effects that particles formed due to corrosion of metallic containment internals by the coolant medium have on the strainer pressure drop. The focus of this presentation is on the numerical models that are used to predict the transport of MWFA by CFD simulations. A number of pseudo-continuous dispersed phases of spherical wetted agglomerates can represent the MWFA. The size, density, the relative viscosity of the fluid-fiber agglomerate mixture and the turbulent dispersion all affect how the fiber agglomerates are transported. In the cases described here, the size is kept constant while the density is modified. This definition affects both the terminal velocity and volume fraction of the dispersed phases. Only one of the single effect experimental scenarios is described here that are used in validation of the numerical models. The scenario examines the suspension and horizontal transport of the fiber agglomerates in a racetrack type channel. The corresponding experiments will be described in an accompanying presentation (see abstract of Seeliger et al.).

The role of organoclay in promoting co-continuous morphology in high-density poly(ethylene)/poly(amide) 6 blends

The effect of organically modified clay on the morphology, rheology and mechanical properties of high-density polyethylene (HDPE) and polyamide 6 (PA6) blends (HDPE/PA6 = 75/25 parts) is studied. Virgin and filled blends were prepared by melt compounding the constituents using a twin-screw extruder. The influence of the organoclay on the morphology of the hybrid was deeply investigated by means of wide-angle X-ray diffractometry, transmission and scanning electron microscopies and quantitative extraction experiments. It has been found that the organoclay exclusively places inside the more hydrophilic polyamide phase during the melt compounding. The extrusion process promotes the formation of highly elongated and separated organoclay-rich PA6 domains. Despite its low volume fraction, the filled minor phase eventually merges once the extruded pellets are melted again, giving rise to a co-continuous microstructure. Remarkably, such a morphology persists for long time in the melt state. A possible compatibilizing action related to the organoclay has been investigated by comparing the morphology of the hybrid blend with that of a blend compatibilized using an ethylene–acrylic acid (EAA) copolymer as a compatibilizer precursor. The former remains phase separated, indicating that the filler does not promote the enhancement of the interfacial adhesion. The macroscopic properties of the hybrid blend were interpreted in the light of its morphology. The melt state dynamics of the materials were probed by means of linear viscoelastic measurements. Many peculiar rheological features of polymer-layered silicate nanocomposites based on single polymer matrix were detected for the hybrid blend. The results have been interpreted proposing the existence of two distinct populations of dynamical species: HDPE not interacting with the filler, and a slower species, constituted by the organoclay-rich polyamide phase, which slackened dynamics stabilize the morphology in the melt state. In the solid state, both the reinforcement effect of the filler and the co-continuous microstructure promote the enhancement of the tensile modulus. Our results demonstrate that adding nanoparticles to polymer blends allows tailoring the final properties of the hybrid, potentially leading to high-performance materials which combine the advantages of polymer blends and the merits of polymer nanocomposites.

Gas-liquid-solid flow modelling in a bubble column

An alternative approach to the modelling of solid-liquid and gas-liquid-solid flows for a 5:1 height to width aspect ratio bubble column is presented here. A modified transport equation for the volume fraction of a dispersed phase has been developed for the investigation of turbulent buoyancy driven flows (Chem. Eng. Proc., in press). In this study, a modified transport equation has been employed for discrete phase motion considering both solid-liquid and gas-liquid-solid flows. The modelling of the three-phase flow in a bubble column was achieved in the following case: injecting a slug of solid particles into the column for 10 s at a velocity of 0.1 m s-1 and then the gas phase flow was initiated with a superficial gas velocity of 0.02 cm s-1. © 2003 Elsevier B.V. All rights reserved.

Gold nanoparticles decorated with oligo(ethylene glycol) thiols:protein resistance and colloidal stability

The interactions between proteins and gold colloids functionalized with protein-resistant oligo(ethylene glycol) (OEG) thiol, HS(CH(2))(11) (OCH(2)CH(2))(6)OMe (EG(6)OMe), in aqueous solution have been studied by small-angle X-ray scattering (SAXS) and UV-vis spectroscopy. The mean size, 2R, and the size distribution of the decorated gold colloids have been characterized by SAXS. The monolayer-protected gold colloids have no correlations due to the low volume fraction in solution and are stable in a wide range of temperatures (5-70 degrees C, pH (1.3-12.4), and ionic strength (0-1.0 M). In contrast, protein (bovine serum albumin) solutions with concentrations in the range of 60-200 mg/mL (4.6-14.5 vol show a pronounced correlation peak in SAXS, which results from the repulsive electrostatic interaction between charged proteins. These protein interactions show significant dependence on ionic strength, as would be expected for an electrostatic interaction (Zhang et al. J. Phys. Chem. B 2007, 111, 251). For a mixture of proteins and gold colloids, the protein-protein interaction changes little upon mixing with OEG-decorated gold colloids. In contrast, the colloid-colloid interaction is found to be strongly dependent on the protein concentration and the size of the colloid itself. Adding protein to a colloidal solution results in an attractive depletion interaction between functionalized gold colloids, and above a critical protein concentration, c*, the colloids form aggregates and flocculate. Adding salt to such mixtures enhances the depletion effect and decreases the critical protein concentration. The aggregation is a reversible process (i.e., diluting the solution leads to dissolution of aggregates). The results also indicate that the charge of the OEG self-assembled monolayer at a curved interface has a rather limited effect on the colloidal stabilization and the repulsive interaction with proteins.

Protein interactions studied by SAXS:effect of ionic strength and protein concentration for BSA in aqueous solutions

We have studied a series of samples of bovine serum albumin (BSA) solutions with protein concentration, c, ranging from 2 to 500 mg/mL and ionic strength, I, from 0 to 2 M by small-angle X-ray scattering (SAXS). The scattering intensity distribution was compared to simulations using an oblate ellipsoid form factor with radii of 17 x 42 x 42 A, combined with either a screened Coulomb, repulsive structure factor, S-SC(q), or an attractive square-well structure factor, S-SW(q). At pH = 7, BSA is negatively charged. At low ionic strength, I <0.3 M, the total interaction exhibits a decrease of the repulsive interaction when compared to the salt-free solution, as the net surface charge is screened, and the data can be fitted by assuming an ellipsoid form factor and screened Coulomb interaction. At moderate ionic strength (0.3-0.5 M), the interaction is rather weak, and a hard-sphere structure factor has been used to simulate the data with a higher volume fraction. Upon further increase of the ionic strength (I >= 1.0 M), the overall interaction potential was dominated by an additional attractive potential, and the data could be successfully fitted by an ellipsoid form factor and a square-well potential model. The fit parameters, well depth and well width, indicate that the attractive potential caused by a high salt concentration is weak and long-ranged. Although the long-range, attractive potential dominated the protein interaction, no gelation or precipitation was observed in any of the samples. This is explained by the increase of a short-range, repulsive interaction between protein molecules by forming a hydration layer with increasing salt concentration. The competition between long-range, attractive and short-range, repulsive interactions accounted for the stability of concentrated BSA solution at high ionic strength.