N/A / Thesis / Master of Applied Science (MASc)

O estudo do transporte de partículas tem uma importância fundamental em diversas áreas de pesquisas como, por exemplo, na formação de gelo em uma aeronave pois pode afetar a sua sustentação e estabilidade. Tamanha é a preocupação com a segurança de vôo que diversos estudos têm sido realizados, resultando em códigos computacionais como o LEWICE nos Estados Unidos, TRAJICE no Reino Unido, ONERA na França e CANICE no Canadá. No Brasil, um dos estudo é feito pela EMBRAER em parceria com algumas instituições. O objetivo deste trabalho é desenvolver um algoritmo que possa ser empregado na trajetória das partículas, utilizando uma formulação euleriana que elimina a dificuldade da semeadura de partículas específica da formulação lagrangiana na determinação da fração de volume da partícula. O método empregado é dos volumes finitos em malhas não-estruturadas cuja principal chave está na reconstrução mínimos quadrados de alta ordem com restrição nos contornos. O desenvolvimento do trabalho engloba 3 etapas: definição da geometria e geração das malhas; utilização de um solver para o tratamento do escoamento do ar e obtenção do campo de velocidade; implementação e utilização do esquema numérico com reconstrução mínimos quadrados de alta ordem para simular o cálculo da fração de volume com imposição de condições limites apropriadas no contorno do corpo. Os resultados dos testes realizados mostram que o esquema numérico com reconstrução mínimos quadrados pode ser empregado na resolução de equações que apresentam uma região de descontinuidade, como é o caso da região de sombra, reduzindo a largura da banda de difusão numérica e overshoots. / The particle transport study has a fundamental importance in diverse research area like in the icing accretion on an aircraft because that can affect its sustentation and stability. The concern is so big that many researches have been carried through, resulting in computational codes like the LEWICE in the United States, TRAJICE in the United Kingdom, ONERA in France and CANICE in Canada. In Brazil, one of the study has been made by the EMBRAER with some institutes. The goal of this work is to develop an algorithm that can be used in the particles trajectory study, using an Eulerian method that eliminates the difficulty particle sowing, particular of the Lagrangian method, in the determination of the droplet fraction volume. This is made by the finite volume method on unstructured meshes whose main key is the high order reconstruction with restriction on the boundary. The development of the work involves 3 stages: geometry definition and mesh generation; using code for the treatment of the air flow and obtained flow velocity; use of the high order numerical scheme least square reconstruction to simulate the droplet fraction volume result with imposition of appropriate limit conditions in the body contour. The realized simulations shown that Least Square method can be used in problem resolution that present descontinuos region like is shadow region reducing numerical diffusion and overshoots.

O estudo do transporte de partículas tem uma importância fundamental em diversas áreas de pesquisas como, por exemplo, na formação de gelo em uma aeronave pois pode afetar a sua sustentação e estabilidade. Tamanha é a preocupação com a segurança de vôo que diversos estudos têm sido realizados, resultando em códigos computacionais como o LEWICE nos Estados Unidos, TRAJICE no Reino Unido, ONERA na França e CANICE no Canadá. No Brasil, um dos estudo é feito pela EMBRAER em parceria com algumas instituições. O objetivo deste trabalho é desenvolver um algoritmo que possa ser empregado na trajetória das partículas, utilizando uma formulação euleriana que elimina a dificuldade da semeadura de partículas específica da formulação lagrangiana na determinação da fração de volume da partícula. O método empregado é dos volumes finitos em malhas não-estruturadas cuja principal chave está na reconstrução mínimos quadrados de alta ordem com restrição nos contornos. O desenvolvimento do trabalho engloba 3 etapas: definição da geometria e geração das malhas; utilização de um solver para o tratamento do escoamento do ar e obtenção do campo de velocidade; implementação e utilização do esquema numérico com reconstrução mínimos quadrados de alta ordem para simular o cálculo da fração de volume com imposição de condições limites apropriadas no contorno do corpo. Os resultados dos testes realizados mostram que o esquema numérico com reconstrução mínimos quadrados pode ser empregado na resolução de equações que apresentam uma região de descontinuidade, como é o caso da região de sombra, reduzindo a largura da banda de difusão numérica e overshoots. / The particle transport study has a fundamental importance in diverse research area like in the icing accretion on an aircraft because that can affect its sustentation and stability. The concern is so big that many researches have been carried through, resulting in computational codes like the LEWICE in the United States, TRAJICE in the United Kingdom, ONERA in France and CANICE in Canada. In Brazil, one of the study has been made by the EMBRAER with some institutes. The goal of this work is to develop an algorithm that can be used in the particles trajectory study, using an Eulerian method that eliminates the difficulty particle sowing, particular of the Lagrangian method, in the determination of the droplet fraction volume. This is made by the finite volume method on unstructured meshes whose main key is the high order reconstruction with restriction on the boundary. The development of the work involves 3 stages: geometry definition and mesh generation; using code for the treatment of the air flow and obtained flow velocity; use of the high order numerical scheme least square reconstruction to simulate the droplet fraction volume result with imposition of appropriate limit conditions in the body contour. The realized simulations shown that Least Square method can be used in problem resolution that present descontinuos region like is shadow region reducing numerical diffusion and overshoots.

Une méthode complètement Eulérienne reposant sur un modèle 1-fluide est présentée afinde résoudre les problèmes d’interaction fluide-structure élastique. L’interface entre le fluideet la structure élastique est représentée par une fonction level-set, transportée par le champde vitesse du fluide et résolue par un schéma d’ordre élevé WENO 5. Les déformationsélastiques sont calculées sur la grille eulérienne à l’aide des caractéristiques rétrogrades.Nous utilisons différents modèles d’hyperélasticité, afin de générer puis d’intégrer les forcesélastiques comme terme source des équations de Navier Stokes. Le couplage vitesse/pressionest résolu par une méthode de correction de pression et les équations sont discrétisées parla méthode des volumes finis sur la grille eulérienne. La principale difficulté réside dansles grands déplacements de fluide autour du solide, source d’instabilités numériques. Afind’éviter ces problèmes, nous effectuons périodiquement une redistanciation de la level-setet une extrapolation linéaire des caractéristiques rétrogrades. Dans un premier temps,nous effectuons la vérification et la validation de notre approche à l’aide de plusieurs castests comme celui proposé par Turek. Ensuite, nous appliquons notre méthode à l’étudedu phénomène d’atténuation des vagues par des structures élastiques. Il s’agit d’une desvoies possibles pour réduire l’impact des fortes houles sur notre littoral. De plus dans lalittérature et à notre connaissance, seules des structures élastiques rigides ou élastiquesmais monodimensionnelles ont été utilisées pour réaliser ces études. Nous proposons deplacer des structures élastiques sur les fonds marins et analysons leur capacité d’absorptionde l’énergie produite par les vagues. / A fully Eulerian method is developed to solve the problem of fluid-elastic structure interactionsbased on a 1-fluid method. The interface between the fluid and the elastic structureis captured by a level set function, advected by the fluid velocity and solved with a WENO5 scheme. The elastic deformations are computed in an Eulerian framework thanks to thebackward characteristics. We use the Neo Hookean or Mooney Rivlin hyperelastic modelsand the elastic forces are incorporated as a source term in the incompressible Navier-Stokesequations. The velocity/pressure coupling is solved with a pressure-correction methodand the equations are discretized by finite volume schemes on a Cartesian grid. The maindifficulty resides in that large deformations in the fluid cause numerical instabilities. Inorder to avoid these problems, we use a re-initialization process for the level set and linearextrapolation of the backward characteristics. First, we verify and validate our approachon several test cases, including the benchmark of FSI proposed by Turek. Next, we applythis method to study the wave damping phenomenon which is a mean to reduce thewaves impact on the coastline. So far, to our knowledge, only simulations with rigid orone dimensional elastic structure has been studied in the literature. We propose to placeelastic structures on the seabed and we analyse their capacity to absorb the wave energy