Resolvo equações de Navier-Stokes compressíveis para várias espécies em uma grade estruturada em 3D. Eu obtive uma solução em uma determinada grade (digamos que seja relativamente grossa). Quero agora refinar minha grade e interpolar minha solução anterior em minha nova grade antes de reiniciar minha simulação. Atualmente, temos uma ferramenta de interpolação que cria uma árvore kd das 2 grades e, em seguida, pode usar 2 métodos diferentes para calcular os valores na nova grade:
- média simples
- distância inversa ponderada (IDW)
- mínimos quadrados móveis (MLS)
Quero focar na precisão, porque, como lida com grandes gradientes, não capturá-los corretamente gerará ondas quando reiniciar meu cálculo. Primeiro tentei calcular a média, mas a precisão não era boa o suficiente.
Eu pensei que o método MLS com polinômios de ordem 2 me daria resultados razoáveis, pois é suposto não ser oscilatório. No entanto, quando olho para o meu campo interpolado, vejo mínimos / máximos locais que excedem os valores do meu campo inicial. Isso significa que a implementação do MLS neste programa não está correta? Devo ter cuidado com o tamanho do meu estêncil e a ordem dos polinômios? Que outro método você recomendaria?
Desde já, obrigado !
fonte