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.

Análise numérica do efeito de fatores influentes da reação álcali-agregado no desempenho de estruturas de concreto

The aim of this work is the numerical simulation of the mechanical performance of concrete affected by Alkali-Aggregate Reaction or RAA, reported by Stanton in 1940. The RAA has aroused attention in the context of Civil Engineering from the early 80, when they were reported consequences of his swelling effect in concrete structures, including cracking, failure and loss of serviceability. Despite the availability of experimental results the problem formulation still lacks refinement so that your solution remains doubtful. The numerical simulation is important resource for the assessment of damages in structures caused by the reaction, and their recoveries The tasks of support of this work were performed by means of the finite element approach, about orthotropic non-linear formulation, and, thermodynamic model of deformation by RAA. The results obtained revealed that the swelling effect of RAA induced decline of the mechanical performance of concrete by decreasing the margin of safety prior to the material failure. They showed that the temperature influences, exclusively, the kinetics of the reaction, so that the failure was the more precocious the higher the temperature of the solid mass of concrete

Viabilidade da Utilização de Máquinas de Indução Convencionais como Motores sem Mancais Mecânicos

Electrical Motors transform electrical energy into mechanic energy in a relatively easy way. In some specific applications, there is a need for electrical motors to function with noncontaminated fluids, in high speed systems, under inhospitable conditions, or yet, in local of difficult access and considerable depth. In these cases, the motors with mechanical bearings are not adequate as their wear give rise to maintenance. A possible solution for these problems stems from two different alternatives: motors with magnetic bearings, that increase the length of the machine (not convenient), and the bearingless motors that aggregate compactness. Induction motors have been used more and more in research, as they confer more robustness to bearingless motors compared to other types of machines building with others motors. The research that has already been carried out with bearingless induction motors utilized prototypes that had their structures of stator/rotor modified, that differ most of the times from the conventional induction motors. The goal of this work is to study the viability of the use of conventional induction Motors for the beringless motors applications, pointing out the types of Motors of this category that can be more useful. The study uses the Finite Elements Method (FEM). As a means of validation, a conventional induction motor with squirrel-cage rotor was successfully used for the beringless motor application of the divided winding type, confirming the proposed thesis. The controlling system was implemented in a Digital Signal Processor (DSP)

Otimização de forma aplicando B-splines sob critério integral de tensões

This work proposes a computational methodology to solve problems of optimization in structural design. The application develops, implements and integrates methods for structural analysis, geometric modeling, design sensitivity analysis and optimization. So, the optimum design problem is particularized for plane stress case, with the objective to minimize the structural mass subject to a stress criterion. Notice that, these constraints must be evaluated at a series of discrete points, whose distribution should be dense enough in order to minimize the chance of any significant constraint violation between specified points. Therefore, the local stress constraints are transformed into a global stress measure reducing the computational cost in deriving the optimal shape design. The problem is approximated by Finite Element Method using Lagrangian triangular elements with six nodes, and use a automatic mesh generation with a mesh quality criterion of geometric element. The geometric modeling, i.e., the contour is defined by parametric curves of type B-splines, these curves hold suitable characteristics to implement the Shape Optimization Method, that uses the key points like design variables to determine the solution of minimum problem. A reliable tool for design sensitivity analysis is a prerequisite for performing interactive structural design, synthesis and optimization. General expressions for design sensitivity analysis are derived with respect to key points of B-splines. The method of design sensitivity analysis used is the adjoin approach and the analytical method. The formulation of the optimization problem applies the Augmented Lagrangian Method, which convert an optimization problem constrained problem in an unconstrained. The solution of the Augmented Lagrangian function is achieved by determining the analysis of sensitivity. Therefore, the optimization problem reduces to the solution of a sequence of problems with lateral limits constraints, which is solved by the Memoryless Quasi-Newton Method It is demonstrated by several examples that this new approach of analytical design sensitivity analysis of integrated shape design optimization with a global stress criterion purpose is computationally efficient

Otimização topológica de estruturas termoelásticas tridimensionais

This work presents an optimization technique based on structural topology optimization methods, TOM, designed to solve problems of thermoelasticity 3D. The presented approach is based on the adjoint method of sensitivity analysis unified design and is intended to loosely coupled thermomechanical problems. The technique makes use of analytical expressions of sensitivities, enabling a reduction in the computational cost through the use of a coupled field adjoint equation, defined in terms the of temperature and displacement fields. The TOM used is based on the material aproach. Thus, to make the domain is composed of a continuous distribution of material, enabling the use of classical models in nonlinear programming optimization problem, the microstructure is considered as a porous medium and its constitutive equation is a function only of the homogenized relative density of the material. In this approach, the actual properties of materials with intermediate densities are penalized based on an artificial microstructure model based on the SIMP (Solid Isotropic Material with Penalty). To circumvent problems chessboard and reduce dependence on layout in relation to the final optimal initial mesh, caused by problems of numerical instability, restrictions on components of the gradient of relative densities were applied. The optimization problem is solved by applying the augmented Lagrangian method, the solution being obtained by applying the finite element method of Galerkin, the process of approximation using the finite element Tetra4. This element has the ability to interpolate both the relative density and the displacement components and temperature. As for the definition of the problem, the heat load is assumed in steady state, i.e., the effects of conduction and convection of heat does not vary with time. The mechanical load is assumed static and distributed

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

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.

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

Modelagem Matemática e Computacional de fenômenos eletrocinéticos em meios porosos carregados eletricamente

In this work we present a mathematical and computational modeling of electrokinetic phenomena in electrically charged porous medium. We consider the porous medium composed of three different scales (nanoscopic, microscopic and macroscopic). On the microscopic scale the domain is composed by a porous matrix and a solid phase. The pores are filled with an aqueous phase consisting of ionic solutes fully diluted, and the solid matrix consists of electrically charged particles. Initially we present the mathematical model that governs the electrical double layer in order to quantify the electric potential, electric charge density, ion adsorption and chemical adsorption in nanoscopic scale. Then, we derive the microscopic model, where the adsorption of ions due to the electric double layer and the reactions of protonation/ deprotanaç~ao and zeta potential obtained in modeling nanoscopic arise in microscopic scale through interface conditions in the problem of Stokes and Nerst-Planck equations respectively governing the movement of the aqueous solution and transport of ions. We developed the process of upscaling the problem nano/microscopic using the homogenization technique of periodic structures by deducing the macroscopic model with their respectives cell problems for effective parameters of the macroscopic equations. Considering a clayey porous medium consisting of kaolinite clay plates distributed parallel, we rewrite the macroscopic model in a one-dimensional version. Finally, using a sequential algorithm, we discretize the macroscopic model via the finite element method, along with the interactive method of Picard for the nonlinear terms. Numerical simulations on transient regime with variable pH in one-dimensional case are obtained, aiming computational modeling of the electroremediation process of clay soils contaminated

A BEM formulation based on Reissner's theory to perform simple bending analysis of plates reinforced by rectangular beams

In this work, the plate bending formulation of the boundary element method (BEM), based on the Reissner's hypothesis, is extended to the analysis of plates reinforced by rectangular beams. This composed structure is modelled by a zoned plate, being the beams represented by narrow sub-regions with larger thickness. The integral equations are derived by applying the weighted residual method to each sub-region, and summing them to get the equation for the whole plate. Equilibrium and compatibility conditions are automatically imposed by the integral equations, which treat this composed structure as a single body. In order to decrease the number of degrees of freedom, some approximations are considered for both displacements and tractions along the beam width. The accuracy of the proposed model is illustrated by simple examples whose exact solution are known as well as by more complex examples whose numerical results are compared with a well-known finite element code.

Discrete generator coordinate: Symmetry restoration in finite-dimensional spaces

Group theoretical-based techniques and fundamental results from number theory are used in order to allow for the construction of exact projectors in finite-dimensional spaces. These operators are shown to make use only of discrete variables, which play the role of discrete generator coordinates, and their application in the number symmetry restoration is carried out in a nuclear BCS wave function which explicitly violates that symmetry. © 1999 Published by Elsevier Science B.V. All rights reserved.

A general technique to embed non-uniform discontinuities into standard solid finite elements

A general technique to embed non-uniform displacement discontinuities into standard solid finite elements is presented. The technique is based on the decomposition of the kinematic fields into a component related to the deformation of the solid portion of the element and one related to the rigid-body motion due to a displacement discontinuity. This decomposition simplifies the incorporation of discontinuity interfaces and provides a suitable framework to account for non-uniform discontinuity modes. The present publication addresses two families of finite element formulations: displacement-based and stress hybrid finite element. © 2005 Elsevier Ltd. All rights reserved.

Blade design and forming for fans using finite elements

The necessity of adapting the standardized fan models to conditions of higher temperature has emerged due to the growth of concerning referring to the consequences of the gas expelling after the Mont Blanc tunnel accident in Italy and France, where even though, with 100 fans in operation, 41 people died. However, since then, the defied solutions have pointed to aerodynamic disadvantages or have seemed nonappropriate in these conditions. The objective of this work is to present an alternative to the market standard fans considering a new technology in constructing blades. This new technology introduces the use of the stainless steel AISI 409 due to its good adaptation to temperatures higher than 400°C, particularly exposed to temperatures of gas exhaust from tunnels in fire situation. Furthermore, it presents a very good resistance to corrosion and posterior welding and pressing, due to its alloyed elements. The innovation is centered in the process of a deep drawing of metallic shells and posterior welding, in order to keep the ideal aerodynamic superficies for the fan ideal performance. On the other hand, the finite element method, through the elasto-plastic software COSMOS permitted the verification of the thickness and structural stability of the blade in relation to the aerodynamic efforts established in the project. In addition, it is not advisable the fabrication of blades with variable localized thickness not even, non-uniform ones, due to the verified concentration of tensions and the difficulties observed in the forming. In this way, this study recommends the construction of blades with uniform variations of thickness. © 2007 Springer.

Analysis of mechanical behavior in bending of polymer composites using the finite elements method

The use of composite materials has increased in the recent decades, mainly in the aeronautics and automotives industries. In the present study is elaborated a computational simulation program of the bending test using the finite elements method, in the commercial software ANSYS. This simulation has the objective of analyze the mechanical behavior in bending of two composites with polymeric matrix reinforced with carbon fibers. Also are realized bending tests of the 3 points to obtain the resistances of the materials. Data from simulation and tests are used to make a comparison between two failures criteria, Tsai-Wu and Hashin criterion. Copyright © 2009 SAE International.

Strain rate of the South American lithospheric plate by SIRGAS-CON geodetic observations

A multiyear solution of the SIRGAS-CON network was used to estimate the strain rates of the earth surface from the changing directions of the velocity vectors of 140 geodetic points located in the South American plate. The strain rate was determined by the finite element method using Delaunay triangulation points that formed sub-networks; each sub-network was considered a solid and homogeneous body. The results showed that strain rates vary along the South American plate and are more significant on the western portion of the plate, as expected, since this region is close to the subduction zone of the Nazca plate beneath the South American plate. After using Euler vectors to infer Nazca plate movement and to orient the velocity vectors of the South American plate, it was possible to estimate the convergence and accommodation rates of the Nazca and South American plates, respectively. Strain rate estimates permitted determination of predominant contraction and/or extension regions and to establish that contraction regions coincide with locations with most of the high magnitude seismic events. Some areas with extension and contraction strains were found to the east within the stable South American plate, which may result from different stresses associated with different geological characteristics. These results suggest that major movements detected on the surface near the Nazca plate occur in regions with more heterogeneous geological structures and multiple rupture events. Most seismic events in the South American plate are concentrated in areas with predominant contraction strain rates oriented northeast-southwest; significant amounts of elastic strain can be accumulated on geological structures away from the plate boundary faults; and, behavior of contractions and extensions is similar to what has been found in seismological studies. © 2013 Elsevier Ltd.