Interpolation estimate for a finite-element space with embedded discontinuities

We consider a recently proposed finite-element space that consists of piecewise affine functions with discontinuities across a smooth given interface Γ (a curve in two dimensions, a surface in three dimensions). Contrary to existing extended finite element methodologies, the space is a variant of the standard conforming Formula space that can be implemented element by element. Further, it neither introduces new unknowns nor deteriorates the sparsity structure. It is proved that, for u arbitrary in Formula, the interpolant Formula defined by this new space satisfies Graphic where h is the mesh size, Formula is the domain, Formula, Formula, Formula and standard notation has been adopted for the function spaces. This result proves the good approximation properties of the finite-element space as compared to any space consisting of functions that are continuous across Γ, which would yield an error in the Formula-norm of order Graphic. These properties make this space especially attractive for approximating the pressure in problems with surface tension or other immersed interfaces that lead to discontinuities in the pressure field. Furthermore, the result still holds for interfaces that end within the domain, as happens for example in cracked domains.

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.

Processing of high resolution magic angle spinning spectra of breast cancer cells by the filter diagonalization method

Proton nuclear magnetic resonance (H-1 NMR) spectroscopy for detection of biochemical changes in biological samples is a successful technique. However, the achieved NMR resolution is not sufficiently high when the analysis is performed with intact cells. To improve spectral resolution, high resolution magic angle spinning (HR-MAS) is used and the broad signals are separated by a T-2 filter based on the CPMG pulse sequence. Additionally, HR-MAS experiments with a T-2 filter are preceded by a water suppression procedure. The goal of this work is to demonstrate that the experimental procedures of water suppression and T-2 or diffusing filters are unnecessary steps when the filter diagonalization method (FDM) is used to process the time domain HR-MAS signals. Manipulation of the FDM results, represented as a tabular list of peak positions, widths, amplitudes and phases, allows the removal of water signals without the disturbing overlapping or nearby signals. Additionally, the FDM can also be used for phase correction and noise suppression, and to discriminate between sharp and broad lines. Results demonstrate the applicability of the FDM post-acquisition processing to obtain high quality HR-MAS spectra of heterogeneous biological materials.

On response time reduction of electrothermomechanical MEMS using topology optimization

Electrothermomechanical MEMS are essentially microactuators that operate based on the thermoelastic effect induced by the Joule heating of the structure. They can be easily fabricated and require relatively low excitation voltages. However, the actuation time of an electrothermomechanical microdevice is higher than the actuation times related to electrostatic and piezoelectric actuation principles. Thus, in this research, we propose an optimization framework based on the topology optimization method applied to transient problems, to design electrothermomechanical microactuators for response time reduction. The objective is to maximize the integral of the output displacement of the actuator, which is a function of time. The finite element equations that govern the time response of the actuators are provided. Furthermore, the Solid Isotropic Material with Penalization model and Sequential Linear Programming are employed. Finally, a smoothing filter is implemented to control the solution. Results aiming at two distinct applications suggest the proposed approach can provide more than 50% faster actuators. (C) 2012 Elsevier B.V. All rights reserved.

Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region

Abstract Background The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. Results We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. Conclusion The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

Finite element analysis of stress distribution in anchor teeth in surgically assisted rapid palatal expansion

The treatment of a transverse maxillary deficiency in skeletally mature individuals should include surgically assisted rapid palatal expansion. This study evaluated the distribution of stresses that affect the expander's anchor teeth using finite element analysis when the osteotomy is varied. Five virtual models were built and the surgically assisted rapid palatal expansion was simulated. Results showed tension on the lingual face of the teeth and alveolar bone, and compression on the buccal side of the alveolar bone. The subtotal Le Fort I osteotomy combined with intermaxillary suture osteotomy seemed to reduce the dissipation of tensions. Therefore, subtotal Le Fort I osteotomy without a step in the zygomaticomaxillary buttress, combined with intermaxillary suture osteotomy and pterygomaxillary disjunction may be the osteotomy of choice to reduce tensions on anchor teeth, which tend to move mesiobuccally (premolar) and distobuccally (molar)

Rainwater major and trace element contents in Southeastern Brazil: an assessment of a sugar cane region in dry and wet period

The purpose of this study was to assess the composition of the rainwater in Araraquara City, Brazil, a region strongly influenced by pre-harvest burning of sugar cane crops. Chemical and mineralogical variables were measured in rainwater collected during the harvest, dry period of 2009 and the non-harvest, wet period of 2010. Ca2+ and NH4+ were responsible for 55% of cations and NO3- for 45% of anions in rainwater. Al and Fe along with K were the most abundant among trace elements in both soluble and insoluble fractions. High volume weighted mean concentration (VWM) for most of the analyzed species were observed in the harvest, dry period, mainly due to agricultural activities and meteorological conditions. The chemistry of the Araraquara rainwater and principal component analysis (PCA) quantification clearly indicate the concurrence of a diversity of sources from natural to anthropogenic especially related to agricultural activities.

Morse taper implants at different bone levels: a finite element analysis of stress distribution

AIM: To explore the biomechanical effects of the different implantation bone levels of Morse taper implants, employing a finite element analysis (FEA). METHODS: Dental implants (TitamaxCM) with 4x13 mm and 4x11 mm, and their respective abutments with 3.5 mm height, simulating a screwed premolar metal-ceramic crown, had their design performed using the software AnsysWorkbench 10.0. They were positioned in bone blocks, covered by 2.5 mm thickness of mucosa. The cortical bone was designed with 1.5 mm thickness and the trabecular bone completed the bone block. Four groups were formed: group 11CBL (11 mm implant length on cortical bone level), group 11TBL (11 mm implant length on trabecular bone level), group 13CBL (13mm implant length on cortical bone level) and group 13TBL (13 mm implant length on trabecular bone level). Oblique 200 N loads were applied. Von Mises equivalent stresses in cortical and trabecular bones were evaluated with the same design program. RESULTS: The results were shown qualitatively and quantitatively by standard scales for each type of bone. By the results obtained, it can be suggested that positioning the implant completely in trabecular bone brings harm with respect to the generated stresses. Its implantation in the cortical bone has advantages with respect to better anchoring and locking, reflecting a better dissipation of the stresses along the implant/bone interfaces. In addition, the search for anchoring the implant in its apical region in cortical bone is of great value to improve stabilization and consequently better stress distribution. CONCLUSIONS: The implant position slightly below the bone in relation to the bone crest brings advantages as the best long-term predictability with respect to the expected neck bone loss.

Numeric reconstruction of cytoskeleton with finite element method and topology optimization method

The importance of mechanical aspects related to cell activity and its environment is becoming more evident due to their influence in stem cell differentiation and in the development of diseases such as atherosclerosis. The mechanical tension homeostasis is related to normal tissue behavior and its lack may be related to the formation of cancer, which shows a higher mechanical tension. Due to the complexity of cellular activity, the application of simplified models may elucidate which factors are really essential and which have a marginal effect. The development of a systematic method to reconstruct the elements involved in the perception of mechanical aspects by the cell may accelerate substantially the validation of these models. This work proposes the development of a routine capable of reconstructing the topology of focal adhesions and the actomyosin portion of the cytoskeleton from the displacement field generated by the cell on a flexible substrate. Another way to think of this problem is to develop an algorithm to reconstruct the forces applied by the cell from the measurements of the substrate displacement, which would be characterized as an inverse problem. For these kind of problems, the Topology Optimization Method (TOM) is suitable to find a solution. TOM is consisted of an iterative application of an optimization method and an analysis method to obtain an optimal distribution of material in a fixed domain. One way to experimentally obtain the substrate displacement is through Traction Force Microscopy (TFM), which also provides the forces applied by the cell. Along with systematically generating the distributions of focal adhesion and actin-myosin for the validation of simplified models, the algorithm also represents a complementary and more phenomenological approach to TFM. As a first approximation, actin fibers and flexible substrate are represented through two-dimensional linear Finite Element Method. Actin contraction is modeled as an initial stress of the FEM elements. Focal adhesions connecting actin and substrate are represented by springs. The algorithm was applied to data obtained from experiments regarding cytoskeletal prestress and micropatterning, comparing the numerical results to the experimental ones

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.

Strumenti per la valutazione dello stato conservativo e per la manutenzione programmata dei ponti della Provincia di Bologna

La tesi tratta di strumenti finalizzati alla valutazione dello stato conservativo e di supporto all'attività di manutenzione dei ponti, dai più generali Bridge Management Systems ai Sistemi di Valutazione Numerica della Condizione strutturale. Viene proposto uno strumento originale con cui classificare i ponti attraverso un Indice di Valutazione Complessiva e grazie ad esso stabilire le priorità d'intervento. Si tara lo strumento sul caso pratico di alcuni ponti della Provincia di Bologna. Su un ponte in particolare viene realizzato un approfondimento specifico sulla determinazione approssimata dei periodi propri delle strutture da ponte. Si effettua un confronto dei risultati di alcune modellazioni semplificate in riferimento a modellazioni dettagliate e risultati sperimentali.

Expression of the repressor element-1 silencing transcription factor (REST) is regulated by IGF-I and PKC in human neuroblastoma cells

The repressor element 1-silencing transcription factor (REST) was first identified as a protein that binds to a 21-bp DNA sequence element (known as repressor element 1 (RE1)) resulting in transcriptional repression of the neural-specific genes [Chong et al., 1995; Schoenherr and Anderson, 1995]. The original proposed role for REST was that of a factor responsible for restricting neuronal gene expression to the nervous system by silencing expression of these genes in non-neuronal cells. Although it was initially thought to repress neuronal genes in non-neuronal cells, the role of REST is complex and tissue dependent. In this study I investigated any role played by REST in the induction and patterning of differentiation of SH-SY5Y human neuroblastoma cells exposed to IGF-I. and phorbol 12- myristate 13-acetate (PMA) To down-regulate REST expression we developed an antisense (AS) strategy based on the use of phosphorothioate oligonucleotides (ODNs). In order to evaluate REST mRNA levels, we developed a real-time PCR technique and REST protein levels were evaluated by western blotting. Results showed that nuclear REST is increased in SH-SY5Y neuroblastoma cells cultured in SFM and exposed to IGF-I for 2-days and it then declines in 5-day-treated cells concomitant with a progressive neurite extension. Also the phorbol ester PMA was able to increase nuclear REST levels after 3-days treatment concomitant to neuronal differentiation of neuroblastoma cells, whereas, at later stages, it is down-regulated. Supporting these data, the exposure to PKC inhibitors (GF10923X and Gö6976) and PMA (16nM) reverted the effects observed with PMA alone. REST levels were related to morphological differentiation, expression of growth coneassociated protein 43 (GAP-43; a gene not regulated by REST) and of synapsin I and βIII tubulin (genes regulated by REST), proteins involved in the early stage of neuronal development. We observed that differentiation of SH-SY5Y cells by IGF-I and PMA was accompanied by a significant increase of these neuronal markers, an effect that was concomitant with REST decrease. In order to relate the decreased REST expression with a progressive neurite extension, I investigated any possible involvement of the ubiquitin–proteasome system (UPS), a multienzymatic pathway which degrades polyubiquinated soluble cytoplasmic proteins [Pickart and Cohen, 2004]. For this purpose, SH-SY5Y cells are concomitantly exposed to PMA and the proteasome inhibitor MG132. In SH-SY5Y exposed to PMA and MG 132, we observed an inverse pattern of expression of synapsin I and β- tubulin III, two neuronal differentiation markers regulated by REST. Their cytoplasmic levels are reduced when compared to cells exposed to PMA alone, as a consequence of the increase of REST expression by proteasome inhibitor. The majority of proteasome substrates identified to date are marked for degradation by polyubiquitinylation; however, exceptions to this principle, are well documented [Hoyt and Coffino, 2004]. Interestingly, REST degradation seems to be completely ubiquitin-independent. The expression pattern of REST could be consistent with the theory that, during early neuronal differentiation induced by IGF-I and PKC, it may help to repress the expression of several genes not yet required by the differentiation program and then it declines later. Interestingly, the observation that REST expression is progressively reduced in parallel with cell proliferation seems to indicate that the role of this transcription factor could also be related to cell survival or to counteract apotosis events [Lawinger et al., 2000] although, as shown by AS-ODN experiments, it does not seem to be directly involved in cell proliferation. Therefore, the decline of REST expression is a comparatively later event during maturation of neuroroblasts in vitro. Thus, we propose that REST is regulated by growth factors, like IGF-I, and PKC activators in a time-dependent manner: it is elevated during early steps of neural induction and could contribute to down-regulate genes not yet required by the differentiation program while it declines later for the acquisition of neural phenotypes, concomitantly with a progressive neurite extension. This later decline is regulated by the proteasome system activation in an ubiquitin-indipendent way and adds more evidences to the hypothesis that REST down-regulation contributes to differentiation and arrest of proliferation of neuroblastoma cells. Finally, the glycosylation pattern of the REST protein was analysed, moving from the observation that the molecular weight calculated on REST sequence is about 116 kDa but using western blotting this transcription factor appears to have distinct apparent molecular weight (see Table 1.1): this difference could be explained by post-translational modifications of the proteins, like glycosylation. In fact recently, several studies underlined the importance of O-glycosylation in modulating transcriptional silencing, protein phosphorylation, protein degradation by proteasome and protein–protein interactions [Julenius et al., 2005; Zachara and Hart, 2006]. Deglycosilating analysis showed that REST protein in SH-SY5Y and HEK293 cells is Oglycosylated and not N-glycosylated. Moreover, using several combination of deglycosilating enzymes it is possible to hypothesize the presence of Gal-β(1-3)-GalNAc residues on the endogenous REST, while β(1-4)-linked galactose residues may be present on recombinant REST protein expressed in HEK293 cells. However, the O-glycosylation process produces an immense multiplicity of chemical structures and monosaccharides must be sequentially hydrolyzed by a series of exoglycosidase. Further experiments are needed to characterize all the post-translational modification of the transcription factor REST.

3-D Boundary element-finite element method for the dynamic analysis of piled buildings.

[EN]A boundary element-finite element model is presented for the three-dimensional dynamic analysis of piled buildings in the frequency domain. Piles are modelled as compressible Euler-Bernoulli beams founded on a linear, isotropic, viscoelastic, zoned-homogeneous, unbounded layered soil, while multi-storey buildings are assumed to be comprised of vertical compressible piers and rigid slabs.

Finite element simulation of a local scale air quality model over complex terrain

[EN]In this paper we propose a finite element method approach for modelling the air quality in a local scale over complex terrain. The area of interest is up to tens of kilometres and it includes pollutant sources. The proposed methodology involves the generation of an adaptive tetrahedral mesh, the computation of an ambient wind field, the inclusion of the plume rise effect in the wind field, and the simulation of transport and reaction of pollutants. We apply our methodology to simulate a fictitious pollution episode in La Palma island (Canary Island, Spain)...