Convergencia asintótica de la solución a un pde parabólico a la solución de un pde elíptico

8

Supongamos que tengo el sistema parabólico con condiciones de contorno de Dirichlet y condición inicial u = g ,

ut=(k(x)u)+f,(x,t)Ω×I
u ( x , t ) = h ,
u=g,xΩ
u(x,t)=h,t=0.

Muchas veces en ingeniería, estamos más interesados ​​en el comportamiento asintótico (estado estacionario) de esta PDE que en el comportamiento transitorio. Entonces, a veces descuidamos el término derivado del tiempo y resolvemos el sistema elíptico lugar. La suposición es que durante un tiempo infinito, .lim t u p a r a b o l i c ( x , t ) = u e l l i p t i c ( x , t )

(k(x)u)=f,(x,t)Ω×I
limtuparabolic(x,t)=uelliptic(x,t)

He observado que cuando f0 , este límite es verdadero, pero no estoy seguro de si este es el caso para f arbitrario fo si hay otras condiciones necesarias para garantizar que este límite sea verdadero. ¿Las condiciones de contorno tienen que converger asintóticamente a un valor constante para que la solución parabólica converja a la solución elíptica?

Aunque mi pregunta está formulada en el caso continuo, también tengo curiosidad por saber si las mismas condiciones son verdaderas para el caso discreto. Es decir, suponiendo que use un esquema de diferencia finita estable y consistente para aproximar uparabolic y uelliptic , si debo esperar

limtuparabolicfdm(x,t)=uellipticfdm(x,t)
si uellipticfdm y uparabolicfdm se discretiza en la misma cuadrícula espacial y Δt ?
Paul
fuente

Respuestas:

3

Sí, con las condiciones de contorno de Dirichlet siempre tiene una convergencia exponencial al estado de staedy. Cualquier libro PDE tendrá una prueba. Para una buena explicación desde una perspectiva numérica, vea el capítulo 2 del libro FDM de LeVeque.

David Ketcheson
fuente
Lo mismo debería ser cierto si hubiera condiciones mixtas (neumann y dirichlet), ¿verdad?
Paul
@Paul Sí, el único caso difícil en 1D es cuando ambos límites son Neumann.
David Ketcheson el