Quantitative risk assessment of groundwater quality utilizing GIS technology and coupled groundwater models

The thesis presents a two-dimensional Risk Assessment Method (RAM) where the assessment of risk to the groundwater resources incorporates both the quantification of the probability of the occurrence of contaminant source terms, as well as the assessment of the resultant impacts. The approach emphasizes the need for a greater dependency on the potential pollution sources, rather than the traditional approach where assessment is based mainly on the intrinsic geo-hydrologic parameters. The risk is calculated using Monte Carlo simulation methods whereby random pollution events were generated to the same distribution as historically occurring events or a priori potential probability distribution. Integrated mathematical models then simulate contaminant concentrations at the predefined monitoring points within the aquifer. The spatial and temporal distributions of the concentrations were calculated from repeated realisations, and the number of times when a user defined concentration magnitude was exceeded is quantified as a risk. The method was setup by integrating MODFLOW-2000, MT3DMS and a FORTRAN coded risk model, and automated, using a DOS batch processing file. GIS software was employed in producing the input files and for the presentation of the results. The functionalities of the method, as well as its sensitivities to the model grid sizes, contaminant loading rates, length of stress periods, and the historical frequencies of occurrence of pollution events were evaluated using hypothetical scenarios and a case study. Chloride-related pollution sources were compiled and used as indicative potential contaminant sources for the case study. At any active model cell, if a random generated number is less than the probability of pollution occurrence, then the risk model will generate synthetic contaminant source term as an input into the transport model. The results of the applications of the method are presented in the form of tables, graphs and spatial maps. Varying the model grid sizes indicates no significant effects on the simulated groundwater head. The simulated frequency of daily occurrence of pollution incidents is also independent of the model dimensions. However, the simulated total contaminant mass generated within the aquifer, and the associated volumetric numerical error appear to increase with the increasing grid sizes. Also, the migration of contaminant plume advances faster with the coarse grid sizes as compared to the finer grid sizes. The number of daily contaminant source terms generated and consequently the total mass of contaminant within the aquifer increases in a non linear proportion to the increasing frequency of occurrence of pollution events. The risk of pollution from a number of sources all occurring by chance together was evaluated, and quantitatively presented as risk maps. This capability to combine the risk to a groundwater feature from numerous potential sources of pollution proved to be a great asset to the method, and a large benefit over the contemporary risk and vulnerability methods.

A CFD approach on the effect of particle size on char entrainment in bubbling fluidised bed reactors

The fluid – particle interaction inside a 41.7 mg s-1 fluidised bed reactor is modelled. Three char particles of sizes 500 µm, 250 µm, and 100 µm are injected into the fluidised bed and the momentum transport from the fluidising gas and fluidised sand is modelled. Due to the fluidising conditions and reactor design the char particles will either be entrained from the reactor or remain inside the bubbling bed. The particle size is the factor that differentiates the particle motion inside the reactor and their efficient entrainment out of it. A 3-Dimensional simulation has been performed with a completele revised momentum transport model for bubble three-phase flow according to the literature as an extension to the commercial finite volume code FLUENT 6.2.

3D simulation of the effects of sphericity on char entrainment in fluidised beds

The paper presents a 3-dimensional simulation of the effect of particle shape on char entrainment in a bubbling fluidised bed reactor. Three char particles of 350 μm side length but of different shapes (cube, sphere, and tetrahedron) are injected into the fluidised bed and the momentum transport from the fluidising gas and fluidised sand is modelled. Due to the fluidising conditions, reactor design and particle shape the char particles will either be entrained from the reactor or remain inside the bubbling bed. The sphericity of the particles is the factor that differentiates the particle motion inside the reactor and their efficient entrainment out of it. The simulation has been performed with a completely revised momentum transport model for bubble three-phase flow, taking into account the sphericity factors, and has been applied as an extension to the commercial finite volume code FLUENT 6.3. © 2010 Elsevier B.V.All rights reserved.

Molecular phase space transport in water:non-stationary random walk model

Molecular transport in phase space is crucial for chemical reactions because it defines how pre-reactive molecular configurations are found during the time evolution of the system. Using Molecular Dynamics (MD) simulated atomistic trajectories we test the assumption of the normal diffusion in the phase space for bulk water at ambient conditions by checking the equivalence of the transport to the random walk model. Contrary to common expectations we have found that some statistical features of the transport in the phase space differ from those of the normal diffusion models. This implies a non-random character of the path search process by the reacting complexes in water solutions. Our further numerical experiments show that a significant long period of non-stationarity in the transition probabilities of the segments of molecular trajectories can account for the observed non-uniform filling of the phase space. Surprisingly, the characteristic periods in the model non-stationarity constitute hundreds of nanoseconds, that is much longer time scales compared to typical lifetime of known liquid water molecular structures (several picoseconds).

Experiments to assess the transport of fibrous insulation debris

Loss of coolant accidents (LOCA) in the primary cooling circuit of a nuclear reactor may result in damage to insulation materials that are located near to the leak. The insulation materials released may compromise the operation of the emergency core cooling system (ECCS). Insulation material in the form of mineral wool fibre agglomerates (MWFA) maybe transported to the containment sump strainers mounted at the inlet of the emergency cooling pumps, where the insulation fibres 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 fibre cakes on the strainers, which can lead to a significant increase in the strainer pressure drop and result in cavitation in the ECCS. Thus, knowledge of transport characteristics of the damaged insulation materials in various scenarios is required to help plan for 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 LOCA. The study entails the generation of fibre agglomerates, the determination of their transport properties in single and multi-effect experiments and the long-term effect that corrosion of the containment internals by the coolant has on the strainer pressure drop. The focus of this presentation is on the experiments performed that characterize the horizontal transport of MWFA, whereas the corresponding CFD simulations are described in an accompanying contribution (see abstract of Cartland Glover et al.). The experiments were performed a racetrack type channel that provided a near uniform horizontal flow. The channel is 0.1 wide by 1.2 m high with a straight length of 5 m and two bends of 0.5 m. The measurement techniques include particle imaging (both wide-angle and macro lens), concurrent particle image velocimetry, ultravelocimetry, laser detection sensors to sense the presence of absence of MWFA and pertinent measurements of the MWFA concentration and quiescent settling characteristics. The transport of the MWFA was observed at velocities of 0.1 and 0.25 m s-1 to verify numerical model behaviour in and just beyond expected velocities in the containment sump of a nuclear reactor.

Fibre agglomerate transport in a horizontal flow

An experimental and theoretical study of the transport of mineral wool fibre agglomerates in nuclear power plant containment sumps is being performed. A racetrack channel was devised to provide data for the validation of numerical models, which are intended to model the transport of fibre agglomerates. The racetrack channel provides near uniform and steady conditions that lead to either the sedimentation or suspension of the agglomerates. Various experimental techniques were used to determine the velocity conditions and the distribution of the fibre agglomerates in the channel. The fibre agglomerates are modelled as fluid particles in the Eulerian reference frame. Simulations of pure sedimentation of a known mass and volume of agglomerations show that the transport of the fibre agglomerates can be replicated. The suspension of the fibres is also replicated in the simulations; however, the definition of the fibre agglomerate phase is strongly dependent on the selected density and diameter. Detailed information on the morphology of the fibre agglomerates is lacking for the suspension conditions, as the fibre agglomerates may undergo breakage and erosion. Therefore, ongoing work, which is described here, is being pursued to improve the experimental characterisation of the suspended transport of the fibre agglomerates.

Generic experiments at the sump model "Zittauer Strömungswanne" (ZSW) for the behaviour of mineral wool in the sump and the reactor core

The investigation of insulation debris transport, sedimentation, penetration into the reactor core and head loss build up becomes important to reactor safety research for PWR and BWR, when considering the long-term behaviour of emergency core cooling systems during loss of coolant accidents. Research projects are being performed in cooperation between the University of Applied Sciences Zittau/Görlitz and the Helmholtz-Zentrum Dresden-Rossendorf. The projects include experimental investigations of different processes and phenomena of insulation debris in coolant flow and the development of CFD models. Generic complex experiments serve for building up a data base for the validation of models for single effects and their coupling in CFD codes. This paper includes the description of the experimental facility for complex generic experiments (ZSW), an overview about experimental boundary conditions and results for upstream and down-stream phenomena as well as for the long-time behaviour due to corrosive processes. © Carl Hanser Verlag, München.