Fenômenos multifásicos são frequentemente observados na natureza, tais como nas gotas de chuva ou neve no ar, nos vulcões e tempestades de areia, e em diversas outras situações. Na solução desses problemas que envolvem escoamentos gás-sólidos e granulares são frequentemente utilizadas duas abordagens: a contínua (formulação Euleriana-Euleriana) e a discreta (formulação Euleriana-Lagrangiana). Na abordagem discreta pode-se utilizar dois modelos para descrever o contato entre as partículas: o modelo de esfera rígida e o modelo de esfera suave. Neste trabalho é realizado um estudo detalhado dos modelos de contato, com foco na modelagem de esfera suave, que é baseada em um sistema dinâmico mola-massa-amortecedor. Por meio desse estudo, com a finalidade de aprimorar o modelo de contato não-linear, são propostas duas relações para o mecanismo de contato de partículas. Essas relações são fundamentadas em um modelo dinâmico, com não-linearidades nas partes conservativas e dissipativas, não apresentando descontinuidades entre as acelerações do início e do fim do contato. A metodologia de desenvolvimento da presente pesquisa está dividida em três partes: pesquisa bibliográfica dos modelos de contato; estudo analítico e numérico desses modelos e testes de problemas com a realização de experimentos numéricos, utilizando o código computacional MFIX (Multiphase Flow with Interphase eXchange). As novas aproximações propostas neste trabalho são analisadas e aplicadas em três diferentes problemas: de dinâmica, escoamento gás-sólido e escoamento granular. Os resultados obtidos utilizando as relações são comparados com dados disponíveis na literatura, mostrando-se adequados para os casos investigados neste trabalho. / Multiphase flow are frequently observed in nature, such as rain drops in air or snowfalls, volcanoes and sandstorms, and several other situations. For solving these problems which involve gas-solid and granular flows are often used two models: the continuous model (Eulerian-Eulerian formulation) and the discrete model (Eulerian-Lagrangian formulation). There are two main contact models used in DEM, the hard-sphere model and the soft-sphere model. In this work is carried out a detailed study of contact models, focusing on soft-sphere model, based on a dynamic system modeled as nonlinear mass-spring-damper. In order to improve the nonlinear contact model, in this study it is proposed two new approximate relations for determining the damping coefficient and duration of contact for a specific nonlinear soft-sphere contact model where the contact force is continuous at the start and end of the contact. The methodology of the development of this work is divided into three parts: literature research of the contact models; analytical and numerical study of these models and test problems with numerical experiments, using the open source code MFIX (Multiphase Flow with Interphase eXchange). The proposed relations are analyzed and applied in three different problems: dynamic problem, gas-solid flow and granular flow. All results are compared with literature data showing good agreement for these cases studied in the present work.

Filtered two-fluid models used to perform large scale simulations of gas-solid fluidized flows of industrial risers require closures for filtered parameters such as filtered and residual stresses, and interphase interaction forces mainly effective drag. Closure models for those filtered parameters may be derived by averaging over results of highly resolved simulations with microscopic two-fluid modeling. This work is a contribution in that context. Recent models for filtered parameters have been written as functions of filter size, filtered solid volume fraction, and filtered slip velocity. A recent study showed that macro-scale variables like domain average solid volume fraction and gas Reynolds number also significantly affect the filtered parameters. In the current work, in addition to these filtered and macro-scale variables, the effects of two new variables over the filtered parameters are investigated: filtered solid kinetic energy and sub-grid gas turbulence. It is shown that the filtered solid kinetic energy should be accounted for in the concerning correlations, thereby improving accuracy. Regarding gas turbulence, literature shows it has no significant effects on the motion of high Stokes particles. Extending on literature, this work investigates the sub-grid gas turbulence effects on meso-scale structures formed of high Stokes particles. Results showed that sub-grid gas turbulence has no significant effects on the meso-scale structures and corresponding filtered parameters. The open source code MFIX was used for all simulations. Ranges of dilute concentration of solid and gas Reynolds number typical of riser flow regimes were considered. A modified two-fluid model with microscopic formulation was used. The sub-grid gas turbulence was generated by means of a forcing function procedure which was implemented in the physical space, over the gravitational term of the momentum conservation equation of the gas phase. First, numerical simulations of the gas phase alone were performed, accounting for literature available data, in order to set a turbulent gas field and calibrate the turbulence intensity. Then, the forcing function was introduced in to the two-fluid model and various gas-solid flows were simulated. While the current results show the necessity of accounting for additional variables in the filtered parameter correlation, they also make it clear the necessity of further developments that are required in the search for better accuracy. / Modelos filtrados de dois-fluidos usados em simulações de grandes escalas de escoamentos fluidizados de gás-sólido de risers industriais exigem fechamentos para parâmetros filtrados tais como as tensões filtradas e residuais, e forças interativas interfases, principalmente arrasto efetivo. Modelos de fechamento para estes parâmetros filtrados podem ser gerados a partir de procedimentos de media aplicados sobre resultados de simulações altamente resolvidas com modelagem microscópica de dois-fluidos. Este trabalho é uma contribuição neste contexto. Modelos de fechamento recentes para parâmetros filtrados tem sido formulados em função de tamanho de filtro, fração volumétrica de sólido filtrada, e velocidade de deslizamento filtrada. Estudo recente mostrou que variáveis de macro-escalas como fração volumétrica de sólido e número de Reynolds de gás médios no domínio também afetam significativamente os parâmetros filtrados. No presente trabalho, além dessas variáveis filtradas e de macro-escala, os efeitos de duas novas variáveis sobre os parâmetros filtrados são investigados: energia cinética filtrada do sólido e turbulência submalha do gás. Em relação à energia cinética filtrada do sólido, mostra-se que a sua consideração refina as correlações em questão, contribuindo assim para melhor acuracidade. Com relação à turbulência do gás, a literatura mostra que não tem efeitos significativos no movimento de particulados de elevados números de Stokes. Acrescentando à literatura, este trabalho investiga os efeitos da turbulência sub-malha do gás sobre estruturas de meso-escala formados de particulados de elevados números de Stokes. Os resultados mostraram que a turbulência sub-malha do gás não tem efeitos significativos sobre estruturas de meso-escalas e parâmetros filtrados correspondentes. O código aberto MFIX foi usado para todas as simulações. Faixas de concentração diluída de sólido e número de Reynolds típicos de escoamentos em risers foram considerados. Um modelo modificado de dois fluidos com formulação microscópica foi utilizado. A turbulência sub-malha do gás foi gerada por meio de um procedimento de \'forcing function\' que foi implementado no espaço físico, sobre o termo fonte gravitacional da equação de momentum da fase gás. Primeiramente, simulações numéricas da fase gás foram realizadas separadamente, levando-se em conta dados disponíveis na literatura, a fim de gerar um campo de gás turbulento e calibrar a intensidade de turbulência. Posteriormente, a \'forcing function\' foi introduzida no modelo de dois-fluidos e vários escoamentos de gás-sólido foram simulados. Enquanto os resultados obtidos mostram a necessidade de consideração de variáveis adicionais para correlação de parâmetros filtrados, também deixam claro a necessidade de desenvolvimentos mais aprofundados na busca de melhor acuracidade.

Dada a complexidade hidrodinâmica de escoamentos gás-sólidos em leitos fluidizados muitas pesquisas têm sido realizadas abordando questões relativas a natureza altamente instável desse processo de fluidização. O conhecimento dessa hidrodinâmica permite estabelecer os parâmetros corretos para reação e transporte de massa em um reator de leito fluidizado, permitindo tomada de decisões à respeito do desempenho do reator. Modelos Eulerianos do contínuo, onde as fases coexistentes são tratadas como meios contínuos interpenetrantes, representam a aproximação mais prática para simular a hidrodinâmica de escoamentos gás-sólido em reatores de leito fluidizado. Aplicando-se procedimentos de médias de Euler obtêm-se os modelos hidrodinâmicos A e B desenvolvidos por GIDASPOW [27] do IIT/ANL, representando as equações conservativas e constitutivas de modelos de duas fases separadas para escoamentos gás-sólido. Neste trabalho, aplica-se um modelo Euleriano de duas fases separadas, desenvolvendo simulações numéricas para diferentes dimensões de colunas ascendentes de leito fluidizado circulante. Para resolução das equações discretas no domínio bidimensional, aplicam-se malhas refinadas. Analisam-se os efeitos nos resultados de simulação causados por essas malhas. Apresenta-se uma discussão geral sobre o modelo aplicado em torno dos resultados obtidos. / Due to the hydrodynamic complexity of gas-solid flows in fluidized beds, several researches have been conducted addressing issues concerning the highly unstable nature of this fluidization process. The knownledge of its hydrodynamics allows establishing the correct parameters for the reaction and mass transport in a fluidized bed reactor, improving its performance. Continuous Eulerian models, where the phases are treated as a continuum, represent the most practical approach to simulate the hydrodynamics of gas-solid flows in fluidized bed reactors. By applying Euler average procedures, it is possible to obtain the hydrodynamic models A and B, developed by GIDASPOW [27] at IIT/ANL. These models represent the conservative and constitutive equations for the models of two separate phases for gas-solid flows. In this work, an Eulerian two-fluids model is applied, developing numerical simulations for different risers of circulating fluidized bed. For the resolution of the discrete equations in the two-dimensional domain, refined meshes are used. An analysis of their effects on the results of the simulation is performed. A general discussion about the results obtained by the model applied is also presented.

Reatores de leito fluidizado circulante são caracterizados por escoamentos altamente heterogêneos e instáveis. Devido a isso, operam em regimes pseudo-estacionários ou estatisticamente permanentes, onde os parâmetros do escoamento oscilam em torno de valores médios bem definidos. Neste trabalho estuda-se o comportamento de uma solução numérica para um escoamento típico de reatores de leito fluidizado circulante na região estatisticamente permanente. A simulação numérica foi realizada com o modelo de dois fluidos. A partir dos dados de simulação são estabelecidas médias e auto-correlações para as várias variáveis de interesse, considerando-se diferentes tempos de simulação ou diferentes tempos de observação. A partir destes resultados discute-se a questão do tempo mínimo de simulação necessário para que os resultados temporais médios sejam representativos do processo. Embora os resultados sugiram que 10 segundos de fluidização sejam suficientes para a geração de resultados médios temporais representativos, uma demonstração mais rigorosa permanece por ser executada. / Circulating fluidized bed reactors are characterized by highly heterogenic and instable flow. Due to this flow, the reactors operate under either pseudo-stationary or statistically permanent patterns, where the flowing parameters vary among well defined average values. This paper aims to study the behavior of a numeric solution for a typical flow of circulating fluidized bed reactors in the statistically permanent region. The numeric simulation was performed with the two-fluid model. From the simulation data both average and self-correlations are established for several variables of interest, considering either different simulation time or different observation time. From these results, the minimum simulation time necessary is discussed, so that the average temporal results can be considered representative in the process. Although the results suggest 10 seconds of fluidization are enough to generate the average temporal results, a more rigorous demonstration is supposed to be executed.

In gas-solid fluidized beds, the generation of electrostatic charges due to continuous contacts between fluidizing particles, and the particles and the fluidization vessel wall, is unavoidable. Industrial operations, such as the production of polyethylene, are susceptible to significant operational challenges caused by electrostatics including reactor wall fouling, a problem known as “sheeting”. The formation of particle sheets can require shutdown periods for clean-up which results in significant economic losses. To gain a better understanding of the underlying mechanisms of electrostatic charging in gas-solid fluidized beds, in an attempt to eliminate or minimize this problem, a pilot-scale pressurized gas-solid fluidization system was designed and built, housing an online electrostatic charge measurement technique consisting of two Faraday cups. The system permits the study of the degree of particle wall fouling at pressures and temperatures up to 2600 kPa and 100°C, respectively, and gas velocities up to 1 m/s (covering a range including turbulent flow regime). The system also allowed, for the first time, the measurement of the fluidizing particles’ mass, net charge and size distribution in various regions of the bed, especially those related to the wall coating under the industrially relevant operating conditions of high pressures and gas velocities. Experimental trials were carried out using polyethylene resin received from commercial reactors to investigate the influence of pressure and gas velocity on the bed hydrodynamics and in turn, the degree of bed electrification. Mechanisms for particle charging, migration and adherence to the column wall were proposed. The size distribution of the gas bubbles shifted towards smaller bubbles as the operating pressure was raised. Thus, higher pressures lead to greater mixing within the bulk of the bed and resulted in a higher degree of particle wall fouling. Moreover, the extent of wall fouling increased linearly with the increase in gas velocity and as the bed transitioned to turbulent regime, due to the increase in particle-wall contacts. Bipolar charging was observed especially within the wall coating with smaller particles being negatively charged. Overall, particle-wall contacts generated negatively charged particles resulting in a net negative charge in the bed, whereas particle-particle contacts generated positively and negatively charged particles resulting in no net charge when entrainment was negligible. The formation of the wall layer and its extent was influenced by the gravitational and drag forces balancing the image force and Coulomb forces (created by the net charge of the bed and the metallic column wall as the attraction between oppositely charged particles).

This thesis builds around the idea of a biofuel production process that is comprised of biomass production, biomass gasification, gas cleaning, and fuel production. In this work, we specifically looked into H2S removal as a part of cleaning the producer gas and flocculation of microalgae which is involved in the harvesting of microalgae after biomass production. One of the impurities to remove from the producer gas is hydrogen sulfide which can be removed by using a packed bed of zinc oxide. Despite the regular use, it was only recently shown that during reaction with H2S, nano-size particles of ZnO exhibit void formation and outward growth. In this work, a micro-scale model was introduced to describe the void formation and outward growth. On the macro-scale, the simulations captured pore clogging of pellets due to the outward growth. The pore clogging prevents the full conversion of pellets and consequently leads to shorter breakthrough times of beds. The second problem investigated here deals with the flocculation of microalgae. Microalgae is produced in relatively low concentrations in the incubator liquid medium and during the harvesting, the concentration is increased to an acceptable level. The harvesting process includes a flocculation followed by a filtration or centrifuge unit. During flocculation, microalgae are stimulated to aggregate and form clusters. The experiments showed that the mean size of clusters formed during flocculation increases with time to a maximum and then starts decreasing, resulting in an overshoot in the mean size profile. The size of clusters influence the efficiency of the afterward filtration or centrifuge, thus it is of interest to carefully track the size evolution of clusters, making the studying of overshoot a crucial research topic. In this work, the possible mechanisms behind this overshoot were investigated. / <p>QC 20170330</p>

Foram desenvolvidos estudos de modelagem e simulação numérica de escoamentos bifásicos gás-sólido na coluna ascendente de leitos fluidizados circulantes utilizando um modelo Euleriano de duas fases separadas. O sistema de equações diferenciais parciais conservativas governantes foi obtido através de um procedimento tradicional. Ambas as fases foram assumidas como meio contínuo. Aplicou-se o procedimento de médias estatísticas de Euler, enfatizando a obtenção dos modelos hidrodinâmicos A e B desenvolvidos no IIT/ANL. Realizou-se análise comparativa de correlações para transferência de quantidade de movimento na interface. Discutiu-se a formulação de condições de contorno apropriadas. As equações diferenciais parciais médias foram discretizadas em volumes de controle Eulerianos. As equações de continuidade foram resolvidas implicitamente. As equações de quantidade de movimento foram resolvidas através de um procedimento explícito-implícito. Foram desenvolvidas simulações numéricas para uma coluna ascendente típica de leitos fluidizados circulantes. Desenvolveu-se análise paramétrica da influência de vários aspectos físicos e matemáticos sobre o escoamento. Avaliou-se resultados de simulação através de metodologia de identificação e caracterização de estruturas coerentes. Estudou-se o efeito da função de arrasto na interface sobre os processos dinâmicos que caracterizam estas estruturas coerentes. Foram realizados estudos numéricos de turbulência a partir de resultados de simulação direta. Várias conclusões e recomendações para futuros trabalhos foram propostas com base nas análises realizadas. Foram apresentadas algumas considerações gerais relativas a aspectos críticos na modelagem e simulação com modelo das duas fases separadas. / Studies were carried out on modeling and numerical simulation of gas-solid two-phase flows in the riser of circulating fluidized beds using an Eulerian two-fluids model. The system of conservative partial differential governing equations was derived through a traditional procedure. Both phases were assumed as a continuum. The Euler averaging procedure was applied emphasizing the derivation of the so called hydrodynamic models A and B developed at IIT/ANL. A comparative analysis was performed among correlations for momentum transfer at the interface. The formulation of suitable boundary conditions was discussed. The average partial differential conservative equations were discretized on Eulerian control volumes. The continuity equations were solved implicitly. The momentum equations were solved through an explicit-implicit procedure. Numerical simulation was performed for a typical circulating fluidized bed riser. A parametric analysis was carried out regarding the influence on the flow of various physical and mathematical aspects. Results of simulation were evaluated through a methodology of identification and characterization of coherent structures. The effect of the interface drag function on dynamic features of those coherent structures was addressed. Numerical studies on turbulence were performed from results of direct simulation. Several conclusions and recommendations for future work were put forward on the basis of the performed analyses. Some general considerations were presented regarding critical features of modeling and simulation through Eulerian two-fluids models.

A presente pesquisa tem como objetivo realizar um estudo teórico e desenvolver simulações computacionais envolvendo a dinâmica de sistemas gás-sólido. O foco principal do trabalho é a modelagem do particulado através da análise das forças de contato entre partículas de materiais granulares utilizando modelos contínuos baseados na mecânica dos solos e na teoria cinética dos escoamentos granulares (sistemas grandes com muitas partículas, formulação Euleriana - Volumes Finitos) e modelos discretos baseados nas características físicas dos materiais (sistemas intermediários e número limitado de partículas, formulação Lagrangeana - Método dos Elementos Discretos). Investigam-se os modelos existentes na literatura com intuito de melhorar os modelos contínuos e discretos baseados na interação entre as partículas que caracterizam a dinâmica do particulado em sistemas gás-sólido. Propõe-se uma nova abordagem para a determinação do coeficiente de rigidez da mola baseada em uma equivalência entre os modelos lineares e não-lineares. Utiliza-se o código fonte MFIX para realizar simulações computacionais da dinâmica de sistemas gás-sólido, analisando o processo de fluidização, mistura e segregação de partículas, influência das correlações de arrasto, e análise das forças de contato entre as partículas através do novo método para a determinação do coeficiente de rigidez da mola . Os resultados obtidos são comparados com dados numéricos e experimentais da literatura. / The purpose of the present study is to perform a theoretical study and develop numerical simulations involving dynamic in gas-solid systems. The focus of the work is the modeling of particulate matter using continuous models based on soil mechanics and the kinetic theory of granular flows (large systems with many particles, Eulerian formulation - Finite Volume) and discrete models based on physical characteristics of the particles (intermediate systems and limited number of particles, Lagrangian formulation - Discrete Element Method). It is proposed a new approach to determine the normal spring stiffness coefficient of the linear model through the numerical solution for the overlap between particles in non-linear models. The linear spring stiffness is determined using an equivalence between the linear and the non-linear models. It is used the MFIX computational code to perform numerical simulations of the dynamics of gas-solid systems. It is analyzed the processes of fluidization, mixing and particle segregation and the influence of drag correlations. The proposed approach for normal spring stiffness coefficient is applied in the numerical simulations of two problems: single freely falling particle and bubbling fluidized bed. The results were compared with numerical and experimental data from literature.

Gas-solid fluidized bed reactors are widely applied in chemical and energy industries, and their design and scale-up are virtually empirical, extremely expensive and time consuming. This scenario has motivated the development of alternative theoretical tools, and two-fluid modeling, where gas and particulate are both treated as interpenetrating continuum phases, has appeared as a most promising approach. Owing to the large domains to be resolved in real-scale fluidized bed reactors, only filtered modeling approaches are feasible, and closure models become necessary to recover sub-grid effects that are filtered by the very coarse numerical grids that are imposed owing to computational limitations. Those closure models, which in hydrodynamic formulations account for filtered interphase momentum exchanges and filtered and residual stresses in the phases, can be derived from results of highly resolved simulations (HRS) performed over small size domains under refined numerical grids. One widely practiced approach consists of applying two-fluid modeling under micro-scale defined closures, generally known as microscopic two-fluid modeling. This approach includes microscopic closures for solid phase stresses derived from the kinetic theory of granular flows (KTGF), which accounts for kinetic-collisional effects only, and is adequate to dilute flows. Otherwise, the conventional KTGF does not account for interparticle friction effects, and its application to dense flow conditions is quite questionable. In this work a literature available modified version of KTGF is applied which also accounts for interparticle friction, and highly resolved simulations are performed for dense flow conditions in order to evaluate the effects of friction over relevant filtered parameters (namely effective drag coefficient, filtered and residual stresses). Ranges of domain average solid volume fractions and gas Reynolds numbers are considered (macro-scale conditions) embracing dense gas-solid fluidized flows from suspensions up to pneumatic transport. The MFIX open source code is used in all the simulations, which are performed over 2D periodical domains for a unique monodisperse particulate. The HRS results (i.e. meso-scale flow fields) are filtered over regions compatible with grid sizes in large scale simulations, and the relevant filtered parameters of concern are derived and classified by ranges of other filtered parameters taken as independent variables (filtered solid volume fraction, filtered slip velocity, and filtered kinetic energy of solid velocity fluctuations, which are referred to as markers). Results show that the relevant filtered parameters of concern are well correlated to all of those filtered markers, and also to all of the imposed macro-scale conditions. Otherwise, interparticle friction showed no significant effects over any filtered parameter. It is recognized that this issue clearly requires further investigation notably regarding the suitability of the markers that were assumed for classifying the filtered results. The current work is intended as a contribution for future developments of more accurate closure models for large scale simulations of gas-solid fluidized flows. / Reatores de leito fluidizado de escoamento gás-sólido são largamente utilizados nas indústrias química e de energia, e o seu projeto e escalonamento são virtualmente empíricos, extremamente caros e demorados. Este cenário tem motivado o desenvolvimento de ferramentas teóricas alternativas, e a modelagem de dois fluidos, onde gás e particulado são ambos tratados com fases contínuas interpenetrantes, tem surgido como uma aproximação muito promissora. Devido aos grandes domínios a serem resolvidos em reatores de leito fluidizado de escala real, apenas aproximações de modelagem filtradas são viáveis, e modelos de fechamento tornam-se necessários para recuperar efeitos sub-malha que são filtrados pelas malhas numéricas grosseiras que são impostas devido as limitações computacionais. Estes modelos de fechamento, que em formulações hidrodinâmicas respondem principalmente por trocas de momentum filtradas entre fases e tensões filtradas e residuais nas fases, podem ser obtidos de resultados de simulações altamente resolvidas (SAR) realizadas em domínios de dimensões reduzidas sob malhas numéricas refinadas. Uma aproximação largamente praticada consiste na aplicação de modelagem de dois fluidos sob fechamentos definidos na micro-escala, genericamente conhecida como modelagem microscópica de dois fluidos. Esta aproximação inclui fechamentos microscópicos para tensões da fase sólida obtidos da teoria cinética dos escoamentos granulares (TCEG), que considera apenas efeitos cinéticos-colisionais, e é adequada para escoamentos diluídos. Por outro lado, a TCEG convencional não leva em conta efeitos de fricção interpartículas, e sua aplicação para condições densas de escoamento é bastante questionável. Neste trabalho aplica-se uma versão modificada da TCEG disponível na literatura que também leva em conta fricção interpartículas, e simulações altamente resolvidas são realizadas para condições de escoamentos densos visando avaliar os efeitos da fricção sobre os parâmetros filtrados relevantes (coeficiente de arrasto efetivo, tensões filtradas e residuais). Considera-se faixas de frações volumétricas de sólido e números de Reynolds do gás médios no domínio (condições de macro-escala) abrangendo escoamentos gás-sólido fluidizados densos desde suspensões até transporte pneumático. O código aberto MFIX é utilizado em todas as simulações, que foram executadas sobre domínios periódicos 2D para um único particulado monodisperso. Os resultados das SAR (i.e., campos de escoamento de meso-escala) foram filtrados sobre regiões compatíveis com tamanhos de malha praticados em simulações de grandes escalas, e os parâmetros filtrados relevantes de interesse são calculados e classificados por faixas de outros parâmetros filtrados tomados como variáveis independentes (fração volumétrica de sólido filtrada, velocidade de deslizamento filtrada, e energia cinética das flutuações de velocidade da fase sólida filtrada, que são referidos como marcadores). Os resultados mostram que os parâmetros filtrados relevantes de interesse são bem correlacionados com todos os marcadores, e também com todas as condições de macro-escala impostas. Por outro lado, a fricção interpartículas não mostrou efeitos significativos sobre qualquer parâmetro filtrado. Reconhece-se que este aspecto claramente requer investigações adicionais, notadamente com respeito à adequação dos marcadores que foram considerados para classificação dos resultados filtrados. O trabalho corrente é posto como uma contribuição para o desenvolvimento futuro de modelos de fechamento mais acurados para simulações de grandes escalas de escoamentos gás-sólido fluidizados.

Les processus élémentaires hétérogènes à l’interface gaz-solide présentent un intérêt fondamental dans de nombreux domaines tels que la catalyse hétérogène, la chimie atmosphérique et des milieux interstellaires, la rentrée atmosphérique de véhicules spatiaux ou encore la description des interactions plama-paroi. Cette thèse a pour objet l’étude de la dynamique des processus de collision non réactive de l’azote N2 sur une surface de tungstène W(100) et des processus de recombinaison moléculaire de l’hydrogène H2 sur des surfaces de tungstène W(100) et W(110). Leur dynamique quasi classique est simulée au moyen de surfaces d’énergie potentielle préalablement construites à partir de calculs de théorie de la fonctionnelle de la densité. Un potentiel multi-adsorbats est notamment développé pour tenir compte du taux de couverture de surface afin d’étudier la compétition entre la recombinaison directe, de type Eley-Rideal et la recombinaison par « atomes chauds » après diffusion hyperthermique d’un atome sur la surface. / Heterogeneous elementary processes at the gas-solid interface are ofgreat interest in many domains such as heterogeneous catalysis, atmospheric and interstellar media chemistry, spacecraft atmospheric re-entry and plasma-wall interactions description. This thesis focus on the dynamics of nitrogen, N2, non reactive scattering on a tungsten W(100) surface and hydrogen, H2, recombination processes on tungsten surfaces W(100) and W(110). The quasiclassical dynamics of these processes is simulated using potential energy surfaces based on density functional theory calculations. In particular, a multi-adsorbate potential is developed to include surface coverage in the dynamics simulation in order to scrutinize the interplay between both direct abstraction, the so-called Eley-Rideal recombination,and the Hot-Atom recombination process after hyperthermal diffusion on the surface