<p>Many practical fluid flows involve liquid-particle systems and so there is a need to better understand the mechanism of particle deposition, adhesion, and agglomeration in suspensions, especially in complex geometries with moving boundaries and free surfaces. In this thesis, the nature of the particle-solid interactions and particle-fluid interactions is studied where the above complexities are present, taking into account particle collision, colloidal, and hydrodynamic forces, and two way coupling between the fluid flow and particles. The research is motivated by the industrial examples of: flow of dross particles near the sink roll surface in a galvanizing bath (moving surface), and the flow of particles in slot coating dies (free and moving surfaces). Particle motion and agglomeration play important roles in the example systems chosen for this fundamental 3-D study. Numerical studies of flow of dispersed suspensions makes it possible to understand the effects of flow conditions, particle characteristics, and flow geometry specifications that lead to agglomeration of particles in complex systems, especially where experimental studies are difficult to perform. Often the effects of these conditions are discovered due to process or product failures, rather than through insight into the processing steps.</p> <p>The modeling methodology used in this work is that micron sized spherical particles are tracked in the fluid phase by solution of Newton`s second law of motion for each particle. Fluid phase applies hydrodynamic forces on particles (drag, lift). Body forces, (soft sphere) particle-particle collisions and particle-wall collisions are considered. Particle concentrations are in the dilute regime between 0.01-5%vol. Flow of particles with the fluid phase is a fully coupled formulation in systems with particle concentrations > 1%vol.</p> <p>The thesis is organized around three example problems taken from industry that pose challenging modeling issues. The first involves particle collisions with a moving wall (dross particles in a zinc bath). The second problem includes particle-particle and particle-wall collisions in a turning flow geometry. The third problem, particle dispersion flows in a slot coating die, has the most complexity and includes particle-particle, particle-wall and free surfaces.</p> <p>Dross particle build up on the sink roll inside the zinc bath is an industrial problem that causes significant down time, and where an experimental study of the molten zinc in a bath running at C is difficult to perform. With the aid of computational fluid dynamics, turbulent flow of molten zinc in galvanizing bath is simulated, compared with previous cold model experiments, and coupled with the motion of dross particles around the sink roll. The presence of fixed position hardware and moving sink roll and guide rolls in a bath with dimensions in the orders of meters, and micron sized (20-100 ) dross particles makes this case a complex study. Drag, buoyancy, lift force and soft sphere nonlinear collision is considered in solution of Newton`s law of motion for each particle. Turbulent flow is simulated using a standard model. Simulations show regions on the sink roll where particles are dragged toward the surface of sink roll and have long residence times. These regions have been reported to experience large particle build-ups in the hot-dip galvanizing process.</p> <p>In another study, formation and breakage of agglomerates in a turning flow is studied. Neutrally buoyant particles with concentration of 5%vol are tracked in a fully coupled flow. Particles form agglomerates at the corner, where drag and lift force from the fluid breaks a number of agglomerates. The presence of a moving wall in the turning flow shifts the suspended particle formations toward the inside of channel. Location of particles agglomerates shifts toward the free surface with the presence of free surface at the turning flow.</p> <p>Motion of micron sized spherical particles with 1-4%vol through a slot die coating system is elucidated in a separate study. The system is complex with presence of moving web and free surface. Discrete element method (DEM) for motion of dispersed phase and volume of fluid (VOF) method for solution of continuous phase are integrated in a simulation study. Particles are 2-4 and the flow dimensions of the system are in the order of 100 . Particles experience collision, colloidal and hydrodynamic forces. Coupling between flow of particles and fluid phase is conducted. The results of this study show particle positions on the coating film can be predefined and depends on their initial positions within the feed slot. Particles agglomerate in recirculating regions of the coating gap and follow the streamlines of flow on the moving web. Regions in the coating gap where particles have high residence times (inside the die and near the feed slot edges) have particle agglomerations in the slot die coating system.</p> / Doctor of Philosophy (PhD)

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).

Le présent travail de recherche concerne l'étude de la convection de Rayleigh-Bénard-Marangoni dans un récipient cylindrique doté d'un fond en substrat solide. Le substrat solide est chauffé sur sa face inférieure par un flux de chaleur localisé. L'étude comporte deux parties : La première partie du travail consiste en une modélisation physique du problème associée à des simulations numériques. Les équations de Navier-Stokes et de l'énergie sont résolues en 3D par une méthode de volumes finis. Un transfert de chaleur conjugué solide-liquide est considéré. Des morphologies originales de cellules (type et nombre) sont observées, elles dépendent des conditions géométriques, des nombres adimensionnels qui régissent la physique de l'écoulement (nombre de Prandtl, de Rayleigh et de Marangoni ainsi que du rapport des conductivités thermiques du substrat solide et du fluide). Les transferts de chaleur sont aussi évalués pour chaque cas d'étude. Dans la deuxième partie, nous allons détaillons une étude expérimentale de la convection de Rayleigh-Bénard-Marangoni dans la même configuration que celle étudiée numériquement. Les structures convectives et leurs évolutions sont étudiées à partir d’images relevées par thermographie infra-rouge. Différents modes d'organisation des cellules convectives ont pu être mis en évidence pour ce type de chauffage à flux thermique imposé non uniforme. / The present research work concerns the study of Rayleigh-Bénard-Marangoni convection in a cylindrical container with a solid substrate base. This solid substrate is heated by a localized heat flux on its underside. The study is divided into two parts : The first part of the work consists of a physical modelling of the problem associated with numerical simulations. The Navier-Stokes and energy equations are solved by using a 3D finite volume method. A conjugate solid-liquid heat transfer is considered. Original morphology of cells (type and number) are observed, they are linked to the geometrical conditions, the dimensionless numbers which govern the physical problem (Prandtl, Rayleigh and Marangoni numbers and the ratio of solid substrate to liquid thermal conductivities). The heat transfers are also evaluated in each case. In the second part of the work, we present an experimental study of Rayleigh-Bénard-Marangoni convection in the same configuration as that studied numerically. Convective structures and their evolutions are studied from images recorded by infrared thermography. Different modes of organization of convective cells have been highlighted for this type of heating with imposed non-uniform heat flux.

Plusieurs études tant numériques qu’expérimentales font état de la présence d’instabilités thermiques dans des films liquides chauffés uniformément par le bas pour des conditions aux limites et d’écoulements particuliers. La présence de ces instabilités modifiera les transferts thermiques associés. Le sujet de ce travail de thèse consiste à étudier numériquement les instabilités thermoconvectives d’un écoulement laminaire tridimensionnel de convection mixte dans un canal horizontal à surface libre. Les variations de la tension de surface avec la température (effet Marangoni ou effet thermocapillaire) sont prises en compte. Bien que d’un intérêt certain pour de nombreuses applications industrielles, cette situation a été très peu étudiée d’un point de vue académique dans la configuration considérée ici. Dans cette de configuration plusieurs types de structures thermoconvectives sont susceptibles d’apparaître. Lorsque les forces induites par les courants de convection naturelle, forcée et thermocapillaire sont du même ordre de grandeur, les premiers résultats montrent un développement des instabilités sous forme de rouleaux convectifs longitudinaux stationnaires semblables à ceux rencontrés pour des écoulements de type Poiseuille-Rayleigh-Bénard. A notre connaissance, c’est la première fois que l’écoulement de convection de type Poiseuille-Rayleigh-Bénard associé aux effets Marangoni est étudié. Le nombre et la distribution spatiale des rouleaux convectifs le long du canal dépendent des conditions de l’écoulement. Nous proposons une étude numérique pour ces conditions particulières d’écoulement pouvant conduire à des instabilités avec une évaluation de leur effet sur les transferts de chaleur. Les équations de Navier-Stokes et de l’énergie sont résolues numériquement par la méthode de volumes finis en prenant en compte les effets thermocapillaires. Les résultats présentés concernent l’influence des paramètres contrôlant l’écoulement (nombres de Reynolds, de Rayleigh, de Biot, de Marangoni et le rapport de forme) sur les motifs de l’écoulement et les échanges thermiques. Dans une seconde partie du travail, complémentaire à la première, une analyse de stabilité linéaire de l’écoulement dans un canal ouvert à surface libre d’extension latérale infinie est réalisée en utilisant la méthode spectrale de type collocation Chebyshev pour résoudre un système aux valeurs propres. Les diagrammes de stabilité déterminant les seuils des paramètres conduisant à l’instabilité thermoconvective ont été obtenus et analysés, ainsi que les structures spatiales associées. / Several studies both numerical and experimental have reported the presence of thermal instabilities in liquid films uniformly heated from below for specific boundary conditions and flows. The presence of these instabilities modifies the associated heat transfer. The subject of this PhD thesis is to study numerically the instability of three-dimensional laminar mixed convection within a liquid flowing on a horizontal channel heated uniformly from below. The upper surface is free and assumed to be flat. The variations of the surface tension with the temperature (Marangoni effect or thermocapillary effect) are taken into account. Although of great interest for many industrial applications, this problem has received little attention from an academic point of view. In this configuration, several types of thermoconvective structures may appear. When the strength of the buoyancy, thermocapillary effects and forced convective currents are comparable, the results show the development of instabilities in the form of steady longitudinal convective rolls similar to those encountered in the Poiseuille-Rayleigh-Bénard flow. To our knowledge, this is the first time that the Poiseuille-Rayleigh-Bénard flow associated to the Marangoni effects has been investigated. The number and spatial distribution of the convective rolls along the channel depend on the flow conditions. We propose a numerical study on the flow conditions that could lead to thermal instabilities with an evaluation of their effect on the heat transfer. The coupled Navier-Stokes and energy equations are solved numerically by the finite volume method taking into account the thermocapillary effects. The results presented concern the influence of several control parameters (the Reynolds, Rayleigh, Biot and Marangoni numbers and the aspect ratio of the channel) on the flow patterns and heat transfer characteristics. In the second part of this work, complimentary to the first, a linear stability analysis of a horizontal liquid film flowing in an open channel, with infinite lateral extension and uniform heating from below, is carried out. An eigenvalue problem is obtained in the course of this analysis which is solved numerically using the Chebyshev collocation spectral method. The stability diagrams determining the threshold parameters leading to thermoconvective instabilities were obtained and analyzed as well as the associated spatial patterns.

Les cavités latérales sont des zones mortes à surface libre situées sur le côté d’un écoulement fluvial ou côtier. Les vitesses caractéristiques au sein de la cavité étant beaucoup plus faibles que celles de l’écoulement, une couche de mélange se développe à l’interface entre ces deux régions. Cette couche de mélange peut alors transférer de la quantité de mouvement de l’écoulement vers la cavité et ainsi mettre en mouvement la cavité et peut aussi transférer de la masse entre les deux régions, telle une pollution venant de l’écoulement amont. L’étude de cette thèse a alors consisté à étudier les caractéristiques de la couche de mélange, qui est rendue spécifique par le fait qu’elle se développe entre deux coins géométriques formés par l’intersection entre les parois de la cavité et celles de l’écoulement principal. Nous avons alors pu identifier l’origine et l’alternance des mouvements de fluide dans la direction transverse: de la cavité vers l’écoulement et inversement. Concernant la mise en mouvement de la cavité, le choix a été fait de considérer un écoulement principal fixé et de modifier l’extension de la cavité dans la direction perpendiculaire à l’écoulement, passant ainsi d’une cavité rectangulaire alignée avec l’écoulement principal à une cavité allongée dans le sens opposé. La mesure de champ de vitesse par PIV 2D a alors montré une forte évolution de la forme de l’écoulement à mesure que la géométrie de la cavité évolue : un système avec deux cellules alignées dans le sens de l’écoulement à un système à une seule cellule, puis un système à deux cellules et enfin un système complexe 3D ont ainsi été observés pour une cavité de plus en plus allongée. Ensuite, une modification du dispositif expérimental a permis de mesurer de deux façons différentes le transport de scalaire de l’écoulement principal vers la cavité, de comprendre les processus associés à ce transfert et enfin de quantifier cette capacité de transfert pour différents écoulements principaux et différentes géométries de cavités. Nous avons notamment montré que la géométrie de la cavité a peu d’effet alors que le nombre de Reynolds et la profondeur d’eau normalisée ont un effet majeur sur cette capacité de transfert de masse entre les deux régions. / Lateral cavities are free-surface dead-zones located on the side of a fluvial or coastal main flow. As the typical velocities are much larger in the main flow than in the cavity, a mixing layer appears at the interface between both regions. This mixing layer is able to transfer between the main flow and the cavity momentum which then sets the fluid in the cavity in motion and also passive scalar, such as a pollution coming from upstream. The objective of this work was then to investigate the characteristics of the mixing layer, which specificity comes from the fact that it is constrained between the upstream and downstream geometrical corners. It was possible to observe the origin and alternation of the transversal fluid motions: from the cavity towards the main flow and conversely. Regarding the motion in the cavity, the choice was made to keep a constant main flow and to measure the 2D horizontal velocity field using PIV as the extension of the cavity increases. The flow pattern then passes from a 2-cell patterns aligned in the direction of the main flow to a single-cell pattern, then a 2-cells patterns aligned along the direction perpendicular to the main flow and finally a complex 3D pattern for the widest cavity. Then a modification of the experimental set-up permitted to investigate the passive scalar exchanges from the main stream towards the cavity. It was possible to understand the processes responsible for such transfer and to quantify the transfer capacity. The analysis dimensional revealed that in the present subcritical, smooth simplified geometry cavity, the three parameters possible responsible for the modification of the transfer capacity are the geometrical aspect ratio of the cavity, the Reynolds number of the main flow and finally the normalized water depth. It was then shown that the impact of the cavity geometry remains negligible but that the Reynolds number and the normalized water do impact this passive scalar transfer capacity.

L’objectif de cette thèse est de développer un modèle numérique 3D afin d’étudier le phénomène de laves torrentielles ; écoulements visqueux fortement chargés en matière solide, surgissant en montagne lors d’orages violents. Aujourd’hui, la prévision des zones de vulnérabilité s’appuie sur des outils de calcul 0D, 1D ou 2D. Or ces outils ne peuvent représenter pleinement le comportement à surface libre des écoulements du fait de nombreuses approximations et hypothèses. Ainsi cette thèse met en œuvre un code numérique 3D pour étudier ces écoulements. Ce travail aboutit à la construction d’un modèle 3D à partir de données réelles de terrain. Plusieurs scénarios ont été étudiés et comparés à des résultats issus d’un modèle 2D. Les résultats mettent en évidence les apports non négligeables de la modélisation 3D : zones d’étalement et de dépôt, phénomènes d’encombrement, modélisation fine des écoulements dans les zones chenalisées. Préalablement, le modèle 3D a été validé en comparant les résultats numériques à des données expérimentales issues de la littérature, pour des typologies d’écoulement représentatives de celles observées sur des sites grandeur nature. / The objective of this thesis is to develop a 3D numerical model to assess debris flow. These viscous flows, heavily loaded with solid matter, form when heavy rain occurs in mountains. Today, forecasts of potentially impacted areas are based on 0D, 1D and 2D numerical tools. However, these tools cannot fully represent the free surface behaviour of debris flows due to the approximations and assumptions on which they are based. Thus, this work utilises a 3D numerical code to study this phenomenon. A specific model is built with real field data. Several flow scenarios are studied and compared with a 2D numerical model. The results highlight the significant benefits of a 3D approach by providing information on the fine representation of flow dynamics over the catchment area. The model also predicts the impact of debris flow (overflowing on a road bridge) and the zones of deposition and spreading. It highlights possible congestion phenomena and reproduces flows in the channels by fully accounting for parietal friction, capabilities not provided by 2D models. Prior to this application, the 3D model was evaluated with five sets of experimental data to validate its ability to represent viscoplastic flows. Different types of flows are studied and are representative of those observed on real sites when debris flow occur.

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

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

Este trabalho tem como objetivo desenvolver um método numérico para simular escoamentos incompressíveis, isotérmicos, confinados ou com superfícies livres, de fuidos viscoelásticos governados pelos modelos integrais de Maxwell e K-BKZ (Kaye-Bernstein, Kearsley e Zapas). A técnica numérica apresentada é uma extensão do método GENSMAC (Tomé McKee - J. Comp. Phys., (110), pp 171--186, 1994 ) para a solução das equações de conservação, juntamente com as equações constitutivas integrais de Maxwell e K-BKZ. As equações governantes são resolvidas pelo método de diferenças finitas em uma malha deslocada. O tensor de Finger, B_t\'(t) é calculado com base nas idéias do método de campos de deformação (Peters et al. - J. Non-Newtonian Fluid Mech. (89), de maneira que não há a necessidade de seguir a trajetória da partícula de fuido para descrever a história de deformação da partícula. Uma abordagem diferente para a discretização do instante passado é utilizada e o tensor de Finger e o tensor das tensões são calculados utilizando um método de segunda ordem. A validação do método numérico descrito nesse trabalho foi feita utilizando o escoamento em um canal bidimensional e a solução numérica obtida para a velocidade e para as componentes de tensão com o modelo de Maxwell foram comparadas com as respectivas soluções analíticas no estado estacionário, mostrando excelente concordância. Os resultados numéricos para a simulação do escoamento em uma contração planar 4 : 1 mostraram bons resultados, tanto qualitativos quanto quantitativos, quando comparados com os resultados experimentais de Quinzani et al. ( J. Non-Newtonian Fluid Mech. (52), pp 1?36, 1994 ). Além disso, utilizando os modelos Maxwel e K-BKZ, o escoamento em uma contração planar 4 : 1 foi simulado para vários números de Weissenberg e os resultados obtidos estão de acordo com os encontrados na literatura. Resultados numéricos de escoamentos com superfícies livres modelados pelas equações integrais de Maxwell e K-BKZ são apresentados. Em particular, a simulação numérica do jato oscilante para diferentes números de Weissenberg e diferentes números de Reynolds é apresentada. / The aim of this work is to develop a numerical technique for simulating incompressible, isothermal, free surface (also con¯ned) viscoelastic flows of fuids governed by the integral models of Maxwell and K-BKZ (Kaye-Bernstein, Kearsley and Zapas). The numerical technique described herein is an extension of the GENSMAC method (Tome and McKee, J. Comput. Phys., 110, pp. 171-186, 1994) to the solution of the momentuum and mass conservation equations together with the integral constitutive Maxwell and K-BKZ equations. The governing equations are solved by the finite difference method on a staggered grid using a Marker-and-Cell approach. The fluid is represented by marker particles on the fluid surface only. This provides the visualization and location of the fluid free surface so that the free surface stress conditions can be applied. The Finger tensor Bt0(t) is computed using the ideias of the deformation fields method (Peters et al. J. Non-Newtonian Fluid Mech., 89, pp. 209-228, 2001) so that it is not necessary to track a fluid particle in order to calculate its deformation history. However, in this work modifcations to the deformation fields method are introduced: the past time is discretized using a different formula, the Finger tensor Bt0(x; t) is obtained by a second order method and the stress tensor ? (x; t) is computed by a second order quadrature formula. The numerical method presented in this work is validated by simulating the flow of a Maxwell fluid in a two-dimensional channel and the numerical solutions of the velocity and the stress components are compared with the respective analytic solutions providing a good agreement. Further, the flow through a 4:1 planar contraction of a specific fuid studied experimentally by Quinzani et al. (J. Non-Newtonian Fluid Mech., 52, pp. 1-36, 1994) was simulated and the numerical results were compared qualitatively and quantitatively with the experimental results and very good agreement was obtained. The Maxwell and the K-BKZ models were applied to simulate the 4:1 planar contraction problem using various Weissenberg numbers and the numerical results were in agreement with those published in the literature. Finally, numerical results of free surface flows using the Maxwell and K-BKZ integral constitutive equations are presented. In particular, the numerical simulation of jet buckling using several Weissenberg numbers and various Reynolds numbers are presented

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