Simulations were performed to evaluate solution adaptive meshing strategies for flows with vortices whose axes of rotation are parallel to the bulk fluid motion. Two configurations were investigated: a wing in a wind tunnel and a missile spinning at 30Hz and 60Hz at 0◦ angle of attack with canards deflected 15◦. Feature-based descriptors were used to identify regions of the flow near vortices that are candidate regions for adaptive meshing. Several different adaptive meshing techniques were evaluated. These techniques include refinement around the vortex core, refinement near the vortex extent surface, refinement inside the extent surface, refinement inside and near the extent surface, and mesh regeneration using the vortex extent surface as an embedded surface. Results for the wing case, compared to experimental data, indicate that it is necessary to refine the region within and near the vortex extent surface to accurately recreate physical characteristics and achieve an acceptable solution.

No description available.

L’objectif de cette étude est de simuler l’érosion d’un sol cohésif sous l’effet d’un écoulement incompressible. Le modèle élaboré décrit une vitesse d’érosion interfaciale qui dépend de la contrainte de cisaillement de l’écoulement. La modélisation numérique proposée est une approche eulérienne, où une méthode de pénalisation de domaines est utilisée pour résoudre les équations de Navier-Stokes autour d’un obstacle. L’interface eau/sol est décrite par une fonction Level Set couplée à une loi d’érosion à seuil.L’approximation numérique est basée sur un schéma DDFV (Discrete Duality Finite Volume) autorisant des raffinements locaux sur maillages non-conformes et non-structurés. L’approche par pénalisation a mis en évidence une couche limite d'inconsistance à l'interface fluide/solide lors du calcul de la contrainte de cisaillement. Deux approches sont proposées pour estimer précisément la contrainte de ce problème à frontière libre. La pertinence du modèle à prédire l’érosion interfaciale du sol est confirmée par la présentation de plusieurs résultats de simulation, qui offrent une meilleure évaluation et compréhension des phénomènes d'érosion / This study focuses on the numerical modelling of the interfacial erosion occurring at a cohesive soil undergoing an incompressible flow process. The model assumes that the erosion velocity is driven by a fluid shear stress at the water/soil interface. The numerical modelling is based on the eulerian approach: a penalization procedure is used to compute Navier-Stokes equations around soil obstacle, with a fictitious domain method, in order to avoid body- fitted unstructured meshes. The water/soil interface’s evolution is described by a Level Set function coupled to a threshold erosion law.Because we use adaptive mesh refinement, we develop a Discrete Duality Finite Volume scheme (DDFV), which allows non-conforming and non-structured meshes. The penalization method, used to take into account a free velocity in the soil with non-body-fitted mesh, introduces an inaccurate shear stress at the interface. We propose two approaches to compute accurately the erosion velocity of this free boundary problem. The ability of the model to predict the interfacial erosion of soils is confirmed by presenting several simulations that provide better evaluation and comprehension of erosion phenomena.

No presente trabalho, utilizou-se o Método da Fronteira Imersa, o qual utiliza dois tipos de malhas computacionais: euleriana (utilizada para o fluido) e lagrangiana (utilizada para representar a interface de separação de dois fluidos). O software livre GMSH foi utilizado para representar um sólido por meio da sua superfície externa e também para gerar uma malha triangular, bidimensional e não estruturada para discretizar essa superfície. Essa superfície foi utilizada como condição inicial para a malha lagrangiana (fronteira imersa). Os dados da malha lagrangiana são armazenados em uma estrutura de dados chamada Halfedge, a qual é largamente utilizada em Computação Gráfica para armazenar superfícies fechadas e orientáveis. Uma vez que a malha lagrangiana esteja armazenada nesta estrutura de dados, passa-se a estudar uma hipotética interação dinâmica entre a fronteira imersa e o escoamento do fluido. Esta interação é estudada apenas em um sentido, considera-se apenas a condição de não deslizamento, isto é, a fronteira imersa acompanhará passivamente um campo de velocidades pré-estabelecido (imposto), sem exercer qualquer força ou influência sobre ele. Foi utilizado um campo de distância local com sinal (função indicadora de fluidos) para identificar o interior e o exterior da superfície que representa a interface entre os fluidos. Este campo de distância é atualizado a cada passo no tempo utilizando idéias de Geometria Computacional, o que tornou o custo computacional para calcular esse campo otimal independente da complexidade geométrica da interface. Esta metodologia mostrou-se robusta e produz uma definição nítida das distintas fases dos fluidos em todos os passos no tempo. Para acompanhar e visualizar de forma mais precisa o comportamento dos fluidos na vizinhança da superfície que representa a interface de separação dos fluido, foi utilizado um algoritmo chamado de Refinamento Adaptativo de Malhas para fazer um refinamento dinâmico da malha euleriana na vizinhança da malha lagrangiana. / The scientific motivation of the present work is the mathematical modeling and the computational simulation of multiphase flows. Specifically, the equations of a two-phase flow are written by combining the Immersed Boundary Method with a suitable fluid indicator function. It is assumed that the fluid equations are discretized on an Eulerian mesh covering completely the flow domain and that the interface between the fluid phases is discretized by a non-structured Lagrangian mesh formed by triangles. In this context, employing tools commonly found in Computational Geometry, the computation of the fluid indicator function is efficiently performed on a block-structured Eulerian mesh bearing dynamical refinement patches. Formed by a set of triangles, the Lagrangian mesh, which is initally generated employing the free software GMSH, is stored in a Halfedge data structure, a data structure which is widely used in Computer Graphics to represent bounded, orientable closed surfaces. Once the Lagrangian mesh has been generated, next, one deals with the hipothetical situation of dealing with the one-way dynamical interaction between the immersed boundary and the fluid flow, that is, considering the non-slip condition, only the action of the flow on the interface is studied. No forces arising on the interface affects the flow, the interface passively being advect with the flow under a prescribed, imposed velocity field. In particular, the Navier-Stokes equations are not solved. The fluid indicator function is given by a signed distance function in a vicinity of the immersed boundary. It is employed to identify interior/exterior points with respect to the bounded, closed region which is assumed to contain one of the fluid phases in its interior. The signed distance is update every time step employing Computational Geometry methods with optimal cost. Several examples in three dimensions, showing the efficiency and efficacy in the computation of the fluid indicator function, are given which employ the dynamical adaptive properties of the Eurlerian mesh for a moving interface.

No presente trabalho, utilizou-se o Método da Fronteira Imersa, o qual utiliza dois tipos de malhas computacionais: euleriana (utilizada para o fluido) e lagrangiana (utilizada para representar a interface de separação de dois fluidos). O software livre GMSH foi utilizado para representar um sólido por meio da sua superfície externa e também para gerar uma malha triangular, bidimensional e não estruturada para discretizar essa superfície. Essa superfície foi utilizada como condição inicial para a malha lagrangiana (fronteira imersa). Os dados da malha lagrangiana são armazenados em uma estrutura de dados chamada Halfedge, a qual é largamente utilizada em Computação Gráfica para armazenar superfícies fechadas e orientáveis. Uma vez que a malha lagrangiana esteja armazenada nesta estrutura de dados, passa-se a estudar uma hipotética interação dinâmica entre a fronteira imersa e o escoamento do fluido. Esta interação é estudada apenas em um sentido, considera-se apenas a condição de não deslizamento, isto é, a fronteira imersa acompanhará passivamente um campo de velocidades pré-estabelecido (imposto), sem exercer qualquer força ou influência sobre ele. Foi utilizado um campo de distância local com sinal (função indicadora de fluidos) para identificar o interior e o exterior da superfície que representa a interface entre os fluidos. Este campo de distância é atualizado a cada passo no tempo utilizando idéias de Geometria Computacional, o que tornou o custo computacional para calcular esse campo otimal independente da complexidade geométrica da interface. Esta metodologia mostrou-se robusta e produz uma definição nítida das distintas fases dos fluidos em todos os passos no tempo. Para acompanhar e visualizar de forma mais precisa o comportamento dos fluidos na vizinhança da superfície que representa a interface de separação dos fluido, foi utilizado um algoritmo chamado de Refinamento Adaptativo de Malhas para fazer um refinamento dinâmico da malha euleriana na vizinhança da malha lagrangiana. / The scientific motivation of the present work is the mathematical modeling and the computational simulation of multiphase flows. Specifically, the equations of a two-phase flow are written by combining the Immersed Boundary Method with a suitable fluid indicator function. It is assumed that the fluid equations are discretized on an Eulerian mesh covering completely the flow domain and that the interface between the fluid phases is discretized by a non-structured Lagrangian mesh formed by triangles. In this context, employing tools commonly found in Computational Geometry, the computation of the fluid indicator function is efficiently performed on a block-structured Eulerian mesh bearing dynamical refinement patches. Formed by a set of triangles, the Lagrangian mesh, which is initally generated employing the free software GMSH, is stored in a Halfedge data structure, a data structure which is widely used in Computer Graphics to represent bounded, orientable closed surfaces. Once the Lagrangian mesh has been generated, next, one deals with the hipothetical situation of dealing with the one-way dynamical interaction between the immersed boundary and the fluid flow, that is, considering the non-slip condition, only the action of the flow on the interface is studied. No forces arising on the interface affects the flow, the interface passively being advect with the flow under a prescribed, imposed velocity field. In particular, the Navier-Stokes equations are not solved. The fluid indicator function is given by a signed distance function in a vicinity of the immersed boundary. It is employed to identify interior/exterior points with respect to the bounded, closed region which is assumed to contain one of the fluid phases in its interior. The signed distance is update every time step employing Computational Geometry methods with optimal cost. Several examples in three dimensions, showing the efficiency and efficacy in the computation of the fluid indicator function, are given which employ the dynamical adaptive properties of the Eurlerian mesh for a moving interface.

Este trabalho apresenta um estudo sobre a solução numérica da equação geral de advecção-difusão usando uma metodologia numérica conservativa. Para a discretização espacial, é usado o Método de Volumes Finitos devido à natureza conservativa da equação em questão. O método é configurado de modo a ter suas variáveis centradas em centro de célula e, para as variáveis, como a velocidade, centradas nas faces um método de interpolação de segunda ordem é utilizado para um ajuste numérico ao centro. Embora a implementação computacional tenha sido feita de forma paramétrica de maneira a acomodar outros esquemas numéricos, a discretização temporal dá ênfase ao Método de Crank-Nicolson. Tal método numérico, sendo ele implícito, dá origem a um sistema linear de equações que, aqui, é resolvido empregando-se o Método Multigrid-Multinível. A corretude do código implementado é verificada a partir de testes por soluções manufaturadas, de modo a checar se a ordem de convergência prevista em teoria é alcançada pelos métodos numéricos. Um jato laminar é simulado, com o acoplamento entre a equação de Navier-Stokes e a equação geral de advecção-difusão, em um domínio computacional tridimensional. O jato é uma forma de vericar se o algoritmo de geração de malhas adaptativas funciona corretamente. O módulo produzido neste trabalho é baseado no código computacional AMR3D-P desenvolvido pelos grupos de pesquisa do IME-USP e o MFLab/FEMEC-UFU (Laboratório de Dinâmica de Fluidos da Universidade Federal de Uberlândia). A linguagem FORTRAN é utilizada para o desenvolvimento da metodologia numérica e as simulações foram executadas nos computadores do LabMAP(Laboratório da Matemática Aplicada do IME-USP) e do MFLab/FEMEC-UFU. / This work presents a study about the numerical solution of variable coecients advectiondi usion equation, or simply, general advection-diusion equation using a conservative numerical methodology. The Finite Volume Method is choosen as discretisation of the spatial domain because the conservative nature of the focused equation. This method is set up to have the scalar variable in a cell centered scheme and the vector quantities, such velocity, are face centered and they need a second order interpolation to get adjusted to the cell center. The computational code is parametric, in which, any implicit temporal discretisation can be choosen, but the emphasis relies on Crank-Nicolson method, a well-known second order method. The implicit nature of aforementioned method gives a linear system of equations which is solved here by the Multilevel-Multigrid method. The correctness of the computational code is checked by manufactured solution method used to inspect if the theoretical order of convergence is attained by the numerical methods. A laminar jet is simulated, coupling the Navier-Stokes equation and the general advection-diusion equation in a 3D computational domain. The jet is a good way to check the corectness of adaptative mesh renement algorithm. The module designed here is based in a previous implemented code AMR3D-P designed by IME-USP and MFLab/FEMEC-UFU (Fluid Dynamics Laboratory, Federal University of Uberlândia). The programming language used is FORTRAN and the simulations were run in LabMAP(Applied Mathematics Laboratoy at IME-USP) and MFLab/FEMEC-UFU computers.

Na obtenção de aproximações numéricas para Equações Diferenciais Parciais Elípticas utilizando o Método dos Elementos Finitos (MEF) alguns problemas apresentam valores maiores para o erro somente em algumas determinadas regiões do domínio como, por exemplo, regiões onde existam singularidades na solução contínua do problema. Uma possível alternativa para reduzir o erro cometido nestas regiões é aumentar o número de elementos nos trechos onde o erro cometido foi considerado grande. A questão principal é como identificar essas regiões, dado que a solução do problema contínuo é desconhecida. Neste trabalho iremos apresentar a chamada estimativa residual, que fornece um estimador do erro cometido na aproximação utilizando apenas os valores conhecidos dos contornos e a aproximação obtida sobre uma dada partição de elementos. Vamos discutir a relação entre a estimativa residual e o erro cometido na aproximação, além de utilizar as estimativas na construção de um algoritmo adaptativo para as malhas em estudo. Utilizando o software FreeFem++ serão obtidas aproximações para a Equação de Poisson e para o sistema de equações associado à Elasticidade Linear e por meio do estimador residual será analisado o erro cometido nas aproximações e a necessidade do refinamento adaptativo das malhas. / In obtaining numerical approximations for solutions to Elliptic Partial Differential Equations using the Finite Element Method (FEM) one sees that some problems have higher values for the error only in certain domain regions such as, for example, regions where the solution of the continous problem is singular. A possible alternative to reduce the error in these regions is to increase the number of elements in the partions where the error was considered large. The main issue is how to identify these regions, since the solution of the continuous problem is unknown. In this work we present the so-called residual estimate, which provides an error estimation approach which uses only the known values on the contours and the obtained approximation on a given discretization. We will discuss the relationship between the residual estimate and the error, and how to use the estimate for adaptively refining the mesh. Solutions for the Poisson equation and the Linear elasticity system of equations, and the residual estimates for the analysis of mesh refinement will be computed using the FreeFem++ software.

Many real-life applications such as the shape optimization of technological devices, the identification of parameters in environmental processes and flow control problems lead to optimization problems governed by systems of convection diusion partial dierential equations (PDEs). When convection dominates diusion, the solutions of these PDEs typically exhibit layers on small regions where the solution has large gradients. Hence, it requires special numerical techniques, which take into account the structure of the convection. The integration of discretization and optimization is important for the overall eciency of the solution process. Discontinuous Galerkin (DG) methods became recently as an alternative to the finite dierence, finite volume and continuous finite element methods for solving wave dominated problems like convection diusion equations since they possess higher accuracy. This thesis will focus on analysis and application of DG methods for linear-quadratic convection dominated optimal control problems. Because of the inconsistencies of the standard stabilized methods such as streamline upwind Petrov Galerkin (SUPG) on convection diusion optimal control problems, the discretize-then-optimize and the optimize-then-discretize do not commute. However, the upwind symmetric interior penalty Galerkin (SIPG) method leads to the same discrete optimality systems. The other DG methods such as nonsymmetric interior penalty Galerkin (NIPG) and incomplete interior penalty Galerkin (IIPG) method also yield the same discrete optimality systems when penalization constant is taken large enough. We will study a posteriori error estimates of the upwind SIPG method for the distributed unconstrained and control constrained optimal control problems. In convection dominated optimal control problems with boundary and/or interior layers, the oscillations are propagated downwind and upwind direction in the interior domain, due the opposite sign of convection terms in state and adjoint equations. Hence, we will use residual based a posteriori error estimators to reduce these oscillations around the boundary and/or interior layers. Finally, theoretical analysis will be confirmed by several numerical examples with and without control constraints

The paper is concerned with the finite element resolution of layers appearing in singularly perturbed problems. A special anisotropic grid of Shishkin type is constructed for reaction diffusion problems. Estimates of the finite element error in the energy norm are derived for two methods, namely the standard Galerkin method and a stabilized Galerkin method. The estimates are uniformly valid with respect to the (small) diffusion parameter. One ingredient is a pointwise description of derivatives of the continuous solution. A numerical example supports the result. Another key ingredient for the error analysis is a refined estimate for (higher) derivatives of the interpolation error. The assumptions on admissible anisotropic finite elements are formulated in terms of geometrical conditions for triangles and tetrahedra. The application of these estimates is not restricted to the special problem considered in this paper.

The paper presents a method for generating curved surfaces of 3D finite element meshes by mesh refinement starting with a very coarse grid. This is useful for parallel implementations where the finest meshes should be computed and not read from large files. The paper deals with simple geometries as sphere, cylinder, cone. But the method may be extended to more complicated geometries. (with 45 figures)