O presente trabalho consiste em uma extensão do sistema FreeFlow-2D para simular escoamentos de fluidos não newtonianos bidimensionais com superfí cies livres, onde o fluido é descrito pelos modelos de Cross ou o modelo ``power-law\'\'. O método numérico empregado é o método GENSMAC. As equações governantes são aproximadas pelo método de diferenças finitas em uma malha deslocada e partículas marcadoras são utilizadas para a visualização do escoamento e localização da superfície livre. Resultados numéricos são apresentados. Em particular, a presente implementação é validada comparando-se a solução numérica com uma solução analítica / This work presents an extention of the Freeflow-2D system to non-Newtonian free surface flows. The governing equations are solved by the finite difference method on a staggered grid. Marker particles are used to describe the fluid providing the location and the visualization of the free surface. The methodology employed is based on the GENSMAC method. The fluid is modelled by the Cross and power-law models. Numerical examples are presented. The code is validated by making a comparison between analytical and numerical solutions

No description available.

Στόχος της παρούσας διδακτορικής διατριβής είναι η μελέτη χρονικά μεταβαλλόμενων και ασταθών υπό συνθήκες ροών ιξωδοπλαστικών και ιξωδο-ελαστικών ρευστών λόγω της σημασίας που παρουσιάζουν τα υλικά αυτά σε διάφορες βιομηχανικές διεργασίες. Στα πλαίσια της εργασίας αυτής μελετήθηκε αρχικά το πρόβλημα της συμπίεσης ιξωδοπλαστικού ρευστού μεταξύ δύο δίσκων το οποίο είναι ένα βασικό πείραμα ρεολογικού χαρακτηρισμού μη Νευτωνικών υλικών. Η χρονικά μεταβαλλόμενη προσομοίωση επέτρεψε τον προσδιορισμό των σημαντικών διαφορών μεταξύ των δύο εκδοχών του συγκεκριμένου ρεολογικού πειράματος (δίσκοι κινούμενοι υπό σταθερή ταχύτητα ή υπό σταθερή δύναμη) που ήταν αδύνατος με τα προηγούμενα μοντέλα, ενώ παρουσιάζεται και πλήρης παραμετρική μελέτη της διεργασίας αυτής. Το υπόλοιπο μέρος της παρούσας εργασίας αφιερώθηκε στη μελέτη της διεργασίας εκβολής ιξωδοελαστικών υλικών η οποία συναντάται ευρέως στη βιομηχανία μορφοποίησης πολυμερών. Παρουσιάζεται μια πλήρης παραμετρική ανάλυση για την εκβολή ενός ιξωδοελαστικού ρευστού Phan-Tien & Tanner (PTT) από ένα δακτυλιοειδή αγωγό προκειμένου να εξετάσουμε την επίδραση της γεωμετρίας του αγωγού, και των ρεολογικών ιδιοτήτων του ρευστού. Επιπλέον ένα σημαντικό πρόβλημα που εμφανίζεται κατά την εκβολή από ένα κυλινδρικό, επίπεδο ή δακτυλιοειδή αγωγό είναι η εμφάνιση ασταθειών οι οποίες επηρεάζουν σημαντικά την ποιότητα του τελικού προϊόντος. Προκειμένου να διερευνηθεί ο μηχανισμός που προκαλεί την εμφάνιση αυτών των ασταθειών πραγματοποιήθηκε γραμμική ανάλυση ευστάθειας για τη αξονοσυμμετρική και καρτεσιανή δισδιάστατη ροή επικόλλησης-ολίσθησης (stick-slip flow) ενός ρευστού PTT. Γύρω από τη λύση μόνιμης κατάστασης πραγματοποιείται γραμμική ανάλυση ευστάθειας, και υπολογίζονται οι ιδιοτιμές του γενικευμένου προβλήματος με τη βοήθεια της μεθόδου Arnoldi. Έτσι προσδιορίζεται η επίδραση των ιδιοτήτων του υλικού στην ευστάθεια ή μη της ροής. / The purpose of this dissertation is to perform a study on time-dependent and conditionally unstable flows of viscoplastic and viscoelastic fluids because of their importance in various industrial processes. A typical experiment which is used widely for the rheological characterization of non-Newtonian fluids is squeeze flow which is the subject of the first chapter of this dissertation. The time-dependent simulation permitted the determination of the distinct differences between the two versions of this rheological experiment (disks moving with constant velocity or under constant force) which was impossible with the models used up to now, while a complete parametric analysis of this process is presented. The rest of this work focuses in the important problem of the extrusion process of a viscoelastic fluid, which is frequently encountered in the polymer industry. A complete parametric analysis is presented for the extrusion of a Phan-Thien Tanner (PTT) fluid from an annular die in order to examine the effect of the die geometry and the rheological properties of the fluid on this process. Moreover, it is widely known that during the extrusion from a cylindrical, planar or an annular die various flow instabilities may arise affecting significantly the quality of the final product. In order to investigate the mechanism, which causes the appearance of these instabilities, a linear stability analysis is performed for the cylindrical or planar stick-slip flow of a PTT fluid. The stability analysis is performed around the steady state solution and the eigenvalues of the generalized problem are calculated using the Arnoldi algorithm. With this method the effect of the various rheological properties of the fluid on the stability of the flow is determined.

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

O presente trabalho consiste em uma extensão do sistema FreeFlow-2D para simular escoamentos de fluidos não newtonianos bidimensionais com superfí cies livres, onde o fluido é descrito pelos modelos de Cross ou o modelo ``power-law\'\'. O método numérico empregado é o método GENSMAC. As equações governantes são aproximadas pelo método de diferenças finitas em uma malha deslocada e partículas marcadoras são utilizadas para a visualização do escoamento e localização da superfície livre. Resultados numéricos são apresentados. Em particular, a presente implementação é validada comparando-se a solução numérica com uma solução analítica / This work presents an extention of the Freeflow-2D system to non-Newtonian free surface flows. The governing equations are solved by the finite difference method on a staggered grid. Marker particles are used to describe the fluid providing the location and the visualization of the free surface. The methodology employed is based on the GENSMAC method. The fluid is modelled by the Cross and power-law models. Numerical examples are presented. The code is validated by making a comparison between analytical and numerical solutions

Bien que répandus dans l'environnement et dans l'industrie, les écoulements granulaires possèdent des caractéristiques particulières qui sont différentes de celles des écoulements de liquides ordinaires comme l'eau (fluides newtoniens). L'une de ces caractéristiques est de présenter un seuil d'écoulement, il est donc fréquent qu'une des frontières de l'écoulement soit constituée de grains au repos (frontière érodable). L'échange entre les deux états solide et liquide d'un écoulement granulaire est à l'origine du phénomène d'érosion et de déposition et constitue le cœur de cette thèse. Nous y effectuons à la fois des études expérimentales et des études théoriques en nous basant sur les modèles phénoménologiques concernant le taux d'échange proposés par Bouchaud, Cates, Ravi Prakash et Edwards et modifiés par Boutreux et de Gennes. A l'aide d'un dispositif expérimental permettant de lâcher sans vitesse initiale une colonne de grains sur une pente granulaire, nous avons quantifié l'importance relative des phénomènes d'érosion et de déposition notamment en déterminant l'angle neutre qui définit la pente critique pour laquelle ces deux processus s'équilibrent. Dans un deuxième temps nous avons appliqué le modèle proposé par Boutreux et de Gennes au cas de l'étalement d'une « marche » granulaire (marche du Sinaï). / Granular flows, which are common in nature and industry, have particular characteristics that are different from the ordinary flow of liquid (eg. water, oil, etc.) and are not yet well understood in the literature. This case underlines the lack of knowledge on the exchange rate between solid and liquid states of granular flows which is the origin of the erosion - deposition phenomenon and constitutes the core of this thesis. To address the issue of the exchange rate solid - liquid, this work is based mainly in two subjects. The first subject is to study the processes of erosion and deposition of a flow generated by the grains stocked in a reservoir. By releasing naturally with zero-velocity the grains to generate the granular flow, we quantified the relation between the erosion and deposition for determining the neutral angle which defines the slope where these two processes balance. The second subject is to analyze the fall of grains in Sinai step problem by studying the spread of the march. For both works, we conduct not only the experimental but also theoretical studies based on the phenomenological models of the exchange rate proposed by BCRE (Bouchaud, Cates, Ravi Prakash and Edwards) and BDG (Boutreux and de Gennes) with some modifications.

Les écoulements en milieu urbain sont complexes et à l’heure actuelle estimés à l’aide d’outils informatiques. Pourtant, le manque de données expérimentales sur des géométries urbaines rend la validation et l’encadrement de l’utilisation de ces derniers difficile. Cette thèse présente les résultats obtenus sur un modèle physique d’un quartier urbain. La distribution des hauteurs d’eau ainsi que la répartition des débits en sortie du quartier expérimental sont mesurées. Leur étude a mis en évidence certains comportements caractéristiques des écoulements. Les données expérimentales ont été comparées aux simulations numériques générées avec un code 3D (Ansys-Fluent®) et un outil de recherche (Neptune 2D) mis au point durant cette thèse. Ce dernier résout les équations de Barré de Saint Venant 2D à l’aide d’un schéma EVR-DG, associé à une modification des solveurs de Riemann qui rend le code de calcul well-balanced.Les écarts observés entre Ansys-Fluent® et l’expérimental sont majoritairement en-dessous de 10%. Le code Neptune 2D apparait quant à lui légèrement moins précis : les écarts peuvent atteindre 20 à 30%. Diverses hypothèses sont avancées pour expliquer ces écarts. / This study deals with urban floods. Nowadays, numerical tools are used to simulate those complex flows. Nevertheless, the lack of experimental make the validation of the softwares difficult. This work presents experimental results of an urban flood physical model. The water height and the outflows at every outlet are measured and compared to numerical results. This study highlights some observed specificities of urban flows. In order to simulate those flows, a numerical tool (Neptune 2D) was developed during this PhD to solve the 2D shallow water equations with an EVR-DG scheme. Modifications of the Riemann solvers lead to a wellbalanced scheme. Numerical results were also provided with a 3D software (Ansys-Fluent®). Differences between Ansys-Fluent® and experimental results are mainly under 10%. Neptune2D is less accurate, with differences reaching 20 to 30%. Some hypotheses are discussed to explain those discrepencies.

Este trabalho é dedicado à análise e implementação de esquemas "upwind" de alta ordem modernos e o modelo de turbulência k-epsilon padrão no Freeflow-2D; um ambiente integrado para simulação numérica em diferenças finitas de problemas de escoamentos incompressíveis com superfícies livres. O propósito do estudo é a simulação de escoamentos de fluidos newtonianos incompressíveis, bidimensionais, confinados e/ou com superfícies livres e a altos valores do número de Reynolds. O desempenho do código Freeflow-2D atual é avaliada na simulação do escoamento numa expansão brusca e de um jato livre incidindo perpendicularmente sobre uma superfície rígida impermeável. O código é então aplicado na simulação de um jato planar turbulento em uma porção de fluido com superfície livre e estacionário. Os resultados numéricos obtidos são comparados com dados experimentais, soluções analíticas e soluções numéricas de outros trabalhos. / This work is devoted to the analysis and implementation of modern high-order upwind schemes and the standard k-epsilon turbulence model into the Freeflow-2D; a finite difference integrated environment for the numerical simulation of incompressible free surface flow problems. The purpose of this study is the two-dimensional simulation of high-Reynolds incompressible newtonian confined and/or free surface flows. The performance of the current Freeflow-2D code is assessed by applying it to the simulation of flow over a backward facing step and of an impinging free jet onto an impermeable rigid surface. The code is then applied to a turbulent planar jet into a pool. The numerical results are compared with experimental data, analytical solution, and numerical simulations of other works.

Neste trabalho apresenta-se um novo esquema prático tipo upwind de alta resolução, denominado EPUS (Eight-degree Polynomial Upwind Scheme), para resolver numericamente equações de conservação TVD e é implementado no contexto do método das diferenças finitas. O desempenho do esquema é investigado na resolução de sistemas hiperbólicos de leis de conservação e escoamentos incompressíveis complexos com superfícies livres. Os resultados numéricos mostraram boa concordãncia com outros resultados numéricos e dados experimentais existentes / Is this work a new practical high resolution upwinding scheme, called EPUS (Eight-degree Polynomial Upwind Scheme), for the numerical solution of transient convection-dominated conservation equations is present. The scheme is based on TVD stability criterion and is implemented in the context of the finite difference methodology. The performance of the scheme is investigated by solving hyperbolic systems of conservation laws and complex incompressible flows with free surfaces. The numerical results displayed good agreement with other existing numerical and experimental data

Três esquemas de captura de descontinuidade são apresentados para simular hiperbólicos de leis de conservação e equações de Navier-Stokes incompressíveis, a saber: FDHERPUS (Five Degree Hermite Upwind Scheme); RUS (Rational Upwind Scheme); e CSPUS (Cubic Spline Polynomial Upwind Scheme). Esses esquemas são baseados nos critérios de estabilidade CBC e TVD e implementados nos contextos das metodologias diferenças finitas e volumes finitos. A precisão local dos esquemas é verificada acessando o erro e a taxa de convergência em problemas testes de referência. Um estudo comparativo entre os esquemas estudados (incluido o WENO5) e o esquema bem estabelecido de van Albada, para resolver leis de conservação lineares e não lineares, é também realizado. O esquema de convecção que fornece melhores resultados em leis de conservação hiperbólicas é então examinado na simulação de escoamentos de fluidos newtonianos com superfícies livres móveis de complexidade crescente; resultados satisfatórios têm sido observados em termos do comportamento global / Three shock capturing schemes for numerical solution of hyperbolic conservation laws and incompressible Navier-Stokes equations are presented, namely: FDHERPUS (Five Degree Hermite Polynomial Upwind Scheme); RUS (Rational Upwind Scheme); and CSPUS ( Cubic Spline Polynomial Upwind Scheme). These schemes are based on CBC and TVD stability criteria and implemented in the context of finite volume methodologies. The local observed accuracy of the schemes is verified by assessing the error and convergence rate on benchmark test cases. A comparative study between the schemes (including WENO5) and the well established van. Albada scheme to solve standard linear and nonlinear hyperbolic conservation laws is also accomplished. The scheme that has provided better results in hyperbolic conservation laws is then examined in the simulation of Newtonian moving free surface flows of increasing complexity, satisfactory agreement has been observed in terms of the overall behavior