Papel de las condiciones de contorno (p. Ej., Periódico) en la ecuación de Poisson


8

Dada la ecuación de Poisson 3D y el lado derecho y el dominio, soy libre de imponer cualquier condición de límite (BC) en la función , ¿o tienen que ser de alguna manera consistentes con el lado derecho? En particular, si impongo BC periódica, ¿habrá exactamente una solución para cualquier lado derecho?

2ϕ(x,y,z)=f(x,y,z)
ϕ

Por ejemplo, let: y resuelvo en un cuadro . Ahora cualquier solución debe ser una suma de donde: porque , y es cualquier función armónica (es decir, ). ¿Correcto?

f(x,y,z)=3π2sin(πx)sin(πy)sin(πz)
(0,1)×(0,1)×(0,1)ϕ0+ϕ1
ϕ0(x,y,z)=sin(πx)sin(πy)sin(πz),
2ϕ0=fϕ12ϕ1=0

Si impongo cero Dirichlet BC, entonces es la única solución, porque satisface el BC, satisface la ecuación y la solución debe ser única. ¿Correcto?ϕ0

¿Qué pasa si impongo BC periódica? ¿Significa que habrá alguna función armónica tal que satisfaga el BC periódico y resuelva la ecuación? ¿Qué es esto explícitamente en este caso?ϕ1ϕ0+ϕ1ϕ1


Técnicamente, esta pregunta es más apropiada para las matemáticas.
David Ketcheson el

No lo mencioné en la pregunta explícitamente, pero supongo implícitamente en el contexto de un solucionador numérico basado en elementos finitos (es por eso que estoy interesado en la respuesta), pero la pregunta en sí es general.
Ondřej Čertík

Respuestas:


10

Desde una perspectiva numérica, quizás sea más fácil discutir las discretizaciones directamente.

Para la ecuación de Poisson con condiciones de contorno de Dirichlet homogéneas, existe una solución única para cualquier lado derecho. Una vez discretizada, la ecuación se puede escribir en la forma , donde es la discretización estándar del operador laplaciano 3D con límite de Dirichlet es una discretización estándar de . Como es positivo definitivo, es invertible, y el sistema tendrá una solución única para cualquier , y por lo tanto, cualquier .A b f A b fAx=bAbfAbf

Existen, por supuesto, los problemas habituales con el submuestreo; Si dos valores diferentes de dan lugar a la misma discretización debido al uso de una grilla gruesa o debido a discontinuidades en , puede haber cierta ambigüedad sobre qué sistema se está resolviendo realmente. Pero siempre que la discretización de se comporte bien, existirá una solución única y significativa.f ffff

La situación es un poco más complicada en el caso de condiciones de límite periódicas porque la discretización estándar del operador laplaciano 3D con límite periódico es semidefinida positiva y tiene un núcleo unidimensional comprende soluciones de la forma con constante.x C CKxCC

Debido a que todavía es simétrica en el caso periódico, tenemos que , por lo que no tendrá una solución a menos que , donde es el vector que consiste en todos los 1s. Esto proporciona la condición de consistencia para el lado derecho en forma discreta.Rango A = K A x = b 1b = 0 1ARangeA=KAx=b1b=01

Tenga en cuenta que, analíticamente, hay una forma algo más simple de verlo. Recuerde que, para , donde es nuestro dominio, tenemos Si estipulamos una condición de límite periódica en , entonces el término de límite en el lado derecho desaparece, y nos queda con Entonces, si satisface , se deduce inmediatamente que debemos tener Este es el análogo analítico deΩ Ω Δ ϕϕC2(Ω)Ωϕ Ω Δ ϕ

ΩΔϕdx=ΩϕndS.
ϕϕ Δ ϕ = f Ω f
ΩΔϕdx=0.
ϕΔϕ=f1b = 0 f b
ΩFreX=0.
1si=0 0, ya que ambos expresan el hecho de que el valor promedio de , y por lo tanto , sobre el dominio debe ser cero.Fsi

John, ¡muchas gracias por esta increíble respuesta! Proporciona una idea clara de cómo funciona. Sentí que hay una trampa para la condición periódica, pero no pude señalarla con el dedo. Gracias de nuevo.
Ondřej Čertík

6

Sé que han pasado 2 años ... pero aquí hay un ejemplo de la física para martillarlo en casa. Considere que está resolviendo la ecuación de Poisson: para potenciales electrostáticos (ϕ) con el lado derecho: Y, ρ (la densidad de carga) se especifica en toda la cuadrícula puntos.f = ( - ρ / ε o )

2ϕ=F
F=(-ρ/ /εo)

Declaración: Para la solución periódica, la condición ∫ρdv = 0 (descrita por Ben arriba) establece que el sistema sea neutral en la red. Esta es una buena manera (física) de pensar sobre el valor promedio de f sobre el cuadro periódico.

Ahora probemos esta afirmación. Para condiciones de contorno periódicas, la integral del campo eléctrico debe ser cero sobre la superficie de la caja (siempre puede encontrar pares de puntos en la superficie de la caja que se cancelan entre sí en la integral).

Luego, según el teorema de divergencia (ley de Gauss), la caja debe ser neutral y esto es lo que nos propusimos probar.

SminortereUNA=1ε0 0Qyonortesyoremi=0 0

Un poco de física aquí, pero creo que proporciona un buen refuerzo para la discusión aquí.


0

Es necesario tener cuidado específicamente en el método de generación de soluciones fabricadas para la ecuación de Poisson. Dado que la definición del término fuente tiene que satisfacer tanto la forma fuerte de la PDE como la forma débil de la PDE. La derivación dada anteriormente es básicamente usar la forma débil de la PDE. En otras palabras, existe una relación entre el gradiente de la solución en el límite y el término fuente integral en el dominio.

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.