A note on the method of weighted differences

Linear ordinary differential equations and systems (34A30) Linear boundary value problems for ordinary differential equations (34B05) Singular perturbations for ordinary differential equations (34E15) Singular perturbations in context of PDEs (35B25) Partial differential equations of mixed type and mixed-type systems of partial differential equations (35M99) Numerical methods for initial value problems involving ordinary differential equations (65L05) Numerical solution of boundary value problems involving ordinary differential equations (65L10) Finite difference methods for boundary value problems involving PDEs (65N06)

Review text

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.

Recommendations

Cites work

Cited in

(5)

This page was built for publication: A note on the method of weighted differences

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1097653)