Developing discontinuous Galerkin methods for the resolution of the incompressible Navier-Stokes equations

Informe de investigación elaborado a partir de una estancia en el Laboratorio de Diseño Computacional en Aeroespacial en el Massachusetts Institute of Technology (MIT), Estados Unidos, entre noviembre de 2006 y agosto de 2007. La aerodinámica es una rama de la dinámica de fluidos referida al estudio de los movimientos de los líquidos o gases, cuya meta principal es predecir las fuerzas aerodinámicas en un avión o cualquier tipo de vehículo, incluyendo los automóviles. Las ecuaciones de Navier-Stokes representan un estado dinámico del equilibrio de las fuerzas que actúan en cualquier región dada del fluido. Son uno de los sistemas de ecuaciones más útiles porque describen la física de una gran cantidad de fenómenos como corrientes del océano, flujos alrededor de una superficie de sustentación, etc. En el contexto de una tesis doctoral, se está estudiando un flujo viscoso e incompresible, solucionando las ecuaciones de Navier- Stokes incompresibles de una manera eficiente. Durante la estancia en el MIT, se ha utilizado un método de Galerkin discontinuo para solucionar las ecuaciones de Navier-Stokes incompresibles usando, o bien un parámetro de penalti para asegurar la continuidad de los flujos entre elementos, o bien un método de Galerkin discontinuo compacto. Ambos métodos han dado buenos resultados y varios ejemplos numéricos se han simulado para validar el buen comportamiento de los métodos desarrollados. También se han estudiado elementos particulares, los elementos de Raviart y Thomas, que se podrían utilizar en una formulación mixta para obtener un algoritmo eficiente para solucionar problemas numéricos complejos.

Estudi, disseny i construcció d'una sonda submarina

Treball de recerca realitzat per un alumne d’ensenyament secundari i guardonat amb un Premi CIRIT per fomentar l'esperit científic del Jovent l’any 2008. S’ha construït una sonda submarina a partir dels següents requeriments: havia de ser senzilla i barata, amb mobilitat en totes les direccions i hauria de poder visualitzar el medi per on es desplaçaria. Un cop construïda la sonda es formularen tres hipòtesis relatives a la resistència que oposa l’aigua al desplaçament de la sonda en funció de la profunditat, la temperatura i la viscositat de l’aigua. Després d’haver observat el comportament dels cossos submergits en un fluid, es va dissenyar la sonda submarina d’acord amb el principi d’Arquimedes. Per construir i dissenyar els circuits elèctrics també es va experimentar amb diferents possibilitats; finalment es va optar per un sistema de relès que controlés individualment cada receptor del circuit, motors i càmera submarina. Per contrastar les hipòtesis, mitjançant un dinamòmetre, es van prendre les mesures de la força de resistència de l’aigua al desplaçament de la sonda. La construcció de la sonda ha estat un èxit, ja que és perfectament estanca i té una bona maniobrabilitat, els motors permeten el moviment de la sonda en qualsevol direcció de l’espai i el senyal de la càmera submarina és prou bo per distingir els objectes que passin pel davant de la sonda i pel fons marí. Dues de les hipòtesis han estat contrastades; la tercera no: s’ha comprovat que la força de resistència augmenta amb la profunditat i amb la disminució de temperatura de l’aigua que envolta a la sonda; la hipòtesi de la viscositat no s’ha pogut contrastar.

Experimental investigation of transient natural convection cooling of high prandtl number fluid with variable viscosity

Report for the scientific sojourn at the James Cook University, Australia, between June to December 2007. Free convection in enclosed spaces is found widely in natural and industrial systems. It is a topic of primary interest because in many systems it provides the largest resistance to the heat transfer in comparison with other heat transfer modes. In such systems the convection is driven by a density gradient within the fluid, which, usually, is produced by a temperature difference between the fluid and surrounding walls. In the oil industry, the oil, which has High Prandtl, usually is stored and transported in large tanks at temperatures high enough to keep its viscosity and, thus the pumping requirements, to a reasonable level. A temperature difference between the fluid and the walls of the container may give rise to the unsteady buoyancy force and hence the unsteady natural convection. In the initial period of cooling the natural convection regime dominates over the conduction contribution. As the oil cools down it typically becomes more viscous and this increase of viscosity inhibits the convection. At this point the oil viscosity becomes very large and unloading of the tank becomes very difficult. For this reason it is of primary interest to be able to predict the cooling rate of the oil. The general objective of this work is to develop and validate a simulation tool able to predict the cooling rates of high Prandtl fluid considering the variable viscosity effects.

Existence and uniqueness of viscosity solution for Hamilton-Jacobi equation with discontinuous coefficients dependent on time

The main result is a proof of the existence of a unique viscosity solution for Hamilton-Jacobi equation, where the hamiltonian is discontinuous with respect to variable, usually interpreted as the spatial one. Obtained generalized solution is continuous, but not necessarily differentiable.

Lateral instability in normal viscous fingers

We study a low-amplitude, long-wavelength lateral instability of the Saffman-Taylor finger by means of a phase-field model. We observe such an instability in two situations in which small dynamic perturbations are overimposed to a constant pressure drop. We first study the case in which the perturbation consists of a single oscillatory mode and then a case in which the perturbation consists of temporal noise. In both cases the instability undergoes a process of selection.

Relevance of dynamic wetting in viscous fingering patterns

We demonstrate that wetting effects at moving contact lines have a strong impact in viscous fingering patterns. Experiments in a rotating Hele-Shaw (HS) cell, dry or prewetted, show consistent morphological differences. When the wetting fluid invades a dry region, contact angle dynamics yield a kinetic contribution to the interface pressure drop that scales with capillary number as Ca2¿3 but is significantly larger than the Park-Homsy kinetic correction. Numerical results are in very good agreement with experiments and show that standard HS equations work best for prewetted cells.

Phase-field model of Hele-Shaw flows in the high-viscosity contrast regime

A one-sided phase-field model is proposed to study the dynamics of unstable interfaces of Hele-Shaw flows in the high viscosity contrast regime. The corresponding macroscopic equations are obtained by means of an asymptotic expansion from the phase-field model. Numerical integrations of the phase-field model in a rectangular Hele-Shaw cell reproduce finger competition with the final evolution to a steady-state finger.

Radial displacement of a fluid annulus in a rotating Hele¿Shaw cell

The radial displacement of a fluid annulus in a rotating circular Hele-Shaw cell has been investigated experimentally. It has been found that the flow depends sensitively on the wetting conditions at the outer interface. Displacements in a prewet cell are well described by Darcy's law in a wide range of experimental parameters, with little influence of capillary effects. In a dry cell, however, a more careful analysis of the interface motion is required; the interplay between a gradual loss of fluid at the inner interface, and the dependence of capillary forces at the outer interface on interfacial velocity and dynamic contact angle, result in a constant velocity for the interfaces. The experimental results in this case correlate in the form of an empirical scaling relation between the capillary number Ca and a dimensionless group, related to the ratio of centrifugal to capillary forces, which spans about three orders of magnitude in both quantities. Finally, the relative thickness of the coating film left by the inner interface, alpha i, is obtained as a function of Ca.

Low viscosity contrast fingering in a rotating Hele-Shaw cell

We study the fingering instability of a circular interface between two immiscible liquids in a radial Hele-Shaw cell. The cell rotates around its vertical symmetry axis, and the instability is driven by the density difference between the two fluids. This kind of driving allows studying the interfacial dynamics in the particularly interesting case of an interface separating two liquids of comparable viscosity. An accurate experimental study of the number of fingers emerging from the instability reveals a slight but systematic dependence of the linear dispersion relation on the gap spacing. We show that this result is related to a modification of the interface boundary condition which incorporates stresses originated from normal velocity gradients. The early nonlinear regime shows nearly no competition between the outgrowing fingers, characteristic of low viscosity contrast flows. We perform experiments in a wide range of experimental parameters, under conditions of mass conservation (no injection), and characterize the resulting patterns by data collapses of two characteristic lengths: the radius of gyration of the pattern and the interface stretching. Deep in the nonlinear regime, the fingers which grow radially outwards stretch and become gradually thinner, to a point that the fingers pinch and emit drops. We show that the amount of liquid emitted in the first generation of drops is a constant independent of the experimental parameters. Further on there is a sharp reduction of the amount of liquid centrifugated, punctuated by periods of no observable centrifugation.

Interfacial instabilities of a fluid annulus in a rotating Hele-Shaw cell

We have studied the interfacial instabilities experienced by a liquid annulus as it moves radially in a circular Hele-Shaw cell rotating with angular velocity Omega. The instability of the leading interface (oil displacing air) is driven by the density difference in the presence of centrifugal forcing, while the instability of the trailing interface (air displacing oil) is driven by the large viscosity contrast. A linear stability analysis shows that the stability of the two interfaces is coupled through the pressure field already at a linear level. We have performed experiments in a dry cell and in a cell coated with a thin fluid layer on each plate, and found that the stability depends substantially on the wetting conditions at the leading interface. Our experimental results of the number of fingers resulting from the instability compare well with the predictions obtained through a numerical integration of the coupled equations derived from a linear stability analysis. Deep in the nonlinear regime we observe the emission of liquid droplets through the formation of thin filaments at the tip of outgrowing fingers.

Three-dimensional aspects of fluid flows in channels. II. Effects of meniscus and thin film regimes on viscous fingers

We perform a three-dimensional study of steady state viscous fingers that develop in linear channels. By means of a three-dimensional lattice-Boltzmann scheme that mimics the full macroscopic equations of motion of the fluid momentum and order parameter, we study the effect of the thickness of the channel in two cases. First, for total displacement of the fluids in the channel thickness direction, we find that the steady state finger is effectively two-dimensional and that previous two-dimensional results can be recovered by taking into account the effect of a curved meniscus across the channel thickness as a contribution to surface stresses. Second, when a thin film develops in the channel thickness direction, the finger narrows with increasing channel aspect ratio in agreement with experimental results. The effect of the thin film renders the problem three-dimensional and results deviate from the two-dimensional prediction.

Three-dimensional aspects of fluid flows in channels. I. Meniscus and thin film regimes

We study the forced displacement of a fluid-fluid interface in a three-dimensional channel formed by two parallel solid plates. Using a lattice-Boltzmann method, we study situations in which a slip velocity arises from diffusion effects near the contact line. The difference between the slip and channel velocities determines whether the interface advances as a meniscus or a thin film of fluid is left adhered to the plates. We find that this effect is controlled by the capillary and Péclet numbers. We estimate the crossover from a meniscus to a thin film and find good agreement with numerical results. The penetration regime is examined in the steady state. We find that the occupation fraction of the advancing finger relative to the channel thickness is controlled by the capillary number and the viscosity contrast between the fluids. For high viscosity contrast, lattice-Boltzmann results agree with previous results. For zero viscosity contrast, we observe remarkably narrow fingers. The shape of the finger is found to be universal.

Computational calculations of magnetic relaxation and viscosity in small magnetic grains

In this article we present a phenomenological model which simulates very well the mag¿ netic relaxation behavior experimentally observed in small magnetic grains and single domain particles. In this model, the occurrence of quantum tunneling of magnetization below a certain temperature is taken into account. Experimental results for different materials are presented to illustrate the most important behavior deduced from our model

Three-dimensional aspects of fluid flows in channels. II. Effects of meniscus and thin film regimes on viscous fingers

We perform a three-dimensional study of steady state viscous fingers that develop in linear channels. By means of a three-dimensional lattice-Boltzmann scheme that mimics the full macroscopic equations of motion of the fluid momentum and order parameter, we study the effect of the thickness of the channel in two cases. First, for total displacement of the fluids in the channel thickness direction, we find that the steady state finger is effectively two-dimensional and that previous two-dimensional results can be recovered by taking into account the effect of a curved meniscus across the channel thickness as a contribution to surface stresses. Second, when a thin film develops in the channel thickness direction, the finger narrows with increasing channel aspect ratio in agreement with experimental results. The effect of the thin film renders the problem three-dimensional and results deviate from the two-dimensional prediction.