Le principe maximum / minimum de l'équation thermique est-il maintenu par la discrétisation de Crank-Nicolson?

10

J'utilise le schéma de différence finie Crank-Nicolson pour résoudre une équation de chaleur 1D. Je me demande si le principe maximum / minimum de l'équation de la chaleur (c'est-à-dire que le maximum / minimum se produit à la condition initiale ou aux limites) est également valable pour la solution discrétisée.

Cela est probablement sous-entendu par le fait que Crank-Nicolson est un schéma stable et convergent. Mais il semble que vous pourriez être en mesure de le prouver directement via un argument d'algèbre linéaire en utilisant les matrices créées à partir du pochoir de Crank-Nicolson.

J'apprécierais toute indication sur la littérature à ce sujet. Merci.

— foobarbaz
source

Salut foobarbaz, et bienvenue sur scicomp! Je suppose que le problème que vous résolvez n'a pas de termes source, n'est-ce pas?

— Paul

8

Le principe maximum pour Crank-Nicolson tiendra si

μ ≐ \frac{k}{h^{2}} \leq 1

$\mu \doteq \frac{k}{h^2} \leq 1$

k

$k$

h

$h$

θ

$\theta$

u^{n + 1} = u^{n} + \frac{μ}{2} ((1 - θ) A u^{n} + θ A u^{n + 1})

$u^{n+1} = u^n + \frac{\mu}{2}\left( (1-\theta)Au^n + \theta Au^{n+1}\right)$

A

$A$

0 \leq θ \leq 1

$0 \leq \theta \leq 1$

μ (1 - 2 θ) \leq \frac{1}{2}

$\mu(1-2\theta) \leq \frac{1}{2}$

μ (1 - θ) \leq \frac{1}{2}

$\mu(1-\theta) \leq \frac{1}{2}$

Pour une preuve, voir Solutions numériques d'équations différentielles partielles par KW Morton . En particulier, regardez les sections 2.10 et 2.11 et le théorème 2.2.

$\mu$

$[0,1]$ $u_i^k$ $k$ $i$ $u^k_0 = u^k_2 = 0$ $k$

(1 - \frac{μ}{2} (- 2)) u_{1}^{n + 1} = (1 + \frac{μ}{2} (- 2)) u_{1}^{n},

$\left(1 - \frac{\mu}{2}(-2)\right)u^{n+1}_1 = \left(1 + \frac{\mu}{2}(-2)\right)u^n_1,$

u_{1}^{n + 1} = (\frac{1 - μ}{1 + μ}) u_{1}^{n} .

$u^{n+1}_1 = \left(\frac{1-\mu}{1+\mu}\right)u^n_1.$

$u_1^0 = 1$

u_{1}^{n} = {(\frac{1 - μ}{1 + μ})}^{n},

$u^n_1 = \left(\frac{1-\mu}{1+\mu}\right)^n,$

u_{1}^{n} \leq 1

$u^n_1 \leq 1$

u_{1}^{n} < 0

$u^n_1 < 0$

n

$n$

μ \leq 1

$\mu \leq 1$

μ \leq 1

$\mu \leq 1$

μ

$\mu$

En réponse à la demande de foobarbaz, j'ai ajouté un croquis de la preuve.

\begin{aligned} (1 + 2 θ μ) u_{j}^{n + 1} & = θ μ (u_{j - 1}^{n + 1} + u_{j + 1}^{n + 1}) \\ + (1 - θ) μ (u_{j - 1}^{n} + u_{j + 1}^{n}) \\ + [1 - 2 (1 - θ) μ] u_{j}^{n} \end{aligned}

$\begin{align*} (1+2\theta\mu)u^{n+1}_j &= \theta\mu(u^{n+1}_{j-1} + u^{n+1}_{j+1})\\ &+ (1-\theta)\mu(u^n_{j-1} + u^n_{j+1})\\ &+ [1-2(1-\theta)\mu]u^n_j \end{align*}$

$\mu(1-\theta)\leq \frac{1}{2}$

$u^{n+1}_j$ $u^{n+1}_{j-1}$ $u^{n+1}_{j+1}$ $u^n_{j-1}$ $u^n_{j+1}$ $u^n_j$ $u^{n+1}_j$ $u^{n+1}_j$

\begin{aligned} (1 + 2 θ μ) u_{j}^{n + 1} & > θ μ (u_{j - 1}^{n + 1} + u_{j + 1}^{n + 1}) \\ + (1 - θ) μ (u_{j - 1}^{n} + u_{j + 1}^{n}) \\ + [1 - 2 (1 - θ) μ] u_{j}^{n} \\ = (1 + 2 θ μ) u_{j}^{n + 1} \end{aligned}

$\begin{align*} (1+2\theta\mu)u^{n+1}_j &> \theta\mu(u^{n+1}_{j-1} + u^{n+1}_{j+1})\\ &+ (1-\theta)\mu(u^n_{j-1} + u^n_{j+1})\\ &+ [1-2(1-\theta)\mu]u^n_j\\ &= (1+2\theta\mu)u^{n+1}_j \end{align*}$

$u^{n+1}_j$ $u$

— Ben
source

Merci! Connaissez-vous une autre référence que Morton? Je ne peux pas accéder à ces sections ou au théorème dans l'aperçu du livre Google. J'aimerais comprendre la preuve.

— foobarbaz

@foobarbaz Je n'ai pas d'autre référence à portée de main, mais j'ai ajouté un aperçu de la preuve. Faites-moi savoir si je peux être plus clair.

— Ben

0

La stabilité signifie qu'une perturbation reste limitée dans le temps. Cela ne signifie pas que le principe maximum est satisfait au niveau discret, c'est une question différente. Satisfaire au principe maximum discret est suffisant mais pas nécessaire pour la stabilité.

— chris
source