Escoamentos sem inércia de fluido elasto-viscoplástico, dentro de uma cavidade, são numericamente analisados. As soluções visam compreender a influência dos efeitos viscosos e elásticos na topologia de superfícies de escoamento. Assumindo-se que o colapso da microestrutura do material é instantâneo, o modelo mecânico é constituído pelas equações governantes de massa e momentum para fluidos incompressíveis, associado a uma equação hiperbólica para o tensor tensão extra, baseado na equação do modelo Oldroyd-B (Nassar et al, 2011). A principal característica do modelo é considerar a viscosidade e o tempo de relaxação como função da taxa de deformação, permitindo a pseudoplasticidade de viscosidade e restringindo os efeitos elásticos para as regiões não deformadas do material. As simulações numéricas são realizadas através do método de Galerkin mínimos quadrados a três campos: tensor tensão extra, pressão e velocidade. Os resultados mostram que as superfícies de escoamento do material são fortemente influenciadas pela ação combinada entre os efeitos elásticos e viscosos, estando em conformidade com a recente visualização experimental dos fluxos elasto-viscoplásticos. / Elasto viscoplastic uid ows without inertia, within a cavity, are numerically analyzed. The solutions aim to understand the in uence of viscous and elastic e ects on the topology of yield surfaces. Assuming that the collapse of the material microstructure is instantaneous, the mechanical model consists of the governing equations of mass and momentum for incompressible uids, associated with a hyperbolic equation for the extra stress tensor, based on the equation of the Oldroyd-B model (Nassar et al, 2011). The main feature of model is to consider the viscosity and the relaxation time as a function of shear rate, allowing the shear-thinning of viscosity and restricting elastic e ects for regions not deformed material. Numerical simulations are performed by the method of Galerkin Least Squares to three elds: extra stress tensor, pressure and velocity. The results show that the yield surfaces of material are strongly in uenced by the combined action between the elastic and viscous e ects, complying with the recent experimental visualization of elasto-viscoplastic ows.

Submitted by RODRIGO KOITI ISHIZAKA (rodrigo_ishizaka@hotmail.com) on 2018-04-09T20:10:27Z No. of bitstreams: 1 Dissertação Final PDF.pdf: 10461391 bytes, checksum: a78974d2db3d166540e03d4cee7def37 (MD5) / Rejected by ALESSANDRA KUBA OSHIRO ASSUNÇÃO (alessandra@fct.unesp.br), reason: Solicitamos que realize correções na submissão seguindo as orientações abaixo: - Número de página das referências no sumário está incorreto (consta p.86 quando é p.87); - No sumário consta Apêndices, porém no trabalho não. É finalizado nas referências. Agradecemos a compreensão. on 2018-04-10T12:44:45Z (GMT) / Submitted by RODRIGO KOITI ISHIZAKA (rodrigo_ishizaka@hotmail.com) on 2018-04-10T16:09:35Z No. of bitstreams: 1 Dissertação Final PDF.pdf: 10461372 bytes, checksum: 679a6e28fa80deb9079156b9b6b70c36 (MD5) / Approved for entry into archive by ALESSANDRA KUBA OSHIRO ASSUNÇÃO (alessandra@fct.unesp.br) on 2018-04-10T20:02:48Z (GMT) No. of bitstreams: 1 ishizaka_rk_me_prud.pdf: 10461372 bytes, checksum: 679a6e28fa80deb9079156b9b6b70c36 (MD5) / Made available in DSpace on 2018-04-10T20:02:48Z (GMT). No. of bitstreams: 1 ishizaka_rk_me_prud.pdf: 10461372 bytes, checksum: 679a6e28fa80deb9079156b9b6b70c36 (MD5) Previous issue date: 2018-02-28 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / O presente trabalho de mestrado consiste em apresentar uma modelagem para escoamentos isotérmicos viscoelásticos em que a viscosidade varia de acordo com a pressão e estudos numéricos para alguns problemas bidimensionais, como escoamentos entre placas paralelas e expansão planar 1:4. Para este trabalho é utilizado o modelo viscoelástico Oldroyd-B e uma modelagem linear da viscosidade com relação a pressão. O método numérico desenvolvido é baseado no método da projeção para desacoplar velocidade e pressão nas equações de Navier-Stokes e depois calcula os tensores pela equação constitutiva com a informação da viscosidade variando com a pressão. As equações são discretizadas em uma malha deslocada pelo método de diferenças finitas. Os resultados numéricos são obtidos das simulações de escoamentos em um canal e expansão planar 1:4 bidimensional, cujo foco é destacar algumas diferenças entre o modelo Oldroyd-B Standard e o modelo Oldroyd-B sob influência da viscosidade que varia linearmente com a pressão. Este trabalho de mestrado também tem por objetivo resolver problemas, com fluidos viscoelásticos, em que haja regiões com variações da pressão e estudar estes resultados numéricos comparando-os com o modelo Oldroyd-B Standard. / The present master's work consists in presenting a model for viscoelastic isothermal flows with pressure-dependent viscosity and numerical studies for some two-dimensional problems, such as channel flow and planar expansion 1:4. In this work will be considered the Oldroyd-B viscoelastic model and a linear modeling viscosity with a linear relationship between viscosity and pressure. The numerical method developed is based on the projection method to decouple velocity and pressure in the Navier-Stokes equations and then calculates the tensor from the constitutive equation taking account the pressure-dependent viscosity. The equations are discretized on a staggered grid by using the finite diference method. The numerical results are obtained for two problems: two-dimensional channel and planar expansion 1:4 flows, whose focus is to highlight some diferences between the standard Oldroyd-B model and the Oldroyd-B model under the influence of viscosity that varies linearly with pressure. This master's work aims to solve problems of viscoelastic fluids in which there exist regions where pressure varies and to study these numerical results comparing them with the standard Oldroyd-B model.

Escoamentos sem inércia de fluido elasto-viscoplástico, dentro de uma cavidade, são numericamente analisados. As soluções visam compreender a influência dos efeitos viscosos e elásticos na topologia de superfícies de escoamento. Assumindo-se que o colapso da microestrutura do material é instantâneo, o modelo mecânico é constituído pelas equações governantes de massa e momentum para fluidos incompressíveis, associado a uma equação hiperbólica para o tensor tensão extra, baseado na equação do modelo Oldroyd-B (Nassar et al, 2011). A principal característica do modelo é considerar a viscosidade e o tempo de relaxação como função da taxa de deformação, permitindo a pseudoplasticidade de viscosidade e restringindo os efeitos elásticos para as regiões não deformadas do material. As simulações numéricas são realizadas através do método de Galerkin mínimos quadrados a três campos: tensor tensão extra, pressão e velocidade. Os resultados mostram que as superfícies de escoamento do material são fortemente influenciadas pela ação combinada entre os efeitos elásticos e viscosos, estando em conformidade com a recente visualização experimental dos fluxos elasto-viscoplásticos. / Elasto viscoplastic uid ows without inertia, within a cavity, are numerically analyzed. The solutions aim to understand the in uence of viscous and elastic e ects on the topology of yield surfaces. Assuming that the collapse of the material microstructure is instantaneous, the mechanical model consists of the governing equations of mass and momentum for incompressible uids, associated with a hyperbolic equation for the extra stress tensor, based on the equation of the Oldroyd-B model (Nassar et al, 2011). The main feature of model is to consider the viscosity and the relaxation time as a function of shear rate, allowing the shear-thinning of viscosity and restricting elastic e ects for regions not deformed material. Numerical simulations are performed by the method of Galerkin Least Squares to three elds: extra stress tensor, pressure and velocity. The results show that the yield surfaces of material are strongly in uenced by the combined action between the elastic and viscous e ects, complying with the recent experimental visualization of elasto-viscoplastic ows.

Escoamentos sem inércia de fluido elasto-viscoplástico, dentro de uma cavidade, são numericamente analisados. As soluções visam compreender a influência dos efeitos viscosos e elásticos na topologia de superfícies de escoamento. Assumindo-se que o colapso da microestrutura do material é instantâneo, o modelo mecânico é constituído pelas equações governantes de massa e momentum para fluidos incompressíveis, associado a uma equação hiperbólica para o tensor tensão extra, baseado na equação do modelo Oldroyd-B (Nassar et al, 2011). A principal característica do modelo é considerar a viscosidade e o tempo de relaxação como função da taxa de deformação, permitindo a pseudoplasticidade de viscosidade e restringindo os efeitos elásticos para as regiões não deformadas do material. As simulações numéricas são realizadas através do método de Galerkin mínimos quadrados a três campos: tensor tensão extra, pressão e velocidade. Os resultados mostram que as superfícies de escoamento do material são fortemente influenciadas pela ação combinada entre os efeitos elásticos e viscosos, estando em conformidade com a recente visualização experimental dos fluxos elasto-viscoplásticos. / Elasto viscoplastic uid ows without inertia, within a cavity, are numerically analyzed. The solutions aim to understand the in uence of viscous and elastic e ects on the topology of yield surfaces. Assuming that the collapse of the material microstructure is instantaneous, the mechanical model consists of the governing equations of mass and momentum for incompressible uids, associated with a hyperbolic equation for the extra stress tensor, based on the equation of the Oldroyd-B model (Nassar et al, 2011). The main feature of model is to consider the viscosity and the relaxation time as a function of shear rate, allowing the shear-thinning of viscosity and restricting elastic e ects for regions not deformed material. Numerical simulations are performed by the method of Galerkin Least Squares to three elds: extra stress tensor, pressure and velocity. The results show that the yield surfaces of material are strongly in uenced by the combined action between the elastic and viscous e ects, complying with the recent experimental visualization of elasto-viscoplastic ows.

We prove a global-in-time and large-data existence of a suitable weak solution to a system of partial differential equations describing an unsteady flow of homogeneous incom- pressible viscoelastic rate-type fluid. The material parameters are continuous functions of temperature and, in particular, the dependence of the shear modulus is assumed to be linear. It is shown that studied models obey the fundamental laws of thermodynamics. The key step towards the existence proof is derivation of the balance of entropy. This in- equality is paramount in the analysis and as its consequence, we obtain sufficient a priori estimates, positivity of temperature and also regularity of the elastic deformation. The second part of the thesis deals with the existence analysis for the isothermal case, however using a completely different method, which is of independent interest. 1

[pt] Filmes finos de líquidos estão presentes em uma gama de aplicações industriais, como processos de atomização e revestimento de substrato. O processo de quebra pode ser divido em duas etapas: o estágio de ruptura, e o estágio de retração. O primeiro, movido pelas forças de van der Waals, ocorre quando uma pequena perturbação cresce e provoca o aparecimento de um pequeno furo no filme. O segundo, movido por forças capilares, provoca o crescimento desse furo levando à desintegração do filme de líquido. A estabilidade de uma cortina de líquido depende das características da perturbação, da espessura do filme e das propriedades do fluido. Análises experimentais mostraram que uma cortina super fina pode ser obtida pela utilização de fluidos viscoelásticos. Os mecanismos físicos associados à esta estabilidade, contudo, não são totalmente compreendidos. Este trabalho apresenta um estudo numérico e teórico dos efeitos das propriedades viscoelásticas na estabilidade de uma cortina de fluido, englobando ambos os estágio do processo. As análises numéricas foram desenvolvidas através da expansão assintótica das variáveis do escoamento com aplicação de um esquema de integração no tempo totalmente implícito. A partir da análise teórica da dinâmica de ruptura foi possível obter um critério de estabilidade linear para perturbações planares e axissimétricas em fluidos Newtonianos e não-Newtonianos. O tempo de ruptura e a velocidade de retração do filme foram calculados numericamente como função das propriedades viscoelásticas do líquido. Resultados mostraram que as forças elásticas atuam de forma a dificultar o processo de quebra e retração. Análises da evolução da espessura mostraram que as propriedades reológicas do fluído também interferem no formato que o filme de fluido assume durante o processo de retração. Para regimes de baixa viscosidade, as forças elásticas atuaram evitando a formação de ondas capilares observadas em fluidos Newtonianos. / [en] Thin free liquid sheets are ubiquitous in many industrial processes, such as atomization and curtain coating. Liquid sheets are susceptible to instabilities at the interface, which can grow, triggering a breakup process. This process can be divided into two different stages: the rupture stage and retraction. The first, driven by van der Waals force, occurs when a small instability grows until it pinches-off the sheet. The second, driven by capillary forces, induces the growth of the hole caused by the pinch-off, leading to the full disintegration of the liquid sheet. The stability of a liquid sheet depends on disturbance characteristics, sheet thickness, and fluid properties. Experimental analyses have shown that thinner stable liquid curtain can be obtained with viscoelastic liquids. The underlyning physical mechanisms associated with increased stability are, however, not fully understood. This work presents a theoretical and numerical analysis of the effect of viscoelasticity on the stability of a thin liquid sheet during both stages of the breakup process. We first analyze the rupture dynamics, deriving linear stability criteria for both planar and axisymmetric perturbations of Newtonian and Oldroyd-B liquids. The time evolution of planar and axisymmetric perturbations in an Oldroyd-B liquid sheet is evaluated using the asymptotic expansion of the flow variables and a fully-implicit time integration scheme. The rupture time and retraction velocity are calculated as a function of the viscoelastic properties. The results show that the liquid rheological behavior does not influence the linear stability criterion. Nevertheless, it has a strong effect on the growth rate of the disturbance and retraction velocity, increasing, thus, the breakup time. The results show that elastic forces act to hinder the rupture and retraction stages. Analysis of the temporal evolution of the thickness profile reveals that liquid rheological behavior also affects the shape of the liquid sheet. For low viscosity regime, the elastic forces damp the capillary waves that arise during the retraction of Newtonian sheets.

Recently, a new kind of turbulence has been discovered in the flow of concentrated polymer melts and solutions. These flows, known as purely elastic flows, become unstable when the elastic forces are stronger than the viscous forces. This contrasts with Newtonian turbulence, a more familiar regime where the fluid inertia dominates. While there is little understanding of purely elastic turbulence, there is a well-established dynamical systems approach to the transition from laminar flow to Newtonian turbulence. In this project, I apply this approach to purely elastic flows. Laminar flows are characterised by ordered, locally-parallel streamlines of fluid, with only diffusive mixing perpendicular to the flow direction. In contrast, turbulent flows are in a state of continuous instability: tiny differences in the location of fluid elements upstream make a large difference to their later locations downstream. The emerging understanding of the transition from a laminar to turbulent flow is in terms of exact coherent structures (ECS) — patterns of the flow that occur near to the transition to turbulence. The problem I address in this thesis is how to predict when a purely elastic flow will become unstable and when it will transition to turbulence. I consider a variety of flows and examine the purely elastic instabilities that arise. This prepares the ground for the identification of a three-dimensional steady state solution to the equations, corresponding to an exact coherent structure. I have organised my research primarily around obtaining a purely elastic exact coherent structure, however, solving this problem requires a very accurate prediction of the exact solution to the equations of motion. In Chapter 2 I start from a Newtonian ECS (travelling wave solutions in two-dimensional flow) and attempt to connect it to the purely elastic regime. Although I found no such connection, the results corroborate other evidence on the effect of elasticity on travelling waves in Poiseuille flow. The Newtonian plane Couette ECS is sustained by the Kelvin-Helmholtz instability. I discover a purely elastic counterpart of this mechanism in Chapter 3, and explore the non-linear evolution of this instability in Chapter 4. In Chapter 5 I turn to a slightly different problem, a (previously unexplained) instability in a purely elastic oscillatory shear flow. My numerical analysis supports the experimental evidence for instability of this flow, and relates it to the instability described in Chapter 3. In Chapter 6 I discover a self-sustaining flow, and discuss how it may lead to a purely elastic 3D exact coherent structure.

Cette thèse est consacrée à la vérification de l'analogie entre l'instabilité viscoélastique (VEI) et l'instabilité magnéto-rotationnel (MRI) dans un écoulement képlérien, afin de mieux comprendre le transport du moment dans les disques d'accrétion. Le discriminant de Rayleigh élasto-rotationnel est établi pour clarifier le rôle de l'élasticité dans le VEI. L'analyse de stabilité linéaire (LSA) avec le modèle d’Oldroyd-B est effectuée pour prédire les paramètres critiques des modes viscoélastiques. Il fait apparaître également l'influence de l'élasticité, la viscosité polymérique et d'autres paramètres de contrôle pour le VEI. Des expériences bien contrôlées avec des solutions aqueuses de polyoxyéthylène (POE) et de polyéthylène glycol (PEG) sont effectuées. Nous avons observé le mode stationnaire axisymétrique supercritique avec des solutions de faible élasticité et modes désordonnés sous-critiques avec des solutions de grande élasticité. Les formes et les valeurs critiques de ces modes sont en bon accord avec les prédictions théoriques de LSA. Selon l'analogie, le mode axisymétrique stationnaire est probablement l'analogue de MRI standard, tandis que le mode désordonné est probable que l'analogue de MRI hélicoïdale. La thèse contient aussi des résultats théoriques expérimentaux sur quatre autres régimes de rotation et un cas de limite d'élasticité infinie. / This thesis is devoted to the verification of the analogy between the viscoelastic instability (VEI) and the magnetorotational instability (MRI) in a Keplerian flow, in order to get better understanding of the momentum transportation in accretion disks.The elasto-rotational Rayleigh discriminant is deduced to clarify the role of the elasticity in the VEI. The linear stability analysis (LSA) with Oldroyd-B model is performed to predict critical parameters of viscoelastic modes, and it reveals the influence of the elasticity, polymer viscosity on the VEI. Experiments with well controlled aqueous solutions of polyoxyethylene (POE) and polyethylene glycol (PEG) are conducted. We have observed supercritical stationary axisymmetric mode with solutions of small elasticity and subcritical disordered modes with solutions of large elasticity. Both the flow patterns and the critical values of these modes are in good agreement with the LSA predictions. According to the analogy, the stationary axisymmetric mode is likely the analog of the standard MRI while the disordered mode is likely the analog of the helical MRI. The thesis contains also theoretical and experimental results with four other rotation regimes and the limit case of infinite elasticity.

企業評價對於投資決策有重大的影響，不論是發行上市、或是機構投資人選擇投資標的、乃至於併購或是清算，企業評價都是一切的基礎。再加上近來各界對於節能產業的重視，發光二極體封裝產業正如日中天的高度成長，如何能夠正確地衡量此產業的企業價值，實是機構投資人或是一般大眾關心的課題。此外，實務界長久以來詬病證管會所採用的承銷價格公式，乃是結合不同評價模式的方式來評斷發行股票公司之正確股票價值，但實證結果往往發現此公式會造成股價被低估的現象。 本研究以台灣地區共六家發光二極體封裝產業上市櫃公司為例，以其民國八十七年至九十四年的財務數據和資料，以五年為一階段，利用七種不同的評價模式：三階段成長現金流量折現法、三階段成長本益比法、三階段成長股價淨值比法、三階段股價銷售額比法、市場比較本益比法、市場比較股價淨值比法以及市場比較股價銷售額比法，預期九十二年初至九十五年初之理論實際股價，並與實際的市價作一比較，利用THEIL所提出的THEIL’S U值來比較不同評價模型與實際市價差距的績效，以選出最適合發光二極體封裝產業之企業評價模式。 本研究更進一步探討長久以來被實務界所詬病的綜合評價模式（結合不同的評價模式），試著經過第一階段實證結果的篩選，利用簡單權重結合本產業最佳和次佳的企業評價模式，以得到一個評價績效更勝於最佳評價模式的綜合評價法。 實證結果顯示，發光二極體最佳評價模型乃為市場比較股價銷售額比法(THEIL’S U=0.3515)，而三階段成長現金流量折現法，則適用於產業較穩定的情況下。突破性的發現則為，利用THEIL’S U值來比較評價績效而選出的最佳和次佳模型，在分別給予簡單權重（ex：50％：50％、60％：40％等）的情況下所得到的綜合評價法，其THEIL’S U值（<0.3515）比當初單一最佳評價法--市場比較股價銷售額比法(0.3515)還要來得低，顯示綜合評價法的有效性的確存在，並值得各界參考。此外，亦發現給予最佳評價法較大權重時，更可以進一步提昇綜合評價法之績效。此結果反駁了實務界長久以來對於綜合評價法的不信任，也給予證管會一個修正承銷價格公式的方向。跨類型的評價法結合並不是不可行，但是需要第一階段各個評價法的評價績效驗證，讓較佳的評價模式彼此結合以產生資訊互補的效果。

Les polymères jouent un rôle important sur la réduction de la traînée et la modification de la structure des écoulements. Nous avons utilisé le modèle Oldroyd-B pour étudier l’effet de la viscoélasticité de solutions de polymères dilués sur des écoulements en deux dimensions autour d’obstacles dans un canal. Les obstacles sont pris en compte par la méthode de pénalisation volumique et des condition aux limites artificielles sans réflexion sont imposées à la sortie du canal. La discrétisation est effectuée par des schémas performants en différences finies et la résolution par une méthode multigrille. Les simulations numériques sont effectuées pour une large gamme de nombres de Reynolds et de nombres de Weissenberg. Les caractéristiques détaillées des écoulements viscoélastiques sont analysées et comparées entre elles et celles de l’écoulement du fluide sans polymère. En particulier les jeux de paramètres conduisant à une augmentation ou à une baisse de la traînée. Enfin les effets du polymère sur des écoulements turbulents sont aussi analysés. / Polymer plays an important role on the drag reduction and modification of the structure of flow. In this thesis, Oldroyd B model is employed to study the effectof viscoelasticity for polymer solutions diluted by two-dimensional direct numerical simulation. The obstacles are taken into account by penalization method. The artificial boundary condition is imposed without any reflection on the channel outlet. Flow pasta cylinder is investigated in detailed by present numerical methods. The numerical codes are valid for predicting the drag force and capturing the important character of viscoelastic flow by comparing with experimental and other numerical results. The drag map of the cylinder is obtained at a wide range of Reynolds number and Weissenberg number space. The detailed characteristics of viscoelastic flow are reported in the thesis. The effects of polymer on two-dimensional turbulent flow are also discussed by grid turbulent flow.