L’objectif de cette thèse est de développer un modèle numérique 3D afin d’étudier le phénomène de laves torrentielles ; écoulements visqueux fortement chargés en matière solide, surgissant en montagne lors d’orages violents. Aujourd’hui, la prévision des zones de vulnérabilité s’appuie sur des outils de calcul 0D, 1D ou 2D. Or ces outils ne peuvent représenter pleinement le comportement à surface libre des écoulements du fait de nombreuses approximations et hypothèses. Ainsi cette thèse met en œuvre un code numérique 3D pour étudier ces écoulements. Ce travail aboutit à la construction d’un modèle 3D à partir de données réelles de terrain. Plusieurs scénarios ont été étudiés et comparés à des résultats issus d’un modèle 2D. Les résultats mettent en évidence les apports non négligeables de la modélisation 3D : zones d’étalement et de dépôt, phénomènes d’encombrement, modélisation fine des écoulements dans les zones chenalisées. Préalablement, le modèle 3D a été validé en comparant les résultats numériques à des données expérimentales issues de la littérature, pour des typologies d’écoulement représentatives de celles observées sur des sites grandeur nature. / The objective of this thesis is to develop a 3D numerical model to assess debris flow. These viscous flows, heavily loaded with solid matter, form when heavy rain occurs in mountains. Today, forecasts of potentially impacted areas are based on 0D, 1D and 2D numerical tools. However, these tools cannot fully represent the free surface behaviour of debris flows due to the approximations and assumptions on which they are based. Thus, this work utilises a 3D numerical code to study this phenomenon. A specific model is built with real field data. Several flow scenarios are studied and compared with a 2D numerical model. The results highlight the significant benefits of a 3D approach by providing information on the fine representation of flow dynamics over the catchment area. The model also predicts the impact of debris flow (overflowing on a road bridge) and the zones of deposition and spreading. It highlights possible congestion phenomena and reproduces flows in the channels by fully accounting for parietal friction, capabilities not provided by 2D models. Prior to this application, the 3D model was evaluated with five sets of experimental data to validate its ability to represent viscoplastic flows. Different types of flows are studied and are representative of those observed on real sites when debris flow occur.

Esta dissertação trata da solução numérica de problemas em dinâmica dos fluidos usando dois novos esquemas upwind de alta resolução, denominados FDPUS-C1 (Five-Degree Polynomial Upwind Scheme of \' C POT. 1\' Class) e SDPUS-C1 (Six-Degree Polynomial Upwind Scheme of \'C POT.1\' Class), para a discretização de termos convectivos lineares e não-lineares. Os esquemas são baseados nos critérios de estabilidade TVD (Total Variation Diminishing) e CBC (Convection Boundedness Criterion) e são implementados, nos contextos das metodologias de diferenças finitas e volumes finitos, no ambiente de simulação Freeflow (an integrated simulation system for Free surface Flow) para escoamentos imcompressíveis 2D, 2D-1/2 e 3D, ou no código bem conhecido CLAWPACK ( Conservation LAW PACKage) para problemaw compressíveis 1D e 2D. Vários testes computacionais são feitos com o objetivo de verificar e validar os métodos numéricos contra esquemas upwind populares. Os novos esqumas são então aplicados na resolução de uma gama ampla de problemas em CFD (Computational Fluids Dynamics), tais como propagação de ondas de choque e escoamentos incompressíveis envolvendo superfícies livres móveis. Em particular, os resultados numéricos para leis de conservação hiperbólicas 2D e equações de Navier-Stokes incompressíveis 2D, 2D-1/2 e 3D demosntram que esses novos esquemas convectivos tipo upwind polinomiais funcionam muito bem / This dissertation deals with the numerical solution of fluid dynamics problems using two new high resolution upwind schemes,. namely FDPUS-C1 and SDPUS-C1, for the discretization of the linear and non-linear convection terms. The Schemes are based on TVD and DBC stability criteria and are implemented in the context of the finite difference and finite volume methodologies, either into the Freeflow code for 2D, 2D-1/2 and 3D incompressible flows or in the well-known CLAWPACK code for 1D and 2D compressible flows. Several computational tests are performed to verify and validate the numerical methods against other popularly used upwind schemes. The new schemes are then applied to solve a wide range of problems in CFD, such as shock wave propagation and incompressible fluid flows involving moving free msurfaces. In particular, the numerical results for 2D hyperbolic conservation laws and 2D, 2D-1/2 and 3D incompressible Navier-Stokes eqautions show that new polynomial upwind convection schemes perform very well

Un des enjeux de la Directive Cadre sur l'Eau (2000) est de rétablir la migration piscicole aval à travers notamment une baisse de la mortalité des poissons au niveau des turbines des centrales hydroélectriques. Les grilles ichtyocompatibles constituent une des solutions possibles grâce à un arrêt physique des poissons couplé à un guidage vers un exutoire de contournement. Dans ce contexte, nous avons mené une étude expérimentale et numérique afin de caractériser différents types de grilles en termes de pertes de charge et de courantologie. Cela a permis de déterminer d'une part les pertes d'énergie liées à la proximité des barreaux pour diverses configurations et d'autre part les vitesses que les poissons sont susceptibles de rencontrer le long des grilles jusqu'à l'exutoire. Ces deux sujets complémentaires ont été abordés via l'utilisation de deux installations hydrauliques couplées à plusieurs systèmes de mesure adaptés à nos modèles de grille. Des mesures de hauteur d'eau ont permis d'établir des nouvelles lois de pertes de charge tandis que des mesures de vitesse avec des systèmes de vélocimétrie acoustique par effet Doppler (ADV) et de vélocimétrie par images de particules (PIV) ont permis de déterminer les angles pour lesquels les grilles satisfont les critères d'ichtyocompatibilité. Certains points de l'étude, liés à la caractérisation des exutoires notamment, ont été résolus grâce à des simulations numériques avec le logiciel OpenFOAM. L'ensemble des résultats obtenus a permis de dimensionner et de comparer les différentes solutions de grilles dans un cas de prise d'eau réel, soulignant ainsi les avantages et problématiques inhérents à chacune des solutions. / The Water Framework Directive (2000) particularly aimed at restoring the fish downstream migration and at preventing fish from being injured or killed in hydropower plant turbines. Fish-friendly trashracks, which physically stop fish and guide them towards bypasses, are one of the possible solutions.In this context, we have carried out an experimental and numerical study focused on the characterization of head losses and velocity distributions for various kinds of trashrack. This helped to quantify the energy loss due to the narrow bar spacing and to assess the magnitude of velocities that fish may face in different trashrack configurations. Two hydraulic installations coupled with several measurement devices allowed the study of both topics on model trashracks. Water depth measurements resulted in new head loss equations while velocity measurements with acoustic Doppler velocimetry (ADV) and particle image velocimetry (PIV) systems allowed the determination of fish-friendly angles in different trashrack configurations. Some other points of the thesis, concerning for instance bypasses, have been investigated with numerical simulations using OpenFOAM software. All the results allowed the conception and the comparison of various trashrack solutions for a real water intake case, highlighting the advantages and drawbacks relative to each of these solutions.

Esta dissertação trata da solução numérica de problemas em dinâmica dos fluidos usando dois novos esquemas upwind de alta resolução, denominados FDPUS-C1 (Five-Degree Polynomial Upwind Scheme of \' C POT. 1\' Class) e SDPUS-C1 (Six-Degree Polynomial Upwind Scheme of \'C POT.1\' Class), para a discretização de termos convectivos lineares e não-lineares. Os esquemas são baseados nos critérios de estabilidade TVD (Total Variation Diminishing) e CBC (Convection Boundedness Criterion) e são implementados, nos contextos das metodologias de diferenças finitas e volumes finitos, no ambiente de simulação Freeflow (an integrated simulation system for Free surface Flow) para escoamentos imcompressíveis 2D, 2D-1/2 e 3D, ou no código bem conhecido CLAWPACK ( Conservation LAW PACKage) para problemaw compressíveis 1D e 2D. Vários testes computacionais são feitos com o objetivo de verificar e validar os métodos numéricos contra esquemas upwind populares. Os novos esqumas são então aplicados na resolução de uma gama ampla de problemas em CFD (Computational Fluids Dynamics), tais como propagação de ondas de choque e escoamentos incompressíveis envolvendo superfícies livres móveis. Em particular, os resultados numéricos para leis de conservação hiperbólicas 2D e equações de Navier-Stokes incompressíveis 2D, 2D-1/2 e 3D demosntram que esses novos esquemas convectivos tipo upwind polinomiais funcionam muito bem / This dissertation deals with the numerical solution of fluid dynamics problems using two new high resolution upwind schemes,. namely FDPUS-C1 and SDPUS-C1, for the discretization of the linear and non-linear convection terms. The Schemes are based on TVD and DBC stability criteria and are implemented in the context of the finite difference and finite volume methodologies, either into the Freeflow code for 2D, 2D-1/2 and 3D incompressible flows or in the well-known CLAWPACK code for 1D and 2D compressible flows. Several computational tests are performed to verify and validate the numerical methods against other popularly used upwind schemes. The new schemes are then applied to solve a wide range of problems in CFD, such as shock wave propagation and incompressible fluid flows involving moving free msurfaces. In particular, the numerical results for 2D hyperbolic conservation laws and 2D, 2D-1/2 and 3D incompressible Navier-Stokes eqautions show that new polynomial upwind convection schemes perform very well

Aujourd’hui, les zones côtières concentrent souvent des villes densément peuplées où le développement économique est associé à une rapide urbanisation. Dans de nombreux endroits, les cours d’eau sont intégrés en zone urbaine présentant à la fois des ressources et des risques potentiels. La concurrence croissante entre les activités économiques et l'espace essentiel aux rivières et aux milieux naturels nécessite une gestion intégrée basée sur des outils fiables capables de fournir des informations hydrologiques. La basse vallée du Var, dans le sud de la France, est spécifiquement exposée à une augmentation de l’activité humaine, et les impacts sociaux sur l’environnement deviennent de plus en plus intenses. Par conséquent, le cycle de l’eau est déséquilibré et induit des problèmes liés à l’eau. Pour traiter ces problèmes, les collectivités territoriales ont besoin d’outils d’aide à la décision qui simulent le comportement du bassin versant. Tel est l’objectif du projet AquaVar qui repose sur trois modèles numériques : Mike SHE pour le bassin versant du Var, Mike 21FM pour la rivière, et Feflow pour l’aquifère. Ces travaux de recherches se concentrent sur le modèle des écoulements à surface libre et son intégration au sein de l’outil d’aide à la décision.Les modèles hydrauliques 2D sont fréquemment utilisés car ils fournissent une vision précise des phénomènes physiques en rivière et de l’hydrodynamique durant des événements extrêmes (inondation, sécheresse et pollution accidentelle). Le modèle des écoulements à surface libre 2D est conçu avec le logiciel Mike 21FM en utilisant la méthodologie classique. Il est calibré et validé pour des conditions spécifiques de la basse vallée du Var : hautes eaux au mois de novembre, hautes eaux au printemps dues à la fonte des neiges et basses eaux en été. Le modèle hydraulique peut non seulement être utilisé pour simuler des scénarios d’événements extrêmes, mais il est également capable d’interagir avec le modèle des écoulements souterrains développé avec Feflow. Une interface spécifique est développée pour connecter le Var et son aquifère dans la basse vallée. Elle nécessite cependant une phase de validation afin de précisément simuler les échanges nappe-rivière.De toute évidence, le modèle Mike 21FM est le cœur de l’outil d’aide à la décision. L'hydrodynamique du fleuve dépend de la transformation des précipitations en débit (modélisée par Mike SHE) et du volume d'eau échangé avec l'aquifère (modélisé par Feflow). La dernière partie est consacrée à la construction de l'outil d’aide à la décision avec l'intégration du modèle hydraulique 2D. Ce système de modèles est implémenté dans une interface web et accessible par différents utilisateurs. Cet outil peut être facilement mis à jour en partageant les données de chaque partenaire. Dans le futur, cet outil pourrait intégrer le petit cycle de l'eau (réseaux d’assainissement et d’alimentation en eau potable). Les autorités locales pourraient donc prévoir et contrôler leurs ressources en eau et les risques potentiels. / Nowadays, the coastal areas concentrate frequently densely populated cities where economic development is associated to a fast urbanization. In many locations rivers are integrated within the urban areas and present both resources and potential risks. The growing competition between the economic activities and the essential space for the rivers requests an efficient management based on reliable tools able to provide hydrological information. Typically, the Lower Var valley, in south of France, faces to an increase of human activities, and the social impacts on the natural environment have become more and more intensive. Therefore, the water cycle is unbalanced implying water problems. To deal with this issue, local authorities needs a Decision Support System (DSS) tool to simulate the behaviour of water system. This is the objective of the AquaVar research which includes three numerical models: Mike SHE for the Var catchment, Mike 21FM for the Lower Var River, and Feflow for the aquifer. Here, the research is focused on the 2D free surface flow model and its integration in the DSS tool.The 2D hydraulic models represent a meaningful approach that can provide an accurate view on the physical processes within the river and on the hydrodynamics during the extreme events (inundation, drought and accidental pollution). The 2D free surface flow model is designed with Mike 21FM software using common methodology. It is calibrated and validated for the three specific weather conditions in the Lower Var valley: floods which occur in November, spring floods due to snow melting and droughts in summer. Not only the 2D hydraulic model allows to simulate scenarios of inundation and accidental pollutions, but it is also able to exchange water volume with the groundwater model developed with Feflow. A specific interface is developed to connect river and aquifer in the Lower Var valley. It required validation cases to accurately simulate the river-aquifer exchanges.Obviously, Mike 21FM is finally the core of the DSS tool because it is the centre of the modelling system. The hydrodynamics of the river depends on precipitations converted in flow by Mike SHE model, and exchanges water volume with the aquifer modelled by Feflow. The last part of this research is dedicated to the construction of the DSS tool and the integration of the 2D hydraulic model. The modelling system is implemented in a web interface adapted for different types of users. This DSS tool can be easily updated by sharing data from all the stakeholders. In the future, this tool could integrate sewage network and drinking water supply system to consider all the water cycle. Hence, local authorities could forecast and control the water resources and the potential risks.

Neste trabalho são apresentados os resultados do desenvolvimento e teste de esquemas upwind de alta resolução para o controle da difusão numérica em leis de conservação gerais e problemas em dinâmica dos fluidos. Em particular, são derivados dois novos esquemas: o ALUS (Adaptive Linear Upwind Scheme) e o TOPUS (Third-Order Polynomial Upwind Scheme). Esses esquemas são testados no transporte de escalares, em equações 1D tipo convecção-difusão, em sistemas hiperbólicos 1D, nas equações de Euler 2D da dinâmica dos gases e nas equações de Navier-Stokes incompressíveis 2D/3D. Os esquemas são então associados a uma modelagem algébrica não linear para a simulação de problemas de escoamentos incompressíveis turbulentos 2D com/sem superfícies livres / In this work, results of the development and testing of high-resolution upwind schemes for controlling of the numerical diffusion for general conservation laws and fluid dynamics problems are presented. In particular, two new high-resolution upwind schemes are derived, namely, the ALUS (Adaptive Linear Upwind Scheme) and the TOPUS (Third-Order Polynomial Upwind Scheme). These schemes are tested in scalar transport, 1D convection-diffusion equations, 1D hyperbolic systems, 2D Euler equations of the gas dynamics, and in 2D/3D incompressible Navier-Stokes equations. The schemes are then combined with a nonlinear Reynolds stress algebraic equation model for the simulation of 2D incompressible turbulent flows with/without free surfaces

Uma parte significativa dos escoamentos encontrados em aplicações tecnológicas é caracterizada por envolver altos números de Reynolds, principalmente aqueles em regime turbulento e com superfície livre. Obter soluções numéricas representativas para essa classe de problemas é extremamente difícil, devido à natureza não-linear das equações diferenciais parciais envolvidas nos modelos. Conseqüentemente, o tema tem sido uma das principais preocupações da comunidade científica moderna em dinâmica de fluidos computacional. Aproximações de primeira ordem para os termos convectivos são as mais adequadas para amortecer oscilações que estão associadas às aproximações de alta ordem não-limitadas. Todavia, elas introduzem dissipação artificial nas representações discretas comprometendo os resultados numéricos. Para minimizar esse efeito não-físico e, ao mesmo tempo, conseguir aproximações incondicionalmente estáveis, é indispensável adotar uma estratégia que combine aproximações de primeira ordem com as de ordem mais alta e que leve em conta a propagação de informações físicas. Os resultados dessa composição são os esquemas "upwind" limitados de alta ordem. Em geral, espera-se que esses esquemas sejam apropriados para a representação das derivadas convectivas nos modelos de turbulência kappa-varepsilon. No contexto de diferenças finitas, a presente tese dedica-se à solução numérica das equações de Navier-Stokes no regime de números de Reynolds elevados. Em particular, ela contém uma análise de algoritmos monotônicos e antidifusivos e modelos de turbulência kappa-varepsilon para a simulação de escoamentos incompressíveis envolvendo superfícies livres. Esquemas de convecção são implementados nos códigos GENSMAC para proporcionar um tratamento robusto dos termos convectivos nas equações de transporte. Duas versões do modelo kappa-varepsilon de turbulência são implementadas nos códigos GENSMAC, para problems bidimensionais e com simetria radial, para descrever os efeitos da turbulência sobre o escoamento médio. Resultados numéricos de escoamentos com simetria radial são comparados com resultados experimentais e analíticos. Simulações numéricas de problemas tridimensionais complexos são apresentadas para avaliar o desempenho de esquemas "upwind". Finalmente, os modelos de turbulência kappa-varepsilon são utilizados para a simulação de escoamentos confinados e com superfícies livres. / A considerable part of fluid flows encountered in technological applications is characterised by involving high-Reynolds numbers, especially those in turbulent regime and with free-surface. It is extremely difficult to obtain representative numerical solutions for this class of problems, due to the non-linear nature of the partial differential equations involved in the models. Consequently, this subject has been one of main concerns in the modern computational fluid dynamics community. First-order approximation to the convective terms is one of the most appropriate to smooth out oscilations/instabilities which are associated with high-order unlimited approximation. However, it introduces numerical dissipation in the discrete representation jeopardizing the numerical results. In order to minimize this non-physical effect and, at the same time, to obtain unconditionally stable approximation, it is essential to adopt a strategy that combines first and high-order approximations and takes into account the propagation of physical information. The results of this composition are the high-order bounded upwind techniques. In general, it is expected that these algorithms are satisfactory for the representation of the convective derivatives in the kappa-varepsilon turbulence model. In the context of finite-difference, the present thesis deals with the numerical solution of the Navier-Stokes equations at high-Reynolds number regimes. In particular, it contains an analysis of monotonic and anti-difusive convection schemes and kappa-varepsilon turbulence models for the simulation of free-surface fluid flows. Upwinding methods are implemented into the GENSMAC codes to provide a robust treatment of the convective terms in the transport equations. Two versions of the K-Epsilon turbulence model are implemented into the two-dimensional and axisymmetric GENSMAC codes, in order to describe the turbulent effects on the average flow. Numerical results of axisymmetric flows are compared with experimental and analytical results. Numerical simulations of complex three-dimensional problems are presented to assess the performance of high-order bounded upwind schemes. Finally, the K-Epsilon turbulence models are employed in the simulation of confined and free-surface flows.

In this work we review 12 years of developments in the field of residual based discretizations for hyperbolic problems and their application to the solution of the shallow water equations. Fundamental concepts related to the topic are recalled and he construction of second and higher order schemes for steady problems is presented. The generalization to time dependent problems by means of multi-step implicit time integration, space-time, and genuinely explicit techniques is thoroughly discussed. Finally, the issues of C-property, super consistency, and wetting/drying are analyzed in this framework showing the power of the residual based approach.

Uma parte significativa dos escoamentos encontrados em aplicações tecnológicas é caracterizada por envolver altos números de Reynolds, principalmente aqueles em regime turbulento e com superfície livre. Obter soluções numéricas representativas para essa classe de problemas é extremamente difícil, devido à natureza não-linear das equações diferenciais parciais envolvidas nos modelos. Conseqüentemente, o tema tem sido uma das principais preocupações da comunidade científica moderna em dinâmica de fluidos computacional. Aproximações de primeira ordem para os termos convectivos são as mais adequadas para amortecer oscilações que estão associadas às aproximações de alta ordem não-limitadas. Todavia, elas introduzem dissipação artificial nas representações discretas comprometendo os resultados numéricos. Para minimizar esse efeito não-físico e, ao mesmo tempo, conseguir aproximações incondicionalmente estáveis, é indispensável adotar uma estratégia que combine aproximações de primeira ordem com as de ordem mais alta e que leve em conta a propagação de informações físicas. Os resultados dessa composição são os esquemas "upwind" limitados de alta ordem. Em geral, espera-se que esses esquemas sejam apropriados para a representação das derivadas convectivas nos modelos de turbulência kappa-varepsilon. No contexto de diferenças finitas, a presente tese dedica-se à solução numérica das equações de Navier-Stokes no regime de números de Reynolds elevados. Em particular, ela contém uma análise de algoritmos monotônicos e antidifusivos e modelos de turbulência kappa-varepsilon para a simulação de escoamentos incompressíveis envolvendo superfícies livres. Esquemas de convecção são implementados nos códigos GENSMAC para proporcionar um tratamento robusto dos termos convectivos nas equações de transporte. Duas versões do modelo kappa-varepsilon de turbulência são implementadas nos códigos GENSMAC, para problems bidimensionais e com simetria radial, para descrever os efeitos da turbulência sobre o escoamento médio. Resultados numéricos de escoamentos com simetria radial são comparados com resultados experimentais e analíticos. Simulações numéricas de problemas tridimensionais complexos são apresentadas para avaliar o desempenho de esquemas "upwind". Finalmente, os modelos de turbulência kappa-varepsilon são utilizados para a simulação de escoamentos confinados e com superfícies livres. / A considerable part of fluid flows encountered in technological applications is characterised by involving high-Reynolds numbers, especially those in turbulent regime and with free-surface. It is extremely difficult to obtain representative numerical solutions for this class of problems, due to the non-linear nature of the partial differential equations involved in the models. Consequently, this subject has been one of main concerns in the modern computational fluid dynamics community. First-order approximation to the convective terms is one of the most appropriate to smooth out oscilations/instabilities which are associated with high-order unlimited approximation. However, it introduces numerical dissipation in the discrete representation jeopardizing the numerical results. In order to minimize this non-physical effect and, at the same time, to obtain unconditionally stable approximation, it is essential to adopt a strategy that combines first and high-order approximations and takes into account the propagation of physical information. The results of this composition are the high-order bounded upwind techniques. In general, it is expected that these algorithms are satisfactory for the representation of the convective derivatives in the kappa-varepsilon turbulence model. In the context of finite-difference, the present thesis deals with the numerical solution of the Navier-Stokes equations at high-Reynolds number regimes. In particular, it contains an analysis of monotonic and anti-difusive convection schemes and kappa-varepsilon turbulence models for the simulation of free-surface fluid flows. Upwinding methods are implemented into the GENSMAC codes to provide a robust treatment of the convective terms in the transport equations. Two versions of the K-Epsilon turbulence model are implemented into the two-dimensional and axisymmetric GENSMAC codes, in order to describe the turbulent effects on the average flow. Numerical results of axisymmetric flows are compared with experimental and analytical results. Numerical simulations of complex three-dimensional problems are presented to assess the performance of high-order bounded upwind schemes. Finally, the K-Epsilon turbulence models are employed in the simulation of confined and free-surface flows.

Neste trabalho são apresentados os resultados do desenvolvimento e teste de esquemas upwind de alta resolução para o controle da difusão numérica em leis de conservação gerais e problemas em dinâmica dos fluidos. Em particular, são derivados dois novos esquemas: o ALUS (Adaptive Linear Upwind Scheme) e o TOPUS (Third-Order Polynomial Upwind Scheme). Esses esquemas são testados no transporte de escalares, em equações 1D tipo convecção-difusão, em sistemas hiperbólicos 1D, nas equações de Euler 2D da dinâmica dos gases e nas equações de Navier-Stokes incompressíveis 2D/3D. Os esquemas são então associados a uma modelagem algébrica não linear para a simulação de problemas de escoamentos incompressíveis turbulentos 2D com/sem superfícies livres / In this work, results of the development and testing of high-resolution upwind schemes for controlling of the numerical diffusion for general conservation laws and fluid dynamics problems are presented. In particular, two new high-resolution upwind schemes are derived, namely, the ALUS (Adaptive Linear Upwind Scheme) and the TOPUS (Third-Order Polynomial Upwind Scheme). These schemes are tested in scalar transport, 1D convection-diffusion equations, 1D hyperbolic systems, 2D Euler equations of the gas dynamics, and in 2D/3D incompressible Navier-Stokes equations. The schemes are then combined with a nonlinear Reynolds stress algebraic equation model for the simulation of 2D incompressible turbulent flows with/without free surfaces