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.

Noether-type discrete conserved quantities arising from a finite element approximation of a variational problem

In this work, we prove a weak Noether-type Theorem for a class of variational problems that admit broken extremals. We use this result to prove discrete Noether-type conservation laws for a conforming finite element discretisation of a model elliptic problem. In addition, we study how well the finite element scheme satisfies the continuous conservation laws arising from the application of Noether’s first theorem (1918). We summarise extensive numerical tests, illustrating the conservation of the discrete Noether law using the p-Laplacian as an example and derive a geometric-based adaptive algorithm where an appropriate Noether quantity is the goal functional.

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.

Numeric simulation of pollutant dispersion by a control-volume based on finite element method

The dispersion of pollutants in the environment is an issue of great interest as it directly affects air quality, mainly in large cities. Experimental and numerical tools have been used to predict the behavior of pollutant species dispersion in the atmosphere. A software has been developed based on the control-volume based on the finite element method in order to obtain two-dimensional simulations of Navier-Stokes equations and heat or mass transportation in regions with obstacles, varying position of the pollutant source. Numeric results of some applications were obtained and, whenever possible, compared with literature results showing satisfactory accordance. Copyright (C) 2010 John Wiley & Sons, Ltd.

Stress Analysis in Simulation Models With or Without Implant Threads Representation

Purpose: This study aimed to evaluate the influence of implants with or without threads representation on the outcome of a two-dimensional finite element (FE) analysis. Materials and Methods: Two-dimensional FE models that reproduced a frontal section of edentulous mandibular posterior bone were constructed using a standard crown/implant/screw system representation. To evaluate the effect of implant threads, two models were created: a model in which the implant threads were accurately simulated (precise model) and a model in which implants with a smooth surface (press-fit implant) were used (simplified model). An evaluation was performed on ANSYS software, in which a load of 133 N was applied at a 30-degree angulation and 2 mm off-axis from the long axis of the implant on the models, The Von Mises stresses were measured. Results: The precise model (1.45 MPa) showed higher maximum stress values than the simplified model (1.2 MPa). Whereas in the cortical bone, the stress values differed by about 36% (292.95 MPa for the precise model and 401.14 MPa for the simplified model), in trabecular bone (19.35 MPa and 20.35 MPa, respectively), the stress distribution and stress values were similar. Stress concentrations occurred around the implant neck and the implant apex. Conclusions: Considering implant and cortical bone analysis, remarkable differences in stress values were found between the models. Although the models showed different absolute stress values, the stress distribution was similar. INT J ORAL MAXILLOFAC IMPLANTS 2009;24:1040-1044

Etch and rinse versus self-etching adhesives systems: Tridimensional micromechanical analysis of dentin/adhesive interface

The purpose of this study was to evaluate stress distribution in the hybrid layer produced by two adhesive systems using three-dimensional finite element analysis (FEA). Four FEA models (M) were developed: Mc, a representation of a dentin specimen (41 x 41 x 82 mu m) restored with composite resin, exhibiting the adhesive layer, hybrid layer (HL), resin tags, peritubular dentin, and intertubular dentin to simulate the etch-and-rinse adhesive system; Mr, similar to Mc, with lateral branches of the adhesive; Ma, similar to Mc, however without resin tags and obliterated tubule orifice, to simulate the environment for the self-etching adhesive system; Mat, similar to Ma, with tags. A numerical simulation was performed to obtain the maximum principal stress (sigma(max)). The highest sigma(max) in the HL was observed for the etch-and-rinse adhesive system. The lateral branches increased the sigma(max) in the HL. The resin tags had a little influence on stress distribution with the self-etching system. (C) 2012 Elsevier Ltd. All rights reserved.

Short Implant to Support Maxillary Restorations: Bone Stress Analysis Using Regular and Switching Platform

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

Regular and Platform Switching: Bone Stress Analysis Varying Implant Type

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

Methods Used for Assessing Stresses in Buccomaxillary Prostheses: Photoelasticity, Finite Element Technique, and Extensometry

The authors describe a literature revision on assessing stresses in buccomaxillary prostheses photoelasticity, finite element technique, and extensometry. They describe the techniques and the importance for use of each method in buccomaxillary prostheses with implants and the need of accomplishing more studies in this scarce literary area.

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.

A finite element free convection model for the side wall heated cavity

This paper presents a finite element numerical solution of free convection in a cavity with side walls maintained at constant but different temperatures. The predictions from the model and the method of solution were validated by comparison with the 'bench mark' solution and Vahl Davis' results and good agreement was found. The present model was used to obtain additional results over a wide range of Rayleigh number (10(3)-10(6)) and L/H ratios varying from 0.1 to 1.0. The predicted stream function patterns, temperature and velocity profiles as well as the mean Nusselt number were presented and discussed. (C) 2000 Elsevier B.V. Ltd. All rights reserved.

Non-linear boundary element analysis of floor slabs reinforced with rectangular beams

In this work, a numerical model to perform non-linear analysis of building floor structures is proposed. The presented model is derived from the Kirchhoff-s plate bending formulation of the boundary element method (BENI) for zoned domains, in which the plate stiffness is modified by the presence of membrane effects. In this model, no approximation of the generalized forces along the interface is required and the compatibility and equilibrium conditions along interfaces are imposed at the integral equation level. In order to reduce the number of degrees of freedom, the Navier Bernoulli hypothesis is assumed to simplify the strain field for the thin sub-regions (rectangular beams). The non-linear formulation is obtained from the linear formulation by incorporating initial internal force fields, which are approximated by using the well-known cell sub-division. Then, the non-linear solution of algebraic equations is obtained by using the concept of the consistent tangent operator. The Von Mises criterion is adopted to govern the elasto-plastic material behaviour checked at points along the plate thickness and along the rectangular beam element axes. The numerical representations are accurately obtained by either computing analytically the element integrals or performing the numerical integration accurately using an appropriate sub-elementation scheme. (C) 2007 Elsevier Ltd. All rights reserved.