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.

Flowback of solids through distribution plates of gas fluidized beds and associated phenomena

This work is concerned with a study of certain phenomena related to the performance and design of distributors in gas fluidized beds with particular regard to flowback of solid particles. The work to be described is divided into two parts. I. In Part one, a review of published material pertaining to distribution plates, including details from the patent specifications, has been prepared. After a chapter on the determination of the incipient fluidizing velocity, the following aspects of multi-orifice distributor plates in gas fluidized beds have been studied: (i) The effect of the distributor on bubble formation related to the way in which even distribution of bubbles on the top surface of the fluidized bed is obtained, e.g. the desirable pressure drop ratio ?PD/?PB for the even distribution of gas across the bed. Ratios of distributor pressure drop ?PD to bed pressure drop at which stable fluidization occurs show reasonable agreement with industrial practice. There is evidence that larger diameter beds tend to be less stable than smaller diameter beds when these are operated with shallow beds. Experiments show that in the presence of the bed the distributor pressure drop is reduced relative to the pressure drop without the bed, and this pressure drop in the former condition is regarded as the appropriate parameter for the design of the distributor. (ii) Experimental measurements of bubble distribution at the surface has been used to indicate maldistribution within the bed. Maldistribution is more likely at low gas flow rates and with distributors having large fractional free area characteristics (i.e. with distributors having low pressure drops). Bubble sizes obtained from this study, as well as those of others, have been successfully correlated. The correlation produced implies the existence of a bubble at the surface of an orifice and its growth by the addition of excess gas from the fluidized bed. (iii) For a given solid system, the amount of defluidized particles stagnating on the distributor plate is influenced by the orifice spacing, bed diameter and gas flow rate, but independent of the initial bed height and the way the orifices are arranged on the distributor plate. II. In Part two, solids flowback through single and multi-orifice distributors in two-dimensional and cylindrical beds of solids fluidized with air has been investigated. Distributors equipped with long cylindrical nozzles have also been included in the study. An equation for the prediction of free flowback of solids through multi-orifice distributors has been derived. Under fluidized conditions two regimes of flowback have been differentiated, namely Jumping and weeping. Data in the weeping regime have been successfully correlated. The limiting gas velocity through the distributor orifices at which flowback is completely excluded is found to be indepnndent of bed height, but a function of distributor design and physical properties of gas and solid used. A criterion for the prediction of this velocity has been established. The decisive advantage of increasing the distributor thickness or using nozzles to minimize solids flowback in fluidized beds has been observed and the opportunity taken to explore this poorly studied subject area. It has been noted, probably for the first time, that with long nozzles, there exists a critical nozzle length above which uncontrollable downflow of solids occurs. A theoretical model for predicting the critical length of a bundle of nozzles in terms of gas velocity through the nozzles has been set up. Theoretical calculations compared favourably with experiments.

CFD-modeling of insulation debris transport phenomena in water flow

The investigation of insulation debris generation, transport and sedimentation becomes important with regard to reactor safety research for PWR and BWR, when considering the long-term behavior of emergency core cooling systems during all types of loss of coolant accidents (LOCA). The insulation debris released near the break during a LOCA incident consists of a mixture of disparate particle population that varies with size, shape, consistency and other properties. Some fractions of the released insulation debris can be transported into the reactor sump, where it may perturb/impinge on the emergency core cooling systems. Open questions of generic interest are the sedimentation of the insulation debris in a water pool, its possible re-suspension and transport in the sump water flow and the particle load on strainers and corresponding pressure drop. A joint research project on such questions is being performed in cooperation between the University of Applied Sciences Zittau/Görlitz and the Forschungszentrum Dresden-Rossendorf. The project deals with the experimental investigation of particle transport phenomena in coolant flow and the development of CFD models for its description. While the experiments are performed at the University at Zittau/Görlitz, the theoretical modeling efforts are concentrated at Forschungszentrum Dresden-Rossendorf. In the current paper the basic concepts for CFD-modeling are described and feasibility studies including the conceptual design of the experiments are presented. © 2009 Elsevier B.V. All rights reserved.

Studies of the shellside performance of shell-and-tube heat exchangers

Accurate prediction of shellside pressure drop in a baffled shell-and-tube heat exchanger is very difficult because of the complicated shellside geometry. Ideally, all the shellside fluid should be alternately deflected across the tube bundle as it traverses from inlet to outlet. In practice, up to 60% of the shellside fluid may bypass the tube bundle or leak through the baffles. This short-circuiting of the main flow reduces the efficiency of the exchanger. Of the various shellside methods, it is shown that only the multi-stream methods, which attempt to obtain the shellside flow distribution, predict the pressure drop with any degree of accuracy, the various predictions ranging from -30% to +70%, generally overpredicting. It is shown that the inaccuracies are mainly due to the manner in which baffle leakage is modelled. The present multi-stream methods do not allow for interactions of the various flowstreams, and yet it is shown that three main effects are identified, a) there is a strong interaction between the main cross flow and the baffle leakage streams, enhancing the crossflow pressure drop, b) there is a further short-circuit not considered previously i.e. leakage in the window, and c) the crossflow does not penetrate as far, on average, as previously supposed. Models are developed for each of these three effects, along with a new windowflow pressure drop model, and it is shown that the effect of baffle leakage in the window is the most significant. These models developed to allow for various interactions, lead to an improved multi-stream method, named the "STREAM-INTERACTION" method. The overall method is shown to be consistently more accurate than previous methods, with virtually all the available shellside data being predicted to within ±30% and over 60% being within ±20%. The method is, thus, strongly recommended for use as a design method.

Coalescence of secondary dispersions in packed beds

A study has been made of the coalescence of secondary dispersions in beds of monosized glass ballotini. The variables investigated were superficial velocity, bed depth, ballotini size and dispersed phase concentration. Equipment was designed to generate a toluene ln water dispersion with phase ratios from 0.1 - 1.0 v/v % and whose mean drop size was determined using a Coulter Counter. The coalesced drops were sized by photography and the mean diameter of the effluent drops was determined using a Malvern Particle Size Analyser. Previous models describing single phase flow in porous media are reviewed and it was found that the experimental data obtained in this study is best represented by the Carman-Kozeny equations. Relative permeability correlations were used to predict the saturation profiles across the bed from measured two phase pressure drop data. Theoretical comparison of drop capture mechanisms indicated that direct and indirect interception are predominant. The total capture efficiency for the bed can also be evaluated using Spielman and Fitzpatrick's correlation.The resulting equation is used to predict the initial, local drop capture rate in a coalescer. A mathematical description of the saturation profiles is formulated and verified by the saturation profiles obtained by relative permeability. Based on the Carman-Kozeny equation, an expression is derived analytically to .predict the two phase pressure drop using the parameters which characterise the saturation profiles. By specifying the local saturation at the inlet face for a given velocity and phase ratio, good agreement between experimental pressure drop data and the model predictions was obtained. An attempt to predict the exit drop size has been made using an analogy for flow through non cylindrical channels.

Mechanisms of secondary dispersions separation in particulate beds

The mechanisms by which drops of secondary liquid dispersion ie. <100μ m, are collected, coalesced and transferred have been studied in particulate beds of different sizes and heights of glass ballotini. The apparatus facilitated different coalescer cell arrangements. The liquid-liquid system was toluene/de-ionised water. The inlet drop size distribution was measured by microscopy and using the Malvern Particle Size analyser; the outlet dispersion was sized by photography. The effect of packed height and packing size upon critical velocity, pressure drop and coalescence efficiency have been investigated. Single and two phase flow pressure drops across the packing were correlated by modified Blake-Kozeny equations. Two phase pressure drop was correlated by two equations, one for large ballotini sizes (267μm - 367μm), the other for small ballotini sizes (93μm- 147.5μm). The packings were efficient coalescers up to critical velocities of 3 x 10-2 m/s to 5 x 10-2 m/s. The saturation was measured across the bed using relative permeability and a mathematical model developed which related this profile to measured pressure drops. Filter coefficients for the range of packing studied were found to be accurately predicted from a modified queueing drop model. 

Transfer coefficients over inclined tube banks

The thesis presents experimental results for shell-side transfer coefficients and pressure drops across four different tube banks, using small-scale models, with yawed tubes, as found in many types of heat exchangers, boilers and nuclear reactors. The tube banks investigated have a staggered tube layout on a rotated square pitch, with a 1.25 pitch-to-diameter ratio. The angle of attack was varied between 45o and 90o. An extensive range of Reynolds number, i.e. 0.5. to 12,600, covering so-called laminar, transition and turbulent flows, was investigated. A diffusion-controlled electrochemical mass transfer technique has been employed to measure mass transfer coefficients. The heat transfer coefficients may be then readily obtained from the mass transfer values by applying the well-established Chilton-Colburn analogy. The results for the normal tube bank, which forms the base case for the study on inclined tube banks, show close agreement with previous work. The transfer coefficients and pressure drops of the inclined tube banks are compared with results from the ideal normal tube bank to examine the effect of inclination angle on heat transfer and pressure drop variations. The variation of the transfer coefficients row-by-row and the entrance and exit effects have also been investigated. An auxilary investigation has been carried out on the role of natural convection. A preliminary correlation of transfer coefficients and pressure drops against the variation in the yaw angle has been attempted. The results are discussed in the light of the few existing theoretical treatments and experimental data for these situations, and recommendations made for future work.

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.

Fin spacing interaction in fluidized bed heat exchangers

The work presented in this thesis is concerned with the heat transfer performance of a single horizontal bare tube and a variety of finned tubes immersed in a shallow air fluidized bed. Results of experimental investigations with the bare tube indicate that the tube position in the bed influences its performance narticularly where fine bed materials are used. In some cases the maximum heat transfer is obtained with the tube in the particle cloud just above the dense phase fluidized bed - a phenomenon that has not been previously observed. This was attributed to the unusual particle circulation in shallow beds. The data is also presented in dimensionless correlations which may be useful for design purposes. A close approximation to the bare tube data can be obtained by using thetransient heating of a spherical robe and this provides a valuable way of accumulating a lot of data very rapidly. The experimental data on finned tubes shows that a fin spacing less than twenty times the average particle diameter can cause a significant reduction in heat transfer due to the interaction which takes place between the particles and the surface of the fins. Furthermore, evidence is provided to show that particle shape plays an important part in the interaction with spherical particles being superior to angular particles at low fin spacing/particle diameter ratio. The finned tube data is less sensitive to tube position in the bed than bare tubes and the best performance is when the tube is positioned at the distributor.A reduction in bed depth decreases the thermal performance of the finned tube but in many practical installations the reduction in pressure drop might more than comnensate for the reduced heat flux. Information is also provided on the theoretical uerformance of fins and the effect of the root contact area between the fins and the tube was investigated.

Constant shear along a boundary of a flume

'I'he accurate rreasurement of bed shear stress has been extremely difficult due to its changing values until white propunded a theory which would give constant shear along the bed of a flume. In this investigation a flume has been designed according to White's theory and by two separate methods proven to give constant shearing force along the bed. The first method applied the Hydrogen Bubble Technique to obtain accurate values of velocity thus allowing the velocity profile to be plotted and the momentum at the various test sections to be calculated. The use of a 16 mm Beaulieu movie camera allowed the exact velocity profiles created by the hydrogen bubbles to be recorded whilst an analysing projector gave the means of calculating the exact velocities at the various test sections. Simultaneously Preston's technique of measuring skin friction using Pitot tubes was applied. Twc banks of open ended water manometer were used for recording the static and velocity head pressure drop along the flume. This tvpe of manometer eliminated air locks in the tubes and was found to be sufficiently accurate. Readings of pressure and velocity were taken for various types and diameters of bed material both natural sands and glass spheres and the results tabulated. Graphs of particle Reynolds Number against bed shear stress were plotted and gave a linear relationship which dropped off at high values of Reynolds number. It was found that bed movement occurred instantaneously along the bed of the flume once critical velocity had been reached. On completion of this test a roof curve inappropriate to the bed material was used and then the test repeated. The bed shearing stress was now no longer constant and yet bed movement started instantaneously along the bed of the flume, showing that there are more parameters than critical shear stress to bed movement. It is concluded from the two separate methods applied that the bed shear stress is constant along the bed of the flume.

Upgrading of fast pyrolysis oils by hot filtration

A hot filtration unit downstream of a 1kg/h fluidised bed fast pyrolysis reactor was designed and built. The filter unit operates at 450oC and consists of 1 exchangeable filter candle with reverse pulse cleaning system. Hot filtration experiments up to 7 hours were performed with beech wood as feedstock. It was possible to produce fast pyrolysis oils with a solid content below 0.01 wt%. The additional residence time of the pyrolysis vapours and secondary vapour cracking on the filter cake caused an increase of non-condensable gases at the expense of organic liquid yield. The oils produced with hot filtration showed superior quality properties regarding viscosity than standard pyrolysis oils. The oils were analysed by rotational viscosimetry and gel permeation chromatography before and after accelerated aging. During filtration the separated particulates accumulate on the candle surface and build up the filter cake. The filter cake leads to an increase in pressure drop between the raw gas and the clean gas side of the filter candle. At a certain pressure drop the filter cake has to be removed by reverse pulse cleaning to regenerate the pressure drop. The experiments showed that successful pressure drop recovery was possible during the initial filtration cycles, thereafter further cycles showed minor pressure drop recovery and therefore a steady increase in differential pressure. Filtration with pre-coating the candle to form an additional layer between the filter candle and cake resulted in total removal of the dust cake.

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

CFD-modeling and Experiments of insulation debris transport phenomena in water flow

The investigation of insulation debris generation, transport and sedimentation becomes important with regard to reactor safety research for PWR and BWR, when considering the long-term behavior of emergency core cooling systems during all types of loss of coolant accidents (LOCA). The insulation debris released near the break during a LOCA incident consists of a mixture of disparate particle population that varies with size, shape, consistency and other properties. Some fractions of the released insulation debris can be transported into the reactor sump, where it may perturb/impinge on the emergency core cooling systems. Open questions of generic interest are the sedimentation of the insulation debris in a water pool, its possible re-suspension and transport in the sump water flow and the particle load on strainers and corresponding pressure drop. A joint research project on such questions is being performed in cooperation between the University of Applied Sciences Zittau/Görlitz and the Forschungszentrum Dresden-Rossendorf. The project deals with the experimental investigation of particle transport phenomena in coolant flow and the development of CFD models for its description. While the experiments are performed at the University at Zittau/Görlitz, the theoretical modeling efforts are concentrated at Forschungszentrum Dresden-Rossendorf. In the current paper the basic concepts for CFD modeling are described and feasibility studies including the conceptual design of the experiments are presented.

Numerical and experimental investigations for insulation particle transport phenomena in water flow

The investigation of insulation debris generation, transport and sedimentation becomes more important with regard to reactor safety research for pressurized and boiling water reactors, when considering the long-term behaviour of emergency core coolant systems during all types of loss of coolant accidents (LOCA). The insulation debris released near the break during a LOCA incident consists of a mixture of a disparate particle population that varies with size, shape, consistency and other properties. Some fractions of the released insulation debris can be transported into the reactor sump, where it may perturb or impinge on the emergency core cooling systems. Open questions of generic interest are for example the particle load on strainers and corresponding pressure-drop, the sedimentation of the insulation debris in a water pool, its possible re-suspension and transport in the sump water flow. A joint research project on such questions is being performed in cooperation with the University of Applied Science Zittau/Görlitz and the Forschungszentrum Dresden-Rossendorf. The project deals with the experimental investigation and the development of computational fluid dynamic (CFD) models for the description of particle transport phenomena in coolant flow. While the experiments are performed at the University Zittau/Görlitz, the theoretical work is concentrated at Forschungszentrum Dresden-Rossendorf. In the present paper, the basic concepts for computational fluid dynamic (CFD) modelling are described and experimental results are presented. Further experiments are designed and feasibility studies were performed.