An experimental research was performed to study the effects of surface roughness on the characteristics of separated and reattached turbulent flows in an open channel. A backward facing step was used to induce flow separation. The rough surfaces comprised wire mesh grit-80 and sand grains of average diameter 1.5 mm. In each experiment, the Reynolds number based on the step height and freestream velocity of approach flow was fixed at 3240 and the Reynolds number based on the approach flow depth and freestream velocity was kept constant at 25130. Particle image velocimetry (PIV) technique was used to measure the flow velocity. The results showed that roughness effects on the mean and turbulent quantities are evident only in the recovery region. Moreover, roughness effects on the flow dynamics are dependent on the specific roughness element.

An experimental research was performed to study the effects of surface roughness on the characteristics of separated and reattached turbulent flows in an open channel. A backward facing step was used to induce flow separation. The rough surfaces comprised wire mesh grit-80 and sand grains of average diameter 1.5 mm. In each experiment, the Reynolds number based on the step height and freestream velocity of approach flow was fixed at 3240 and the Reynolds number based on the approach flow depth and freestream velocity was kept constant at 25130. Particle image velocimetry (PIV) technique was used to measure the flow velocity. The results showed that roughness effects on the mean and turbulent quantities are evident only in the recovery region. Moreover, roughness effects on the flow dynamics are dependent on the specific roughness element.

Ce travail a pour but de comparer les approches de simulation numérique des grandes échelles explicite (LES) et implicite (ILES) pour un jet turbulent non-réactif ou réactif d’hydrogène à grande vitesse dans un co-courant d’air, typique d’un super-statoréacteur. La résolution des calculs va de 32 × 32 × 128 à 256 × 256 × 1024, à l’aide d’un schéma WENO d’ordre 5. Les LES explicites emploient les modèles sous-maille de Smagorinsky et de Fonction de Structure Sélective, associés au transport moléculaire. Les LES implicites sont réalisées avec et sans modèle de transport moléculaire, en résolvant les équations de Navier- Stokes ou d’Euler. Dans le cas non-réactif, le modèle de Smagorinsky est trop dissipatif. Le modèle de Fonction de Structure Sélective améliore les résultats, sans faire mieux que l’approche ILES quelle que soit la résolution. Dans le cas réactif, une coupure physique visqueuse est indispensable pour fixer une épaisseur à la flamme, et assurer la convergence en maillage de l’approche ILES. On montre aussi que les résultats LES/ILES sont moins sensibles aux conditions d’injection que ceux de l’approche RANS. Le premier chapitre est une introduction générale au contexte de l’étude. Au second chapitre, on rappelle les équations générales pour un écoulement réactif et on détaille les modèles thermodynamique et de transport retenus. Au troisième chapitre, les équations de la LES et les modèles sous-maille sont présentés. On examine également quelques propriétés du schéma numérique. Le chapitre 4 est consacré à la méthode numérique et au code de calcul. Enfin, on présente les cas-tests et on discute les résultats au chapitre 5. / This work is intended to compare Large Eddy Simulation and Implicit Large Eddy Simulation (LES and ILES) for a turbulent, non-reacting or reacting high speed H2 jet in co-flowing air, typical of scramjet engines. Numerical simulations are performed at resolutions ranging from 32 × 32 × 128 to 256 × 256 × 1024, using a 5th order WENO scheme. Physical LES are carried out with the Smagorinsky and the Selective Structure Function models associated to molecular diffusion. Implicit LES are performed with and without molecular diffusion, by solving either the Navier-Stokes or the Euler equations. In the nonreacting case, the Smagorinsky model is too dissipative. The Selective Structure Function leads to better results, but does not show any superiority compared to ILES, whatever the grid resolution. In the reacting case, a molecular viscous cut-off in the simulation is mandatory to set a physical width for the reaction zone in the ILES approach, hence to achieve grid-convergence. It is also found that LES/ILES are less sensitive to the inlet conditions than the RANS approach. The first chapter is an introduction to the context of this study. In the second chapter, the governing equations for multispecies reacting flows are presented, with emphasis on the thermodynamic and transport models. In the third chapter, physical LES equations and explicit sub-grid modeling strategies are detailed. Some properties of the numerical scheme are also investigated. In chapter four, the numerical scheme and some aspects of the solver are explained. Finally, non-reacting and reacting numerical experiments are presented and the results are discussed.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The rice paddy is responsible for 12-16% of the anthropogenic emission of CH4 in the world. Brazil is one of most important rice paddy producer. The goal of this study is to understand the dynamics of methane exchange between rice paddy and the atmosphere in the south of Brazil. Two estimation techniques are used: eddy covariance and soil chamber. Data were obtained in a rice paddy farm in Cachoeira do Sul municipality. The total emission of methane in 2015-2016 rice paddy growing season (November 14, 2015 to April 16 2016) was 35.7 gCH4m−2. The solar radiation is the weather forcing that showed a higher correlation with methane emissions. / A produção de arroz irrigado por inundação é responsável por 12-26% de toda a emissão antropogênica de CH4 do mundo. O Brasil é um dos principais produtores de arroz irrigado por inundação. O objetivo deste estudo é compreender a dinâmica das trocas de metano entre o arroz irrigado por alagamento e a atmosfera no sul do Brasil. São utilizadas duas técnicas de estimativa de fluxos superficiais: Eddy Covariance e câmara de solo. Os dados foram obtidos sobre uma lavoura de arroz irrigado por inundação no município de Cachoeira do Sul. A emissão total de metano no arroz irrigado na safra de 2015-2016 (14 de novembro de 2015 a 16 de abril de 2016) foi de 35,7 gCH4m��2. A radiação solar é o forçante meteorológico que mostrou maior correlação com o fluxo de metano.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / This thesis quantifies the net change of CO2 (NEE) in an irrigated rice crop cultivation in RS (Brazil). The research is part of the project CT-HIDRO, which presents as main purpose the description of surface conditions for different ecosystems in the country. Data were gathered during the cultivation period, from October 2003 to March 2004. During this time, the average air temperature and the incident radiation were suitable for the development of the culture in such an area. The turbulent flows, calculated by the eddy covariance method, overestimated the NEE value, probably due to the low estimate of breathing tax, because, during some night-times, sensor registered little turbulence. Moreover, other factor could have interfered in the result of the flows, for instance, the time in which the flows are calculated. The Webb correction was applied in the curve of NEE, and a value smaller than 1% was found. Such result didn't modify the previous result significantly. The behavior of CO2 flow during the whole analyzed period showed a maximum of -0,47mgCO2 m−2 s −1 in January 2004, what demonstrates the largest absorption of CO2 from the atmosphere by the plant during the day and a smaller liberation at night, due to the development phase in which the plant is. In this case, the Webb correction was also little emphasized, not interfering in the obtained result. The curves that demonstrate the three phases of development of the rice presented a larger absorption of CO 2 during the reproductive phase. The behavior of the flow of CO 2 showed a negative flow during the day, when the plant absorbs CO2 from the atmosphere, and a positive flow during the night, when the inverse process occurs, what was waited for. The latent heat flow presented larger diary average values from January to February, always with positive values, consequence of the rice crop irrigation. Two days with high values of flow of latent heat were obtained, what motivated the analysis of the energy balance. The closing of the energy balance was not reached in both days. A possible reason for that could be the low estimate of flows when the eddy covariance method was applied. / Este trabalho quantifica a troca líquida de CO2 (NEE) em uma área de cultivo de arroz irrigado no RS, como parte do projeto CT-HIDRO, projeto este que tem como finalidade descrever condições de superfície para diferentes ecossistemas do país. Os dados utilizados foram coletados durante o período de cultivo, entre outubro de 2003 e março de 2004. Neste período a temperatura média do ar e a radiação incidente são adequadas para o desenvolvimento da cultura nesta região. Os fluxos turbulentos, calculados pelo método de covariância de vórtices (eddy covariance), superestimaram o valor da NEE, provavelmente devido, a subestimativa da taxa de respiração, pois em alguns períodos noturnos os sensores registraram pouca turbulência. Além disso, outro fator que pode interferir no resultado dos fluxos, como por exemplo, o intervalo em que os fluxos são calculados. Foi aplicada a correção de Webb na curva de NEE e encontrou-se um valor menor que 1%, fato este que não modificou significativamente o resultado encontrado. O comportamento do fluxo de CO 2 durante todo o período analisado, mostrou um pico de -0,47mgCO2 m−2 s −1 , em janeiro de 2004, o que demonstra a maior captura de CO2 da atmosfera pela planta durante o dia e uma menor liberação à noite, devido a fase de desenvolvimento em que a planta se encontra. Neste caso, a correção de Webb também foi pouco acentuada, não interferindo no resultado obtido. As curvas que demonstram as três fases de desenvolvimento do arroz mostraram claramente uma maior absorção de CO 2 durante a fase reprodutiva. O comportamento do fluxo de CO 2 mostra um fluxo negativo durante à dia, quando a planta absorve CO2 da atmosfera e um fluxo positivo durante à noite quando ocorre o processo inverso, o que era esperado. O fluxo de calor latente apresentou maiores valores médios diários entre os meses de janeiro e fevereiro, apresentando sempre valores positivos, conseqüência da irrigação do campo de arroz. Obtivemos dois dias com valores elevados de fluxo de calor latente, o que fez com que analisássemos o balanço de energia. O fechamento do balanço de energia não foi alcançado em ambos os dias. Uma possível razão para isto pode ser o devido à subestimativa dos fluxos quando aplicado o método de covariância de vórtices.

In this work the predictive capability of a number of Reynolds stress transport(RST) models was first tested in a range of non-equilibrium homogeneous flows, comparisons being drawn with existing direct numerical simulation (DNS) results and physical measurements. The cases considered include both shear and normally strained flows, in some cases with a constant applied strain rate, and in others where this varied with time. Models were generally found to perform well in homogeneous shear at low shear rates, but their performance increasingly deteriorated at higher shear rates. This was attributed mainly to weaknesses in the pressure-strain rate models, leading to over-prediction of the shear stress component of the stress anisotropy tensor at high shear rates. Performance in irrotational homogeneous strains was generally good, and was more consistent over a much wider range of strain rates. In the experimental plane strain and axisymmetric contraction cases, with time-varying strain rates, there was evidence of an accelerated dissipation rate generation. Significant improvement was achieved through the use of an alternative dissipation rate generation term, Pε , in these cases, suggesting a possible route for future modelling investigation. Subsequently, the models were also tested in the inhomogeneous case of pulsating channel flow over a wide range of frequencies, the reference for these cases being the LES of Scotti and Piomelli (2001). A particularly challenging feature in this problem set was the partial laminarisation and re-transition that occurred cyclically at low and, to a lesser extent, intermediate frequencies. None of the models tested were able to reproduce correctly all of the observed flow features, and none returned consistently superior results in all the cases examined. Finally, models were tested in the case of a plane jet interacting with a rectangular dead-end enclosure. Two geometric configurations are examined, corresponding a steady regime, and an intrinsically unsteady regime in which periodic flow oscillations are experimentally observed (Mataoui et al., 2003). In the steady case generally similar flow patterns were returned by the models tested, with some differences arising in the degree of downward deflection of the impinging jet, which in turn affected the level of turbulence energy developing in the lower part of the cavity. In the unsteady case, only two of the models tested, a two-equation k-ε model and an advanced RST model, correctly returned purely periodic solutions. The other two RST models, based on linear pressure-strain rate terms, returned unsteady flow patterns that exhibited complex oscillations with significant cycle-to-cycle variations. Unfortunately, the limited availability of reliable experimental data did not allow a detailed quantitative examination of model performance.

Pas de résumé / In order to give a general statistical description of turbulence, one tries to identify universal statistical features, common to a wider class of turbulent flows. In 1988, Kraichnan and Panda discovered one such possibly universal feature, namely, the fact that Navier-Stokes turbulence tends to reduce the strength of the nonlinearity in its governing equations. In the flrst part of the manuscript we consider the strength of the nonlinear term and, more precisely, of its fluctuations in isotropic turbulence. In order to measure this strength, we compare to the case of a flow fleld with the same energy distribution where the modes are statistically independent, as is the case in Gaussian noise. It is shown that the turbulent flow self-organizes towards a state in which the nonlinearity is reduced, and it is discussed what the implications of this reduction are. Also, in two dimensions it is illustrated how this self-organization manifests itself through the appearance of well-deflned vortical flow structures. In the second part of the manuscript, we investigate the dynamics of the Hasegawa- Wakatani model, a model relevant in the study of magnetically conflned fusion plasmas. The two-dimensional version of this model is considered, which includes some key features of the turbulent dynamics of a tokamak-edge. We consider the limit of the model in which the nonlinearity is reduced with respect to the linear forces. For this weakly nonlinear, wave dominated regime, analytical predictions suggest the presence of a feedback loop in which energy is transferred to highly anisotropic zonal flows by nonlocal interactions. We confirm these predictions and we demonstrate a strong suppression of the turbulent radial particle flux. In wall bounded geometry, the same mechanism is observed and here also the flux is eflciently reduced by the turbulence-zonal flow interaction.

Les travaux numériques et expérimentaux de cette thèse contribuent à une meilleure compréhension de la dynamique de grosses particules dans un écoulement turbulent. Un premier volet m’a permis de quantifier leur mouvement relatif au fluide, ainsi que leur influence locale sur l’écoulement turbulent. Dans un second volet, j'ai trouvé que l'effet collectif des particules est d'atténuer l’amplitude des fluctuations turbulentes. En revanche, celles-ci n’ont pas d’influence sur les propriétés statistiques fines de l’écoulement. Aussi, ces mesures suggèrent l’existence d’une transition de phase dans les grandes échelles de l’écoulement au-delà d’un seuil critique du nombre de particules. / The numerical and experimental work of this thesis contribute to a better understanding of the dynamics of large particles in a turbulent flow. The first part allowed me to quantify their relative motion to the flow and their local influence on the surrounding flow. In a second part, I found that the collective effect of particles is to reduce the amplitude of turbulent fluctuations. In revanche, they have no influence on the fine statistical properties of the flow. Also, these measurements suggest the existence of a phase transition in the larger scales of the flow beyond a critical threshold of the number of particles.

Yes / A set of shallow-water equations (SWEs) based on a k-epsilon Reynold stress model is established to simulate the turbulent flows over a complex roughness bed. The fundamental equations are discretized by the second-order finite-difference method (FDM), in which spatial and temporal discretization are conducted by staggered-grid and leap-frog schemes, respectively. The turbulent model in this study stems from the standard k-epsilon model, but is enhanced by replacing the conventional vertical production with a more rigorous and precise generation derived from the energy spectrum and turbulence scales. To verify its effectiveness, the model is applied to compute the turbulence in complex flow surroundings (including a rough bed) in an abrupt bend and in a natural waterway. The comparison of the model results against experimental data and other numerical results shows the robustness and accuracy of the present model in describing hydrodynamic characteristics, especially turbulence features on the complex roughness bottom. / National Key Research and Development Program of China (Grant No: 2016YFE0122500, 2013CB036401 and 2013CB036402), China Postdoctoral Science Foundation (Grant No: 2016M591184) and Programme of Introducing Talents of Discipline to Universities (Grant No: BC2018038) / Research Development Fund Publication Prize Award winner, June 2019.

Settling and consolidation of suspended clay particles are significant issue in many fields such as geotechnical engineering, coastal and hydraulic engineering, and environmental engineering. A comprehensive literature review was conducted on the settling, consolidation and erosion of mixed soil material (cohesive and non-cohesive). Soil beds formed by sedimentation process of loose particles will be either show a segregated or homogeneous in structure, depending on the depositional environment. These sediments initially undergo self-weight consolidation and may be eroded under high flow rate. A number of studies have recently investigated the characteristic of consolidated clay bed in stagnant water. Hence, consolidation parameters were determined using a well-known vertical settling column consolidation test setup. However, limited research studies are available for deposition and consolidation of a mixture of sediment (clay, silt and sand) under flow conditions which are more representative of what happens in nature. A long flume and pump were used to create different turbulent conditions and simulate the natural process... The results for deposition and consolidation of different mixtures under stagnant and turbulence conditions were analyzed and compered in term of compressibility, permeability as well as shear strength. The results of this experimental research program indicated that the flow rate, initial concentration, height of settling and composition of sediment are all important factors that could affect the final bed dry unit weight. Two non-intrusive techniques were applied for measuring the dry unit weight at settling and consolidation stages. Impact echo technique has never been applied to measure the dry unit weight of self-weight consolidation along the vertical stratification of cohesive and non-cohesive particles. Also, a novel conductance sensor has been developed to improve the efficiency of this technique. The limitations of using these techniques will be highlighted in this study.