📚 الرئيسية

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

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

Mathématiques Appliques

التمرين 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=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 = 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 :

(1+k)n(1+k)T/k(1+k)^n \leq (1+k)^{T/k}

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

UneTU0\|U^n\|_\infty \leq e^T \|U^0\|_\infty

التمرين 2

تمرين 2

L'état des contraintes en un point d'un solide est donné par :

σ=[133313331]\sigma = \begin{bmatrix} 1 & 3 & 3 \\ 3 & 1 & 3 \\ 3 & 3 & 1 \end{bmatrix}
  1. Calculer les contraintes principales (σ1σ2σ3)(\sigma_1 \geq \sigma_2 \geq \sigma_3) du tenseur σ\sigma.

  2. Déterminer le vecteur normal n\vec{n} au plan PP qui coupe les axes aux points A(1,0,0)A(1,0,0), B(0,1,0)B(0,1,0) et C(0,0,2)C(0,0,2).

  3. Calculer les composantes normale σn\sigma_n et tangentielle σr\sigma_r du vecteur contrainte agissant sur le plan PP.

  4. Déterminer les directions principales du tenseur σ\sigma.

التمرين 3

تمرين 3

Soit Ω\Omega un ouvert quelconque de RN\mathbb{R}^N et fL2(Ω)f \in L^2(\Omega). On considère le problème :

{Δu+u=f,dans Ω,u=0,sur Ω.(2)\begin{cases} \Delta u + u = f, & \text{dans } \Omega, \\ u = 0, & \text{sur } \partial\Omega. \end{cases} \tag{2}
  1. Trouver la formulation variationnelle du problème (2)(2) à la forme :
{uH01(Ω),a(u,v)=L(v),vH01(Ω).(3)\begin{cases} u \in H^1_0(\Omega), \\ a(u,v) = L(v), \quad \forall v \in H^1_0(\Omega). \end{cases} \tag{3}
  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}