Dissertations / Theses on the topic 'Hele-Shaw flow; Stokes flow'

The Mullins-Sekerka problem, also called two-sided Hele-Shaw flow, arises in modeling a binary material with two stable concentration phases. A coarsening process occurs, and large particles grow while smaller particles eventually dissolve. Single particles become spherical. This process is described by evolving harmonic functions within the two phases with the moving interface driven by the jump in the normal derivatives of the harmonic functions at the interface. The harmonic functions are continuous across the interface, taking on values equal to the mean curvature of the interface. This dissertation reformulates the three-dimensional problem as one on the two-dimensional interface by using boundary integrals. A semi-implicit scheme to solve the free boundary problem numerically is implemented. Numerical analysis tasks include discretizing surfaces, overcoming node bunching, and dealing with topology change in a toroidal particle. A particle (node)-cluster technique is developed with the aim of alleviating excessive run time caused by filling the dense matrix used in solving a system of linear equations.

La microfluidique, domaine de recherche qui a émergé il y a juste 20 ans, a permis de réduire les dimensions des dispositifs d’analyse biologique ouvrant la porte au concept de « laboratoire-sur-puce » (lab-on-chip). Les succès de cette approche sont déjà nombreux, depuis l’analyse du génome en passant par la réduction du coût des analyses médicales. L’utilisation de gouttes comme enceinte réactionnelle au sein de ces dispositifs est une évolution récente qui permet de réduire encore le volume des échantillons biologiques et d’augmenter la vitesse d’analyse en parallélisant les mesures.Notre équipe développe des capteurs acoustiques dédiés à la détection d’analytes biologiques en milieu liquide. Ce type de capteur a pour principal défaut de ne permettre qu’une mesure contraignant au remplacement de l’interface de biodétection pour une réutilisation éventuelle du capteur. Dès lors, ils utilisent majoritairement une chambre de détection tout ou partie jetable, même si quelques travaux de recherche ont pu montrer la régénération d’un capteur par traitement chimique.Nous proposons ici de s’affranchir des étapes lourdes de remplacement ou de traitement de l’interface de détection qui conduisent entre chaque mesure au démontage du dispositif de détection. Nous employons dans ce cas les gouttes non plus comme enceinte réactionnelle mais comme interface de détection mobile. Elles ont ainsi le potentiel d’être générées et fonctionnalisées directement dans le dispositif pour détecter un analyte spécifique et peuvent être simplement évacuées afin de régénérer l’interface pour effectuer une mesure différente.Les travaux présentés dans cette thèse visent plus particulièrement la capture sur gouttes fonctionnalisées dans ce type de capteur innovant. Ils exposent le développement, incluant la fabrication et la caractérisation, de ces dispositifs microfluidiques ainsi que le montage d’un banc de test expérimental dédié. Ce sujet est suivi de deux projets ancillaires de développement de dispositifs microfluidiques liés aux capteurs acoustiques et à l’utilisation de gouttes. Le premier vise à homogénéiser les vitesses d’écoulements dans une chambre réactionnelle. Le second exploite les propriétés de génération de gouttes pour réaliser un condensateur à capacité variable
Since two decades the research on microfluidics systems allowed creating devices for biological detection with regular improvement in compactness, functionality integration and quantity of biological sample, leading to the concept of lab-on-chip. This approach has resulted in dramatic changes in the biomedical field, for example, opening the possibility to perform genomic analysis or improving the medical analysis cost. Using droplet as reaction chamber is a recent evolution that leads to a decrease in biological sample volume and an increase in analysis speed by multiplexing.Our team develops acoustical sensors dedicated to detect biomarker of interest in liquids. The principal weakness of theses sensors lies in their need for replacement of the biodetection interface for performing a new measurement. Accordingly, they use a detection chamber partially or totally disposable. However, few research works showed reusability of sensor by regenerating the bioreceptor layer on the detection interface by chemical treatment.We are proposing to avoid the replacement or the chemical treatment of the detection interface that requires dismounting the device between measurements. We are using here droplets, not as reaction chambers but as movable detection interface. They can be generated and configured directly inside the device to detect a specific biomarker. Then, droplets can be easily evacuated and replaced through the device, which allows to chain measurement of various configurations without dismounting it.The research work conducted in this thesis focuses on the fluidic aspects of this innovative sensor. They show development, including realization and characterization, of theses microfluidic devices and its dedicated characterization setup. This project is followed by two ancillary works about development of microfluidic devices for acoustical sensors and droplets systems. The first one is aiming at the homogenization of the flow velocity inside a reaction chamber. The second one is exploiting property of droplet generation for the realization of a variable capacitance capacitor

La digitation visqueuse est une instabilité hydrodynamique apparaissant lorsque, dans un milieu poreux, un fluide moins visqueux déplace un fluide plus visqueux. L'objectif de notre thèse est l'étude expérimentale des propriétés des motifs de digitation lorsque l'échantillon de fluide visqueux est de taille finie et lorsqu'une réaction chimique modifie la viscosité dans un milieu poreux modèle, en l'occurrence une cellule de Hele-Shaw. En particulier, notre étude a permis de quantifier la contribution de dispersion et de la digitation visqueuse, l'étalement dans l'espace d'échantillons de taille finie en fonction des paramètres expérimentaux (contraste de viscosité, vitesse de déplacement et taille de l'échantillon). Pour les fluides réactifs, nous analysons la digitation induite par une réaction A + B C dont le produit C est plus visqueux que les réactifs A et B, ceux-ci ayant la même viscosité. Nous mettons en évidence l'effet des concentrations en réactifs, du choix du fluide vecteur et du débit d'injection sur le motif de digitation.
Doctorat en Sciences
info:eu-repo/semantics/nonPublished

Submitted by Fabio Sobreira Campos da Costa (fabio.sobreira@ufpe.br) on 2016-12-12T14:44:32Z No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Dissertação - Pedro Henrique Amorim Anjos.pdf: 9970150 bytes, checksum: 29bca81f576fbdecd909c0ab24d68ce5 (MD5)
Made available in DSpace on 2016-12-12T14:44:32Z (GMT). No. of bitstreams: 2 license_rdf: 1232 bytes, checksum: 66e71c371cc565284e70f40736c94386 (MD5) Dissertação - Pedro Henrique Amorim Anjos.pdf: 9970150 bytes, checksum: 29bca81f576fbdecd909c0ab24d68ce5 (MD5) Previous issue date: 2015-02-05
CNPQ
Consideramos o movimento da interface entre dois fluidos viscosos, imiscíveis e incompressíveis confinados em uma célula de Hele-Shaw radial. Quando o fluido de baixa viscosidade é injetado no centro da célula e desloca radialmente o fluido de alta viscosidade, a instabilidade de Saffman-Taylor entra em ação levando ao crescimento de estruturas em formato de dedos. Dependendo da natureza dos fluidos envolvidos, eles podem molhar as paredes da célula de Hele-Shaw, deixando para trás um filme molhante de espessura finita. Diante disso, investigamos a influência da camada de fluido molhante, deixada pelo fluido deslocado, na dinâmica linear e fracamente não-linear do sistema. Mais especificamente, examinamos como a instabilidade da interface e os mecanismos de formação de padrões (bifurcação e competição de dedos) são afetados pela presença da camada de filme molhante no limite de baixo número de capilaridade. Estudamos também uma variação da instabilidade de Saffman-Taylor usual induzida pela injeção: o problema do levantamento em célula de Hele-Shaw. Nossos resultados analíticos indicam que o molhamento tem um impacto significativo nos padrões não-lineares resultantes. Ele é responsável por uma atenuação na variação do comprimento dos dedos, enquanto que induz o aparecimento de estruturas que apresentam dedos pequenos e largos, provenientes do fluido não molhante, alternados por dedos pequenos e finos pertencentes ao fluido molhante. Durante o processo de levantamento, medidas adquiridas pelo teste de adesão quantificam a força adesiva de um fluido viscoso confinado entre placas paralelas. Experimentos e intensas simulações numéricas indicam que a formação de dedos afeta a força de levantamento, promovendo uma diminuição na magnitude da mesma. Concluímos este trabalho propondo um modelo analítico que forneça a força de adesão considerando não só o efeito de dedos da interface, mas também a ação do molhamento e de estresses viscosos normais.
We consider the interfacial motion between two immiscible, incompressible, viscous fluids in the confined geometry of a radial Hele-Shaw cell. When the less viscous fluid is injected at the center and drives radially the more viscous fluid, the Saffman-Taylor instability takes place leading to the growth of fingerlike shapes. Depending on the nature of the fluids involved, they can wet the walls of the Hele-Shaw cell plates, leaving behind a film of finite thickness. In this framework, we investigate the influence of a thin wetting film trailing behind the displaced fluid on the linear and weakly nonlinear dynamics of the system. More specifically, we examine how the interface instability and the pattern formation mechanisms of finger tip-splitting and finger competition are affected by the presence of such a film in the low capillary number limit. We also examined a variant of the usual injection-driven Saffman-Taylor instability: the lifting Hele-Shaw cell problem. Our analytical results indicate that wettability has a significant impact on the resulting nonlinear patterns. It restrains finger length variability while inducing the development of structures presenting short, blunt penetrating fingers of the nonwetting fluid, alternated by short, sharp fingers of the wetting fluid. During the lifting process, probe-tack measurements evaluate the adhesion strength of viscous fluids confined between parallel plates. Existing meticulous experiments and intensive numerical simulations indicate that fingering formation affects the lifting force, making it to decrease in intensity. We conclude this work by proposing an analytical model that computes the lifting adhesion force by taking into account not only the effect of interfacial fingering, but also the action of wetting, and viscous normal stresses.

La propagation d'un front de réaction-diffusion est étudiée expérimentalement dans un écoulement cellulaire multiéchelle. Le front est produit par réaction autocatalytique en solution. L'écoulement est réalisé en géométrie de Hele-Shaw par électroconvection, son caractère multiéchelle étant réalisé par l'action combinée de deux nappes d'aimants d'échelles différentes. La géométrie du front et sa vitesse moyenne de propagation sont déterminées pour une large gamme d'intensité des vortex de chaque échelle. Elles sont confortées par une simulation numérique de l'avancée du domaine brulé dans le domaine frais. L'effet de la nature multiéchelle de l'écoulement sur la vitesse moyenne du front est compris par une méthode de renormalisation dont la validation est fournie par l'obtention d'un courbe maitresse pour l'ensemble des données
The propagation of a reaction-diffusion front is experimentally studied in a multi-scale cellular flow. The front is produced by an autocatalytic chemical reaction in an aqueous solution. The flow is generated by electroconvection and its multi-scale nature is induced by overlaying magnets of different scales. This enables an independent tune of the flow intensity at each scale. The geometry and the mean velocity of the front have been determined over a large range of scale intensities. These features are confirmed by a numerical simulation based on a burnt and fresh domain dynamics, the burnt domain expanding across the fresh one. The effect of the multi-scale nature of the flow on the mean front velocity is recovered by a renormalisation method validated by a collapse of the data onto a single curve

Integrated theoretical and experimental analysis is carried out in this work to investigate the fundamental failure mechanisms and flow patterns involved in the process of fluid injection into dense granular media. The experimental work is conducted with aqueous glycerin solutions, utilizing a novel setup based on a Hele-Shaw cell filled with dense dry sand. The two dimensional nature of the setup allows direct visualization and imaging analysis of the real-time fluid and grain kinematics. The experimental results reveal that the fluid flow patterns show a transition from simple radial flow to a ramified morphology while the granular media behaviors change from that of rigid porous media to localized failure that lead to development of fluid channels. Based on the failure/flow patterns, four distinct failure/flow regimes can be identified, namely, (i) a simple radial flow regime, (ii) an infiltration-dominated regime, (iii) a grain displacement-dominated regime, and (iv) a viscous fingering-dominated regime. These distinct failure/flow regimes emerge as a result of competition among various energy dissipation mechanisms, namely, viscous dissipation through infiltration, dissipation due to grain displacements, and viscous dissipation through flow in thin channels and can be classified based on the characteristic times associated with fluid injection, hydromechanical coupling and viscoelastoplasticity. The injection process is also analyzed numerically using the discrete element method (DEM) coupled with two fluid flow scheme, a fixed coarse grid scheme based on computational fluid dynamics (CFD) and a pore network modeling scheme. The numerical results from the two complementary methods reproduce phenomena consistent with the experimental observations and justify the concept of associating the displacement regimes with the partition among energy dissipation mechanisms. The research in this work, though fundamental in nature, will have direct impacts on many engineering problems in civil, environmental and petroleum engineering such as ground improvement, environmental remediation and reservoir stimulation.

Dans le chapitre 1, on considère le modèle micropolaire de Navier-Stokes avec conditions de bords de type Dirichlet non homogènes en dimension deux. On donnera un résultat d'existence d'une solution faible en utilisant le théorème du point fixe de Leray-Schauder, puis on prouvera l'unicité de la solution faible du problème sous certaines hypothèses. On établiera une justification mathématique de l’équation de Reynolds généralisé à partir de ce modèle là. On étudiera ensuite la forme de l'équation de Reynolds suivant le choix de la viscosité et des données initiales. Dans le chapitre 2, nous considérons le modèle de Hele-Shaw généralisé dans une cellule laminaire, qui consiste à injecter du fluide, avec un débit non constant w 0, à travers un trou de frontière 1, situé sur l'une des deux surfaces ; et à tenir compte que l'une des surfaces a une géométrie quelconque et animée d'un mouvement relatif vertical. En introduisant un changement de variable de type Baiocchi, le problème initial se ramène à l'étude d'une inéquation variationnelle avec terme de Volterra. L'existence d'une solution pour cette dernière est donnée par le théorème du point fixe de Banach. Des résultats de régularité en espace pour la solution seront prouvés en introduisant un problème pénalisé et en utilisant la méthode de Rothe (semi-discrétisation en temps), puis on montrera que la dérivée par rapport à t de la solution de l'inéquation variationnelle est dans l#(0, t, h#2()), ce dernier résultat nous permet de revenir au problème initial. Dans le chapitre 3, on considère un problème de Stefan à deux phases avec convection. Le problème est gouverné par un système couple non linéaire, comprenant la loi de Darcy pour un fluide non newtonien et l'équation d'équilibre d'énergie avec second membre dans l#1. Pour prouver l'existence de solutions du problème faible on introduira une famille de solutions approchées (#, p#), > 0, définies sur le domaine entier , en insérant une fonction de pénalité convenable dans l'équation de pression. On considère ensuite séparement les problèmes en # et p#, respectivement, et en utilisant le principe de point fixe de Schauder, on montre l'existence de couples solutions (#, p#) du problème approché, pour tout > 0. En faisant tendre vers zéro, on montre que les solutions du problème approché convergent vers une limite (, p) qui est une solution faible du problème variationnel. On montre aussi que la fonction est continue d'où le domaine où > 0 est un ensemble ouvert, et l'interface des deux phases est définie a posteriori comme l'ensemble de niveau = 0. On établira, enfin, quelques relations entre les solutions faibles et classiques, dans le cas d’une courbe assez régulière

This thesis is in the field of non-linear partial differential equations (PDE), focusing on problems which show some type of phase-transition. A single phase Hele-Shaw flow models a Newtoninan fluid which is being injected in the space between two narrowly separated parallel planes. The time evolution of the space that the fluid occupies can be modelled by a semi-linear PDE. This is a problem within the field of free boundary problems. In the multi-phase problem we consider the time-evolution of a system of phases which interact according to the principle that the joint boundary which emerges when two phases meet is fixed for all future times. The problem is handled by introducing a parameterized equation which is regularized and penalized. The penalization is non-local in time and tracks the history of the system, penalizing the joint support of two different phases in space-time. The main result in the first paper is the existence theory of a weak solution to the parameterized equations in a Bochner space using the implicit function theorem. The family of solutions to the parameterized problem is uniformly bounded allowing us to extract a weakly convergent subsequence for the case when the penalization tends to infinity. The second problem deals with a parameterized highly oscillatory quasi-linear elliptic equation in divergence form. As the regularization parameter tends to zero the equation gets a jump in the conductivity which occur at the level set of a locally periodic function, the obstacle. As the oscillations in the problem data increases the solution to the equation experiences high frequency jumps in the conductivity, resulting in the corresponding solutions showing an effective global behaviour. The global behavior is related to the so called homogenized solution. We show that the parameterized equation has a weak solution in a Sobolev space and derive bounds on the solutions used in the analysis for the case when the regularization is lost. Surprisingly, the limiting problem in this case includes an extra term describing the interaction between the solution and the obstacle, not appearing in the case when obstacle is the zero level-set. The oscillatory nature of the problem makes standard numerical algorithms computationally expensive, since the global domain needs to be resolved on the micro scale. We develop a multi scale method for this problem based on the heterogeneous multiscale method (HMM) framework and using a finite element (FE) approach to capture the macroscopic variations of the solutions at a significantly lower cost. We numerically investigate the effect of the obstacle on the homogenized solution, finding empirical proof that certain choices of obstacles make the limiting problem have a form structurally different from that of the parameterized problem.

QC 20180222

Cette thèse expérimentale s’inscrit dans le domaine des sciences séparatives et se base sur la technique de SPLITT (SPLIT-flow Thin fractionation). Son objectif consiste en l’étude des mécanismes qui sont à l’origine de la séparation, en continu et sans membrane, d’objets de taille micrométrique dans des mini-séparateurs fluidiques (Step-SPLITT). Les expériences menées, en laboratoire et lors de vols paraboliques, ont révélé le couplage complexe comme l’influence des effets hydrodynamiques et du champ gravitationnel sur la migration transverse des espèces en écoulement. Des visualisations tridimensionnelles par holographie digitale ont corroboré nos résultats et dévoilé des comportements inattendus. Les capacités séparatives des Step-SPLITT ont rendu possible l’analyse et la séparation d’objets biologiques et biomimétiques. Enfin, cette étude complétée par une modélisation tridimensionnelle de l’écoulement nous a permis de mettre au point un nouveau prototype de séparateur.

This experimental thesis belongs to the field of separative sciences and is based on the SPLITT technique (SPLIT-flow Thin fractionation). The objective is to study the mechanisms that are at the origin of continuous and membraneless separation of micron-size species in mini fluidic separators (Step-SPLITT). Experiments undertaken in laboratory and during parabolic flights revealed the complex coupling of the hydrodynamic effects and the gravitational field influencing the transverse migration of the flowing species. Three-dimensional visualizations performed by digital holography confirmed our results and disclosed unexpected behaviours. The separation capacities of Step-SPLITT made the analysis and the separation of biological and biomimetic species possible. In addition this study in conjunction with a three-dimensional flow modelling enabled us to develop a new prototype of separator.
Doctorat en sciences appliquées
info:eu-repo/semantics/nonPublished

This thesis discusses the radial version of the Hele-Shaw problem. Different from the channel version, traveling-wave solutions do not exist in this version. Under algebraic potentials, in the case that the droplets expand, in finite time, cusps will appear on the boundary and classical solutions may not exist afterwards. Physicists have suggested that for (2p+1,2)-cusps, that near cusp singularities of Hele-Shaw flow, after scaling X, Y by some powers of time t respectively, the main part of Y(X, t) is a one-parameter family and does not depend on time t. They have also suggested that the solutions of the Hele-Shaw problem are connected with dispersionless KdV (dKdV) hierarchy. In this study, we rigorously proved that this is the case for (3,2)-cusps when the droplets are simply connected and the external potentials are algebraic. We gave exact solutions and showed that the main parts of the exact solutions are some special solutions of the dispersionless string equation. More over, borrowed from the physical paper$\cite{Teo}$ with a little more details, we showed the arguments of how these special solutions are related to dKdV hierarchy.

碩士
國立交通大學
機械工程系所
106
The LOCA coating machine uses stress to force two plates where OCR material is placed between close-fitting. The glue optimal pattern spreads close to its base plate and no glue should excess the plate boundary. The processes mimic the problems of compressing Hele-Shaw flows. Accordingly, we assume that the inner flow pattern formed from the pulled plate is same as the original figure which upper plate pushed down. In this study, we use computer simulation method to find out that how the viscosity ratio, the initial aspect ratio, the lifting plate velocity and the surface tension induce the optimization of glue pattern. As the result of interfacial instability, less viscous fluid displaces another fluid of higher viscosity while the height increasing. Four sides of the rectangle will rapidly shrink inward. On the contrary, the pattern will appear four fingers on the corner of the rectangle.

博士
國立中央大學
土木工程研究所
97
This study explores the 2-D transient seepage flow and steady flow from cracks into different porous media, including both coarse material and fine material. The Crack Apparatus of Hele-Shaw Model is chosen to analyze the configuration of wetting front in the coarse material. Due to the effect of gravity on the flow condition, the wetting front in the coarse material is more elongated than that of the fine material. The limited width of wetting front is approaching a constant value in the coarse material. On the other hand, the effect of capillary causes the wetting front in the fine sand expands continuously as time increases. In this paper, a Green-Ampt type of conceptual model is proposed to describe the propagating process of the wetting front in different materials. The time-dependent wetting front and the corresponding wetted area are obtained in this study, which are beneficial for examining the leaky contaminant transport from cracks into different porous media or for the analysis of the slope stability problems. The second topic of this study is development and validation of an analytical solution for the flow through horizontal and sloping bed. The parameter, B, shown in the analytical solution of horizontal bed increases with the increasing ratio of width to upstream level. The Hele–Shaw model is employed to examine the phreatic surface and exit point of seepage flow by adjusting gap width and fluids density. Aluminum powders were used to as tracking particles inside the gap and the flow visualization was performed. Both horizontal and sloping bed, the phreatic (free) surfaces and exit points are in good agreement between experimental data and the analytical solution. The assumption of a parabolic piezometric head was verified by experimental data. In the horizontal bed, the ratio of down-upstream water depth is large than 0.9 and the ratio of upstream depth -width is less than 0.1, the solution based on Dupuit model is close to the analytical solution.

碩士
國立成功大學
航空太空工程學系碩博士班
100
In this study, the interfacial instabilities of miscible magnetic fluids and diesel were theoretically and experimentally studied in a rotating Hele-Shaw cell. First, the dimensional analysis was used to understand the relation magnetic force, viscous force, centrifugal force, pressure force and Coriolis force. In addition, it could assist in experimental design. Two kinds of the flow field were studied, the open Hele-Shaw cell and the closed Hele-Shaw cell. The labyrinthine fingering phenomena of the miscible magnetic fluid interface on a perpendicular magnetic field were investigated. First, the interface of miscible magnetic fluids and diesel became more unstable when the rotating speed and the consequent centrifugal force were increased without magnetic field applied. Three ferrofluids were used in this study. To compare the results, the gyration radius will be grown fast with ferrofluid EMG905. The first reason was that the initial susceptibility and particle concentration for EMG905 were higher than that for the other two ferrofluids. So the diffusion and magnetic effects for EMG905 will be the most obvious one. The second reason was that the density difference between EMG905 and the diesel was the highest, so the Rayleigh-Taylor instability will increase for the centrifugal force. In the first part, the subject focused on the labyrinthine fingering phenomena of the interface in an open Hele-Shaw cell. Experimental results showed that the gyration radius will be grown fast with increasing the rotating speed in the same magnetic field. Under the same rotating speed condition, the gyration radius will be grown fast with higher magnetic field in the first stage. And in the later stage, there will be some differences between the developments of gyration radius with different magnetic fields. For lower magnetic field, the growth will be demonstrated by the centrifugal force. For higher magnetic field, the viscosity will be increased because of the perpendicular magnetic field, and the viscous force will stabilize the interface. The second part of this study focused on the experiment in a closed Hele-Shaw cell. For the closed flow field, it can be regard as a solid-body rotation and there will produce an opposite pressure force, which can help to stabilize the interface. Experimental results showed that the growth of gyration radius will be fast in an open Hele-Shaw cell. By doing so, the appearance of the opposite pressure force can be proved. At last, the experimental results also showed that the growth of gyration radius will be fast in an open Hele-Shaw cell with higher gap. So, the three-dimensional effect played an important role in miscible magnetic fluids and diesel.

We study experimentally the gravity-driven bubbly flow inside a large Hele-Shaw cell. The bubbles and foam were created by a series of upside-down overturns of the half-filled cell about its horizontal axis. When the liquid flows down it entraps a large number of bubbles, which remain stable as the liquid contains surfactant molecules. The total number and sizes of these bubbles slowly asymptote to a steady state after dozens of overturns. It takes longer to reach this asymptote when the viscosity of the liquid is larger. The bubbles also become more monodisperse with more cell over-turns. The number and distribution of the bubbles in turn affects the average motion of the liquid phase, which is characterized by the downwards motion of the liquid center of mass. We use high-resolution 6k video-camera to track the trajectories of thousands of bubbles. This required the development of software codes to identify individual bubbles and follow them between video frames. Successful thresholding algorithm required a machine-learning component, which was integrated into the program. This program also needed to account for possible splitting or coalescence of adjacent bubbles. The program can also find the velocities along the trajectories. In this way we can find the vertical velocity of bubbles as a function of their sizes. The smaller bubbles are sometimes observed to move downwards against their buoyancy. This occurs when the viscous stress from the surrounding liquid phase overcomes the upwards buoyancy force. Bubbles with similar sizes were often found to be stacking together and having worm-like rising movement that is faster than their individual rising velocity. The occurrence of the bubble stacking was dependent on the distance between the bubbles, their sizes and their wakes. Clusters of tiny bubbles that are much smaller than the gap of the Hele-Shaw cell were observed to form layers which can severely hinder the overall liquid motion.

博士
國立交通大學
機械工程系所
102
Viscous fingering is an interfacial fluid flow instability that occurs when less viscous fluid displaces another more viscous one in a Hele-Shaw cell or porous media, leading to the formation of finger-like pattern at the interface of both fluids. The interfacial evolution of multiphase flows will severely impact on the quality of production and efficiency in a variety of practical application of industrial process. Most frequent example of this instability is that of oil recovery for which viscous fingering takes place when an aqueous solution displaces more viscous oil in underground reservoirs, leading to the formation of nontrivial fingerlike structure and reduce the efficiency of the displacement process. Another particularly interesting variation of the classic radial flow is the investigation of fingering instabilities in Hele-Shaw cells presenting variable gap spacing. This is also a very important issue in many industrial areas including adhesion, lubrication, and colloidal hydrodynamics. In this dissertation, we carried out the highly accurate simulation to investigate the interfacial evolution in two scenarios－radial injection-driven miscible flow and lifting radial Hele-Shaw flow, both with the monotonic and nonmonotonic viscosity profile. So, the thesis consists of two parts: Part 1 focus on radial injection-driven miscible flow in a Hele-Shaw cell and covers three major topics. To begin with, we perform numerical experiments in a wide range to study the dispersion relation on both the Péclet number and the parameters of the viscosity profile. A monotonic viscosity-concentration relation of exponential type (concave) by other scholars is assumed, and a linear and reverse (convex) monotonic viscosity profiles and nonmonotonic one are also discussed. Results of this study show that as the overall viscosity contrast held constant, nonmonotonic viscosity profile lead to a more stable flow than that of monotonic one, and there are no significant differences in different viscosity profiles. However, if the nonmonotonic viscosity profile crosses the convex monotonic viscosity profile, the nonmonotonic feature enhances the prominence of interfacial instability. Then, a great variety of morphological behaviors is systematically introduced. In general, the nonmonotonic feature enhances the prominence of interfacial instability. Formation of dual vortex pairs and “reverse fingering”, where the fingers spread farther in the backward than in the forward direction are observed, which are not present in monotonic viscosity profile. Finally, we have carried out a parameter study to understand the effects of nonmonotonicity on the stability of the injection flow. In part 2, discussions start with the investigation of the influence of lifting scenario and the perturbation set. Contrast to the injection-driven miscible flow in radial Hele-Shaw cells which leads to the formation of morphing flow phenomenon of finger tip-splitting and side-branch events are plentiful if the injection rate is constant with time. More complicated flow are present for time-dependent gap flow which results in different kinds of patterns, and leads to intricate morphologies if the cell’s gap width grows exponentially with time. Recent studies show that the growing of intricate patterns due to lifting can be ntrolled by properly adjusting the time-dependent gap width. Moreover, we found the exponential lifting case will cause the flow more unstable than the variant lifting situation. We also deduce higher Péclet number and viscous contrast (A in monotonic viscosity profile and μm in nonmonotonic one) demonstrate more vigorous fingering. The sensitivity of the system to changes in the initial conditions and perturbation set is also discussed. Next, the effects of four viscosity profiles as stated in part 1 have been investigated. Unlike injection flow,the stability of three monotonic viscosity profiles are always in the series of concave, linear and convex. However, as injection flow, if the nonmonotonic viscosity profile crosses the convex curve will enhances the prominence of interfacial instability. Finally, we have carried out a parameter study to understand the effects of nonmonotonicity viscosity profile on the stability of the lifting flow.