We present an Embedded Boundary with Adaptive Mesh Refinement technique for solving the compressible Navier Stokes equations in arbitrary complex domains; followed by a numerical studies for the effect of circular cylinders on the transient dynamics of the Richtmyer-Meshkov Instability. A PDE multidimensional extrapolation approach is used to reconstruct the solution in the ghost-fluid regions and imposing boundary conditions on the fluid-solid interface, coupled with a multi-dimensional algebraic interpolation for freshly cleared cells. The Navier Stokes equations are numerically solved by the second order multidimensional upwind method. Block-structured AMR, implemented with the Chombo framework, is utilized to reduce the computational cost while keeping high resolution mesh around the Embedded Boundary and regions of high gradient solutions. The versatility of the method is demonstrated via several numerical examples, in both static and moving geometry, ranging from low Mach number nearly incompressible to supersonic flows. Our simulation results are extensively verified against other numerical results and validated against available experimental results where applicable. The effects on the transient dynamics of the Richtmyer-Meshkov instability due to small scale perturbations introduced on the shock-wave or the material interface by a single or set of solid circular cylinders were computationally investigated using the developed technique. First, we discuss the RMI initiated on a flat interface by a rippled shock-wave that is disturbed by a single circular cylinder. Then, we study the effect of introducing a number of circular cylinders on the interface. The arrangement of the cylinders set mimic (in a two dimensional domain) the presence of the solid supporting grid wires used in the formation of the material interface in the experimental setup. We analyzed their effects on the mixing layer growth and the mixedness level, and qualitatively demonstrate the cylinders' perturbation effects on the mixing layer structure. We modeled the cylinders' influence based on their diameters; and showed the model ability to predict the variation of the mixing layer growth for different flow parameters.

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.

O objetivo desse trabalho consiste no estudo de métodos de identificação de parâmetros em problemas envolvendo a equação de advecção-difusão 2D. Essa equação é resolvida utilizando o método dos volumes finitos, sendo empregada métodos de reconstrução de alta ordem em malhas não-estruturadas de triângulos para calcular os fluxos nas faces dos volumes de controle. Como ferramenta de busca dos parâmetros é empregada a técnica baseadas em gradientes, sendo os mesmos calculados utilizando processos baseados em métodos adjuntos. / The aim of this work concern to study parameter identification methods on problems involving the advection-diffusion equation in two dimensions. This equation is solved employing the finite volume methods, and high-order reconstruction methods, on triangle unstructured meshes to solve the fluxes across the faces of control volumes. As parameter searching tool is employed technicals based on gradients. The gradients are solved using processes based on adjoint methods.

O objetivo desse trabalho consiste no estudo de métodos de identificação de parâmetros em problemas envolvendo a equação de advecção-difusão 2D. Essa equação é resolvida utilizando o método dos volumes finitos, sendo empregada métodos de reconstrução de alta ordem em malhas não-estruturadas de triângulos para calcular os fluxos nas faces dos volumes de controle. Como ferramenta de busca dos parâmetros é empregada a técnica baseadas em gradientes, sendo os mesmos calculados utilizando processos baseados em métodos adjuntos. / The aim of this work concern to study parameter identification methods on problems involving the advection-diffusion equation in two dimensions. This equation is solved employing the finite volume methods, and high-order reconstruction methods, on triangle unstructured meshes to solve the fluxes across the faces of control volumes. As parameter searching tool is employed technicals based on gradients. The gradients are solved using processes based on adjoint methods.