A finite-volume shallow layer method for the MHD instabilities in an aluminium production cell

An electrolytic cell for Aluminium production contains molten metal and molten electrolyte, which are subject to high dc-currents and magnetic fields. Lorentz forces arising from the cross product of current and magnetic field may amplify natural gravity waves at the interface between the two fluids, leading to short circuits in extreme cases. The external magnetic field and current distribution in the production cell is computed through a detailed finite element analysis at Torino Polytechnic. The results are then used to compute the magnetohydrodynamic and thermal effects in the aluminium/electrolyte bath. Each cell has lateral dimensions of 6m x 2m, whilst the bath depth is only 30cm. the electrically resistive electrolyte path, which is critical in the operation of the cell, has layer depth of only a few centimetres below each carbon anode. Because the shallow dimensions of the liquid layer a finite-volume shallow-layer technique has been used at Greenwich to compute the resulting flow-field and interface perturbations. The information obtained from this method, i.e. depth averaged velocities and aluminium/electrolyte interface position is then embedded in the three-dimensional finite volume code PHYSICA and will be used to compute the heat transfer and phase change in the cell.

Modélisation des machines par circuits couplés pour l’observation des courants de barres au rotor

Ce travail présente une modélisation rapide d’ordre élévé capable de modéliser une configuration rotorique en cage complète ou en grille, de reproduire les courants de barre et tenir compte des harmoniques d’espace. Le modèle utilise une approche combinée d’éléments finis avec les circuits-couplés. En effet, le calcul des inductances est réalisé avec les éléments finis, ce qui confère une précision avancée au modèle. Cette méthode offre un gain important en temps de calcul sur les éléments finis pour des simulations transitoires. Deux outils de simulation sont développés, un dans le domaine du temps pour des résolutions dynamiques et un autre dans le domaine des phaseurs dont une application sur des tests de réponse en fréquence à l’arrêt (SSFR) est également présentée. La méthode de construction du modèle est décrite en détail de même que la procédure de modélisation de la cage du rotor. Le modèle est validé par l’étude de machines synchrones: une machine de laboratoire de 5.4 KVA et un grand alternateur de 109 MVA dont les mesures expérimentales sont comparées aux résultats de simulation du modèle pour des essais tels que des tests à vide, des courts-circuits triphasés, biphasés et un test en charge.

Coupled CFD Shape Optimization for aerodynamic profiles

The present document deals with the optimization of shape of aerodynamic profiles -- The objective is to reduce the drag coefficient on a given profile without penalising the lift coefficient -- A set of control points defining the geometry are passed and parameterized as a B-Spline curve -- These points are modified automatically by means of CFD analysis -- A given shape is defined by an user and a valid volumetric CFD domain is constructed from this planar data and a set of user-defined parameters -- The construction process involves the usage of 2D and 3D meshing algorithms that were coupled into own- code -- The volume of air surrounding the airfoil and mesh quality are also parametrically defined -- Some standard NACA profiles were used by obtaining first its control points in order to test the algorithm -- Navier-Stokes equations were solved for turbulent, steady-state ow of compressible uids using the k-epsilon model and SIMPLE algorithm -- In order to obtain data for the optimization process an utility to extract drag and lift data from the CFD simulation was added -- After a simulation is run drag and lift data are passed to the optimization process -- A gradient-based method using the steepest descent was implemented in order to define the magnitude and direction of the displacement of each control point -- The control points and other parameters defined as the design variables are iteratively modified in order to achieve an optimum -- Preliminary results on conceptual examples show a decrease in drag and a change in geometry that obeys to aerodynamic behavior principles

Análise da fluência em vigas de concreto armado simplesmente apoiadas usando o método dos elementos finitos

Reinforced concrete creep is a phenomenon of great importance. Despite being appointed as the main cause of several pathologies, its effects are yet considered in a simplified way by the structural designers. In addition to studying the phenomenon in reinforced concrete structures and its current account used in the structural analysis, this paper compares creep strains at simply supported reinforced concrete beams in analytical and in experimental forms with the finite element method (FEM) simulation results. The strains and deflections obtained through the analytical form were calculated with the Brazilian code NBR 6118 (2014) recommendations and the simplified method from CEB-FIP 90 and the experimental results were extracted from tests available in the literature. Finite element simulations are performed using ANSYS Workbench software, using its 3D SOLID 186 elements and the structure symmetry. Analyzes of convergence using 2D PLANE 183 elements are held as well. At the end, it is concluded that FEM analyses are quantitative and qualitative efficient for the estimation of this non-linearity and that the method utilized to obtain the creep coefficients values is sufficiently accurate.

Avaliação do potencial passível de conversão em energia elétrica das correntes na plataforma continental sul brasileira

O contínuo crescimento da população mundial aumenta a demanda e a competição por energia, colocando grande esforço sobre as fontes de energia não renováveis existentes. Devido a isso, políticas globais para geração de energias renováveis e menos poluentes estão sendo fortalecidas, além de promoverem o desenvolvimento de novas tecnologias. Várias formas de conversão de energia foram desenvolvidas no decorrer dos anos, com destaque para os conversores de energia das correntes a base de turbinas, que demonstram alta capacidade de conversão energética e já se encontram em funcionamento. O modelo tridimensional TELEMAC3D foi utilizado para a investigação dos processos hidrodinâmicos. Este modelo foi acoplado ao módulo de conversão de energia para as análises nos locais de maior viabilidade e conversão energética na Plataforma Continental do Sul do Brasil. A região de estudo demonstrou possuir duas regiões com alto potencial para a exploração da energias das correntes marinhas, entretanto a região mais viável para a instalação de conversores de corrente é a região norte delimitada entre o Farol da Conceição e o Farol da Solidão, podendo atingir potência média de 10kW=Dia, e alcançando valores integrados de 3:5MW=Ano. Através de uma análise da sazonalidade foram observados, durante a primavera os períodos mais energéticos em ambas as regiões estudadas. As maiores intensidades de conversão de energia foram estimadas com variabilidade temporal de 16 dias, demonstrando alta correlação com eventos associados à passagem de frentes meteorológicas na região. O sítio da região norte, com a presença de barreiras que representam a forma dos conversores, se destaca mantendo boa conversão durante os eventos de ótimo potencial energético. Esta melhora se deve ao efeito de intensificação do campo de correntes associado à presença da estrutura física que otimiza a eficiência do sítio. Não foram observadas diferenças significativas no padrão de variabilidade temporal das simulações estudadas, indicando que a presença das barreiras não induz grandes alterações no padrão temporal da conversão de energia nas escalas temporais analisadas neste trabalho. Os eventos de alta geração de energia foram relacionados a incidência de fortes ventos de quadrante sul e norte, indicando que pelo formato e disposição dos conversores, ventos de sudoeste e norte podem favorecer ótimos eventos de conversão de energia. As simulações dos sítios de conversão demonstraram alta capacidade de geração energética, com quatro eventos de extrema geração de energia. Entretanto, o sítio da região norte demonstrou eficiência superior a 59,39 GWh ao ano, equivalendo a 0.22% do consumo energético do estado do Rio Grande do Sul no ano de 2010.

Hybrid mathematical and informational modeling of beam-to-column connections

The analysis of steel and composite frames has traditionally been carried out by idealizing beam-to-column connections as either rigid or pinned. Although some advanced analysis methods have been proposed to account for semi-rigid connections, the performance of these methods strongly depends on the proper modeling of connection behavior. The primary challenge of modeling beam-to-column connections is their inelastic response and continuously varying stiffness, strength, and ductility. In this dissertation, two distinct approaches—mathematical models and informational models—are proposed to account for the complex hysteretic behavior of beam-to-column connections. The performance of the two approaches is examined and is then followed by a discussion of their merits and deficiencies. To capitalize on the merits of both mathematical and informational representations, a new approach, a hybrid modeling framework, is developed and demonstrated through modeling beam-to-column connections. Component-based modeling is a compromise spanning two extremes in the field of mathematical modeling: simplified global models and finite element models. In the component-based modeling of angle connections, the five critical components of excessive deformation are identified. Constitutive relationships of angles, column panel zones, and contact between angles and column flanges, are derived by using only material and geometric properties and theoretical mechanics considerations. Those of slip and bolt hole ovalization are simplified by empirically-suggested mathematical representation and expert opinions. A mathematical model is then assembled as a macro-element by combining rigid bars and springs that represent the constitutive relationship of components. Lastly, the moment-rotation curves of the mathematical models are compared with those of experimental tests. In the case of a top-and-seat angle connection with double web angles, a pinched hysteretic response is predicted quite well by complete mechanical models, which take advantage of only material and geometric properties. On the other hand, to exhibit the highly pinched behavior of a top-and-seat angle connection without web angles, a mathematical model requires components of slip and bolt hole ovalization, which are more amenable to informational modeling. An alternative method is informational modeling, which constitutes a fundamental shift from mathematical equations to data that contain the required information about underlying mechanics. The information is extracted from observed data and stored in neural networks. Two different training data sets, analytically-generated and experimental data, are tested to examine the performance of informational models. Both informational models show acceptable agreement with the moment-rotation curves of the experiments. Adding a degradation parameter improves the informational models when modeling highly pinched hysteretic behavior. However, informational models cannot represent the contribution of individual components and therefore do not provide an insight into the underlying mechanics of components. In this study, a new hybrid modeling framework is proposed. In the hybrid framework, a conventional mathematical model is complemented by the informational methods. The basic premise of the proposed hybrid methodology is that not all features of system response are amenable to mathematical modeling, hence considering informational alternatives. This may be because (i) the underlying theory is not available or not sufficiently developed, or (ii) the existing theory is too complex and therefore not suitable for modeling within building frame analysis. The role of informational methods is to model aspects that the mathematical model leaves out. Autoprogressive algorithm and self-learning simulation extract the missing aspects from a system response. In a hybrid framework, experimental data is an integral part of modeling, rather than being used strictly for validation processes. The potential of the hybrid methodology is illustrated through modeling complex hysteretic behavior of beam-to-column connections. Mechanics-based components of deformation such as angles, flange-plates, and column panel zone, are idealized to a mathematical model by using a complete mechanical approach. Although the mathematical model represents envelope curves in terms of initial stiffness and yielding strength, it is not capable of capturing the pinching effects. Pinching is caused mainly by separation between angles and column flanges as well as slip between angles/flange-plates and beam flanges. These components of deformation are suitable for informational modeling. Finally, the moment-rotation curves of the hybrid models are validated with those of the experimental tests. The comparison shows that the hybrid models are capable of representing the highly pinched hysteretic behavior of beam-to-column connections. In addition, the developed hybrid model is successfully used to predict the behavior of a newly-designed connection.

Numerical investigations on the strain-adaptive bone remodelling in the periprosthetic femur: Influence of the boundary conditions

Background: There are several numerical investigations on bone remodelling after total hip arthroplasty (THA) on the basis of the finite element analysis (FEA). For such computations certain boundary conditions have to be defined. The authors chose a maximum of three static load situations, usually taken from the gait cycle because this is the most frequent dynamic activity of a patient after THA. Materials and methods: The numerical study presented here investigates whether it is useful to consider only one static load situation of the gait cycle in the FE calculation of the bone remodelling. For this purpose, 5 different loading cases were examined in order to determine their influence on the change in the physiological load distribution within the femur and on the resulting strain-adaptive bone remodelling. First, four different static loading cases at 25%, 45%, 65% and 85% of the gait cycle, respectively, and then the whole gait cycle in a loading regime were examined in order to regard all the different loadings of the cycle in the simulation. Results: The computed evolution of the apparent bone density (ABD) and the calculated mass losses in the periprosthetic femur show that the simulation results are highly dependent on the chosen boundary conditions. Conclusion: These numerical investigations prove that a static load situation is insufficient for representing the whole gait cycle. This causes severe deviations in the FE calculation of the bone remodelling. However, accompanying clinical examinations are necessary to calibrate the bone adaptation law and thus to validate the FE calculations.

Algoritmos de otimização e o problema inverso aplicados à identificação de danos estruturais com o uso do MEF

Dissertação (mestrado)—Universidade de Brasília, Faculdade de Tecnologia, Departamento de Engenharia Civil e Ambiental, 2015.

Validación de la simulación numérica del flujo bifásico hidrodinámico en sistemas de lecho fluido

Este trabajo tiene como objetivo la mejora en la validación de la simulación numérica del flujo bifásico característico del transporte de lecho fluido, mediante la formulación y desarrollo de un modelo numérico combinado Volúmenes Finitos - Elementos Finitos. Para ello se simula numéricamente el flujo de mezcla sólido-gas en una Cámara de Lecho Fluido, bajo implementación en código COMSOL, cuyos resultados son mejores comparativamente a un modelo basado en el método de Elementos Discretos implementado en código abierto MFIX. El problema fundamental de la modelización matemática del fenómeno de lecho fluido es la irregularidad del dominio, el acoplamiento de las variables en espacio y tiempo y, la no linealidad. En esta investigación se reformula apropiadamente las ecuaciones conservativas del fenómeno, tales que permitan obtener un problema variacional equivalente y solucionable numéricamente. Entonces; se define una ecuación de estado en función de la presión hidrodinámica y la fracción volumétrica de sólidos, quedando desacoplado el sistema en tres sub-problemas, garantizando así la existencia de solución del problema general. Una vez aproximados numéricamente ambos modelos, se comparan los resultados de donde se observa que el modelo materia del presente artículo, verifica de forma más eficaz las condiciones de mezcla óptima, reflejada en la calidad del burbujeo y velocidad de mezcla.

Effect of different feed-rate in bone drilling: experimental and numerical study

The behaviour of bone tissue during drilling has been subject of recent studies due to its great importance. Because of thermal nature of the bone drilling, high temperatures and thermal mechanical stresses are developed during drilling that affect the process quality. However, there is still a lack information with regard to the distribution of mechanical and thermal stresses during bone drilling. The present paper describes a sequentially coupled thermal-stress analysis to assess the mechanical and thermal stress distribution during bone drilling. A three-dimensional thermo-mechanical model was developed using the ANSYS/LSDYNA finite element code under different drilling conditions. The model incorporates the dynamic characteristics of drilling process, as well as the thermo-mechanical properties of the involved materials. Experimental tests with polyurethane foam materials were also carried out. It was concluded that the use of higher feed-rates lead to a decrease of normal stresses and strains in the foam materials. The experimental and numerical results were compared and showed good agreement. The proposed numerical model could be used to predict the better drilling parameters and minimize the bone injuries.

Otimização e análise mecânica de pastas geopoliméricas para uso em poços sujeitos à injeção cíclica de vapor

Oil wells subjected to cyclic steam injection present important challenges for the development of well cementing systems, mainly due to tensile stresses caused by thermal gradients during its useful life. Cement sheath failures in wells using conventional high compressive strength systems lead to the use of cement systems that are more flexible and/or ductile, with emphasis on Portland cement systems with latex addition. Recent research efforts have presented geopolymeric systems as alternatives. These cementing systems are based on alkaline activation of amorphous aluminosilicates such as metakaolin or fly ash and display advantageous properties such as high compressive strength, fast setting and thermal stability. Basic geopolymeric formulations can be found in the literature, which meet basic oil industry specifications such as rheology, compressive strength and thickening time. In this work, new geopolymeric formulations were developed, based on metakaolin, potassium silicate, potassium hydroxide, silica fume and mineral fiber, using the state of the art in chemical composition, mixture modeling and additivation to optimize the most relevant properties for oil well cementing. Starting from molar ratios considered ideal in the literature (SiO2/Al2O3 = 3.8 e K2O/Al2O3 = 1.0), a study of dry mixtures was performed,based on the compressive packing model, resulting in an optimal volume of 6% for the added solid material. This material (silica fume and mineral fiber) works both as an additional silica source (in the case of silica fume) and as mechanical reinforcement, especially in the case of mineral fiber, which incremented the tensile strength. The first triaxial mechanical study of this class of materials was performed. For comparison, a mechanical study of conventional latex-based cementing systems was also carried out. Regardless of differences in the failure mode (brittle for geopolymers, ductile for latex-based systems), the superior uniaxial compressive strength (37 MPa for the geopolymeric slurry P5 versus 18 MPa for the conventional slurry P2), similar triaxial behavior (friction angle 21° for P5 and P2) and lower stifness (in the elastic region 5.1 GPa for P5 versus 6.8 GPa for P2) of the geopolymeric systems allowed them to withstand a similar amount of mechanical energy (155 kJ/m3 for P5 versus 208 kJ/m3 for P2), noting that geopolymers work in the elastic regime, without the microcracking present in the case of latex-based systems. Therefore, the geopolymers studied on this work must be designed for application in the elastic region to avoid brittle failure. Finally, the tensile strength of geopolymers is originally poor (1.3 MPa for the geopolymeric slurry P3) due to its brittle structure. However, after additivation with mineral fiber, the tensile strength became equivalent to that of latex-based systems (2.3 MPa for P5 and 2.1 MPa for P2). The technical viability of conventional and proposed formulations was evaluated for the whole well life, including stresses due to cyclic steam injection. This analysis was performed using finite element-based simulation software. It was verified that conventional slurries are viable up to 204ºF (400ºC) and geopolymeric slurries are viable above 500ºF (260ºC)

Estudo da aplicação de molas magnéticas na recolha da energia de fontes renováveis

Dissertação de Mestrado, Engenharia Elétrica e Eletrónica, Especialização em Sistemas de Energia e Controlo, Instituto Superior de Engenharia, Universidade do Algarve, 2016

Estudo comparativo de pás para aerogeradores de grande porte fabricadas em materiais compósitos reforçadas com fibra de carbono ou fibra de vidro

The research and development of wind turbine blades are essential to keep pace with worldwide growth in the renewable energy sector. Although currently blades are typically produced using glass fiber reinforced composite materials, the tendency for larger size blades, particularly for offshore applications, has increased the interest on carbon fiber reinforced composites because of the potential for increased stiffness and weight reduction. In this study a model of blade designed for large generators (5 MW) was studied on a small scale. A numerical simulation was performed to determine the aerodynamic loading using a Computational Fluid Dynamics (CFD) software. Two blades were then designed and manufactured using epoxy matrix composites: one reinforced with glass fibers and the other with carbon fibers. For the structural calculations, maximum stress failure criterion was adopted. The blades were manufactured by Vacuum Assisted Resin Transfer Molding (VARTM), typical for this type of component. A weight comparison of the two blades was performed and the weight of the carbon fiber blade was approximately 45% of the weight of the fiberglass reinforced blade. Static bending tests were carried out on the blades for various percentages of the design load and deflections measurements were compared with the values obtained from finite element simulations. A good agreement was observed between the measured and calculated deflections. In summary, the results of this study confirm that the low density combined with high mechanical properties of carbon fibers are particularly attractive for the production of large size wind turbine blades

Análise e Prevenção de avarias de um veio de manivelas de um motor Pielstick

Na engenharia mecânica há cada vez mais necessidade de utilizar e de prever o comportamento das máquinas térmicas, mais propriamente dos motores de combustão interna, em especial na área da manutenção e da prevenção de falha num dos componentes vitais de um motor a 4 tempos: o veio de manivelas. Esta situação já tem sido bastante observada na indústria mecânica naval, nomeadamente na Marinha Portuguesa e, devido ao seu elevado grau de importância no desempenho de qualquer motor, decidiu-se focar o trabalho desta tese no estudo dos motores a diesel S.E.M.T Pielstick das unidades navais da Marinha Portuguesa, mais especificamente das corvetas da classe “João Coutinho” e da classe “ Baptista de Andrade”, devido ao historial de ocorrência de falhas no veio de manivelas nesta classe de navios e em outras da Marinha Portuguesa. Para efetuar este estudo, utilizaram-se todos os dados relativos ao historial de ocorrências de falhas destes motores, bem como todos os dados disponíveis do fabricante destes motores, por forma a reproduzir da forma mais fiável possível um modelo tridimensional do veio de manivelas no programa de modelação informática CAD Solidworks®, e possibilitar a análise cinemática do veio de manivelas. Desta forma, foi possível simular as condições de funcionamento do motor, assim como analisar e determinar a causa de falha do veio de manivelas, visando prolongar a vida útil dos veios de manivelas, contribuindo não só para menores custos de manutenção mas também para o aumento da operacionalidade destes navios.

Electro-optical characterization of IC compatible microcantilevers

The aim of this work is to simulate and optically characterize the piezoelectric performance of complementary metal oxide semiconductor (CMOS) compatible microcantilevers based on aluminium nitride (AlN) and manufactured at room temperature. This study should facilitate the integration of piezoelectric micro-electro-mechanical systems (MEMS) such as microcantilevers, in CMOS technology. Besides compatibility with standard integrated circuit manufacturing procedures, low temperature processing also translates into higher throughput and, as a consequence, lower manufacturing costs. Thus, the use of the piezoelectric properties of AlN manufactured by reactive sputtering at room temperature is an important step towards the integration of this type of devices within future CMOS technology standards. To assess the reliability of our fabrication process, we have manufactured arrays of free-standing microcantilever beams of variable dimension and studied their piezoelectric performance. The characterization of the first out-of-plane modes of AlN-actuated piezoelectric microcantilevers has been carried out using two optical techniques: laser Doppler vibrometry (LDV) and white light interferometry (WLI). In order to actuate the cantilevers, a periodic chirp signal in certain frequency ranges was applied between the device electrodes. The nature of the different vibration modes detected has been studied and compared with that obtained by a finite element model based simulation (COMSOL Multiphysics), showing flexural as well as torsional modes. The correspondence between theoretical and experimental data is reasonably good, probing the viability of this high throughput and CMOS compatible fabrication process. To complete the study, X-ray diffraction as well as d33 piezoelectric coefficient measurements were also carried out.