Computational fluid dynamics investigation of a centrifugal blood pump

In the development of a ventricular assist device, computational fluid dynamics (CFD) analysis is an efficient tool to obtain the best design before making the final prototype. In this study, different designs of a centrifugal blood pump were developed to investigate flow characteristics and performance. This study assumed the blood flow as being an incompressible homogeneous Newtonian fluid. A constant velocity was applied at the inlet; no slip boundary conditions were applied at device wall; and pressure boundary conditions were applied at the outlet. The CFD code used in this work was based on the finite volume method. In the future, the results of CFD analysis can be compared with flow visualization and hemolysis tests.

The influence of immobilized biomass and particle size on the fluid dynamics of an anaerobic fluidized bed reactor

Fluid dynamic analysis is an important branch of several chemical engineering related areas, such as drying processes and chemical reactors. However, aspects concerning fluid dynamics in wastewater treatment bioreactors still require further investigation, as they highly influence process efficiency. Therefore, it is essential to evaluate the influence of biofilm on the reactor fluid dynamic behavior, through the analysis of a few important parameters, such as minimum fluidization velocity, bed expansion and porosity, and particle terminal velocity. The main objective of the present work was to investigate the fluid dynamics of an anaerobic fluidized bed reactor, having activated carbon particles as support media for biomass immobilization. Reactor performance was tested using synthetic residual water, which was prepared using the solution employed in BOD determination. The results showed that the presence of immobilized biomass increased particle density and altered the main fluid dynamic parameters investigated.

Flow Visualization in Arteriovenous Fistula and Aneurysm Using Computational Fluid Dynamics

The arteriovenous fistula (AVF) is characterized by enhanced blood flow and is the most widely used vascular access for chronic haemodialysis (Sivanesan et al., 1998). A large proportion of the AVF late failures are related to local haemodynamics (Sivanesan et al., 1999a). As in AVF, blood flow dynamics plays an important role in growth, rupture, and surgical treatment of aneurysm. Several techniques have been used to study the flow patterns in simplified models of vascular anastomose and aneurysm. In the present investigation, Computational Fluid Dynamics (CFD) is used to analyze the flow patterns in AVF and aneurysm through the velocity waveform obtained from experimental surgeries in dogs (Galego et al., 2000), as well as intra-operative blood flow recordings of patients with radiocephalic AVF ( Sivanesan et al., 1999b) and physiological pulses (Aires, 1991), respectively. The flow patterns in AVF for dog and patient surgeries data are qualitatively similar. Perturbation, recirculation and separation zones appeared during cardiac cycle, and these were intensified in the diastole phase for the AVF and aneurysm models. The values of wall shear stress presented in this investigation of AVF and aneurysm models oscillated in the range that can both cause damage to endothelial cells and develop atherosclerosis.

DEVELOPMENT OF STATIC MIXERS FOR MISCIBLE FLUIDS IN LAMINAR FLOW WITH THE USE OF COMPUTATIONAL FLUID DYNAMICS (CFD)

Static mixers with improved performance were developed from CFD simulations in a stepwise approach. The relevant geometric features of simple mixer designs and the corresponding mixing mechanisms-laminar shear, elongational flow, and distributive mixing-were identified first. This information was used to formulate guidelines for the development of new geometries. The solid elements of the static mixer should: (a) provide restrictions to the flow; (b) deflect the flow; (c) be sequentially rotated around the flow direction to provide symmetry; (d) extend from the center of the pipe to the vicinity of the walls to avoid short-circuiting; and (e) distribute and remix the flow. Based on these guidelines, two improved mixer designs were developed: the DS A-I mixer has a good mixing efficiency and an acceptable pressure drop; the Fins 35 degrees mixer is more efficient and compact, but requires a larger pressure drop. Their performance indicates that their use is possible on industrial applications.

Thermal Characteristics of Polybia scutellaris Nests (Hymenoptera: Vespidae) Using Computational Fluid Dynamics: A Possible Adaptation to Tropical Climates

Polybia scutellaris constructs huge nests characterized by numerous spinal projections on the surface. We investigated the thermal characteristics of P scutellaris nests in order to determine whether their nest temperature is homeothermically maintained and whether the spines play a role in the thermoregulation of the nests. In order to examine these hypotheses, we measured the nest temperature in a active nest and in an abandoned nest. The temperature in the active nest was almost stable at 27 degrees C, whereas that of the abandoned nest varied with changes in the ambient temperature, suggesting that nest temperature was maintained by the thermogenesis of colony individuals. In order to predict the thermal properties of the spines, a numerical simulation was employed. To construct a 3D-model of a P scutellaris nest, the nest architecture was simplified into an outer envelope and the surface spines, for both of which the initial temperature was set at 27 degrees C. The physical properties of the simulated nest were regarded to be those of wood since the nest of this species is constructed from plant materials. When the model was exposed to cool air (12 degrees C), the temperature was lower in the models with more spines. On the other hand, when the nest was heated (42 degrees C), the temperature increase was smaller in models with more spines. It is suggested that the spines act as a heat radiator, not as an insulator, against the changes in ambient temperature.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov-Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics. (C) 2007 Elsevier B.V. All rights reserved.

Investigation of a method for predicting AUV derivatives

The paper addresses the problem of autonomous underwater vehicle (AUV) modelling and parameter estimation as a means to predict the dynamic performance of underwater vehicles and thus provide solid guidelines during their design phase. The use of analytical and semi-empirical (ASE) methods to estimate the hydrodynamic derivatives of a popular class of AUVs is discussed. A comparison is done with the results obtained by using computational fluid dynamics to evaluate the bare hull lift force distribution around a fully submerged body. An application is made to the estimation of the hydrodynamic derivatives of the MAYA AUV, an autonomous underwater vehicle developed under a joint Indian-Portuguese project. The estimates obtained were used to predict the turning diameter of the vehicle during sea trials. (C) 2008 Elsevier Ltd. All rights reserved.

Inverse aerodynamic design applications using the MGM hybrid formulation

The well-known modified Garabedian-Mcfadden (MGM) method is an attractive alternative for aerodynamic inverse design, for its simplicity and effectiveness (P. Garabedian and G. Mcfadden, Design of supercritical swept wings, AIAA J. 20(3) (1982), 289-291; J.B. Malone, J. Vadyak, and L.N. Sankar, Inverse aerodynamic design method for aircraft components, J. Aircraft 24(2) (1987), 8-9; Santos, A hybrid optimization method for aerodynamic design of lifting surfaces, PhD Thesis, Georgia Institute of Technology, 1993). Owing to these characteristics, the method has been the subject of several authors over the years (G.S. Dulikravich and D.P. Baker, Aerodynamic shape inverse design using a Fourier series method, in AIAA paper 99-0185, AIAA Aerospace Sciences Meeting, Reno, NV, January 1999; D.H. Silva and L.N. Sankar, An inverse method for the design of transonic wings, in 1992 Aerospace Design Conference, No. 92-1025 in proceedings, AIAA, Irvine, CA, February 1992, 1-11; W. Bartelheimer, An Improved Integral Equation Method for the Design of Transonic Airfoils and Wings, AIAA Inc., 1995). More recently, a hybrid formulation and a multi-point algorithm were developed on the basis of the original MGM. This article discusses applications of those latest developments for airfoil and wing design. The test cases focus on wing-body aerodynamic interference and shock wave removal applications. The DLR-F6 geometry is picked as the baseline for the analysis.

Stability and in vitro penetration study of rutin incorporated in a cosmetic emulsion through an alternative model biomembrane

Rutin is employed as antioxidant and to prevent the capillary fragility and, when incorporated in cosmetic emulsions, it must target the action site. In vitro cutaneous penetration studies through human skin is the ideal situation, however, there are difficulties to obtain and to maintain this tissue viability. Among the membrane models, shed snake skin presents itself as pure stratum corneum, providing barrier function similar to human and it is obtained without the animal sacrifice. The objectives of this research were the development and stability evaluation of a cosmetic emulsion containing rutin and propylene glycol (penetration enhancer) and the evaluation or rutin in vitro cutaneous penetration and retention from the emulsion, employing an alternative model biomembrane. Emulsion was developed with rutin and propylene glycol, both at 5.0% w/w. Active substance presented on the formulation was quantified by a validated spectrophotometric method at 361.0 nm. Rutin Rutin cutaneous penetration and retention was performed in vertical diffusion cells with shed snake skin of Crotalus durissus, as alternative model biomembrane, and distilled water and ethanol 99.5% (1:1), as receptor fluid. The experiment was conducted for six hours, at 37.0 +/- 0.5 degrees C with constant stirring of 300 rpm. Spectrophotometry at 410.0 nm, previously validated, determined the active substance after cutaneous penetration/ retention. Emulsion did not promote rutin cutaneous penetration through C. durissus skin, retaining 0.931 +/- 0.0391 mu g rutin/mg shed snake skin. The referred formulation was chemically stable for 30 days after stored at 25.0 +/- 2.0 degrees C, 5.0 +/- 0.5 degrees C and 45.0 +/- 0.5 degrees C. In conclusion, it has not been verified the active cutaneous penetration through the model biomembrane, but only its retention on the Crotalus durissus stratum corneum, condition considered stable for 30 days.

CR3 (CD11b/CD18) activation of nasal neutrophils: a measure of upper airway endotoxin exposure

Inhaled endotoxin (lipopolysaccharide, LPS) initiates an inflammatory response and leads to the expression of CR3 (CD11b/CD18) receptors on polymorphonuclear leukocytes (PMNs). We determined if PMN activation in nasal lavage fluid (NLF) is a possible biomarker of occupational endotoxin exposure. Seven subjects exposed to endotoxin provided NLF samples that were split into three aliquots (negative control - 1 M nicotinamide; sham; positive control - 11 eta g of exogenous LPS) and PMN activation was measured using a chemiluminometer. Differences in mean PMN activation were apparent, negative control: 548 +/- 15.65 RLU 100 mu l(-1); sham: 11469 +/- 2582 RLU 100 mu l(-1); positive control: 42026 +/- 16659 RLU 100 mu l (n = 7; p < 0.05). This technique shows promise as a diagnostic method for measuring upper airway LPS exposure.

Resonant Wave Interactions in the Equatorial Waveguide

Weakly nonlinear interactions among equatorial waves have been explored in this paper using the adiabatic version of the equatorial beta-plane primitive equations in isobaric coordinates. Assuming rigid lid vertical boundary conditions, the conditions imposed at the surface and at the top of the troposphere were expanded in a Taylor series around two isobaric surfaces in an approach similar to that used in the theory of surface-gravity waves in deep water and capillary-gravity waves. By adopting the asymptotic method of multiple time scales, the equatorial Rossby, mixed Rossby-gravity, inertio-gravity, and Kelvin waves, as well as their vertical structures, were obtained as leading-order solutions. These waves were shown to interact resonantly in a triad configuration at the O(epsilon) approximation. The resonant triads whose wave components satisfy a resonance condition for their vertical structures were found to have the most significant interactions, although this condition is not excluding, unlike the resonant conditions for the zonal wavenumbers and meridional modes. Thus, the analysis has focused on such resonant triads. In general, it was found that for these resonant triads satisfying the resonance condition in the vertical direction, the wave with the highest absolute frequency always acts as an energy source (or sink) for the remaining triad components, as usually occurs in several other physical problems in fluid dynamics. In addition, the zonally symmetric geostrophic modes act as catalyst modes for the energy exchanges between two dispersive waves in a resonant triad. The integration of the reduced asymptotic equations for a single resonant triad shows that, for the initial mode amplitudes characterizing realistic magnitudes of atmospheric flow perturbations, the modes in general exchange energy on low-frequency (intraseasonal and/or even longer) time scales, with the interaction period being dependent upon the initial mode amplitudes. Potential future applications of the present theory to the real atmosphere with the inclusion of diabatic forcing, dissipation, and a more realistic background state are also discussed.

Propriedades dinâmicas de fluidos por simulação computacional: métodos híbridos atomístico-contínuo

Computational methods for the calculation of dynamical properties of fluids might consider the system as a continuum or as an assembly of molecules. Molecular dynamics (MD) simulation includes molecular resolution, whereas computational fluid dynamics (CFD) considers the fluid as a continuum. This work provides a review of hybrid methods MD/CFD recently proposed in the literature. Theoretical foundations, basic approaches of computational methods, and dynamical properties typically calculated by MD and CFD are first presented in order to appreciate the similarities and differences between these two methods. Then, methods for coupling MD and CFD, and applications of hybrid simulations MD/CFD, are presented.

Padronização do eletrorretinograma de campo total em cães

Eletrorretinograma (ERG) é o meio diagnóstico objetivo e não-invasivo para avaliar a função retiniana e detectar precocemente, em várias espécies, lesões nas suas camadas mais externas. As indicações mais comuns para ERG em cães são: avaliação pré-cirúrgica de pacientes com catarata, caracterização de distúrbios que causam cegueira, além de servir como importante modelo para o estudo da distrofia retiniana que acomete o homem. Vários são os fatores que podem alterar o ERG tais como: eletrorretinógrafo, fonte de estimulação luminosa, tipo do eletrodo, tempo de adaptação ao escuro, tamanho pupilar, opacidade de meios e protocolo de sedação ou anestesia; além da espécie, raça e idade. Objetivou-se com este estudo padronizar o ERG para cães submetidos à sedação, seguindo o protocolo da International Society for Clinical Electrophysiology of Vision (ISCEV), utilizando Ganzfeld e eletrodos Burian Allen. Foram realizados 233 eletrorretinogramas em cães, 147 fêmeas e 86 machos, com idades entre um e 14 anos. Dos 233 cães examinados, 100 apresentavam catarata em diferentes estágios de maturação, 72 eram diabéticos e apresentavam catarata madura ou hipermadura, 26 apresentaram eletrorretinograma compatível com degeneração retiniana progressiva, três apresentaram eletrorretinograma compatível com síndrome da degeneração retiniana adquirida subitamente e 32 não apresentaram lesão retiniana capaz de atenuar as respostas do ERG, sendo considerados normais quanto à função retiniana. A sedação foi capaz de produzir boa imobilização do paciente sem rotacionar o bulbo ocular, permitindo adequada estimulação retiniana bilateralmente, com auxílio do Ganzfeld. O sistema eletrodiagnóstico Veris registrou com sucesso e simultaneamente de ambos os olhos, as cinco respostas preconizadas pela ISCEV. Como o ERG de campo total tornou-se exame fundamental na rotina oftalmológica, sua padronização é indispensável quando se objetiva comparar resultados de laboratórios distintos. A confiabilidade e reprodutibilidade deste protocolo foi demonstrada com a obtenção de registros de ótima qualidade utilizando protocolo padrão da ISCEV, eletrorretinógrafo Veris, Ganzfeld e eletrodos Burian Allen nos cães submetidos à sedação.

Comparative study of the hygienic behavior of Carniolan and Africanized honey bees directed towards grouped versus isolated dead brood cells

In Apis mellifera, hygienic behavior involves recognition and removal of sick, damaged or dead brood from capped cells. We investigated whether bees react in the same way to grouped versus isolated damaged capped brood cells. Three colonies of wild-type Africanized honey bees and three colonies of Carniolan honey bees were used for this investigation. Capped worker brood cells aged 12 to 14 days old were perforated with the pin-killing method. After making holes in the brood cells, the combs were placed back into the hives; 24 h later the number of cleaned cells was recorded in areas with pin-killed and control brood cells. Four repetitions were made in each colony. Isolated cells were more frequently cleaned than grouped cells, though variance analysis showed no significant difference (P = 0.1421). Carniolan bees also were somewhat, though not significantly more hygienic than Africanized honey bees (P = 0.0840). We conclude that honey bees can detect and remove both isolated and grouped dead brood. The tendency towards greater hygienic efficiency directed towards grouped pin-killed brood may be a consequence of a greater concentration of volatiles emanating from the wounds in the dead pupae.