Neste trabalho, o escoamento incompressível ao redor de um cilindro fixo é estudado numericamente através do método de elementos finitos. Foram realizadas simulações bidimensionais no domínio do tempo, com números de Reynolds variando entre 100 e 600, utilizando-se, para tanto, malhas não-estruturadas com elementos triangulares. Pretende-se aqui analisar a eficácia da solução das equações de Navier-Stokes utilizando o método das penalidades, meio pelo qual o acoplamento pressão-velocidade foi tratado. Avalia-se a convergência da solução para diferentes valores do fator de penalidade e sugere-se um método para estimá-lo. Analisa-se, ainda, a sensibilidade da resposta à utilização da matriz de inércia nos formatos consistente e concentrada. Por fim, é realizada a comparação dos coeficientes de arrasto médio, flutuação do coeficiente de sustentação e número de Strouhal obtidos neste trabalho com resultados de outras publicações. / In this work, the incompressible flow around a stationary cylinder is investigated by using the Finite Element Method. Two-dimensional simulations in time domain have been carried out, with Reynolds number varying from 100 to 600, using non-structured meshes with triangular elements. The aim of this work is to analyze the efficiency of Penalty Methods, which is the way that the velocity-pressure coupling problem is treated here, in Navier-Stokes equations solution. The solution convergence from different values of penalty parameter is evaluated and it is suggested a method to estimate it. In addition, it is studied the sensibilty of response when using the mass matrix in consistent or lumped format. At last, a comparison between average drag coefficient, fluctuating lift and Strouhal number obtained here and those found in other publications is shown.

O trabalho de doutorado teve como objetivo desenvolver um modulador óptico baseado em filme magnetostrictivo para aplicação como magnetômetro. A modelagem e simulação do dispositivo foi realizada utilizando software com cálculo por método dos elementos finitos (MEF) e teve como finalidade auxiliar iterativamente os processos de projeto e fabricação do modulador. A originalidade da proposta baseia-se na construção de um guia de onda em óptica integrada recoberto por um filme magnetostrictivo para permitir a modulação , via efeito elasto-óptico, da onda guiada pela aplicação de campos magnéticos externos. O campo magnético aplicado provoca a deformação o material magnetostrictivo que induz uma modificação no perfil de esforço aplicado ao substrato. O substrato tem suas propriedades ópticas alteradas devido ao efeito elasto-óptico, o que provoca mudanças nas propriedades da luz transmitida. O trabalho tem seu início com o estudo e a caracterização de filmes magnetostrictivos de Tb25F275 e Tb23Co77 depositados por sputtering sobre substratos de silício. Uma técnica para preparação das amostras e medição da magnetostricção foi estabelecida e os coeficientes de magnetostrição dos filmes foram determinados a partir das medições diretas dos deslocamentos das amostras, em função dos campos magnéticos aplicados, utilizando a técnica de Atomic Force Microscopy (AFM). Os resultados experimentais obtidos permitiram a realização de simulações por MEF para verificação dos modos de guiamento da luz gerados pelo perfil de esforços induzidos termicamente no processo de deposição do filme magnetostrictivo sobre substrato de B12GeO4 (BGO). Foi modelado e simulado também o efeito da aplicação do campo magnético sobre o guia óptico obtido inicialmente pelo efeito de esforço térmico. No resultado das simulações foi possível verificar as alterações do índice de refração efetivo e da intensidade óptica do modo guiado em função de campos magnéticos aplicados ao modulador. Ao final do trabalho realizaram-se a fabricação de alguns protótipos. Os resultados das caracterizações dos moduladores construídos permitirão, no futuro, ajustar os modelos de simulação elaborados. / The doctoral work aims are the development and simulation of an optical modulator based on the effect of magnetostriction for application as magnetometer. The multiphysics simulations were performed using the Finite Elements Method (FEM). In the manufacturing process of optical modulator integrated optics techniques were applied. The originality of the proposal is based on the construction of an integrated optical waveguide covered by a magnetostrictive film to allow the modulation of the guided wave, through the elasto-optical effect, by the application of external magnetic fields. The applied magnetic field causes deformation of the magnetostrictive material that induces a modification of the stress profile produced in substrate. The substrate has its optical properties altered due to the elasto-optical effect, which causes changes in the properties of transmitted light. The work begins with the study and characterization of TbFe and TbCo2 magnetostrictive films deposited by sputtering on silicon substrates. A method for sample preparation and measurement of magnetostriction was established. The magnetostrictive coefficient of the films was determined from the direct measurement of displacements of samples by AFM technique for magnetic fields applied. The experimental results obtained allowed to perform MEF simulations to verify the light guided modes generated by the profile of thermally induced stress created by deposition process of magnetostrictive film on B12GeO4 (BGO) substrate. It was also modeled and simulated the effect of the application of magnetic field on the optical guide obtained initially by the effect of thermal stress. In simulation results, it was possible to verify the changes of effective refractive index and optical intensity of guided modes as functions of magnetic fields applied to the modulator. At the end of the work, some prototypes were fabricated. The results of characterizations of the built modulators will allow, in the future, adjustments in simulation models.

<p>Föreliggande avhandling omfattar en fallstudie av utmattningsrisken för järnvägsbroarna mellan Stockholm Central och Söder Mälarstrand. Ett stort antal utmattningssprickor i broarna över Söderström och Söder Mälarstrand har sedan länge varit kända. I samband med en bärighetsutredning av sträckan har ett stort antal anslutningspunkter i samtliga broar identifierats som kritiska avseende utmattningsrisk. Sträckan är den mest trafikerade i Sverige och de förenklade beräkningsmetoder som anges i Banverkets föreskrift för bärighetsberäkningar är i vissa fall inte tillämpliga.</p><p>De formella bärighetsberäkningarna har visat på flerfaldiga överskridanden i utnyttjandegrad gällande utmattning. I de punkter som beräkningarna visar störst risk för utmattning har inga sprickor påträffats, trots omfattande förbandsbesiktningar. Beräkningarna baseras på typiserade spänningskollektiv och givet antal spännings­växlingar, oberoende av den aktuella trafikmängden. Enligt bärighetsnormen ges möjlighet att beakta den verkliga trafikmängden baserat på historiskt bruttotonnage och typiserade trafiklaster. Dessa beräkningar har utförts och visar på ännu större risk för utmattning, jämfört med den förenklade metoden. I syfte att på ett ännu mer nyanserat sätt uppskatta utmattningsrisken har en metod använts, baserat på uppskattning av den faktiska trafikmängden och dess fördelningar. Beräkningarna baseras på trafikdata tillgänglig från Banverket. Metoden är tillämplig på andra järnvägsbroar på andra sträckor.</p><p>Utöver teoretiska beräkningar har töjningsmätningar utförts. Under 2006 utförde dåvarande Carl Bro AB töjningsmätningar på bro över Söder Mälarstrand och viadukt söder om Söder Mälarstrand, avseende några enskilda tågöverfarter. Under 2008 utförde KTH Brobyggnad omfattande töjningsmätningar på bro över Söderström, avseende all trafik under en månads tid.</p> / <p>The present thesis comprises a case study of the risk of fatigue of the railway bridges chaining between Stockholm Central Station and the district Söder Mälarstrand. A large number of fatigue-related cracks in the bridges at Söderström and Söder Mälarstrand have been known for a long time. During a capacity assessment of the current bridges, a large number of connections have been identified as critical concerning fatigue resistance. The route is the most frequent in all of Sweden and the simplified methods of fatigue assessment defined by Banverket may not always be applicable.</p><p>A conventional capacity assessment has shown numerous exceeds in fatigue resistance, using the stated safety margins. No fatigue cracks have been identified at the locations showing the largest theoretical risk of fatigue, in spite of extensive investigations. The conventional calculations are based on a uniform stress collective and a fixed number of stress cycles, independent of the actual traffic volume. According to the regulations stated by Banverket, the assessment may optionally be performed using historical data of the gross tonnage and standardised traffic loading. Such calculations have been undertaken and show even greater risk of fatigue, compared to the conventional assessment. To estimate the risk of fatigue in more detail, a method has been used, based on estimations of the real traffic volume and its distributions. The analysis is based on available data of the traffic volume and may be applied to other railway bridges on other locations.</p><p>Besides theoretical analyses, field measurements have been performed. In 2006, former Carl Bro AB carried out strain gauge measurements on the bridge passing Söder Mälarstrand and the viaduct south of Söder Mälarstrand. The measurements comprised a small amount of individual train passages. During 2008, the division of Structural Design and Bridges at KTH performed an extensive field measurement programme on the bridge passing Söderström. Continuous measurements collecting data of all traffic during a period of one month was performed.</p>

In the last years multilevel preconditioners like BPX became more and more popular for solving second-order elliptic finite element discretizations by iterative methods. P. Oswald has adapted these methods for discretizations of the fourth order biharmonic problem by rectangular conforming Bogner-Fox-Schmidt elements and nonconforming Adini elements and has derived optimal estimates for the condition numbers of the preconditioned linear systems. In this paper we generalize the results from Oswald to the construction of BPX and Multilevel Diagonal Scaling (MDS-BPX) preconditioners for the elasticity problem of thin smooth shells of arbitrary forms where we use Koiter's equations of equilibrium for an homogeneous and isotropic thin shell, clamped on a part of its boundary and loaded by a resultant on its middle surface. We use the two discretizations mentioned above and the preconditioned conjugate gradient method as iterative method. The parallelization concept is based on a non-overlapping domain decomposition data structure. We describe the implementations of the multilevel preconditioners. Finally, we show numerical results for some classes of shells like plates, cylinders, and hyperboloids.

Coupled systems are ubiquitous in modern engineering and science. Such systems can encompass fluid dynamics, structural mechanics, chemical species transport and electrostatic effects among other components, all of which can be coupled in many different ways. In addition, such models are usually multiscale, making their numerical simulation challenging, and necessitating the use of adaptive modeling techniques. The multiscale, multiphysics models of electrosomotic flow (EOF) constitute a particularly challenging coupled flow system. A special feature of such models is that the coupling between the electric physics and hydrodynamics is via the boundary. Numerical simulations of coupled systems are typically targeted towards specific Quantities of Interest (QoIs). Adjoint-based approaches offer the possibility of QoI targeted adaptive mesh refinement and efficient parameter sensitivity analysis. The formulation of appropriate adjoint problems for EOF models is particularly challenging, due to the coupling of physics via the boundary as opposed to the interior of the domain. The well-posedness of the adjoint problem for such models is also non-trivial. One contribution of this dissertation is the derivation of an appropriate adjoint problem for slip EOF models, and the development of penalty-based, adjoint-consistent variational formulations of these models. We demonstrate the use of these formulations in the simulation of EOF flows in straight and T-shaped microchannels, in conjunction with goal-oriented mesh refinement and adjoint sensitivity analysis. Complex computational models may exhibit uncertain behavior due to various reasons, ranging from uncertainty in experimentally measured model parameters to imperfections in device geometry. The last decade has seen a growing interest in the field of Uncertainty Quantification (UQ), which seeks to determine the effect of input uncertainties on the system QoIs. Monte Carlo methods remain a popular computational approach for UQ due to their ease of use and "embarassingly parallel" nature. However, a major drawback of such methods is their slow convergence rate. The second contribution of this work is the introduction of a new Monte Carlo method which utilizes local sensitivity information to build accurate surrogate models. This new method, called the Local Sensitivity Derivative Enhanced Monte Carlo (LSDEMC) method can converge at a faster rate than plain Monte Carlo, especially for problems with a low to moderate number of uncertain parameters. Adjoint-based sensitivity analysis methods enable the computation of sensitivity derivatives at virtually no extra cost after the forward solve. Thus, the LSDEMC method, in conjuction with adjoint sensitivity derivative techniques can offer a robust and efficient alternative for UQ of complex systems. The efficiency of Monte Carlo methods can be further enhanced by using stratified sampling schemes such as Latin Hypercube Sampling (LHS). However, the non-incremental nature of LHS has been identified as one of the main obstacles in its application to certain classes of complex physical systems. Current incremental LHS strategies restrict the user to at least doubling the size of an existing LHS set to retain the convergence properties of LHS. The third contribution of this research is the development of a new Hierachical LHS algorithm, that creates designs which can be used to perform LHS studies in a more flexibly incremental setting, taking a step towards adaptive LHS methods. / text

No description available.

Made available in DSpace on 2014-06-11T19:25:21Z (GMT). No. of bitstreams: 0 Previous issue date: 2004-01-23Bitstream added on 2014-06-13T19:53:02Z : No. of bitstreams: 1 waidemam_l_me_ilha.pdf: 2247815 bytes, checksum: aab885504632d6169da23c5ac9796141 (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / Este trabalho tem como objetivo principal analisar o comportamento dinâmico de estruturas laminares planas com carregamento perpendicular ao plano médio, em particular as placas delgadas, utilizando-se, para isso, a teoria clássica de flexão de placas e a discretização estrutural feita com os elementos finitos triangulares e retangulares trabalhando em conjunto e em separado. Na dedução das matrizes de rigidez e de massas dos elementos finitos em questão utiliza-se a formulação com parâmetros generalizados e com coordenadas homogêneas, cujas funções aproximadoras contêm nove e doze monômios, respectivamente, extraídos do polinômio algébrico cúbico em “x” e “y”. Para a consideração do amortecimento utiliza-se o Método de Rayleigh e para a integração numérica ao longo do tempo utiliza-se o Método de Newmark, via algoritmo previsor / corretor. Ao final deste são elaborados vários exemplos elucidativos visando uma análise quantitativa e qualitativa dos resultados obtidos. / In this work the dynamic behavior of plane laminate structures, with load applied perpendicularly to the middle plan, has been analyzed. The classic theory of bending plates and structural subdivision - done with triangular and rectangular finite elements working together as well as in separate – are used to study thin plates. The formulation employing generalized parameters and homogeneous coordinates, using approximating functions containing nine and twelve terms starting from the cubic algebraic polynomial in Cartesian coordinates x and y, is used to obtain the stiffness and mass matrices for the triangular and rectangular finite element, respectively. The Rayleigh Method is used to take into account the structural dumping while the Newmark Method is used to perform the numeric integration in the time, by means of predictor / corrector scheme. Additionally, several elucidating examples are elaborated in order to analyze the final results.

Made available in DSpace on 2014-06-11T19:25:27Z (GMT). No. of bitstreams: 0 Previous issue date: 2006-02-03Bitstream added on 2014-06-13T19:12:11Z : No. of bitstreams: 1 bagagli_rp_me_bauru.pdf: 1886943 bytes, checksum: a693fc27655c63ad68c8b6f960ca899f (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / Com a crise energética recente, houve uma nova conscientização da necessidade de utilização mais racional da energia. Desta feita, uma série de pesquisas com fontes alternativas de energia, que vinham sendo preteridos em função da impressão que a crise energética do início da década de 1970 havia passado, tem ganho nova força. Dentre todas as alternativas para aproveitamento de energia solar, uma das mais utilizadas são os chamados sistemas domésticos de aquecimento de água. Este tipo de sistema, entretanto, ainda é complexo, constituído de uma série de dutos e conexões entre coletor e tanque armazenador, que contribuem para o elevado custo destes dispositivos. Uma alternativa para otimizar o custo final é o uso dos sistemas solar integrados coletor-tanque (ICS). Neste trabalho, foi avaliado o processo de movimentação natural do fluido em uma das geometrias mais comuns de sistemas deste tipo: a trapezoidal. Foi aplicada a condição de fluxo de calor constante na face inclinada para avaliação... / The recent energy crisis has developed a new conscience for necessity of rational energy use. Several works treating about renewable energy was stopped in past based in a false idea that the 1970's energy crisis was finished. Nowadays, these works have been retaken with the large use of solar energy in the Solar Domestic Hot Water Systems. However, this device is quite complex and has several components like pipes and fittings coupling solar collectors and storage tanks. This characteristic makes it an expensive system and bring difficulties for his large use. An alternative to turn it cheaper is the construction of a device with solar collector and storage tank integrated in one single component (ICS). In this work was done an evaluation of free convection process in a common geometry of this device: the trapezoidal shape. For this analysis, a constant heat flux condition was applied to the inclined face for evaluation of free convection process. Numerical results... (Complete abstract click electronic access below)

In this dissertation, we study the stationary heat equation in a heterogeneous tridimensional plate, using a "dimension reduction" techinique called hierarchical modelling and we generate model the original problem in a two-dimensional domain. To estimate the error modelling, we develop an asymptotic expansion for the original problem solution and for the aproximate solution. Comparing both solutions with their own asymptotic expansions, we obtain an estimative of the error modelling. We perform some computational experiments, using the Residual Free Bubbles (RFB) Method and the Multiscale Finite Element Method for the diffusion problem and for the diffusion-reaction problem in a two-dimensional domain, with small parameters. Finally, we extend the numeric solutions found the original tridimensional problem. / Neste trabalho, estudamos a equação do calor estacionária em uma placa heterogênea tridimensional. Para a modelagem deste problema, utilizamos uma técnica de redução de dimensão conhecida por Modelagem Hierárquica. Desta forma, geramos um modelo para o problema original em um domínio bidimensional. Com o objetivo de estimar o erro de modelagem, desenvolvemos a expansão assintótica da solução do problema original e da solução aproximada. Comparando as soluções com suas respectivas expansões assintóticas, obtemos uma estimativa para o erro de modelagem. Realizamos alguns experimentos computacionais, desenvolvendo o método Residual Free Bubbles (RFB) e o método de Elementos Finitos Multiescala (MEFM) para o problema de difusão e para o problema de difusão-reação em um domínio bidimensional, com parâmetros pequenos. Com base nestes experimentos, encontramos algumas soluções numéricas para o problema da placa tridimensional.

Orientador: Rogério de Oliveira Rodrigues / Banca: Wilson Sergio Venturini / Banca: Mônica Pinto Barbosa / Resumo: Este trabalho tem como objetivo principal analisar o comportamento dinâmico de estruturas laminares planas com carregamento perpendicular ao plano médio, em particular as placas delgadas, utilizando-se, para isso, a teoria clássica de flexão de placas e a discretização estrutural feita com os elementos finitos triangulares e retangulares trabalhando em conjunto e em separado. Na dedução das matrizes de rigidez e de massas dos elementos finitos em questão utiliza-se a formulação com parâmetros generalizados e com coordenadas homogêneas, cujas funções aproximadoras contêm nove e doze monômios, respectivamente, extraídos do polinômio algébrico cúbico em "x" e "y". Para a consideração do amortecimento utiliza-se o Método de Rayleigh e para a integração numérica ao longo do tempo utiliza-se o Método de Newmark, via algoritmo previsor / corretor. Ao final deste são elaborados vários exemplos elucidativos visando uma análise quantitativa e qualitativa dos resultados obtidos. / Abstract: In this work the dynamic behavior of plane laminate structures, with load applied perpendicularly to the middle plan, has been analyzed. The classic theory of bending plates and structural subdivision - done with triangular and rectangular finite elements working together as well as in separate - are used to study thin plates. The formulation employing generalized parameters and homogeneous coordinates, using approximating functions containing nine and twelve terms starting from the cubic algebraic polynomial in Cartesian coordinates "x" and "y", is used to obtain the stiffness and mass matrices for the triangular and rectangular finite element, respectively. The Rayleigh Method is used to take into account the structural dumping while the Newmark Method is used to perform the numeric integration in the time, by means of predictor / corrector scheme. Additionally, several elucidating examples are elaborated in order to analyze the final results. / Mestre