Siguiendo mi pregunta anterior , estoy tratando de aplicar condiciones de contorno a esta malla de volumen finito no uniforme,
Me gustaría aplicar una condición de límite de tipo Robin a las lhs del dominio ( , de modo que,
donde es el valor límite; son coeficientes definidos en el límite, advección y difusión respectivamente; , es la derivada de evaluada en el límite es la variable por la que estamos resolviendo.
Posibles enfoques
Puedo pensar en dos formas de implementar esta condición de límite en la malla de volumen finito anterior:
Un enfoque de célula fantasma.
Escriba como una diferencia finita, incluida una celda fantasma.
A. Luego use interpolación lineal con puntos y para encontrar el valor intermedio, .
B. Alternativamente encuentre promediando sobre las celdas,
En cualquier caso, la dependencia de la célula fantasma puede eliminarse de la forma habitual (mediante sustitución en la ecuación de volumen finito).
Un enfoque de extrapolación.
Ajuste una función lineal (o cuadrática) a utilizando los valores en los puntos ( ). Esto proporcionará el valor en . La función lineal (o cuadrática) se puede diferenciar para encontrar una expresión para el valor de la derivada, , en el límite. Este enfoque no utiliza una célula fantasma.
Preguntas
- ¿Qué enfoque de los tres (1A, 1B o 2) es "estándar" o recomendaría?
- ¿Qué enfoque introduce el error más pequeño o es el más estable?
- Creo que puedo implementar el enfoque de célula fantasma yo mismo, sin embargo, ¿cómo se puede implementar el enfoque de extrapolación? ¿Este enfoque tiene un nombre?
- ¿Hay alguna diferencia de estabilidad entre ajustar una función lineal o una ecuación cuadrática?
Ecuación específica
Deseo aplicar este límite a la ecuación de advección-difusión (en forma de conservación) con un término fuente no lineal,
Discretización esta ecuación en lo anterior con malla usando el -method da,
Sin embargo, para el punto límite ( ) prefiero usar un esquema totalmente implícito ( θ = 1 ) para reducir la complejidad,
Observe el punto fantasma , esto se eliminará aplicando la condición de contorno.
Los coeficientes tienen las definiciones,
Todas las variables " " se definen como en el diagrama anterior. Finalmente, Δ t , que es el paso de tiempo ( NB este es un simplificado caso con constantes a y d coeficientes, en la práctica los " r coeficientes" son un poco más complicado por esta razón).
fuente
Respuestas:
Esto es más bien una observación general sobre FVM que una respuesta a las preguntas concretas. Y el mensaje es que no debería existir la necesidad de una discretización ad hoc de las condiciones de contorno.
A diferencia de los métodos FE o FD, donde el punto de partida es una respuesta discreta para la solución, el enfoque FVM deja la solución intacta (al principio) pero promedia una segmentación del dominio. La discretización de la solución entra en juego solo cuando el sistema de ecuaciones de equilibrio obtenido se convierte en un sistema de ecuaciones algebraicas al aproximar los flujos a través de las interfaces.
En este sentido, en vista de las condiciones de contorno, aconsejo mantener la forma continua de la solución el mayor tiempo posible e introducir las aproximaciones discretas solo al final.
Digamos, la ecuación mantiene en todo el dominio. Luego se mantiene en el subdominio [ 0 , h 1 ) , y una integración en el espacio da ∫ h 1 0 u t d x
Pero ahora, para convertir esto en una ecuación algebraica, normalmente se supone que en la celda la función u es constante en el espacio, es decir, u ( t , x ) | C i = u i ( t ) . Así, habiendo asociado u ( x i ) ≈ u i , se puede expresar u x | h i en las celdas a través del cociente de diferencia en u i y u i + 1Cyo tu u ( t , x ) |Cyo= uyo( t ) u ( xyo) ≈ uyo tuXEl |hyo tuyo tui + 1 tu
DadoT |x = 0= gre uno puede introducir una célula fantasma y la condición de que un interpolante entre tu0 0 y tu1 es igual a solre en el huésped.
DadotuXEl |x = 0= gnorte uno puede introducir una celda fantasma y la condición de que una aproximación a la derivada entre tu0 0 y tu1 partidos solnorte en el huésped
Si se prescribe el flujo en sí:( - a u + dtuX) |x = 0= gR , no hay necesidad de una discretización.
Sin embargo, no estoy seguro de qué hacer en caso de que haya Robin tipo bc que no coincidan directamente con el flujo. Esto necesitará cierta regularización debido a la discontinuidad de los parámetros de advección y difusión.
===> Algunos pensamientos personales sobre FVM <===
fuente