Comenzando con la ecuación de advección en forma de conservación.
donde es una velocidad que depende del espacio, yu es una concentración de una especie que se conserva.
Discretizar el flujo (donde el flujo , se define en los bordes de las celdas entre los puntos de malla) da, u t = 1
Usando un viento de primer orden, aproximamos los flujos como,
Lo que da, ut=1
Si fue constante, esto se reducirá al esquema familiar de viento arriba, es decir, u t = a.
Mi pregunta es, ¿cómo podemos tratar los coeficientes no constantes de la ecuación de advección? La velocidad se define en los centros celulares, por lo que un enfoque simple sería el siguiente,
Este es mi enfoque preferido porque es muy simple de implementar.
Sin embargo, también podríamos usar (supongo) un esquema de promedio para definir la velocidad en los bordes de la celda,
En el libro de LeVeque dice:
Pero en realidad no elabora demasiado después de eso. ¿Qué es un enfoque común?
Estoy resolviendo un problema de conservación (estoy usando la ecuación de advección como una ecuación de continuidad), por lo que quiero asegurarme de que después de aplicar la discretización se conserve la propiedad de conservación. ¡Me gustaría evitar sorpresas ocultas con respecto a estos coeficientes variables! ¿Alguien tiene algunos comentarios generales y orientación?
Actualización A continuación hay dos respuestas realmente buenas y solo podría elegir una :(