This work is a numerical study of Burgers' equation with Robin's boundary conditions. The goal is to determine the behavior of the solutions in the limiting cases of Dirichlet and Neumann boundary conditions. We develop and test two separate finite element and Galerkin schemes. The Galerkin/Conservation method is shown to give better results and is then used to compute the response as the Robin's Boundary conditions approach both the Dirichlet and Neumann boundary conditions. Burgers' equation is treated as a perturbation of the linear heat equation with the appropriate realistic constants. The goal is to determine if the use of the Robin's boundary conditions to approximate Dirichlet and Neumann boundary conditions affords any advantage over schemes that employ only "exact" Dirichlet or Neumann boundary conditions. Our finite element results indicate that solutions using appropriate Robin's boundary conditions approach the same solutions obtained by "exact" Dirichlet or Neumann boundary conditions. This allows us to obtain realistic solutions in some cases where the other schemes had previously failed. / Master of Science

Cette étude concerne la simulation de systèmes thermiques comportant plusieurs domaines à l’aide d’une méthode d’ordre réduit adaptée au mono-corps. Les modèles modaux réduits sont construits à partir de modes de branche. Ceux-ci sont déterminés à l’aide d’un modèle détaillé et permettent la prise en compte de non-linéarités. La réduction est effectuée par la technique de l’amalgame modal découplé qui permet de façon automatique et immédiate d’obtenir une base réduite qui ne comporte qu’un faible nombre de modes à partir de la base initiale. Le couplage entre les différents domaines s’effectue par l’intermédiaire d’une résistance thermique de contact. On montre la nécessité de faire intervenir un terme supplémentaire de pénalisation de saut de flux qui vient améliorer les résultats. Les simulations numériques effectuées sur des cas tests (microprocesseur et radiateur en 2D, bloc métallique avec cartouches chauffantes en 3D) montrent la pertinence de la méthode. / This study concerns the simulation of thermal systems with multiple fields with a reduced-order method suited to a single body. Reduced models are constructed from modal branch eigenmodes. They are determined using a detailed model and allow the inclusion of non-linearities. The reduction is carried out by the simplified amalgam method which allows an automatic and immediate way to obtain a reduced basis which contains only a small number of modes from the original basis. The coupling between the different areas is carried out through a thermal contact resistance. It shows the need to involve an additional flux jump penalty term to improve the results. The numerical simulations carried out on test cases (microprocessor and radiator in 2D, metal block with hot cartridges in 3D) show the relevance of the method.

<p>Techniques for a posteriori error estimation for finite element approximations of an elliptic partial differential equation are studied.This extends previous work on localized error control in finite element methods for linear elasticity.The methods are then applied to the problem of homogenization of periodic structures. In particular, error estimates for the effective elastic properties are obtained. The usefulness of these estimates is twofold.First, adaptive methods using mesh refinements based on the estimates can be constructed.Secondly, one of the estimates can give reasonable measure of the magnitude ofthe error. Numerical examples of this are given.</p>

Techniques for a posteriori error estimation for finite element approximations of an elliptic partial differential equation are studied.This extends previous work on localized error control in finite element methods for linear elasticity.The methods are then applied to the problem of homogenization of periodic structures. In particular, error estimates for the effective elastic properties are obtained. The usefulness of these estimates is twofold.First, adaptive methods using mesh refinements based on the estimates can be constructed.Secondly, one of the estimates can give reasonable measure of the magnitude ofthe error. Numerical examples of this are given.

No description available.

[pt] Métodos estabilizados de elementos finitos são propostos e analisados para problemas de fluidos, com particular ênfase nas equações de Navier-Stokes incomprenssível. Após a apresentação da mecânica dos escoamentos dos fluidos, introduzimos no Capítulo 3, no contexto de problema de Stokes, as dificuldadas numéricas associadas ao método de Galerkin em problemas de fluidos e simulamos em sucesso alguns escoamentos lentos através de formulações finitos para estabilizadas. No capítulo 4, propomos uma nova formulação de elementos finitos para a equação da energia, mais precisamente para o modelo da advecção-difusão do calor. Graças a um novo desenho do parâmetro de estabilidade T, o qual permite adicionar difusão às regiões advectivas e difusivas-dominadas do escoamento de maneira diferemciada, obtivemos um bom desempenho novo método mesmo em situações de altíssimo número de Péclet (10(2) menor que Pe menor que 10 e (6)), conforme ilustram os testes numéricos realizados. Coletando as experiências adquiridas com modelos lineares de Stokes e da advencção-difusão, nos foi possível propor, analisar o erro e testar dois novos métodos estabilizados para o problema de Navier-Stokes transiente. Construídos de maneira a herdar as boas características de estabilidade dos métodos propostos apresentam bom desempenho em escoamentos fortemente advectivos, bem como não necessitam atender a priori à condição de Baduska-Brezzi. Através de um algoritmo preditor/ multi-corretor de integração do termo inercial da equação de movimento, estes ,métodos foram capazes de de simular de maneira precisa escoamentos de interesse em Mecânica(400 menor que Re< menor que 500), captando escoamentos secundários, tais como recirculações de fluido. / [en] Stabilized methods for fluid problems are proposed and analysed with particular emphasis to the incompressible Navier-Stokes equations. We Begin in Chapter 2 introducing the balance equations of fluid Mechanics. Next. In Chapter 3, we discuss the numerical difficulties of the Galerkin method in fluids(in the contexto f the Stokes problem) and performance some succeful simulations of creeping flows, employing stabilized formulations. In Chapter 4, we propose a new finite element formulation for the energy equation, or more preciselly for the advective-diffusive model. Taking advantage of new design of the stability parameter T, which permits to add diffusion to advective and diffusive regions of the flow in a different way, we success to obtain a good performance of the new method in flows with very high Péclet numbers (10(2) lass than Pe lessa than 10(6)), as illustred at numerical testes performed. By collecting the Stokes and advective-diffusive experiences,it was possible to propose, analyse and test two new stabilized methods for the transient Navier-Stokes problem. These methods were built in a way to heritage the good characteristics showed by the stabilized methods introduced for the Stokes and adventive-diffusive models. The new methods propoposed have a good performance in high advective flows, besides there is no need to satisfy the Babuska-Brezzi condition. Employing a predictor/multi-corretor algorithm, we were able to simulate accruratly some useful flows(400 less than Re less than 500), such as fluid recirculations.

Finite element models have been successfully used to analyze adhesive bonds in actual structures, but this takes a considerable amount of time and a high computational cost. The objective of this study is to develop a simple and cost-effective finite element model for adhesively bonded joints which could be used in industry. Stress and durability analyses of crack patch geometries are possible applications of this finite element model. For example, the lifetime of aging aircraft can be economically extended by the application of patches bonded over the flaws located in the wings or the fuselage. Special two and three- dimensional adhesive elements have been developed for stress and displacement analyses in adhesively bonded joints. Both the 2-D and 3-D elements are used to model the whole adhesive system: adherends and adhesive layer. In the 2-D elements, adherends are represented by Bernoulli beam elements with axial deformation and the adhesive layer by plane stress or plane strain elements. The nodes of the plane stress-strain elements that lie in the adherend-adhesive interface are rigidly linked with the nodes of the beam elements. The 3-D elements consist of shell elements that represent the adherends and solid brick elements to model the adhesive. This technique results in smaller models with faster convergence than ordinary finite element models. The resulting mesh can represent arbitrary geometries of the adhesive layer and include cracks. Since large displacements are often observed in adhesively bonded joints, geometric nonlinearity is modeled. 2-D and 3-D stress analyses of single lap joints are presented. Important 3-D effects can be appreciated. Fracture mechanics parameters are computed for both cases. A stress analysis of a crack patch geometry is presented. A numerical simulation of the debonding of the patch is also included. / Ph. D.

This thesis studied the finite element modeling of plane wave acoustic phenomena in ducts. The study looked into finite element factors such as shape functions, mesh refinement, and element distortion. The study concluded that the higher order shape function eight-node quadrilateral element gave considerably better results than lower order shape function four-node quadrilateral element. The eight-node element converged much faster to the analytical solution than the four-node element. The average error, taking all the cases in consideration, for the four-node element was around 30 % for a mesh refinement of about 14 elements per wavelength at 100 Hz frequency. The eight-node element in the other hand had average absolute errors of less than 1% under the same conditions. This section also found that the eight-node element was substantially more resistant to solution deterioration due to element distortion than the four-node element. For example distorting the four-node element up to 60* degrees usually increased errors very rapidly to above 100 % errors. The eight-node element on the other hand usually produced errors of less than 5 % for the same level of distortion. The study showed that the type of boundary condition used had a significant effect on the solution accuracy. The study demonstrated that the effect of the natural boundary conditions was more global. Meeting this kind of boundary condition through mesh convergence produced accurate results throughout the duct. / Master of Science

No description available.

In this thesis a constitutive model is developed for the numerical prediction of UD composite material behaviour under impact loading. Impact induced loading usually results in three dimensional stress states which significantly influences the failure behaviour. The heterogeneous nature of composite materials, in particular, results in a complex failure behaviour which manifests itself in various failure modes. Predicting the onset and evolution of these failure modes requires the use of physically based three dimensional theories for the prediction of the onset of damage and subsequent damage evolution. Furthermore, the use of polymeric matrices in continuous fibre reinforced composites results in a distinct directional strain rate dependent material behaviour which needs to be incorporated in constitutive models for the numerical simulation of impact events. The developed constitutive model relies on the prediction of the onset of damage evolution by the use of physically based three dimensional stress based failure criteria. A special feature of the proposed model is the identification of potential fracture planes. Numerically efficient algorithms for finding such planes are developed thus enabling the implementation into an explicit FE environment which was prohibitive so far. Damage evolution is simulated by degrading the tractions which are acting on the failure mode dependent fracture planes. The damage evolution and consequent energy dissipation is thereby driven by physically based dissipation potentials which consider only stresses which contribute to damage growth. The well known mesh dependent energy dissipation in Continuum Damage Mechanics is reduced by the introduction of an element size dependent parameter into the constitutive equations. An experimental program is conducted to investigate the compressive behaviour of composites. The focus of the study is on the rate dependent failure behaviour. The experiments are designed such that the failure mechanisms can be studied at varying strain rates with identical boundary conditions. This allows for direct conclusions about the strain rate dependent material behaviour. Novel optical measurement techniques are applied across all investigated strain rates thus ensuring an improved observation of the failure modes. The proposed constitutive model is finally verified by modelling of three point beam bending experiments which were performed quasi-statically and at impact velocities. The experimental technique for beam bending at impact loading was therefore improved thus yielding significantly more accurate experimental data.