Computer simulations of soft particles in flow

Understanding the dynamics of blood cells is a crucial element to discover biological mechanisms, to develop new efficient drugs, design sophisticated microfluidic devices, for diagnostics. In this work, we focus on the dynamics of red blood cells in microvascular flow. Microvascular blood flow resistance has a strong impact on cardiovascular function and tissue perfusion. The flow resistance in microcirculation is governed by flow behavior of blood through a complex network of vessels, where the distribution of red blood cells across vessel cross-sections may be significantly distorted at vessel bifurcations and junctions. We investigate the development of blood flow and its resistance starting from a dispersed configuration of red blood cells in simulations for different hematocrits, flow rates, vessel diameters, and aggregation interactions between red blood cells. Initially dispersed red blood cells migrate toward the vessel center leading to the formation of a cell-free layer near the wall and to a decrease of the flow resistance. The development of cell-free layer appears to be nearly universal when scaled with a characteristic shear rate of the flow, which allows an estimation of the length of a vessel required for full flow development, $l_c \approx 25D$, with vessel diameter $D$. Thus, the potential effect of red blood cell dispersion at vessel bifurcations and junctions on the flow resistance may be significant in vessels which are shorter or comparable to the length $l_c$. The presence of aggregation interactions between red blood cells lead in general to a reduction of blood flow resistance. The development of the cell-free layer thickness looks similar for both cases with and without aggregation interactions. Although, attractive interactions result in a larger cell-free layer plateau values. However, because the aggregation forces are short-ranged at high enough shear rates ($\bar{\dot{\gamma}} \gtrsim 50~\text{s}^{-1}$) aggregation of red blood cells does not bring a significant change to the blood flow properties. Also, we develop a simple theoretical model which is able to describe the converged cell-free-layer thickness with respect to flow rate assuming steady-state flow. The model is based on the balance between a lift force on red blood cells due to cell-wall hydrodynamic interactions and shear-induced effective pressure due to cell-cell interactions in flow. We expect that these results can also be used to better understand the flow behavior of other suspensions of deformable particles such as vesicles, capsules, and cells. Finally, we investigate segregation phenomena in blood as a two-component suspension under Poiseuille flow, consisting of red blood cells and target cells. The spatial distribution of particles in blood flow is very important. For example, in case of nanoparticle drug delivery, the particles need to come closer to microvessel walls, in order to adhere and bring the drug to a target position within the microvasculature. Here we consider that segregation can be described as a competition between shear-induced diffusion and the lift force that pushes every soft particle in a flow away from the wall. In order to investigate the segregation, on one hand, we have 2D DPD simulations of red blood cells and target cell of different sizes, on the other hand the Fokker-Planck equation for steady state. For the equation we measure force profile, particle distribution and diffusion constant across the channel. We compare simulation results with those from the Fokker-Planck equation and find a very good correspondence between the two approaches. Moreover, we investigate the diffusion behavior of target particles for different hematocrit values and shear rates. Our simulation results indicate that diffusion constant increases with increasing hematocrit and depends linearly on shear rate. The third part of the study describes development of a simulation model of complex vascular geometries. The development of the model is important to reproduce vascular systems of small pieces of tissues which might be gotten from MRI or microscope images. The simulation model of the complex vascular systems might be divided into three parts: modeling the geometry, developing in- and outflow boundary conditions, and simulation domain decomposition for an efficient computation. We have found that for the in- and outflow boundary conditions it is better to use the SDPD fluid than DPD one because of the density fluctuations along the channel of the latter. During the flow in a straight channel, it is difficult to control the density of the DPD fluid. However, the SDPD fluid has not that shortcoming even in more complex channels with many branches and in- and outflows because the force acting on particles is calculated also depending on the local density of the fluid.

Sinusite maxilar odontogénica

O seio maxilar é o seio paranasal mais susceptível a invasões bacterianas, tanto pelo óstio nasal, como pela cavidade oral. As sinusites maxilares têm como causas mais frequentes, as infecções víricas, rinites alérgicas ou não alérgicas, variações anatómicas, diabetes mellitus, fumar, nadar, mergulhar, escalar a altas atitudes, e as infecções e tratamentos dentários. A pesquisa bibliográfica, foi realizada sem quaisquer limitações temporais, com restrição linguística a Português, Espanhol e Inglês, sendo excluídos os artigos de outros idiomas; em vários livros e revistas, assim como artigos científicos obtidos, entre Maio e Julho de 2015, nos motores de busca Pubmed, ScienceDirect, Scielo, Elsevier e B-on. A sinusite maxilar odontogénica é uma doença infecto-inflamatória, habitualmente associada à ruptura da membrana de Schneider e a processos infecciosos dentários crónicos. Causa hiperplasia e hipertrofia da mucosa, o que origina sinais e sintomas próprios, assim como mudanças radiográficas perceptíveis. Existem diferentes etiologias de causa odontogénica: cárie, doença periodontal, quistos odontogénicos e iatrogenia – tratamento endodôntico não cirúrgico, cirurgia endodôntica, comunicações oro-antrais, implantes dentários, elevação do seio maxilar, cirurgia pré-protética e cirurgia ortognática – sendo que a iatrogenia é a mais comum (cerca de 56%). Esta patologia afecta com mais frequência indivíduos dos 42,7 aos 51, 7 anos, e preferencialmente a região molar, seguida dos pré-molares e em alguns casos, caninos. Os organismos que dominam na fase aguda e crónica, são sensivelmente os mesmos, mas em número diferente, e existe uma conexão entre a flora comensal periapical e a flora patogénica em caso de sinusite maxilar odontogénica. O diagnóstico é essencialmente clínico, no entanto existem diferentes exames complementares para confirmarem ou formarem o diagnóstico. Pela grande acessibilidade ao método radiográfico, torna-se fundamental que o médico dentista saiba diferencial as diversas patologias que afectam o seio maxilar. O tratamento abrange a eliminação da causa dentária e o tratamento farmacológico, da infecção, essencialmente à base de antibióticos, e da dor se esta existir. E o tratamento cirúrgico, que contempla a punção-lavagem sinusal, antrostomia intranasal, técnica de Caldwell-Luc e cirurgia sinusal endoscópica. Concluindo, o médico dentista deve ter um amplo conhecimento sobre esta patologia para que a possa reconhecer, tratar ou preveni-la.

An experimental investigation by towing tank on VIV of a long flexible cylinder for deepwater riser application

In this study, the dynamic response of a vertical flexible cylinder vibrating at low mode numbers with combined x-y motion was investigated in a towing tank. The uniform flow was simulated by towing the flexible cylinder along the tank in still water; therefore, the turbulence intensity of the free flow was negligible in obtaining more reliable results. A lower branch of dominant frequencies with micro vibration amplitude was found in both cross-flow and in-line directions. This justifiable discrepancy was likely caused by an initial lock-in. The maximum attainable amplitude, modal analysis and x-y trajectory in cross-flow and in-line directions are reported here and compared with previous literature, along with some good agreements and different observations that were obtained from the study. Drag and lift coefficients are also evaluated by making use of a generalized integral transform technique approach, yielding an alternative method to study fluid force acting upon a flexible cylinder.

Coupled CFD Shape Optimization for aerodynamic profiles

The present document deals with the optimization of shape of aerodynamic profiles -- The objective is to reduce the drag coefficient on a given profile without penalising the lift coefficient -- A set of control points defining the geometry are passed and parameterized as a B-Spline curve -- These points are modified automatically by means of CFD analysis -- A given shape is defined by an user and a valid volumetric CFD domain is constructed from this planar data and a set of user-defined parameters -- The construction process involves the usage of 2D and 3D meshing algorithms that were coupled into own- code -- The volume of air surrounding the airfoil and mesh quality are also parametrically defined -- Some standard NACA profiles were used by obtaining first its control points in order to test the algorithm -- Navier-Stokes equations were solved for turbulent, steady-state ow of compressible uids using the k-epsilon model and SIMPLE algorithm -- In order to obtain data for the optimization process an utility to extract drag and lift data from the CFD simulation was added -- After a simulation is run drag and lift data are passed to the optimization process -- A gradient-based method using the steepest descent was implemented in order to define the magnitude and direction of the displacement of each control point -- The control points and other parameters defined as the design variables are iteratively modified in order to achieve an optimum -- Preliminary results on conceptual examples show a decrease in drag and a change in geometry that obeys to aerodynamic behavior principles

A aplicação de plasma rico em plaquetas na regeneração dos tecidos periodontais

O periodonto é uma unidade biológica e funcional, constituída pela gengiva, pelo cemento, pelo ligamento periodontal e pelo osso alveolar. O seu processo de cicatrização envolve mecanismos fisiológicos complexos que requerem a ação dos fatores de crescimento, péptidos oriundos da desgranulação das plaquetas. Neste sentido surge o Plasma Rico em Plaquetas como um produto autólogo, obtido a partir da centrifugação do sangue do próprio paciente e que visa melhorar a cicatrização dos tecidos após procedimentos enquadrados na Medicina Dentária. Esta revisão bibliográfica baseou-se numa pesquisa realizada na base de dados MEDLINE, via pubmed. Foram utilizadas as palavras-chave “plasma rich in growth factors”, “platelet-rich plasma”, “oral surgery”, “dental implants”, “sinus lift”, “third molar surgery” e “bone regeneration”. Após leitura de 40 artigos, foram descartados 9 pela sua pouca relevância no contexto. O objetivo é avaliar a efetividade da aplicação de plasma rico em plaquetas na regeneração dos tecidos periodontais em situações clínicas como alvéolos pós-extracionais, cirurgias de implantes, cirurgias de elevação do seio maxilar e na regeneração óssea. A efetividade em tecidos moles parece ser consensual. A efetividade no tecido ósseo é alvo de contradição entre os diversos autores, concluindo-se que há necessidade de mais estudos randomizados e controlados para que se possa responder à questão com evidência científica suficiente.

SISAL - Um sistema supervisório para elevação artificial de petróleo

MEDEIROS, Adelardo A. D. et al. SISAL - Um Sistema Supervisório para Elevação Artificial de Petróleo. In: Rio Oil and Gas Expo Conference, 2006, Rio de Janeiro, RJ. Anais... Rio de Janeiro, 2006.

Effect of Geometric Uncertainties on the Aerodynamic Characteristic of Offshore Wind Turbine Blades

Offshore wind turbines operate in a complex unsteady flow environment which causes unsteady aerodynamic loads. The unsteady flow environment is characterized by a high degree of uncertainty. In addition, geometry variations and material imperfections also cause uncertainties in the design process. Probabilistic design methods consider these uncertainties in order to reach acceptable reliability and safety levels for offshore wind turbines. Variations of the rotor blade geometry influence the aerodynamic loads which also affect the reliability of other wind turbine components. Therefore, the present paper is dealing with geometric uncertainties of the rotor blades. These can arise from manufacturing tolerances and operational wear of the blades. First, the effect of geometry variations of wind turbine airfoils on the lift and drag coefficients are investigated using a Latin hypercube sampling. Then, the resulting effects on the performance and the blade loads of an offshore wind turbine are analyzed. The variations of the airfoil geometry lead to a significant scatter of the lift and drag coefficients which also affects the damage-equivalent flapwise bending moments. In contrast to that, the effects on the power and the annual energy production are almost negligible with regard to the assumptions made.

Desenvolvimento de um método para correção de curvas de desempenho em bombas centrífugas submersas operando com fluidos viscosos

Electrical Submersible Pump (ESP) is used as an artificial lift technique. However, pumping viscous oil is generally associated with low Reynolds number flows. This condition leads to a performance degradation respect to the performance expected from the regular operation with water that most of the centrifugal pumps are originally designed for. These issues are considered in this investigation through a numerical study of the flow in two different multistage, semi-axial type ESPs. This investigation is carried out numerically using a Computational Fluid Dynamics (CFD) package, where the transient RANS equations are solved numerically. The turbulence is modeled using the SST model. Head curves for several operating conditions are compared with manufacturer’s curves and experimental data for a three-stage ESP, showing good agreement for a wide range of fluid viscosities and rotational speeds. Dimensionless numbers (n, n, n e Re) are used to investigate performance degradation of the ESPs. In addition, flow phenomena through the impellers of the ESPs are investigated using flow field from numerical results. Results show that performance degradation is directly related to rotational Reynolds number, Re. In addition, it was verified that performance degradation occurs for constant normalized specific speedn, which shows that performance degradation occurs similarly for different centrifugal pumps. Moreover, experimental data and numerical results agreed with a correlation from literature between head and flow correction factors proposed by Stepanoff (1967). A definition of modified Reynolds number was proposed and relates the head correction factor to viscosity. A correlation between head correction factor and the modified Reynolds number was proposed, which agreed well with numerical and experimental data. Then, a method to predict performance degradation based on the previous correlations was proposed. This method was compared with others from literature. In general, results and conclusions from this work can also be useful to bring more information about the flow of highly viscous fluids in pumps, especially in semi-axial, multistage ESPs.

Aerodynamic behavior of an airfoil with morphing trailing edge for wind turbine applications

The length of wind turbine rotor blades has been increased during the last decades. Higher stresses arise especially at the blade root because of the longer lever arm. One way to reduce unsteady blade-root stresses caused by turbulence, gusts, or wind shear is to actively control the lift in the blade tip region. One promising method involves airfoils with morphing trailing edges to control the lift and consequently the loads acting on the blade. In the present study, the steady and unsteady behavior of an airfoil with a morphing trailing edge is investigated. Two-dimensional Reynolds-Averaged Navier-Stokes (RANS) simulations are performed for a typical thin wind turbine airfoil with a morphing trailing edge. Steady-state simulations are used to design optimal geometry, size, and deflection angles of the morphing trailing edge. The resulting steady aerodynamic coefficients are then analyzed at different angles of attack in order to determine the effectiveness of the morphing trailing edge. In order to investigate the unsteady aerodynamic behavior of the optimal morphing trailing edge, time-resolved RANS-simulations are performed using a deformable grid. In order to analyze the phase shift between the variable trailing edge deflection and the dynamic lift coefficient, the trailing edge is deflected at four different reduced frequencies for each different angle of attack. As expected, a phase shift between the deflection and the lift occurs. While deflecting the trailing edge at angles of attack near stall, additionally an overshoot above and beyond the steady lift coefficient is observed and evaluated.

Sistema supervisor inteligente para procesos de producción de petróleo

Maximizar la producción de pozos de crudo pesado y extra pesado es el principal beneficio que se desea obtener de los sistemas de control que están corrientemente operativos en empresas de petróleo. Dada la naturaleza compleja y cambiante con el tiempo de los métodos existentes de levantamiento artificial para extracción de crudo, se dificulta el cumplimiento de las especificaciones pre establecidas para el procesamiento del crudo por parte de los lazos de control regulatorios. Tomando esto en cuenta, en éste trabajo se propone un sistema de supervisión inteligente que permite detectar cambios en las condiciones de operación del proceso productivo y realizar ajustes automáticos de sus consignas. Además, el sistema supervisor propuesto tiene la capacidad de detectar fallas en los sensores involucrados en los lazos de control, garantizando de esta manera una operación confiable del proceso. La propuesta fue probada en un pozo de petróleo real obteniéndose resultados que superaron las expectativas iniciales.

An Optimal Order Interior Penalty Discontinuous Galerkin Discretization of the Compressible Navier-Stokes Equations

In this article we propose a new symmetric version of the interior penalty discontinuous Galerkin finite element method for the numerical approximation of the compressible Navier-Stokes equations. Here, particular emphasis is devoted to the construction of an optimal numerical method for the evaluation of certain target functionals of practical interest, such as the lift and drag coefficients of a body immersed in a viscous fluid. With this in mind, the key ingredients in the construction of the method include: (i) An adjoint consistent imposition of the boundary conditions; (ii) An adjoint consistent reformulation of the underlying target functional of practical interest; (iii) Design of appropriate interior-penalty stabilization terms. Numerical experiments presented within this article clearly indicate the optimality of the proposed method when the error is measured in terms of both the L_2-norm, as well as for certain target functionals. Computational comparisons with other discontinuous Galerkin schemes proposed in the literature, including the second scheme of Bassi & Rebay, cf. [11], the standard SIPG method outlined in [25], and an NIPG variant of the new scheme will be undertaken.

On MDS convolutional codes over Z_p^r

Maximum distance separable (MDS) convolutional codes are characterized through the property that the free distance meets the generalized Singleton bound. The existence of free MDS convolutional codes over Zpr was recently discovered in Oued and Sole (IEEE Trans Inf Theory 59(11):7305–7313, 2013) via the Hensel lift of a cyclic code. In this paper we further investigate this important class of convolutional codes over Zpr from a new perspective. We introduce the notions of p-standard form and r-optimal parameters to derive a novel upper bound of Singleton type on the free distance. Moreover, we present a constructive method for building general (non necessarily free) MDS convolutional codes over Zpr for any given set of parameters.

Comparison of two methods for measuring the pubococcygeal line from sagittal-plane magnetic resonance imaging

Aims The pubococcygeal line (PCL) is an important reference line for determining measures of pelvic organ support on sagittal-plane magnetic resonance imaging (MRI); however, there is no consensus on where to place the posterior point of the PCL. As coccyx movement produced during pelvic floor muscle (PFM) contractions may affect other measures, optimal placement of the posterior point is important. This study compared two methods for measuring the PCL, with different posterior points, on T2-weighted sagittal MRI to determine the effect of coccygeal movement on measures of pelvic organ support in older women. Methods MRI of the pelvis was performed in the midsagittal plane, at rest and during PFM contractions, on 47 community-dwelling women 60 and over. The first PCL was measured to the tip of the coccyx (PCLtip) and the second to the sacrococcygeal joint (PCLjnt). Four measures of pelvic organ support were made using each PCL as the reference line: urethrovesical junction height, uterovaginal junction height, M-line and levator plate angle. Results During the PFM contraction the PCLtip shortened and lifted (P < 0.001); the PCLjnt did not change (P > 0.05). The changes in the four measures of pelvic organ support were smaller when measured relative to the PCLtip as compared to those to the PCLjnt (P < 0.001). Conclusions Coccyx movement affected the length and position of the PCLtip, which resulted in underestimates of the pelvic-organ lift produced by the PFM contraction. Therefore, we recommend that the PCL be measured to the sacrococcygeal joint and not to the tip of the coccyx

Construção, caracterização e aplicação de eletrodos para descarga por arco de cátodo oco

The system built to characterize electrodes and, consequently, deposited fine films are constituted by a hollow cathode that works to discharges and low pressures (approximately 10-3 to 5 mbar), a source DC (0 to 1200 V), a cylindrical camera of closed borossilicato for flanges of stainless steel with an association of vacuum bombs mechanical and spread. In the upper flange it is connected the system of hollow cathode, which possesses an entrance of gas and two entrances for its refrigeration, the same is electrically isolated of the rest of the equipment and it is polarized negatively. In front of the system of hollow cathode there is a movable sample in stainless steel with possibility of moving in the horizontal and vertical. In the vertical, the sample can vary its distance between 0 and 70 mm and, in the horizontal, can leave completely from the front of the hollow cathode. The sample and also the cathode hollow are equipped with cromel-alumel termopares with simultaneous reading of the temperatures during the time of treatment. In this work copper electrodes, bronze, titanium, iron, stainless steel, powder of titanium, powder of titanium and silício, glass and ceramic were used. The electrodes were investigated relating their geometry change and behavior of the plasma of the cavity of hollow cathode and channel of the gas. As the cavity of hollow cathode, the analyzed aspects were the diameter and depth. With the channel of the gas, we verified the diameter. In the two situations, we investigated parameters as flow of the gas, pressure, current and applied tension in the electrode, temperature, loss of mass of the electrode with relationship at the time of use. The flow of gas investigated in the electrodes it was fastened in a work strip from 15 to 6 sccm, the constant pressure of work was among 2.7 to 8 x 10-2 mbar. The applied current was among a strip of work from 0,8 to 0,4 A, and their respective tensions were in a strip from 400 to 220 V. Fixing the value of the current, it was possible to lift the curve of the behavior of the tension with the time of use. That curves esteem in that time of use of the electrode to its efficiency is maximum. The temperatures of the electrodes were in the dependence of that curves showing a maximum temperature when the tension was maximum, yet the measured temperatures in the samples showed to be sensitive the variation of the temperature in the electrodes. An accompaniment of the loss of mass of the electrode relating to its time of use showed that the electrodes that appeared the spherical cavities lost more mass in comparison with the electrodes in that didn't appear. That phenomenon is only seen for pressures of 10-2 mbar, in these conditions a plasma column is formed inside of the channel of the gas and in certain points it is concentrated in form of spheres. Those spherical cavities develop inside of the channel of the gas spreading during the whole extension of the channel of the gas. The used electrodes were cut after they could not be more used, however among those electrodes, films that were deposited in alternate times and the electrodes that were used to deposit films in same times, those films were deposited in the glass substrata, alumina, stainless steel 420, stainless steel 316, silício and steel M2. As the eletros used to deposit films in alternate time as the ones that they were used to deposit in same times, the behavior of the thickness of the film obeyed the curve of the tension with relationship the time of use of the electrode, that is, when the tension was maximum, the thickness of the film was also maximum and when the tension was minimum, the thickness was minimum and in the case where the value of the tension was constant, the thickness of the film tends to be constant. The fine films that were produced they had applications with nano stick, bio-compatibility, cellular growth, inhibition of bacterias, cut tool, metallic leagues, brasagem, pineapple fiber and ornamental. In those films it was investigated the thickness, the adherence and the uniformity characterized by sweeping electronic microscopy. Another technique developed to assist the production and characterization of the films produced in that work was the caloteste. It uses a sphere and abrasive to mark the sample with a cap impression, with that cap form it is possible to calculate the thickness of the film. Through the time of life of the cathode, it was possible to evaluate the rate of waste of its material for the different work conditions. Values of waste rate up to 3,2 x 10-6 g/s were verified. For a distance of the substratum of 11 mm, the deposited film was limited to a circular area of 22 mm diameter mm for high pressures and a circular area of 75 mm for pressure strip. The obtained films presented thickness around 2,1 µm, showing that the discharge of arch of hollow cathode in argon obeys a curve characteristic of the tension with the time of life of the eletrodo. The deposition rate obtained in this system it is of approximately 0,18 µm/min

Aerodynamic parametric analysis of an unconventional joined-wing aircraft configuration

In this paper, the experimental results of an unconventional joined-wing aircraft configuration are presented. The test model uses two different wings, forward and rear, both joined in tandem and forming diamond shapes both in plant and front views. The wings are joined in such a way that it is possible to change the rear wing dihedral angle values and the rear wing sweep angle values in 25 different positions that modify the relative distance and the relative height between the wings. To measure the system aerodynamic coefficients itis necessary to perform wind tunnel tests. The datapresented corresponds to the lift, drag and induced drag aerodynamic coefficients, as well as the aerodynamic efficiency and the parameter for minimum required power, from the calculated values of the lift and drag time series measured by a 6-axis force and torque sensor. The results show the influence on the aerodynamic coefficients of the rear wing sweep and dihedral angles parameters. As a main result, it can be concluded that, in general terms, the lift and induced drag aerodynamic coefficients values decrease as both the distance and height between the wings increase, on the other hand, the total drag aerodynamic coefficient decreases if the distance between the wings increases, but nevertheless shows a slight tendency to increase if the height of the rear wing increases, whereas the aerodynamic efficiency and the parameter for minimum required power increase if the distance between the wings increases