Écriture de la matrice des différences finies de l'équation de Poisson avec les conditions aux limites de Neumann

Je suis intéressé à résoudre l'équation de Poisson en utilisant l'approche des différences finies. Je voudrais mieux comprendre comment écrire l'équation matricielle avec les conditions aux limites de Neumann. Quelqu'un pourrait-il revoir ce qui suit, est-ce exact?

La matrice des différences finies

L'équation de Poisson,

\frac{\partial^{2} u (X)}{\partial X^{2}} = ré (X)

$\frac{\partial^2u(x)}{\partial x^2} = d(x)$

peut être approximée par une équation matricielle à différence finie,

\frac{1}{(Δ X)^{2}} M ∙ \hat{u} = \hat{ré}

$\frac{1}{(\Delta x)^2} \textbf{M}\bullet \hat u = \hat d$

où est un matrice et et sont vecteurs (colonne), $\textbf{M}$ $n \times n$ $\hat u$ $\hat d$ $1 \times n$

Matrice de différences finies de l'équation de Poisson

Ajout d'une condition aux limites Neumann

Une condition aux limites de Neumann impose un flux connu à la frontière (ici, nous l'appliquons sur le côté gauche où la frontière est à ), $x=0$

écrivant cette condition aux limites comme une différence finie centrée,

\frac{\partial u (X = 0)}{\partial X} = σ

$\frac{\partial u(x=0)}{\partial x} = \sigma$

Erreur dans l'équation. NB. À l'origine, j'ai fait une erreur ici, j'ai signé une erreur et je n'ai pas divisé par 2. Ce qui suit a été corrigé.

\frac{u_{2} - u_{0}}{2 Δ X} = σ

$\frac{u_2 - u_0}{2\Delta x} = \sigma$

Notez l'introduction d'un point de maillage en dehors du domaine d'origine ( ). Ce terme peut être éliminé en introduisant la deuxième équation, $u_0$

\frac{u_{0} - 2 u_{1} + u_{2}}{(Δ X)^{2}} = {ré}_{1}

$\frac{u_0 - 2u_1 + u_2}{(\Delta x)^2} = d_1$

L'équation ne contient plus d'informations en raison de l'introduction du nouveau point de maillage. Il nous permet d'écrire la dérivée double de comme frontière en termes de utilisant une différence finie centrée. $u_1$ $u_0$

La partie dont je ne suis pas sûr

La combinaison de ces deux équations peut être éliminée. Pour montrer le fonctionnement, réorganisons d'abord l'inconnu, $u_0$

u_{0} = - 2 σ Δ X + u_{2} u_{0} = (Δ X)^{2} {ré}_{1} + 2 u_{1} - u_{2}

$u_0 = -2\sigma\Delta x + u_2 \\ u_0 = (\Delta x)^2 d_1 + 2 u_1 - u_2$

Ensuite, ils sont définis égaux et réorganisés dans la forme,

\frac{u_{2} - u_{1}}{(Δ X)^{2}} = \frac{{ré}_{1}}{2} + \frac{σ}{Δ X}

$\frac{u_2 - u_1}{(\Delta x)^2} = \frac{d_1}{2} + \frac{\sigma}{\Delta x}$

$u$ $(\Delta x)^2$

Enfin, en utilisant cette équation comme première ligne de la matrice,

Équation de Poisson avec une condition aux limites de Neumann sur le côté gauche (corrigée)

Quelques dernières réflexions,

Cette matrice finale est-elle correcte?
Aurais-je pu utiliser une meilleure approche?
Existe-t-il un moyen standard d'écrire cette matrice?

— boyfarrell
source

2 Δ x

$2 \Delta x$

u_{0} = - σ Δ x + u_{2}

$u_0 = - \sigma \Delta x + u_2$

Ceci est très bien réglé dans le texte de différence finie de LeVeque , chapitre 2.

— David Ketcheson

Ces problèmes sont également bien expliqués dans scientificpython.net/1/post/2013/01/…

— Evgeni Sergeev

pourriez-vous s'il vous plaît voir ce scicomp.stackexchange.com/questions/14306/…

— usumdelphini

Réponses:

Je pense que vous êtes sur la bonne voie. Si vous corrigez vos erreurs, cela ressemblera beaucoup à http://www.math.toronto.edu/mpugh/Teaching/Mat1062/notes2.pdf .

— vanCompute
source

Merci pour les commentaires et corrections. J'ai mis à jour ma question avec des corrections. Je pense que c'est correct maintenant.

— boyfarrell

$u_0$

Prenez du recul et réfléchissez au problème pendant une seconde. La spécification d'une équation de Laplace indique fondamentalement que chaque point est la moyenne de ses voisins. Ceci est généralement visualisé comme une feuille de caoutchouc et m'aide à penser à ces choses. (Poisson est similaire avec des points plus ou moins extensibles)

Lorsque vous spécifiez la valeur de la surface de la solution sur les bords les plus externes, vous "épinglez" la feuille dans l'espace à ces points. Lorsque vous spécifiez la feuille par sa dérivée sur les bords, il existe un certain nombre de solutions qui remplissent l'équation de translation de la feuille dans l'espace tout en conservant la même forme réelle et donc des dérivés.

$u_0 = 0$

— meawoppl
source

Donc, généralement l'équation de Poisson est résolue avec au moins une condition aux limites de Dirichlet, de sorte qu'une solution unique puisse être trouvée? Je suppose qu'il est logique que les conditions aux limites de Neumann n'aient de sens que lorsque la source et les puits sont inclus, sinon il existe un nombre infini de solutions. Cependant, si je prends l'équation de diffusion à la place, des conditions aux limites de Neumann sont parfois nécessaires pour la physique correcte (par exemple, aucun flux de la quantité à travers une frontière lorsque du / dx = 0). C'est ce qui m'intéresse vraiment. La méthode ci-dessus est-elle la bonne approche pour appliquer les BC Neumann?

— boyfarrell

Vous ne pouvez pas appliquer les BC Neumann sur tous les côtés de votre papier. Si vous le faites, vous n'aurez pas de solution unique. Il doit être fixé au moins sur un côté.

— vanCompute

@meawoppl: Comment spécifier le point fixe tout en faisant une résolution directe de matrice?

— jvriesem

En règle générale, il suffit d'attribuer un point à une constante, en définissant un seul terme d'affilée sur 1, le reste zéro et une valeur sur le RHS qui correspond au plan de solution que vous souhaitez voir.

— meawoppl