The two-dimensional depth-averaged shallow water equations have attracted considerable attentions as a practical way to solve flows with free surface. Compared to three-dimensional Navier-Stokes equations, the shallow water equations give essentially the same results at much lower cost. Solving the shallow water equations by the Godunov-type finite volume method is a newly emerging area. The Godunov-type finite volume method is good at capturing the discontinuous fronts in numerical solutions. This makes the method suitable for solving the system of shallow water equations. In this dissertation, both the shallow water equations and the Godunov-type finite volume method are described in detail. A new surface flow routing method is proposed in the dissertation. The method does not limit the shallow water equations to open channels but extends the shallow water equations to the whole domain. Results show that the new routing method is a promising method for prediction of watershed runoff. The method is also applied to turbulence modeling of free surface flow. The κ - ε turbulence model is incorporated into the system of shallow water equations. The outcomes prove that the turbulence modeling is necessary for calculation of free surface flow. At last part of the dissertation, a total load sediment transport model is described and the model is tested against 1D and 2D laboratory experiments. In summary, the proposed numerical method shows good potential in solving free surface flow problems. And future development will be focusing on river meandering simulation, non-equilibrium sediment transport and surface flow - subsurface flow interaction.

This study is part of a project that was conducted by Florida International University (FIU), which designed, built, and characterized a constructed wetland (CW) wastewater treatment system in El Salvador, Central America. This study focuses on the detailed performance of a field-scale CW in the tropics, consisting of a subsurface flow wetland (SSF) and a surface flow wetland (SF). The SSF had a total area of 151.2 sq. m., filled with gravel and planted with Phragmites, Thalia, and Brachiaria, along three independent cells operating with a flow in parallel and receiving the influent domestic wastewater from a facultative lagoon and preceded the SF. The SF was a shallow open basin containing alternating regions of open water (195 sq. m. total) with submerged aquatic plants (Elodea) and regions of shallower water (605 sq. m. total) filled with emergent wetland vegetation (Typha, Thalia, and Cyperus). The design, construction, startup and operation phases of the field scale CW in that tropical setting are thus described with detail, including both SSF and SF characterization of influent and effluents. The SSF average hydraulic detention times during the wet and dry season were 2 days ± 0.9 days and 4 days ± 0.4 days, respectively; and the SF average hydraulic detention times during the wet and dry season were 20 days ± 11.1 days and 77 days ± 19.5 days, respectively. Brachiaria/Cyperus presented better results during the wet season with average BOD5 residuals of 36% ± 25%, and Thalia/Thalia during the dry season with average BOD5 residuals of 33% ± 22%. Phragmites/Typha presented better results during both seasons with average TSS residuals of 2% ± 3% (wet season), and 2% ± 2% (dry season). Residuals are also presented for COD, Oil & Grease, Total Fecal Coliform, Total Phosphorus, and Total Nitrogen. In addition, an assessment of the empirical models used in the design of the system is completed, having the EPA methodology as the preferred for BOD5 removal and three methodologies for TSS removal under tropical climate conditions. A comparison of the differences in treatment associated with each one of the selected plants and their combinations is also discussed. In summary, results strongly suggest that the CW system can effectively reduce contaminants in wastewater to levels that are comparable with the objective levels (i.e., secondary treatment levels).

This thesis describes a system for dealing with free surface fluid simulations, and the components needed in order to construct such a system. It builds upon recent research, but in a computer graphics context the amount of available literature is limited and difficult to implement. Because of this, the text aims at providing a solid foundation of the mathematics needed, at explaining in greater detail the steps needed to solve the problem, and lastly at improving some aspects of the animation process as it has been described in earlier works. The aim of the system itself is to provide visually plausible renditions of animated fluids in three dimensions in a manner that allows it to be usable in a visual effects production context. The novel features described include a generalized interaction layer providing greater control to artists, a new way of dealing with moving objects that interact with the fluid and a method for adding source and drain capabilities.

Le rotomoulage réactif est un procédé de fabrication de pièces creuses en polymère où la synthèse du matériau intervient pendant la mise en œuvre. Cette méthode présente plusieurs avantages comparée à la méthode traditionnelle utilisant des poudres thermoplastiques : réduction du temps de cycle, utilisation possible de matériaux techniques, et baisse de la consommation d'énergie et du coût des matières premières. Cependant le rotomoulage réactif est plus complexe à mettre en œuvre car la polymérisation provoque un changement important et rapide de la viscosité. Une des solutions pour optimiser ce procédé est de simuler l'écoulement du système réactif pendant la mise en œuvre.Pour ce travail nous avons utilisé un polyuréthane thermodurcissable. Des analyses thermiques et rhéologiques ont permis d'étudier les phénomènes de gélification et de vitrification du matériau et le diagramme Temps-Température-Transformation a été établi. Le comportement rhéocinétique du système a également été modélisé.Le procédé a été simulé en utilisant un code de calcul basé sur la méthode « Smoothed Particle Hydrodynamics » (SPH). Ce code a été développé par notre équipe et plusieurs améliorations ont été apportées au cours de cette étude. Pour effectuer des simulations plus réalistes en utilisant un plus grand nombre de particules, la première amélioration a consisté à accélérer la résolution des calculs. Ensuite l'évolution de la viscosité a été prise en compte grâce à l'utilisation d'un modèle rhéocinétique et une nouvelle condition limite a été développée pour simuler l'adhésion du polymère sur la paroi du moule. Les modifications nécessaires à la simulation d'écoulements 3D ont également été apportées au code SPH. / Reactive rotational molding is a process to manufacture hollow plastic parts where synthesis occurs during the shaping. This method has several advantages compared to traditional rotomolding using thermoplastic powders: shorter cycle time, possible use of high performance materials, and decrease of energy consumption and raw materials costs. However reactive rotational molding is more complex to implement mainly because of the important and quick change of viscosity occurring during polymerization. One of the solutions to optimize this process is to simulate the reactive system flow during processing.In this work we used thermoset polyurethane as reactive system. Thanks to thermal and rheological analysis, gelation and vitrification phenomena were studied and Time-Temperature-Transformation diagram was established. Material chemiorheological behavior was also modeled.The process has been simulated using a solver based on Smoothed Particle Hydrodynamics (SPH) method. This solver was developed in our research team and several improvements have been added during this study. To be able to simulate realistic flows with a high number of particles, the first improvement was to accelerate the resolution of calculations. Then the change of viscosity has been taken into account using a chemiorheological model and a new boundary condition was developed to simulate adhesion of polymer on the mold wall. To be able to simulate 3D flows, the needed modifications have been also added to the SPH solver.

This project deals with the CFD modelling of a free surface flow. The aim is to develop and validate a fast and accurate numerical model for stratified two-phase flows. Volume of Fluid (VOF) multiphase model is employed. The purpose is to use the developed numerical model for the design of an element within a compact nuclear reactor.Unsteady Reynolds Averaged Navier-Stokes (RANS) simulations are conducted. Two free surface test cases are simulated to verify and ensure robustness of the model: a dam break and a vertical cylindrical obstacle set in a channel. From there, an optimization is performed in order to find the best compromise between accuracy and rapidity with the solver. The proper set of parameter models is found by carrying out extensive sensitivity studies and compare the solutions with available measurements.The obtained numerical results show a reasonable good agreement with the experimental data for the dam-break. Significant time savings are achieved thanks to the implemented optimization process while maintaining accuracy. The optimized model is then applied to the second test case and comparisons with experimental measurements are carried out. The same physical behavior of the flow as in experiments is captured with the simulations. The differences found between the simulation data and experiments are partly due to the difficulty to monitor experimentally with a high accuracy the highly non uniform regions within the flow.

La rencontre entre une couche limite laminaire se développant sur une plaque plane sous une surface libre et un obstacle rectangulaire émergent mène à l’apparition d’une zone de recirculation complexe appelée tourbillon en fer à cheval à l’amont de l’obstacle. Cette zone de recirculation est composée d’un certain nombre de vortex s’enroulant autour de l’obstacle et pouvant présenter des dynamiques variées. Le but de ce travail est de réaliser une étude expérimentale paramétrique du tourbillon en fer à cheval laminaire et faiblement turbulent. Les mesures sont réalisées par PIV et par trajectographies et des méthodes de détection basés sur la décomposition modale (POD) et la détection des points critiques sont utilisées afin de faciliter l’étude de la dynamique tourbillonnaire. Ce manuscrit se propose de : (i) Déterminer l’évolution des caractéristiques du tourbillon en fer à cheval en fonction des paramètres adimensionnels de l’écoulement. (ii) Mettre en place une typologie des différentes dynamiques exhibées. (iii) Déterminer les mécanismes à l’origine de l’apparition de ces différentes dynamiques à l’aide d’un modèle conceptuel. (iv) Et finalement, étudier l’influence de la longueur de l’obstacle sur le tourbillon en fer à cheval. / A laminar boundary layer developing under a free-surface and interacting with an emerging obstacle lead to the creation of a complex recirculation zone called horseshoe vortex upstream from the obstacle. This recirculation zone is composed of a given number of vortices wrapping around the obstacle with various kind of dynamics. The aim of this work is to perform an experimental parametric study of the laminar, and slightly turbulent horseshoe vortex. Measurements are obtained by PIV and trajectographies and detection method based on modal decomposition (POD) and critical points are used to ease the study of the vortex dynamics. This manuscript intends to : (i) Explain the horseshoe vortex main characteristics evolution with the non-dimensional parameters of the flow. (ii) Define a typology of the different observed vortex dynamics. (iii) Identify the physical mechanisms at the origin of those dynamics through the use of a conceptual model. (iv) And finally, show the influence of the obstacle length on the horseshoe vortex.

The marine current turbine (MCT) is an exciting proposition for the extraction of renewable tidal and marine current power. However, the numerical prediction of the performance of the MCT is difficult due to its complex geometry, the surrounding turbulent flow and the free surface. The main purpose of this research is to develop a computational tool for the simulation of a MCT in turbulent flow and in this thesis, the author has modified a 3D Large Eddy Simulation (LES) numerical code to simulate a three blade MCT under a variety of operating conditions based on the Immersed Boundary Method (IBM) and the Conservative Level Set Method (CLS). The interaction between the solid structure and surrounding fluid is modelled by the immersed boundary method, which the author modified to handle the complex geometrical conditions. The conservative free surface (CLS) scheme was implemented in the original Cgles code to capture the free surface effect. A series of simulations of turbulent flow in an open channel with different slope conditions were conducted using the modified free surface code. Supercritical flow with Froude number up to 1.94 was simulated and a decrease of the integral constant in the law of the wall has been noticed which matches well with the experimental data. Further simulations of the marine current turbine in turbulent flow have been carried out for different operating conditions and good match with experimental data was observed for all flow conditions. The effect of waves on the performance of the turbine was also investigated and it has been noticed that this existence will increase the power performance of the turbine due to the increase of free stream velocity.

Submitted by Debora de Oliveira Medeiros null (deboraomedeiros@gmail.com) on 2017-06-13T19:08:30Z No. of bitstreams: 1 dissertacao_arquivar.pdf: 4524135 bytes, checksum: 33a027db23ff40bebc026b9090443d53 (MD5) / Approved for entry into archive by Luiz Galeffi (luizgaleffi@gmail.com) on 2017-06-19T12:58:15Z (GMT) No. of bitstreams: 1 medeiros_do_me_prud.pdf: 4524135 bytes, checksum: 33a027db23ff40bebc026b9090443d53 (MD5) / Made available in DSpace on 2017-06-19T12:58:15Z (GMT). No. of bitstreams: 1 medeiros_do_me_prud.pdf: 4524135 bytes, checksum: 33a027db23ff40bebc026b9090443d53 (MD5) Previous issue date: 2017-04-24 / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / Neste trabalho é apresentado um estudo das equações de Navier Stokes incompressível com superfície livre e métodos de projeção com uma formulação recente denominada laplaciano superficial. Esta formulação tem como finalidade uma melhor descrição da força de tensão superficial e grandezas, como curvatura e vetores tangencial e normal, descritas sobre a interface. Assim, uma condição de tensão normal alternativa é definida, e esta quando discretizada implicitamente e combinada com o método de projeção, descrevem a formulação laplaciano superficial, cuja sua solução é utilizada como condições de contorno para resolver o sistema de equações que descrevem o escoamento. A nova formulação destina-se a resolver um sistema tridiagonal de equações gerado sobre a interface, e usar a solução deste sistema de equações como uma condição de contorno na superfície livre para o sistema linear da correção da pressão no interior do domínio, que é resultante da aplicação do método de projeção. A nova equação que define a condição de tensão normal conta com grandezas definidas na malha euleriana, no contexto Marker-And-Cell (MAC), que devem ser projetadas sobre a malha lagrangeana e também considera a curvatura e os vetores tangente e normal na sua descrição, sendo importante um estudo detalhado de geometria diferencial. Finalmente, variações da formulação laplaciano superficial com tensão superficial para diferenças finitas são aplicadas para resolver os testes numéricos da oscilação da gota e da gota apoiada que possuem solução de referência, além da simulação de um problema com movimento de interface (fountain flow). Nestes testes, concluímos que a variação mais precisa e estável é aquela que aplica uma discretização da equação da tensão normal utilizando médias de valores alocados na malha MAC. / This work presents a study of the Navier-Stokes equations incompressible with free surface and a projection methods with a recent formulation defined as surface laplacian. The purpose of this formulation is to improve the description of the force of stress tension and quantities, as curvature and tangent and normal vectors, present at the interface. Thus, an alternative normal stress condition is defined, and when this is discretized implicitly and combined with the employ of projection method, describe the surface laplacian formulation, whose solution is used as boundary condition to solve the system of equations describing the flow. The new formulation is intended to solve the tridiagonal system of equations generated at the interface, and to use this solution as a boundary condition at free surface for the linear system of the pressure correction inside of domain, which results of the application of the projection method. The new equation used to define the normal stress condition considers quantities defined in the Eulerian mesh, in the Marker-And-Cell context (MAC), that should be projected on the Lagrangian mesh. In addition, in the new equation, it is also considered the influence of the curvature and normal and tangential vectors in your description, so that a detailed study of differential geometry for this computation is important. Finally, variations of the surface laplacian formulation for finite differences are applied for solving numerical tests of the drop oscillation and sessile drop which have reference solutions, beyond of the simulation of a problem with free surface moving (fountain flow). In these tests, we concluded that the most accurate and stable variation is the one that applies a discretization of the normal stress equation using the mean of values in the MAC mesh. / FAPESP: 2015/01243-0

Submitted by (lucia.rodrigues@ufrpe.br) on 2016-07-11T14:45:22Z No. of bitstreams: 1 Flavio Adriano Marques.pdf: 3199315 bytes, checksum: a208fb8314e3ffa0d451678350503edb (MD5) / Made available in DSpace on 2016-07-11T14:45:22Z (GMT). No. of bitstreams: 1 Flavio Adriano Marques.pdf: 3199315 bytes, checksum: a208fb8314e3ffa0d451678350503edb (MD5) Previous issue date: 2013-02-21 / The Brazilian region known as Brejo Paraibano, in recent decades has been undergoing a process of replacing its native coverage by crops, mainly for pasture for extensive livestock farm, which has favored the occurrence of soil degradation in the form of water erosion. Given this context, this study aimed to quantify the rates of interrill erosion under simulated rainfall and evaluate the roughness parameters vegetable emerged to laminar flow in a pasture area hilly relief. The experiment was conducted in a completely randomized design with 4 treatments, 15%, 25%, 35% and 45% slope in a Ultisol under pasture, 5 replications, totaling 20 plots. The mean flow velocity and infiltration rates of pasture varied significantly with increasing slope contributing to increased erosion rates up to 35%. In the condition of 45% slope was not the greatest soil losses have occurred because of the removal of the horizon, in which the erosion process is acting on the argilic B horizon. The vegetation drag coefficient to pasture expressed elevation when there was a decrease of turbulent flow between plants. In the condition height flow increase, there was a greater energy flow between the mass of water and the plants structure of Brachiaria decumbens that resulted in decreasing on the force drag tension of plants. / A microrregião do Brejo Paraibano nas últimas décadas vem sofrendo um processo de substituição de sua cobertura nativa por culturas agrícolas, sobretudo por pastagens para exploração de pecuária extensiva, o que vem favorecendo a ocorrência de degradação do solo na forma de erosão hídrica. Diante desse contexto, este estudo teve como proposta quantificar as taxas de erosão em entressulcos sob chuva simulada e avaliar os parâmetros da rugosidade vegetal emersa ao escoamento em área de pastagem em relevo movimentado. O experimento foi conduzido em delineamento inteiramente casualizado nas condições de declive: 15%, 25%, 35% e 45%, em um Argissolo Vermelho-Amarelo sob pastagem, em 5 repetições, totalizando 20 parcelas experimentais. A velocidade média do escoamento e as taxas de infiltração na pastagem variaram significativamente com a elevação do declive contribuindo para elevação das taxas de erosão ate o declive de 35%. Na condição de 45% de declive não ocorreram as maiores perdas de solo em virtude de ter ocorrido a remoção do horizonte A, estando o processo erosivo atuando sobre o horizonte B textural. O coeficiente de arraste vegetal da pastagem expressou elevação quando ocorreu diminuição da turbulência do escoamento entre as plantas. Com o aumento da lâmina de escoamento, houve um maior fluxo de energia entre a massa de água e a estrutura vegetal da gramínea Brachiaria decumbens, que se refletiram na diminuição dos valores de tensão de arraste na planta.

Les rivières forment un système dynamique complexe soumis à des variations de grande ampleur, en effet les précipitations sont jugées comme la cause fondamentale de ces fluctuations. Au cours du temps, la morphologie des rivières évolue sous l’influence de plusieurs paramètres, en particulier les crues, les ouvrages hydrauliques, le transport sédimentaire. Le développement de la végétation dans le lit de la rivière et sur les berges peut affecter les conditions hydrodynamiques et le comportement des cours d'eau, ainsi l'impact de la végétation est une question cruciale pour la gestion des réseaux d'irrigation et des flux naturels. En réduisant la vitesse, la présence de végétation peut augmenter les dépôts de sédiments et modifier le risque d'inondation du fait des effets combinés de l'augmentation de la rugosité et de la diminution de la zone d'écoulement du chenal principal du fleuve. Ces aspects sont mis en avant par l’application de simulations numériques à des cas réels : cas de la Medjerda et du Canal Medjerda Cap Bon en Tunisie. La végétation est une caractéristique commune des eaux côtières et fluviales naturelles, interagissant à la fois avec le débit d'eau et le transport de sédiments. Cependant, les processus physiques qui régissent ces interactions sont encore mal compris, ce qui rend difficile la prévision du transport des sédiments et de la morpho dynamique. L’enjeu de cette thèse est d’améliorer la connaissance des processus physiques régissant les interactions entre végétation, écoulement et transport sédimentaire. Le but final est de pouvoir améliorer la gestion des hydro-systèmes artificiels ou naturels. Ce travail impliquera deux approches complémentaires d’expériences et de modélisation analytique et numérique. Dans un premier temps, on s’attachera à mieux caractériser les processus physiques d’interaction entre végétation et écoulement. Pour cela des expériences sur différents canaux fourniront l’hydrodynamique au-dessus de végétations modèles. On mettra l’accent sur le développement des méthodes expérimentales spécifiques à l’étude de couche limite au-dessus de macro-rugosités. Ces résultats seront dans un second temps analysés à partir de modèles analytiques qui permettent d’avoir les relations hauteur-débits nécessaires pour la gestion. Les caractéristiques et les performances de plusieurs modèles seront évaluées en regard des différents types de végétations. Dans un troisième temps, des expériences avec des sédiments préciseront l’influence de la végétation sur la modification du transport solide. La réduction des contraintes sur les lits engendre une adaptation nécessaire des lois de transport classique. Un modèle d’ajustement de ces lois sera proposé. / Rivers form a complex dynamic system subject to wide variations, in fact precipitation is considered as the fundamental cause of these fluctuations. Over time, the morphology of rivers evolves under the influence of several parameters, especially floods, hydraulic structures, sediment transport. The development of vegetation in the river bed and on the banks can affect the hydrodynamic conditions and the behavior of a watercourse, so the impact of vegetation on sediment transport is a crucial issue for the management of irrigation networks and natural flows. By reducing velocity, the presence of vegetation can increase sediment deposition and modify the risk of flooding due to the combined effects of increasing roughness and decreasing of flowing area of the river main channel. These aspects are highlighted by the application of numerical simulations to real cases: the case of the Medjerda River and the Channel of Medjerda Cap Bon in Tunisia. Vegetation is a common feature of natural coastal and riverine waters, interacting with both water flow and sediment transport. However, the physical processes governing these interactions are still poorly understood, making it difficult to predict sediment transport and morpho dynamics. The aim of this thesis is to improve the knowledge of the physical processes governing the interactions between vegetation, flow and sediment transport on the one hand, and to select the appropriate model that will be applied in the real case of rivers. The ultimate purpose is to improve the management of artificial or natural hydro-systems. This work will involve two complementary approaches to experiments and analytical and numerical modeling. At first, we will focus on better characterizing the physical processes of interaction between vegetation and flow. For this, experiments on different channels will provide hydrodynamics over model vegetation’s. Emphasis will be placed on the development of experimental methods specific to boundary layer studies over macro-roughness. These results will be analyzed in a second time from analytical models that allow the stage-discharge relationships required for management. The characteristics and performances of several models will be evaluated with regard to different types of vegetation. Thirdly, experiments with sediments will specify the influence of vegetation on the modification of solid transport. The reduction of the constraints on the beds generates a necessary adaptation of the classical transport laws. An adjustment model for these laws will be proposed.