La mise à l'échelle variable est-elle essentielle lors de la résolution numérique de certains problèmes PDE?

Dans la simulation de semi-conducteur, il est courant que les équations soient mises à l'échelle afin qu'elles aient des valeurs normalisées. Par exemple, dans des cas extrêmes, la densité électronique dans les semi-conducteurs peut varier sur 18 ordres de grandeur, et le champ électrique peut changer de manière galbée, sur 6 (ou plus) ordres de grandeur.

Cependant, les journaux ne donnent jamais vraiment de raison de le faire. Personnellement, je suis content de traiter des équations en unités réelles, y a-t-il un avantage numérique à le faire, est-ce impossible autrement? J'ai pensé avec une double précision qu'il y aurait suffisamment de chiffres pour faire face à ces fluctuations.

---

Les deux réponses sont très utiles, merci beaucoup!

La résolution d'un PDE (linéaire) consiste à discrétiser l'équation pour donner un système linéaire, qui est ensuite résolu par un solveur linéaire dont la convergence (taux) dépend du nombre de condition de la matrice. La mise à l'échelle des variables réduit souvent ce nombre de conditions, améliorant ainsi la convergence. (Cela revient essentiellement à appliquer un préconditionneur diagonal, voir Nicholas Higham's Accuracy and Stability of Numerical Algorithms .)

La résolution des EDP non linéaires nécessite en outre une méthode pour résoudre des équations non linéaires telles que la méthode de Newton, où la mise à l'échelle peut également influencer la convergence.

Étant donné que tout normaliser nécessite généralement très peu d'efforts, c'est presque toujours une bonne idée.

$$-\varepsilon \Delta u + u = 0 \text{ on } \Omega,\\ u = 1 \text{ on } \partial\Omega.$$

Cela dit, il n'y a pas de mise à l'échelle des variables ou des domaines qui supprime cette difficulté.

$u_{\alpha}$

$$-\alpha^2\Delta u = f_\alpha \text{ on } \alpha\Omega$$ $\alpha$$u_1$ $$-\Delta u = f \text{ on } \Omega.$$ $u_{\alpha}(\mathbf{x}):=u_1(\mathbf{x}/\alpha)$$\alpha$ . Cela montre que la mise à l'échelle du paramètre$\alpha$n'influence pas vraiment le comportement de la solution - un fait qui est assez évident ici, mais qui peut être assez difficile à voir dans d'autres cas. Vous pouvez construire des exemples assez complexes avec plus de paramètres de la même manière, par exemple, Navier-Stokes et sa non-dimensionnalisation dans la formulation du nombre de Reynolds.

Le traitement des nombres à virgule flottante peut être délicat en ce qui concerne la soustraction de petits nombres de plus grands nombres, ainsi que de nombreux autres aspects. Je recommanderais de lire les articles de blog de John D. Cooks à leur sujet, tels que

Anatomie d'un nombre à virgule flottante

ainsi que d'Oracle

Ce que tout informaticien devrait savoir sur l'arithématique en virgule flottante

Certains algorithmes numériques pour la minimisation ou la maximisation nécessitent également une normalisation pour la stabilité numérique.