Estudo numérico da sedimentação em correntes de turbidez com evolução do relevo de fundo

Abstract

Correntes de densidade são fluxos gravitacionais gerados pela diferença de densidade entre dois fluidos. Correntes de turbidez fazem parte de uma sub-classificação das correntes de densidade, na qual o fluido mais denso tem, na sua composição, partículas em suspensão. Muitos trabalhos numéricos já estudaram a dinâmica das correntes de turbidez, mas, nenhum dos encontrados aplicou mudanças de relevo concomitantes com a simulação, causadas pela sedimentação das próprias partículas da corrente e nem alterou o relevo após a passagem de cada evento em um domínio tridimensional. O presente trabalho pretende analisar a alteração no relevo de fundo causada por uma corrente de turbidez. No código Incompact3d, as equações de Navier-Stokes, Continuidade e Transporte e Difusão são resolvidas em uma malha cartesiana tridimensional. A condição inicial adotada é a de Lock-Exchange. As simulações realizadas utilizaram Simulação Numérica Direta (DNS). O código utiliza um esquema compacto centrado de sexta ordem, em diferenças finitas, para o esquema espacial, e Adams-Bashfort de terceira ordem para o esquema temporal. A validação do código foi realizada comparando-se com trabalhos experimental e numéricos. A análise das diferentes proporções granulométricas mostrou que quanto maior é a quantidade de material grosso na condição inicial, maior será seu depósito para um dado tempo. Em consequência, mais relevante se torna a consideração da alteração do relevo de fundo. Além disso, quanto maior o fator de compactação do sedimento, maior será o erro de não considerar a atualização de fundo. Os resultados também apontaram que os erros médios ao não considerar a atualização do fundo são da ordem de 4% da massa de depósito em 20 tempos adimensionais, para os parâmetros utilizados. Ao se propagar uma corrente de turbidez sobre o depósito de outra, os erros se mostram menores.

Gravity currents are gravitational fluxes triggered by density difference between two fluids. A sub-classification of those are turbidity currents, in which the denser fluid is composed by the lighter fluid plus suspended particles. Many papers had shown turbidity currents dynamics, although none of the papers found had applied changes in the simulated topography due to deposit during the own simulation, neither they had altered a 3D domain topography after each flux, applying the changes caused by the previous current. The present dissertation aims to analyse the turbidity current dynamics alteration caused by the influence of its own deposit, altering the topography during the very simulation. The analysis is conducted in a polidispersed turbidity current. The Incompact3d code solves Navier-Stokes, continuity and transport-diffusion equation, in a tridimensional cartesian mesh. Lock-exchange was chosen to be the initial condition. Direct Numerical Simulations (DNS) are performed. Sixth order compact finite-difference schemes are used on the spatial domain, while third order Adams-Bashfort is applied for the temporal evaluation. Comparisons with numerical and experimental papers were performed for code verification. Results showed the coarser the particles on the starting lock-exchange, the higher its deposit is, and the more the terrain will be altered. Nevertheless, the bigger the compacting factor, the bigger the error of not considering bathymetry alteration. Results also point that the average errors of not considering the update are in order of 4% on the mass deposit, after 20 dimensionless times, for the used parameters. When a current propagates over the deposit of a previous one, these errors are smaller.