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.

delta-Superderivations of Semisimple Finite-Dimensional Jordan Superalgebras

In the paper, a complete description of the delta-derivations and the delta-superderivations of semisimple finite-dimensional Jordan superalgebras over an algebraically closed field of characteristic p not equal 2 is given. In particular, new examples of nontrivial (1/2)-derivations and odd (1/2)-superderivations are given that are not operators of right multiplication by an element of the superalgebra.

Causal amplitudes in the Schwinger model at finite temperature

We show, in the imaginary time formalism, that the temperature dependent parts of all the retarded (advanced) amplitudes vanish in the Schwinger model. We trace this behavior to the CPT invariance of the theory and give a physical interpretation of this result in terms of forward scattering amplitudes of on-shell thermal particles.

Exercise training improves relaxation response and SOD-1 expression in aortic and mesenteric rings from high caloric diet-fed rats

Abstract Background Obesity has been associated with a variety of disease such as type II diabetes mellitus, arterial hypertension and atherosclerosis. Evidences have shown that exercise training promotes beneficial effects on these disorders, but the underlying mechanisms are not fully understood. The aim of this study was to investigate whether physical preconditioning prevents the deleterious effect of high caloric diet in vascular reactivity of rat aortic and mesenteric rings. Methods Male Wistar rats were divided into sedentary (SD); trained (TR); sedentary diet (SDD) and trained diet (TRD) groups. Run training (RT) was performed in sessions of 60 min, 5 days/week for 12 weeks (70–80% VO2max). Triglycerides, glucose, insulin and nitrite/nitrate concentrations (NOx-) were measured. Concentration-response curves to acetylcholine (ACh) and sodium nitroprusside (SNP) were obtained. Expression of Cu/Zn superoxide dismutase (SOD-1) was assessed by Western blotting. Results High caloric diet increased triglycerides concentration (SDD: 216 ± 25 mg/dl) and exercise training restored to the baseline value (TRD: 89 ± 9 mg/dl). Physical preconditioning significantly reduced insulin levels in both groups (TR: 0.54 ± 0.1 and TRD: 1.24 ± 0.3 ng/ml) as compared to sedentary animals (SD: 0.87 ± 0.1 and SDD: 2.57 ± 0.3 ng/ml). On the other hand, glucose concentration was slightly increased by high caloric diet, and RT did not modify this parameter (SD: 126 ± 6; TR: 140 ± 8; SDD: 156 ± 8 and TRD 153 ± 9 mg/dl). Neither high caloric diet nor RT modified NOx- levels (SD: 27 ± 4; TR: 28 ± 6; SDD: 27 ± 3 and TRD: 30 ± 2 μM). Functional assays showed that high caloric diet impaired the relaxing response to ACh in mesenteric (about 13%), but not in aortic rings. RT improved the relaxing responses to ACh either in aortic (28%, for TR and 16%, to TRD groups) or mesenteric rings (10%, for TR and 17%, to TRD groups) that was accompanied by up-regulation of SOD-1 expression and reduction in triglycerides levels. Conclusion The improvement in endothelial function by physical preconditioning in mesenteric and aortic arteries from high caloric fed-rats was directly related to an increase in NO bioavailability to the smooth muscle mostly due to SOD-1 up regulation.

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)

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.