[pt] Soluções poliméricas diluídas podem ser usadas na indústria de Óleo e Gás para injeção em reservatórios de petróleo como agentes de controle de mobilidade em técnicas de recuperação avançada de óleo, com o objetivo de aumentar o fator de recuperação do campo e postergar a indesejada alta taxa de produção de água. No entanto, o fenômeno macroscópico esperado, aumento da eficiência de varrido, é influenciado diretamente por fatores microscópicos que ocorrem na escala de poro. Devido às suas grandes cadeias moleculares, ao avançar através do meio poroso, soluções poliméricas de alto peso molecular dificultam o avanço da água alterando a razão de mobilidades água/óleo. Porém, o estudo de escoamento de soluções poliméricas ainda é um desafio devido ao comportamento não Newtoniano da relação vazão x diferença de pressão e à sensibilidade a diversos fatores como: tipo de polímero utilizado, concentração, razão de viscosidades, temperatura, salinidade e taxa de cisalhamento. O presente trabalho analisa a implementação de dois modelos constitutivos, Oldroyd-B e FENE-CR, aplicados ao método de solução de elementos finitos, EF, utilizando a formulação dos resíduos ponderados de Galerkin para modelar o escoamento de soluções poliméricas diluídas através de capilares com constrição e entender os fenômenos viscosos e elásticos envolvidos. Os resultados obtidos fornecem uma descrição mais detalhada da aplicação dos modelos testados e do escoamento de soluções poliméricas diluídas em gargantas de poros. / [en] Diluted polymeric solutions can be used in Oil & Gas industry for reservoir injection as mobility control agents in enhanced oil recovery technics, which goal is to increase the oil field recovery factor and postpone the high water production rates. However, the macroscopic expected phenomena, increase of the areal sweep, is directly influenced by microscopic factors that occurs at the pore scale. Due to their huge molecular chains, by advancing through the porous media, high molecular weight polymer solutions hinder the water flow changing the water/oil mobility ratio. Nevertheless, the comprehension of polymeric solutions flow is still a great challenge due to the Non-Newtonian behavior of the flow rate x pressure gradient relation and the sensibility to several variables such as: type of polymer used, concentration, viscosity ratio, temperature, salinity and shear rate. The present work analyzes the implementation of two constitutive models, Oldroyd-B and FENE-CR, coupled with the Finite Element method using the Galerkin weighted residual formulation to model the flow of diluted polymeric solutions through constricted capillaries and understand the viscous and elastic phenomena involved. The obtained results provide a more detailed portrait of the constitutive models and the flow of diluted polymeric solutions through pore throats.

本研究では懸濁系の磁気粘弾性流体の流動特性と圧力特性について四要素モデルを用いて考察を行う．実験では，磁気粘弾性流体がオリフィスを流れる際の圧力損失について，スタートアップの圧力変化と磁場印加時の圧力の過渡的な挙動を測定した．また実験結果の検証として数値解析も行った．その結果，磁気粘弾性流体の流れは，ケルビン力，磁気スピンの作用，物性を決定づける4要素モデルによって数値的に表現されることを明らかにした． / The flow and pressure characteristics of magneto-viscoelastic fluids in suspension systems are discussed in this study using a four-element model. In the experiment, the pressure drop of the magneto-viscoelastic fluid flowing through the orifice was measured in terms of pressure change at startup and transient pressure behavior when a magnetic field is applied. Numerical analysis was also performed to verify the experimental results. The results show that the flow of a magneto-viscoelastic fluid can be numerically described by a four-element model that determines the Kelvin force, the action of magnetic spin, and physical properties. / 博士(工学) / Doctor of Philosophy in Engineering / 同志社大学 / Doshisha University

[pt] Filmes finos de líquido estão presentes em uma variedade de sistemas e aplicações. Estamos interessados em filmes compostos por duas camadas, que são comuns no processo de revestimento por cortina. No revestimento por cortina, o líquido cai de uma matriz formando uma cortina formada por um filme fino antes de molhar o substrato em movimento. Um dos limites mais importantes do processo é a ruptura da cortina, que define um limite inferior para a vazão do líquido de revestimento. Consequentemente, este limite inferior da vazão define a espessura mínima viável do filme depositado. Evidências experimentais mostraram que o uso de uma cortina compostas por duas camadas, com uma das camadas sendo mais fina e viscoelástica, pode atrasar a ruptura da cortina para taxas de fluxo mais baixas. A quebra de filmes líquidos de duas camadas, compostas por um líquido newtoniano e um viscoelástico, é estudado por meio da resolução das equações diferenciais que descrevem a evolução da configuração do filme até seu rompimento. O efeito de diferentes parâmetros no tempo de ruptura é determinado. Os resultados mostram o mesmo comportamento observado experimentalmente, a fina camada de líquido viscoelástico retarda o rompimento, estabilizando o filme líquido. / [en] Thin liquid sheets are present in a variety of systems and applications. Here, we are interested in double-layered sheets, which are common in the curtain coating process. In curtain coating, the liquid falls from a die forming a thin curtain before wetting the moving substrate. One of the most important process limits is the curtain breakup, which sets a lower limit for the coating liquid flow rate. Consequently, this flow rate lower limit defines the feasible minimum deposited film thickness. Experimental evidence have shown that using a two-layer curtain, with a viscoelastic thin layer, may delay the curtain breakup to lower flow ratios. The breakup of two-layer liquid sheets, composed of a Newtonian and a viscoelastic liquid, is studied by solving the differential equations that describe the evolution of the liquid sheet configuration until breakup. The effect of different parameters on the breakup time is determined. The results show the same behavior observed experimentally, thin viscoelastic liquid layer delays the breakup, stabilizing the liquid sheet.

Surface coating is widely used in microelectronic industry to produce thin films over surfaces with uneven topography. Such processes are used in fabricating integrated circuits, storage devices, such as magnetic disks, memory devices and optical disks as well as for manufacturing adhesives, magnetic tapes, magazines which can produce thicker films over patterns of similar depth and width at higher speeds. Other applications of film flow over uneven surfaces come from specific designs of surfaces of heat-exchangers and the surfaces of various structured packings used to improve heat and mass transfer operations. The one-dimensional, gravity-driven film-flow of a linear or exponential PTT liquid, flowing either on the outer or on the inner surface of a vertical cylinder or over a planar wall is analyzed. Numerical solution of the governing equations is generally possible. Analytical solutions are derived only for: (i) linear PTT model in cylindrical and planar geometries in the absence of solvent and the affinity parameter set at zero; (ii) linear or exponential PTT model in a planar geometry in the absence of solvent and the affinity parameter the affinity parameter obtains nonzero values; (iii) exponential PTT model in planar geometry in the absence of solvent and the affinity parameter set at zero. Then, the two-dimensional, steady flow of a viscoelastic film over a periodic topography under the action of a body force is studied. It is examined the interplay of elastic, viscous, inertia and capillary forces on the film thickness and planarization efficiency over steep topographical changes of the substrate. The code is validated by verifying that in isolated topographies the periodicity conditions result in fully developed viscoelastic film flow at the inflow/outflow boundaries and that its predictions for Newtonian fluids over 2D topography under creeping flow conditions coincide with those of previous works. Finally, the steady film-flow of a Newtonian fluid has been studied over a trench examining the various types of inclusions that can be formed. It can be distinguished three possible flow configurations when (a) the triple contact points are ‘pinned’ at the lips of the cavity, (b) the triple contact points are at the left side and the bottom of the cavity so that the cavity is not filled with liquid only around its left concave corner and (c) the two triple contact points are at the two sides of the cavity so that its bottom remains empty. / Η μελέτη ρευστών υμένων κατέχει σημαντική θέση στη σύγχρονη επιστήμη και τεχνολογία και συναντάται ευρύτατα σε βιομηχανικές διεργασίες (π.χ. κατασκευή ηλεκτρονικών εξαρτημάτων όπως είναι μικροεπεξεργαστές, ολοκληρωμένα κυκλώματα και συσκευές μνήμης, κατασκευή ψηφιακών αποθηκευτικών μέσων όπως μαγνητικοί δίσκοι, δισκέτες, κασέτες και οπτικοί δίσκοι, διεργασίες φωτολιθογραφικής επικάλυψης και εκτύπωσης όπως επικαλυπτικές ουσίες, μαγνητικές ταινίες και περιοδικά κ.α.). Επιπλέον, οι ρευστοί υμένες χρησιμοποιούνται για την βελτιστοποίηση διεργασιών μεταφοράς θερμότητας και μάζας (εναλλάκτες θερμότητας). Πρώτα αναλύεται η μονοδιάστατη ροή λόγω βαρύτητας ενός υμένα ιξωδοελαστικού ρευστού που ακολουθεί το καταστατικό μοντέλο Phan-Thien and Tanner (PTT) σε γραμμική ή εκθετική μορφή. Το ρευστό ρέει είτε στην εξωτερική είτε στην εσωτερική επιφάνεια ενός κατακόρυφου κυλίνδρου ή σε επίπεδο τοίχωμα. Η αριθμητική επίλυση των διεπουσών εξισώσεων είναι πάντοτε δυνατή. Αναλυτικές εκφράσεις εξάγονται μόνο για (i) γραμμικό μοντέλο PTT σε κυλινδρικές και επίπεδες γεωμετρίες απουσία διαλύτη και η παράμετρος συγγένειας είναι μηδενική, (ii) γραμμικό ή εκθετικό μοντέλο PTT σε επίπεδη γεωμετρία απουσία διαλύτη με μη μηδενική παράμετρο συγγένειας και (iii) εκθετικό μοντέλο PTT σε επίπεδη γεωμετρία απουσία διαλύτη και με μηδενική παράμετρο συγγένειας. Στη συνέχεια μελετάται η δισδιάστατη μόνιμη ροή ενός ιξωδοελαστικού υμένα κατά μήκος περιοδικής τοπογραφίας υπό την επίδραση είτε της βαρύτητας είτε της φυγοκέντρου δύναμης. Εξετάζεται τόσο η επίδραση των ιξωδοελαστικών ιδιοτήτων υπό συνθήκες έρπουσας ροής όσο και των τριχοειδών και αδρανειακών δυνάμεων καθώς και της γεωμετρίας του υποστρώματος στο πάχος του υμένα και στη δυνατότητα εξομάλυνσής του παρά τις απότομες αλλαγές του υποστρώματος. Ο κώδικας ελέγχθηκε μέσω της σύγκρισης των προβλέψεών του με προηγούμενες εργασίες που αφορούν σε Νευτωνικά και σε ιξωδοελαστικά ρευστά. Τέλος, εξετάζεται η δισδιάστατη μόνιμη ροή ενός Νευτωνικού υμένα πάνω από παρόμοια μεταβαλλόμενο υπόστρωμα, ενώ επιτρέπεται να σχηματίζονται εγκλείσματα αέρα μεταξύ του υμένα και του υποστρώματος. Διακρίνονται οι εξής δυνατότητες σχηματισμού εγκλεισμάτων αέρα ανάλογα με τη θέση της κάτω επιφάνειας του υμένα: (α) εμφανίζονται σταθερά σημεία επαφής στις άνω (κυρτές) γωνίες του υποστρώματος, (β) σχηματίζεται έγκλεισμα γύρω από την πρώτη κοίλη γωνία του υποστρώματος και (γ) εμφανίζονται δύο τριπλά σημεία επαφής με τις παράπλευρες επιφάνειες του υποστρώματος.

<p>Particle suspensions in viscoelastic fluids (e.g., polymeric fluids, liquid crystalline solutions, gels) are ubiquitous in industrial processes and in biology. In such fluids, particles often acquire lift forces that push them to preferential streamlines in the flow domain. This lift force depends greatly on the fluid’s rheology, and plays a vital role in many applications such as particle separations in microfluidic devices, particle rinsing on silicon wafers, and particle resuspension in enhanced oil recovery. Previous studies have provided understanding on how fluid rheology affects the motion of spherical particles in simple viscoelastic fluid flows such as shear flows. However, the combined effect of more complex flow profiles and particle shape is still under-explored. The main contribution of this thesis is to: (a) provide understanding on the migration and rotation dynamics of an arbitrary-shaped particle in complex flows of a viscoelastic fluid, and (b) develop guidelines for designing such suspensions for general applications.</p> <p><br></p> <p>In the first part of the thesis, we develop theories based on the second-order fluid (SOF) constitutive model to provide solutions for the polymeric force and torque on an arbitrary-shaped solid particle under a general quadratic flow field. When the first and second normal stress coefficients satisfy  <strong>Ψ</strong>₁  = −2 <strong>Ψ</strong> ₂ (corotational limit), the fluid viscoelasticity modifies only the fluid pressure and we provide exact solutions to the polymer force and torque on the particle. For a general SOF with  <strong>Ψ</strong> ₁ ≠  −2 <strong>Ψ</strong> ₂, fluid viscoelasticity modifies the shear stresses, and we provide a procedure for numerical solutions. General scaling laws are also identified to quantify the polymeric lift based on different particle shapes and orientation. We find that the particle migration speed is directly proportional to the length the particle spans in the shear gradient direction (L_sg), and that polymeric torques lead to unique orientation behavior under flow.</p> <p><br></p> <p>Secondly, we investigate the migration and rotational behavior of prolate and oblate spheroids in various viscoelastic, pressure-driven flows. In a 2-D slit flow, fluid viscoelasticity causes prolate particles to transition to a log-rolling motion where the particles orient perpendicular to the flow-flow gradient plane. This behavior leads to a slower overall migration speed (i.e., lift) of prolate particles towards the flow centerline compared to spherical particles of the same volume. In a circular tube flow, prolate particles align their long axis along the flow direction due to the extra polymer torque generated by the velocity curvature in all radial directions. Again, this effect causes prolate particles to migrate slower to the flow centerline than spheres of the same volume. For oblate particles, we quantify their long-time orientation and find that they migrate slower than spheres of the same volume, but exhibit larger migration speeds than prolate particles. Lastly, we examine the effect of normal stress ratio ? <strong>α</strong>  = <strong>Ψ</strong> ₂ /<strong>Ψ</strong>₁on the particle motion and find that this parameter only quantitatively impacts the particle migration velocity but has negligible effect on the rotational dynamics. We therefore can utilize the exact solution derived under the corotational limit (?<strong>α</strong> = −1/2) for a quick and reasonable prediction on the particle dynamics.</p> <p><br></p> <p>We next experimentally investigate the migration behavior of spheroidal particles in microfluidic systems and draw comparisons to our theoretical predictions. A dilute suspension of prolate/oblate microparticles in a density-matched 8% aqueous polyvinylpyrrolidone (PVP) solution is used as the model suspension system. Using brightfield microscopy, we qualitatively confirm our theoretical predictions for flow Deborah numbers 0 < De < 0.1 – i.e., that spherical particles show faster migration speed than prolate and oblate particles of the same volume in tube flows.</p> <p><br></p> <p>We finally design a holographic imaging method to capture the 3-D position and orientation of dynamic microparticles in microfluidic flow. We adopt in-line holography setup and propose a straightforward hologram reconstruction method to extract the 3-D position and orientation of a non-spherical particle. The method utilizes image moment to locate the particle and localize the detection region. We detect the particle position in the depth direction by quantifying the image sharpness at different depth position, and uses principal component analysis (PCA) to detect the orientation of the particle. For a semi-transparent particle that produces complex diffraction patterns, a mask based on the image moment information can be utilized during the image sharpness process to better resolve the particle position.</p> <p><br></p> <p>In the last part of this thesis, we conclude our work and discuss the future research perspectives. We also comment on the possible application of current work to various fields of research and industrial processes.</p> <p><br></p>