Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

Simulation of needle-type corona electrodes by the finite element method

This paper describes a software tool, called LEVSOFT, suitable for the electric field simulations of corona electrodes by the Finite Element Method (FEM). Special attention was paid to the user friendly construction of geometries with corners and sharp points, and to the fast generation of highly refined triangular meshes and field maps. The execution of self-adaptive meshes was also implemented. These customized features make the code attractive for the simulation of needle-type corona electrodes. Some case examples involving needle type electrodes are presented.

Study of the phenomena of fatigue using the finite element method

About 99% of mechanical failures are consequence of the phenomena of fatigue, which consists on the progressive weakening of the resistant section of a mechanical component due to the growing of cracks caused by fluctuating loadings. A broad diversity of factors influences the fatigue life of a mechanical component, like the surface finishing, scale factors, among others, but none is as significantly as the presence of geometric severities. Stress concentrators are places where fatigue cracks have a greater probability to occur, and so on, the intuit of this work is to develop a consistent and trustfully methodology to determine the theoretical stress concentration factor of mechanical components. Copyright © 2007 SAE International.

Influence of voids in the hybrid layer based on self-etching adhesive systems: a 3-D FE analysis

The presence of porosities at the dentin/adhesive interface has been observed with the use of new generation dentin bonding systems. These porosities tend to contradict the concept that etching and hybridization processes occur equally and simultaneously. Therefore, the aim of this study was to evaluate the micromechanical behavior of the hybrid layer (HL) with voids based on a self-etching adhesive system using 3-D finite element (FE) analysis. Three FE models (Mr) were built: Mr, dentin specimen (41x41x82 μm) with a regular and perfect (i.e. pore-free) HL based on a self-etching adhesive system, restored with composite resin; Mp, similar to M, but containing 25% (v/v) voids in the HL; Mpp, similar to Mr, but containing 50% (v/v) voids in the HL. A tensile load (0.03N) was applied on top of the composite resin. The stress field was obtained by using Ansys Workbench 10.0. The nodes of the base of the specimen were constrained in the x, y and z axes. The maximum principal stress (σmax) was obtained for all structures at the dentin/adhesive interface. The Mpp showed the highest peak of σmax in the HL (32.2 MPa), followed by Mp (30 MPa) and Mr (28.4 MPa). The stress concentration in the peritubular dentin was high in all models (120 MPa). All other structures positioned far from voids showed similar increase of stress. Voids incorporated into the HL raised the σmax in this region by 13.5%. This behavior might be responsible for lower bond strengths of self-etching and single-bottle adhesives, as reported in the literature.

Influence of high insertion torque on implant placement - an anisotropic bone stress analysis

The aim of this study was to evaluate the influence of the high values of insertion torques on the stress and strain distribution in cortical and cancellous bones. Based on tomography imaging, a representative mathematical model of a partial maxilla was built using Mimics 11.11 and Solid Works 2010 softwares. Six models were built and each of them received an implant with one of the following insertion torques: 30, 40, 50, 60, 70 or 80 Ncm on the external hexagon. The cortical and cancellous bones were considered anisotropic. The bone/implant interface was considered perfectly bonded. The numerical analysis was carried out using Ansys Workbench 10.0. The convergence of analysis (6%) drove the mesh refinement. Maximum principal stress (σ max) and maximum principal strain (ε max) were obtained for cortical and cancellous bones around to implant. Pearson's correlation test was used to determine the correlation between insertion torque and stress concentration in the periimplant bone tissue, considering the significance level at 5%. The increase in the insertion torque generated an increase in the σ max and ε max values for cortical and cancellous bone. The σmax was smaller for the cancellous bone, with greater stress variation among the insertion torques. The ε max was higher in the cancellous bone in comparison to the cortical bone. According to the methodology used and the limits of this study, it can be concluded that higher insertion torques increased tensile and compressive stress concentrations in the periimplant bone tissue.

Regular and switching platform: Bone stress analysis with varying implant diameter

The aim of this study was to evaluate stress distribution of the peri-implant bone by simulating the biomechanical influence of implants with different diameters of regular or platform switched connections by means of 3-dimensional finite element analysis. Five mathematical models of an implant-supported central incisor were created by varying the diameter (5.5 and 4.5 mm, internal hexagon) and abutment platform (regular and platform switched). For the cortical bone, the highest stress values (rmax and rvm) were observed in situation R1, followed by situations S1, R2, S3, and S2. For the trabecular bone, the highest stress values (rmax) were observed in situation S3, followed by situations R1, S1, R2, and S2. The influence of platform switching was more evident for cortical bone than for trabecular bone and was mainly seen in large platform diameter reduction.

Free tools and strategies for the generation of 3D finite element meshes: Modeling of the cardiac structures

The Finite Element Method is a well-known technique, being extensively applied in different areas. Studies using the Finite Element Method (FEM) are targeted to improve cardiac ablation procedures. For such simulations, the finite element meshes should consider the size and histological features of the target structures. However, it is possible to verify that some methods or tools used to generate meshes of human body structures are still limited, due to nondetailed models, nontrivial preprocessing, or mainly limitation in the use condition. In this paper, alternatives are demonstrated to solid modeling and automatic generation of highly refined tetrahedral meshes, with quality compatible with other studies focused on mesh generation. The innovations presented here are strategies to integrate Open Source Software (OSS). The chosen techniques and strategies are presented and discussed, considering cardiac structures as a first application context. © 2013 E. Pavarino et al.

Use of Stress Analysis Methods to Evaluate the Biomechanics of Oral Rehabilitation With Implants

Because the biomechanical behavior of dental implants is different from that of natural tooth, clinical problems may occur. The mechanism of stress distribution and load transfer to the implant/bone interface is a critical issue affecting the success rate of implants. Therefore, the aim of this study was to conduct a brief literature review of the available stress analysis methods to study implant-supported prosthesis loading and to discuss their contributions in the biomechanical evaluation of oral rehabilitation with implants. Several studies have used experimental, analytical, and computational models by means of finite element models (FEM), photoelasticity, strain gauges and associations of these methods to evaluate the biomechanical behavior of dental implants. The FEM has been used to evaluate new components, configurations, materials, and shapes of implants. The greatest advantage of the photoelastic method is the ability to visualize the stresses in complex structures, such as oral structures, and to observe the stress patterns in the whole model, allowing the researcher to localize and quantify the stress magnitude. Strain gauges can be used to assess in vivo and in vitro stress in prostheses, implants, and teeth. Some authors use the strain gauge technique with photoelasticity or FEM techniques. These methodologies can be widely applied in dentistry, mainly in the research field. Therefore, they can guide further research and clinical studies by predicting some disadvantages and streamlining clinical time.

Potential Application in Energy Harvesting of Intermodal Energy Exchange in a Frame: FEM Analysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Scientific use of the finite element method in Orthodontics

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The reduced modulus in the analysis of sharp instrumented indentation tests

In the analysis of instrumented indentation data, it is common practice to incorporate the combined moduli of the indenter (E-i) and the specimen (E) in the so-called reduced modulus (E-r) to account for indenter deformation. Although indenter systems with rigid or elastic tips are considered as equivalent if E-r is the same, the validity of this practice has been questioned over the years. The present work uses systematic finite element simulations to examine the role of the elastic deformation of the indenter tip in instrumented indentation measurements and the validity of the concept of the reduced modulus in conical and pyramidal (Berkovich) indentations. It is found that the apical angle increases as a result of the indenter deformation, which influences in the analysis of the results. Based upon the inaccuracies introduced by the reduced modulus approximation in the analysis of the unloading segment of instrumented indentation applied load (P)-penetration depth (delta) curves, a detailed examination is then conducted on the role of indenter deformation upon the dimensionless functions describing the loading stages of such curves. Consequences of the present results in the extraction of the uniaxial stress-strain characteristics of the indented material through such dimensional analyses are finally illustrated. It is found that large overestimations in the assessment of the strain hardening behavior result by neglecting tip compliance. Guidelines are given in the paper to reduce such overestimations.

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.

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.

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)...

Adaptive finite element simulation of stack pollutant emissions over complex terrains

[EN]A three-dimensional finite element model for the pollutant dispersion is presented. In these environmental processes over a complex terrain, a mesh generator capable of adapting itself to the topographic characteristics is essential. The first stage of the model consists on the construction of an adaptive tetrahedral mesh of a rectangular region bounded in its lower part by the terrain and in its upper part by a horizontal plane. Once the mesh is constructed, an adaptive local refinement of tetrahedra is used in order to capture the plume rise. Wind measurements are used to compute an interpolated wind field, that is modified by using a mass-consistent model and perturbing its vertical component to introduce the plume rise effect...