Uso de iteración de punto fijo para desacoplar un sistema de pde


12

Supongamos que tengo un problema de valor límite:

du

d2udx2+dvdx=f in Ω
u=h en Ω
dudx+d2vdx2=g in Ω
u=h in Ω

Mi objetivo es descomponer la solución de este problema acoplado en una secuencia de PDE no acoplados. Para desacoplar el sistema, estoy aplicando una iteración de punto fijo sobre una secuencia de aproximaciones modo que(uk,vk)

du k - 1

d2ukdx2+dvk1dx=f
duk1dx+d2vkdx2=g

Teóricamente, esto me permitiría resolver ambas ecuaciones como PDE puramente elípticas. Sin embargo, nunca he visto iteraciones de punto fijo aplicadas a PDE de esta manera. He visto iteraciones de punto fijo aplicadas a las ecuaciones numéricamente discretizadas (método de diferencias finitas, método de elementos finitos, etc.), pero nunca a las ecuaciones continuas directamente.

¿Estoy violando algún principio matemático descarado al hacer esto? ¿Es esto matemáticamente válido? ¿Podría resolver el PDE acoplado como una secuencia de PDE no acoplados utilizando la iteración de punto fijo aplicada al problema de la variable CONTINUO, en lugar del problema de la variable DISCRETA?

En este punto, no estoy realmente preocupado sobre si es práctico usar este método, sino más bien si es teóricamente plausible. Cualquier comentario sería muy apreciado!


66
En la literatura de PDE hiperbólica, los métodos de división fraccional de pasos y operadores hacen lo que usted describe anteriormente.
Geoff Oxberry

(uk,vk)(uk,pk)

@BillBarth: ¡Sí! Acabo de corregirlo.
Paul

@ GeoffOxberry: Me parece que la división del operador es muy diferente en carácter.
anónimo

@Paul: Puedo pensar en al menos otro problema en el que los "PDE acoplados" se resuelven a través de una iteración de punto fijo (y no solo formulados como problemas de punto fijo): descomposición del dominio, véase, por ejemplo, el método de Neumann-Dirichlet. (la diferencia aquí es que tiene dos PDE pero viven en dominios diferentes y el acoplamiento es solo a través de una interfaz).
anónimo

Respuestas:


11

C(Ω)×C(Ω)

d2ukdx2+dvk1dx=fd2vkdx2+duk1dx=g
(más condiciones de contorno).

Está claro que si esta secuencia converge, será una solución de su conjunto original de PDE.

xkxk+1u0v0

(ukvk)(u^kv^k)q(uk1vk1)(u^k1v^k1)
|q|<1(uk1,vk1)(u^k1,v^k1)

Esta lógica funciona tanto en el espacio continuo como en el discreto.


3
¿No debería ? |q|<1
Paul
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.