Resuelvo ecuaciones de Navier-Stokes comprimibles de varias especies en una cuadrícula estructurada en 3D. He obtenido una solución en una cuadrícula dada (digamos una relativamente gruesa). Ahora quiero refinar mi grilla e interpolar mi solución anterior en mi nueva grilla antes de reiniciar mi simulación. Actualmente, tenemos una herramienta de interpolación que construye un árbol kd de las 2 cuadrículas y luego puede usar 2 métodos diferentes para calcular los valores en la nueva cuadrícula:
- promedio simple
- distancia inversa ponderada (IDW)
- mover mínimos cuadrados (MLS)
Quiero centrarme en la precisión porque, dado que trato con grandes gradientes, no capturarlos correctamente generará ondas cuando reinicie mi cálculo. Primero probé el promedio simple pero la precisión no fue lo suficientemente buena.
Pensé que el método MLS con polinomios de orden 2 me daría resultados razonables ya que se supone que no es oscilatorio. Sin embargo, cuando miro mi campo interpolado, veo mínimos / máximos locales que superan los valores de mi campo inicial. ¿Significa esto que la implementación de MLS en este programa no es correcta? ¿Debo tener cuidado con el tamaño de mi plantilla y el orden de los polinomios? ¿Qué otro método recomendarías?
Gracias por adelantado !
fuente