A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Synthesis of YAP phase by a polymeric method and phase progression mechanisms

The phase formation kinetics of YAP (YAlO(3)) synthesized through the polymeric precursor method was investigated by thermal analysis, X-ray diffraction and FT-IR spectroscopy. We demonstrated that the YAP synthesis is highly dependent on the heat and mass transport during all stages of the synthesis route. In the first stages, during the preparation of amorphous precursor, ""hot spots"" need to be suppressed to avoid the occurrence of chemical inhomogeneities. Very high heating rates combined with small amorphous particles are advantageous in the last stage during the formation of crystalline phase. We were able to synthesize nanosized particles of YAP single phase at temperatures around 1100 A degrees C for future preparation of phosphors or ceramics for optics.

Structured minimal-memory inexact quasi-Newton method and secant preconditioners for augmented Lagrangian optimization

Augmented Lagrangian methods for large-scale optimization usually require efficient algorithms for minimization with box constraints. On the other hand, active-set box-constraint methods employ unconstrained optimization algorithms for minimization inside the faces of the box. Several approaches may be employed for computing internal search directions in the large-scale case. In this paper a minimal-memory quasi-Newton approach with secant preconditioners is proposed, taking into account the structure of Augmented Lagrangians that come from the popular Powell-Hestenes-Rockafellar scheme. A combined algorithm, that uses the quasi-Newton formula or a truncated-Newton procedure, depending on the presence of active constraints in the penalty-Lagrangian function, is also suggested. Numerical experiments using the Cute collection are presented.

Modeling and simulation of the anode in direct ethanol fuels cells

Mathematical modeling has been extensively applied to the study and development of fuel cells. In this work, the objective is to characterize a mechanistic model for the anode of a direct ethanol fuel cell and perform appropriate simulations. The software Comsol Multiphysics (R) (and the Chemical Engineering Module) was used in this work. The software Comsol Multiphysics (R) is an interactive environment for modeling scientific and engineering applications using partial differential equations (PDEs). Based on the finite element method, it provides speed and accuracy for several applications. The mechanistic model developed here can supply details of the physical system, such as the concentration profiles of the components within the anode and the coverage of the adsorbed species on the electrode surface. Also, the anode overpotential-current relationship can be obtained. To validate the anode model presented in this paper, experimental data obtained with a single fuel cell operating with an ethanol solution at the anode were used. (C) 2008 Elsevier B.V. All rights reserved.

Movements and mechanical stresses in gas-filled flat plate solar collectors

Sealed gas filled flat plate solar collectors will have stresses in the material since volume and pressure varies in the gas when the temperature changes. Several geometries were analyzed and it could be seen that it is possible reducing the stresses and improve the safety factor of the weakest point in the construction by using larger area and/or reducing the distance between glass and absorber and/or change width and height relationship so the tubes are getting longer. Further it could be shown that the safety factor won't always get improved with reinforcements. It is so because when an already strong part of the collector gets reinforced it will expose weaker parts for higher stresses. The finite element method was used for finding out the stresses.

Fundações superficiais para linhas de transmissão em argilas moles saturadas

As fundações empregadas para linhas de transmissão de alta tensão em depósitos de argilas moles saturadas são geralmente estacas. No Rio Grande do Sul, projetos alternativos com fundações superficiais foram empregados em forma experimental para torres de suspensão, as quais apresentaram bom desempenho com menor custo. Esta pesquisa visa apontar critérios racionais da mecânica dos solos que permitam projetar, a partir do conhecimento do perfil de resistência não drenada, fundações superficiais em depósitos de argilas moles saturadas com maior confiabilidade. Para a determinaçao do perfil de resistência não drenada foi desenvolvido um cone manual com leituras elétricas, o qual foi testado apresentando bons resultados. A vantagem do uso deste equipamento reside no fato de que pode ser transportado manualmente a locais de difícil acesso, onde são implantadas as torres de alta tensão. Os demais parâmetros de solos do local de estudo, situado no município de Canoas no Rio Grande do Sul, foram obtidos a partir de ensaios de compressão triaxial do tipo adensado isotropicamente não drenado (CIU) com medida de poro pressão e de ensaios de adensamento. Foram analisadas as fundações superficiais das torres de suspensão projetadas, empregando teorias analíticas convencionais, baseadas no conhecimento da resistência não drenada, e o método de elementos finitos. A partir deste estudo, é sugerida uma metodologia para projetos de fundações superficiais de torres de suspensão em argilas moles saturadas.

Análise de estruturas de superfície delgadas axissimétricas via método dos elementos finitos com utilização de elemento retilíneo

The present work deals with the linear analysis of bi-dimensional axisymmetric structures, through development and implementation of a Finite Element Method code. The structures are initially studied alone and afterwards compatibilized into coupled structures, that is, assemblages, including tanks and pressure vessels. Examples are analysed and, in order to prove accuracy, the results were compared with those furnished by the analytical solutions

Análise não-linear de treliças pelo método dos elementos finitos posicional

This work presents the positional nonlinear geometric formulation for trusses using different strain measures. The positional formulation presents an alternative approach for nonlinear problems. This formulation considers nodal positions as variables of the nonlinear system instead of displacements (widely found in literature). The work also describes the arc-length method used for tracing equilibrium paths with snap-through and snap-back. Numerical applications for trusses already established in the literature and comparisons with other studies are provided to prove the accuracy of the proposed formulation

Estudo do efeito de substratos metamateriais em parâmetros de antenas de microfita

Microstrip antennas are subject matter in several research fields due to its numerous advantages. The discovery, at 1999, of a new class of materials called metamaterials - usually composed of metallic elements immersed in a dielectric medium, have attracted the attention of the scientific community, due to its electromagnetic properties, especially the ability to use in planar structures, such as microstrip, without interfering with their traditional geometry. The aim of this paper is to analyze the effects of one and bidimensional metamaterial substrates in microstrip antennas, with different configurations of resonance rings, SRR, in the dielectric layer. Fractal geometry is applied to these rings, in seeking to verify a multiband behavior and to reduce the resonance frequency of the antennas. The results are then given by commercial software Ansoft HFSS, used for precise analysis of the electromagnetic behavior of antennas by Finite Element Method (FEM). To reach it, this essay will first perform a literature study on fractal geometry and its generative process. This paper also presents an analysis of microstrip antennas, with emphasis on addressing different types of substrates as part of its electric and magnetic anisotropic behavior. It s performed too an approach on metamaterials and their unique properties

Método de Otimização Topológica em Estruturas Tridimensionais

The topology optimization problem characterize and determine the optimum distribution of material into the domain. In other words, after the definition of the boundary conditions in a pre-established domain, the problem is how to distribute the material to solve the minimization problem. The objective of this work is to propose a competitive formulation for optimum structural topologies determination in 3D problems and able to provide high-resolution layouts. The procedure combines the Galerkin Finite Elements Method with the optimization method, looking for the best material distribution along the fixed domain of project. The layout topology optimization method is based on the material approach, proposed by Bendsoe & Kikuchi (1988), and considers a homogenized constitutive equation that depends only on the relative density of the material. The finite element used for the approach is a four nodes tetrahedron with a selective integration scheme, which interpolate not only the components of the displacement field but also the relative density field. The proposed procedure consists in the solution of a sequence of layout optimization problems applied to compliance minimization problems and mass minimization problems under local stress constraint. The microstructure used in this procedure was the SIMP (Solid Isotropic Material with Penalty). The approach reduces considerably the computational cost, showing to be efficient and robust. The results provided a well defined structural layout, with a sharpness distribution of the material and a boundary condition definition. The layout quality was proporcional to the medium size of the element and a considerable reduction of the project variables was observed due to the tetrahedrycal element

CENSUS HISTOGRAMS: A SIMPLE FEATURE EXTRACTION and MATCHING APPROACH FOR FACE RECOGNITION

Most face recognition approaches require a prior training where a given distribution of faces is assumed to further predict the identity of test faces. Such an approach may experience difficulty in identifying faces belonging to distributions different from the one provided during the training. A face recognition technique that performs well regardless of training is, therefore, interesting to consider as a basis of more sophisticated methods. In this work, the Census Transform is applied to describe the faces. Based on a scanning window which extracts local histograms of Census Features, we present a method that directly matches face samples. With this simple technique, 97.2% of the faces in the FERET fa/fb test were correctly recognized. Despite being an easy test set, we have found no other approaches in literature regarding straight comparisons of faces with such a performance. Also, a window for further improvement is presented. Among other techniques, we demonstrate how the use of SVMs over the Census Histogram representation can increase the recognition performance.

2D evaluation of crack openings using smeared and embedded crack models

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

NUMERICAL SIMULATION of THE BEHAVIOR of THE BEHAVIOR of REINFORCED CONCRETE BEAMS and COLUMNS BY FEM CAASTEM 2000 PROGRAM

The objective of this paper is the numerical study of the behavior of reinforced concrete beams and columns by non-linear numerical simulations. The numerical analysis is based on the finite element method implemented in CASTEM 2000. This program uses the constitutive elastoplastic perfect model for the steel, the Drucker-Prager model for the concrete and the Newton-Raphson for the solution of non-linear systems. This work concentrates on the determination of equilibrium curves to the beams and force-strain curves to the columns. The numeric responses are confronted with experimental results found in the literature in order to check there liability of the numerical analyses.

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

Passivity Versus Unilateral Angular Misfit: Evaluation of Stress Distribution on Implant-Supported Single Crowns: Three-Dimensional Finite Element Analysis

The aim of this study was to evaluate the effect of unilateral angular misfit of 100 Km on stress distribution of implant-supported single crowns with ceramic veneering and gold framework by three-dimensional finite element analysis. Two three-dimensional models representing a maxillary section of premolar region were constructed: group 1 (control)-crown completely adapted to the implant and group 2-crown with unilateral angular misfit of 100 Km. A vertical force of 100 N was applied on 2 centric points of the crown. The von Mises stress was used as an analysis criterion. The stress values and distribution in the main maps (204.4 MPa for group 1 and 205.0 MPa for group 2) and in the other structures (aesthetic veneering, framework, retention screw, implant, and bone tissue) were similar for both groups. The highest stress values were observed between the first and second threads of the retention screw. Considering the bone tissue, the highest stress values were exhibited in the peri-implant cortical bone. The unilateral angular misfit of 100 Km did not influence the stress distribution on the implant-supported prosthesis under static loading.