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

A consequence of a loss of coolant accident is that the local insulation material is damaged and maybe transported to the containment sump where it can penetrate and/or block the sump strainers. An experimental and theoretical study, which examines the transport of mineral wool fibers via single and multi-effect experiments is being performed. This paper focuses on the experiments and simulations performed for validation of numerical models of sedimentation and resuspension of mineral wool fiber agglomerates in a racetrack type channel. Three velocity conditions are used to test the response of two dispersed phase fiber agglomerates to two drag correlations and to two turbulent dispersion coefficients. The Eulerian multiphase flow model is applied with either one or two dispersed phases.

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 application of contact mechanics to the numerical simulation of particulate material

Particulate solids are complex redundant systems which consist of discrete particles. The interactions between the particles are complex and have been the subject of many theoretical and experimental investigations. Invetigations of particulate material have been restricted by the lack of quantitative information on the mechanisms occurring within an assembly. Laboratory experimentation is limited as information on the internal behaviour can only be inferred from measurements on the assembly boundary, or the use of intrusive measuring devices. In addition comparisons between test data are uncertain due to the difficulty in reproducing exact replicas of physical systems. Nevertheless, theoretical and technological advances require more detailed material information. However, numerical simulation affords access to information on every particle and hence the micro-mechanical behaviour within an assembly, and can replicate desired systems. To use a computer program to numerically simulate material behaviour accurately it is necessary to incorporte realistic interaction laws. This research programme used the finite difference simulation program `BALL', developed by Cundall (1971), which employed linear spring force-displacement laws. It was thus necessary to incorporate more realistic interaction laws. Therefore, this research programme was primarily concerned with the implementation of the normal force-displacement law of Hertz (1882) and the tangential force-displacement laws of Mindlin and Deresiewicz (1953). Within this thesis the contact mechanics theories employed in the program are developed and the adaptations which were necessary to incorporate these laws are detailed. Verification of the new contact force-displacement laws was achieved by simulating a quasi-static oblique contact and single particle oblique impact. Applications of the program to the simulation of large assemblies of particles is given, and the problems in undertaking quasi-static shear tests along with the results from two successful shear tests are described.

A method of fundamental solutions for two-dimensional heat conduction

We investigate an application of the method of fundamental solutions (MFS) to heat conduction in two-dimensional bodies, where the thermal diffusivity is piecewise constant. We extend the MFS proposed in Johansson and Lesnic [A method of fundamental solutions for transient heat conduction, Eng. Anal. Bound. Elem. 32 (2008), pp. 697–703] for one-dimensional heat conduction with the sources placed outside the space domain of interest, to the two-dimensional setting. Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate results can be obtained efficiently with small computational cost.

Relaxation procedures for an iterative MFS algorithm for the stable reconstruction of elastic fields from Cauchy data in two-dimensional isotropic linear elasticity

We investigate two numerical procedures for the Cauchy problem in linear elasticity, involving the relaxation of either the given boundary displacements (Dirichlet data) or the prescribed boundary tractions (Neumann data) on the over-specified boundary, in the alternating iterative algorithm of Kozlov et al. (1991). The two mixed direct (well-posed) problems associated with each iteration are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method, while the optimal value of the regularization parameter is chosen via the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases the iterative procedure at the point where the accumulation of noise becomes dominant and the errors in predicting the exact solutions increase, is also presented. The MFS-based iterative algorithms with relaxation are tested for Cauchy problems for isotropic linear elastic materials in various geometries to confirm the numerical convergence, stability, accuracy and computational efficiency of the proposed method.

Universal profile of the vortex condensate in two-dimensional turbulence

An inverse turbulent cascade in a restricted two-dimensional periodic domain creates a condensate—a pair of coherent system-size vortices. We perform extensive numerical simulations of this system and carry out theoretical analysis based on momentum and energy exchanges between the turbulence and the vortices. We show that the vortices have a universal internal structure independent of the type of small-scale dissipation, small-scale forcing, and boundary conditions. The theory predicts not only the vortex inner region profile, but also the amplitude, which both perfectly agree with the numerical data.

Lichen competition: two-dimensional warfare in slow motion

Lichens, a symbiotic association between a filamentous fungus and an alga, are often dominant in stressful environments such as the surfaces of rock and tree bark. Under these conditions, lichens experience extremes of temperature, moisture supply, and low availability of nutrients. As a consequence, lichens sequester a high proportion of their carbon production for stress resistance rather than for growth. Hence, as a group lichens are particularly slow growing organisms with many species growing at less than 2mm per year and some at less than 0.5mm per year. Whether or not competition occurs between lichen thalli in these communities is controversial. This article discusses the evidence that competition occurs between lichens on rock and tree bark and assesses whether competitive effects are likely to be important in structuring these communities.

Two-dimensional optical power distribution of side-out-coupled radiation from tilted FBGs in multimode fibre

Two-dimensional optical power distribution of the out-coupled radiation from UV-inscribed tilted fibre Bragg gratings (TFBGs) in multimode fibre (MMF) has been side-detected with high spatial/spectral resolution, showing a near-identical radiation mode profile to that measured from the fibre-end detection method. A comparative investigation of the radiation characteristics of TFBGs fabricated in singlemode fibre (SMF) and MMF clearly indicates that the radiation out-coupling is stronger and spatially more confined in MMF. The unique spatial-to-spectral encoding property of the coupling mechanism offers potential application in low-cost WDM devices.

A three dimensional numerical investigation of hillslope flow processes

Physically based distributed models of catchment hydrology are likely to be made available as engineering tools in the near future. Although these models are based on theoretically acceptable equations of continuity, there are still limitations in the present modelling strategy. Of interest to this thesis are the current modelling assumptions made concerning the effects of soil spatial variability, including formations producing distinct zones of preferential flow. The thesis contains a review of current physically based modelling strategies and a field based assessment of soil spatial variability. In order to investigate the effects of soil nonuniformity a fully three dimensional model of variability saturated flow in porous media is developed. The model is based on a Galerkin finite element approximation to Richards equation. Accessibility to a vector processor permits numerical solutions on grids containing several thousand node points. The model is applied to a single hillslope segment under various degrees of soil spatial variability. Such variability is introduced by generating random fields of saturated hydraulic conductivity using the turning bands method. Similar experiments are performed under conditions of preferred soil moisture movement. The results show that the influence of soil variability on subsurface flow may be less significant than suggested in the literature, due to the integrating effects of three dimensional flow. Under conditions of widespread infiltration excess runoff, the results indicate a greater significance of soil nonuniformity. The recognition of zones of preferential flow is also shown to be an important factor in accurate rainfall-runoff modelling. Using the results of various fields of soil variability, experiments are carried out to assess the validity of the commonly used concept of `effective parameters'. The results of these experiments suggest that such a concept may be valid in modelling subsurface flow. However, the effective parameter is observed to be event dependent when the dominating mechanism is infiltration excess runoff.

Comparison of BER estimation methods in numerical simulation of 40 Gbit/s RZ-DPSK transmission

This study analyzes the validity of different Q-factor models in the BER estimation in RZ-DPSK transmission at 40 Gb/s channel rate. The impact of the duty cycle of the carrier pulses on the accuracy of the BER estimates through the different models has also been studied.

Comparative analysis of BER estimation methods in numerical simulation of 40-Gb/s RZ-DPSK transmission with in-line SOAs

Applying direct error counting, we compare the accuracy and evaluate the validity of different available numerical approaches to the estimation of the bit-error rate (BER) in 40-Gb/s return-to-zero differential phase-shift-keying transmission. As a particular example, we consider a system with in-line semiconductor optical amplifiers. We demonstrate that none of the existing models has an absolute superiority over the others. We also reveal the impact of the duty cycle on the accuracy of the BER estimates through the differently introduced Q-factors.

Comparison of BER estimation methods in numerical simulation of 40 Gbit/s RZ-DPSK transmission

This study analyzes the validity of different Q-factor models in the BER estimation in RZ-DPSK transmission at 40 Gb/s channel rate. The impact of the duty cycle of the carrier pulses on the accuracy of the BER estimates through the different models has also been studied.

Comparative analysis of BER estimation methods in numerical simulation of 40-Gb/s RZ-DPSK transmission with in-line SOAs

Applying direct error counting, we compare the accuracy and evaluate the validity of different available numerical approaches to the estimation of the bit-error rate (BER) in 40-Gb/s return-to-zero differential phase-shift-keying transmission. As a particular example, we consider a system with in-line semiconductor optical amplifiers. We demonstrate that none of the existing models has an absolute superiority over the others. We also reveal the impact of the duty cycle on the accuracy of the BER estimates through the differently introduced Q-factors. © 2007 IEEE.

Two-dimensional optical power distribution of side-out-coupled radiation from tilted FBGs in multimode fibre

Two-dimensional optical power distribution of the out-coupled radiation from UV-inscribed tilted fibre Bragg gratings (TFBGs) in multimode fibre (MMF) has been side-detected with high spatial/spectral resolution, showing a near-identical radiation mode profile to that measured from the fibre-end detection method. A comparative investigation of the radiation characteristics of TFBGs fabricated in singlemode fibre (SMF) and MMF clearly indicates that the radiation out-coupling is stronger and spatially more confined in MMF. The unique spatial-to-spectral encoding property of the coupling mechanism offers potential application in low-cost WDM devices.

Two-dimensional viscoelastic discrete triangular system with negative-stiffness components

Potential applications of high-damping and high-stiffness composites have motivated extensive research on the effects of negative-stiffness inclusions on the overall properties of composites. Recent theoretical advances have been based on the Hashin-Shtrikman composite models, one-dimensional discrete viscoelastic systems and a two-dimensional nested triangular viscoelastic network. In this paper, we further analyze the two-dimensional triangular structure containing pre-selected negative-stiffness components to study its underlying deformation mechanisms and stability. Major new findings are structure-deformation evolution with respect to the magnitude of negative stiffness under shear loading and the phenomena related to dissipation-induced destabilization and inertia-induced stabilization, according to Lyapunov stability analysis. The evolution shows strong correlations between stiffness anomalies and deformation modes. Our stability results reveal that stable damping peaks, i.e. stably extreme effective damping properties, are achievable under hydrostatic loading when the inertia is greater than a critical value. Moreover, destabilization induced by elemental damping is observed with the critical inertia. Regardless of elemental damping, when the inertia is less than the critical value, a weaker system instability is identified.