Continuous milli-scale and micro-scale structures such as FlowPlate® microreactors have emerged as a promising element of process intensification due to their inherently effective rates of mass and heat transfer. These microfluidic devices have proven to be a preferred solution in place of energy-intensive batch processes for certain pathways of fine chemical and pharmaceutical synthesis, most notably fast reactions taking place on the scale of milliseconds to seconds. Computational fluid dynamics (CFD) has become an increasingly valuable tool in the field of microreactor design and optimization for its ability to locally map complex fluid flow patterns and resolve microscopic scales of reactive mixing that are challenging to characterize experimentally. The primary objective of this research was thus to develop and validate a mathematical model for the simulation of chaotic flow and homogeneous mixing in continuous microreactors. The model needed to be versatile enough to handle transition between flow regimes within a given reactor as well as the coexistence of both chaotic and laminar flow patterns in the micromixing elements that comprise said reactors. This was successfully achieved through the implementation of a k-ω SST (shear-stress transport) turbulence model that accounts for the impact of small-scale temporal and spatial fluctuations generated in the micromixer geometries studied herein; be it a liquid-liquid mixer (LLM), a serpentine (SZ) or a tangential (TG) mixer. In a first CFD study, the computational predictions were validated based on excellent agreement with experimental pressure loss (R^2 > 0.997) and residence time distribution (RTD) data (R^2 > 0.97) in several LL microreactors at Reynolds numbers ranging from 210 to 2140. Furthermore, the local velocity distribution and streamlines were mapped across the 3D domain of these reactors and it was discovered, based on the emergence of advective recirculation zones and turbulent dispersion, that a drastic change in flow behaviour occurred in these mixing elements at a Reynolds number of about 640. The interspacing of LLM elements with straight microchannels proved to be a suitable approach to modulating pressure loss while concurrently maintaining the chaotic secondary flow patterns generated from the mixers. In a second CFD study, the impact of micromixer geometry on the local velocity fields and advective transport performance was investigated both from a macromixing and micromixing perspective. Like the LLM, the SZ and TG mixers conferred chaotic secondary flow patterns at characteristic Reynolds numbers between 500 and 1000. As such, it was concluded that it would be ideal to operate these mixers at water flow rates of at least 30 ml/min. Contour plots of the velocity magnitude coupled with the computation of RTD showed that the SZ virtually mimics a plug-flow profile over a volume of 77 mm3 or greater at 50 g/min. The RTD of the LLM and TG resembles that of a mixed flow pattern given that approximately 65-80% of their fluid volume is occupied by recirculation zones. As such, it required 65 LLMs in series (3105 mm3) and 80 TGs (1142 mm3) to approach the same pattern as 10 SZs (77 mm3) from a macromixing perspective. Micromixing time distributions (MTD) were also characterized by locally computing the decay time of small-scale segregation (t_SSS) as a function of flow rate, wherein higher flow rates generated lower characteristic mixing times. The TG and LLM conferred the broadest range of mixing times, spanning nearly four orders of magnitude in the range of [0.02 ms, 10 ms], whereas the SZ generated a much narrower MTD ranging between [0.024 ms, 0.69 ms]. Finally, the impact of geometry and flow conditions on reaction yield was assessed by characterizing the extent of a finite-rate reaction relative to an infinitely fast reaction taking place in parallel. The calculated yield for the competitive-parallel reaction scheme showed that the second Damköhler number (Dall) computed based on the mean tSSS provides useful information about whether the process will be limited by the intrinsic rate of reaction or by the rate of mass transfer, even though the reaction process is controlled by a combination of the RTD as well as loss of LSS and SSS. It was concluded that the change in MTD as a function of power dissipation should coincide with the reaction yield response, and that any deviation in that relationship is because of macroscopic blending of reactants in the entrance region.

Nesta tese é apresentado um esquema novo de alta resolução upwind (denominado TDPUS-C3) para reconstrução de fluxos numéricos para leis de conservação não lineares e problemas relacionados em DFC. O esquema é baseado nos critérios de estabilidade CBC e TVD e desenvolvido utilizando condições de diferenciabilidade \'C POT. 3\'. Além disso, é realiozada a implementação da associação do esquema TDPLUS-C3 com a modelagem de turbulência RNG \'\\kappa - \\epsilon\'. O propósito é obter soluções numéricas de sistemas hiperbólicos de leis de conservação para dinâmica dos gases e equações de Navier-Stokes para escoamento incompreensível de fluidos newtonianos e não newtonianos (viscoelásticos). Fazendo o uso do esquema TDPUS-C3, a precisão global dos métodos numéricos é verificada acessando o erro em problemas teste (benchmark) 1D e 2D. Um estudo comparativo entre os resultados do esquema TDPUS-C3 e os esquemas upwind convencionais para leis de conservação hiperbólicas complexas é também realizado. A Associação das modelagens numéricas (upwinding mais RNG \'\\kappa - \\epsilon\') é , então, examinada na simulação de escoamentos turbulentos de fluidos newtonianos envolvendo superfícies livres móveis, usando a metodologia URANS. No geral, em termos do comportamento global, concordância satisfatória é observada / In this thesis, a new high-resolution upwind scheme (named TDPUS-C3) for reconstruction of numerical fluxes for nonlinear conservation laws and related CFD problems in presented. The scheme is based on CBC and TVD stability criteria and developed by employing differentiability condictions (\'C POT. 3\'). In additon, the implementation of an association of the TDPUS-C3 scheme with the RNG \'\\kappa - \\epsilon\' turbulence modelling is also performed. The purpose is to obtain numerical solutions of systems of hyperbolic conservation laws for gas dynamics and Navier-Stokes equations for incompressible flow of Newtonian and non-Newtonian (viscoelstic) fluids. By using the TDPUS-C3 scheme, the global accuracy of the numerical methods is verified by assessing the error on 1D and 2D benchmark test cases. A comparative study between the TDPUS-C3 scheme and convectional upwind schemes to solve standard and complex hyperbolic conservation laws is also accomplished. The association of the numerical modelling (upwinding plus RNG \'\\kappa - epsilon\') is then examined in the simulation of turbulent Newtonian fluid flows involving moving free surfaces, by using URANS methodology. Overall, satisfactory agreement is found in terms of the overall behaviour

There is currently a great deal of interest within the aviation community in the design of small, slow-flying but manoeuvrable uninhabited vehicles for reconnaissance, surveillance, and search and rescue operations in urban environments. Inspired by observation of birds, insects, fish and cetaceans, flapping wings are being actively studied in the hope that they may provide greater propulsive efficiencies than propellers and rotors at low Reynolds numbers for such Micro-Air Vehicles (MAVs). Researchers have posited the Strouhal number (combining flapping frequency, amplitude and forward speed) as the parameter controlling flapping wing aerodynamics in cruising flight, although there is conflicting evidence. This thesis explores the effect of flapping frequency and amplitude on forces and wake structures, as well as physical mechanisms leading to optimum propulsive efficiency. Two-dimensional rigid airfoils are considered at Reynolds number 2,000 ??? 40,000. A compressible Navier-Stokes simulation is combined with numerical and analytical potential flow techniques to isolate and evaluate the effect of viscosity, leading and trailing edge vortex separation, and wake vortex dynamics. The wake structures of a plunging airfoil are shown to be sensitive to the flapping frequency independent of the Strouhal number. For a given frequency, the wake of the airfoil exhibits ???vortex lock-in??? as the amplitude of motion is increased, in a manner analogous to an oscillating circular cylinder. This is caused by interaction between the flapping frequency and the ???bluff-body??? vortex shedding frequency apparent even for streamlined airfoils at low Reynolds number. The thrust and propulsive efficiency of a plunging airfoil are also shown to be sensitive to the flapping frequency independent of Strouhal number. This dependence is the result of vortex shedding from the leading edge, and an interaction between the flapping frequency and the time for vortex formation, separation and convection over the airfoil surface. The observed propulsive efficiency peak for a pitching and plunging airfoil is shown to be the result of leading edge vortex shedding at low flapping frequencies (low Strouhal numbers), and high power requirements at large flapping amplitudes (high Strouhal numbers). The efficiency peak is governed by flapping frequency and amplitude separately, rather than the Strouhal number directly.

Wenn auf dem Foto oben der Westwind Zephyr und in seinen Armen die Morgenbrise Aura nicht Venus an die Gestaden Zyperns treibt, sondern stattdessen den Geburtstagskuchen ausbläst , dann ist sein Atem das, was in der Strömungsmechanik als „laminare Strömung“ bezeichnet wird. Eine Strömung, deren Stromlinien parallel zueinander verlaufen und deren Einzelelemente, hier die Luftmoleküle, einen gleichgerichteten Weg verfolgen. „Turbulent“ ist hingegen der von den Kerzen aufsteigende Rauch über dem Kuchen, der „Richtung Osten“ hinweggeblasen wird. Diese Turbulenz von Flüssigkeiten und Gasen ist allgegenwärtig in unserer Welt, ob für unser Auge direkt sichtbar oder unsichtbar: die Luft, die tief in unsere Lunge eingesogen wird bis hin zu den Lungenbläschen, die Spuren der Milch beim Umrühren in einer Kaffeetasse, der Rauch, der von einem Schornstein aufsteigt, das Wasser rund um die großen und kleinen Kiesel in einem Bach, der Wind, der den Kirchturm und die Hausecke umwirbelt, das heiße Plasma, das in Feuerfackeln von der Sonnenoberfläche ins Weltall hinauslodert, oder die großen Wolken kosmischen Staubs, die sich in Strudeln und Wirbeln zu Galaxien oder Sternen verdichten. „Turbulent“ ist also eine Strömung, deren Stromlinien sich zu überkreuzen scheinen und deren Einzelelemente keinen gleichgerichteten Weg verfolgen. Stattdessen existieren vielfältigste Formen und Muster von miteinander verschränkten Wirbeln auf allen Größenskalen. Aufgrund dieser Komplexität in Formen und Skalen gehört die Beschreibung und Vorhersage von Turbulenz schon seit Jahrhunderten zu den großen Rätseln in Physik und Mathematik. Da turbulente Strömungen gleichwohl derart zentral sind für viele Bereiche menschlichen Lebens und Handelns, werden Grundlagen- und Anwendungsforschung mit Nachdruck vorangetrieben. Die vorliegende Arbeit umfaßt gleich drei Anwendungsfälle von Turbulenzforschung, und es darf als bezeichnend für die Allgegenwart der Turbulenz angesehen werden, daß sich diese drei Anwendungen in solch unterschiedlichen Größenskalen abspielen. Die Windabkühlung von Rindern, die bis in den Sub-Millimeter-Maßstab im Bereich von Fell und Hautoberfläche hinein betrachtet werden muß, die Sturmge-fährdung von Wäldern, für die Größen zwischen einem halben Meter an den Bäumen und mehreren Kilometern in der Landschaft relevant sind, und zu guter Letzt das turbulente Geschehen in kos-mischen Gaswolken und Galaxienhaufen, das sich im Größenbereich von vielen Millionen Lichtjahren abspielt. Nicht nur in den Techniken der Modellierung, sondern auch in der physikalischen Wirklich-keit sind diese Phänomene trotz der gewaltigen Größenunterschiede eng verwandt. In diesem Sinne: vom Kosmos zur Kuh.

There is currently a great deal of interest within the aviation community in the design of small, slow-flying but manoeuvrable uninhabited vehicles for reconnaissance, surveillance, and search and rescue operations in urban environments. Inspired by observation of birds, insects, fish and cetaceans, flapping wings are being actively studied in the hope that they may provide greater propulsive efficiencies than propellers and rotors at low Reynolds numbers for such Micro-Air Vehicles (MAVs). Researchers have posited the Strouhal number (combining flapping frequency, amplitude and forward speed) as the parameter controlling flapping wing aerodynamics in cruising flight, although there is conflicting evidence. This thesis explores the effect of flapping frequency and amplitude on forces and wake structures, as well as physical mechanisms leading to optimum propulsive efficiency. Two-dimensional rigid airfoils are considered at Reynolds number 2,000 ??? 40,000. A compressible Navier-Stokes simulation is combined with numerical and analytical potential flow techniques to isolate and evaluate the effect of viscosity, leading and trailing edge vortex separation, and wake vortex dynamics. The wake structures of a plunging airfoil are shown to be sensitive to the flapping frequency independent of the Strouhal number. For a given frequency, the wake of the airfoil exhibits ???vortex lock-in??? as the amplitude of motion is increased, in a manner analogous to an oscillating circular cylinder. This is caused by interaction between the flapping frequency and the ???bluff-body??? vortex shedding frequency apparent even for streamlined airfoils at low Reynolds number. The thrust and propulsive efficiency of a plunging airfoil are also shown to be sensitive to the flapping frequency independent of Strouhal number. This dependence is the result of vortex shedding from the leading edge, and an interaction between the flapping frequency and the time for vortex formation, separation and convection over the airfoil surface. The observed propulsive efficiency peak for a pitching and plunging airfoil is shown to be the result of leading edge vortex shedding at low flapping frequencies (low Strouhal numbers), and high power requirements at large flapping amplitudes (high Strouhal numbers). The efficiency peak is governed by flapping frequency and amplitude separately, rather than the Strouhal number directly.

The production and transportation of hydrocarbon fluids in multiphase pipelines could be severely hindered by particulate solids deposit such as produced sand particles which accompany hydrocarbon production. Knowledge of the flow characteristics of solid particles in fluids transported in pipelines is important in order to accurately predict solid particles deposition in pipelines. This research thesis presents the development of a three-dimensional (3D) Computational Fluids Dynamics (CFD) modelling technique for the prediction of liquid-solids multiphase flow in pipes, with special emphasis on the flow in V-inclined pipe bends. The Euler-Euler (two-fluid) multiphase modelling methodology has been adopted and the multiphase model equations and closure models describing the liquid-solids flow have been implemented and calculated using the finite volume method in a CFD code software. The liquid phase turbulence has been modelled using a two-equation k−ε turbulence model which contains additional terms to account for the effects of the solid-particles phase on the multiphase turbulence structure. The developed CFD numerical framework has been verified for the relevant forces and all the possible interaction mechanisms of the liquid-solids multiphase flow by investigating four different numerical frameworks, in order to determine the optimum numerical framework that captures the underlying physics and covers the interaction mechanisms that lead to sand deposition and the range of sand transport flow regimes in pipes. The flow of liquid-sand in pipe has been studied extensively and the numerical results of sand concentration distribution across pipe and other flow properties are in good agreement with published experimental data on validation. The numerical framework has been employed to investigate the multiphase flow in V-inclined pipe bends of ±4o−6o, seemingly small inclined bend angles. The predicted results which include the sand segregation, deposition velocity and flow turbulence modulation in the pipe bend show that the seemingly small pipe bends have significant effect on the flow differently from that of horizontal pipes. The pipe bend causes abrupt local change in the multiphase flow characteristic and formation of stationary sand deposit in the pipe at a relatively high flow velocity. The threshold velocity to keep sand entrained in liquid in pipe bends is significantly higher than that required for flow horizontal pipes. A critical implication of this is that the correlations for predicting sand deposition in pipelines must account for the effect of pipe bend on flow characteristics in order to provide accurate predictions of the critical sand transport velocity (MTV) in subsea petroleum flowlines, which V-inclined pipe bends are inevitable due to seabed topology.

Nesta tese é apresentado um esquema novo de alta resolução upwind (denominado TDPUS-C3) para reconstrução de fluxos numéricos para leis de conservação não lineares e problemas relacionados em DFC. O esquema é baseado nos critérios de estabilidade CBC e TVD e desenvolvido utilizando condições de diferenciabilidade \'C POT. 3\'. Além disso, é realiozada a implementação da associação do esquema TDPLUS-C3 com a modelagem de turbulência RNG \'\\kappa - \\epsilon\'. O propósito é obter soluções numéricas de sistemas hiperbólicos de leis de conservação para dinâmica dos gases e equações de Navier-Stokes para escoamento incompreensível de fluidos newtonianos e não newtonianos (viscoelásticos). Fazendo o uso do esquema TDPUS-C3, a precisão global dos métodos numéricos é verificada acessando o erro em problemas teste (benchmark) 1D e 2D. Um estudo comparativo entre os resultados do esquema TDPUS-C3 e os esquemas upwind convencionais para leis de conservação hiperbólicas complexas é também realizado. A Associação das modelagens numéricas (upwinding mais RNG \'\\kappa - \\epsilon\') é , então, examinada na simulação de escoamentos turbulentos de fluidos newtonianos envolvendo superfícies livres móveis, usando a metodologia URANS. No geral, em termos do comportamento global, concordância satisfatória é observada / In this thesis, a new high-resolution upwind scheme (named TDPUS-C3) for reconstruction of numerical fluxes for nonlinear conservation laws and related CFD problems in presented. The scheme is based on CBC and TVD stability criteria and developed by employing differentiability condictions (\'C POT. 3\'). In additon, the implementation of an association of the TDPUS-C3 scheme with the RNG \'\\kappa - \\epsilon\' turbulence modelling is also performed. The purpose is to obtain numerical solutions of systems of hyperbolic conservation laws for gas dynamics and Navier-Stokes equations for incompressible flow of Newtonian and non-Newtonian (viscoelstic) fluids. By using the TDPUS-C3 scheme, the global accuracy of the numerical methods is verified by assessing the error on 1D and 2D benchmark test cases. A comparative study between the TDPUS-C3 scheme and convectional upwind schemes to solve standard and complex hyperbolic conservation laws is also accomplished. The association of the numerical modelling (upwinding plus RNG \'\\kappa - epsilon\') is then examined in the simulation of turbulent Newtonian fluid flows involving moving free surfaces, by using URANS methodology. Overall, satisfactory agreement is found in terms of the overall behaviour

The turbine inlet guide vanes of a gas-turbine are subjected to extreme hot gas temperatures which increases the risk of mechanical failure and overall reduces the component lifespan. Hence, it is of great interest for gas-turbine manufacturers to establish methods for accurately estimating the temperature distribution along the vane surface. Due to the three-dimensional nature of turbine flow, it is of interest to establish Computational Fluid Dynamics (CFD) methodology which capture these three-dimensional effects. This thesis is one in a collection of theses conducted at Siemens Energy AB on the subject. Previous studies have investigated and validated the implementation of RANS simulations on non-cooled turbine vanes and endwalls. In this study, the focus is on studying a film cooled vane and establishing one RANS as well as one hybrid modelling strategy for heat transfer coefficient (HTC) predictions. The HTC prediction capabilities are compared and validated against experimental data presented in the doctoral thesis by Fabrizio Fontaneto on the LS94 vane at Von Karman Institute. The chosen RANS modelling method was the Shear Stress Transport (SST) k-ω turbulence model, with γ-Reθ transition modelling, based on the findings by Enico (2021) and Daugulis (2022). The model proved capable in estimating the HTC well on mainly the suction side of the vane. The pressure side HTC was largely under-predicted, a common issue with the SST model also seen in the previous theses as well as the hybrid simulations. The strength of the SST k-ω turbulence model, with γ-Reθ transition modelling, is in accurately capturing the HTC magnitude, most likely due to the well-predicted turbulence intensity decay at the inlet. However, it largely under-predicts the HTC along the suction side film-coolant layer, implying that it may be over-estimating the film-cooling capabilities. The hybrid model chosen was the Scale Resolving Hybrid (SRH) model, with underlying RANS SST k-ω. Compared to RANS, hybrid results were under-estimated, seemingly offset from the experimental data by a constant 200 units along the entire vane midspan. This is likely due to the inaccurate turbulence intensity presented in the SRH simulations, which decays quickly along the inlet compared to RANS and experimental data. Yet still, the hybrid model showed potential in capturing certain results not seen with RANS, such as the secondary flow effects by the vane endwalls, as well as arguably capturing the general HTC trend at midspan seen in the experimental data. Additionally, the section of severely under-predicted HTC by the suction side film-coolant seen with RANS is not present in the hybrid results. Although the hybrid model has proven promising in many aspects, in its current state it is not a viable method for HTC predictions due to its general under-prediction of HTC. Largely, the authors suspect this is due to the undesirably coarse mesh around the cooling holes, which leads to RANS computation in regions where SRH is desired. Thus, improvements would need to be made to the model, where, for example, implementing a zonal hybrid RANS-LES model would be an option. Considering the hybrid model in its current state, RANS is the preferred method, especially when considering the greater computational cost and the labor associated with hybrid simulations which were experienced during this study. In conclusion, it is evident that the correct capture of inlet turbulence intensity decay as well as suitable mesh refinement by the cooling holes are crucial for obtaining the correct magnitudes of HTC, and thus, the capture of it should be of utmost priority in future work within the field.

L'objectif de la présente étude est d'analyser la modélisation de la turbulence de simulations en moyenne de Reynolds (RANS) dans le cadre d'écoulements de type turbomachines, en utilisant des simulations aux grandes échelles (SGE) comme référence. L'étude porte sur deux cas test: un décollement de coin dans une grille d'aubes rectiligne, et un écoulement de jeu pour un aubage isolé dans un jet. Deux lois de comportement, la loi de comportement de Boussinesq et la loi de comportement quadratique (quadratic constitutive relation ou QCR), sont analysées, avec deux versions du modèle de turbulence k-omega de Wilcox. Les lois de comportement étudiées reposent sur deux hypothèses: une hypothèse d'alignement entre le tenseur de Reynolds et un tenseur construit à partir de l'écoulement moyen, et une hypothèse sur la viscosité turbulente. L'hypothèse d'alignement est étudiée à partir de la SGE, pour laquelle les deux tenseurs sont indépendamment connus, en utilisant un indicateur construit sur le produit scalaire des tenseurs. Les résultats sont présentés sous forme d'une fonction de répartition de la valeur de l'indicateur pour le domaine complet, puis pour trois sous-domaines d'intérêt: l'entrée, une région où l'écoulement interagit fortement avec les parois, et une région où l'écoulement est fortement tourbillonnaire. L'hypothèse d'alignement n'est que rarement valide pour la loi de comportement de Boussinesq. Pour la QCR, les résultats sont meilleurs en entrée, comparé à la loi de Boussinesq. Il ne sont cependant pas meilleurs pour les régions où l'écoulement est plus tourbillonnaire. Une amélioration de la loi de comportement est nécessaire pour pouvoir faire progresser la modélisation turbulente en RANS. En revanche, l'utilisation de l'énergie cinétique turbulente et du taux de dissipation spécifique semble correcte pour estimer la valeur de la viscosité turbulente. L'analyse de la modélisation de l'équation d'énergie cinétique turbulente (ECT) est réalisée au travers d'une comparaison terme à terme avec l'équation d'ECT résolue par la SGE. Les résultats SGE présentent une turbulence qui n'est pas à l'équilibre : la production et la dissipation ne sont pas superposées, et le terme de transport est important. Pour le RANS, la turbulence est à l'équilibre : la production et la dissipation sont superposées, et le terme de transport est de faible intensité. Un modèle de turbulence qui prend en compte le déséquilibre est nécessaire pour améliorer ce point. En dernier lieu, une nouvelle formulation hybride RANS/SGE est proposée, fondée sur la distance à la paroi en unités de paroi. La formulation est validée dans un canal bi-périodique et un premier essai est réalisé sur le cas de décollement de coin, mais d'autres analyses sont nécessaires avant que cette formulation ne soit fonctionnelle. / The present study aims at analysing turbulence modelling in Reynolds-averaged Navier-Stokes (RANS) simulations, in the context of turbomachinery flows, using large-eddy simulations (LES) as references. Two test cases are considered: a corner separation (CS) flow in a linear compressor cascade, and a tip-leakage (TL) flow of a single blade in a jet. Two constitutive relations, the Boussinesq constitutive relation and the quadratic constitutive relation (QCR), are investigated, with two versions of Wilcox's $k-\omega$ turbulence model. The studied constitutive relations rely on two hypotheses: an alignment hypothesis between the Reynolds stress tensor and a mean flow tensor, and an hypothesis on the turbulent viscosity. The alignment hypothesis is investigated using LES, where both the tensors are known independently, with an indicator built on the inner product of the tensors. The results are presented as probability density functions of the indicator value for the entire domain first, and then for three specific areas of interest: the inlet area, similar to a boundary-layer flow, an area of strong interaction between the flow and the walls (CS: passage area, TL: tip clearance) and an area of highly vortical flow (CS: separation wake, TL: tip-leakage vortex). The alignment hypothesis is rarely verified in any area for the Boussinesq constitutive relation. For the QCR, the results are improved for the inlet areas compared to the Boussinesq constitutive relation, but no significant improvement is found in the highly vortical regions. An improvement of the constitutive relation is needed in order to improve the RANS turbulence modelling. In contrast, the use of the turbulent kinetic energy and the specific dissipation rate appears quite correct to estimate the turbulent viscosity. The modelling of the RANS turbulent kinetic energy (TKE) budget equation is investigated through a term to term comparison with the resolved LES TKE budget equation. The LES presents a turbulence that is not at equilibrium, with the production and the dissipation not superimposed, and an important amount of transport. This differs from the RANS models, at equilibrium: the production and the dissipation are superimposed, with a small amount of transport. The development of a non-equilibrium turbulence model for RANS simulations could improve this aspect of turbulence modelling. Finally, a new hybrid RANS-LES formulation, based on the wall distance in wall units, is also proposed. It is validated on a bi-periodical channel flow, and a first attempt is made on the corner separation case, but further investigations are still needed for the model to be fully operational.

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