El efecto de desacoplar un sistema acoplado de PDE


8

Hice una pregunta algo similar anteriormente, pero tal vez podría haber sido demasiado específica para que alguien realmente respondiera. Aquí hay un poco más general de una pregunta con la que estoy luchando. Considere el siguiente sistema:

u 2

(D1(u2)u1)=f1(u2)
u2t+f2(u1,u2)(D2(u2)u2)=0

suponiendo un conjunto general de BC: D iu in = u i , N ,

ui=ui,D,onΓD
Diuin=ui,N,onΓN

Usando DGFEM para discretizaciones espaciales y Euler hacia atrás para la derivada del tiempo. Nos desacoplamos así:

  1. Resolver para usando u k 2 : - ( D 1 ( u k 2 ) u k + 1 1 ) = f 1 ( u k 2 )u1k+1u2k

    (D1(u2k)u1k+1)=f1(u2k)
  2. Resolver para usando u k + 1 1 : u 2u2k+1u1k+1

    u2t+f2(u1k+1,u2k,u2k+1)(D2(u2k+1)u2k+1)=0

El método de Newton se usa para manejar la no linealidad.

f2(u1,u2)=u2(u1+u2u2)

entonces

f2(u1k+1,u2k,u2k+1)=u2k+1(u1k+1+u2ku2k)

f2u1k+1u2ku2k+1f2u2ku1k+1+u2ku2k

f2=f2(u1,exactk+1,u2k+1,u2,exactk+1)=u2k+1(u1,exactk+1+u2,exactk+1u2,exactk+1)

¿Es este rendimiento subóptimo en el primer escenario esperado dado el desacoplamiento? Sé que el desacoplamiento limitará el paso de tiempo, pero estoy bastante seguro de que estoy tomando el paso de tiempo lo suficientemente pequeño.

f2u2u2u1k+1

No había podido encontrar mucha literatura sobre el efecto que tendría el desacoplamiento en la convergencia, por lo que si alguien pudiera orientarme en la dirección correcta u ofrecer algún consejo, sería genial.


3
No los desacoplas, estás aplicando un método de división.
David Ketcheson

Gracias por la corrección. No me he olvidado de esta pregunta, solo sigo trabajando en ello. Actualizaré si lo hago funcionar ...
Justin Dong

Respuestas:


3

Esto es definitivamente lo que se conoce como un esquema de división y es muy popular en la óptica no lineal, así como en las simulaciones cuánticas, un ejemplo es la ecuación de dirac, pero también puede hacerlo en PDES no lineales '.

u1k+1u2ku1k+1u2k+1

u1u2u1u2Lu1f=f2(u2)G(x;u2)LG=δ(xs)u1

u1=ΩG((xs);u2)fds

u2k+1u1k+1u2k+1

Sin embargo, definitivamente hay flexibilidad en este desacoplamiento o división que podría ser un buen esquema en sus ecuaciones. Se puede encontrar una buena introducción a la división en el contexto que le interesa en el Capítulo 13 'División y sus primos' de Boyd 'Chebyshev y métodos espectrales de Fourier'.

Espero que esto ayude


1

Es un poco difícil de decir sin saber con precisión qué PDE estás resolviendo; Dicho esto, algunos lugares para mirar son

  1. Para los sistemas lineales acoplados , los métodos Uzawa y lagrangianos aumentados son ejemplos de esquemas generales de división. Existen resultados conocidos sobre la convergencia de estos métodos basados ​​en los espectros de los bloques y el complemento de Schur del sistema de bloques acoplados (ver Bramble et al, sobre la resolución de problemas de punto de silla de montar).
  2. Existe una gran cantidad de material sobre diferentes métodos de división para PDE. Aunque están formulados en gran medida para problemas lineales, existe un amplio uso y teoría para su extensión a ecuaciones no lineales, por ejemplo, al menos para las PDE no lineales que rigen los tipos de convección-difusión y flujo incompresible, los esquemas de división / proyección tienen una historia bastante rica .

u1u2u1u2

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.