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

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.

This thesis consists of two papers dealing with a stabilization method for free-surface flows. The method was initially developed to stabilize mantle-convection simulations, but is in this work extended to ice-sheet- and glacier modeling. The objective of this thesis is to assess the method when applied glaciological simulations, with regards to stability and accuracy. It is shown that the method works well and increases stable time-step sizes substantially both for ice-sheet- and glacier simulations, without loss of accuracy. The increased stability properties might be useful for performing long-term simulations and increasing sea-level-rise predictions on a centennial time scale. / Denna avhandling består av två artiklar som inom ramen för ismodellering undersöker en stabiliseringsmetod för flöden med en fri yta. Metoden framtogs först för stabilisering av simuleringar av mantelkonvektion, men har i den här avhandlingen anpassats till ismodellering. Avhandlingens mål har varit att utvärdera metoden med avseende på stabiltet och noggrannhet. Från de utförda studierna framkommer det att metoden ökar längden på stabila tidssteg avsevärt, utan att nämnvärt påverka noggrannheten hos islösaren. De ökade stabilitetsegenskaperna hos metoden kan exempelvis innebära ökad nogrannheten i fastställandet av framtida havsnivåhöjning genom möjliggörandet av långtidssimuleringar på en tidsskala av flera hundra år.

No description available.

Convection in a layer of fluid with a free surface due to a combination of thermocapillary stresses and buoyancy is a classic problem of fluid mechanics. It has attracted increasing attentions recently due to its relevance for two-phase cooling. Many of the modern thermal management technologies exploit the large latent heats associated with phase change at the interface of volatile liquids, allowing compact devices to handle very high heat fluxes. To enhance phase change, such cooling devices usually employ a sealed cavity from which almost all noncondensable gases, such as air, have been evacuated. Heating one end of the cavity, and cooling the other, establishes a horizontal temperature gradient that drives the flow of the coolant. Although such flows have been studied extensively at atmospheric conditions, our fundamental understanding of the heat and mass transport for volatile fluids at reduced pressures remains limited. A comprehensive and quantitative numerical model of two-phase buoyancy-thermocapillary convection of confined volatile fluids subject to a horizontal temperature gradient has been developed, implemented, and validated against experiments as a part of this thesis research. Unlike previous simplified models used in the field, this new model incorporates a complete description of the momentum, mass, and heat transport in both the liquid and the gas phase, as well as phase change across the entire liquid-gas interface. Numerical simulations were used to improve our fundamental understanding of the importance of various physical effects (buoyancy, thermocapillary stresses, wetting properties of the liquid, etc.) on confined two-phase flows. In particular, the effect of noncondensables (air) was investigated by varying their average concentration from that corresponding to ambient conditions to zero, in which case the gas phase becomes a pure vapor. It was found that the composition of the gas phase has a crucial impact on heat and mass transport as well as on the flow stability. A simplified theoretical description of the flow and its stability was developed and used to explain many features of the numerical solutions and experimental observations that were not well understood previously. In particular, an analytical solution for the base return flow in the liquid layer was extended to the gas phase, justifying the previous ad-hoc assumption of the linear interfacial temperature profile. Linear stability analysis of this two-layer solution was also performed. It was found that as the concentration of noncondensables decreases, the instability responsible for the emergence of a convective pattern is delayed, which is mainly due to the enhancement of phase change. Finally, a simplified transport model was developed for heat pipes with wicks or microchannels that gives a closed-form analytical prediction for the heat transfer coefficient and the optimal size of the pores of the wick (or the width of the microchannels).