Innovació tecnològica a la construcció subterrània

Report for the scientific sojourn carried out in the International Center for Numerical Methods in Engineering (CIMNE) –state agency – from February until November 2007. The work within the project Technology innovation in underground construction can be grouped into the following tasks: development of the software for modelling underground excavation based on the discrete element method - the numerical algorithms have been implemented in the computer programs and applied to simulation of excavation using roadheaders and TBM-s -; coupling of the discrete element method with the finite element method; development of the numerical model of rock cutting taking into account of wear of rock cutting tools -this work considers a very important factor influencing effectiveness of underground works -.

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

Development of an HTR-10 model in the Serpent reactor physics code

The use of exact coordinates of pebbles and fuel particles of pebble bed reactor modelling becoming possible in Monte Carlo reactor physics calculations is an important development step. This allows exact modelling of pebble bed reactors with realistic pebble beds without the placing of pebbles in regular lattices. In this study the multiplication coefficient of the HTR-10 pebble bed reactor is calculated with the Serpent reactor physics code and, using this multiplication coefficient, the amount of pebbles required for the critical load of the reactor. The multiplication coefficient is calculated using pebble beds produced with the discrete element method and three different material libraries in order to compare the results. The received results are lower than those from measured at the experimental reactor and somewhat lower than those gained with other codes in earlier studies.

Mathematical Modelling of the Dynamics of Granular Materials in a Rotating Cylinder

The current study is aimed at the development of a theoretical simulation tool based on Discrete Element Method (DEM) to 'interpret granular dynamics of solid bed in the cross section of the horizontal rotating cylinder at the microscopic level and subsequently apply this model to establish the transition behaviour, mixing and segregation.The simulation of the granular motion developed in this work is based on solving Newton's equation of motion for each particle in the granular bed subjected to the collisional forces, external forces and boundary forces. At every instant of time, the forces are tracked and the positions velocities and accelarations of each partcle is The software code for this simulation is written in VISUAL FORTRAN 90 After checking the validity of the code with special tests, it is used to investigate the transition behaviour of granular solids motion in the cross section of a rotating cylinder for various rotational speeds and fill fraction.This work is hence directed towards a theoretical investigation based on Discrete Element Method (DEM) of the motion of granular solids in the radial direction of the horizontal cylinder to elucidate the relationship between the operating parameters of the rotating cylinder geometry and physical properties ofthe granular solid.The operating parameters of the rotating cylinder include the various rotational velocities of the cylinder and volumetric fill. The physical properties of the granular solids include particle sizes, densities, stiffness coefficients, and coefficient of friction Further the work highlights the fundamental basis for the important phenomena of the system namely; (i) the different modes of solids motion observed in a transverse crosssection of the rotating cylinder for various rotational speeds, (ii) the radial mixing of the granular solid in terms of active layer depth (iii) rate coefficient of mixing as well as the transition behaviour in terms of the bed turnover time and rotational speed and (iv) the segregation mechanisms resulting from differences in the size and density of particles.The transition behaviour involving its six different modes of motion of the granular solid bed is quantified in terms of Froude number and the results obtained are validated with experimental and theoretical results reported in the literature The transition from slumping to rolling mode is quantified using the bed turnover time and a linear relationship is established between the bed turn over time and the inverse of the rotational speed of the cylinder as predicted by Davidson et al. [2000]. The effect of the rotational speed, fill fraction and coefficient of friction on the dynamic angle of repose are presented and discussed. The variation of active layer depth with respect to fill fraction and rotational speed have been investigated. The results obtained through simulation are compared with the experimental results reported by Van Puyvelde et. at. [2000] and Ding et at. [2002].The theoretical model has been further extended, to study the rmxmg and segregation in the transverse direction for different particle sizes and their size ratios. The effect of fill fraction and rotational speed on the transverse mixing behaviour is presented in the form of a mixing index and mixing kinetics curve. The segregation pattern obtained by the simulation of the granular solid bed with respect to the rotational speed of the cylinder is presented both in graphical and numerical forms. The segregation behaviour of the granular solid bed with respect to particle size, density and volume fraction of particle size has been investigated. Several important macro parameters characterising segregation such as mixing index, percolation index and segregation index have been derived from the simulation tool based on first principles developed in this work.

Análise da forma da partícula na aplicação do Método dos Elementos Discretos (DEM)

O presente trabalho analisa o efeito da forma da partícula por meio do método dos Elementos Discretos (DEM). São analisadas quatro partículas com formas e esfericidades diferentes – uma esférica e três não-esféricas. O ângulo de repouso é o principal instrumento de avaliação dos escoamentos. A análise inicia com a calibração da partícula esférica usando dados da literatura, em seguida, estes parâmetros calibrados são empregados para a simulação dos escoamentos das partículas não-esféricas. São feitas comparações dos esforços computacionais e, estas informações, são usadas para verificar as vantagens que a partícula esférica oferece sobre as outras três formas. Nesse cenário, procedimentos foram desenvolvidos para ajudar no processo de calibração dos ângulos de repousos, tendo como base, o conhecimento da sensibilidade de alguns parâmetros DEM. Os resultados mostram a influência das partículas não-esféricas e, principalmente, que é possível obter com a partícula esférica, escoamentos similares aos das partículas não-esféricas, com ganho computacional.

Interação CFD-DEM em fluidização: aplicação para o setor mineral

A fluidização de partículas é amplamente utilizada na indústria, principalmente devido às altas taxas de transferência de calor e massa entre as fases. O acoplamento entre a Dinâmica dos Fluidos Computacional (CFD – Computational Fluid Dynamics) e o Método dos Elementos Discretos (DEM – Discrete Element Method) tem se tornado atrativo para a simulação de fluidização, já que nesse caso o movimento das partículas é analisado de forma mais direta do que em outros tipos de abordagens. O grande problema do acoplamento CFD-DEM é a alta exigência computacional para rastrear todas as partículas do sistema, o que leva ao uso de estratégias de redução do tempo de simulação que em caso de utilização incorreta podem comprometer os resultados. O presente trabalho trata da aplicação do acoplamento CFD-DEM na análise de fluidização de alumina, que é um problema importante para o setor mineral. Foram analisados diversos parâmetros capazes de influenciar os resultados e o tempo de simulação como os passos de tempo, os modelos de arrasto, a distribuição granulométrica das partículas, a constante de rigidez, a utilização de partículas representativas com tamanho maior que o das partículas reais, etc. O modelo de força de interação DEM utilizado foi o modelo de mola e amortecedor lineares (LSD – Linear Spring Dashpot). Todas as simulações foram realizadas com o software ANSYS FLUENT 14.5 e os resultados obtidos foram comparados com dados experimentais e da literatura. Tais resultados permitiram comprovar a capacidade do modelo linear LSD em predizer o comportamento global de leitos de alumina e reduzir o tempo de simulação, desde que os parâmetros do modelo sejam definidos de forma adequada.

Simulation von Förderprozessen bei Vibrationsförderanlagen

Die erzielbare Fördergeschwindigkeit bei Vibrationsförderern hängt maßgeblich von der Bewegungsfunktion des Förderorganes ab. Für die gezielte Simulation dieser Anlagen mittels der diskreten Elemente Methode (DEM) ist es notwendig die geometrisch vernetzen Förderorgannachbildungen mit praxisrelevanten Bewegungsfunktionen zu beaufschlagen. Der Artikel beschreibt die Einbindung dieser Bewegungsfunktionen in die quellenoffene DEM-Software LIGGGHTS. Während des Simulationsprozesses wird eine Bewegung vernetzter CAD-Modelle durch trigonometrische Reihen ermöglicht.

Simulation der peristaltischen Förderung von Stückgütern als Schüttgut

Im Beitrag wird ein neuartiges Förderprinzip zur federnden Aufnahme und zum Transport von massenhaft anfallenden Paketstrukturen vorgestellt. Das Förderprinzip beruht auf einem flächigen Tragmittel in Form eines veränderbaren, elastischen Verbundes von kleinskaligen Fördermodulen. Das konzipierte Transportprinzip mit peristaltischen Eigenschaften soll entstehende Staus der Pakete schnell auflösen und eine dedizierte Steuerung von Teilmengen zulassen, um den erforderlichen Durchsatz innerhalb eines Materialflusssystems zu erreichen. Diese Lösung ermöglicht eine sinnvolle Verknüpfung von Wirkprinzipien der Schüttgut- und Stückgutförderung zur Aufnahme und Fortbewegung von Pakete als Schüttgut. Die Grundfunktionalität des Förderkonzepts wird durch die numerische Simulation auf Basis der Diskrete Elemente Methode sowie der Mehrkörpersimulation überprüft.

Determination of the mechanical properties of corn grains and olive fruits required in DEM simulations

Discrete element method (DEM) is a numerical technique widely used for simulating the mechanical behavior of granular materials involved in many food and agricultural industry processes. Additionally, this technique is also a powerful tool to understand many complex phenomena related to the mechanics of granular materials. However, to make use of the potential of this technique it is necessary to develop DEM models capable of representing accurately the reality. For that, among some other questions, it is essential that the values of the microscopic material properties used to define the numerical model are accurately determined.

Determinación del coeficiente de rozamiento grano-pared en silos de pared corrugada mediante el método de los elementos discretos

El comportamiento de los materiales granulares almacenados en silos se ve afectado por varios parámetros, tanto aquellos característicos del material como de la geometría del silo. La determinación del coeficiente de rozamiento pared-partícula es uno de los parámetros de mayor importancia, siendo habituales en su determinación el uso de ensayos de corte directo. En el presente trabajo se estudia dicho coeficiente para el supuesto de una pared corrugada, el cual, teóricamente, debe representar a un valor efectivo que depende tanto del coeficiente de rozamiento grano-pared (para el caso de una pared lisa) y el ángulo de rozamiento interno del material. La determinación del rozamiento efectivo se ha realizado a través de la simulación por elementos discretos de un ensayo de corte sobre una pared corrugada. Los valores obtenidos han sido comparados con las prescripciones expuestas en la normativa vigente. La potencialidad del método de los elementos discretos permite el estudio de diversas configuraciones geométricas de la pared corrugada de los silos sin necesidad de realizar ensayos de laboratorio. Esto permitirá, en trabajos sucesivos, investigar la influencia de muy diversos factores en el valor del rozamiento efectivo grano-pared en este tipo de paredes. The mechanical behaviour of granular materials stored in silos is affected by numerous parameters, some of them being related to the characteristic of the stored materials and others to the geometry of the silo. The determination of the particle-wall friction coefficient (or wall friction) is of great importance and direct shear tests are usually conducted in order to obtain its value. In this work this variable is analysed for the case of a corrugated wall. This value is expected to be an effective value between the particle-wall friction coefficient (obtained for a flat wall) and the internal friction coefficient of the material under study. The effective wall friction determination has been carried out by using a discrete element model to simulate a direct shear test on a corrugated wall. The values obtained have been compared with prescriptions given in the current standards. The potential of the discrete element method allows different geometries of the silo corrugated wall to be considered without the necessity of developing laboratory tests. In future works it will allow the influence of numerous parameters on the effective wall friction in corrugated walls to be studied.

Numerical and experimental investigation of structural stiffness influence on ratcheting convection cell in granular soils under cyclic loading

Lateral cyclic loaded structures in granular soils can lead to an accumulation of irreversible strains by changing their mechanical response (densification) and forming a closed convective cell in the upper layer of the bedding. In the present thesis the convective cell dimension, formation and grain migration inside this closed volume have been studied and presented in relation to structural stiffness and different loads. This relation was experimentally investigated by applying a cyclic lateral force to a scaled flexible vertical element embedded in dry granular soil. The model was monitored with a camera in order to derive the displacement field by means of the PIV technique. Modelling large soil deformation turns out to be difficult, using mesh-based methods. Consequently, a mesh-free approach (DEM) was chosen in order to investigate the granular flow with the aim of extracting interesting micromechanical information. In both the numerical and experimental analyses the effect of different loading magnitudes and different dimensions of the vertical element were considered. The main results regarded the different development, shape and dimensions of the convection cell and the surface settlements. Moreover, the Discrete Element Method has proven to give satisfactory results in the modelling of large deformation phenomena such as the ratcheting convective cell.

Using DEM to model ore breakage within a pilot scale SAG mill

The best accepted method for design of autogenous and semi-autogenous (AG/SAG) mills is to carry out pilot scale test work using a 1.8 m diameter by 0.6 m long pilot scale test mill. The load in such a mill typically contains 250,000-450,000 particles larger than 6 mm, allowing correct representation of more than 90% of the charge in Discrete Element Method (DEM) simulations. Most AG/SAG mills use discharge grate slots which are 15 mm or more in width. The mass in each size fraction usually decreases rapidly below grate size. This scale of DEM model is now within the possible range of standard workstations running an efficient DEM code. This paper describes various ways of extracting collision data front the DEM model and translating it into breakage estimates. Account is taken of the different breakage mechanisms (impact and abrasion) and of the specific impact histories of the particles in order to assess the breakage rates for various size fractions in the mills. At some future time, the integration of smoothed particle hydrodynamics with DEM will allow for the inclusion of slurry within the pilot mill simulation. (C) 2004 Elsevier Ltd. All rights reserved.

Effect of granule morphology on breakage behaviour during compression

In the area of dry particle breakage, Discrete Element Method (DEM) simulations have been widely used to analyse the sensitivity of various physical parameters to the behaviour of agglomerates during breakage. This paper looks at the effect of agglomerate shape and structure on the mechanisms and extent of breakage of dry agglomerates under compressive load using DEM simulations. In the simulations, a spherical-shaped agglomerate produced within the DEM code is compared with an irregularly shaped agglomerate, whose structure is that of an actual granule that was characterised with X-ray microtomography (muCT). Both agglomerates have identical particle size distribution, coordination number and surface energy values, with only the agglomerate shape and structure differing between the two. The work here details the breakage behaviour with a number of traditional DEM output parameters (i.e., contact/cluster distributions) with showing vastly different behaviour between the two agglomerates. (C) 2004 Elsevier B.V. All rights reserved.

A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics

The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to Moore's law, i.e., to double every 18-24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures.

Influence of charge size distribution on net-power draw of tumbling mill based on DEM modelling

Modelling and optimization of the power draw of large SAG/AG mills is important due to the large power draw which modern mills require (5-10 MW). The cost of grinding is the single biggest cost within the entire process of mineral extraction. Traditionally, modelling of the mill power draw has been done using empirical models. Although these models are reliable, they cannot model mills and operating conditions which are not within the model database boundaries. Also, due to its static nature, the impact of the changing conditions within the mill on the power draw cannot be determined using such models. Despite advances in computing power, discrete element method (DEM) modelling of large mills with many thousands of particles could be a time consuming task. The speed of computation is determined principally by two parameters: number of particles involved and material properties. The computational time step is determined by the size of the smallest particle present in the model and material properties (stiffness). In the case of small particles, the computational time step will be short, whilst in the case of large particles; the computation time step will be larger. Hence, from the point of view of time required for modelling (which usually corresponds to time required for 3-4 mill revolutions), it will be advantageous that the smallest particles in the model are not unnecessarily too small. The objective of this work is to compare the net power draw of the mill whose charge is characterised by different size distributions, while preserving the constant mass of the charge and mill speed. (C) 2004 Elsevier Ltd. All rights reserved.