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
suponiendo un conjunto general de BC: D i ∇ u i ⋅ n = u i , N ,
Usando DGFEM para discretizaciones espaciales y Euler hacia atrás para la derivada del tiempo. Nos desacoplamos así:
Resolver para usando u k 2 : - ∇ ⋅ ( D 1 ( u k 2 ) ∇ u k + 1 1 ) = ∇ ⋅ f 1 ( u k 2 )
Resolver para usando u k + 1 1 : ∂ u 2
El método de Newton se usa para manejar la no linealidad.
entonces
¿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.
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.