Análise quantitativa de dispersão de vazamentos de substâncias inflamáveis e/ou tóxicas em ambientes com barreiras ou semi confinados.

Com o atual desenvolvimento industrial e tecnológico da sociedade, a presença de substâncias inflamáveis e/ou tóxicas aumentou significativamente em um grande número de atividades. A possível dispersão de gases perigosos em instalações de armazenamento ou em operações de transporte representam uma grande ameaça à saúde e ao meio ambiente. Portanto, a caracterização de uma nuvem inflamável e/ou tóxica é um ponto crítico na análise quantitativa de riscos. O objetivo principal desta tese foi fornecer novas perspectivas que pudessem auxiliar analistas de risco envolvidos na análise de dispersões em cenários complexos, por exemplo, cenários com barreiras ou semi-confinados. A revisão bibliográfica mostrou que, tradicionalmente, modelos empíricos e integrais são usados na análise de dispersão de substâncias tóxicas / inflamáveis, fornecendo estimativas rápidas e geralmente confiáveis ao descrever cenários simples (por exemplo, dispersão em ambientes sem obstruções sobre terreno plano). No entanto, recentemente, o uso de ferramentas de CFD para simular dispersões aumentou de forma significativa. Estas ferramentas permitem modelar cenários mais complexos, como os que ocorrem em espaços semi-confinados ou com a presença de barreiras físicas. Entre todas as ferramentas CFD disponíveis, consta na bibliografia que o software FLACS® tem bom desempenho na simulação destes cenários. Porém, como outras ferramentas similares, ainda precisa ser totalmente validado. Após a revisão bibliográfica sobre testes de campo já executados ao longo dos anos, alguns testes foram selecionados para realização de um exame preliminar de desempenho da ferramenta CFD utilizado neste estudo. Foram investigadas as possíveis fontes de incertezas em termos de capacidade de reprodutibilidade, de dependência de malha e análise de sensibilidade das variáveis de entrada e parâmetros de simulação. Os principais resultados desta fase foram moldados como princípios práticos a serem utilizados por analistas de risco ao realizar análise de dispersão com a presença de barreiras utilizando ferramentas CFD. Embora a revisão bibliográfica tenha mostrado alguns dados experimentais disponíveis na literatura, nenhuma das fontes encontradas incluem estudos detalhados sobre como realizar simulações de CFD precisas nem fornecem indicadores precisos de desempenho. Portanto, novos testes de campo foram realizados a fim de oferecer novos dados para estudos de validação mais abrangentes. Testes de campo de dispersão de nuvem de propano (com e sem a presença de barreiras obstruindo o fluxo) foram realizados no campo de treinamento da empresa Can Padró Segurança e Proteção (em Barcelona). Quatro testes foram realizados, consistindo em liberações de propano com vazões de até 0,5 kg/s, com duração de 40 segundos em uma área de descarga de 700 m2. Os testes de campo contribuíram para a reavaliação dos pontos críticos mapeados durante as primeiras fases deste estudo e forneceram dados experimentais para serem utilizados pela comunidade internacional no estudo de dispersão e validação de modelos. Simulações feitas utilizando-se a ferramenta CFD foram comparadas com os dados experimentais obtidos nos testes de campo. Em termos gerais, o simulador mostrou bom desempenho em relação às taxas de concentração da nuvem. O simulador reproduziu com sucesso a geometria complexa e seus efeitos sobre a dispersão da nuvem, mostrando claramente o efeito da barreira na distribuição das concentrações. No entanto, as simulações não foram capazes de representar toda a dinâmica da dispersão no que concerne aos efeitos da variação do vento, uma vez que as nuvens simuladas diluíram mais rapidamente do que nuvens experimentais.

Desenvolvimento da modelagem de turbulência e interação fluido-estrutura para as vibrações induzidas por vórtices de cilindro rígido.

Esta tese apresenta o desenvolvimento e aplicação de modelos de turbulência, transição laminar-turbulenta e de interações fluido-estrutura ao escoamento externo em cilindro rígido estacionário e em vibrações induzidas por vórtices. Tais desenvolvimentos foram realizados no código ReFRESCO, baseado em técnicas de dinâmica de fluidos computacional (CFD). Realizou-se um estudo quanto ao desempenho do modelo k- SST em extensa faixa de números de Reynolds, segundo o qual se identificaram as deficiências de modelagem para este escoamento. A modelagem adaptativa das escalas (SAS) e o modelo de transição por correlações locais (LCTM), ambos combinados ao SST, melhoraram a aderência aos resultados experimentais para este escoamento, em uma contribuição original deste trabalho. A aplicação de técnicas de verificação e validação possibilitou a estimação de incertezas e erros para os modelos e números de Reynolds e também de identificada como outra contribuição deste trabalho. A combinação da modelagem em SST, SAS e LCTM com movimentos impostos de realizada para números de Reynolds moderados, diferentes frequências e amplitudes de vibração, algo que poucas publicações abordam em detalhes. Com relação aos movimentos livres, este trabalho traz contribuições com a aplicação dos modelos SST e SAS ao estudo de vibrações induzidas por vórtices em dois graus de liberdade, baixa razão de massa e números de Reynolds moderados, mais altos do que normalmente observados na literatura. Por fim, a investigação da importância relativa de efeitos da turbulência aos casos de movimentos livres e impostos, com relação ao caso de cilindro estacionário, comprovou a conjetura formulada na parte inicial deste trabalho, no que tange à escolha do modelo de turbulência em determinadas aplicações. Tal escolha mostrou-se menos decisiva no caso do cilindro em movimento imposto e ainda menos nos movimentos livres, em comparação ao caso estacionário, uma vez que a resposta em movimentos do corpo filtra grande parte dos efeitos turbulentos de ordem superior. Esta observação mostra-se relevante, uma vez que pode permitir simplificações na modelagem e aplicação de ferramentas de CFD em uma classe importante de projetos de engenharia.

Otimização acústica e análise numérica do escoamento ao redor de um conjunto cilindro-placa separadora.

RESUMO Simulações de aeroacústica computacional demandam uma quantidade considerável de tempo, o que torna complicada a realização de estudos paramétricos. O presente trabalho propõe uma metodologia viável para otimização aeroacústica. Através da análise numérica utilizando dinâmica dos fluidos computacional, foi estudada a aplicação de uma placa separadora desacoplada como método de controle passivo da esteira turbulenta de um cilindro e avaliou-se a irradiação de ruído causado pela interação do escoamento com ambos os corpos, empregando ferramentas de aeroacústica computacional baseadas no método de Ffowcs-Williams e Hawkings. Algumas abordagens distintas de metodologias de otimização de projeto foram aplicadas neste problema, com o objetivo de chegar a uma configuração otimizada que permita a redução do nível sonoro ao longe. Assim, utilizando uma ferramenta de otimização multidisciplinar, pode-se avaliar a capacidade de modelos heurísticos e a grande vantagem do emprego de algoritmos baseados em método de superfície de resposta quando aplicados em um problema não linear, pois requerem a avaliação de um menor número de alternativas para se obter um ponto ótimo. Além disso, foi possível identificar e agrupar os resultados em 5 clusters baseados em seus parâmetros geométricos, nível de pressão sonora global e o valor quadrático médio do coeficiente de arrasto, confirmando a eficiência da aplicação de placas separadoras longas desacopladas posicionadas próximas ao cilindro na estabilização da esteira turbulenta, enquanto que o posicionamento de placas acima de um espaçamento crítico aumentou o nível de pressão acústica irradiado devido à formação de vórtices no espaço entre o cilindro e a placa separadora.

Coupled Biomechanical Response of the Cornea Assessed by Non-Contact Tonometry. A Simulation Study

The mechanical response of the cornea subjected to a non-contact air-jet tonometry diagnostic test represents an interplay between its geometry, the corneal material behavior and the loading. The objective is to study this interplay to better understand and interpret the results obtained with a non-contact tonometry test. A patient-specific finite element model of a healthy eye, accounting for the load free configuration, was used. The corneal tissue was modeled as an anisotropic hyperelastic material with two preferential directions. Three different sets of parameters within the human experimental range obtained from inflation tests were considered. The influence of the IOP was studied by considering four pressure levels (10–28 mmHg) whereas the influence of corneal thickness was studied by inducing a uniform variation (300–600 microns). A Computer Fluid Dynamics (CFD) air-jet simulation determined pressure loading exerted on the anterior corneal surface. The maximum apex displacement showed a linear variation with IOP for all materials examined. On the contrary, the maximum apex displacement followed a cubic relation with corneal thickness. In addition, a significant sensitivity of the apical displacement to the corneal stiffness was also obtained. Explanation to this behavior was found in the fact that the cornea experiences bending when subjected to an air-puff loading, causing the anterior surface to work in compression whereas the posterior surface works in tension. Hence, collagen fibers located at the anterior surface do not contribute to load bearing. Non-contact tonometry devices give useful information that could be misleading since the corneal deformation is the result of the interaction between the mechanical properties, IOP, and geometry. Therefore, a non-contact tonometry test is not sufficient to evaluate their individual contribution and a complete in-vivo characterization would require more than one test to independently determine the membrane and bending corneal behavior.

Aceleración de corrientes eólicas de reducida velocidad para la generación eléctrica doméstica

En los últimos años, los sistemas que utilizan como fuente recursos renovables se han posicionado como una interesante alternativa para la producción de energía. Entre las fuentes disponibles, la energía eólica viene configurándose como una de las fuentes de energía renovable con mayor crecimiento en los últimos años. En este trabajo se propone el aprovechamiento de las corrientes eólicas circulantes sobre la superficie terrestre para la producción de electricidad con el fin de abastecer buena parte de la demanda en viviendas aisladas, pequeñas instalaciones agropecuarias, equipamiento de servicio ubicado en lugares remotos, etc. Por lo general estas brisas tienen una baja densidad energética, por ello proponemos una interfaz mecánica que concentre las masas de aire, acelerando su circulación y alcanzando importantes incrementos en la velocidad de impulsión. La primera parte se centra en la elaboración de un procedimiento de caracterización a partir de la metodología científica con el cual modelar una estructura concentradora de flujo eólico válida para un aerogenerador de eje vertical. Este método trata el diseño de un elemento acelerador capaz de optimizar el aprovechamiento de estas brisas con independencia de la dirección de éstas. Su diseño viene dado por la resolución de un conjunto de objetivos fundamentales, dotando al sistema de unas prestaciones particulares en relación a su arquitectura y operatividad. Estos objetivos son los siguientes: - Operatividad ante cualquier dirección eólica adoptada - Incremento del rendimiento potencial de la turbina de eje vertical - Minimizar el desarrollo de efectos turbulentos alrededor del sistema integrado - Capacidad resolutiva ante la presencia de fuertes vientos - Estabilidad estructural - Compatibilidad ante instalaciones propias del volumen arquitectónico - Control global del rendimiento del sistema. La segunda parte aborda el modelado del prototipo y el análisis de su comportamiento mediante simulaciones en el ámbito de la fluidodinámica computacional. El resultado es un prototipo caracterizado por una arquitectura capaz de sectorizar la entrada de viento en diferentes tramos inyectando el flujo eólico estratégicamente. La incorporación de la interfaz sobre el rotor aumenta la superficie de captación eólica, facilitando su entrada a través de las diferentes aberturas y llevando a cabo su concentración según avanza por los tramo de circulación. Una vez finalizado dicho avance, el flujo es inyectado en un rango angular de nido por la elevada fuerza de sustentación capaz de generarse gracias a la incidencia del flujo eólico, aprovechando las particularidades que ofrece la rotación de este tipo de rotores. El resultado de la inyección sectorizada es el desarrollo de una circulación interior vorticial que incide permanentemente en el rango de sustentación característico del perfil aerodinámico que define la geometría de la pala de rotación. Ello provoca que se alcance un funcionamiento nominal a velocidades reducidas. En este proceso se incluyen las acciones necesarias para dar una respuesta eficiente a cualquier tipo e solicitación eólica. En presencia de velocidades de relativa importancia, la interfaz concentradora adapta su arquitectura con el fin de regular la entrada de flujo, retrasando la activación de los dispositivos reguladores propios de la turbina eólica. En presencia de vientos importantes, la interfaz dispone de los mecanismos necesarios para proceder al cierre de las aberturas procediendo a la parada del rotor. La validación de los prototipos elaborados se ha llevado a cabo mediante simulación computacional CFD. Los resultados confirman la consecución de un funcionamiento nominal a velocidades más reducidas, una mayor superficie de captación y periodo de tiempo de funcionamiento efectivo en comparación a las turbinas convencionales. Para el caso práctico modelado los resultados mejoran en más de 2.5 veces la potencia generada y multiplican por cuatro la fuerza ejercida sobre las palas de la turbina. La elaboración de un método preciso para el diseño de este tipo de estructuras concentradoras posibilita que se pueda alcanzar un diseño en función a las necesidades del usuario, o las condiciones eólicas de un emplazamiento dado. A ello hay que unir su compatibilidad de uso y montaje con sistemas de captación solar, conformando un sistema híbrido capaz de aprovechar tanto la energía solar como eólica para el abastecimiento autónomo. Esta característica incrementa el ratio de zonas geográficas donde se puede llevar a cabo su implantación. Las últimas páginas están reservadas a esbozar las líneas futuras de desarrollo y evolución, tanto en términos de e ciencia productiva como en su incorporación a nuevos volúmenes arquitectónicos y estructuras civiles en general.

Análise de balística externa de um projétil de calibre 155 mm

O objetivo desta Dissertação é propor uma geometria para um projétil de combate a incêndios. É utilizado o código de volume finito Star - CCM+R para calcular forças e momentos que atuam no projétil e é criado um código no MathematicaR para definir as caraterísticas da trajetória. O PRODASV 3 R é utilizado de modo a verificar e validar os modelos implementados. Inicialmente, o procedimento adotado consistiu em criar uma metodologia completa para uma munição existente que permitisse ser verificada e validada com a literatura disponível. Foi realizada uma extensa caraterização aerodinâmica através de CFD, a fim de se estabelecer o comportamento das forças e momentos que atuam nos projéteis e assim fornecer os dados necessários para o modelo 6 DOF. As forças e momentos dependem do ângulo de ataque. Concluiu-se que a cinta de travamento, a cinta de guiamento e o efeito de Coriolis eram desprezáveis quando comparados com os resultados obtidos. Na fase seguinte foram testadas seis geometrias diferentes, variando-se o comprimento e o nariz. Definiu-se para este projétil um a velocidade inicial de 100m/s, 10kg, concluindo-se que o comprimento total seria de 697,1mm e perfil semiesférico, com capacidade para transportar 7,5 dm3 no interior. Foram realizados testes com vento para avaliar a estabilidade em voo, e o software e processo de análise desenvolvido neste trabalho estão disponíveis, sendo um contributo para o projeto FIREND.

Metodologie di calcolo automatico di schiere palettate

La progettazione e la modellazione delle geometrie di corpi complessi come le schiere palettate delle turbomacchine, da sempre impegna il tecnico che, dapprima su carta, e successivamente in forma digitale, deve scontrarsi con le difficoltà sia analitiche di risoluzione di sistemi di equazioni differenziali, che geometriche a causa della doppia curvatura dei profili stessi. L’avvento dei calcolatori ha inevitabilmente giocato un ruolo fondamentale nella rapida evoluzione di tecniche di modellazione, calcolo e rappresentazione, per aiutare il progettista a risolvere completamente il problema, o almeno riscontrare risultati approssimativamente corretti, al fine di ridurre i tempi di realizzazione e i costi dell’impresa. Si vuole dunque cercare di descrivere le fasi che la progettazione oggi richiede, sfruttando quello che i software moderni mettono a disposizione, con l’obiettivo di mostrare uno dei molteplici percorsi che il progettista oggi può seguire per riuscire nel suo scopo.

Axisymmetric shock wave interaction with a cone: a benchmark test

Results of the benchmark test are presented of comparing numerical schemes solving shock wave of M-s = 2.38 in nitrogen and argon interacting with a 43 degrees semi-apex angle cone and corresponding experiments. The benchmark test was announced in Shock Waves Vol. 12, No. 4, in which we tried to clarify the effects of viscosity and heat conductivity on shock reflection in conical flows. This paper summarizes results of ten numerical and two experimental applications. State of the art in studies regarding the shock/cone interaction is clarified.

A parametric study of droplet deformation through a microfluidic contraction: Low viscosity Newtonian droplets

Numerical simulations are conducted to investigate how a droplet of Newtonian liquid. entrained in a higher viscosity Newtonian liquid, behaves when passing through an axisymmetric microfluidic contraction. Simulations are performed using a transient Volume of Fluid finite volume algorithm, and cover ranges of Reynolds and Weber numbers relevant to microfluidic flows. Results are presented for a droplet to surrounding fluid viscosity ratio of 0.001. In contrast to behaviour at higher viscosity ratios obtained previously by the authors, shear and interfacial tension driven instabilities often develop along the droplet Surface. leading to complex shape development, and in some instances, droplet breakup. (c) 2006 Elsevier Ltd. All rights reserved.

Numerical simulation of magnetite segregation in a dense medium cyclone

Numerical simulations of turbulent driven flow in a dense medium cyclone with magnetite medium have been conducted using Fluent. The predicted air core shape and diameter were found to be close to the experimental results measured by gamma ray tomography. It is possible that the Large eddy simulation (LES) turbulence model with Mixture multi-phase model can be used to predict the air/slurry interface accurately although the LES may need a finer grid. Multi-phase simulations (air/water/medium) are showing appropriate medium segregation effects but are over-predicting the level of segregation compared to that measured by gamma-ray tomography in particular with over prediction of medium concentrations near the wall. Further, investigated the accurate prediction of axial segregation of magnetite using the LES turbulence model together with the multi-phase mixture model and viscosity corrections according to the feed particle loading factor. Addition of lift forces and viscosity correction improved the predictions especially near the wall. Predicted density profiles are very close to gamma ray tomography data showing a clear density drop near the wall. The effect of size distribution of the magnetite has been fully studied. It is interesting to note that the ultra-fine magnetite sizes (i.e. 2 and 7 mu m) are distributed uniformly throughout the cyclone. As the size of magnetite increases, more segregation of magnetite occurs close to the wall. The cut-density (d(50)) of the magnetite segregation is 32 gm, which is expected with superfine magnetite feed size distribution. At higher feed densities the agreement between the [Dungilson, 1999; Wood, J.C., 1990. A performance model for coal-washing dense medium cyclones, Ph.D. Thesis, JKMRC, University of Queensland] correlations and the CFD are reasonably good, but the overflow density is lower than the model predictions. It is believed that the excessive underflow volumetric flow rates are responsible for under prediction of the overflow density. (c) 2006 Elsevier Ltd. All rights reserved.

Study of flow behaviour in a three-product cyclone using computational fluid dynamics

Simplicity in design and minimal floor space requirements render the hydrocyclone the preferred classifier in mineral processing plants. Empirical models have been developed for design and process optimisation but due to the complexity of the flow behaviour in the hydrocyclone these do not provide information on the internal separation mechanisms. To study the interaction of design variables, the flow behaviour needs to be considered, especially when modelling the new three-product cyclone. Computational fluid dynamics (CFD) was used to model the three-product cyclone, in particular the influence of the dual vortex finder arrangement on flow behaviour. From experimental work performed on the UG2 platinum ore, significant differences in the classification performance of the three-product cyclone were noticed with variations in the inner vortex finder length. Because of this simulations were performed for a range of inner vortex finder lengths. Simulations were also conducted on a conventional hydrocyclone of the same size to enable a direct comparison of the flow behaviour between the two cyclone designs. Significantly, high velocities were observed for the three-product cyclone with an inner vortex finder extended deep into the conical section of the cyclone. CFD studies revealed that in the three-product cyclone, a cylindrical shaped air-core is observed similar to conventional hydrocyclones. A constant diameter air-core was observed throughout the inner vortex finder length, while no air-core was present in the annulus. (c) 2006 Elsevier Ltd. All rights reserved.

Numerical analysis of gas and micro-particle interactions in a hand-held shock-tube device

A unique hand-held gene gun is employed for ballistically delivering biomolecules to key cells in the skin and mucosa in the treatment of the major diseases. One of these types of devices, called the Contoured Shock Tube (CST), delivers powdered micro-particles to the skin with a narrow and highly controllable velocity distribution and a nominally uniform spatial distribution. In this paper, we apply a numerical approach to gain new insights in to the behavior of the CST prototype device. The drag correlations proposed by Henderson (1976), Igra and Takayama (1993) and Kurian and Das (1997) were applied to predict the micro-particle transport in a numerically simulated gas flow. Simulated pressure histories agree well with the corresponding static and Pitot pressure measurements, validating the CFD approach. The calculated velocity distributions show a good agreement, with the best prediction from Igra & Takayama correlation (maximum discrepancy of 5%). Key features of the gas dynamics and gas-particle interaction are discussed. Statistic analyses show a tight free-jet particle velocity distribution is achieved (570 +/- 14.7 m/s) for polystyrene particles (39 +/- 1 mu m), representative of a drug payload.

Expansion tube operating conditions for studying nonequilibrium radiation relevant to Titan Aerocapture

The predictions of nonequilibrium radiation in the shock layer for a Titan aerocapture aeroshell vary significantly amongst Computational Fluid Dynamics (CFD) analyses and are limited by the physical models of the nonequilibrium flow processes. Of particular interest are nonequilibrium processes associated with the CN molecule which is a strong radiator. It is necessary to have experimental data for these radiating shock layers which will allow for validation of the CFD models. This paper describes the development of a test flow condition for subscale aeroshell models in a superorbital expansion tunnel. We discuss the need for a Titan gas condition that closely simulates the atmospheric composition and present experimental data of the free stream test flow conditions. Furthermore, we present finite-rate CFD calculations of the facility to estimate the remaining free stream conditions, which cannot be directly measured during experiments.

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.

Sump simulations

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 in the containment sump. Two dispersed phases were conditions to determine the influence of entrained air from a jet on the transport of fibre agglomerates through the sump. The strainer model of A. Grahn was implemented to observe the impact that the accumulation of the fibres have on the pressure drop across the strainers. The geometry considered is similar to the containment sump configurations found in Nuclear Power Plants.