This work investigates through numerical simulation, the operation of a big bore direct fuel injection spark ignition engine, at part load under stratified operation. It evaluates fuel-air mixture preparation and combustion process with the adoption of detailed chemical kinetics mechanisms for both Isooctane and Ethanol, and applying adaptive mesh refinement to capture the turbulent flame brush. The investigation is split in 3 main parts. In the first part, with Isooctane as fuel, the impact of in-cylinder turbulence level induced by squish has shown that the attempt to isolate the squish ratio, maintaining the bowl shape, for the evaluated cases have led to a scenario not more appropriate for flame initiation and propagation for 2 of the 3 geometries. But the observations made during this initial stage have led to the proposal of a fourth geometry to improve the mixture formation and combustion process. As it was seen the combustion process was about 11.5 deg faster with the new piston bowl proposed. In the second part, still with Isooctane and maintaining the new proposed piston, evaluates the influence of two types of hollow cone fuel injectors, an inwardly and an outwardly opening types, where maintain fixed spark timing, the end of injection is varied and compared among the two cases, while targeting for the same gross IMEP output. The main results are that the outwardly opening injector case resulted in better fuel-air mixture preparation, even with a late end of injection. This led to higher combustion efficiency and lower unburned hydrocarbon, CO and soot emissions, while increasing NOx emissions. The 10-90% MFB burn duration is higher for the outwardly opening injector case. In the last part the outwardly opening spray injector from the previous part, but using Ethanol as fuel has shown that to attain the same IMEP level the injected fuel mass is increased with Ethanol, and with its higher latent heat of vaporization, the time required to have an ignitable fuel-air mixture more than doubled that for the Isooctane case. Another important effect of these is the excessive increase of THC emissions. The overall combustion duration was faster for the Ethanol, mainly as the last part of the combustion was almost twice as fast as for the Isooctane case. It may be a consequence of a more homogeneous fuel-air mixture cloud as the fuel has more time to diffuse as the EOI is more advanced. The in-cylinder charge cooling effect of Ethanol led to a reduction in the in-cylinder temperature, leading to a reduction in NOx formation. CO emissions was also lowered, which is maybe attributed to either the reduced chemical dissociation with the lower temperatures, or reduced fuel rich regions. The reduced fuel rich regions also explain the reason for lower soot emissions.

O escoamento estratificado óleo-água é bastante comum na indústria do petróleo, especialmente em poços direcionas offshore, oleodutos e gasodutos. Entretanto, existem poucos trabalhos na literatura sobre a natureza da estrutura ondulatória observada no escoamento em dutos ascendentes ou descendentes a partir da horizontal. O objetivo deste trabalho foi estudar as propriedades geométricas e cinemáticas da onda interfacial, i.e, forma média, comprimento, amplitude e celeridade, e assim contribuir para a compreensão do papel da onda interfacial na dinâmica do escoamento estratificado. Um software baseado em plataforma Labview® possibilitou a automação para obtenção dos dados das ondas interfaciais extraídas de imagens de vídeos de alta resolução. Além das propriedades das ondas, também se coletaram valores de fração volumétrica in situ e de gradiente de pressão bifásico para cinco ângulos de inclinação (-20°, -10°, 0°, 10° e 20°) em diferentes pares de vazões de óleo e água. Desenvolveu-se um modelo fenomenológico considerando os termos ondulatórios do escoamento para o cálculo da fração volumétrica in situ e do gradiente de pressão bifásico, sendo comparado com modelos disponíveis na literatura e dados experimentais. A concordância do modelo proposto com os dados coletados neste trabalho se mostrou muito boa, o que sugere um avanço em comparação ao existente na literatura. O estudo da equação da onda de perturbação interfacial para o escoamento estratificado óleo-água indicou que a natureza da onda observada é cinemática e não dinâmica; e baseado na equação da celeridade da onda cinemática pode-se confrontar a celeridade experimental com a teórica, revelando boa concordância. / The oil-water stratified flow is quite common in the oil industry, especially in offshore directional wells and pipelines. However, there are few studies on the physics of the wavy structure observed in upward and downward stratified flow. The goal of this work was to study the geometric and kinematic properties of interfacial waves, i.e., the average shape, wavelength, amplitude and celerity. A homemade Labview®-based software enabled the automatic acquisition of data extracted from frames obtained via high resolution video recording. In situ volume fraction and two-phase pressure gradient data for five inclination angles (-20°, -10°, 0°, 10° and 20°) at several pairs of oil and water flow rates were also collected. A phenomenological model that takes into account the wavy structure is proposed to calculate volume fractions and two-phase pressure gradient and it was compared with available models from the literature and experimental data. The good agreement of the proposed model with the data collected in this study is promising and suggests that it may provide better predictions in comparison with models from the literature. The study of the interfacial perturbation wave equation for stratified flow indicates that the observed waves nature is kinematic and not dynamic; and based on the kinematic wave velocity equation we could compare the experimental celerity with the theoretical one, with good agreement.

Embora o escoamento estratificado óleo-água seja comum na indústria do petróleo, existem poucos trabalhos na literatura sobre a estrutura ondulatória encontrada no escoamento em dutos ascendentes ou descendentes. O objetivo deste trabalho é entender e caracterizar as estruturas interfaciais ondulatórias em escoamentos estratificados óleo e água em dutos levemente inclinados e comparar com os resultados retirados da modelagem proposta. Uma primeira tentativa de modelar a onda interfacial óleo-água através da equação da energia bifásica unidimensional para regime permanente é apresentada e comparada aos resultados obtidos experimentalmente. Valores de comprimento de onda foram coletados para quatro ângulos de inclinação a partir da horizontal (-5º, -2º, 2º e 5º) e diversos pares de vazão de óleo e água. Os dados foram extraídos de frames de vídeos de alta resolução através de uma técnica manual. Um software baseado em LabView foi desenvolvido para possibilitar a automação da obtenção do comprimento de onda e amplitude. Os comprimentos de onda obtidos com a técnica automática foram comparados com os dados obtidos através da técnica manual. O mensuramento é bem satisfatório e sugere que a ferramenta experimental proposta possa ser aplicada para o estudo de qualquer padrão de escoamento óleo-água, onde uma estrutura ondulatória possa ser identificada. A validade da modelagem da onda interfacial proposta foi avaliada através da comparação entre os resultados teóricos e dados experimentais. A concordância observada é encorajadora. / Even though the oil-water stratified flow pattern has a common occurrence in the upstream oil industry, quite often in directional wells, there are only a few works in the literature dealing with its interfacial wavy structure. This work has the mean goal of comprehending a characterizing the interfacial wave structure in inclined water-oil stratified pipe flow, comparing it with the proposed model. A model, based on the steady-state one-dimensional two-phase flow energy equation is presented and compared to our experimental data. Wave lengths and amplitude data were collected at four inclinations from the horizontal (-5º, -2º, 2º and 5º) and several oil and water flow rates. The data were extracted form high-resolution video images through a manual technique. A LabView based software was developed in order to obtain the lengths and amplitude automatically. Wave length obtained from it were compared to the manual technique data. The satisfactory measurement suggests the applicability of the experimental tool to any water-oil pipe flow pattern, where a wave structure can be identified. The validity by the comparison between theoretical and experimental data. The agreement observed is promising.

Embora o escoamento estratificado óleo-água seja comum na indústria do petróleo, existem poucos trabalhos na literatura sobre a estrutura ondulatória encontrada no escoamento em dutos ascendentes ou descendentes. O objetivo deste trabalho é entender e caracterizar as estruturas interfaciais ondulatórias em escoamentos estratificados óleo e água em dutos levemente inclinados e comparar com os resultados retirados da modelagem proposta. Uma primeira tentativa de modelar a onda interfacial óleo-água através da equação da energia bifásica unidimensional para regime permanente é apresentada e comparada aos resultados obtidos experimentalmente. Valores de comprimento de onda foram coletados para quatro ângulos de inclinação a partir da horizontal (-5º, -2º, 2º e 5º) e diversos pares de vazão de óleo e água. Os dados foram extraídos de frames de vídeos de alta resolução através de uma técnica manual. Um software baseado em LabView foi desenvolvido para possibilitar a automação da obtenção do comprimento de onda e amplitude. Os comprimentos de onda obtidos com a técnica automática foram comparados com os dados obtidos através da técnica manual. O mensuramento é bem satisfatório e sugere que a ferramenta experimental proposta possa ser aplicada para o estudo de qualquer padrão de escoamento óleo-água, onde uma estrutura ondulatória possa ser identificada. A validade da modelagem da onda interfacial proposta foi avaliada através da comparação entre os resultados teóricos e dados experimentais. A concordância observada é encorajadora. / Even though the oil-water stratified flow pattern has a common occurrence in the upstream oil industry, quite often in directional wells, there are only a few works in the literature dealing with its interfacial wavy structure. This work has the mean goal of comprehending a characterizing the interfacial wave structure in inclined water-oil stratified pipe flow, comparing it with the proposed model. A model, based on the steady-state one-dimensional two-phase flow energy equation is presented and compared to our experimental data. Wave lengths and amplitude data were collected at four inclinations from the horizontal (-5º, -2º, 2º and 5º) and several oil and water flow rates. The data were extracted form high-resolution video images through a manual technique. A LabView based software was developed in order to obtain the lengths and amplitude automatically. Wave length obtained from it were compared to the manual technique data. The satisfactory measurement suggests the applicability of the experimental tool to any water-oil pipe flow pattern, where a wave structure can be identified. The validity by the comparison between theoretical and experimental data. The agreement observed is promising.

O escoamento estratificado óleo-água é bastante comum na indústria do petróleo, especialmente em poços direcionas offshore, oleodutos e gasodutos. Entretanto, existem poucos trabalhos na literatura sobre a natureza da estrutura ondulatória observada no escoamento em dutos ascendentes ou descendentes a partir da horizontal. O objetivo deste trabalho foi estudar as propriedades geométricas e cinemáticas da onda interfacial, i.e, forma média, comprimento, amplitude e celeridade, e assim contribuir para a compreensão do papel da onda interfacial na dinâmica do escoamento estratificado. Um software baseado em plataforma Labview® possibilitou a automação para obtenção dos dados das ondas interfaciais extraídas de imagens de vídeos de alta resolução. Além das propriedades das ondas, também se coletaram valores de fração volumétrica in situ e de gradiente de pressão bifásico para cinco ângulos de inclinação (-20°, -10°, 0°, 10° e 20°) em diferentes pares de vazões de óleo e água. Desenvolveu-se um modelo fenomenológico considerando os termos ondulatórios do escoamento para o cálculo da fração volumétrica in situ e do gradiente de pressão bifásico, sendo comparado com modelos disponíveis na literatura e dados experimentais. A concordância do modelo proposto com os dados coletados neste trabalho se mostrou muito boa, o que sugere um avanço em comparação ao existente na literatura. O estudo da equação da onda de perturbação interfacial para o escoamento estratificado óleo-água indicou que a natureza da onda observada é cinemática e não dinâmica; e baseado na equação da celeridade da onda cinemática pode-se confrontar a celeridade experimental com a teórica, revelando boa concordância. / The oil-water stratified flow is quite common in the oil industry, especially in offshore directional wells and pipelines. However, there are few studies on the physics of the wavy structure observed in upward and downward stratified flow. The goal of this work was to study the geometric and kinematic properties of interfacial waves, i.e., the average shape, wavelength, amplitude and celerity. A homemade Labview®-based software enabled the automatic acquisition of data extracted from frames obtained via high resolution video recording. In situ volume fraction and two-phase pressure gradient data for five inclination angles (-20°, -10°, 0°, 10° and 20°) at several pairs of oil and water flow rates were also collected. A phenomenological model that takes into account the wavy structure is proposed to calculate volume fractions and two-phase pressure gradient and it was compared with available models from the literature and experimental data. The good agreement of the proposed model with the data collected in this study is promising and suggests that it may provide better predictions in comparison with models from the literature. The study of the interfacial perturbation wave equation for stratified flow indicates that the observed waves nature is kinematic and not dynamic; and based on the kinematic wave velocity equation we could compare the experimental celerity with the theoretical one, with good agreement.

Escoamentos bifásicos estão presentes em diversos processos naturais e industriais, como na indústria de petróleo. Podem apresentar-se em diferentes configurações topológicas, ou, padrões de escoamento, entre eles o escoamento estratificado ondulado e o estratificado com mistura na interface. Os escoamentos bifásicos estratificados óleo-água têm sido utilizados como uma forma conveniente de evitar a formação de emulsões de água em óleo em oleodutos e possuem uma ocorrência comum em poços de petróleo direcionais. Quando a onda interfacial ultrapassa determinado limite geométrico e cinemático, surge o fenômeno do entranhamento de gotas, representado por misturas entre as fases junto à interface que promovem um aumento na queda de pressão. Modelos têm sido apresentados pela literatura na tentativa de descrever o fenômeno do entranhamento de gotas. Neste trabalho é apresentada uma nova proposta de modelagem matemática unidimensional para o entranhamento de gotas com o objetivo de melhorar a previsão dos parâmetros envolvidos, em especial, da fração volumétrica de óleo e da queda de pressão bifásica. Também foi utilizada simulação numérica computacional, CFD (Computational Fluid Dynamics), com o uso de software comercial para obtenção dos valores dos parâmetros do escoamento estratificado ondulado óleo-água (fração volumétrica de óleo, queda de pressão, amplitude e comprimento da onda interfacial). Os resultados da modelagem fenomenológica para entranhamento e os de CFD foram comparados com bancos de dados experimentais. Os resultados em CFD mostram concordância com os resultados experimentais, tanto na análise qualitativa das propriedades geométricas das ondas interfaciais, quanto na comparação direta com os dados para fração volumétrica e queda de pressão. Os resultados numéricos da modelagem fenomenológica para fatores de entranhamento apresentam boa concordância com dados da literatura. / Two-phase flows are present in many natural and industrial processes, such as in the oil industry. They may be found arranged in several flow patterns, including the wavy stratified flow and the stratified with mixing at the interface. The stratified oil-water flow has been used as a convenient way to avoid the formation of emulsions of water in oil and have a common occurrence in directional oil wells. When the interfacial wave exceeds a certain geometric and kinematic limit the phenomenon of drop entrainment arises at the interface, causing an increase of pressure drop. Models have been presented in the literature in an attempt to describe the phenomenon of drop entrainment. In this work, we present a new method for the one-dimensional mathematical modeling of entrainment in order to improve the prediction of oil volume fraction and pressure drop. It was also used a commercial computational fluid dynamics tool (CFD) to obtain the values of flow parameters of wavy stratified oil-water flow, such as oil volume fraction, pressure drop, amplitude and length of the interfacial waves. The results of the phenomenological modeling for entrainment and CFD were compared with experimental databases. The CFD results are in agreement with the experimental results in both the qualitative analysis of the geometric properties of the interfacial waves and in direct comparison with oil-volumetric-fraction and pressure-drop data. The numerical results of the phenomenological model for entrainment factors are in agreement with data from the literature.

Escoamentos bifásicos estão presentes em diversos processos naturais e industriais, como na indústria de petróleo. Podem apresentar-se em diferentes configurações topológicas, ou, padrões de escoamento, entre eles o escoamento estratificado ondulado e o estratificado com mistura na interface. Os escoamentos bifásicos estratificados óleo-água têm sido utilizados como uma forma conveniente de evitar a formação de emulsões de água em óleo em oleodutos e possuem uma ocorrência comum em poços de petróleo direcionais. Quando a onda interfacial ultrapassa determinado limite geométrico e cinemático, surge o fenômeno do entranhamento de gotas, representado por misturas entre as fases junto à interface que promovem um aumento na queda de pressão. Modelos têm sido apresentados pela literatura na tentativa de descrever o fenômeno do entranhamento de gotas. Neste trabalho é apresentada uma nova proposta de modelagem matemática unidimensional para o entranhamento de gotas com o objetivo de melhorar a previsão dos parâmetros envolvidos, em especial, da fração volumétrica de óleo e da queda de pressão bifásica. Também foi utilizada simulação numérica computacional, CFD (Computational Fluid Dynamics), com o uso de software comercial para obtenção dos valores dos parâmetros do escoamento estratificado ondulado óleo-água (fração volumétrica de óleo, queda de pressão, amplitude e comprimento da onda interfacial). Os resultados da modelagem fenomenológica para entranhamento e os de CFD foram comparados com bancos de dados experimentais. Os resultados em CFD mostram concordância com os resultados experimentais, tanto na análise qualitativa das propriedades geométricas das ondas interfaciais, quanto na comparação direta com os dados para fração volumétrica e queda de pressão. Os resultados numéricos da modelagem fenomenológica para fatores de entranhamento apresentam boa concordância com dados da literatura. / Two-phase flows are present in many natural and industrial processes, such as in the oil industry. They may be found arranged in several flow patterns, including the wavy stratified flow and the stratified with mixing at the interface. The stratified oil-water flow has been used as a convenient way to avoid the formation of emulsions of water in oil and have a common occurrence in directional oil wells. When the interfacial wave exceeds a certain geometric and kinematic limit the phenomenon of drop entrainment arises at the interface, causing an increase of pressure drop. Models have been presented in the literature in an attempt to describe the phenomenon of drop entrainment. In this work, we present a new method for the one-dimensional mathematical modeling of entrainment in order to improve the prediction of oil volume fraction and pressure drop. It was also used a commercial computational fluid dynamics tool (CFD) to obtain the values of flow parameters of wavy stratified oil-water flow, such as oil volume fraction, pressure drop, amplitude and length of the interfacial waves. The results of the phenomenological modeling for entrainment and CFD were compared with experimental databases. The CFD results are in agreement with the experimental results in both the qualitative analysis of the geometric properties of the interfacial waves and in direct comparison with oil-volumetric-fraction and pressure-drop data. The numerical results of the phenomenological model for entrainment factors are in agreement with data from the literature.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada. / To study and forecast the environmental impacts that a hydroelectric power plant reservoir may produce, a numerical simulation that takes into account the effects of stratification and scalar mixing (such as temperature, salinity or water-soluble substance) is required. This work proposes a methodology for the study of the environmental fluid flow phenomena, mainly for flows in which the knowledge of the interaction between stratification and mixing can give important notions of the phenomena that occur. For this, a numerical simulation tool for 3D environmental flow is developed. A tetrahedral mesh generator of the reservoir based on the terrain topology and an algebraic turbulence model based on the Richardson number are the main tools developed. The main difficulty in tetrahedral mesh generation of a reservoir is nonuniform distribution of the points related to the huge ratio between the horizontal and vertical scales of the reservoir. In this type of point distributions, conventional tetrahedron mesh generation algorithm may become unstable. For this reason, a unstructured tetrahedral mesh generator is developed and the methodology used to obtain conforming elements is described. Triangular surface mesh generation using the Delaunay triangulation and the construction of the tetrahedra from the triangular surface mesh are the main steps to the mesh generator. The hydrodynamic simulation of reservoirs with a turbulence model provides a useful tool that is computationally viable for engineering purposes. Furthermore, the turbulence model based on the Richardson number takes into account the effects of interaction between turbulence and stratification. The algebraic model is the simplest among the various models of turbulence, but provides realistic results with the fitting of a small amount of parameters. Eddy-Viscosity/Diffusivity models for stratified turbulent flows models are incorporated. Using the Finite Element Method (FEM) approximation the Reynolds-averaged Navier-Stokes (RANS) and mean scalar transport equations are approximated. The convective terms are discretized employing the Semi-Lagrangian method, and the spatial discretization is based on the Galerkin method. The computational results are compared with the results available in the literature. Finally, the simulation of the flow in a branch of a reservoir is presented.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Uma simulação numérica que leva em conta os efeitos de estratificação e mistura escalar (como a temperatura, salinidade ou substância solúvel em água) é necessária para estudar e prever os impactos ambientais que um reservatório de usina hidrelétrica pode produzir. Este trabalho sugere uma metodologia para o estudo de escoamentos ambientais, principalmente aqueles em que o conhecimento da interação entre a estratificação e mistura pode dar noções importantes dos fenômenos que ocorrem. Por esta razão, ferramentas de simulação numérica 3D de escoamento ambiental são desenvolvidas. Um gerador de malha de tetraedros do reservatório e o modelo de turbulência algébrico baseado no número de Richardson são as principais ferramentas desenvolvidas. A principal dificuldade na geração de uma malha de tetraedros de um reservatório é a distribuição não uniforme dos pontos relacionada com a relação desproporcional entre as escalas horizontais e verticais do reservatório. Neste tipo de distribuição de pontos, o algoritmo convencional de geração de malha de tetraedros pode tornar-se instável. Por esta razão, um gerador de malha não estruturada de tetraedros é desenvolvido e a metodologia utilizada para obter elementos conformes é descrita. A geração de malha superficial de triângulos utilizando a triangulação Delaunay e a construção do tetraedros a partir da malha triangular são os principais passos para o gerador de malha. A simulação hidrodinâmica com o modelo de turbulência fornece uma ferramenta útil e computacionalmente viável para fins de engenharia. Além disso, o modelo de turbulência baseado no número de Richardson leva em conta os efeitos da interação entre turbulência e estratificação. O modelo algébrico é o mais simples entre os diversos modelos de turbulência. Mas, fornece resultados realistas com o ajuste de uma pequena quantidade de parâmetros. São incorporados os modelos de viscosidade/difusividade turbulenta para escoamento estratificado. Na aproximação das equações médias de Reynolds e transporte de escalar é utilizando o Método dos Elementos Finitos. Os termos convectivos são aproximados utilizando o método semi-Lagrangeano, e a aproximação espacial é baseada no método de Galerkin. Os resultados computacionais são comparados com os resultados disponíveis na literatura. E, finalmente, a simulação de escoamento em um braço de reservatório é apresentada. / To study and forecast the environmental impacts that a hydroelectric power plant reservoir may produce, a numerical simulation that takes into account the effects of stratification and scalar mixing (such as temperature, salinity or water-soluble substance) is required. This work proposes a methodology for the study of the environmental fluid flow phenomena, mainly for flows in which the knowledge of the interaction between stratification and mixing can give important notions of the phenomena that occur. For this, a numerical simulation tool for 3D environmental flow is developed. A tetrahedral mesh generator of the reservoir based on the terrain topology and an algebraic turbulence model based on the Richardson number are the main tools developed. The main difficulty in tetrahedral mesh generation of a reservoir is nonuniform distribution of the points related to the huge ratio between the horizontal and vertical scales of the reservoir. In this type of point distributions, conventional tetrahedron mesh generation algorithm may become unstable. For this reason, a unstructured tetrahedral mesh generator is developed and the methodology used to obtain conforming elements is described. Triangular surface mesh generation using the Delaunay triangulation and the construction of the tetrahedra from the triangular surface mesh are the main steps to the mesh generator. The hydrodynamic simulation of reservoirs with a turbulence model provides a useful tool that is computationally viable for engineering purposes. Furthermore, the turbulence model based on the Richardson number takes into account the effects of interaction between turbulence and stratification. The algebraic model is the simplest among the various models of turbulence, but provides realistic results with the fitting of a small amount of parameters. Eddy-Viscosity/Diffusivity models for stratified turbulent flows models are incorporated. Using the Finite Element Method (FEM) approximation the Reynolds-averaged Navier-Stokes (RANS) and mean scalar transport equations are approximated. The convective terms are discretized employing the Semi-Lagrangian method, and the spatial discretization is based on the Galerkin method. The computational results are compared with the results available in the literature. Finally, the simulation of the flow in a branch of a reservoir is presented.