Aerodynamic study of three-dimensional larynx models using finite element methods

The airflow velocities and pressures are calculated from a three-dimensional model of the human larynx by using the finite element method. The laryngeal airflow is assumed to be incompressible, isothermal, steady, and created by fixed pressure drops. The influence of different laryngeal profiles (convergent, parallel, and divergent), glottal area, and dimensions of false vocal folds in the airflow are investigated. The results indicate that vertical and horizontal phase differences in the laryngeal tissue movements are influenced by the nonlinear pressure distribution across the glottal channel, and the glottal entrance shape influences the air pressure distribution inside the glottis. Additionally, the false vocal folds increase the glottal duct pressure drop by creating a new constricted channel in the larynx, and alter the airflow vortexes formed after the true vocal folds. (C) 2007 Elsevier Ltd. All rights reserved.

Efficient Modeling of Thin Wires in a Lossy Medium by Finite Elements Applied to Grounding Systems

A procedure is proposed to accurately model thin wires in lossy media by finite element analysis. It is based on the determination of a suitable element width in the vicinity of the wire, which strongly depends on the wire radius to yield accurate results. The approach is well adapted to the analysis of grounding systems. The numerical results of the application of finite element analysis with the suitably chosen element width are compared with both analytical results and those computed by a commercial package for the analysis of grounding systems, showing very good agreement.

Shell curvature as an initial deformation: A geometrically exact finite element approach

An alternative approach for the analysis of arbitrarily curved shells is developed in this paper based on the idea of initial deformations. By `alternative` we mean that neither differential geometry nor the concept of degeneration is invoked here to describe the shell surface. We begin with a flat reference configuration for the shell mid-surface, after which the initial (curved) geometry is mapped as a stress-free deformation from the plane position. The actual motion of the shell takes place only after this initial mapping. In contrast to classical works in the literature, this strategy enables the use of only orthogonal frames within the theory and therefore objects such as Christoffel symbols, the second fundamental form or three-dimensional degenerated solids do not enter the formulation. Furthermore, the issue of physical components of tensors does not appear. Another important aspect (but not exclusive of our scheme) is the possibility to describe exactly the initial geometry. The model is kinematically exact, encompasses finite strains in a totally consistent manner and is here discretized under the light of the finite element method (although implementation via mesh-free techniques is also possible). Assessment is made by means of several numerical simulations. Copyright (C) 2009 John Wiley & Sons, Ltd.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids. (C) 2010 Elsevier B.V. All rights reserved.

Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

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.

A 2-D Delaunay refinement algorithm using an initial prerefinement from the boundary mesh

An improvement to the quality bidimensional Delaunay mesh generation algorithm, which combines the mesh refinement algorithms strategy of Ruppert and Shewchuk is proposed in this research. The developed technique uses diametral lenses criterion, introduced by L. P. Chew, with the purpose of eliminating the extremely obtuse triangles in the boundary mesh. This method splits the boundary segment and obtains an initial prerefinement, and thus reducing the number of necessary iterations to generate a high quality sequential triangulation. Moreover, it decreases the intensity of the communication and synchronization between subdomains in parallel mesh refinement.

An EFG method for the nonlinear analysis of plates undergoing arbitrarily large deformations

The applicability of a meshfree approximation method, namely the EFG method, on fully geometrically exact analysis of plates is investigated. Based on a unified nonlinear theory of plates, which allows for arbitrarily large rotations and displacements, a Galerkin approximation via MLS functions is settled. A hybrid method of analysis is proposed, where the solution is obtained by the independent approximation of the generalized internal displacement fields and the generalized boundary tractions. A consistent linearization procedure is performed, resulting in a semi-definite generalized tangent stiffness matrix which, for hyperelastic materials and conservative loadings, is always symmetric (even for configurations far from the generalized equilibrium trajectory). Besides the total Lagrangian formulation, an updated version is also presented, which enables the treatment of rotations beyond the parameterization limit. An extension of the arc-length method that includes the generalized domain displacement fields, the generalized boundary tractions and the load parameter in the constraint equation of the hyper-ellipsis is proposed to solve the resulting nonlinear problem. Extending the hybrid-displacement formulation, a multi-region decomposition is proposed to handle complex geometries. A criterium for the classification of the equilibrium`s stability, based on the Bordered-Hessian matrix analysis, is suggested. Several numerical examples are presented, illustrating the effectiveness of the method. Differently from the standard finite element methods (FEM), the resulting solutions are (arbitrary) smooth generalized displacement and stress fields. (c) 2007 Elsevier Ltd. All rights reserved.

Numerical prediction of gas concentrations and fluctuations above a triangular hill within a turbulent boundary layer

The objective of the present work is to propose a numerical and statistical approach, using computational fluid dynamics, for the study of the atmospheric pollutant dispersion. Modifications in the standard k-epsilon turbulence model and additional equations for the calculation of the variance of concentration are introduced to enhance the prediction of the flow field and scalar quantities. The flow field, the mean concentration and the variance of a flow over a two-dimensional triangular hill, with a finite-size point pollutant source, are calculated by a finite volume code and compared with published experimental results. A modified low Reynolds k-epsilon turbulence model was employed in this work, using the constant of the k-epsilon model C(mu)=0.03 to take into account the inactive atmospheric turbulence. The numerical results for the velocity profiles and the position of the reattachment point are in good agreement with the experimental results. The results for the mean and the variance of the concentration are also in good agreement with experimental results from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Topology optimization for designing strain-gauge load cells

Load cells are used extensively in engineering fields. This paper describes a novel structural optimization method for single- and multi-axis load cell structures. First, we briefly explain the topology optimization method that uses the solid isotropic material with penalization (SIMP) method. Next, we clarify the mechanical requirements and design specifications of the single- and multi-axis load cell structures, which are formulated as an objective function. In the case of multi-axis load cell structures, a methodology based on singular value decomposition is used. The sensitivities of the objective function with respect to the design variables are then formulated. On the basis of these formulations, an optimization algorithm is constructed using finite element methods and the method of moving asymptotes (MMA). Finally, we examine the characteristics of the optimization formulations and the resultant optimal configurations. We confirm the usefulness of our proposed methodology for the optimization of single- and multi-axis load cell structures.

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.

Memory with emotional content, brain amygdala and Alzheimer`s disease

Objectives - A highly adaptive aspect of human memory is the enhancement of explicit, consciously accessible memory by emotional stimuli. We studied the performance of Alzheimer`s disease (AD) patients and elderly controls using a memory battery with emotional content, and we correlated these results with the amygdala and hippocampus volume. Methods - Twenty controls and 20 early AD patients were subjected to the International Affective Picture System (IAPS) and to magnetic resonance imaging-based volumetric measurements of the medial temporal lobe structures. Results - The results show that excluding control group subjects with 5 or more years of schooling, both groups showed improvement with pleasant or unpleasant figures for the IAPS in an immediate free recall test. Likewise, in a delayed free recall test, both the controls and the AD group showed improvement for pleasant pictures, when education factor was not controlled. The AD group showed improvement in the immediate and delayed free recall test proportional to the medial temporal lobe structures, with no significant clinical correlation between affective valence and amygdala volume. Conclusion - AD patients can correctly identify emotions, at least at this early stage, but this does not improve their memory performance.

Brain structure and symptom dimension relationships in obsessive-compulsive disorder: A voxel-based morphometry study

Background: Obsessive-compulsive disorder (OCD) is a clinically heterogenous disorder characterized by temporally stable symptom dimensions. Past inconsistent results from structural neuroimaging studies of OCD may have resulted from the effects of these specific symptom dimensions as well as other socio-demographic and clinical variables upon gray matter (GM) volume. Methods: GM volume was measured in 25 adult OCD patients and 20 adult healthy controls using voxel-based morphometry (VBM), controlling for age and total brain GM volume. Univariate and multivariate regression analyses were carried out between regions of GM difference and age, age of onset, medication load, OCD severity, depression severity, and separate symptom dimension scores. Results: Significant GM volumetric differences in OCD patients relative to controls were found in dorsal cortical regions, including bilateral BA6, BA46, BA9 and right BA8 (controls > patients), and bilateral midbrain (patients > controls). Stepwise regression analyses revealed highly significant relationships between greater total OCD symptom severity and smaller GM volumes in dorsal cortical regions and larger GM volumes in bilateral midbrain. Greater age was independently associated with smaller GM volumes in right BA6, left BA9, left BA46 and larger GM volumes in right midbrain. Greater washing symptom severity was independently associated with smaller GM volume in right BA6, while there was a trend association between greater hoarding symptom severity and lower GM volume in left BA6. Limitations: The sample was relatively small to examine the relationship between symptom scores and GM volumes. Multiple patients were taking medication and had comorbid disorders. Conclusions: These analyses suggest dorsal prefrontal cortical and bilateral midbrain GM abnormalities in OCD that appear to be primarily driven by the effects of total OCD symptom severity. The results regarding the relationship between GM volumes and symptom dimension scores require examination in larger samples. (C) 2008 Elsevier B.V. All rights reserved.

Effects of changing blood viscosity and heart rate on vena contracta width as evaluated by color Doppler flow mapping. An in vitro study with a pulsatile flow model

Background: There is only limited knowledge on how the quantification of valvular regurgitation by color Doppler is affected by changing blood viscosity. This study was designed to evaluate the effect of changing blood viscosity on the vena contracta width using an in vitro model of valvular insufficiency capable of providing ample variation in the rate and stroke volume. Methods: We constructed a pulsatile flow model filled with human blood at varying hematocrit (15%, 35%, and 55%) and corresponding blood viscosity (blood/water viscosity: 2.6, 4.8, 9.1) levels in which jets were driven through a known orifice (7 mm(2)) into a 110 mL compliant receiving chamber (compliance: 2.2 mL/mmHg) by a pulsatile pump. In addition, we used variable pump stroke volumes (5, 7.5, and 10 mL) and rates (40, 60, and 80 ppm). Vena contracta region was imaged using a 3.5 MHz transducer. Pressure and volume in the flow model were kept constant during each experimental condition, as well as ultrasound settings. Results: Blood viscosity variation in the experimental range did not induce significant changes in vena contracta dimensions. Also, vena contracta width did not change from normal to low hematocrit and viscosity levels. A very modest increase only in vena contracta dimension was observed at very high level of blood viscosity when hematocrit was set to 55% . Pump rate, in the evaluated range, did not influence vena contracta width. These results in controlled experimental settings suggest that the vena contracta is an accurate quantitative method for quantifying valvular regurgitation even when this condition is associated with anemia, a frequent finding in patients with valvular heart disease.

Mapping the cytoskeletal prestress

Park CY, Tambe D, Alencar AM, Trepat X, Zhou EH, Millet E, Butler JP, Fredberg JJ. Mapping the cytoskeletal prestress. Am J Physiol Cell Physiol 298: C1245-C1252, 2010. First published February 17, 2010; doi: 10.1152/ajpcell.00417.2009.-Cell mechanical properties on a whole cell basis have been widely studied, whereas local intracellular variations have been less well characterized and are poorly understood. To fill this gap, here we provide detailed intracellular maps of regional cytoskeleton (CSK) stiffness, loss tangent, and rate of structural rearrangements, as well as their relationships to the underlying regional F-actin density and the local cytoskeletal prestress. In the human airway smooth muscle cell, we used micropatterning to minimize geometric variation. We measured the local cell stiffness and loss tangent with optical magnetic twisting cytometry and the local rate of CSK remodeling with spontaneous displacements of a CSK-bound bead. We also measured traction distributions with traction microscopy and cell geometry with atomic force microscopy. On the basis of these experimental observations, we used finite element methods to map for the first time the regional distribution of intracellular prestress. Compared with the cell center or edges, cell corners were systematically stiffer and more fluidlike and supported higher traction forces, and at the same time had slower remodeling dynamics. Local remodeling dynamics had a close inverse relationship with local cell stiffness. The principal finding, however, is that systematic regional variations of CSK stiffness correlated only poorly with regional F-actin density but strongly and linearly with the regional prestress. Taken together, these findings in the intact cell comprise the most comprehensive characterization to date of regional variations of cytoskeletal mechanical properties and their determinants.