Desenvolvimento de um simulador numérico empregando o método Smoothed Particle Hydrodynamics para a resolução de escoamentos incompressíveis. Implementação computacional em paralelo (CUDA)

Neste trabalho, foi desenvolvido um simulador numérico baseado no método livre de malhas Smoothed Particle Hydrodynamics (SPH) para a resolução de escoamentos de fluidos newtonianos incompressíveis. Diferentemente da maioria das versões existentes deste método, o código numérico faz uso de uma técnica iterativa na determinação do campo de pressões. Este procedimento emprega a forma diferencial de uma equação de estado para um fluido compressível e a equação da continuidade a fim de que a correção da pressão seja determinada. Uma versão paralelizada do simulador numérico foi implementada usando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA) da NVIDIA Corporation. Foram simulados três problemas, o problema unidimensional do escoamento de Couette e os problemas bidimensionais do escoamento no interior de uma Cavidade (Shear Driven Cavity Problem) e da Quebra de Barragem (Dambreak).

Simulação de escoamentos incompressíveis empregando o método Smoothed Particle Hydrodynamics utilizando algoritmos iterativos na determinação do campo de pressões

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas.

Aerodynamic effects in ink-jet printing on a moving web

The airflow between the fast-moving substrate and stationary print heads in a web print press may cause print quality issues in high-speed, roll-to-roll printing applications. We have studied the interactions between ink drops and the airflow in the gap between the printhead and substrate, by using an experimental flow channel and high-speed imaging. The results show: 1) the gap airflow is well approximated by a standard Couette flow profile; 2) the effect of gap airflow on the flight paths of main drops and satellites is negligible; and 3) the interaction between the gap airflow and the wakes from the printed ink drops should be investigated as the primary source of aerodynamically- related print quality issues. ©2012 Society for Imaging Science and Technology.

Influence de l’agrégation érythrocytaire sur la migration axiale de microparticules simulant des plaquettes sanguines

Lors du phénomène d’hémostase primaire ou de thrombose vasculaire, les plaquettes sanguines doivent adhérer aux parois afin de remplir leur fonction réparatrice ou pathologique. Pour ce faire, certains facteurs rhéologiques et hémodynamiques tels que l’hématocrite, le taux de cisaillement local et les contraintes de cisaillement pariétal, entrent en jeu afin d’exclure les plaquettes sanguines de l’écoulement principal et de les transporter vers le site endommagé ou enflammé. Cette exclusion pourrait aussi être influencée par l’agrégation de globules rouges qui est un phénomène naturel présent dans tout le système cardiovasculaire selon les conditions d’écoulement. La dérive de ces agrégats de globules rouges vers le centre des vaisseaux provoque la formation de réseaux d’agrégats dont la taille et la complexité varient en fonction de l’hématocrite et des conditions de cisaillement présentes. Il en résulte un écoulement bi-phasique avec un écoulement central composé d’agrégats de globules rouges avoisinés par une région moins dense en particules où l’on peut trouver des globules rouges singuliers, des petits rouleaux de globules rouges et une importante concentration en plaquettes et globules blancs. De ce fait, il est raisonnable de penser que plus la taille des agrégats qui occupent le centre du vaisseau augmente, plus il y aura de plaquettes expulsées vers les parois vasculaires. L'objectif du projet est de quantifier, in vitro, la migration des plaquettes sanguines en fonction du niveau d’agrégation érythrocytaire présent, en faisant varier l’hématocrite, le taux de cisaillement et en promouvant l’agrégation par l’ajout d’agents tels que le dextran à poids moléculaire élevé. Cependant, le comportement non Newtonien du sang dans un écoulement tubulaire peut être vu comme un facteur confondant à cause de son impact sur l’organisation spatiale des agrégats de globules rouges. De ce fait, les études ont été réalisées dans un appareil permettant de moduler, de façon homogène, la taille et la structure de ces agrégats et de quantifier ainsi leur effet sur la migration axiale des plaquettes. Du sang de porc anti coagulé a été ajusté à différents taux d’hématocrite et insérer dans un appareil à écoulement de Couette, à température ambiante. Les plaquettes sanguines, difficilement isolables in vitro sans en activer certains ligands membranaires, ont été remplacées par des fantômes en polystyrène ayant un revêtement de biotine. La quantification de la migration de ces fantômes de plaquettes a été réalisée grâce à l’utilisation de membranes biologiques fixées sur les parois internes de l’entrefer du rhéomètre de Couette. Ces membranes ont un revêtement de streptavidine assurant une très forte affinité d’adhésion avec les microparticules biotynilées. À 40% d’hématocrite, à un cisaillement de 2 s-1, 566 ± 53 microparticules ont été comptées pour un protocole préétabli avec du sang non agrégeant, comparativement à 1077 ± 229 pour du sang normal et 1568 ± 131 pour du sang hyper agrégeant. Les résultats obtenus suggèrent une nette participation de l’agrégation érythrocytaire sur le transport des fantômes de plaquettes puisque l’adhésion de ces derniers à la paroi du rhéomètre de Couette augmente de façon quasi exponentielle selon le niveau d’agrégation présent.

Spectroscopic signatures of an Fmoc–tetrapeptide, Fmoc and fluorene

The hydrophobic Fmoc [N-(fluorenyl)-9-methoxycarbonyl] unit is often conjugated to peptides to confer amphiphilicity on the molecules, and to introduce the potential for aromatic stacking interactions. In this paper we report spectroscopic data for fluorene, Fmoc and Fmoc conjugated to a small peptide (GRDS in this case) to enable interpretation of spectroscopic data collected on complex structures that can be created with Fmoc–peptide molecules. We report absorbance data (ultra violet and infrared), circular dichroism, linear dichroism (both film and Couette flow orientation) and Raman spectra. Orientations of the Fmoc chromophores in the linear dichroism studies are deduced.

Contributions to the theory of magnetorotational instability and waves in a rotating plasma

The one-fluid magnetohydrodynamic (MHD) theory of magnetorotational instability (MRI) in an ideal plasma is presented. The theory predicts the possibility of MRI for arbitrary 0, where 0 is the ratio of the plasma pressure to the magnetic field pressure. The kinetic theory of MRI in a collisionless plasma is developed. It is demonstrated that as in the ideal MHD, MRI can occur in such a plasma for arbitrary P. The mechanism of MRI is discussed; it is shown that the instability appears because of a perturbed parallel electric field. The electrodynamic description of MRI is formulated under the assumption that the dispersion relation is expressed in terms of the permittivity tensor; general properties of this tensor are analyzed. It is shown to be separated into the nonrotational and rotational parts. With this in mind, the first step for incorporation of MRI into the general theory of plasma instabilities is taken. The rotation effects on Alfven waves are considered.

Desenvolvimento de software utilizando a técnica sph (smoothed particle hydrodynamics) na geração de ondas de submersão

Pós-graduação em Engenharia Elétrica - FEIS

Medidores inteligentes: o primeiro passo em direção as redes inteligentes

Smart grids are the focus of major study today because of the necessity of modernization in electrical systems and reduction of greenhouse gas emissions that increases global warming. Reaching the best deployment method, you must first of all know the current electrical system and how to use them for the benefit of this new technology. Preparing the action plan we should be aware of the main points of smart grids in each step of the electricity system - generation, transmission and distribution. Analyzed these topics, this work will focus on the first step in the implementation of the smart grids: the smart meters, tool which is already being implemented in Brazil. The main characteristics and applications of these devices, as well as their communication structure with the core distributors will be showed during the paper. Finally, we present a case study which will be discussed and analyzed based in the results obtained with the implementation of smart meters in the city of Vancouver, Canada, where we have a considerable savings already in the first year, with fully paying the initial investment and still have a profit

Cross-correlation-aided transport in stochastically driven accretion flows

The origin of linear instability resulting in rotating sheared accretion flows has remained a controversial subject for a long time. While some explanations of such non-normal transient growth of disturbances in the Rayleigh stable limit were available for magnetized accretion flows, similar instabilities in the absence of magnetic perturbations remained unexplained. This dichotomy was resolved in two recent publications by Chattopadhyay and co-workers [Mukhopadhyay and Chattopadhyay, J. Phys. A 46, 035501 (2013)1751-811310.1088/1751-8113/46/3/035501; Nath, Phys. Rev. E 88, 013010 (2013)PLEEE81539-375510.1103/PhysRevE.88.013010] where it was shown that such instabilities, especially for nonmagnetized accretion flows, were introduced through interaction of the inherent stochastic noise in the system (even a "cold" accretion flow at 3000 K is too "hot" in the statistical parlance and is capable of inducing strong thermal modes) with the underlying Taylor-Couette flow profiles. Both studies, however, excluded the additional energy influx (or efflux) that could result from nonzero cross correlation of a noise perturbing the velocity flow, say, with the noise that is driving the vorticity flow (or equivalently the magnetic field and magnetic vorticity flow dynamics). Through the introduction of such a time symmetry violating effect, in this article we show that nonzero noise cross correlations essentially renormalize the strength of temporal correlations. Apart from an overall boost in the energy rate (both for spatial and temporal correlations, and hence in the ensemble averaged energy spectra), this results in mutual competition in growth rates of affected variables often resulting in suppression of oscillating Alfven waves at small times while leading to faster saturations at relatively longer time scales. The effects are seen to be more pronounced with magnetic field fluxes where the noise cross correlation magnifies the strength of the field concerned. Another remarkable feature noted specifically for the autocorrelation functions is the removal of energy degeneracy in the temporal profiles of fast growing non-normal modes leading to faster saturation with minimum oscillations. These results, including those presented in the previous two publications, now convincingly explain subcritical transition to turbulence in the linear limit for all possible situations that could now serve as the benchmark for nonlinear stability studies in Keplerian accretion disks.

Couette Gas Flow Between Cylinders with Different Temperature

Петър Господинов, Добри Данков, Владимир Русинов, Стефан Стефанов - Изследвано е стационарно течение на Кует на разреден газ в случая на въртене на вътрешния цилиндър и неподвижен външен цилиндър чрез използване на DSMC метод и числено решение на уравненията на Навие–Стокс за относително малка (дозвукова) скорост на въртене. Изследвани са различни случаи при промяна на температурата на въртящият се цилиндър и числото на Кнудсен. Целта на изследването е да се установи влиянието на малки скорости на въртене върху макрохарактеристиките – плътността, скоростта и температурата на газа. Установено е добро съвпадение на резултатите получени по двата метода за Kn = 0.02. Получените резултати са важни при решаването на неравнинни, задачи от микрофлуидиката с отчитане на ефектите на кривината. Ключови думи: механика на флуидите, кинетична теория, разреден газ, DSMC.

A dilation-driven vortex flow in sheared granular materials explains a rheometric anomaly

Granular flows occur widely in nature and industry, yet a continuum description that captures their important features is yet not at hand. Recent experiments on granular materials sheared in a cylindrical Couette device revealed a puzzling anomaly, wherein all components of the stress rise nearly exponentially with depth. Here we show, using particle dynamics simulations and imaging experiments, that the stress anomaly arises from a remarkable vortex flow. For the entire range of fill heights explored, we observe a single toroidal vortex that spans the entire Couette cell and whose sense is opposite to the uppermost Taylor vortex in a fluid. We show that the vortex is driven by a combination of shear-induced dilation, a phenomenon that has no analogue in fluids, and gravity flow. Dilatancy is an important feature of granular mechanics, but not adequately incorporated in existing models.

I. Studies on bacteriophage DNA's and minicircular DNA's in M. lysodeikticus and E. coli 15. II. Flow dichroism of DNA solutions

Part I

Chapter 1.....A physicochemical study of the DNA molecules from the three bacteriophages, N1, N5, and N6, which infect the bacterium, M. lysodeikticus, has been made. The molecular weights, as measured by both electron microscopy and sedimentation velocity, are 23 x 10⁶ for N5 DNA and 31 x 10⁶ for N1 and N6 DNA's. All three DNA's are capable of thermally reversible cyclization. N1 and N6 DNA's have identical or very similar base sequences as judged by membrane filter hybridization and by electron microscope heteroduplex studies. They have identical or similar cohesive ends. These results are in accord with the close biological relation between N1 and N6 phages. N5 DNA is not closely related to N1 or N6 DNA. The denaturation T_m of all three DNA's is the same and corresponds to a (GC) content of 70%. However, the buoyant densities in CsCl of Nl and N6 DNA's are lower than expected, corresponding to predicted GC contents of 64 and 67%. The buoyant densities in Cs₂SO₄ are also somewhat anomalous. The buoyant density anomalies are probably due to the presence of odd bases. However, direct base composition analysis of N1 DNA by anion exchange chromatography confirms a GC content of 70%, and, in the elution system used, no peaks due to odd bases are present.

Chapter 2.....A covalently closed circular DNA form has been observed as an intracellular form during both productive and abortive infection processes in M. lysodeikticus. This species has been isolated by the method of CsC1-ethidium bromide centrifugation and examined with an electron microscope.

Chapter 3.....A minute circular DNA has been discovered as a homogeneous population in M. lysodeikticus. Its length and molecular weight as determined by electron microscopy are 0.445 μ and 0.88 x 10⁶ daltons respectively. There is about one minicircle per bacterium.

Chapter 4.....Several strains of E. coli 15 harbor a prophage. Viral growth can be induced by exposing the host to mitomycin C or to uv irradiation. The coliphage 15 particles from E. coli 15 and E, coli 15 T^- appear as normal phage with head and tail structure; the particles from E. coli 15 TAU are tailless. The complete particles exert a colicinogenic activity on E.coli 15 and 15 T^-, the tailless particles do not. No host for a productive viral infection has been found and the phage may be defective. The properties of the DNA of the virus have been studied, mainly by electron microscopy. After induction but before lysis, a closed circular DNA with a contour length of about 11.9 μ is found in the bacterium; the mature phage DNA is a linear duplex and 7.5% longer than the intracellular circular form. This suggests the hypothesis that the mature phage DNA is terminally repetitious and circularly permuted. The hypothesis was confirmed by observing that denaturation and renaturation of the mature phage DNA produce circular duplexes with two single-stranded branches corresponding to the terminal repetition. The contour length of the mature phage DNA was measured relative to φX RFII DNA and λ DNA; the calculated molecular weight is 27 x 10⁶. The length of the single-stranded terminal repetition was compared to the length of φX 174 DNA under conditions where single-stranded DNA is seen in an extended form in electron micrographs. The length of the terminal repetition is found to be 7.4% of the length of the nonrepetitious part of the coliphage 15 DNA. The number of base pairs in the terminal repetition is variable in different molecules, with a fractional standard deviation of 0.18 of the average number in the terminal repetition. A new phenomenon termed "branch migration" has been discovered in renatured circular molecules; it results in forked branches, with two emerging single strands, at the position of the terminal repetition. The distribution of branch separations between the two terminal repetitions in the population of renatured circular molecules was studied. The observed distribution suggests that there is an excluded volume effect in the renaturation of a population of circularly permuted molecules such that strands with close beginning points preferentially renature with each other. This selective renaturation and the phenomenon of branch migration both affect the distribution of branch separations; the observed distribution does not contradict the hypothesis of a random distribution of beginning points around the chromosome.

Chapter 5....Some physicochemical studies on the minicircular DNA species in E. coli 15 (0.670 μ, 1.47 x 10⁶ daltons) have been made. Electron microscopic observations showed multimeric forms of the minicircle which amount to 5% of total DNA species and also showed presumably replicating forms of the minicircle. A renaturation kinetic study showed that the minicircle is a unique DNA species in its size and base sequence. A study on the minicircle replication has been made under condition in which host DNA synthesis is synchronized. Despite experimental uncertainties involved, it seems that the minicircle replication is random and the number of the minicircles increases continuously throughout a generation of the host, regardless of host DNA synchronization.

Part II

The flow dichroism of dilute DNA solutions (A₂₆₀≈0.1) has been studied in a Couette-type apparatus with the outer cylinder rotating and with the light path parallel to the cylinder axis. Shear gradients in the range of 5-160 sec.^-1 were studied. The DNA samples were whole, "half," and "quarter" molecules of T4 bacteriophage DNA, and linear and circular λb₂b₅c DNA. For the linear molecules, the fractional flow dichroism is a linear function of molecular weight. The dichroism for linear A DNA is about 1.8 that of the circular molecule. For a given DNA, the dichroism is an approximately linear function of shear gradient, but with a slight upward curvature at low values of G, and some trend toward saturation at larger values of G. The fractional dichroism increases as the supporting electrolyte concentration decreases.

Simulations of Two-Dimensional Foam under Couette Shear

null Sponsorship: Financial support is acknowledged from the University of Wales Aberystwyth Senate Fund, the Ulysses France-Ireland Exchange Scheme and EPSRC (EP/D014956/1, EP/D048397/1, EP/D071127/1).

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).