Le but de cette thèse est de décrire en détail la dynamique linéaire et non linéaire d'une couche limite attachée sur une plaque plane à bas nombre de Reynolds. La dynamique linéaire, pilotée par les interactions entre les vecteurs propres non-orthogonaux, est étudiée à travers deux méthodes différentes d'instabilité globale : une analyse globale aux vecteurs propres et une optimisation directe-adjointe. Dans ces analyses globales, aucune structure spatiale n'est imposée à priori pour la perturbation, les effets convectifs dus au fort non parallélisme de l'écoulement sont pris en compte. Pour le cas de la couche limite décollée, le déclenchement des instationnarités a été clarifié : i) pour une forte amplification des perturbations de nature convective et bidimensionnelle ; ii) pour des effets de non normalité longitudinale engendrant le phénomène du flapping ; iii) pour une forte sensibilité vis-à-vis d'un forçage harmonique ; iv) pour un monde tridimensionnel globalement instable. Pour une couche limite attachée, le but a été d'identifier les perturbations localisées caractérisées par des fréquences multiples dans les directions longitudinales et transversales de l'écoulement, induisant une amplification de l'énergie des perturbations. Pour évaluer les effets les effets de la non linéarité dans les mécanismes d'instabilité identifiés par les analyses de stabilité globale, des simulations numériques directes ont été réalisées pour les écoulements de couche limite attachées ou décollées, bidimensionnelles et tridimensionnelles. La dynamique des perturbations permettant une plus rapide vers la turbulence a été étudiée. Différents scénarios de transition ont été observés, les différents mécanismes de transition ont été analysés.

Les situations dans lesquelles une onde de choc interagit avec une couche limite sont nombreuses dans les industries aéronautiques et spatiales. Sous certaines conditions (nombre de Mach élevé, grand angle de choc…), ces interactions entrainent un décollement de la couche limite. Des études antérieures ont montré que la zone de recirculation et le choc réfléchi sont tous deux soumis à un mouvement d'oscillation longitudinale à basse fréquence connu sous le nom d’instabilité de l’interaction onde de choc / couche limite (IOCCL). Ce phénomène appelé soumet les structures à des chargement oscillants à basse fréquence qui peuvent endommager les structures.L’objectif du travail de thèse est de réaliser des simulations instationnaires de l’IOCCL afin de contribuer à une meilleure compréhension de l’instabilité de l’IOCCL et des mécanismes physiques sous-jacents.Pour effectuer cette étude, une approche numérique originale est utilisée. Un schéma « One step » volume fini qui couple l’espace et le temps, repose sur une discrétisation des flux convectifs par le schéma OSMP développé jusqu’à l’ordre 7 en temps et en espace. Les flux visqueux sont discrétisés en utilisant un schéma aux différences finies centré standard. Une contrainte de préservation de la monotonie (MP) est utilisée pour la capture de choc. La validation de cette approche démontre sa capacité à calculer les écoulements turbulents et la grande efficacité de la procédure MP pour capturer les ondes de choc sans dégrader la solution pour un surcoût négligeable. Il est également montré que l’ordre le plus élevé du schéma OSMP testé représente le meilleur compromis précision / temps de calcul. De plus un ordre de discrétisation des flux visqueux supérieur à 2 semble avoir une influence négligeable sur la solution pour les nombres de Reynolds relativement élevés considérés.En simulant un cas d’IOCCL 3D avec une couche limite incidente laminaire, l’influence des structures turbulentes de la couche limite sur l’instabilité de l’IOCCL est supprimée. Dans ce cas, l’unique cause d’IOCCL suspectée est liée à la dynamique de la zone de recirculation. Les résultats montrent que seul le choc de rattachement oscille aux fréquences caractéristiques de la respiration basse fréquence du bulbe de recirculation. Le point de séparation ainsi que le choc réfléchi ont une position fixe. Cela montre que dans cette configuration, l’instabilité de l’IOCCL n’a pas été reproduite.Afin de reproduire l’instabilité de l’IOCCL, la simulation de l’interaction entre une onde de choc et une couche limite turbulente est réalisée. Une méthode de turbulence synthétique (Synthetic Eddy Method - SEM) est développée et utilisée à l’entrée du domaine de calcul pour initier une couche limite turbulente à moindre coût. L’analyse des résultats est effectuée en utilisant notamment la méthode snapshot-POD (Proper Orthogonal Decomposition). Pour cette simulation, l’instabilité de l’IOCCL a été reproduite. Les résultats suggèrent que la dynamique du bulbe de recirculation est dominée par une respiration à moyenne fréquence. Ces cycles successifs de remplissage / vidange de la zone séparée sont irréguliers dans le temps avec une taille maximale du bulbe de recirculation variant d’un cycle à l’autre. Ce comportement du bulbe de recirculation traduit une modulation basse fréquence des amplitudes des oscillations des points de séparation et de recollement et donc une respiration basse fréquence de la zone séparée. Ces résultats suggèrent que l’instabilité de l’IOCCL est liée à cette dynamique basse fréquence du bulbe de recirculation, les oscillations du pied du choc réfléchi étant en phase avec le point de séparation. / Situations where an incident shock wave impinges upon a boundary layer are common in the aeronautical and spatial industries. Under certain circumstances (High Mach number, large shock angle...), the interaction between an incident shock wave and a boundary layer may create an unsteady separation bubble. This bubble, as well as the subsequent reflected shock wave, are known to oscillate in a low-frequency streamwise motion. This phenomenon, called the unsteadiness of the shock wave boundary layer interaction (SWBLI), subjects structures to oscillating loads that can lead to damages for the solid structure integrity.The aim of the present work is the unsteady numerical simulation of (SWBLI) in order to contribute to a better understanding of the SWBLI unsteadiness and the physical mechanism causing these low frequency oscillations of the interaction zone.To perform this study, an original numerical approach is used. The one step Finite Volume approach relies on the discretization of the convective fluxes of the Navier Stokes equations using the OSMP scheme developed up to the 7-th order both in space and time, the viscous fluxes being discretized using a standard centered Finite-Difference scheme. A Monotonicity-Preserving (MP) constraint is employed as a shock capturing procedure. The validation of this approach demonstrates the correct accuracy of the OSMP scheme to predict turbulent features and the great efficiency of the MP procedure to capture discontinuities without spoiling the solution and with an almost negligible additional cost. It is also shown that the use of the highest order tested of the OSMP scheme is relevant in term of simulation time and accuracy compromise. Moreover, an order of accuracy higher than 2-nd order for approximating the diffusive fluxes seems to have a negligible influence on the solution for such relatively high Reynolds numbers.By simulating the 3D unsteady interaction between a laminar boundary layer and an incident shock wave, we suppress the suspected influence of the large turbulent structures of the boundary layer on the SWBLI unsteadiness, the only remaining suspected cause of unsteadiness being the dynamics of the separation bubble. Results show that only the reattachment point oscillates at low frequencies characteristic of the breathing of the separation bubble. The separation point of the recirculation bubble and the foot of the reflected shock wave have a fixed location along the flat plate with respect to time. It shows that, in this configuration, the SWBLI unsteadiness is not observed.In order to reproduce and analyse the SWBLI unsteadiness, the simulation of a shock wave turbulent boundary layer interaction (SWTBLI) is performed. A Synthetic Eddy Method (SEM), adapted to compressible flows, has been developed and used at the inlet of the simulation domain for initiating the turbulent boundary layer without prohibitive additional computational costs. Analyses of the results are performed using, among others, the snapshot Proper Orthogonal Decomposition (POD) technique. For this simulation, the SWBLI unsteadiness has been observed. Results suggest that the dominant flapping mode of the recirculation bubble occurs at medium frequency. These cycles of successive enlargement and shrinkage of the separated zone are shown to be irregular in time, the maximum size of the recirculation bubble being submitted to discrepancies between successive cycles. This behaviour of the separation bubble is responsible for a low frequency temporal modulation of the amplitude of the separation and reattachment point motions and thus for the low frequency breathing of the separation bubble. These results tend to suggest that the SWBLI unsteadiness is related to this low frequency dynamics of the recirculation bubble; the oscillations of the reflected shocks foot being in phase with the motion of the separation point.

Nous cherchons à modéliser numériquement un jet flottant air-hélium dans une cavité avec deux ouvertures à partir de simulations aux grandes échelles (LES) et de simulations numériques directes (DNS). La configuration considérée est basée sur une étude expérimentale menée au CEA de Saclay reproduisant une fuite d'hydrogène en environnement confiné. La dimension de la cavité a été choisie pour permettre une transition laminaire-turbulent intervenant environ à la mi-hauteur de la cavité. Cette étude porte principalement sur trois points majeurs : l'influence des conditions aux limites sur le développement du jet et son interaction avec l'environnement extérieur, la validité du modèle numérique qui est analysée en comparant la distribution de vitesse obtenue numériquement aux mesures expérimentales (PIV) et, enfin, la compréhension de la distribution air-hélium et le phénomène de stratification qui s'établit à l'intérieur de la cavité. Nous observons dans un premier temps que des conditions limites de pression constante appliquées directement au ras des évents conduisent à une sous-estimation du débit volumique d'air entrant dans la cavité et donc à une surestimation de la masse de l'hélium à l'intérieur de la cavité, ce qui n'est pas acceptable dans un contexte d'évaluation du risque hydrogène. En revanche, la prise en compte, dans le domaine de calcul, d'une région extérieure à la cavité prédit correctement le flux d'air entrant. Les résultats numériques sont alors en bon accord avec les données PIV. Il a été montré que les prédictions de la DNS, par rapport à la LES, concordent mieux avec les mesures de vitesse par PIV. Le champ de concentration prédit numériquement présente une couche homogène en haut de la cavité, dont la concentration est en accord avec le modèle théorique de Linden et al. 1990. Cependant, sa position et son épaisseur ne correspondent pas au modèle. Ceci est principalement dû aux interactions directes entre le jet flottant et, d'une part, avec les limites solides de la cavité et d'autre part, avec l'environnement extérieur. L'analyse statistique concernant la production de la flottabilité de l'énergie cinétique turbulente (TKE) a permis d'identifier les limites du jet flottant. / We present numerical results from large eddy simulations (LES) and coarse direct numerical simulations (DNS) of an air-helium buoyant jet rising in a two vented cavity. The geometrical configuration mimics the helium release experimental set-up studied at CEA Saclay in the framework of security assessment of hydrogen-based systems with an indoor usage. The dimension of the enclosure was chosen to ensure a laminar-turbulent transition occurring at about the middle height of the cavity. This study focuses mainly on three key points : the influence of the boundary conditions on the jet development and its interaction with the exterior environment, the validity of the numerical model which is analyzed by comparing the numerical velocity distribution versus the measured particle image velocimetry (PIV) ones, and finally understanding the distribution of air-helium and the stratification phenomenon that takes place inside the cavity. We observe at first that applying constant pressure outlet boundary conditions directly at the vent surfaces underestimates the volumetric flow rate of air entering the enclosure and thus overestimate the helium mass inside the cavity. On the contrary, modelling an exterior region in the computational domain better predicts the air flow-rate entrance and numerical results matches better with the experimental PIV data. It has been figured out that the coarse DNS predictions match better with the velocities PIV measurements, compared to the LES. Numerical prediction of the helium field depicts a homogeneous layer formed at the top of the cavity, with a concentration in good agreement with the theoretical model of Linden et al. 1990. However, the position and the thickness of the layer do not correspond to the theory. This is mainly due to the direct interactions between the buoyant jet and both the solid boundaries of the cavity and the exterior environment. Statistical analysis regarding the buoyancy production of the turbulent kinetic energy (TKE) served to identify the limits of the buoyant jet.

La marée interne contribue au maintien de la circulation méridienne de renversement. Il existe, à l'heure actuelle, une controverse sur la nature exacte des mécanismes pilotant cette circulation. Une meilleure quantification des apports énergétiques associés aux ondes internes permettrait d'apporter quelques clés de compréhension de ce mécanisme complexe. Dans cette thèse, différents régimes d'ondes internes topographiques inspirés par des congurations océaniques sont étudiés an d'évaluer quantitativement les transferts énergétiques associés. L'utilisation complémentaire des outils numériques et expérimentaux permet de détailler la dynamique de ces régimes de manière exhaustive. La première partie de la thèse porte sur la génération d'ondes internes linéaires par l'oscillation d'un mont Gaussien dans un fluide linéairement stratifié. L'approche choisie s'appuie sur un jeu d'expériences de laboratoire pour lesquelles la pente relative du rayon d'onde interne par rapport à la pente maximale du mont varie. Nous montrons qu'un maximum énergétique est atteint dans le régime critique pour lequel les pentes du rayon et du mont sont similaires. Dans la suite de la thèse, la dynamique d'ondes internes de forte amplitude se propageant dans des régions de fort gradient de densité, comme la pycnocline océanique, est étudiée. Nous utilisons dans un premier temps le modèle numérique Symphonie-NH pour décrire leur génération et leur dynamique, sur une configuration académique bidimensionnelle. Tout d'abord, la génération primaire d'ondes internes interfaciales est étudiée. On s'intéresse en particulier à des régimes fortement non-linéaires pour lesquelles des ondes solitaires sont observées. Elle sont induites par l'interaction directe entre la marée barotrope et la topographie et est observée dans des régimes de pycnocline de forte intensité dans l'océan, comme en mer de Sulu. La structure des ondes internes solitaires est étudiée avec des modèles analytiques simples comme l'équation KdV. En modifiant la forme du mont, un contrôle topographique important de la génération primaire d'ondes internes solitaires est observé. Un paramètre adimensionné est proposé pour décrire ce contrôle. Ensuite, la génération secondaire d'ondes internes solitaires induites par l'interaction de rayons d'ondes internes émis sur une topographie avec une pycnocline d'intensité modérée, comme dans le Golfe de Gascogne, est étudiée. Des simulations numériques directes sont effectuées pour décrire la dynamique des ondes internes solitaires, et leur atténuation par radiation d'énergie dans la couche inférieure. L'évolution de la structure des modes normaux en fonction de l'intensité de la pycnocline, et le rôle joué par la forme du mont sont détaillés. Des expériences sont menées pour étudier la génération primaire et la génération secondaire d'ondes internes solitaires dans le grand canal du CNRM-GAME. Une configuration expérimentale utilisant un mont sinusoïdal oscillant dans la couche inférieure, stratifée ou non, d'un fluide bicouche est adoptée. Cette configuration, inspirée des simulations numériques précédentes, permet d'explorer une gamme plus large de régimes d'ondes interfaciales. Des mesures de déplacement interfacial avec des sondes à ultrasons d'une part, et avec des mesures optiques d'autre part, permettent de discuter la dynamique, et la structure tridimensionnelle de ces ondes. La structure des ondes internes solitaires dans le cas de la génération primaire apparaît plus stable que pour la génération secondaire. Dans ce deuxième cas, des structures transverses régulières sont mesurées.

Cette thèse porte sur la caractérisation numérique et la modélisation stochastique de l'accélération du fluide pour l'écoulement en canal à grand nombre de Reynolds. La motivation concerne l'observation et l'analyse des effets de l'intermittence liés aux interactions à longue portée à travers le canal. Dans la première partie, l'accélération est étudiée par simulation numérique directe pour trois différents nombres de Reynolds (180, 590 et 1000). La lognormalité de la norme de l'accélération est observée quelle que soit la distance à la paroi. Un profil universel de la norme de l'accélération est également recherché par analyse dimensionnelle. La seconde partie présente une modélisation stochastique de l'accélération basée sur la décomposition norme/orientation. Le modèle stochastique pour la norme s'appuie sur un processus de fragmentation afin de représenter les interactions à longue portée à travers le canal. Pour l'orientation, l'évolution vers l'isotropie lorsque la distance à la paroi augmente (observée par la DNS) est reproduite grâce à un modèle de marche aléatoire sur une sphère. Ces modèles ont été appliqués à l'approche LES-SSAM (Stochastic Subgrid Acceleration Model) introduite par Sabel'nikov, Chtab et Gorokhovski. Nos calculs montrent que les estimations de la vitesse moyenne, du spectre d'énergie, des contraintes de l'écoulement et de la non-gaussianité des statistiques de l'accélération peuvent être améliorées de façon significative par rapport à la LES classique. L'intérêt de l'approche LES-SSAM, donnant un accès vers la structure intermittente de sous-maille, est illustré dans la dernière partie, par l'étude du transport de particules inertielles ponctuelles par l'écoulement de canal. Cette étude commence par l'analyse par DNS de l'influence des structures de paroi sur la dynamique des particules

Les contacteurs centrifuges, basés sur les écoulements de Taylor-Couette, sont bien adaptés pour la mise en œuvre de réactions chimiques ou biochimiques, y compris en milieu polyphasique. Ils possèdent particulièrement plusieurs propriétés favorables à la mise en œuvre des opérations d'extraction liquide-liquide. Un dispositif expérimental a été conçu avec cette idée en tête. Il est constitué de deux cylindres concentriques avec le cylindre intérieur entraîné en rotation et l'externe fixe. L’écoulement de Taylor-Couette se produit dans l’espace annulaire entre eux. Il présente la particularité d’évoluer vers la turbulence par apparition successive d’instabilités. La dispersion axiale ainsi que le mélange, sont extrêmement sensibles à ces structures d’écoulement, ce qui rend difficile la modélisation du couplage entre l’hydrodynamique et le transfert de matière. Ce point particulier a été étudié expérimentalement et numériquement. L’écoulement et le mélange ont été caractérisés par des mesures simultanées de PIV (Vélocimétrie par Imagerie de Particules) et PLIF (Fluorescence Induite par Laser). Les champs de concentration PLIF ont permis d’identifier les différents mécanismes de transport intra et inter-vortex. Pour les régimes ondulatoires (WVF et MWVF), le mélange intra-vortex est contrôlé par l’advection chaotique, directement lié aux caractéristiques du champ de vitesse, qui confère aux vortex une capacité plus importante à convecter et à étirer les filets de fluide. En revanche, l’apparition des vagues brisent les frontières qui séparent les vortex ce qui favorise le transport inter-vortex. La combinaison de ces deux mécanismes contrôle principalement la dispersion axiale. Nous avons également mis en évidence le comportement non monotone des propriétés de mélange en fonction de l’histoire de l’écoulement. Notamment l’état d’onde (la longueur d’onde axiale et l’amplitude de la vague). Nous avons calculé le coefficient de dispersion axiale Dx à l’aide des mesures de distribution de temps de séjour (DTS) et de suivi Lagrangien de particules (DNS). Les deux résultats numériques et expérimentaux ont confirmé l’effet significatif des structures de l’écoulement et de l’histoire sur la dispersion axiale. / Taylor-Couette flows between two concentric cylinders have great potential applications in chemical engineering. They are particularly convenient for two-phase small scale devices enabling solvent extraction operations. An experimental device was designed with this idea in mind. It consists of two concentric cylinders with the inner one rotating and the outer one fixed. Taylor-Couette flows take place in the annular gap between them, and are known to evolve towards turbulence through a sequence of successive instabilities. Macroscopic quantities, such as axial dispersion and mixing index, are extremely sensitive to these flow structures, which may lead to flawed modelling of the coupling between hydrodynamics and mass transfer. This particular point has been studied both experimentally and numerically. The flow and mixing have been characterized by means of flow visualization and simultaneous PIV (Particle Imaging Velocimetry) and PLIF (Planar Laser Induced Fluorescence) measurements. PLIF visualizations showed clear evidences of different transport mechanisms including « intravortex mixing » and « inter-vortex mixing ». Under WVF and MWVF regimes, intra-vortex mixing is controlled by chaotic advection, due to the 3D nature of the flow, while inter-vortex transport occurs due to the presence of waves between neighbouring vortices. The combination of these two mechanisms results in enhanced axial dispersion. We showed that hysteresis may occur between consecutive regimes depending on flow history and this may have a significant effect on mixing for a given Reynolds number. The axial dispersion coefficient Dx evolution along the successive flow states was investigated thanks to dye Residence Time Distribution measurements (RTD) and particle tracking (DNS). Both experimental and numerical results have confirmed the significant effect of the flow structure and history on axial dispersion. Our study confirmed that the commonly used 1-parameter chemical engineering models (e.g. the « well-mixed stirred tanks in serie » model) are not valid for Taylor-Couette reactors modelling : two parameters are at least required for an efficient description of mixing in Taylor-Couette flows.

Résumé : La mesure des profils transitoires et de la vitesse de solidification sont deux données importantes pour le contrôle de procédés industriels impliquant un changement de phase. Dans le cas de l’électrolyse de l’aluminium, ce processus de solidification assure la protection du système et influe sur la performance énergétique du procédé de fabrication. Malheureusement, ces données se révèlent, dans la plupart des cas, difficilement accessibles. Ce travail de thèse porte sur le développement de nouveaux outils permettant l’étude et la caractérisation de la solidification de matériaux à changement de phase et à haute température. L’objectif est de développer un système de mesure du front de solidification de matériaux à changement de phase non destructif et ne perturbant pas le milieu de mesure, tout en assurant une précision et une réponse suffisamment rapide pour exploiter de nouvelles stratégies de contrôle dans les cuves d’électrolyse. Ce travail couple une étude expérimentale fondamentale de la solidification de la cryolithe avec une modélisation numérique de phénomène de changement de phase solide-liquide dans des conditions proches du fonctionnement de cuves d’électrolyse. // Abstract : Measurement of transient solidification fronts and of solidification rate are two important data for controlling industrial processes involving a solid-liquid phase change. In the case of aluminium electrolysis, this solidification process protects the system and affects the energy performance of the manufacturing process. Unfortunately, these data are not easy to obtain in most cases. This thesis focuses on the development of new tools for the study and on the solidification characterization of phase change materials at high temperature. The goal is to develop a nondestructive solidification front measurement system for phase change materials without disturbing the measurement medium, while ensuring accuracy and a fast enough response time to exploit new control strategies in electrolysis cells. This work couples a fundamental experimental study of the cryolite solidification with numerical modeling of solid-liquid phase change phenomenon under conditions close to those during normal operation of electrolytic cells.

Cette thèse porte sur la caractérisation numérique et la modélisation stochastique de l’accélération du fluide pour l’écoulement en canal à grand nombre de Reynolds. La motivation concerne l’observation et l’analyse des effets de l’intermittence liés aux interactions à longue portée à travers le canal. Dans la première partie, l’accélération est étudiée par simulation numérique directe pour trois différents nombres de Reynolds (180, 590 et 1000). La lognormalité de la norme de l’accélération est observée quelle que soit la distance à la paroi. Un profil universel de la norme de l’accélération est également recherché par analyse dimensionnelle. La seconde partie présente une modélisation stochastique de l’accélération basée sur la décomposition norme/orientation. Le modèle stochastique pour la norme s’appuie sur un processus de fragmentation afin de représenter les interactions à longue portée à travers le canal. Pour l’orientation, l’évolution vers l’isotropie lorsque la distance à la paroi augmente (observée par la DNS) est reproduite grâce à un modèle de marche aléatoire sur une sphère. Ces modèles ont été appliqués à l’approche LES-SSAM (Stochastic Subgrid Acceleration Model) introduite par Sabel’nikov, Chtab et Gorokhovski. Nos calculs montrent que les estimations de la vitesse moyenne, du spectre d’énergie, des contraintes de l’écoulement et de la non-gaussianité des statistiques de l’accélération peuvent être améliorées de façon significative par rapport à la LES classique. L’intérêt de l’approche LES-SSAM, donnant un accès vers la structure intermittente de sous-maille, est illustré dans la dernière partie, par l’étude du transport de particules inertielles ponctuelles par l’écoulement de canal. Cette étude commence par l’analyse par DNS de l’influence des structures de paroi sur la dynamique des particules / The main objective of this thesis is to observe numerically and to analyze the effects of intermittency in a high Reynolds number turbulent channel flow. To this end, the thesis is focused on characterization and stochastic modelling of the fluid acceleration in such a flow, with emphasis on long-range interactions across the channel. In the first part, the acceleration is studied using DNS for three Reynolds numbers (180, 590 et1000). It is observed that the norm of acceleration is log-normal whatever the wall distance is. The universal form of scaling law for the acceleration is proposed by dimensional analysis. In the second part, the acceleration is simulated stochastically, assuming the norm/orientation decomposition. The stochastic model for the norm is based on the fragmentation process in order to represent the long-range interactions across the channel. The orientation is simulated by random walk on a sphere in order to reproduce the relaxation towards isotropy with increasing the wall distance. This was observed preliminary in our DNS. These models were applied in the framework of LES-SSAM approach (Stochastic Subgrid AccelerationModel), which was introduced by Sabel’nikov, Chtab and Gorokhovski and assessed in the case of the box turbulence. Our computations showed that the mean velocity, the energy spectra, the viscous and turbulent stresses, as well as the non-gaussianity of acceleration statistics can be considerably improved in comparison with standard LES. The advantage of the LES-SSAM approach, which accounts for intermittency on subgrid scales, is demonstrated in the last part of this thesis. Here the transport of inertial point wise particles was studied by DNS and by LES-SSAM. The influence of wall structures on the particle’s dynamics is analyzed.