As análises de erros são conduzidas antes de qualquer projeto a ser desenvolvido. A necessidade do conhecimento do comportamento do erro numérico em malhas estruturadas e não-estruturadas surge com o aumento do uso destas malhas nos métodos de discretização. Desta forma, o objetivo deste trabalho foi criar uma metodologia para analisar os erros de discretização gerados através do truncamento na Série de Taylor, aplicados às equações de Poisson e de Advecção-Difusão estacionárias uni e bidimensionais, utilizando-se o Método de Volumes Finitos em malhas do tipo Voronoi. A escolha dessas equações se dá devido a sua grande utilização em testes de novos modelos matemáticos e função de interpolação. Foram usados os esquemas Central Difference Scheme (CDS) e Upwind Difference Scheme(UDS) nos termos advectivos. Verificou-se a influência do tipo de condição de contorno e a posição do ponto gerador do volume na solução numérica. Os resultados analíticos foram confrontados com resultados experimentais para dois tipos de malhas de Voronoi, uma malha cartesiana e outra triangular comprovando a influência da forma do volume finito na solução numérica obtida. Foi percebido no estudo que a discretização usando o esquema CDS tem erros menores do que a discretização usando o esquema UDS conforme literatura. Também se percebe a diferença nos erros em volumes vizinhos nas malhas triangulares o que faz com que não se tenha uma uniformidade nos gráficos dos erros estudados. Percebeu-se que as malhas cartesianas com nó no centróide do volume tem menor erro de discretização do que malhas triangulares. Mas o uso deste tipo de malha depende da geometria do problema estudado / The analyses of errors are conducted before any project to be developed. The necessity of studying the behavior of the numerical error on structured and unstructured grids comes up with the increasing use of these methods of discretization meshes. Thus, the objective was to create a methodology to analyze the errors generated by discretization of the truncation in the Taylor series, applied to the equations of Poisson and Advection-Diffusion stationary and uni and bi-dimensional, using the Finite Volume Method on Voronoi mesh. The choice of these equations is due to its wide use in testing new mathematical models and interpolation function. The schemes used were the Central Difference Scheme (CDS) and the Upwind Difference Scheme (UDS) in the advective terms. There was the influence of boundary condition and position of the generator in the numerical solution of the volume. The analytical results were compared with experimental results for two types of Voronoi meshes, a Cartesian mesh and a triangular shape showing the influence of finite volume in the numerical solution obtained. It was perceived that the discretization in the study using the CDS scheme has smaller errors than the discretization scheme using the UDS as literature. Also notice the difference in the errors in neighboring volumes in triangular meshes which means that there has been no uniformity in the graphs of errors studied. It was noticed that the Cartesian meshes with node at the centroid of the volume is smaller than discretization error triangular meshes. But using this type of meshes depends on the geometry of the problem studied

Moldagem por injeção é um dos mais importantes processos industriais para produção de produtos plásticos finos. Esse processo é dividido essencialmente em quatro estágios: plastificação, preenchimento, empacotamento e resfriamento. O escoamento de um fluido caracterizado por alta viscosidade em uma cavidade estreita é um problema tipicamente encontrado em processos de moldagem por injeção.Neste caso, o escoamento pode ser descrito por uma formulação conhecida como aproximação de Hele-Shaw. Tal formulação pode ser derivada das equações de conservação tridimensionais usando um número de suposições a respeito do polímero injetado e da geometria da cavidade do molde, juntamente com a integração e o acoplamento das equações da conservação da quantidade de movimento e da continuidade. Essa formulação, referindo às limitações da geometria do molde como sendo canais estreitos e quase sem curvatura, é comumente denominada formulação 2 1/2D. Neste trabalho, é apresentada uma técnica para a simulação da fase de preenchimento de um processo de moldagem por injeção, usando essa formulação 2 1/2D, com um método de volumes finitos e malhas não estruturadas. O modelo de Cross modificado com dependência da temperatura de Arrhenius é empregado para descrever a viscosidade do polímero fundido. O campo de distribuição de temperatura é tridimensional e é resolvido usando um esquema semi-Lagrangeano baseado em volumes finitos. As malhas não estruturadas utilizadas são geradas por triangulação de Delaunay e o método numérico implementado usa a estrutura de dados topológica SHE - Singular Handle Edge, que é capaz de lidar com condições de contorno e singularidades, aspectos comumente encontrados em simulações numéricas de escoamento de fluidos. / Injection molding is one of the most important industrial processes for the manufacturing of thin plastic products. This process can be divided into four stages: plastic melting, filling, packing and cooling phases. The flow of a fluid characterized by high viscosity in a narrow gap is a problem typically found in injection molding processes. In this case, the flow can be described by a formulation known as Hele-Shaw approach. Such formulation can be btained from the three-dimensional conservation equation using a number of assumptions regarding the injected polymer and the geometry of the mold, together with the integration and the coupling of the momentum and continuity equations. This approach, referring to limitations of the mould geometry to narrow, weakly curved channels, is usually called 2 1/2D approach. In this work a technique for the simulation of the filling stage of the injection molding process, using this 2 1/2D approach, with a finite volume method and unstructured meshes, is presented. The modified-Cross model with Arrhenius temperature dependence is employed to describe the viscosity of the melt. The temperature field is 3D and it is solved using a semi-Lagrangian scheme based on the finite volume method. The employed unstructured meshes are generated by Delaunay triangulation and the implemented numerical method uses the topological data structure SHE - Singular Handle Edge, capable to deal with boundary conditions and singularities, aspects commonly found in numerical simulation of fluid flow.

Thickened tailings production and disposal continue to grow in importance in the mining industry. In particular, the transport of oil sands tailings is of interest in this study. These tailings must be in a homogeneous state (non-segregating) during pipeline flow and subsequent discharge. Tailings are often transported in an open channel or flume. Slurries containing both clay and coarse sand particles typically exhibit non-Newtonian rheological behaviour. The prediction of the flow behaviour of these slurries is complicated by the limited research activity in this area. As a result, the underlying mechanisms of solids transport in these slurries are not well understood. To address this deficiency, experimental studies were conducted with kaolin clay slurries containing coarse sand in an open circular channel.<p> A numerical model has been developed to predict the behaviour of coarse solid particles in laminar, open channel, non-Newtonian flows. The model involves the simultaneous solution of the Navier-Stokes equations and a scalar concentration equation describing the behaviour of coarse particles within the flow. The model uses the theory of shear-induced particle diffusion (Phillips et al., 1992) to provide a number of relationships to describe the diffusive flux of coarse particles within laminar flows. A sedimentation flux has been developed and incorporated into the Phillips et al. (1992) model to account for gravitational flux of particles within the flow. Previous researchers (Gillies et al., 1999) have shown that this is a significant mechanism of particle migration.<p> The momentum and concentration partial differential equations have been solved numerically by applying the finite volume method. The differential equations are non-linear, stiff and tightly coupled which requires a novel means of analysis. Specific no-flux, no-slip and no-shear boundary conditions have been applied to the channel walls and free surface to produce simulated velocity and concentration distributions. The results show that the model is capable of predicting coarse particle settling in laminar, non-Newtonian, open channel flows. The results of the numerical simulations have been compared to the experimental results obtained in this study, as well as the experimental results of previous studies in the literature.

This thesis studies advection modeling for heat transfer and fluid flow problems using a new Non--Inverted Skew Upwind Scheme (called NISUS). Variants of the new scheme are formulated and developed with 8-noded hexahedral elements using the Finite Element Method (FEM)and rectangular elements based on a Finite Volume Method (FVM). A new method of mass weighting to predict convective fluxes of each scalar from the nodal point values is developed. Due to an explicit representation in terms of nodal variables, local inversion of the upwind coefficient matrix is not needed. Also, this thesis evaluates two variants of the new scheme (i.e., 3-node / 3-point and 4-node / 8-point formulations) within a 3--D FEM and a third variant within a 2--D FVM. The 3--D FEM variants are applied to a variety of test problems involving the transport of a scalar variable, while the 2--D FVM variant is applied to fluid flow problems including natural convection in an enclosure and micro--channel flow simulations. The promising performance of NISUS, as compared with exact and previous solutions, is demonstrated both in terms of accuracy and stability. Furthermore, a new data storage format called Compressed Banded Data (CBD) is developed for sparse banded matrices generated by the control volume finite element method (CVFEM). The platform of the new CBD structure permits dynamic switching between various solvers, without any procedural change in the implementation of existing simulation software. The performance of different Krylov techniques with an ILU(0) preconditioner is observed and compared in three test problems with a direct solver. / October 2006

Electric breakdown in power transformers is preceded by pre-breakdown events such as streamers. The understanding of these phenomena is important in order to optimize liquid insulation systems. Earlier works have derived a model that describes streamers in transformer oil and utilized a finite element method to produce numerical solutions. This research investigates the consequences of changing the numerical method to a finite volume-based solver implemented in OpenFOAM. Using a standardized needle-sphere geometry, a number of oil and voltage combinations were simulated, and the results are for the most part similar to those produced by the previous method. In cases with differing results the change is attributed to the more stable numerical performance of the OpenFOAM solver. A proof of concept for the extension of the simulation from a two-dimensional axial symmetry to three dimensions is also presented. / Elektriska genomslag i högspänningstransformatorer föregås av bildandet av elektriskt ledande kanaler som kallas streamers. En god förståelse av detta fenomen är viktigt vid konstruktionen av oljebaserad elektrisk isolation. Tidigare forskning i ämnet har tagit fram en modell för fortplantningen av streamers. Denna modell har sedan lösts numeriskt av ett beräkningsverktyg baserat på finita elementmetoden. I denna uppsats undersöks konsekvenserna av att byta metod till finita volymsmetoden genom att implementera en lösare i OpenFOAM. En standardiserad nål-sfär-geometri har ställts upp och ett flertal kombinationer av oljor och spänningar har simulerats. De flesta resultaten visar god överensstämmande med tidigare forskning medan resultat som avviker har tillskrivits de goda numeriska egenskaperna hos OpenFOAM-lösaren. En ny typ av simulering har även genomförts där simulationen utökas från en tvådimensionell axisymmetrisk geometri til tre dimensioner.

In this thesis we address several analytical and numerical problems related with the general relativistic study of black hole space-times and boson stars. We have developed a new centered finite volume method based on the flux splitting approach. The techniques for dealing with the singularity, steep gradients and apparent horizon location, are studied in the context of a single Schwarzschild black hole, in both spherically symmetric and full 3D simulations. We present an extended study of gauge instabilities related with a class of singularity avoiding slicing conditions and show that, contrary to previous claims, these instabilities are not generic for evolved gauge conditions. We developed an alternative to the current space coordinate conditions, based on a generalized Almost Killing Equation. We performed a general relativistic study regarding the long term stability of Mixed-State Boson Stars configurations and showed that they are suitable candidates for dark matter models. / En esta tesis abordamos varios problemas analíticos y numéricos relacionados con el estudio de agujeros negros relativistas y modelos de materia oscura. Hemos desarrollado un nuevo método de volúmenes finitos centrados basado en el enfoque de la división de flujo. Discutimos las técnicas para tratar con la singularidad, los gradientes abruptos y la localización del horizonte aparente en el contexto de un solo agujero negro de Schwarzschild, en simulaciones tanto con simetría esférica como completamente tridimensionales. Hemos extendido el estudio de una familia de condiciones de foliaciones evitadoras de singularidad y mostrado que ciertas inestabilidades no son genéricas para condiciones de gauge dinámicas. Desarrollamos una alternativa a las prescripciones actuales basada en una Almost Killing Equation generalizada. Hemos realizado también un estudio con respecto a la estabilidad a largo plazo de configuraciones de Mixed-State Boson Stars, el cual sugiere que estas podrían ser candidatas apropiadas para modelos de materia oscura.

The linear elastodynamic two-dimensional anti-plane stress problem, where deformations occur in only one direction is considered for one sided non-planar faults. Fault dynamics are modeled using purely velocity dependent friction laws, and applied on boundaries with complex geometry. Summation-by-parts operators and energy estimates are used to couple a high-order finite difference method with an unstructured finite volume method. The unstructured finite volume method is used near the fault and the high-order finite difference method further away from the fault where no complex geometry is present. Boundary conditions are imposed weakly on characteristic form using the simultaneous approximation term technique, allowing explicit time integration to be used. Numerical computations are performed to verify the accuracy and time stability, of the method.

This master thesis was provided by Scania AB. The objective of this thesis was to modify an application in the free Computational Fluid Dynamics software OpenFOAM to be able to handle spray and wall film modeling of a Urea Water Solution together with Conjugate Heat Transfer. The basic purpose is to widen the knowledge of the vaporization process of a Urea Water Solution in the exhaust gas after treatment system for a diesel engine by using OpenFOAM. First, urea has been modeled as a very viscous liquid at low temperature to mimic the solidication process of urea. Second, the development of the new application has been done. At last, test simulations of a simple test case are performed with the new application. The results are then compared with simplied hand calculations to verify a correct behavior of certain exposed source terms. The new application is working properly for the test case but to ensure the reliability, the results need to be compared with another Computational Fluid Dynamics software or more preferable, real experiments. For more advanced geometries, the continued development presented last in this thesis is highly recommended to follow.

Thickened tailings production and disposal continue to grow in importance in the mining industry. In particular, the transport of oil sands tailings is of interest in this study. These tailings must be in a homogeneous state (non-segregating) during pipeline flow and subsequent discharge. Tailings are often transported in an open channel or flume. Slurries containing both clay and coarse sand particles typically exhibit non-Newtonian rheological behaviour. The prediction of the flow behaviour of these slurries is complicated by the limited research activity in this area. As a result, the underlying mechanisms of solids transport in these slurries are not well understood. To address this deficiency, experimental studies were conducted with kaolin clay slurries containing coarse sand in an open circular channel.<p> A numerical model has been developed to predict the behaviour of coarse solid particles in laminar, open channel, non-Newtonian flows. The model involves the simultaneous solution of the Navier-Stokes equations and a scalar concentration equation describing the behaviour of coarse particles within the flow. The model uses the theory of shear-induced particle diffusion (Phillips et al., 1992) to provide a number of relationships to describe the diffusive flux of coarse particles within laminar flows. A sedimentation flux has been developed and incorporated into the Phillips et al. (1992) model to account for gravitational flux of particles within the flow. Previous researchers (Gillies et al., 1999) have shown that this is a significant mechanism of particle migration.<p> The momentum and concentration partial differential equations have been solved numerically by applying the finite volume method. The differential equations are non-linear, stiff and tightly coupled which requires a novel means of analysis. Specific no-flux, no-slip and no-shear boundary conditions have been applied to the channel walls and free surface to produce simulated velocity and concentration distributions. The results show that the model is capable of predicting coarse particle settling in laminar, non-Newtonian, open channel flows. The results of the numerical simulations have been compared to the experimental results obtained in this study, as well as the experimental results of previous studies in the literature.

A numerical method (SIMPLE DIRK Method) for unsteady incompressible viscous flow simulation is presented. The proposed method can be used to achieve arbitrarily high order of accuracy in time-discretization which is otherwise limited to second order in majority of the currently used simulation techniques. A special class of implicit Runge-Kutta methods is used for time discretization in conjunction with finite volume based SIMPLE algorithm. The algorithm was tested by solving for velocity field in a lid-driven square cavity. In the test case calculations, power law scheme was used in spatial discretization and time discretization was performed using a second-order implicit Runge-Kutta method. Time evolution of velocity profile along the cavity centerline was obtained from the proposed method and compared with that obtained from a commercial computational fluid dynamics software program, FLUENT 6.2.16. Also, steady state solution from the present method was compared with the numerical solution of Ghia, Ghia, and Shin and that of Erturk, Corke, and Goökçöl. Good agreement of the solution of the proposed method with the solutions of FLUENT; Ghia, Ghia, and Shin; and Erturk, Corke, and Goökçöl establishes the feasibility of the proposed method.