La demande croissante pour alléger les moteurs d’avions et diminuer les émissions polluantes de la propulsion aéronautique réclame à rendre plus compact le système de compression des moteurs, qui représente environ 40%-50% de la masse totale. Or, à taux de compression global égal, la réduction du nombre d’étage exige d’un point de vue aérodynamique une augmentation de la charge des aubes de compresseur par étage. La charge d’aube est aujourd’hui limitée car elle induit différents mécanismes de pertes tridimensionnelles très pénalisant. L’un des plus importants est le décollement de coin qui se forme à la jonction entre l’extrados de l’aube et le moyeu ou le carter. Bien que des travaux existent sur les mécanismes et paramètres intervenant dans le décollement de coin, il est encore difficile de proposer une méthode de contrôle efficace. Cela est principalement dû à deux raisons : (i) le manque de compréhension fine des mécanismes physiques, (ii) l’utilisation pour la conception de modèles de turbulence classiques de type RANS (Reynolds-averaged Navier-Stokes) qui ne sont pas capables de prédire précisément le décollement de coin, car ils ne peuvent pas décrire correctement les mécanismes de transport turbulent. Des simulations de type RANS et LES (large-eddy simulation = simulation des grandes échelles) sont présentées dans cette thèse sur une configuration de grille d’aubes de compresseur, et comparées avec les données expérimentales obtenues au LMFA (issues de travaux séparés). L’approche RANS surestime globalement le décollement de coin. Une amélioration significative est obtenue par la méthode LES, en particulier pour le coefficient de pression statique sur l’aube et les pertes de pression totale. Ces résultats montrent que la zone de décollement de coin, qui est la source principale des pertes, génère des tourbillons de grande échelle associés à de forts niveaux d’énergie. Les histogrammes bimodaux de la vitesse tangentielle qui ont été observés expérimentalement semblent confirmés par les résultats LES. En ce qui concerne les amplitudes des fluctuations de vitesse tangentielle, les résultats expérimentaux et ceux de la LES mettent en évidence deux pics sur certains profils perpendiculaires aux parois. Enfin, grâce à l’approche LES, les bilans de l’énergie cinétique turbulente sont calculés et analysés. Ils décrivent l’équilibre entre les termes de production, de dissipation et de transport. Une des perspectives de cette analyse est d’aider à améliorer la modélisation de la turbulence en approche RANS. / The increasing demand to reduce the mass of aircraft jet engines and emissions of aircraft propulsion requires to make the compression system of engines more compact, since this component accounts for about 40%-50% of the total mass. However, at a given overall pressure ratio, decreasing the number of stages will raise the compressor blade loading per stage. The blade loading is extremely restricted by different three-dimensional flow loss mechanisms. One of them is the corner separation that forms between the blade suction side and the hub or shroud. Although some works previously investigated the mechanisms and the parameters of corner separation, it is still difficult to propose an effective control method of the corner separation. That is mainly due to two reasons: (i) the lack of knowledge of the physical mechanisms, (ii) the nowadays classical RANS (Reynolds-averaged Navier-Stokes) turbulence models are not capable to accurately predict the corner separation, since they cannot correctly describe the turbulent transport mechanisms. RANS (Reynolds-averaged Navier-Stokes) and LES (large-eddy simulation) simulations are here presented on a compressor cascade configuration, in comparison with experimental data obtained at LMFA (from separate works). The RANS approach globally over-estimates the corner separation, whereas a significant improvement is achieved with the LES, especially for the blade surface static pressure coefficient and the total pressure losses. The corner separation region, which is the main source of the total pressure losses, is shown to generate large-scale energy-containing eddies. The bimodal histograms of the streamwise velocity that were observed experimentally seem to be confirmed by the LES results. Concerning the streamwise velocity fluctuations (RMS), both the experiment and the LES show some profiles with two peaks. Finally, thanks to the LES approach, the turbulent kinetic energy budget, which represents the balance between the production, dissipation and transport terms, are computed and analyzed. This may help the improvement of RANS turbulence modeling.

Afin d’élaborer des stratégies de contrôle des écoulements en temps réel, il est nécessaire d’avoir recours à des modèles d’ordre réduit (ROMs), car la résolution des équations complètes est trop coûteuse en temps de calcul (des jours, des semaines) et en espace mémoire. Dans cette thèse, les modèles réduits ont été construits avec la méthode POD (Proper Orthogonal Decomposition). Une méthode de projection basée sur la minimisation des résidus, initiée par les travaux de Leblond et al. [134] a été proposée. Dans certaines configurations, la précision des résultats est significativement augmentée, par rapport à une projection de Galerkin classique. Dans un second temps, un algorithme d’optimisation non-linéaire, à direction de descente basée sur la méthode des équations adjointes, a été couplé avec des modèles réduits utilisant des bases POD. Deux méthodes de construction de base POD ont été employées : soit avec un paramètre (un nombre de Reynolds,. . . ), soit avec plusieurs paramètres (plusieurs nombres de Reynolds, . . . ). Les ROMs obtenus ont été utilisés pour contrôler la dispersion d’un polluant dans une cavité ventilée puis pour contrôler le champ de température dans une cavité entraînée différentiellement chauffée. Le contrôle est réalisé en temps quasi-réel et les résultats obtenus sont plutôt satisfaisants. Néanmoins, ces méthodes sont encore trop coûteuses en espace mémoire pour être aujourd’hui embarquées dans les boîtiers de contrôle utilisés dans le bâtiment. Une autre stratégie de contrôle, s’appuyant sur les contrôleurs actuels, a ainsi été développée. Celle-ci permet d’obtenir la température (ainsi que la vitesse) dans la zone d’occupation du bâtiment, en utilisant une décomposition des champs par POD et un algorithme d’optimisation de Levenberg-Marquardt. Elle a été validée sur une cavité différentiellement chauffée, puis appliquée sur une cavité ventilée 3D, proche d’un cas réel. / In order to control flows in real-time, it is necessary to resort to reduced-order models (ROMs) because the classical method of simulations is too expensive in CPU time (several days, weeks) and memory storage. In this thesis, the ROMs have been built with the POD (Proper Orthogonal Decomposition) technique. First, a projection method based on the minimization of the equations residuals and established starting from the works of Leblond et al. [134] have been developed. In some cases, the results accuracy is significantly increased. Secondly, a direct descent optimization algorithm based on adjoint-equations has been coupled with POD/ROMs. Two construction methods of POD bases has been employed: either with simulations for only one parameter (one Reynolds number, . . . ), or with simulations for several parameters (several Reynolds numbers,. . . ). The obtained ROMs have been applied in order to control the pollutant dispersion and then to control the temperature field in a lid-driven cavity heated by the left. The control is realized in quasi-real time and the results are rather satisfying. Nevertheless, these methods are still too expensive in memory storage to be embedded in the current controllers. Thus, another control strategy has been proposed, using POD and an optimization algorithm (Levenberg-Marquardt). This one enables to obtain the temperature (and the velocity) in the occupation zone of the building and has been validated on the lid-driven cavity heated by the left and applied on a 3D-ventilated cavity, similar to a real case.

Au cours des dernières décennies, des progrès remarquables ont été réalisés au niveau de la simulation d’écoulements sanguins dans des modèles anatomiques réalistes construits à partir de données d'imagerie médicale 3D en vue de simulation hémodynamique et physiologique 3D à grande échelle. Alors que les modèles anatomiques précis sont d'une importance primordiale pour simuler le flux sanguin, des conditions aux limites réalistes sont également importantes surtout lorsqu’il s’agit de calculer des champs de vitesse et de pression. La première cible de cette thèse était d'étudier l'analyse de convergence des inconnus pour différents types de conditions aux limites permettant un cadre flexible par rapport au type de données d'entrée (vitesse, pression, débit, ...). Afin de faire face au grand coût informatique associé, nécessitant un calcul haute performance, nous nous sommes intéressés à comparer les performances de deux préconditionneurs par blocs; le preconditionneur LSC (Least-Squared Commutator et le preconditionneur PCD (Pressure Convection Diffusion). Dans le cadre de cette thèse, nous avons implémenté ce dernier dans la bibliothèque Feel++. Dans le but de traiter l'interaction fluide-structure, nous nous sommes focalisés sur l'approximation de la force exercée par le fluide sur la structure, un champ essentiel intervenant dans la condition de continuité pour assurer le couplage du modèle de fluide avec le modèle de structure. Enfin, afin de valider nos choix numériques, deux cas tests ont été réalisés et une comparaison avec les données expérimentales et numériques a été établie et validée (le benchmark FDA et le benchmark Phantom). / Towards a large scale 3D computational model of physiological hemodynamics, remarkable progress has been made in simulating blood flow in realistic anatomical models constructed from three-dimensional medical imaging data in the past few decades. When accurate anatomic models are of primary importance in simulating blood flow, realistic boundary conditions are equally important in computing velocity and pressure fields. Thus, the first target of this thesis was to investigate the convergence analysis of the unknown fields for various types of boundary conditions allowing for a flexible framework with respect to the type of input data (velocity, pressure, flow rate, ...). In order to deal with the associated large computational cost, requiring high performance computing, we were interested in comparing the performance of two block preconditioners; the least-squared commutator preconditioner and the pressure convection diffusion preconditioner. We implemented the latter, in the context of this thesis, in the Feel++ library. With the purpose of handling the fluid-structure interaction, we focused of the approximation of the force exerted by the fluid on the structure, a field that is essential while setting the continuity condition to ensure the coupling of the fluid model with the structure model. Finally, in order to assess our numerical choices, two benchmarks (the FDA benchmark and the Phantom benchmark) were carried out, and a comparison with respect to experimental and numerical data was established and validated.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal.

Flow of granular materials is usually initiated when the shear stress is large enough and exceeds certain critical value. This can result in the presence of the dead-zones in which the flow itself does not take place. Motions of such materials are frequently described by Bingham model. Flows of granular fluids are frequently connected with the presence of free surface. In the thesis Bingham model is incorporated into a more general framework of Bingham-Korteweg fluids, which is a suitable way how to transfer free- boundary problems into the problems on fixed domains. A part of the thesis concerns mathematical analysis of interesting relevant problems for incompressible fluids. 1

Ce travail a pour objectif de générer des modèles vasculaires et de simuler des écoulements sanguins réalistes à l'intérieur de ces modèles. La première étape consiste à segmenter/reconstruire le volume 3D du réseau vasculaire. Une fois de tels volumes vasculaires segmentés et maillés, il est alors possible de simuler des écoulements sanguins à l'intérieur de ceux-ci. Pour la segmentation, nous utilisons une approche variationnelle. Nous proposons un premier modèle qui inclut un a priori de tubularité dans les modèles de débruitage ROF et TV-L1. Néanmoins, bien que ces modèles permettent de réhausser les vaisseaux, ils ne permettent pas de les segmenter. C'est pourquoi nous proposons un deuxième modèle amélioré qui inclut à la fois un a priori de tubularité et de direction dans le modèle de segmentation de Chan-Vese. Les résultats sont présentés sur des images synthétiques 2D, ainsi que sur des images rétiniennes. En ce qui concerne la simulation, nous nous intéressons d'abord au réseau veineux cérébral, encore peu étudié. Les équations de la dynamique des fluides qui régissent les écoulements sanguins dans notre géométrie sont alors les équations de Navier-Stokes. Pour résoudre ces équations, la méthode classique des caractéristiques est comparée avec un schéma d'ordre plus élevé. Ces deux schémas sont validés sur des solutions analytiques avant d'être appliqués aux cas réalistes du réseau veineux cérébral premièrement, puis du polygone artériel de Willis. / The aim of this work is to generate vascular models and simulate blood flows inside these models. A first step consists of segmenting/reconstructing the 3D volume of the vascular network. Once such volumes are segmented and meshed, it is then possible to simulate blood flows. For segmentation purposes, we use a variational approach. We first propose a model that embeds a vesselness prior in the denoising models ROF and TV-L1. Although these models can enhance vessels, they are not designed for segmentation. Then, we propose a second, improved model that includes both vesselness and direction priors in the Chan-Vese segmentation model. The results are presented on 2D synthetic images, as well as retinal images. In the second part, devoted to simulation, we first focus on the cerebral venous network, that has not been intensively studied. The equations governing blood flows inside our geometry are the Navier-Stokes equations. For their resolution, the classical method of characteristics is compared with a high-order scheme. Both schemes are validated on analytical solutions before their application on the realistic cases of the cerebral venous network, and the arterial polygon of Willis.

Turbulence is ubiquitous in the flows of fluids and plasmas. This thesis is devoted to studies of the statistical properties of turbulence in the three-dimensional (3D) Hall magnetohydrodynamic (Hall-MHD) equations, the two-dimensional (2D) MHD equations, the one-dimensional (1D) hyperviscous Burgers equation, and the 3D Navier-Stokes equations. Chapter 1 contains a brief introduction to statistically homogeneous and isotropic turbulence. This is followed by an over-view of the equations we study in the subsequent chapters, the motivation for the studies and a summary of problems we investigate in chapters 2-6. In Chapter 2 we present our study of Hall-MHD turbulence [1]. We show that a shell-model version of the 3D Hall-MHD equations provides a natural theoretical model for investigating the multiscaling behaviors of velocity and magnetic structure functions. We carry out extensive numerical studies of this shell model, obtain the scaling exponents for its structure functions, in both the low-k and high-k power-law ranges of 3D Hall-MHD, and find that the extended-self-similarity procedure is helpful in extracting the multiscaling nature of structure functions in the high-k regime, which otherwise appears to display simple scaling. Our results shed light on intriguing solar-wind measurements. In Chapter 3 we present our study of the inverse-cascade regime in two-dimensional magnetohydrodynamic turbulence [2]. We present a detailed direct numerical simulation (DNS) of statistically steady, homogeneous, isotropic, two-dimensional magnetohydrodynamic (2D MHD) turbulence. Our study concentrates on the inverse cascade of the magnetic vector potential. We examine the dependence of the statistical properties of such turbulence on dissipation and friction coefficients. We extend earlier work significantly by calculating fluid and magnetic spectra, probability distribution functions (PDFs) of the velocity, magnetic, vorticity, current, stream-function, and magnetic-vector-potential fields and their increments. We quantify the deviations of these PDFs from Gaussian ones by computing their flatnesses and hyperflatnesses. We also present PDFs of the Okubo-Weiss parameter, which distinguishes between vortical and extensional flow regions, and its magnetic analog. We show that the hyperflatnesses of PDFs of the increments of the stream-function and the magnetic vector potential exhibit significant scale dependence and we examine the implication of this for the multiscaling of structure functions. We compare our results with those of earlier studies. In Chapter 4 we compare the statistical properties of 2D MHD turbulence for two different energy injection scales. We present systematic DNSs of statistically steady 2D MHD turbulence. Our two DNSs are distinguished by kinj, the wave number at which we inject energy into the system. In our first DNS (run R1), kinj = 2 and, in the second (run R2) kinj = 250. We show that various statistical properties of the turbulent states in the runs R1 and R2 are strikingly different The nature of energy spectrum, probability distribution functions, and topological structures are compared for the two runs R1 and R2 are found to be strikingly different. In Chapter 5 we study the hyperviscous Burgers equation for very high α, order of hyperviscosity [3]. We show, by using direct numerical simulations and theory, how, by increasing α in equations of hydrodynamics, there is a transition from a dissipative to a conservative system. This remarkable result, already conjectured for the asymptotic case α →∞ [U. Frisch et al., Phys. Rev. Lett. 101, 144501 (2008)], is now shown to be true for any large, but finite, value of α greater than a crossover value α crossover. We thus provide a self-consistent picture of how dissipative systems, under certain conditions, start behaving like conservative systems, and hence elucidate the subtle connection between equilibrium statistical mechanics and out-of-equilibrium turbulent flows. In Chapter 6 we show how to use asymptotic-extrapolation and Richardson extrapolation methods to extract the exponents ξ p that characterize the dependence of the order-p moments of the velocity gradients on the Reynolds number Re. To use these extrapolation methods we must have high-precision data for such moments. We obtain these high-precision data by carrying out the most extensive, quadruple precision, pseudospectral DNSs of the Navier-Stokes equation.

Cette thèse a pour but l'étude de la mise en œuvre de commandes par asservissement visuel pour le contrôle actif d'un écoulement de Poiseuille. D'un point de vue général, le contrôle d'écoulements vise à modifier ou à maintenir l'état de l'écoulement, malgré une éventuelle perturbation extérieure. Une des situations d'intérêt concerne par exemple la transition vers la turbulence où l'écoulement peut devenir turbulent avec la croissance de sa densité d'énergie cinétique. La réduction de la traînée est également une application potentielle dans des problèmes d'ingénierie. Un des buts applicatifs de cette thèse cherchera ainsi à minimiser à la fois la densité d'énergie cinétique et la traînée. Des modèles numériques peuvent être utilisés pour générer un modèle d'état des équations aux dérivées partielles d'un écoulement de Poiseuille. Le modèle d'état considéré dans cette thèse s'appuie sur une représentation spectrale afin de transformer les équations aux dérivées partielles originelles en un système d'équations différentielles ordinaires. Le vecteur d'état rassemble dans notre cas la vitesse et la vorticité. Les signaux de commande dépendent eux de conditions aux limites de type Dirichlet non homogènes qui correspondent à des actions de soufflage/aspiration. Le nombre de degrés de liberté commandé du problème correspond à la dimension du signal de commande. La densité d'énergie cinétique et la traînée sont modélisées en fonction du vecteur d'état et du signal de commande. Dans cette thèse nous avons plus particulièrement considéré un asservissement visuel partitionné. Celui-ci est appliqué au modèle d'état de l'écoulement avec deux degrés de liberté afin de minimiser simultanément la densité d'énergie cinétique et la traînée. La traînée, contrairement à l'énergie cinétique, diminue de façon monotone en fonction du temps. Une augmentation du nombre de degrés de liberté permet d'améliorer la décroissance de la densité d'énergie cinétique. Lorsque le nombre de degré de liberté correspond à la dimension du vecteur d'état, et en s'appuyant sur une commande par asservissement visuel, nous montrons que la densité d'énergie cinétique décroit de façon monotone au cours du temps. Le modèle d'état de l'écoulement de Poiseuille vit dans un espace de très grande dimension. Par conséquent, il est nécessaire d'un point de vue pratique de réduire la dimension du contrôleur. Nous démontrons que la loi de commande s'appuyant sur un modèle réduit peut être appliquée au système complet. Dans ce cas la densité d'énergie cinétique décroit presque de façon monotone au cours du temps en utilisant une commande par asservissement visuel à deux degrés de liberté. / The visual servoing control approach is formulated for the flow control of the plane Poiseuille flow. Generally, the flow control can lead the flow from its current state to a desired state. In transition to turbulence, the growth of kinetic energy density can lead the flow to turbulence. Moreover, the drag reduction is a potential application in the engineering applications. Therefore, this thesis aims to minimize the kinetic energy density and the skin friction drag. The governing equations of the plane Poiseuille flow are modeled to a standard form in the automatic control. More precisely, the partial differential equations of the plane Poiseuille flow are transformed to a state space representation by using the spectral method. The streamwise and spanwise directions are discretized based on the Fourier series while the wall-normal direction is discretized based on the Chebyshev polynomials. The state vector involves the wall-normal velocity and vorticity. The control signals depend on the inhomogeneous Dirichlet boundary conditions which correspond to blowing/suction boundary control. The number of independent control signals is called the number of the degree of freedom. Moreover, the skin-friction drag and the kinetic energy density are modeled as a function of the state vector. The goal is to minimize both the skin-friction drag and the kinetic energy density by appropriate methods. The partitioned visual servoing control is used to minimize, simultaneously, the skin-friction drag and the kinetic energy density with two degrees of freedom. As a result, the behavior of the skin-friction drag monotonically decreases in time. However, the behavior of the kinetic energy density does not monotonically decrease in time, the similar results from the other methods such as: PID and LQR controls. Therefore, the number of the degree of freedom increases, which leads to the improvement of the kinetic energy density. In addition, when the number of the degree of freedom equals the number of state vector, the kinetic energy density monotonically decreases in time by using the visual servoing control. The dimension of linearized plane Poiseuille flow is large, therefore, we need to reduce the order of controller. We demonstrate that the control law based on a mode reduction can be applied for the full system. Moreover, the kinetic energy density almost will monotonically decreases in time even using two degrees of freedom when the visual servoing control is designed based on the model order reduction.

Dans cette thèse, nous présentons notre travail de recherche dans le domaine du calcul haute performance en mécanique des fluides. Avec la demande croissante de simulations à haute résolution, il est devenu important de développer des solveurs numériques pouvant tirer parti des architectures récentes comprenant des processeurs multi-cœurs et des accélérateurs. Nous nous proposons dans cette thèse de développer un solveur efficace pour la résolution sur architectures hétérogènes CPU/GPU des équations de Navier-Stokes (NS) relatives aux écoulements 3D de fluides incompressibles.Tout d'abord nous présentons un aperçu de la mécanique des fluides avec les équations de NS pour fluides incompressibles et nous présentons les méthodes numériques existantes. Nous décrivons ensuite le modèle mathématique, et la méthode numérique choisie qui repose sur une technique de prédiction-projection incrémentale.Nous obtenons une distribution équilibrée de la charge de calcul en utilisant une méthode de décomposition de domaines. Une parallélisation à deux niveaux combinée avec de la vectorisation SIMD est utilisée dans notre implémentation pour exploiter au mieux les capacités des machines multi-cœurs. Des expérimentations numériques sur différentes architectures parallèles montrent que notre solveur NS obtient des performances satisfaisantes et un bon passage à l'échelle.Pour améliorer encore la performance de notre solveur NS, nous intégrons le calcul sur GPU pour accélérer les tâches les plus coûteuses en temps de calcul. Le solveur qui en résulte peut être configuré et exécuté sur diverses architectures hétérogènes en spécifiant le nombre de processus MPI, de threads, et de GPUs.Nous incluons également dans ce manuscrit des résultats de simulations numériques pour des benchmarks conçus à partir de cas tests physiques réels. Les résultats obtenus par notre solveur sont comparés avec des résultats de référence. Notre solveur a vocation à être intégré dans une future bibliothèque de mécanique des fluides pour le calcul sur architectures parallèles CPU/GPU. / In this PhD thesis, we present our research in the domain of high performance software for computational fluid dynamics (CFD). With the increasing demand of high-resolution simulations, there is a need of numerical solvers that can fully take advantage of current manycore accelerated parallel architectures. In this thesis we focus more specifically on developing an efficient parallel solver for 3D incompressible Navier-Stokes (NS) equations on heterogeneous CPU/GPU architectures. We first present an overview of the CFD domain along with the NS equations for incompressible fluid flows and existing numerical methods. We describe the mathematical model and the numerical method that we chose, based on an incremental prediction-projection method.A balanced distribution of the computational workload is obtained by using a domain decomposition method. A two-level parallelization combined with SIMD vectorization is used in our implementation to take advantage of the current distributed multicore machines. Numerical experiments on various parallel architectures show that this solver provides satisfying performance and good scalability.In order to further improve the performance of the NS solver, we integrate GPU computing to accelerate the most time-consuming tasks. The resulting solver can be configured for running on various heterogeneous architectures by specifying explicitly the numbers of MPI processes, threads and GPUs. This thesis manuscript also includes simulation results for two benchmarks designed from real physical cases. The computed solutions are compared with existing reference results. The code developed in this work will be the base for a future CFD library for parallel CPU/GPU computations.

Este trabalho de pesquisa é dedicado ao desenvolvimento, análise e implementação de um novo esquema upwind de alta resolução (denominada PFDPUS) para a aproximação de termos convectivos em leis de conservação e problemas relacionados em mecânica dos fluídos. Usando variáveis normalizadas de Leonard, o equema PFDPUS é baseado em uma função polinomial por partes que satisfaz os critérios de estabilidade CBC e TVD. O desempenho do esquema PEDPUS é investigado na solução das equações de advecção de escalares, Burgers, Euler e MHD. O novo esquema é então aplicado para simular escoamentos incompressíveis envolvendo superfícies livres móveis. Para tanto, o esquema PFDPUS é implementado dentro do software CLAWPACK para problemas compressíveis, e no código Freeflow para poblemas incompressíveis. Os resultados numéricos são comparados com dados experimentais, numéricos e analíticos / This work is dedicated to the development, analysis and implementation of a new high-resolution upwind scheme (called PFDPUS) for approximation of convective terms in conservation laws and related fluid mechanics problems. By using the normalized variables of Leonard, the PFDPUS scheme is based on a piecewise polynomical function that satisfies the CBC and TVD stability criteria. The performance of the PFDPUS scheme is assessed by solving advection of scalars, Burgers, Euler and MHD equations. Then the new scheme is applied to simulate incompressible flows involving moving free surfaces. The PFDPUS scheme is implemented into the CLAWPACK software for compressible problems, and in the Freeflow code for incompressible problems. The numerical results are compared with experimental, numerical and analytical data