📚 الرئيسية

مسابقة دكتوراه 2021Université Echahid Hamma Lakhdar d'El Oued

مسابقة تخصص · EDP · المدة: 3سا

MathématiquesAppliques

التمرين 1

تمرين 1

On s'intéresse à la discrétisation du problème suivant :

(P):{ut2ux2u(x,t)=0,(x,t)[0,1]×R+,u(0,t)=u(1,t)=0,tR+,u(x,0)=u0(x),x[0,1].(\mathcal{P}) : \begin{cases} \dfrac{\partial u}{\partial t} - \dfrac{\partial^2 u}{\partial x^2} - u(x,t) = 0, & (x,t) \in [0,1] \times \mathbb{R}_+, \\ u(0,t) = u(1,t) = 0, & t \in \mathbb{R}^*_+, \\ u(x,0) = u_0(x), & x \in [0,1]. \end{cases}

Pour :

xi=ih,h=1N+1,i=0,1,,N+1,tn=nk, nN.x_i = ih, \quad h = \frac{1}{N+1}, \quad i = 0,1,\ldots,N+1, \qquad t_n = nk,\ n \in \mathbb{N}.

On considère le schéma suivant :

{uin+1uinkui+1n+12uin+1+ui1n+1h2uin+1=0,i=1,,N, nN,u0n+1=uN+1n+1=0,nN,ui0=u0(xi),i=1,,N.\begin{cases} \dfrac{u_i^{n+1} - u_i^n}{k} - \dfrac{u_{i+1}^{n+1} - 2u_i^{n+1} + u_{i-1}^{n+1}}{h^2} - u_i^{n+1} = 0, & i = 1,\ldots,N,\ n \in \mathbb{N}, \\[6pt] u_0^{n+1} = u_{N+1}^{n+1} = 0, & n \in \mathbb{N}, \\[4pt] u_i^0 = u_0(x_i), & i = 1,\ldots,N. \end{cases}
  1. Montrer que le schéma demande, à chaque pas de temps, la résolution du système linéaire :
AUn+1=b(1)AU^{n+1} = b \tag{1}

ARN,NA \in \mathbb{R}^{N,N}, bRNb \in \mathbb{R}^N à déterminer.

  1. Prouver que le système (1)(1) admet une solution unique.

  2. Pour nNn \in \mathbb{N}, on pose Un=supi=1,,Nuin\|U^n\|_\infty = \displaystyle\sup_{i=1,\ldots,N} |u_i^n|.

    3.1. Montrer que :

11+k1+βk,\frac{1}{1+k} \leq 1 + \beta k,

pour β=11α\beta = \dfrac{1}{1-\alpha}, k]0,α[k \in ]0,\alpha[ et α]0,1[\alpha \in ]0,1[.

3.2. Sous quelle condition sur kk et hh le schéma est-il stable ? En déduire que :

UneβTU0\|U^n\|_\infty \leq e^{\beta T} \|U^0\|_\infty

التمرين 2

تمرين 2

La répartition des contraintes est donnée par le tenseur suivant dans le repère (O,e1,e2,e3)(O, e_1, e_2, e_3) :

σ(x1,x2,x3)=[x1+x2σ12(x1,x2)0σ12(x1,x2)x12x2300x2]\sigma(x_1, x_2, x_3) = \begin{bmatrix} x_1 + x_2 & \sigma_{12}(x_1,x_2) & 0 \\ \sigma_{12}(x_1,x_2) & x_1 - 2x_2 & 3 \\ 0 & 0 & x_2 \end{bmatrix}

La contrainte agissant au point M(0,1)M(0,1) sur un plan vertical de normale inclinée de 45°45° par rapport à l'axe Ox1Ox_1 est une contrainte de cisaillement pur τ\tau.

Déterminer σ12(x1,x2)\sigma_{12}(x_1, x_2) et donner la valeur de τ\tau.

التمرين 3

تمرين 3

Soit Ω\Omega un ouvert borné régulier de classe C1C^1, fL2(Ω)f \in L^2(\Omega) et gg une trace sur Ω\partial\Omega d'une fonction de H1(Ω)H^1(\Omega). On considère le problème :

{Δu=f,dans Ω,u+un=g,sur Ω.(2)\begin{cases} -\Delta u = f, & \text{dans } \Omega, \\ u + \dfrac{\partial u}{\partial n} = g, & \text{sur } \partial\Omega. \end{cases} \tag{2}
  1. Montrer que toute solution du problème (2)(2) est une solution d'un problème de la forme :
{Trouver uV,a(u,v)=L(v),vV(3)\begin{cases} \text{Trouver } u \in V, \\ a(u,v) = L(v), \quad \forall v \in V \end{cases} \tag{3}

VV est un espace de Hilbert, aa et LL sont des formes bilinéaire et linéaire respectivement.

  1. Montrer qu'il existe une solution unique du problème (3)(3).

  2. Déduire une autre formulation variationnelle du problème (2)(2) à la forme d'optimisation :

{uV,J(u)=minvVJ(v).\begin{cases} u \in V, \\ J(u) = \displaystyle\min_{v \in V} J(v). \end{cases}