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.

Adaptivity and a posteriori error control for bifurcation problems I: the Bratu problem

This article is concerned with the numerical detection of bifurcation points of nonlinear partial differential equations as some parameter of interest is varied. In particular, we study in detail the numerical approximation of the Bratu problem, based on exploiting the symmetric version of the interior penalty discontinuous Galerkin finite element method. A framework for a posteriori control of the discretization error in the computed critical parameter value is developed based upon the application of the dual weighted residual (DWR) approach. Numerical experiments are presented to highlight the practical performance of the proposed a posteriori error estimator.

A class of domain decomposition preconditioners for hp-discontinuous Galerkin finite element methods

In this article we address the question of efficiently solving the algebraic linear system of equations arising from the discretization of a symmetric, elliptic boundary value problem using hp-version discontinuous Galerkin finite element methods. In particular, we introduce a class of domain decomposition preconditioners based on the Schwarz framework, and prove bounds on the condition number of the resulting iteration operators. Numerical results confirming the theoretical estimates are also presented.

High-order hp-adaptive discontinuous Galerkin finite element methods for compressible fluid flows

This article is concerned with the construction of general isotropic and anisotropic adaptive strategies, as well as hp-mesh refinement techniques, in combination with dual-weighted-residual a posteriori error indicators for the discontinuous Galerkin finite element discretization of compressible fluid flow problems.

Influence of implant-abutment angulations and crown material on stress distribution on central incisor: a 3D FEA

Aim: To investigate the effect of implant-abutment angulation and crown material on stress distribution of central incisors. Finite element method was used to simulate the clinical situation of a maxillary right central incisor restored by two different implant-abutment angulations, 15° and 25°, using two different crown materials (IPS E-Max CAD and zirconia). Methods: Two 3D finite element models were specially prepared for this research simulating the abutment angulations. Commercial engineering CAD/CAM package was used to model crown, implant abutment complex and bone (cortical and spongy) in 3D. Linear static analysis was performed by applying a 178 N oblique load. The obtained results were compared with former experimental results. Results: Implant Von Mises stress level was negligibly changed with increasing abutment angulation. The abutment with higher angulation is mechanically weaker and expected to fail at lower loading in comparison with the steeper one. Similarly, screw used with abutment angulation of 25° will fail at lower (about one-third) load value the failure load of similar screw used with abutment angulated by 15°. Conclusions: Bone (cortical and spongy) is insensitive to crown material. Increasing abutment angulation from 15° to 25°, increases stress on cortical bone by about 20% and reduces it by about 12% on spongy bone. Crown fracture resistance is dramatically reduced by increasing abutment angulation. Zirconia crown showed better performance than E-Max one.

Avances en el proyecto “Determinación de la forma exacta de la tierra a través de MEMS en un Nano-Satélite”

El presente trabajo muestra el avance que se tiene en el proyecto denominado “Diseño y Construcción de un Nano-Satélite”, que se lleva a cabo en la ESIME U.P.-Ticoman del Instituto Politécnico Nacional. Este proyecto tiene como objetivo principal introducir a los alumnos en el que hacer espacial, se trata de un satélite con un peso de un kilogramo y es un cubo de 10 cm de lado y cuya carga útil tiene MEMS (sistema de micro-electromecánicos) que harán las veces de péndulo electrónico, lo anterior y los armónicos gravitacionales permitirán determinar la forma de la tierra. Se presenta el fundamento matemático, los componentes del satélite, el diseño preliminar de la estructura, la manufactura de la estructura, los resultados de un ensaye estructural destructivo y el análisis estático por elementos finitos.

High Frequency Class DE Converter Using a Multilayer Coreless PCB Transformer

In modern power electronics equipment, it is desirable to design a low profile, high power density, and fast dynamic response converter. Increases in switching frequency reduce the size of the passive components such as transformers, inductors, and capacitors which results in compact size and less requirement for the energy storage. In addition, the fast dynamic response can be achieved by operating at high frequency. However, achieving high frequency operation while keeping the efficiency high, requires new advanced devices, higher performance magnetic components, and new circuit topology. These are required to absorb and utilize the parasitic components and also to mitigate the frequency dependent losses including switching loss, gating loss, and magnetic loss. Required performance improvements can be achieved through the use of Radio Frequency (RF) design techniques. To reduce switching losses, resonant converter topologies like resonant RF amplifiers (inverters) combined with a rectifier are the effective solution to maintain high efficiency at high switching frequencies through using the techniques such as device parasitic absorption, Zero Voltage Switching (ZVS), Zero Current Switching (ZCS), and a resonant gating. Gallium Nitride (GaN) device technologies are being broadly used in RF amplifiers due to their lower on- resistance and device capacitances compared with silicon (Si) devices. Therefore, this kind of semiconductor is well suited for high frequency power converters. The major problems involved with high frequency magnetics are skin and proximity effects, increased core and copper losses, unbalanced magnetic flux distribution generating localized hot spots, and reduced coupling coefficient. In order to eliminate the magnetic core losses which play a crucial role at higher operating frequencies, a coreless PCB transformer can be used. Compared to the conventional wire-wound transformer, a planar PCB transformer in which the windings are laid on the Printed Board Circuit (PCB) has a low profile structure, excellent thermal characteristics, and ease of manufacturing. Therefore, the work in this thesis demonstrates the design and analysis of an isolated low profile class DE resonant converter operating at 10 MHz switching frequency with a nominal output of 150 W. The power stage consists of a class DE inverter using GaN devices along with a sinusoidal gate drive circuit on the primary side and a class DE rectifier on the secondary side. For obtaining the stringent height converter, isolation is provided by a 10-layered coreless PCB transformer of 1:20 turn’s ratio. It is designed and optimized using 3D Finite Element Method (FEM) tools and radio frequency (RF) circuit design software. Simulation and experimental results are presented for a 10-layered coreless PCB transformer operating in 10 MHz.

Modelling and simulations of ductile iron solidification-induced variations in mechanical behaviour on component and microstructural level

The mechanical behaviour and performance of a ductile iron component is highly dependent on the local variations in solidification conditions during the casting process. Here we show a framework which combine a previously developed closed chain of simulations for cast components with a micro-scale Finite Element Method (FEM) simulation of the behaviour and performance of the microstructure. A casting process simulation, including modelling of solidification and mechanical material characterization, provides the basis for a macro-scale FEM analysis of the component. A critical region is identified to which the micro-scale FEM simulation of a representative microstructure, generated using X-ray tomography, is applied. The mechanical behaviour of the different microstructural phases are determined using a surrogate model based optimisation routine and experimental data. It is discussed that the approach enables a link between solidification- and microstructure-models and simulations of as well component as microstructural behaviour, and can contribute with new understanding regarding the behaviour and performance of different microstructural phases and morphologies in industrial ductile iron components in service.

Models for reliability estimation of HVDC cable systems

The work carried out in this thesis aims at: - studying – in both simulative and experimental methods – the effect of electrical transients (i.e., Voltage Polarity Reversals VPRs, Temporary OverVoltages TOVs, and Superimposed Switching Impulses SSIs) on the aging phenomena in HVDC extruded cable insulations. Dielectric spectroscopy, conductivity measurements, Fourier Transform Infra-Red FTIR spectroscopy, and space charge measurements show variation in the insulating properties of the aged Cross-Linked Polyethylene XLPE specimens compared to non-aged ones. Scission in XLPE bonds and formation of aging chemical bonds is also noticed in aged insulations due to possible oxidation reactions. The aged materials show more ability to accumulate space charges compared to non-aged ones. An increase in both DC electrical conductivity and imaginary permittivity has been also noticed. - The development of life-based geometric design of HVDC cables in a detailed parametric analysis of all parameters that affect the design. Furthermore, the effect of both electrical and thermal transients on the design is also investigated. - The intrinsic thermal instability in HVDC cables and the effect of insulation characteristics on the thermal stability using a temperature and field iterative loop (using numerical methods – Finite Difference Method FDM). The dielectric loss coefficient is also calculated for DC cables and found to be less than that in AC cables. This emphasizes that the intrinsic thermal instability is critical in HVDC cables. - Fitting electrical conductivity models to the experimental measurements using both models found in the literature and modified models to find the best fit by considering the synergistic effect between field and temperature coefficients of electrical conductivity.

A Trans-dimensional inversion algorithm to model deformation sources with unconstrained shape in finite element domains

Ground deformation provides valuable insights on subsurface processes with pattens reflecting the characteristics of the source at depth. In active volcanic sites displacements can be observed in unrest phases; therefore, a correct interpretation is essential to assess the hazard potential. Inverse modeling is employed to obtain quantitative estimates of parameters describing the source. However, despite the robustness of the available approaches, a realistic imaging of these reservoirs is still challenging. While analytical models return quick but simplistic results, assuming an isotropic and elastic crust, more sophisticated numerical models, accounting for the effects of topographic loads, crust inelasticity and structural discontinuities, require much higher computational effort and information about the crust rheology may be challenging to infer. All these approaches are based on a-priori source shape constraints, influencing the solution reliability. In this thesis, we present a new approach aimed at overcoming the aforementioned limitations, modeling sources free of a-priori shape constraints with the advantages of FEM simulations, but with a cost-efficient procedure. The source is represented as an assembly of elementary units, consisting in cubic elements of a regular FE mesh loaded with a unitary stress tensors. The surface response due to each of the six stress tensor components is computed and linearly combined to obtain the total displacement field. In this way, the source can assume potentially any shape. Our tests prove the equivalence of the deformation fields due to our assembly and that of corresponding cavities with uniform boundary pressure. Our ability to simulate pressurized cavities in a continuum domain permits to pre-compute surface responses, avoiding remeshing. A Bayesian trans-dimensional inversion algorithm implementing this strategy is developed. 3D Voronoi cells are used to sample the model domain, selecting the elementary units contributing to the source solution and those remaining inactive as part of the crust.

Planning preoperatorio basato sulla ricostruzione di un avatar muscolo-scheletrico customizzato

L’obbiettivo dell’elaborato che segue è quello di fornire a medici e chirurghi uno strumento di valutazione dei trattamenti di cura, per pazienti affetti da malattie del sistema muscolo scheletrico e della deambulazione. Lo studio incomincia dall’analisi del movimento del paziente tramite strumenti di gate analysis. Impiegando come modello la ricostruzione dell’apparato muscolo scheletrico del paziente, tramite sistemi di modellazione 3Dda TAC (Tomografia Assiale Computerizzata) o risonanza magnetica, realizzando un modello personalizzato che verrà analizzato tramite il software OpenSim. Tramite lo studio del movimento andremo successivamente ad identificare le forze articolari agenti sul ginocchio e grazie alla modellazione FEM (Finite Element Method) verrà realizzato uno strumento di valutazione della strategia operatoria ottimale per il paziente, simulando l’operazione del ginocchio, che agisce sulle forze muscolari variando le forze agenti sull’articolazione e come queste variazioni modificano l’andamento delle tensioni nell’articolazione.

Analisi vibroacustica di un sistema di scarico per monoposto Formula SAE

Lo scopo del presente lavoro è quello di determinare una metodologia di analisi vibroacustica di validità generale applicabile a tutti i casi nei quali la forzante sia di tipo vibrazionale. Nello specifico si è analizzato il comportamento strutturale e acustico di un impianto di scarico di una monoposto Formula SAE. Ricorrendo all’analisi FEM (Finite Element Method) è possibile determinare e quantificare gli effetti dannosi causati dalle vibrazioni già nella fase di prototipazione permettendo una sostanziale riduzione dei tempi e costi. La determinazione del comportamento strutturale del modello alle vibrazioni è iniziata dall’analisi modale, grazie alla quale sono state determinate le frequenze naturali e i modi propri dell’impianto di scarico. Successivamente, l’analisi FRF (Frequency Response Function) ha permesso di conoscere la risposta del nostro sistema ad una forzante imposta mettendo in luce le diverse criticità strutturali. Con il presupposto di ottenere delle condizioni di carico che fossero il più vicine possibili alle normali condizioni operative si è impostata un’analisi PSD (Power Spectral Density). Per concludere la prima parte dell’analisi si è reso necessario indagare anche il comportamento a fatica vibrazionale, valutando in questo modo le zone soggette a vita finita e quindi le prime a cedere in fase di esercizio. La parte finale è stata dedicata all’analisi del rumore. Dall’analisi FRF si è determinata la SPL (Sound Pressure Level) ottenendo come output un valore di pressione sonora prodotto dall’effetto della propagazione delle onde di pressione generate dalla vibrazione strutturale dell’impianto di scarico. Infine, l’analisi di Transmission Loss ha permesso di valutare l’efficacia della geometria del silenziatore sulla riduzione del livello acustico generato dal transito dei gas di scarico alle diverse frequenze.

Keratoconus prediction using a finite element model of the cornea with local biomechanical properties

PURPOSE: The ability to predict and understand which biomechanical properties of the cornea are responsible for the stability or progression of keratoconus may be an important clinical and surgical tool for the eye-care professional. We have developed a finite element model of the cornea, that tries to predicts keratoconus-like behavior and its evolution based on material properties of the corneal tissue. METHODS: Corneal material properties were modeled using bibliographic data and corneal topography was based on literature values from a schematic eye model. Commercial software was used to simulate mechanical and surface properties when the cornea was subject to different local parameters, such as elasticity. RESULTS: The simulation has shown that, depending on the corneal initial surface shape, changes in local material properties and also different intraocular pressures values induce a localized protuberance and increase in curvature when compared to the remaining portion of the cornea. CONCLUSIONS: This technique provides a quantitative and accurate approach to the problem of understanding the biomechanical nature of keratoconus. The implemented model has shown that changes in local material properties of the cornea and intraocular pressure are intrinsically related to keratoconus pathology and its shape/curvature.