الرئيسية

مسابقة دكتوراه 2015Université Kasdi Merbah - Ouargla — الموضوع 01

مسابقة تخصص · Analyse Numérique & Optimisation

FB_IMG_1527867318026.pdf, concours du 17 octobre 2015, Analyse Numérique, variante 2

التمرين 1

Différences finies, stabilité et convergence d'un problème aux limites

#différences finies#stabilité#convergence#matrice monotone

Soit ff une fonction régulière. On souhaite approcher la solution du problème

{u(x)+u(x)1+x=f(x),0<x<1,u(0)=u(1)=0,\begin{cases} -u''(x)+\dfrac{u'(x)}{1+x}=f(x),&0<x<1,\\ u(0)=u(1)=0, \end{cases}

par la méthode des différences finies. On considère une discrétisation (x0,,xn+1)(x_0,\ldots,x_{n+1}) de [0,1][0,1] à pas constant h=1n+1h=\frac1{n+1} et les formules centrées pour uu' et uu''. On note uiu_i la valeur approchée de u(xi)u(x_i).

  1. Écrire le problème approché sous la forme Au=bA\vec u=\vec b, en explicitant la matrice A=(aij)1i,jnA=(a_{ij})_{1\le i,j\le n} et le vecteur b=(bi)1in\vec b=(b_i)_{1\le i\le n}.

On admet que AA est monotone, c'est-à-dire que A1=(cij)A^{-1}=(c_{ij}) existe et cij0c_{ij}\ge0.

  1. On considère
ϕ(x)=23(x2+2x)ln2(1+x)22ln(1+x),x[0,1].\phi(x)=\frac23(x^2+2x)\ln 2-\frac{(1+x)^2}{2}\ln(1+x),\qquad x\in[0,1].

a. Calculer ϕ(0)\phi(0), ϕ(1)\phi(1) et

ϕ(x)+ϕ(x)1+x.-\phi''(x)+\frac{\phi'(x)}{1+x}.

b. Montrer qu'il existe une constante CC, indépendante de hh, telle que

max1inϕ(xi+1)2ϕ(xi)+ϕ(xi1)h2+ϕ(xi+1)ϕ(xi1)2h(1+ih)1Ch2.\max_{1\le i\le n}\left| -\frac{\phi(x_{i+1})-2\phi(x_i)+\phi(x_{i-1})}{h^2} +\frac{\phi(x_{i+1})-\phi(x_{i-1})}{2h(1+ih)}-1 \right|\le Ch^2.

c. En posant ϕ=(ϕ(xi))1in\vec\phi=(\phi(x_i))_{1\le i\le n}, déduire que

min1in(Aϕ)i1Ch2.\min_{1\le i\le n}(A\vec\phi)_i\ge1-Ch^2.
  1. Montrer qu'il existe une constante MM, indépendante de hh, telle que
A1M,A1=max1inj=1ncij.\|A^{-1}\|_\infty\le M, \qquad \|A^{-1}\|_\infty=\max_{1\le i\le n}\sum_{j=1}^n|c_{ij}|.

التمرين 2

Problème de transmission et approximation de Galerkin

#éléments finis#Galerkin#Sobolev#problème de transmission

Soient I=(1,1)I=(-1,1), I=(1,0)I_-=(-1,0), I+=(0,1)I_+=(0,1) et fL2(I)f\in L^2(I). On note V=H01(I)V=H_0^1(I).

Pour tout λ>0\lambda>0, on considère le problème de transmission

{λ(uλ)=fdans (1,0),(uλ+)=fdans (0,1),uλ+(0)=uλ(0),duλ+dx(0)=λduλdx(0),uλ+(1)=uλ(1)=0.\begin{cases} -\lambda(u_\lambda^-)''=f&\text{dans }(-1,0),\\ -(u_\lambda^+)''=f&\text{dans }(0,1),\\ u_\lambda^+(0)=u_\lambda^-(0),\\ \dfrac{du_\lambda^+}{dx}(0)=\lambda\dfrac{du_\lambda^-}{dx}(0),\\ u_\lambda^+(1)=u_\lambda^-(-1)=0. \end{cases}
  1. Trouver le problème variationnel correspondant.
  2. Montrer que le problème variationnel est bien posé.
  3. Montrer l'estimation de Céa
uλuhH1(I)CinfvhVhuλvhH1(I),\|u_\lambda-u_h\|_{H^1(I)} \le C\inf_{v_h\in V_h}\|u_\lambda-v_h\|_{H^1(I)},

CC est indépendante de hh et uhu_h est la solution approchée obtenue par la méthode de Galerkin P1P_1. 4. Montrer que si λ1\lambda\ne1, alors uλH2(I)u_\lambda\notin H^2(I). 5. Est-il raisonnable d'espérer une estimation de la forme

uλuhH1(I)Ch,\|u_\lambda-u_h\|_{H^1(I)}\le C'h,

CC' est indépendante de hh ?