Finite element methods for the Stokes system based on a Zienkiewicz type N-simplex

Hermite interpolation is increasingly showing to be a powerful numerical solution tool, as applied to different kinds of second order boundary value problems. In this work we present two Hermite finite element methods to solve viscous incompressible flows problems, in both two- and three-dimension space. In the two-dimensional case we use the Zienkiewicz triangle to represent the velocity field, and in the three-dimensional case an extension of this element to tetrahedra, still called a Zienkiewicz element. Taking as a model the Stokes system, the pressure is approximated with continuous functions, either piecewise linear or piecewise quadratic, according to the version of the Zienkiewicz element in use, that is, with either incomplete or complete cubics. The methods employ both the standard Galerkin or the Petrov–Galerkin formulation first proposed in Hughes et al. (1986) [18], based on the addition of a balance of force term. A priori error analyses point to optimal convergence rates for the PG approach, and for the Galerkin formulation too, at least in some particular cases. From the point of view of both accuracy and the global number of degrees of freedom, the new methods are shown to have a favorable cost-benefit ratio, as compared to velocity Lagrange finite elements of the same order, especially if the Galerkin approach is employed.

Optimal reactive power flow via the modified barrier Lagrangian function approach

A new approach called the Modified Barrier Lagrangian Function (MBLF) to solve the Optimal Reactive Power Flow problem is presented. In this approach, the inequality constraints are treated by the Modified Barrier Function (MBF) method, which has a finite convergence property: i.e. the optimal solution in the MBF method can actually be in the bound of the feasible set. Hence, the inequality constraints can be precisely equal to zero. Another property of the MBF method is that the barrier parameter does not need to be driven to zero to attain the solution. Therefore, the conditioning of the involved Hessian matrix is greatly enhanced. In order to show this, a comparative analysis of the numeric conditioning of the Hessian matrix of the MBLF approach, by the decomposition in singular values, is carried out. The feasibility of the proposed approach is also demonstrated with comparative tests to Interior Point Method (IPM) using various IEEE test systems and two networks derived from Brazilian generation/transmission system. The results show that the MBLF method is computationally more attractive than the IPM in terms of speed, number of iterations and numerical conditioning. (C) 2011 Elsevier B.V. All rights reserved.

GENERALIZED FINITE ELEMENT METHOD ON NONCONVENTIONAL HYBRID-MIXED FORMULATION

The generalized finite element method (GFEM) is applied to a nonconventional hybrid-mixed stress formulation (HMSF) for plane analysis. In the HMSF, three approximation fields are involved: stresses and displacements in the domain and displacement fields on the static boundary. The GFEM-HMSF shape functions are then generated by the product of a partition of unity associated to each field and the polynomials enrichment functions. In principle, the enrichment can be conducted independently over each of the HMSF approximation fields. However, stability and convergence features of the resulting numerical method can be affected mainly by spurious modes generated when enrichment is arbitrarily applied to the displacement fields. With the aim to efficiently explore the enrichment possibilities, an extension to GFEM-HMSF of the conventional Zienkiewicz-Patch-Test is proposed as a necessary condition to ensure numerical stability. Finally, once the extended Patch-Test is satisfied, some numerical analyses focusing on the selective enrichment over distorted meshes formed by bilinear quadrilateral finite elements are presented, thus showing the performance of the GFEM-HMSF combination.

Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.

Analysis of elastic functionally graded materials under large displacements via high-order tetrahedral elements

The purpose of this study is to present a position based tetrahedral finite element method of any order to accurately predict the mechanical behavior of solids constituted by functionally graded elastic materials and subjected to large displacements. The application of high-order elements makes it possible to overcome the volumetric and shear locking that appears in usual homogeneous isotropic situations or even in non-homogeneous cases developing small or large displacements. The use of parallel processing to improve the computational efficiency, allows employing high-order elements instead of low-order ones with reduced integration techniques or strain enhancements. The Green-Lagrange strain is adopted and the constitutive relation is the functionally graded Saint Venant-Kirchhoff law. The equilibrium is achieved by the minimum total potential energy principle. Examples of large displacement problems are presented and results confirm the locking free behavior of high-order elements for non-homogeneous materials. (C) 2011 Elsevier B.V. All rights reserved.

A generalized finite element method with global-local enrichment functions for confined plasticity problems

The main feature of partition of unity methods such as the generalized or extended finite element method is their ability of utilizing a priori knowledge about the solution of a problem in the form of enrichment functions. However, analytical derivation of enrichment functions with good approximation properties is mostly limited to two-dimensional linear problems. This paper presents a procedure to numerically generate proper enrichment functions for three-dimensional problems with confined plasticity where plastic evolution is gradual. This procedure involves the solution of boundary value problems around local regions exhibiting nonlinear behavior and the enrichment of the global solution space with the local solutions through the partition of unity method framework. This approach can produce accurate nonlinear solutions with a reduced computational cost compared to standard finite element methods since computationally intensive nonlinear iterations can be performed on coarse global meshes after the creation of enrichment functions properly describing localized nonlinear behavior. Several three-dimensional nonlinear problems based on the rate-independent J (2) plasticity theory with isotropic hardening are solved using the proposed procedure to demonstrate its robustness, accuracy and computational efficiency.

Relationship between global structural parameters and Enzyme Commission hierarchy: Implications for function prediction

In protein databases there is a substantial number of proteins structurally determined but without function annotation. Understanding the relationship between function and structure can be useful to predict function on a large scale. We have analyzed the similarities in global physicochemical parameters for a set of enzymes which were classified according to the four Enzyme Commission (EC) hierarchical levels. Using relevance theory we introduced a distance between proteins in the space of physicochemical characteristics. This was done by minimizing a cost function of the metric tensor built to reflect the EC classification system. Using an unsupervised clustering method on a set of 1025 enzymes, we obtained no relevant clustering formation compatible with EC classification. The distance distributions between enzymes from the same EC group and from different EC groups were compared by histograms. Such analysis was also performed using sequence alignment similarity as a distance. Our results suggest that global structure parameters are not sufficient to segregate enzymes according to EC hierarchy. This indicates that features essential for function are rather local than global. Consequently, methods for predicting function based on global attributes should not obtain high accuracy in main EC classes prediction without relying on similarities between enzymes from training and validation datasets. Furthermore, these results are consistent with a substantial number of studies suggesting that function evolves fundamentally by recruitment, i.e., a same protein motif or fold can be used to perform different enzymatic functions and a few specific amino acids (AAs) are actually responsible for enzyme activity. These essential amino acids should belong to active sites and an effective method for predicting function should be able to recognize them. (C) 2012 Elsevier Ltd. All rights reserved.

Contribution to assessing the stiffness reduction of structural elements in the global stability analysis of precast concrete multi-storey buildings

This study deals with the reduction of the stiffness in precast concrete structural elements of multi-storey buildings to analyze global stability. Having reviewed the technical literature, this paper present indications of stiffness reduction in different codes, standards, and recommendations and compare these to the values found in the present study. The structural model analyzed in this study was constructed with finite elements using ANSYS® software. Physical Non-Linearity (PNL) was considered in relation to the diagrams M x N x 1/r, and Geometric Non-Linearity (GNL) was calculated following the Newton-Raphson method. Using a typical precast concrete structure with multiple floors and a semi-rigid beam-to-column connection, expressions for a stiffness reduction coefficient are presented. The main conclusions of the study are as follows: the reduction coefficients obtained from the diagram M x N x 1/r differ from standards that use a simplified consideration of PNL; the stiffness reduction coefficient for columns in the arrangements analyzed were approximately 0.5 to 0.6; and the variation of values found for stiffness reduction coefficient in concrete beams, which were subjected to the effects of creep with linear coefficients from 0 to 3, ranged from 0.45 to 0.2 for positive bending moments and 0.3 to 0.2 for negative bending moments.

Influence of chest wall on lung mechanics during inspiration

RATIONALE: The interaction between lungs and chest wall influences lung volume, that determines lung history during respiration cycle. In this study, the influence of chest wall mechanics on respiratory system is assessed by the evaluation of inspiration pressure-volume curve (PV curve) under three different situations: closed-chest, open-chest and isolated lung. The PV curve parameters in each situation allow us to further understand the role played by different chest wall elements in the respiratory function. Methods: Twenty-four male Wistar rats (236 ± 29 g) were used. The animals were weighted and then anesthetized with xylazine 2% (O,SmL/kg) and ketamine 10% (0,9mL/kg), exsanguinated and later tracheostomies with a metallic cannula (14 gauge).The cannula was connected to an automatic small animal insufflator. This setup was connected to a pressure transducer (32 samples/s). The 24 animals were randomly separated in three groups:(i) closed chest,(ii) open chest and (iii) isolated lung. The rats were insufflated with 20mL quasi-statically (constant speed of 0,1mUs). lnsufflated volume and measured pressure data were kept and PV curves were obtained for all animals. The PV curves were fitted (non-linear least squares) against the sigmoid equation (1) to obtain the sigmoid equation parameters (a,b,c,d). Elastance measurements were obtained from linear regression of pressure/volume measurements in a 0,8s interval before and after the calculated point. Results: The parameters a,b and c showed no significant change, but the parameter d showed a significant variation among the three groups. The initial elastance also varied between open and closed chest, indicating the need of a higher pressure for the lung expansion, as can be seen in Table 1. Conclusion: A supporting effect of the chest wall was observed at the initial moments of inspiration, observed as a higher initial elastance in open chest situations than in closed chest situations (p=0,00001). The similar initial elastance for the isolated lung and closed chest may be explained by the specific method used for the isolated lung experiment. As the isolated lung is supported by the trachea vertically, the weight of the tissue may have a similar effect of the residual negative pressure in the thorax, responsible for maintaining the residual volume.

Influence of chest wall on lung mechanics during inspiration.

RATIONALE: The interaction between lungs and chest wall influences lung volume, that determines lung history during respiration cycle. In this study, the influence of chest wall mechanics on respiratory system is assessed by the evaluation of inspiration pressure-volume curve (PV curve) under three different situations: closed-chest, open-chest and isolated lung. The PV curve parameters in each situation allow us to further understand the role played by different chest wall elements in the respiratory function. Methods: Twenty-four male Wistar rats (236 ± 29 g) were used. The animals were weighted and then anesthetized with xylazine 2% (0,5mL/kg) and ketamine 10% (0,9mL/kg), exsanguinated and later tracheostomized with a metallic cannula (14 gauge). The cannula was connected to an automatic small animal insufflator. This setup was connected to a pressure transducer (32 samples/s). The 24 animals were randomly separated in three groups: (i) closed chest, (ii) open chest and (iii) isolated lung. The rats were insufflated with 20mL quasi-statically (constant speed of 0,1mL/s). Insufflated volume and measured pressure data were kept and PV curves were obtained for all animals. The PV curves were fitted (non-linear least squares) against the sigmoid equation (1) to obtain the sigmoid equation parameters (a,b,c,d). Elastance measurements were obtained from linear regression of pressure/volume measurements in a 0,8s interval before and after the calculated point. Results: The parameters a, b and c showed no significant change, but the parameter d showed a significant variation among the three groups. The initial elastance also varied between open and closed chest, indicating the need of a higher pressure for the lung expansion, as can be seen in Table 1. Table 1: Mean and Standard Deviation of parameters obtained for each protocol. Protocol: Closed Chest – a (mL) -0.35±0.33; b (mL) 13.93±0.89; c (cm H2O) 21.28±2.37; d (cm H2O) 6.17±0.84; r²** (%) 99.4±0.14; Initial Elastance* (cm H2)/mL) 12.72±6.66; Weight (g) 232.33±5.72. Open Chest - a (mL) 0.01±0.28; b (mL) 14.79±0.54; c (cm H2O) 19.47±1.41; d (cm H2O) 3.50±0.28; r²** (%) 98.8±0.34; Initial Elastance* (cm H2)/mL) 28.68±2.36; Weight (g) 217.33±7.97. Isolated Lung - a (mL) -0.09±0.46; b (mL) 14.22±0.75; c (cm H2O) 21.76±1.43; d (cm H2O) 4.24±0.50; r²** (%) 98.9±0.19; Initial Elastance* (cm H2)/mL) 7.13±8.85; Weight (g) 224.33±16.66. * Elastance measures in the 0-0,1 mL range. ** Goodness of sigmoid fit versus measured data Conclusion: A supporting effect of the chest wall was observed at the initial moments of inspiration, observed as a higher initial elastance in open chest situations than in closed chest situations (p=0,00001). The similar initial elastance for the isolated lung and closed chest may be explained by the specific method used for the isolated lung experiment. As the isolated lung is supported by the trachea vertically, the weight of the tissue may have a similar effect of the residual negative pressure in the thorax, responsible for maintaining the residual volume.

Analysis, manufacture and characterization of Ni/Cu functionally graded structures

In this work, an experimental and numerical analysis and characterization of functionally graded structures (FGSs) is developed. Nickel (Ni) and copper (Cu) materials are used as basic materials in the numerical modeling and experimental characterization. For modeling, a MATLAB finite element code is developed, which allows simulation of harmonic and modal analysis considering the graded finite element formulation. For experimental characterization, Ni-Cu FGSs are manufactured by using spark plasma sintering technique. Hardness and Young's modulus are found by using microindentation and ultrasonic measurements, respectively. The effective gradation of Ni/Cu FGS is addressed by means of optical microscopy, energy dispersive spectrometry, scanning electron microscopy and hardness testing. For the purpose of comparing modeling and experimental results, the hardness curve, along the gradation direction, is used for identifying the gradation profile; accordingly, the experimental hardness curve is used for approximating the Young's modulus variation and the graded finite element modeling is used for verification. For the first two resonance frequency values, a difference smaller than 1% between simulated and experimental results is obtained. (C) 2012 Elsevier Ltd. All rights reserved.

Myocardial Deformation by Speckle Tracking in Severe Dilated Cardiomyopathy

Background: The high and increasing prevalence of Dilated Cardiomyopathy (DCM) represents a serious public health issue. Novel technologies have been used aiming to improve diagnosis and the therapeutic approach. In this context, speckle tracking echocardiography (STE) uses natural myocardial markers to analyze the systolic deformation of the left ventricle (LV). Objective: Measure the longitudinal transmural global strain (GS) of the LV through STE in patients with severe DCM, comparing the results with normal individuals and with echocardiographic parameters established for the analysis of LV systolic function, in order to validate the method in this population. Methods: Seventy-one patients with severe DCM (53 +/- 12 years, 72% men) and 20 controls (30 +/- 8 years, 45% men) were studied. The following variables were studied: LV volumes and ejection fraction calculated by two and three-dimensional echocardiography, Doppler parameters, Tissue Doppler Imaging systolic and diastolic LV velocities and GS obtained by STE. Results: Compared with controls, LV volumes were higher in the DCM group; however, LVEF and peak E-wave velocity were lower in the latter. The myocardial performance index was higher in the patient group. Tissue Doppler myocardial velocities (S', e', a') were significantly lower and E/e' ratio was higher in the DCM group. GS was decreased in the DCM group (-5.5% +/- 2.3%) when compared with controls (-14.0% +/- 1.8%). Conclusion: In this study, GS was significantly lower in patients with severe DCM, bringing new perspectives for therapeutic approaches in this specific population. (Arq Bras Cardiol 2012;99(3):834-842)

Structural properties and energetics of diffuse Rb-87 clusters in three-dimension

A correlated two-body basis function is used to describe the three-dimensional bosonic clusters interacting via two-body van der Waals potential. We calculate the ground state and the zero orbital angular momentum excited states for Rb-N clusters with up to N = 40. We solve the many-particle Schrodinger equation by potential harmonics expansion method, which keeps all possible two-body correlations in the calculation and determines the lowest effective many-body potential. We study energetics and structural properties for such diffuse clusters both at dimer and tuned scattering length. The motivation of the present study is to investigate the possibility of formation of N-body clusters interacting through the van der Waals interaction. We also compare the system with the well studied He, Ne, and Ar clusters. We also calculate correlation properties and observe the generalised Tjon line for large cluster. We test the validity of the shape-independent potential in the calculation of the ground state energy of such diffuse cluster. These are the first such calculations reported for Rb clusters. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4730972]

Atomically resolved study of the morphology change of InAs/GaAs quantum dot layers induced by rapid thermal annealing

The optoelectronic properties of InAs/GaAs quantum dots can be tuned by rapid thermal annealing. In this study, the morphology change of InAs/GaAs quantum dots layers induced by rapid thermal annealing was investigated at the atomic-scale by cross-sectional scanning tunneling microscopy. Finite elements calculations that model the outward relaxation of the cleaved surface were used to determine the indium composition profile of the wetting layer and the quantum dots prior and post rapid thermal annealing. The results show that the wetting layer is broadened upon annealing. This broadening could be modeled by assuming a random walk of indium atoms. Furthermore, we show that the stronger strain gradient at the location of the quantum dots enhances the intermixing. Photoluminescence measurements show a blueshift and narrowing of the photoluminescence peak. Temperature dependent photoluminescence measurements show a lower activation energy for the annealed sample. These results are in agreement with the observed change in morphology. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4770371]

Lung function monitoring in patients with duchenne muscular dystrophy on steroid therapy

Background Duchenne muscular dystrophy (DMD) is a sex-linked inherited muscle disease characterized by a progressive loss in muscle strength and respiratory muscle involvement. After 12 years of age, lung function declines at a rate of 6 % to 10.7 % per year in patients with DMD. Steroid therapy has been proposed to delay the loss of motor function and also the respiratory involvement. Method In 21 patients with DMD aged between seven and 16 years, the forced vital capacity (FVC) and the forced expiratory volume in one second (FEV1) were evaluated at three different times during a period of two years. Results We observed in this period of evaluation the maintenance of the FVC and the FEV1 in this group of patients independently of chronological age, age at onset of steroid therapy, and walking capacity. Conclusion The steroid therapy has the potential to stabilize or delay the loss of lung function in DMD patients even if they are non-ambulant or older than 10 years, and in those in whom the medication was started after 7 years of age.