Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Os escoamentos altamente convectivos representam um desafio na simulação pelo método de elementos finitos. Com a solução de elementos finitos de Galerkin para escoamentos incompressíveis, a matriz associada ao termo convectivo é não simétrica, e portanto, a propiedade de aproximação ótima é perdida. Na prática as soluções apresentam oscilações espúrias. Muitos métodos foram desenvolvidos com o fim de resolver esse problema. Neste trabalho apresentamos um método semi- Lagrangeano, o qual é implicitamente um método do tipo upwind, que portanto resolve o problema anterior, e comparamos o desempenho do método na solução das equações de convecção-difusão e Navier-Stokes incompressível com o Streamline Upwind Petrov Galerkin (SUPG), um método estabilizador de reconhecido desempenho. No SUPG, as funções de forma e de teste são tomadas em espaços diferentes, criando um efeito tal que as oscilações espúrias são drasticamente atenuadas. O método semi-Lagrangeano é um método de fator de integração, no qual o fator é um operador de convecção que se desloca para um sistema de coordenadas móveis no fluido, mas restabelece o sistema de coordenadas Lagrangeanas depois de cada passo de tempo. Isto prevê estabilidade e a possibilidade de utilizar passos de tempo maiores.Existem muitos trabalhos na literatura analisando métodos estabilizadores, mas não assim com o método semi-Lagrangeano, o que representa a contribuição principal deste trabalho: reconhecer as virtudes e as fraquezas do método semi-Lagrangeano em escoamentos dominados pelo fenômeno de convecção. / Convection dominated flows represent a challenge for finite element method simulation. Many methods have been developed to address this problem. In this work we compare the performance of two methods in the solution of the convectiondiffusion and Navier-Stokes equations on environmental flow problems: the Streamline Upwind Petrov Galerkin (SUPG) and the semi-Lagrangian method. In Galerkin finite element methods for fluid flows, the matrix associated with the convective term is non-symmetric, and as a result, the best approximation property is lost. In practice, solutions are often corrupted by espurious oscillations. In this work, we present a semi- Lagrangian method, which is implicitly an upwind method, therefore solving the spurious oscillations problem, and a comparison between this semi-Lagrangian method and the Streamline Upwind Petrov Galerkin (SUPG), an stabilizing method of recognized performance. The SUPG method takes the interpolation and the weighting functions in different spaces, creating an effect so that the spurious oscillations are drastically attenuated. The semi-Lagrangean method is a integration factor method, in which the factor is an operator that shifts to a coordinate system that moves with the fluid, but it resets the Lagrangian coordinate system after each time step. This provides stability and the possibility to take bigger time steps. There are many works in the literature analyzing stabilized methods, but they do not analyze the semi-Lagrangian method, which represents the main contribution of this work: to recognize the strengths and weaknesses of the semi-Lagrangian method in convection dominated flows.

Least Square Kinetic Upwind Method (LSKUM) belongs to the class of mesh-less method that solves compressible Euler equations of gas dynamics. LSKUM is kinetic theory based upwind scheme that operates on any cloud of points. Euler equations are derived from Boltzmann equation (of kinetic theory of gases) after taking suitable moments. The basic update scheme is formulated at Boltzmann level and mapped to Euler level by suitable moments. Mesh-less solvers need only cloud of points to solve the governing equations. For a complex configuration, with such a solver, one can generate a separate cloud of points around each component, which adequately resolves the geometric features, and then combine all the individual clouds to get one set of points on which the solver directly operates. An obvious advantage of this approach is that any incremental changes in geometry will require only regeneration of the small cloud of points where changes have occurred. Additionally blanking and de-blanking strategy along with overlay point cloud can be adapted in some applications like store separation to avoid regeneration of points. Naturally, the mesh-less solvers have advantage in tackling complex geometries and moving components over solvers that need grids. Conventionally, higher order accuracy for space derivative term is achieved by two step defect correction formula which is computationally expensive. The present solver uses low dissipation single step modified CIR (MCIR) scheme which is similar to first order LSKUM formulation and provides spatial accuracy closer to second order. The maximum time step taken to march solution in time is limited by stability criteria in case of explicit time integration procedure. Because of this, explicit scheme takes a large number of iterations to achieve convergence. The popular explicit time integration schemes like four stages Runge-Kutta (RK4) are slow in convergence due to this reason. The above problem can be overcome by using the implicit time integration procedure. The implicit schemes are unconditionally stable i.e. very large time steps can be used to accelerate the convergence. Also it offers superior robustness. The implicit Lower-Upper Symmetric Gauss-Seidel (LU-SGS) scheme is very attractive due to its low numerical complexity, moderate memory requirement and unconditional stability for linear wave equation. Also this scheme is more efficient than explicit counterparts and can be implemented easily on parallel computers. It is based on the factorization of the implicit operator into three parts namely lower triangular matrix, upper triangular matrix and diagonal terms. The use of LU-SGS results in a matrix free implicit framework which is very economical as against other expensive procedures which necessarily involve matrix inversion. With implementation of the implicit LU-SGS scheme larger time steps can be used which in turn will reduce the computational time substantially. LU-SGS has been used widely for many Finite Volume Method based solvers. The split flux Jacobian formulation as proposed by Jameson is most widely used to make implicit procedure diagonally dominant. But this procedure when applied to mesh-less solvers leads to block diagonal matrix which again requires expensive inversion. In the present work LU-SGS procedure is adopted for mesh-less approach to retain diagonal dominancy and implemented in 2-D and 3-D solvers in matrix free framework. In order to assess the efficacy of the implicit procedure, both explicit and implicit 2-D solvers are tested on NACA 0012 airfoil for various flow conditions in subsonic and transonic regime. To study the performance of the solvers on different point distributions two types of the cloud of points, one unstructured distribution (4074 points) and another structured distribution (9600 points) have been used. The computed 2-D results are validated against NASA experimental data and AGARD test case. The density residual and lift coefficient convergence history is presented in detail. The maximum speed up obtained by use of implicit procedure as compared to explicit one is close to 6 and 14 for unstructured and structured point distributions respectively. The transonic flow over ONERA M6 wing is a classic test case for CFD validation because of simple geometry and complex flow. It has sweep angle of 30° and 15.6° at leading edge and trailing edge respectively. The taper ratio and aspect ratio of the wing are 0.562 and 3.8 respectively. At M∞=0.84 and α=3.06° lambda shock appear on the upper surface of the wing. 3¬D explicit and implicit solvers are tested on ONERA M6 wing. The computed pressure coefficients are compared with experiments at section of 20%, 44%, 65%, 80%, 90% and 95% of span length. The computed results are found to match very well with experiments. The speed up obtained from implicit procedure is over 7 for ONERA M6 wing. The determination of the aerodynamic characteristics of a wing with the control surface deflection is one of the most important and challenging task in aircraft design and development. Many military aircraft use some form of the delta wing. To demonstrate the effectiveness of 3-D solver in handling control surfaces and small gaps, implicit 3-D code is used to compute flow past clipped delta wing with aileron deflection of 6° at M∞ = 0.9 and α = 1° and 3°. The leading edge backward sweep is 50.4°. The aileron is hinged from 56.5% semi-span to 82.9% of semi-span and at 80% of the local chord from leading edge. The computed results are validated with NASA experiments

Este trabalho é dedicado à análise e implementação de esquemas "upwind" de alta ordem modernos e o modelo de turbulência k-epsilon padrão no Freeflow-2D; um ambiente integrado para simulação numérica em diferenças finitas de problemas de escoamentos incompressíveis com superfícies livres. O propósito do estudo é a simulação de escoamentos de fluidos newtonianos incompressíveis, bidimensionais, confinados e/ou com superfícies livres e a altos valores do número de Reynolds. O desempenho do código Freeflow-2D atual é avaliada na simulação do escoamento numa expansão brusca e de um jato livre incidindo perpendicularmente sobre uma superfície rígida impermeável. O código é então aplicado na simulação de um jato planar turbulento em uma porção de fluido com superfície livre e estacionário. Os resultados numéricos obtidos são comparados com dados experimentais, soluções analíticas e soluções numéricas de outros trabalhos. / This work is devoted to the analysis and implementation of modern high-order upwind schemes and the standard k-epsilon turbulence model into the Freeflow-2D; a finite difference integrated environment for the numerical simulation of incompressible free surface flow problems. The purpose of this study is the two-dimensional simulation of high-Reynolds incompressible newtonian confined and/or free surface flows. The performance of the current Freeflow-2D code is assessed by applying it to the simulation of flow over a backward facing step and of an impinging free jet onto an impermeable rigid surface. The code is then applied to a turbulent planar jet into a pool. The numerical results are compared with experimental data, analytical solution, and numerical simulations of other works.

Este trabalho de pesquisa é dedicado ao desenvolvimento, análise e implementação de um novo esquema upwind de alta resolução (denominada PFDPUS) para a aproximação de termos convectivos em leis de conservação e problemas relacionados em mecânica dos fluídos. Usando variáveis normalizadas de Leonard, o equema PFDPUS é baseado em uma função polinomial por partes que satisfaz os critérios de estabilidade CBC e TVD. O desempenho do esquema PEDPUS é investigado na solução das equações de advecção de escalares, Burgers, Euler e MHD. O novo esquema é então aplicado para simular escoamentos incompressíveis envolvendo superfícies livres móveis. Para tanto, o esquema PFDPUS é implementado dentro do software CLAWPACK para problemas compressíveis, e no código Freeflow para poblemas incompressíveis. Os resultados numéricos são comparados com dados experimentais, numéricos e analíticos / This work is dedicated to the development, analysis and implementation of a new high-resolution upwind scheme (called PFDPUS) for approximation of convective terms in conservation laws and related fluid mechanics problems. By using the normalized variables of Leonard, the PFDPUS scheme is based on a piecewise polynomical function that satisfies the CBC and TVD stability criteria. The performance of the PFDPUS scheme is assessed by solving advection of scalars, Burgers, Euler and MHD equations. Then the new scheme is applied to simulate incompressible flows involving moving free surfaces. The PFDPUS scheme is implemented into the CLAWPACK software for compressible problems, and in the Freeflow code for incompressible problems. The numerical results are compared with experimental, numerical and analytical data

In this work, we deal with the finite element method Streamline Upwind/Petrov-Galerkin (SUPG) and use it to solve boundary value problem for the stationary convection-diffusion equation with dominant convection with Dirichlet boundary condition on the whole boundary of bounded polyhedral computational domain of dimension 1 and 2, respectively. We consider a quadratic Lagrangian finite elements on the line segments and triangles, respectively. The core of the work is a proposition of choice of stabilizing parameter of SUPG method as an elementwise affine function in outflow boundary layer and as an elementwise constant function in the rest of the computational domain. We show that this choice gives a more accurate solution than the choice of the stabilization parameter as a constant in each element. 1

Computational fluid dynamics (CFD) is a significant tool routinely used indesign and optimization in aerospace industry. Often cases with unsteadyflows must be computed, and the long compute times of standard methods hasmotivated the present work on new implicit methods to replace the standardexplicit schemes. The implementation and numerical experiments were donewith the Swedish national flow solver Edge, developed by FOI,universities, and collaboration partners.The work is concentrated on a Lower-Upper Symmetric Gauss-Seidel (LU-SGS)type of time stepping. For the very anisotropic grids needed forReynolds-Averaged Navier-Stokes (RANS) computations of turbulent boundary layers,LU-SGS is combined with a line-implicit technique.  The inviscid flux Jacobians which contribute to the diagonalblocks of the system matrix are based on a flux splitting method with upwind type dissipation giving  control over diagonal dominance and artificial dissipation.The method is  controlled by several parameters, and comprehensivenumerical experiments were carried out to identify their influence andinteraction so that close to optimal values can be suggested. As an example,the optimal number of iterations carried out in a time-step increases with increased resolution of the computational grid.The numbering of the unknowns is important, and the numberings produced by mesh generators of Delaunay- and advancing front-type wereamong the best.The solver has been parallelized with the Message Passing Interface (MPI) for runs on multi-processor hardware,and its performance scales with the number of processors at least asefficiently as the explicit methods. The new method saves typicallybetween 50 and 80 percent of the runtime, depending on the case, andthe largest computations have reached 110M grid nodes. Theclassical multigrid acceleration for 3D RANS simulations was foundineffective in the cases tested in combination with the LU-SGS solverusing optimal parameters. Finally, preliminary time-accurate simulations for unsteady flows have shown promising results. / <p>QC 20141201</p>

An unstructured explicit, Reynolds averaged Navier-Stokes solver is developed to operate on inviscid flows, laminar flows and turbulent flows and one equation Spalart-Allmaras turbulence modeling is implemented to the solver. A finite volume formulation, which is cell-center based, is used for numerical discretization of Navier-Stokes equations in conservative form. This formulation is combined with one-step, explicit time marching upwind numerical scheme that is the first order accurate in space. Turbulent viscosity is calculated by using one equation Spalart-Allmaras turbulence transport equation. In order to increase the convergence of the solver local time stepping technique is applied. Eight test cases are used to validate the developed solver,for inviscid flows, laminar flows and turbulent flows. All flow regimes are tested on NACA-0012 airfoil. The results of NACA-0012 are compared with the numerical and experimental data.

No description available.

No

This thesis describes several approaches for stabilizing POD-ROMs (that is, reduced order models based on basis functions derived from the proper orthogonal decomposition) for both the CDR (convection-diffusion-reaction) equation and the NSEs (Navier-Stokes equations). Stabilization is necessary because standard POD-ROMs of convection-dominated problems usually display numerical instabilities. The first stabilized ROM investigated is a streamline-upwind Petrov-Galerkin ROM (SUPG-ROM). I prove error estimates for the SUPG-ROM and derive optimal scalings for the stabilization parameter. I test the SUPG-ROM with the optimal parameter in the numerical simulation of a convection-dominated CDR problem. The SUPG-ROM yields more accurate results than the standard Galerkin ROM (G-ROM) by eliminating the inherent numerical artifacts (noise) in the data and dampening spurious oscillations. I next propose two regularized ROMs (Reg-ROMs) based on ideas from large eddy simulation and turbulence theory: the Leray ROM (L-ROM) and the evolve-then-filter ROM (EF-ROM). Both Reg-ROMs use explicit POD spatial filtering to regularize (smooth) some of the terms in the standard G-ROM. I propose two different POD spatial filters: one based on the POD projection and a novel POD differential filter. These two new Reg-ROMs and the two spatial filters are investigated in the numerical simulation of the three-dimensional flow past a circular cylinder problem at Re = 100. The numerical results show that EF-ROM-DF is the most accurate Reg-ROM and filter combination and the differential filter generally yields better results than the projection filter. The Reg-ROMs perform significantly better than the standard G-ROM and decrease the CPU time (compared against the direct numerical simulation) by orders of magnitude (from about four days to four minutes). / Ph. D.