O estudo de escoamentos de fluidos não-Newtonianos através de expansões desperta um grande interesse em pesquisadores nas diversas áreas da engenharia, devido a sua ampla aplicação em indústrias e no meio acadêmico. O objetivo principal desta Dissertação é simular problemas de escoamentos envolvendo fluidos viscoplásticos através de expansões axissimétricas abruptas. O modelo mecânico empregado é baseado nas equações de conservação de massa e de momentum para escoamentos isocóricos acoplados com a equação constitutiva de um Fluido Newtoniano Generalizada (GNL), com a função de viscosidade de Herschel-Bulkley regularizada pela equação de Papanastasiou. O modelo mecânico é aproximado por um modelo estabilizado de elementos finitos, denominado método Galerkin Mínimos-Quadrados, ou Galerkin Least-squares (GLS). Esse método (GLS) é usado a fim superar as dificuldades numéricas do modelo de Galerkin clássico: a condição de Babuška-Brezzi e a instabilidade inerente em regiões advectivas do escoamento. O método é construído adicionando termos de malha-dependentes a fim aumentar a estabilidade da formulação de Galerkin clássica sem danificar sua consistência. A formulação GLS é aplicada para estudar a influência do índice power-law, da tensão limite de cisalhamento e razão de aspecto na dinâmica do escoamento de fluidos de Herschel-Bulkley através de expansões axissimétricas abruptas de razão de aspecto 1:2 e 1:4. Os problemas que envolvem números de Reynolds desprezíveis, para uma escala do número de Herschel-Bulkley entre 0 e 100 e índice de comportamento entre 0,2 e 1,0 são apresentados. Os resultados são fisicamente detalhados e estão de acordo com a literatura. / The study of non-Newtonian fluid flows in expansions is of great interest for researchers in the several branches of engineering, due to their wide application both in industry and academy. The objective of this Dissertation is to simulate flow problems involving a viscoplastic fluid through an axisymmetric abrupt expansion. The mechanical model employed is based on the mass and momentum conservative equations for isochoric flows coupled with the Generalized Newtonian Liquid (GNL) constitutive equation, with the Papanastasiou-regularized Herschel-Bulkley viscosity function. The mechanical model is approximated by a stabilized finite element scheme, namely the Galerkin Least-squares method. This method (GLS) is used in order to overcome the numerical difficulties of the classical Galerkin method: the Babuška- Brezzi condition and the inherent instability in advective flow regions. The method is built adding mesh-dependent terms in order to increase the stability of the classical Galerkin formulation without damaging its consistency. The GLS formulation is applied to study the influence of power-law index, yield stress and aspect reason in the flow dynamics of Herschel- Bulkley fluids through an axisymmetric abrupt expansions of aspect reason 1:2 and 1:4. Problems involving negligible Reynolds numbers, for a Herschel-Bulkley number range between 0 and 100 and e power-law index range between 0.2 and 1.0 are presented. The results are physically comprehensive and are in accordance with the literature.

A grande maioria dos fluidos encontrados na natureza se comportam como fluidos não- Newtonianos o que torna o seu estudo muito importante para diversas áreas da engenharia. Este trabalho tem como objetivo simular o problema específico de escoamentos de fluidos viscoplásticos em expansões abruptas planares com razão de aspecto de 1:4. O problema em questão se mostra interessante pois em muitos sistemas industriais são apresentados a geometria proposta para estudo associada a fluidos viscoplásticos. Foi empregado o modelo mecânico multi-campos (mult-field) baseado nas equações de conservação de massa e balanço de momentum para escoamentos isocóricos acoplados com a equação constitutiva de um Fluido Newtoniano Generalizada (GNL), associada à função de viscosidade de Herschel-Bulkley regularizada através da equação de Papanastasiou. O modelo mecânico é aproximado por um modelo estabilizado de elementos finitos, denominado método Galerkin Mínimos-Quadrados (GLS). A fim de se pesquisar os fenômenos reológicos ali presentes é feito o estudo da influência do índice de power-law na topologia de um escoamento creeping flow (Re@0) para uma vasta faixa de números de Herschel-Bulkley variando entre 0.1 e 100. Os resultados mostraram-se satisfátorios, apresentando uma forte influencia do número de Herschel-Bulkley e do índice de power law na topologia e na dinâmica do escoamento. Obteve-se uma validação do trabalho comparando-se os resultados obtidos nesta dissertação com os obtidos em artigo de grande credibilidade e aceito por toda a academia científica. / The majority of fluids found in nature behave like non-Newtonian fluids what makes their study of great importance to various areas of engineering. The present work aims to simulate the specific problem of the viscoplastic fluids flows through in planar abrupt expansion with the aspect ratio of 1:4. The problem in question is interesting because in many industrial systems are the proposal to study geometry associated with viscoplastic fluid. It employed the multi-field mechanical model based on equations of conservation of mass and momentum balance for the isochoric flow coupled with the constitutive equation of the Generalized Newtonian Liquids (GNL), associated with the function of viscosity of Herschel-Bulkley regularized by the equation of Papanastasiou. The mechanical model is approximated by a stabilized finite element model, called Galerkin Least-Squares method (GLS). In order to research the rheological phenomena present is done studying the influence of power-law index in the topology of the creeping flow (Re@0) for a wide range of numbers of Herschel-Bulkley ranging between 0.1 and 100. The results have proved satisfactory, showing a strong influence on the number of Herschel-Bulkley and power law index of the topology and dynamics of the flow. Got to be a validation of the work comparing the results with those obtained in this thesis in an article of great credibility and accepted throughout the scientific academy.

[pt] O escoamento de bolhas de gás em fluidos não-newtonianos é muito importante em diversas aplicações, tais como biorreatores, processamento de alimentos e operações de perfuração e cimentação de poços de petróleo, contudo poucos estudos abordam este tema. No presente trabalho é analisado o escoamento ascendente de bolhas de gás em fluido viscoplástico, modelado aqui como fluido Herschel-Bulkley. Utilizou-se uma abordagem bidimensional, aproximada por placas paralela. O escoamento multifásico foi numericamente simulado fazendo uso do método de volumes finitos VOF (volume of fluid), o qual resolve as equações de conservação de massa e momento aclopado a uma equação para a fração volumétrica dos fluidos. A influência de fatores tais como tensão limite de escoamento, dimensão da bolha, número de bolhas e distância entre as bolhas escoando em um fluido viscoplástico foram investigadas. Os resultados indicaram que a tensão limite de escoamento tem grande impacto na velocidade de deslocamento da bolha. No caso de mais de uma bolha escoando foi também observado que o deslocamento de uma bolha altera a velocidade de ascenção das outras, fazendo com que elas se unam, e a medida em que a distância entre as bolhas aumenta a interferência é eliminada e as mesmas escoam como bolhas individuais. Além disso, foi verificado que, dependendo do tamanho da bolha escoando, a parede interfere na sua velocidade de ascenção. Por fim foi analisada a mudança no formato da bolha, podendo-se observar que em número de Reynolds menores a mesma apresenta um formato esférico e a medida em que este parâmetro aumenta, a bolha vai se deformando e adquire uma forma elipsoidal. Os resultados qualitativos do presente estudo numérico foram comparados com alguns trabalhos experimentais encontrados na literatura e corresponderam relativamente bem. / [en] The gas bubble displacement in non-Newtonian fluids is important in many applications, including bioreactors, food processing and drilling and cementing of oil wells however a few studies have investigated this issue. The motion of gas bubbles in a viscoplastic fluid, modeled as a Herschel-Bulkley fluid, is analysed in the present work. A bidimensional approach was used, approximated by parallel plates. The multiphase flow was numerically simulated using the finite volume technique VOF (volume of fluid), which solves the conservation equations of mass and momentum coupled to an equation for the volume fraction of the fluids. Is was investigated the influence of factors such as fluid yield stress, bubble size, number of bubbles rising in the viscoplastic fluid and also the distance between them. The results indicated that the yield stress has great impact on the bubble rising velocity. In the case of multiple bubbles flowing it was also observed that the displacement of one bubble influences the rising velocity of the others, causing them to join together. As the distance between the bubbles increase the interference is eliminated and the bubbles flow separately. Furthermore, it was found that there was an wall interference depending on the size of the bubble. Lastly, it was analysed the change in the shape of the bubble. It could be observed that for lower Reynolds number the bubbles presented an spherical shape and as this parameter increases it turns into an elipsoidal shape. The qualitative results of the present study agreed relatively well with experimental works found in the literature.

O presente trabalho destina-se ao estudo numérico da geometria de tubos de seção elíptica que facilite a transferência de calor adimensional e diminua a queda de pressão adimensional (Δ̃) sofrida pelo escoamento. O método aplicado é o Design Construtal, que visa determinar a geometria que apresentará a menor resistência ao escoamento, ou seja, busca-se determinar a razão de aspecto da elipse (=⁄) que favorece a transferência de calor e diminui a queda de pressão do escoamento. O fluido empregado neste estudo apresenta características de viscoplasticidade. A relação entre a tensão cisalhante e a taxa de deformação obedece ao modelo de Herschel-Bulkley modificado. Considera-se que o escoamento é incompressível, laminar, bidimensional, externo e ocorre em regime permanente. A solução numérica do problema proposto é realizada com um código comercial baseado no método dos volumes finitos. É investigada a influência do índice de potência, , sobre a seção elíptica que facilita o escoamento e, para isso, este índice é variado de 0,4 a 1. A influência dos números de Reynolds (√), Herschel-Bulkley (√) e Prandtl (√) sobre o comportamento do escoamento também é avaliada. √ é variado de 1 a 40, √ é variado de 1 a 100 e √ é variado de 0,1 a 100 Os resultados mostram que, para um escoamento com √=1, √=1 e √=1, o aumento do índice de potência influencia negativamente na transferência de calor adimensional e a seção elíptica, que maximiza esta transferência de calor adimensional, tende a ser mais alongada na direção do escoamento. Já e influenciam positivamente na transferência de calor adimensional. Para um escoamento com √=1, √=1, =0,4 conclui-se que com o aumento de a razão de aspecto ótima (q,opt), do ponto de vista térmico, diminui. Quando é considerado um escoamento com √=1, √=1, =0,4 conclui-se que q,opt diminui com o aumento de , ou seja, a elipse torna-se mais alongada no sentido do escoamento. A variação de √ em um escoamento com √=1, √=1, =0,4 mostra que o aumento deste parâmetro acarreta em aumento da taxa de transferência adimensional e de Δ̃. / The present work is aimed at the numerical study of the geometry of elliptic section tubes that facilitates the dimensionless heat transfer and decreases the dimensionless pressure drop (Δ̃) suffered by the flow. The applied method is the Construtal Design, which aims to determine the geometry that will present the least resistance to the flow, that is, to determine the aspect ratio of the ellipse (=⁄) that favors heat transfer and decreases the flow pressure drop. The fluid used in this study has viscoplasticity characteristics. The relationship between shear stress and strain rate follows the modified Herschel-Bulkley model. It is considered that the flow is incompressible, laminar, two-dimensional, external and occurs in steady state. The numerical solution of the proposed problem is carried out with a commercial code based on the finite volume method. The influence of the power index, n, on the elliptical section facilitating the flow is investigated, and for this, the index is varied from 0.4 to 1. The influence of the Reynolds number (√), Herschel-Bulkley number (√) and Prandtl number (√) on the flow behavior is also evaluated √ is varied from 1 to 40, √ is varied from 1 to 100 and √ is varied from 0.1 to 100. The results show that for a flow with √=1, √=1 and √=1, the increase of the power index negatively influences the dimensionless heat transfer and the elliptic section, which maximizes this dimensionless heat transfer, tends to be more elongated in the direction of flow. Already √ and √ influence positively the dimensionless heat transfer. For a flow with √=1, √=1, =0.4 it is concluded that with the increase of √ the optimum aspect ratio (q,opt), from the thermal point of view, decreases. When a flow is considered with √=1, √=1, =0.4 it is concluded that q,opt decreases with the increase of √, that is, ellipse becomes more elongated in the flow direction. The variation of √ in a flow with √=1, √=1, =0.4 shows that the increase of this parameter causes an increase of the dimensionless transfer rate and Δ̃.

Cette thèse expérimentale et numérique s’intéresse à l’écoulement à inertie négligeable d’un fluide newtonien et d’un fluide à seuil autour d’une plaque ou d’un cylindre en présence d’une paroi. Les efforts générés et les structures d’écoulement ont été particulièrement étudiés. Pour le fluide à seuil, l’étude est centrée sur le domaine des effets de seuil grands par rapport aux effets visqueux. Pour l’approche expérimentale, un dispositif a été conçu permettant de mesurer simultanément les forces de traînée et de portance en fonction de la vitesse et de la distance obstacle/paroi. Les fluides étudiés sont un sirop de glucose Newtonien et un gel de Carbopol considéré comme un fluide à seuil modèle. Ces matériaux ont été finement caractérisés du point de vue rhéologique. Des simulations numériques utilisant le modèle viscoplastique anélastique d’Herschel-Bulkley régularisé par le modèle de Papanastasiou ont été réalisées. Les résultats obtenus avec le fluide newtonien ont permis de valider l’installation expérimentale et l’approche numérique. Les effets de la contrainte seuil, de la distance paroi/obstacle, de la rhéofluidification, de l’angle d’inclinaison de la plaque sur les coefficients de traînée et de portance ont été examinés en détail. Le régime de lubrification a été étudié pour le cylindre expérimentalement et théoriquement. La morphologie de l’écoulement (zones rigides) et les grandeurs locales ont été montrées. Des solutions analytiques de la traînée et de la portance ont été proposées. Les comparaisons entre les résultats expérimentaux et numériques ainsi qu’avec les solutions issues de la plasticité de la mécanique des sols sont analysées. Les écarts sont discutés en termes d’influence de l’élasticité et de plasticité du fluide à seuil. / This experimental and numerical thesis deals with the creeping flow of a Newtonian fluid and a yield stress fluid around a plate or a cylinder in the presence of a wall. The generated forces and the flow structures have been particularly studied. The study focuses for the yield stress fluid on large yield stress effects compared to the viscous effects. For the experimental approach, a device has been designed to measure simultaneously the drag and lift forces as a function of velocity and obstacle / wall distance. The fluids studied are a Newtonian glucose syrup and a Carbopol gel considered as a yield stress fluid model. These materials have been finely characterized from the rheological point of view. Numerical simulations using the Herschel-Bulkley anelastic viscoplastic model regularized by the Papanastasiou model were carried out. The results obtained with the Newtonian fluid allowed to validate the experimental set-up and the numerical approach. The effects of the yield stress, the wall / obstacle distance, the power law index and the inclined plate on the drag and lift coefficients were carefully examined. The lubrication regime has been studied for the cylinder both experimentally and theoretically. The morphology of the flow (rigid zones) and the local field quantities have been shown. Analytical solutions of drag and lift have been proposed. Comparisons between the experimental and numerical results as well as the solutions resulting from the plasticity of soil mechanics are analyzed. Deviations are discussed in terms of the elasticity and plasticity influence of the yield stress fluid.

Dans le présent travail, nous étudions l'impact sur un obstacle rigide d'un écoulement transitoire à surface libre de fluide viscoplastique. Cette étude est conduite numériquement à l'aide de la méthode SPH (Smoothed Particle Hydrodynamics), en y intégrant le modèle rhéologique de Herschel-Bulkley. Le code employé est adapté à nos besoins et validé sur des cas test classiques. Les caractéristiques locales de l'écoulement à proximité de l'obstacle sont analysées et deux régimes d'impact sont mis en évidence en fonction de la pente d'écoulement. L'étude des pressions exercées sur l'obstacle, conduite spatialement et temporellement en fonction de ces régimes d'impact, nous permet de mettre en évidence les rôles respectifs des composantes gravitationnelle et cinétique de la pression. Nos résultats sont comparés systématiquement à des résultats expérimentaux issus de travaux précédents et sont cohérents avec ces derniers. Une étude comparative de nos écoulements de fluide viscoplastique avec des écoulements de matériau granulaires de propriétés similaires nous conduit à mettre en évidence des caractéristiques communes entre les deux matériaux. / In this work, we study the impact of a transient free-surface flow of viscoplastic fluid on a rigid obstacle. This study is conducted numerically using the SPH (Smoothed Particle Hydrodynamics) method, and the Herschel-Bulkley rheological model. The SPH code is adapted to our needs and validated on classic benchmarks. The local characteristics of the flow near the obstacle are analysed and two impact regimes are highlighted depending on the slope angle. By studying of the pressure exerted on the obstacle, both spatially and temporally, with regards to these impact regimes, we evidence the respective roles of the gravitational and kinetic components of the pressure. Our results are systematically compared with experimental data from a previous work and are shown to be consistent. A comparative study conducted on both our viscoplastic flows and flows of granular material of similar properties highlights common characteristics of the two materials.

Les écoulements rampants, recirculants et instationnaires d'un fluide viscoplastique autour d'un cylindre ont été étudiés.Numériquement, les morphologies des écoulements, la localisation des zones rigides, les champs de contraintes et pression autour du cylindre ainsi que le coefficient de traînée, ont été déterminés sur un large domaine des nombres de Reynolds et d'Oldroyd.Expérimentalement, les fluides étudiés sont des gels de polymère Carbopol®. Le comportement élastoviscoplastique de ces gels a été modélisé par une loi d'Herschel-Bulkley adaptée. Le montage expérimental conçu et réalisé a été validé par l'étude de l'écoulement d'un fluide newtonien autour d'un cylindre et la mise en place d'une procédure adaptée pour les fluides à seuil.On a pu constater l'influence des conditions d'interface avec l'apparition d'une morphologie de lâchers de tourbillons simultanés et symétriques.

Dans le présent travail, nous étudions l'impact sur un obstacle rigide d'un écoulement transitoire à surface libre de fluide viscoplastique. Cette étude est conduite numériquement à l'aide de la méthode SPH (Smoothed Particle Hydrodynamics), en y intégrant le modèle rhéologique de Herschel-Bulkley. Le code employé est adapté à nos besoins et validé sur des cas test classiques. Les caractéristiques locales de l'écoulement à proximité de l'obstacle sont analysées et deux régimes d'impact sont mis en évidence en fonction de la pente d'écoulement. L'étude des pressions exercées sur l'obstacle, conduite spatialement et temporellement en fonction de ces régimes d'impact, nous permet de mettre en évidence les rôles respectifs des composantes gravitationnelle et cinétique de la pression. Nos résultats sont comparés systématiquement à des résultats expérimentaux issus de travaux précédents et sont cohérents avec ces derniers. Une étude comparative de nos écoulements de fluide viscoplastique avec des écoulements de matériau granulaires de propriétés similaires nous conduit à mettre en évidence des caractéristiques communes entre les deux matériaux. / In this work, we study the impact of a transient free-surface flow of viscoplastic fluid on a rigid obstacle. This study is conducted numerically using the SPH (Smoothed Particle Hydrodynamics) method, and the Herschel-Bulkley rheological model. The SPH code is adapted to our needs and validated on classic benchmarks. The local characteristics of the flow near the obstacle are analysed and two impact regimes are highlighted depending on the slope angle. By studying of the pressure exerted on the obstacle, both spatially and temporally, with regards to these impact regimes, we evidence the respective roles of the gravitational and kinetic components of the pressure. Our results are systematically compared with experimental data from a previous work and are shown to be consistent. A comparative study conducted on both our viscoplastic flows and flows of granular material of similar properties highlights common characteristics of the two materials.

Cette thèse est une étude expérimentale de la formation et de la maturation d’un anneau tourbillonnaire. Elle porte sur les écoulements faiblement inertiels (Reynolds : 5 à 500) en fluide newtonien puis non newtonien. Les anneaux sont générés par un système cylindre-piston. Ils sont analysés par visualisation et par vélocimétrie par images de particules (PIV). La dynamique en fluide newtonien à faible nombre de Reynolds se révèle plus complexe que celle à grands Reynolds avec l’apparition d’un anneau secondaire contrarotatif. Les résultats obtenus en fluide rhéofluidifiant montrent l’influence de l’indice de comportement ainsi que les zones de comportement rhéofluidifiant et newtonien. Les explorations en fluides viscoplastique et viscoélastique montrent la formation d’un ou plusieurs anneaux secondaires contrarotatifs, qui diffèrent dans leur formation et leur dynamique des anneaux observés en fluide newtonien et que l’on associe aux propriétés physiques intrinsèques du fluide. / The present work is an experimental study of the generation and the maturation of vortex rings, in order to characterize their structure and their global dynamics for small to moderate Reynolds numbers (between 5 and 500) in Newtonian and non-Newtonian fluids. The experimental set-up consists of a vertical cylindrical piston-tube system with the lower part immersed in a filled tank. Measurement campaigns have been carried out using dye visualization and Particle Image Velocimetry (PIV). A first part of the work is focussed on Newtonian fluid and allows the dynamics at low Reynolds numbers to be investigated qualitatively and quantitatively. This dynamics turns out to be more complex than the one classically observed at high Reynolds numbers, and is characterized by the production of a counter-rotating secondary vortex ring. The results obtained for shear thinning fluids show the influence of the power-law index on the development and the propagation of the ring. The computation of the shear rate field allows the results to be analyzed in terms of shear thinning and Newtonian regions. Finally a preliminary investigation for viscoplastic and viscoelastic fluids has been performed. In both cases, it is shown that one (for viscoplastic fluids) or several (for viscoelastic fluids) counter-rotating secondary vortex rings are generated, a phenomenon that can be associated with the intrinsic physical properties of the fluid. All these results provide several perspectives of study in the field of vortex rings dynamics in the weakly inertial regimes.

O presente trabalho destina-se ao estudo numérico da geometria de tubos de seção elíptica que facilite a transferência de calor adimensional e diminua a queda de pressão adimensional (Δ̃) sofrida pelo escoamento. O método aplicado é o Design Construtal, que visa determinar a geometria que apresentará a menor resistência ao escoamento, ou seja, busca-se determinar a razão de aspecto da elipse (=⁄) que favorece a transferência de calor e diminui a queda de pressão do escoamento. O fluido empregado neste estudo apresenta características de viscoplasticidade. A relação entre a tensão cisalhante e a taxa de deformação obedece ao modelo de Herschel-Bulkley modificado. Considera-se que o escoamento é incompressível, laminar, bidimensional, externo e ocorre em regime permanente. A solução numérica do problema proposto é realizada com um código comercial baseado no método dos volumes finitos. É investigada a influência do índice de potência, , sobre a seção elíptica que facilita o escoamento e, para isso, este índice é variado de 0,4 a 1. A influência dos números de Reynolds (√), Herschel-Bulkley (√) e Prandtl (√) sobre o comportamento do escoamento também é avaliada. √ é variado de 1 a 40, √ é variado de 1 a 100 e √ é variado de 0,1 a 100 Os resultados mostram que, para um escoamento com √=1, √=1 e √=1, o aumento do índice de potência influencia negativamente na transferência de calor adimensional e a seção elíptica, que maximiza esta transferência de calor adimensional, tende a ser mais alongada na direção do escoamento. Já e influenciam positivamente na transferência de calor adimensional. Para um escoamento com √=1, √=1, =0,4 conclui-se que com o aumento de a razão de aspecto ótima (q,opt), do ponto de vista térmico, diminui. Quando é considerado um escoamento com √=1, √=1, =0,4 conclui-se que q,opt diminui com o aumento de , ou seja, a elipse torna-se mais alongada no sentido do escoamento. A variação de √ em um escoamento com √=1, √=1, =0,4 mostra que o aumento deste parâmetro acarreta em aumento da taxa de transferência adimensional e de Δ̃. / The present work is aimed at the numerical study of the geometry of elliptic section tubes that facilitates the dimensionless heat transfer and decreases the dimensionless pressure drop (Δ̃) suffered by the flow. The applied method is the Construtal Design, which aims to determine the geometry that will present the least resistance to the flow, that is, to determine the aspect ratio of the ellipse (=⁄) that favors heat transfer and decreases the flow pressure drop. The fluid used in this study has viscoplasticity characteristics. The relationship between shear stress and strain rate follows the modified Herschel-Bulkley model. It is considered that the flow is incompressible, laminar, two-dimensional, external and occurs in steady state. The numerical solution of the proposed problem is carried out with a commercial code based on the finite volume method. The influence of the power index, n, on the elliptical section facilitating the flow is investigated, and for this, the index is varied from 0.4 to 1. The influence of the Reynolds number (√), Herschel-Bulkley number (√) and Prandtl number (√) on the flow behavior is also evaluated √ is varied from 1 to 40, √ is varied from 1 to 100 and √ is varied from 0.1 to 100. The results show that for a flow with √=1, √=1 and √=1, the increase of the power index negatively influences the dimensionless heat transfer and the elliptic section, which maximizes this dimensionless heat transfer, tends to be more elongated in the direction of flow. Already √ and √ influence positively the dimensionless heat transfer. For a flow with √=1, √=1, =0.4 it is concluded that with the increase of √ the optimum aspect ratio (q,opt), from the thermal point of view, decreases. When a flow is considered with √=1, √=1, =0.4 it is concluded that q,opt decreases with the increase of √, that is, ellipse becomes more elongated in the flow direction. The variation of √ in a flow with √=1, √=1, =0.4 shows that the increase of this parameter causes an increase of the dimensionless transfer rate and Δ̃.