Sensores capacitivos são amplamente utilizados em vários sistemas de medida, tanto na indústria como em laboratórios de pesquisa. Estes dispositivos medem variadas grandezas físicas, como deslocamento, força, pressão, densidade e concentração de substâncias em escoamentos bifásicos, entre outros. Sua aplicação vêm sendo utilizada há algum tempo de forma relativamente elementar, necessitando ainda de estudos específicos para entender a fundo o fenômeno de transferência de carga. A sonda proposta neste trabalho visa melhorar a aquisição do sinal capacitivo e obter uma medida precisa de fração volumétrica em condições bifásicas e trifásicas. Portanto, é necessário um circuito transdutor que transmita de forma precisa a capacitância registrada entre os eletrodos dos sensores até uma central de aquisição. Neste trabalho buscamos desenvolver um modelo de sonda para medir frações volumétricas in situ. Para tanto, apresentamos as etapas de calibração dos sensores, através de válvulas de fechamento rápido, e validamos a técnica sob determinadas condições e variados padrões de escoamento. Dois sensores de geometrias distintas (anéis e hélices) foram utilizados na composição da sonda capacitiva não intrusiva. Propõe-se a solução do sistema linear para medição direta de fração volumétrica em escoamento trifásico. / Capacitive sensors are widely employed in many measurement systems in both industry and laboratory activities. These devices measure several physical quantities such as displacement, forces, pressure, density and concentration of substances in two-phase flows. Their application in industry is still insipient, in part because specific studies are still necessary to properly understand the phenomenon of charge transfer. The goal of this research is the development of a non-intrusive double-sensor probe to on-line measurement of in-situ volumetric fractions in two and three-phase flows. Therefore, we present the calibration work, through quick-closing valves, and validation techniques under certain conditions and for different flows patterns. A transducer circuit to transmit the capacitance sign accurately from the sensors to the central acquisition was developed. Two sensors with different geometries (rings and helix) comprise the double-sensor capacitive probe. Finally, we propose the solution of a linear system for the direct measurement of in-situ volumetric fraction in three-phase flow.

A economia mundial é fortemente dependente da disponibilidade de óleo, no entanto as reservas de óleo leve, mais atrativas ao mercado, tendem nos próximos anos ao esgotamento, trazendo à luz as reservas de óleo pesado para atender as demandas do mercado. No entanto, os processos associados ao transporte e processamento desses óleos viscosos consomem muita energia, impondo um grande desafio para a indústria do petróleo. Nesse contexto, se transportar óleos viscosos com água através de um escoamento bifásico em padrão Core Annular Flow (CAF) é muito promissor. Nessa prática, uma película anular de água envolve o núcleo que contém o óleo, minimizando o contato entre este e a parede do duto e reduzindo as perdas de energia por atrito. Visou o presente estudo mapear os padrões de escoamento de óleo pesado com água em dutos horizontais e verticais, mensurar os fatores de redução de potência e do diferencial de pressão associados à um trecho reto e à uma válvula de gaveta aberta e determinar o holdup para avaliar a eficiência do CAF. Para atingir esses objetivos, uma unidade de bancada foi construída, consistindo de tanques de armazenamento e de separação da mistura. Interligando os tanques, foram dispostos dutos transparentes com 27 mm de diâmetro interno e 8 m de extensão, contendo duas seções horizontais e uma vertical para análise e diversos acessórios hidráulicos. O óleo utilizado possuía viscosidade de 3200 cP a 22 °C e densidade de 945 kg/m3 . Os resultados experimentais demonstraram haver diversas configurações de fluxo segundo o posicionamento do duto, e que a diferença de densidade dos fluidos descentraliza o óleo nos escoamentos horizontais. Além disso, o trecho vertical intercalado entre os trechos horizontais se comportou como um retificador de fluxo, melhorando os índices do CAF. Foi determinado um fator de redução de perda de pressão máximo de 250 vezes para o trecho reto e de 12 vezes para a válvula de gaveta. O fator global máximo de redução de potência foi mensurado em 2,2 vezes. Concluiu-se que, para um trecho reto, o projeto de instalações hidráulicas para o escoamento bifásico deve considerar 15 % a mais de perda de carga em relação ao escoamento de água pura. Já para a válvula de gaveta, esse fator deve ser de 700 %. / The world economy is strongly dependent on the availability of oil, however, light oil reserves, more market-oriented, tend to deplete in coming years, bringing to light the heavy oil reserves to meet the demands of the market. However, the processes associated with the transportation and production of these viscous oils consume a lot of energy and pose a great challenge for the oil industry. In this context, transporting viscous oils through a liquid-liquid two-phase oil-water flow in a core annular pattern (CAF) is very promising. In this method, an annular water film surrounds the oil-containing core, minimizing its contact with the pipe wall and reducing energy losses by friction. The aim of the present work was to map the flow patterns of a biphasic oil-water flow in horizontal and vertical pipes, to evaluate the overall energy savings provided by the CAF technique, to measure the pressure gradient reduction factor along the pipe and in a gate valve and to determine the holdup as a way of assessing the energy efficiency of the biphasic oil-water flow transport. To achieve these goals, an experimental facility was built and it consisted of cargo and separation tanks. Connecting these tanks, approximately 8 m of 27mm-ID clear transparent PVC pipes, two horizontal and one vertical sections and various hydraulic fittings. It was used in the tests lubricating oil with 3200 cP and 945 kg/m3 at 22°C and distilled water. The experimental results showed that there are several flow configurations according to the pipe positioning, and that the oil and the water density difference decentralizes the oil core in horizontal flows. In addition, the vertical section placed between the two horizontal ones behaved like a flow rectifier, improving the core annular flow energy savings basis. An average pressure gradient reduction factor of 250 times in a straight pipe and of 12 times in a gate valve was determined. The maximum overall power reduction factor was measured as being 2.2 times. It was concluded that the design of a hydraulic installation to transport heavy oil with water in a core annular pattern should consider a pressure drop increase by a factor of 15 % in a straight pipe and by a factor of 700 % in a gate valve based on the monophasic water transport at similar flow rates.

As Vibrações Induzidas por Escoamentos (VIEs) são dos problemas dinâmicos mais importantes no projeto de trocadores de calor de tipo carcaça e tubos. Esta interação fluido - estrutura pode gerar vibrações de grande amplitude dos tubos ou partes internas do equipamento, o que pode provocar o impacto entre componentes e até falha por fadiga. Várias bancadas experimentais têm sido construídas para o estudo desta interação fluido - estrutura, no entanto, alguns mecanismos de vibração, principalmente aqueles relacionados com escoamentos bifásicos, não são completamente compreendidos. Portanto, neste trabalho, o projeto de uma estrutura dinâmica para o estudo das vibrações induzidas por escoamento bifásico é apresentado. Esta estrutura é composta por um sistema de fios de aço tensionados que permitem que a primeira frequência de ressonância do tubo seja sintonizada. O tubo instrumentado foi instalado numa seção de testes de tubos rígidos em configuração triangular normal com uma razão de passo transversal e diâmetro de 1,26. Os resultados experimentais foram contrastados com uma modelagem teórica desenvolvida para esta estrutura dinâmica. Pôde-se verificar que a primeira frequência de ressonância e seu fator de amortecimento coincidem com os valores estimados pelo modelo, enquanto que a segunda frequência de ressonância apresenta valor teórico 6,6% maior do que o experimental. A seguir, um estudo experimental no ar foi desenvolvido para o reconhecimento dos modos de vibração e as frequências de ressonância da estrutura dinâmica. O fator de amortecimento no ar também foi estimado, o método de Kennedy - Pancu combinado com o algoritmo ERA foi usado para este fim. Assim que a primeira frequência de ressonância foi reconhecida e sintonizada, testes em escoamentos de água foram realizados. A influência da velocidade do escoamento monofásico de água na amplitude de vibração e nas características da Densidade Espectral (PSD) de aceleração foi analisada. Finalmente, testes em escoamento bifásico foram feitos, a sensibilidade do tubo ao escoamento bifásico fica evidente quando analisadas as amplitudes nas PSDs. Além disso, a influência da fração de vazio na amplitude de vibração e no fator de amortecimento foi revisada, no entanto, resultados de estudos anteriores não foram conferidos pelo alto amortecimento concentrado introduzido por alguns elementos construtivos utilizados. / Flow - Induced Vibration (FIV) is the most critical dynamic issue in the design of heat exchangers. This fluid - structure phenomenon may generate high amplitude vibration of tubes or structural parts, which may lead to impacts between internal components or even failure due to fatigue. Many test benches have been constructed to study this fluid - structure interactions, however, some vibration mechanisms, mostly those related to multiphase flow, are not yet fully understood. Therefore, in this work, the design of a dynamic structure for the study of multiphase flow induced vibration is presented. This structure is composed by a system of tensioned piano wires that allow the first natural frequency of the instrumented tube to be calibrated. The tube was installed in a rigid bundle configured in a normal triangular pattern with pitch-todiameter ratio of 1.26. Then, an experimental study in air environment was performed aiming at addressing the mode shapes and resonance frequencies of the dynamic structure. The damping ratio in air environment was also estimated, the Kennedy - Pancu method in combination with the Eigensystem Realization Algorithm have been used to that purpose. Experimental results in air were compared with a theoretical model. The first resonance frequency and its damping factor presented a good agreement with the model, while the second resonance frequency deviated some 6,6% from the predicted value. Tests in water and multiphase flow were also performed. The influence of water flow velocity on vibration amplitudes and on the Power Spectrum Densities (PSDs) has been analyzed. Finally, tests for two - phase water/air flow were carried out; the tube sensitivity to multiphase flow is evident when the vibration level in acceleration PSDs is by analyzed. Also, the influence of void fraction on vibration amplitude and damping ratio have been addressed, however, previous studies could not be validated due to the high damping introduced by some constructive elements of the structure.

A economia mundial é fortemente dependente da disponibilidade de óleo, no entanto as reservas de óleo leve, mais atrativas ao mercado, tendem nos próximos anos ao esgotamento, trazendo à luz as reservas de óleo pesado para atender as demandas do mercado. No entanto, os processos associados ao transporte e processamento desses óleos viscosos consomem muita energia, impondo um grande desafio para a indústria do petróleo. Nesse contexto, se transportar óleos viscosos com água através de um escoamento bifásico em padrão Core Annular Flow (CAF) é muito promissor. Nessa prática, uma película anular de água envolve o núcleo que contém o óleo, minimizando o contato entre este e a parede do duto e reduzindo as perdas de energia por atrito. Visou o presente estudo mapear os padrões de escoamento de óleo pesado com água em dutos horizontais e verticais, mensurar os fatores de redução de potência e do diferencial de pressão associados à um trecho reto e à uma válvula de gaveta aberta e determinar o holdup para avaliar a eficiência do CAF. Para atingir esses objetivos, uma unidade de bancada foi construída, consistindo de tanques de armazenamento e de separação da mistura. Interligando os tanques, foram dispostos dutos transparentes com 27 mm de diâmetro interno e 8 m de extensão, contendo duas seções horizontais e uma vertical para análise e diversos acessórios hidráulicos. O óleo utilizado possuía viscosidade de 3200 cP a 22 °C e densidade de 945 kg/m3 . Os resultados experimentais demonstraram haver diversas configurações de fluxo segundo o posicionamento do duto, e que a diferença de densidade dos fluidos descentraliza o óleo nos escoamentos horizontais. Além disso, o trecho vertical intercalado entre os trechos horizontais se comportou como um retificador de fluxo, melhorando os índices do CAF. Foi determinado um fator de redução de perda de pressão máximo de 250 vezes para o trecho reto e de 12 vezes para a válvula de gaveta. O fator global máximo de redução de potência foi mensurado em 2,2 vezes. Concluiu-se que, para um trecho reto, o projeto de instalações hidráulicas para o escoamento bifásico deve considerar 15 % a mais de perda de carga em relação ao escoamento de água pura. Já para a válvula de gaveta, esse fator deve ser de 700 %. / The world economy is strongly dependent on the availability of oil, however, light oil reserves, more market-oriented, tend to deplete in coming years, bringing to light the heavy oil reserves to meet the demands of the market. However, the processes associated with the transportation and production of these viscous oils consume a lot of energy and pose a great challenge for the oil industry. In this context, transporting viscous oils through a liquid-liquid two-phase oil-water flow in a core annular pattern (CAF) is very promising. In this method, an annular water film surrounds the oil-containing core, minimizing its contact with the pipe wall and reducing energy losses by friction. The aim of the present work was to map the flow patterns of a biphasic oil-water flow in horizontal and vertical pipes, to evaluate the overall energy savings provided by the CAF technique, to measure the pressure gradient reduction factor along the pipe and in a gate valve and to determine the holdup as a way of assessing the energy efficiency of the biphasic oil-water flow transport. To achieve these goals, an experimental facility was built and it consisted of cargo and separation tanks. Connecting these tanks, approximately 8 m of 27mm-ID clear transparent PVC pipes, two horizontal and one vertical sections and various hydraulic fittings. It was used in the tests lubricating oil with 3200 cP and 945 kg/m3 at 22°C and distilled water. The experimental results showed that there are several flow configurations according to the pipe positioning, and that the oil and the water density difference decentralizes the oil core in horizontal flows. In addition, the vertical section placed between the two horizontal ones behaved like a flow rectifier, improving the core annular flow energy savings basis. An average pressure gradient reduction factor of 250 times in a straight pipe and of 12 times in a gate valve was determined. The maximum overall power reduction factor was measured as being 2.2 times. It was concluded that the design of a hydraulic installation to transport heavy oil with water in a core annular pattern should consider a pressure drop increase by a factor of 15 % in a straight pipe and by a factor of 700 % in a gate valve based on the monophasic water transport at similar flow rates.

Les simulations numériques à l'échelle du pore sont fréquemment utilisées pour étudier le comportement des écoulements multiphasiques largement rencontrées dans phénomènes naturels et applications industrielles. Dans ce travail, la morphologie de structures liquides et l'action capillaire sont examinées à l'échelle des pores par la méthode de Boltzmann sur réseau (LBM) à plusieurs composants selon le modèle de Shan-Chen. Les résultats numériques obtenus sont en bon accord avec les solutions théoriques. Les simulations numériques sont étendues à microstructures complexes au-delà du régime pendulaire.La LBM a été utilisée pour modéliser l'écoulement multiphasique à travers un milieu poreux idéalisé dans des conditions de drainage primaire quasi-statique. Les simulations LBM ont fourni une excellente description du déplacement de l'interface fluide-fluide à travers les grains. Pendant le drainage, les simulations LBM sont capables de reproduire la déconnexion d'une phase dans le milieu granulaire sous la forme de ponts pendulaires ou structures liquides complexes. Malheureusement, le temps de calcul nécessaire pour ce type de simulations est assez élevé. Afin d’optimiser les ressources de calcul, nous présentons un modèle 2D (modèle Throat-Network) basé sur des solutions analytiques pour décrire l'écoulement biphasique à travers un ensemble de disques dans un temps de calcul très réduit, donc le modèle 2D est susceptible de remplacer les simulations LBM lorsque les ressources de calcul sont limitées. L'approche souligne l'importance de simuler le problème a l'échelle de la gorge du pore pour obtenir les relations volume - pression capillaire locales. Le modèle Throat-Network est un point de départ pour le modèle hybride proposé pour résoudre les problèmes en 3D. Le modèle hybride combine l’efficacité de l’approche réseau de pores et la précision du LBM à l’échelle des pores. Le modèle hybride est basé sur la décomposition de l’échantillon en petits sous-domaines, dans lesquels des simulations LBM sont effectuées pour déterminer les propriétés hydrostatiques principales (pression capillaire d'entrée, courbe de drainage primaire et morphologie du liquide pour chaque gorge du pore). Malgré la réduction significative des temps de calcul obtenus avec le modèle hybride, le temps n’est pas négligeable et les modélisations numériques d'échantillons de grandes tailles ne sont pas réalistes. Les approximations données par les méthodes Incircle et MS-P, qui prédisent les propriétés hydrostatiques, sont comparées à celles de LBM et du modèle hybride. / Numerical simulations at the pore scale are a way to study the behavior of multiphase flows encountered in many natural processes and industrial applications. In this work, liquid morphology and capillary action are examined at the pore-scale by means of the multicomponent Shan-Chen lattice Boltzmann method (LBM). The accuracy of the numerical model is first contrasted with theoretical solutions. The numerical results are extended to complex microstructures beyond the pendular regime.The LBM has been employed to simulate multiphase flow through idealized granular porous media under quasi-static primary drainage conditions. LBM simulations provide an excellent description of the fluid-fluid interface displacement through the grains. Additionally, the receding phase trapped in the granular media in form of pendular bridges or liquid clusters is well captured. Unfortunately, such simulations require a significant computation time. A 2D model (Throat-Network model) based on analytical solutions is proposed to mimic the multiphase flow with very reduced computation cost, therefore, suitable to replace LBM simulations when the computation resources are limited. The approach emphasizes the importance of simulating at the throat scale rather than the pore body scale in order to obtain the local capillary pressure - liquid content relationships. The Throat-Network model is a starting point for the a hybrid model proposed to solve 3D problems. The hybrid model combines the efficiency of the pore-network approach and the accuracy of the LBM at the pore scale to optimize the computational resources. The hybrid model is based on the decomposition of the granular assembly into small subsets, in which LBM simulations are performed to determine the main hydrostatic properties (entry capillary pressure, capillary pressure - liquid content relationship and liquid morphology for each pore throat). Despite the reduction of computation time, it is still not negligible and not affordable for large granular packings. Approximations by the Incircle and the MS-P method, which predict hydrostatic properties, are contrasted with the results provided by LBM and the hybrid model. Relatively accurate predictions are given by the approximations.

The development, analysis, and implementation of efficient methods to solve algebraic systems of equations are main research directions in the field of numerical simulation and are the focus of this thesis. Due to their lesser demands for computer resources, iterative solution methods are the choice to make, when very large scale simulations have to be performed. To improve their efficiency, iterative methods are combined with proper techniques to accelerate convergence. A general technique to do this is to use a so-called preconditioner. Constructing and analysing various preconditioning methods has been an active field of research already for decades. Special attention is devoted to the class of the so-called optimal order preconditioners, that possess both optimal convergence rate and optimal computational complexity. The preconditioning techniques, proposed and studied in this thesis, utilise the block structure of the underlying matrices, and lead to methods that are of optimal order. In the first part of the thesis, we construct an Algebraic MultiLevel Iteration (AMLI) method for systems arising from discretizations of parabolic problems, using Crouzeix-Raviart finite elements. The developed AMLI method is based on an approximated block factorization of the original system matrix, where the partitioning is associated with a sequence of nested discretization meshes. In the second part of the thesis we develop solution methods for the numerical simulation of multiphase flow problems, modelled by the Cahn-Hilliard (C-H) equation. We consider the discrete C-H problem, obtained via finite element discretization in space and implicit schemes in time. We propose techniques to precondition the Jacobian of the discrete nonlinear system, based on its natural two-by-two block structure. The preconditioners are used in the framework of inexact Newton methods. We develop two nonlinear solution algorithms for the Cahn-Hilliard problem. Both lead to efficient optimal order methods. One of the main advantages of the proposed methods is that they are implemented using available software toolboxes for both sequential and distributed execution. The theoretical analysis of the solution methods presented in this thesis is combined with numerical studies that confirm their efficiency.

Surface energy plays a major role in many phenomena that are important in technological and industrial processes, for example in wetting, grain growth and sintering. In this thesis, such surface-energy driven processes are studied by means of the phase-field method. The phase-field method is often used to model mesoscale microstructural evolution in materials. It is a diffuse interface method, i.e., it considers the surface or phase boundary between two bulk phases to have a non-zero width with a gradual variation in physical properties such as energy density, composition and crystalline structure. Neck formation and coarsening are two important diffusion-controlled features in solid-state sintering and are studied using our multiphase phase-field method. Inclusion of Navier-Stokes equation with surface-tension forces and convective phase-field equations into the model, enables simulation of reactive wetting and liquid-phase sintering. Analysis of a spreading liquid on a surface is investigated and is shown to follow the dynamics of a known hydrodynamic theory. Analysis of important capillary phenomena with wetting and motion of two particles connected by a liquid bridge are studied in view of important parameters such as contact angles and volume ratios between the liquid and solid particles. The interaction between solute atoms and migrating grain boundaries affects the rate of recrystallization and grain growth. The phenomena is studied using a phase-field method with a concentration dependent double-well potential over the phase boundary. We will show that with a simple phase-field model it is possible to model the dynamics of grain-boundary segregation to a stationary boundary as well as solute drag on a moving boundary. Another important issue in phase-field modeling has been to develop an effective coupling of the phase-field and CALPHAD methods. Such coulping makes use of CALPHAD's thermodynamic information with Gibbs energy function in the phase-field method. With the appropriate thermodynamic and kinetic information from CALPHAD databases, the phase-field method can predict mictrostructural evolution in multicomponent multiphase alloys. A phase-field model coupled with a TQ-interface available from Thermo-Calc is developed to study spinodal decomposition in FeCr, FeCrNi and TiC-ZrC alloys. / QC 20100622

Coupled Thermo-Hydro-Mechanical (THM) behavior in geological porous media has been a subject of great interest in many geoengineering disciplines. Many attempts have been made to develop numerical prediction capabilities associated with topics such as the movement of pollutant plumes, gas injection, energy storage, geothermal energy extraction, and safety assessment of repositories for radioactive waste and spent nuclear fuel. This thesis presents a new numerical modeling approach and a new computer code for simulating coupled THM behavior in geological porous media in general, and compacted bentonite clays in particular, as buffer materials in underground radioactive waste repositories. New governing equations were derived according to the theory of mixtures, considering interactions among solid-phase deformation, flows of water and gases, heat transport, and phase change of water. For three-dimensional problems, eight governing equations were formulated to describe the coupled THM processes. A new thermal conductivity model was developed to predict the thermal conductivity of geological porous media as composite mixtures. The proposed model considers the combined effects of solid mineral composition, temperature, liquid saturation degree, porosity and pressure on the effective thermal conductivity of the porous media. The predicted results agree well with the experimental data for MX80 bentonite. A new water retention curve model was developed to predict the suction-saturation behavior of the geological porous media, as a function of suction, effective saturated degree, temperature, porosity, pore-gas pressure, and the rate of saturation degree change with time. The model was verified against experimental data of the FEBEX bentonite, with good agreement between measured and calculated results. A new finite element code (ROLG) was developed for modeling fully coupled thermo-hydro-mechanical processes in geological porous media. The new code was validated against several analytical solutions and experiments, and was applied to simulate the large scale in-situ Canister Retrieval Test (CRT) at Äspö Hard Rock Laboratory, SKB, Sweden, with good agreement between measured and predicted results. The results are useful for performance and safety assessments of radioactive waste repositories. / QC20100720 / THERESA

Integrated theoretical and experimental analysis is carried out in this work to investigate the fundamental failure mechanisms and flow patterns involved in the process of fluid injection into dense granular media. The experimental work is conducted with aqueous glycerin solutions, utilizing a novel setup based on a Hele-Shaw cell filled with dense dry sand. The two dimensional nature of the setup allows direct visualization and imaging analysis of the real-time fluid and grain kinematics. The experimental results reveal that the fluid flow patterns show a transition from simple radial flow to a ramified morphology while the granular media behaviors change from that of rigid porous media to localized failure that lead to development of fluid channels. Based on the failure/flow patterns, four distinct failure/flow regimes can be identified, namely, (i) a simple radial flow regime, (ii) an infiltration-dominated regime, (iii) a grain displacement-dominated regime, and (iv) a viscous fingering-dominated regime. These distinct failure/flow regimes emerge as a result of competition among various energy dissipation mechanisms, namely, viscous dissipation through infiltration, dissipation due to grain displacements, and viscous dissipation through flow in thin channels and can be classified based on the characteristic times associated with fluid injection, hydromechanical coupling and viscoelastoplasticity. The injection process is also analyzed numerically using the discrete element method (DEM) coupled with two fluid flow scheme, a fixed coarse grid scheme based on computational fluid dynamics (CFD) and a pore network modeling scheme. The numerical results from the two complementary methods reproduce phenomena consistent with the experimental observations and justify the concept of associating the displacement regimes with the partition among energy dissipation mechanisms. The research in this work, though fundamental in nature, will have direct impacts on many engineering problems in civil, environmental and petroleum engineering such as ground improvement, environmental remediation and reservoir stimulation.

Gas-Partikel-Stroemungen finden sich in weiten Bereichen der Energie- und Verfahrenstechnik. Beispiele fuer haeu- fig anzutreffende Problemstellungen sind der Transport, die Separation oder die Injektion eines Gemisches aus festen Partikeln und einem Traegergas. Fuer die numerische Simulation solcher disperser Mehr- phasenstroemungen hat sich das Lagrange-Verfahren als besonders geeignet erwiesen. Andererseits stellt die An- wendung dieses Berechnungsverfahrens hoechste Anforderun- gen an die Ressourcen der verwendeten Rechner. Dies gilt im besonderen Masse fuer die Simulation von Stroemungen mit einer moderaten bis hohen Partikelbeladung, in denen die Partikel-Partikel-Kollisionen einen grossen Einfluss auf das Stroemungsverhalten haben. Um das grosse Leistungspotential, das heutige massiv par- allele Hochleistungsrechner bieten, effizient zu nutzen, wurden im Rahmen dieser Arbeit parallele Simulationsalgo- rithmen fuer die numerische Berechnung kollisionsbehafte- ter Gas-Partikel-Stroemungen entwickelt. Die Effizienz dieser Algorithmen wurde anhand verschiedener Testfaelle untersucht. Auf der Grundlage der dabei erzielten Ergeb- nisse wurden Vorschlaege fuer weitere Entwicklungsmoeg- lichkeiten erarbeitet. / Gas-particle-flows can be found widely in the field of energy production and process engineering. Examples for applications of such kind of flows are transport, se- paration or injection of a mixture of solid particles and a gaseous phase. The Lagrangian approach has proved to be a suitable means for the numerical simulation of disperse multiphase flows. On the other hand its application requires a large amount of computational power, especially when flows with a mo- derate or high particle loading are computed and particle- particle collisions have a significant influence on the flow. In order to use efficiently the large computational power that parallel computers provide nowadays, parallel algo- rithms for the numerical simulation of gas-particle flows including particle-particle collisions were developed in the cource of this work. The algorithms' efficiency was investigated considering different test cases. On the basis of the results suggestions for further developments were made.