A note on the method of weighted differences (Q1097653)

On introduit une nouvelle approximation avec des différences pour le problème singulier avec de perturbation et on démontre une condition nécessaire de convergence uniforme. On considère le problème initial \[ (1)\quad L^ 1u_{\epsilon}\equiv \epsilon u'(x)+a(x)u(x)=f(x),\quad (x>0),\quad u(0)=\phi, \] où \(\epsilon\) est un petit paramètre, a(x) et f(x) sont des fonctions unies de la variable x; on suppose que \(a(x)\geq \alpha >0\). Pour approcher la valeur de \(u(x_ i)\) on utilise l'expression \((b_ i+c_ i)u_{i+1}+(1-2b_ i)u_ i+(b_ i-c_ i)u_{i-1},\) où \(b_ i\) et \(c_ i\) sont des paramètres qui doivent déterminer. On obtient la suivante équation aux différences \[ (2)\quad \epsilon \frac{(u_{i+1}-u_{i- 1})}{2h}+a_ i\{(b_ i+c_ i)u_{i+1}+(1-2b_ i)u_ i+(b_ i-c_ i)u_{i-1}\}=f_ i. \] On fait l'observation que (2) est compatible avec (1). L'équation aux différences (2) devient \(L_ h^ 1u_ i\equiv 2a_ ib_ iu_{i+1}+a_ i(1-2b_ i)u_ i=f_ i\). On démontre le théorème 1: Supposons que la solution du problème \(L^ 1_ hu_ i=f_ i\) \((i>0)\) \(u_ 0=\phi\), converge uniformement dans \(\epsilon\) vers la solution de l'équation (1). Soient \(\rho =h/\epsilon\) et n, un nombre entier non-négatif, fixés. Dans ces conditions on détermine lim \(b_ n\), lorsque \(h\to 0\). On considère aussi le problème avec des valeurs aux limites \[ (3)\quad L^ 2u_{\epsilon}\equiv \epsilon u''(x)+a(x)u'(x)=f(x),\quad (0<x<1),\quad u(0)=A,\quad u(1)=B,\quad a(x)\geq \alpha >0. \] Pour approcher \(u'(x_ i)\) on utilise une expression, où \(d_ i\) est un paramètre qui doit déterminé. Dans ce cas on obtient la suivante équation aux différences \[ (4)\quad L^ 2_ hu_ i\equiv \epsilon \frac{u_{i+1}-2u_ i+u_{i-2}}{h^ 2}+a_ i\frac{(1+d_ i)(u_{i+1}-u_ i)+(1-d_ i)(u_ i-u_{i-1})}{2h}=f_ i, \] \((i=1,2,...,N-1).\) L'A. fait l'observation que (4) est compatible avec (3). On démontre le théorème 2: Supposons que la solution pour \(L^ 2_ hu_ i=f_ i\quad (i>0,\quad i=1(1)N-1)\quad u_ 0=A,\quad u_ N=B,\) converge uniformément dans \(\epsilon\) vers la solution de (3). Soient \(\rho =h/\epsilon\) et n, un nombre entier non-négatif, fixés. Dans ces conditions on détermine lim \(d_ n\), lorsque \(h\to 0\). Dans la dernière partie du travail on considère le problème elliptique- parabolique singulier avec de perturbation de la forme \[ L\quad 3u_{\epsilon}\equiv \epsilon (\partial^ 2u/\partial y^ 2)+(\partial^ 2u/\partial x^ 2)-a(x,y)(\partial u/\partial y)=f(x,y),\quad (0<x,\quad y<1)\quad u|_{\Gamma}=0, \] où \(\epsilon >0\) est un petit paramètre, \(a(x,y)\geq \alpha >0\). On fait l'observation que si l'on prend \(\epsilon_{i,j}\), dans la manière convenable, alors on obtient le même schéma aux différences comme dans le cas de la méthode d'Ilin.