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.

Simulación de fluidos inmiscibles con el Particle Finite Element Method (PFEM)

Proyecto de investigación realizado a partir de una estancia en el Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC), Argentina, entre febrero y abril del 2007. La simulación numérica de problemas de mezclas mediante el Particle Finite Element Method (PFEM) es el marco de estudio de una futura tesis doctoral. Éste es un método desarrollado conjuntamente por el CIMEC y el Centre Internacional de Mètodos Numèrics en l'Enginyeria (CIMNE-UPC), basado en la resolución de las ecuaciones de Navier-Stokes en formulación Lagrangiana. El mallador ha sido implementado y desarrollado por Dr. Nestor Calvo, investigador del CIMEC. El desarrollo del módulo de cálculo corresponde al trabajo de tesis de la beneficiaria. La correcta interacción entre ambas partes es fundamental para obtener resultados válidos. En esta memoria se explican los principales aspectos del mallador que fueron modificados (criterios de refinamiento geométrico) y los cambios introducidos en el módulo de cálculo (librería PETSc, algoritmo predictor-corrector) durante la estancia en el CIMEC. Por último, se muestran los resultados obtenidos en un problema de dos fluidos inmiscibles con transferencia de calor.

Contact melting of a three-dimensional phase change material on a flat substrate

In this paper a model is developed to describe the three dimensional contact melting process of a cuboid on a heated surface. The mathematical description involves two heat equations (one in the solid and one in the melt), the Navier-Stokes equations for the flow in the melt, a Stefan condition at the phase change interface and a force balance between the weight of the solid and the countering pressure in the melt. In the solid an optimised heat balance integral method is used to approximate the temperature. In the liquid the small aspect ratio allows the Navier-Stokes and heat equations to be simplified considerably so that the liquid pressure may be determined using an igenfunction expansion and finally the problem is reduced to solving three first order ordinary differential equations. Results are presented showing the evolution of the melting process. Further reductions to the system are made to provide simple guidelines concerning the process. Comparison of the solutions with experimental data on the melting of n-octadecane shows excellent agreement.

Discretized integral hydrodynamics

Using an interpolant form for the gradient of a function of position, we write an integral version of the conservation equations for a fluid. In the appropriate limit, these become the usual conservation laws of mass, momentum, and energy. We also discuss the special cases of the Navier-Stokes equations for viscous flow and the Fourier law for thermal conduction in the presence of hydrodynamic fluctuations. By means of a discretization procedure, we show how the integral equations can give rise to the so-called particle dynamics of smoothed particle hydrodynamics and dissipative particle dynamics.

Transient shear banding in time-dependent fluids

We study the dynamics of shear-band formation and evolution using a simple rheological model. The description couples the local structure and viscosity to the applied shear stress. We consider in detail the Couette geometry, where the model is solved iteratively with the Navier-Stokes equation to obtain the time evolution of the local velocity and viscosity fields. It is found that the underlying reason for dynamic effects is the nonhomogeneous shear distribution, which is amplified due to a positive feedback between the flow field and the viscosity response of the shear thinning fluid. This offers a simple explanation for the recent observations of transient shear banding in time-dependent fluids. Extensions to more complicated rheological systems are considered.

Accounting for adsorption and desorption in Lattice Boltzmann simulations

We report a Lattice-Boltzmann scheme that accounts for adsorption and desorption in the calculation of mesoscale dynamical properties of tracers in media of arbitrary complexity. Lattice Boltzmann simulations made it possible to solve numerically the coupled Navier-Stokes equations of fluid dynamics and Nernst-Planck equations of electrokinetics in complex, heterogeneous media. With the moment propagation scheme, it became possible to extract the effective diffusion and dispersion coefficients of tracers, or solutes, of any charge, e.g., in porous media. Nevertheless, the dynamical properties of tracers depend on the tracer-surface affinity, which is not purely electrostatic and also includes a species-specific contribution. In order to capture this important feature, we introduce specific adsorption and desorption processes in a lattice Boltzmann scheme through a modified moment propagation algorithm, in which tracers may adsorb and desorb from surfaces through kinetic reaction rates. The method is validated on exact results for pure diffusion and diffusion-advection in Poiseuille flows in a simple geometry. We finally illustrate the importance of taking such processes into account in the time-dependent diffusion coefficient in a more complex porous medium.

Simulació de vent en temps real per paisatges

La simulació de la realitat és un fenomen que va sorgir fa uns anys per tal de predir esdeveniments sense haver de malbaratar recursos. El problema inicial de la simulació va ser la necessitat de simplificar la realitat a causa de la manca de capacitat dels ordinadors de l’època. Amb aquest projecte volem ajudar, per exemple, a estudis científics sobre la difusió de la contaminació en grans nuclis a causa de l’efecte del vent, càlculs de trajectòries amb forces externes degudes al vent, o incorporar en el món de la multimèdia efectes realistes de vent. El principal objectiu d’aquest projecte és desenvolupar un sistema que permeti realitzar simulacions realistes de vent per un paisatge 2D, i estudiar com el vent és veu afectat per la geometria de l’escena. Un punt important, és que tot ha de ser en temps real. Per aconseguir-ho, utilitzarem tècniques basades en el mètode de Lattice-Boltzmann, el qual consisteix en una xarxa regular que representa el fluid en posicions discretes, i estudiar com flueix. Escollint els paràmetres correctes de la simulació, es pot demostrar que aquest mètode convergeix a les equacions continues de Navier-Stokes, les qual són les més importants per descriure el comportament macroscòpic d’un fluid. Per accelerar tots els càlculs, utilitzarem la capacitat i la potencia de les targes gràfiques, ajustarem l’algorisme per poder-lo utilitzar en paral•lel, tot tenint en compte les restriccions de les GPUs. També haurem de generar un sistema per poder llegir les escenes 2D sobre les que realitzarem la simulació. Finalment, haurem de “pintar” el vent per tal de poder visualitzar el resultat de la simulació

A model of the wind direction signature in the Stokes emission vector from the ocean surface at microwave frequencies

This paper presents a model of the Stokes emission vector from the ocean surface. The ocean surface is described as an ensemble of facets with Cox and Munk's (1954) Gram-Charlier slope distribution. The study discusses the impact of different up-wind and cross-wind rms slopes, skewness, peakedness, foam cover models and atmospheric effects on the azimuthal variation of the Stokes vector, as well as the limitations of the model. Simulation results compare favorably, both in mean value and azimuthal dependence, with SSM/I data at 53° incidence angle and with JPL's WINDRAD measurements at incidence angles from 30° to 65°, and at wind speeds from 2.5 to 11 m/s.

Caraterización y uso de pantallas de cristal líquido para diferentes aplicaciones

Proyecto de investigación realizado a partir de una estancia en el Departamento de Física de la Universidad Autònoma de Barcelona entre octubre de 2007 y febrero del 2008. La investigación realizada se ha centrado en la caracterización y uso en diferentes aplicaciones de pantallas de cristal líquido por reflexión. Las fluctuaciones temporales de la modulación de estos dispositivos origina una despolarización parcial de la luz transmitida. Por ello para su caracterización se ha utilizado el formalismo de Mueller-Stokes. Se ha desarrollado un montaje experimental para la caracterización y se ha propuesto un modelo para predecir tanto su modulación en maplictud come en fase y polarización. También se ha estudiado la variación del comportamiento en función de la longitud de onda utilizada y del ángulo de incidencia. La respuesta de los sistemas ópticos en cuanto a su profundidad de enfoque y su resolución puede ser modificada mediante filtros de transmisión no uniforme. Hemos estudiado la capacidad de las pantallas de cristal líquido para implementar este tipo de filtros. Se han estudiado los efectos de diferentes tipos de filtros tanto teóricamente como experimentalmente. Se han diseñado sistemas ópticos telescópicos utilizando dos pantallas de cristal líquido. En estos sistemas se puede variar el aumento en tiempo real. Debido a la versatilidad que ofrecen las pantallas de cristal líquido se pueden implementar lentes cilíndricas con lo que se pueden variar los aumentos en cada una de las direcciones, formando así un procesador anamórfico programable. Por último se ha propuesto un procesador óptico de polarización. Las pantallas de cristal líquido cambian el estado de polarización de la luz incidente, lo que da lugar a la posibilidad de modular espacialmente el estado de polarización. Se han estudiado teóricamente y experimentalmente la propagación de estos haces.

Could cell membranes produce acoustic streaming? Making the case for synechococcus self-propulsion

Sir James Lighthill proposed in 1992 that acoustic streaming occurs in the inner ear, as part of the cochlear amplifier mechanism. Here we hypothesize that some of the most ancient organisms use acoustic streaming not only for self-propulsion but also to enhance their nutrient uptake. We focus on a motile strain of Synechococcus, a yanobacteria whose mechanism for self-propulsion is not known. Molecular motors could work like piezoelectric transducers acting on the crystalline structure surrounding the outer cell membrane. Our calculations show that a traveling surface acoustic wave (SAW)could account for the observed velocities. These SAW waves will also produce a non-negligible Stokes layer surrounding the cell: motion within this region being essentially chaotic. Therefore, an AS mechanism would be biologically advantageous, enhancing localized diffusion processes and consequently, chemical reactions. We believe that acoustic streaming, produced by nanometer scale membrane vibrations could be widespread in cell biology. Other possible instances are yeast cells and erythrocytes. Flows generated by acoustic streaming may also be produced by silica coated diatoms along their raphe. We note that microelectromechanical (MEMS) acoustic streaming devices were first introduced in the 1990’s. Nature may have preceded this invention by 2.7 Gyr.

Why are slip lengths so large in carbon nanotubes?

The enhanced flow in carbon nanotubes is explained using a mathematical model that includes a depletion layer with reduced viscosity near the wall. In the limit of large tubes the model predicts no noticeable enhancement. For smaller tubes the model predicts enhancement that increases as the radius decreases. An analogy between the reduced viscosity and slip-length models shows that the term slip-length is misleading and that on surfaces which are smooth at the nanoscale it may be thought of as a length-scale associated with the size of the depletion region and viscosity ratio. The model therefore provides a physical interpretation of the classical Navier slip condition and explains why `slip-lengths' may be greater than the tube radius.

Influence of average size and interface passivation on the spectral emission of Si nanocrystals embedded in SiO2

The correlation between the structural (average size and density) and optoelectronic properties [band gap and photoluminescence (PL)] of Si nanocrystals embedded in SiO2 is among the essential factors in understanding their emission mechanism. This correlation has been difficult to establish in the past due to the lack of reliable methods for measuring the size distribution of nanocrystals from electron microscopy, mainly because of the insufficient contrast between Si and SiO2. With this aim, we have recently developed a successful method for imaging Si nanocrystals in SiO2 matrices. This is done by using high-resolution electron microscopy in conjunction with conventional electron microscopy in dark field conditions. Then, by varying the time of annealing in a large time scale we have been able to track the nucleation, pure growth, and ripening stages of the nanocrystal population. The nucleation and pure growth stages are almost completed after a few minutes of annealing time at 1100°C in N2 and afterward the ensemble undergoes an asymptotic ripening process. In contrast, the PL intensity steadily increases and reaches saturation after 3-4 h of annealing at 1100°C. Forming gas postannealing considerably enhances the PL intensity but only for samples annealed previously in less time than that needed for PL saturation. The effects of forming gas are reversible and do not modify the spectral shape of the PL emission. The PL intensity shows at all times an inverse correlation with the amount of Pb paramagnetic centers at the Si-SiO2 nanocrystal-matrix interfaces, which have been measured by electron spin resonance. Consequently, the Pb centers or other centers associated with them are interfacial nonradiative channels for recombination and the emission yield largely depends on the interface passivation. We have correlated as well the average size of the nanocrystals with their optical band gap and PL emission energy. The band gap and emission energy shift to the blue as the nanocrystal size shrinks, in agreement with models based on quantum confinement. As a main result, we have found that the Stokes shift is independent of the average size of nanocrystals and has a constant value of 0.26±0.03 eV, which is almost twice the energy of the Si¿O vibration. This finding suggests that among the possible channels for radiative recombination, the dominant one for Si nanocrystals embedded in SiO2 is a fundamental transition spatially located at the Si¿SiO2 interface with the assistance of a local Si-O vibration.

Brownian motion of multidimensional systems in nonpotential velocity-dependent fields of force

We study the Brownian motion in velocity-dependent fields of force. Our main result is a Smoluchowski equation valid for moderate to high damping constants. We derive that equation by perturbative solution of the Langevin equation and using functional derivative techniques.

Relaxation dynamics in suspensions of ferromagnetic particles

We have studied the relaxation dynamics of a dilute assembly of ferromagnetic particles in suspension. A formalism based on the Smoluchowski equation, describing the evolution of the probability density for the directions of the magnetic moment and of the axis of easy magnetization of the particles, has been developed. We compute the rotational viscosity from a Green-Kubo formula and give an expression for the relaxation time of the particles which comes from the dynamic equations of the correlation functions. Concerning the relaxation time for the particles, our results agree quite well with experiments performed on different samples of ferromagnetic particles for which the magnetic energy, associated with the interaction between the magnetic moments and the external field, or the energy of anisotropy plays a dominant role.

Aggregation of liposomes induced by calcium: A structural and kinetic study

In this work, the calcium-induced aggregation of phosphatidylserine liposomes is probed by means of the analysis of the kinetics of such process as well as the aggregate morphology. This novel characterization of liposome aggregation involves the use of static and dynamic light-scattering techniques to obtain kinetic exponents and fractal dimensions. For salt concentrations larger than 5 mM, a diffusion-limited aggregation regime is observed and the Brownian kernel properly describes the time evolution of the diffusion coefficient. For slow kinetics, a slightly modified multiple contact kernel is required. In any case, a time evolution model based on the numerical resolution of Smoluchowski's equation is proposed in order to establish a theoretical description for the aggregating system. Such a model provides an alternative procedure to determine the dimerization constant, which might supply valuable information about interaction mechanisms between phospholipid vesicles.