A moving-mesh finite element method and its application to the numerical solution of phase-change problems

A distributed Lagrangian moving-mesh finite element method is applied to problems involving changes of phase. The algorithm uses a distributed conservation principle to determine nodal mesh velocities, which are then used to move the nodes. The nodal values are obtained from an ALE (Arbitrary Lagrangian-Eulerian) equation, which represents a generalization of the original algorithm presented in Applied Numerical Mathematics, 54:450--469 (2005). Having described the details of the generalized algorithm it is validated on two test cases from the original paper and is then applied to one-phase and, for the first time, two-phase Stefan problems in one and two space dimensions, paying particular attention to the implementation of the interface boundary conditions. Results are presented to demonstrate the accuracy and the effectiveness of the method, including comparisons against analytical solutions where available.

A finite element-spherical harmonics model for radiative transfer in inhomogeneous clouds. Part II. some applications

EVENT has been used to examine the effects of 3D cloud structure, distribution, and inhomogeneity on the scattering of visible solar radiation and the resulting 3D radiation field. Large eddy simulation and aircraft measurements are used to create realistic cloud fields which are continuous or broken with smooth or uneven tops. The values, patterns and variance in the resulting downwelling and upwelling radiation from incident visible solar radiation at different angles are then examined and compared to measurements. The results from EVENT confirm that 3D cloud structure is important in determining the visible radiation field, and that these results are strongly influenced by the solar zenith angle. The results match those from other models using visible solar radiation, and are supported by aircraft measurements of visible radiation, providing confidence in the new model.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

Stress Analysis in Platform-Switching Implants: A 3-Dimensional Finite Element Study

The aim of this study was to evaluate the influence of the platform-switching technique on stress distribution in implant, abutment, and pen-implant tissues, through a 3-dimensional finite element study. Three 3-dimensional mandibular models were fabricated using the Solid Works 2006 and InVesalius software. Each model was composed of a bone block with one implant 10 mm long and of different diameters (3.75 and 5.00 mm). The UCLA abutments also ranged in diameter from 5.00 mm to 4.1 mm. After obtaining the geometries, the models were transferred to the software FEMAP 10.0 for pre- and postprocessing of finite elements to generate the mesh, loading, and boundary conditions. A total load of 200 N was applied in axial (0 degrees), oblique (45 degrees), and lateral (90) directions. The models were solved by the software NeiNastran 9.0 and transferred to the software FEMAP 10.0 to obtain the results that were visualized through von Mises and maximum principal stress maps. Model A (implants with 3.75 mm/abutment with 4.1 mm) exhibited the highest area of stress concentration with all loadings (axial, oblique, and lateral) for the implant and the abutment. All models presented the stress areas at the abutment level and at the implant/abutment interface. Models B (implant with 5.0 mm/abutment with 5.0 mm) and C (implant with 5.0 mm/abutment with 4.1 mm) presented minor areas of stress concentration and similar distribution pattern. For the cortical bone, low stress concentration was observed in the pen-implant region for models B and C in comparison to model A. The trabecular bone exhibited low stress that was well distributed in models B and C. Model A presented the highest stress concentration. Model B exhibited better stress distribution. There was no significant difference between the large-diameter implants (models B and C).

Effect of Passive Fit Absence in the Prosthesis/Implant/Retaining Screw System: A Two-Dimensional Finite Element Analysis

The misfit between prostheses and implants is a clinical reality, but the level that can be accepted without causing mechanical or biologic problem is not well defined. This study investigates the effect of different levels of unilateral angular misfit prostheses in the prosthesis/implant/retaining screw system and in the surrounding bone using finite element analysis. Four models of a two-dimensional finite element were constructed: group I (control), prosthesis that fit the implant; groups 2 to 4, prostheses with unilateral angular misfit of 50, 100, and 200 mu m, respectively. A load of 133 N was applied with a 30-degree angulation and off-axis at 2 mm from the long axis of the implant at the opposite direction of misfit on the models. Taking into account the increase of the angular misfit, the stress maps showed a gradual increase of prosthesis stress and uniform stress in the implant and trabecular bone. Concerning the displacement, an inclination of the system due to loading and misfit was observed. The decrease of the unilateral contact between prosthesis and implant leads to the displacement of the entire system, and distribution and magnitude alterations of the stress also occurred.

Effect of Metal-Ceramic or All-Ceramic Superstructure Materials on Stress Distribution in a Single Implant-Supported Prosthesis: Three-Dimensional Finite Element Analysis

Purpose: This three-dimensional finite element analysis study evaluated the effect of different material combinations on stress distribution within metal-ceramic and all-ceramic single implant-supported prostheses. Materials and Methods: Three-dimensional finite element models reproducing a segment of the maxilla with a missing left first premolar were created. Five groups were established to represent different superstructure materials: GP, porcelain fused to gold alloy; GR, modified composite resin fused to gold alloy; TP, porcelain fused to titanium; TR, modified composite resin fused to titanium; and ZP, porcelain fused to zirconia. A 100-N vertical force was applied to the contact points of the crowns. All models were fixed in the superior region of bone tissue and in the mesial and distal faces of the maxilla section. Stress maps were generated by processing with finite element software. Results: Stress distribution and stress values of supporting bone were similar for the GP, GR, TP, and ZP models (1,574.3 MPa, 1,574.3 MPa, 1,574.3 MPa, and 1,574.2 MPa, respectively) and different for the TR model (1,838.3 MPa). The ZP model transferred less stress to the retention screw (785 MPa) than the other groups (939 MPa for GP, 961 MPa for GR, 1,010 MPa for TP, and 1,037 MPa for TR). Conclusion: The use of different materials to fabricate a superstructure for a single implant-supported prosthesis did not affect the stress distribution in the supporting bone. The retention screw received less stress when a combination of porcelain and zirconia was used. Int J Oral Maxillofac Implants 2011;26:1202-1209

Cortical Bone Stress Distribution in Mandibles with Different Configurations Restored with Prefabricated Bar-Prosthesis Protocol: A Three-Dimensional Finite-Element Analysis

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

Influence of Implant Design on the Biomechanical Environment of Immediately Placed Implants: Computed Tomography-Based Nonlinear Three-Dimensional Finite Element Analysis

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

Biomechanical Evaluation of Platform Switching in Different Implant Protocols: Computed Tomography-Based Three-Dimensional Finite Element Analysis

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

MATHEMATICA notebook for computing tetrahedral finite element shape functions and matrices for the Helmholtz equation

A MATHEMATICA notebook to compute the elements of the matrices which arise in the solution of the Helmholtz equation by the finite element method (nodal approximation) for tetrahedral elements of any approximation order is presented. The results of the notebook enable a fast computational implementation of finite element codes for high order simplex 3D elements reducing the overheads due to implementation and test of the complex mathematical expressions obtained from the analytical integrations. These matrices can be used in a large number of applications related to physical phenomena described by the Poisson, Laplace and Schrodinger equations with anisotropic physical properties.

Stress distribution in ceramic restorations over natural tooth using finite element analysis. lithium disilicate x aluminum oxide material

Background: Data on stress distribution in tooth-restoration interface with different ceramic restorative materials are limited. The aim of this chapter was to assess the stress distribution in the interface of ceramic restorations with laminate veneer or full-coverage crown with two different materials (lithium dissilicate and densely sintered aluminum oxide) under different loading areas through finite element analysis. Materials and Methods: Six two-dimensional finite element models were fabricated with different restorations on natural tooth: laminate veneer (IPS Empress, IPS Empress Esthetic and Procera AllCeram) or full-coverage crown (IPS e.max Press and Procera AllCeram). Two different loading areas (L) (50N) were also determined: palatal surface at 45° in relation to the long axis of tooth (L1) and perpendicular to the incisal edge (L2). A model with higid natural tooth was used as control. von Mises equivalent stress (σ vM) and maximum principal stress (σ max) were obtained on Ansys software. Results: The presence of ceramic restoration increased σ vM and σ max in the adhesive interface, mainly for the aluminum oxide (Procera AllCeram system) restorations. The full-coverage crowns generated higher stress in the adhesive interface under L1 while the same result was observed for the laminate veneers under L2. Conclusions: Lithium dissilicate and densely sintered aluminum oxide restorations exhibit different behavior due to different mechanical properties and loading conditions. © 2011 Nova Science Publishers, Inc.

Comparison of different designs of implant-retained overdentures and fixed full-arch implant-supported prosthesis on stress distribution in edentulous mandible - A computed tomography-based three-dimensional finite element analysis

A finite element analysis was used to compare the effect of different designs of implant-retained overdentures and fixed full-arch implant-supported prosthesis on stress distribution in edentulous mandible. Four models of an human mandible were constructed. In the OR (O'ring) group, the mandible was restored with an overdenture retained by four unsplinted implants with O'ring attachment; in the BC (bar-clip) -C and BC groups, the mandibles were restored with overdentures retained by four splinted implants with bar-clip anchor associated or not with two distally placed cantilevers, respectively; in the FD (fixed denture) group, the mandible was restored with a fixed full-arch four-implant-supported prosthesis. Models were supported by the masticatory muscles and temporomandibular joints. A 100-N oblique load was applied on the left first molar. Von Mises (σvM), maximum (σmax) and minimum (σmin) principal stresses (in MPa) analyses were obtained. BC-C group exhibited the highest stress values (σvM=398.8, σmax=580.5 and σmin=-455.2) while FD group showed the lowest one (σvM=128.9, σmax=185.9 and σmin=-172.1). Within overdenture groups, the use of unsplinted implants reduced the stress level in the implant/prosthetic components (59.4% for σvM, 66.2% for σmax and 57.7% for σmin versus BC-C group) and supporting tissues (maximum stress reduction of 72% and 79.5% for σmax, and 15.7% and 85.7% for σmin on the cortical and trabecular bones, respectively). Cortical bone exhibited greater stress concentration than the trabecular bone for all groups. The use of fixed implant dentures and removable dentures retained by unsplinted implants to rehabilitate edentulous mandible reduced the stresses in the periimplant bone tissue, mucosa and implant/prosthetic components. © 2013 Elsevier Ltd.

Avaliação da verticalização de segundo molar inferior ancorado em minimplante através da metodologia dos elementos finitos 3D

Pós-graduação em Odontologia - FOA

Influence of Platform and Abutment Angulation on Peri-Implant Bone. A Three-Dimensional Finite Element Stress Analysis

The aim of this study was to evaluate stress distribution on the pen-implant bone, simulating the influence of Nobel Select implants with straight or angulated abutments on regular and switching platform in the anterior maxilla, by means of 3-dimensional finite element analysis. Four mathematical models of a central incisor supported by external hexagon implant (13 mm x 5 mm) were created varying the platform (R, regular or S. switching) and the abutments (S, straight or A, angulated 15 degrees). The models were created by using Mimics 13 and Solid Works 2010 software programs. The numerical analysis was performed using ANSYS Workbench 10.0. Oblique forces (100 N) were applied to the palatine surface of the central incisor. The bone/implant interface was considered perfectly integrated. Maximum (sigma(max)) and minimum (sigma(min)) principal stress values were obtained. For the cortical bone the highest stress values (sigma(max)) were observed in the RA (regular platform and angulated abutment, 51 MPa), followed by SA (platform switching and angulated abutment, 44.8 MPa), RS (regular platform and straight abutment, 38.6 MPa) and SS (platform switching and straight abutment, 36.5 MPa). For the trabecular bone, the highest stress values (sigma(max)) were observed in the RA (6.55 MPa), followed by RS (5.88 MPa), SA (5.60 MPa), and SS (4.82 MPa). The regular platform generated higher stress in the cervical periimplant region on the cortical and trabecular bone than the platform switching, irrespective of the abutment used (straight or angulated).

Scientific use of the finite element method in Orthodontics

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