📚 الرئيسية

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

مسابقة تخصص · Analyse Fonctionnelle · المعامل: 2 · المدة: 2سا

Concours d'accès au doctorat 3ème cycle LMD 2020/2021 — Spécialité Equations Différentielles et Applications — Épreuve : Equations Différentielles et Applications (Variante 1), Université Echahid Hamma Lakhdar d'El Oued, Faculté des Sciences Exactes, Département des Mathématiques — Date 6/03/2021 — Coefficient 02 — Durée 2h.

التمرين 1

Exercice 1 — Problème biharmonique et théorème de Lax-Milgram

#pde#lax-milgram#sobolev-spaces#biharmonic#variational-formulation

Soit ΩRN\Omega \subset \mathbb{R}^N un ouvert (non vide) borné et assez régulier. Étant donné fL2(Ω)f \in L^2(\Omega), on considère le problème

Trouver uH2(Ω)u \in H^2(\Omega), tel que

(Pvar)i,j=1NΩ2u(x)xi22v(x)xj2dx=Ωf(x)v(x)dx,vH2(Ω).(P_{var}) \quad \sum_{i,j=1}^N \int_\Omega \frac{\partial^2 u(x)}{\partial x_i^2} \frac{\partial^2 v(x)}{\partial x_j^2} dx = \int_\Omega f(x) v(x) dx, \quad \forall v \in H^2(\Omega).
  1. (3 pts) Quelles sont les hypothèses du théorème de Lax-Milgram vérifiées par le problème (Pvar)(P_{var}) ?
  2. (3 pts) Peut-on utiliser le théorème de Lax-Milgram pour montrer l'existence et l'unicité de la solution du problème (Pvar)(P_{var}) ? (justifier votre réponse).
الحل

1.

H2(Ω)H^2(\Omega) est un espace de Hilbert. La forme bilinéaire a(u,v)=i,jΩ2uxi22vxj2dxa(u,v) = \sum_{i,j} \int_\Omega \frac{\partial^2 u}{\partial x_i^2} \frac{\partial^2 v}{\partial x_j^2} dx est bien une forme bilinéaire continue sur H2(Ω)H^2(\Omega) (par inégalité de Hölder). La forme linéaire L(v)=Ωf(x)v(x)dxL(v) = \int_\Omega f(x)v(x)dx est continue sur H2(Ω)H^2(\Omega) (par Hölder).

2.

Non, on ne peut pas utiliser le théorème de Lax-Milgram car a(,)a(\cdot,\cdot) n'est pas coercive sur H2(Ω)H^2(\Omega). En effet, pour toute fonction vv constante et non nulle, vH2(Ω)v \in H^2(\Omega) mais a(v,v)=0a(v,v) = 0 car toutes les dérivées secondes sont nulles. Donc il n'existe pas α>0\alpha \gt 0 tel que a(v,v)αvH22a(v,v) \geq \alpha \|v\|_{H^2}^2. La coercivité est une condition nécessaire pour appliquer Lax-Milgram.

Non, car a(,) n’est pas coercive sur H2(Ω)\boxed{\text{Non, car } a(\cdot,\cdot) \text{ n'est pas coercive sur } H^2(\Omega)}

التمرين 2

Exercice 2 — Équation intégrale de Volterra : point fixe de Banach

#integral-equations#fixed-point#banach-contraction#volterra

Soient E=C([a,b],R)E = C([a,b], \mathbb{R}) muni de la norme de la convergence uniforme E\|\cdot\|_E, K:[a,b]×[a,b]RK : [a,b] \times [a,b] \to \mathbb{R} une fonction continue, φE\varphi \in E et λR\lambda \in \mathbb{R}. On considère l'opérateur A:EEA : E \to E défini par

Au(t)=φ(t)+λatK(t,s)u(s)ds.Au(t) = \varphi(t) + \lambda \int_a^t K(t,s)u(s)ds.
  1. (3 pts) Montrer par récurrence sur nn que M>0\exists M \gt 0 :
AnuAnvE(λM(ba))nn!uvEu,vE.\|A^n u - A^n v\|_E \leq \frac{(|\lambda| M(b-a))^n}{n!} \|u - v\|_E \quad \forall u, v \in E.
  1. (3 pts) Déduire que pour tout réel λ\lambda, l'équation u(t)=φ(t)+λatK(t,s)u(s)dsu(t) = \varphi(t) + \lambda \int_a^t K(t,s)u(s)ds admet une unique solution dans EE.
الحل

1.

Posons M=maxat,sbK(t,s)M = \max_{a \leq t,s \leq b} |K(t,s)|.

Pour n=1n = 1 : Au(t)Av(t)=λatK(t,s)(u(s)v(s))dsλM(ta)uvEλM(ba)uvE|Au(t) - Av(t)| = |\lambda| \left|\int_a^t K(t,s)(u(s)-v(s))ds\right| \leq |\lambda| M (t-a) \|u-v\|_E \leq |\lambda| M (b-a) \|u-v\|_E.

Si vrai au rang nn, alors An+1u(t)An+1v(t)λMat(λM(sa))nn!uvEds=(λM(ta))n+1(n+1)!uvE|A^{n+1}u(t) - A^{n+1}v(t)| \leq |\lambda| M \int_a^t \frac{(|\lambda|M(s-a))^n}{n!} \|u-v\|_E ds = \frac{(|\lambda|M(t-a))^{n+1}}{(n+1)!}\|u-v\|_E.

Donc An+1uAn+1vE(λM(ba))n+1(n+1)!uvE\|A^{n+1}u - A^{n+1}v\|_E \leq \frac{(|\lambda|M(b-a))^{n+1}}{(n+1)!}\|u-v\|_E.

2.

Pour tout λR\lambda \in \mathbb{R}, on peut choisir nn assez grand tel que (λM(ba))nn!<1\frac{(|\lambda|M(b-a))^n}{n!} \lt 1. Alors AnA^n est contractant. Par suite AA admet un unique point fixe dans EE (principe de contraction de Banach), qui est la solution de l'équation intégrale.

u(t)=φ(t)+λatK(t,s)u(s)ds admet une unique solution\boxed{u(t) = \varphi(t) + \lambda \int_a^t K(t,s)u(s)ds \text{ admet une unique solution}}

التمرين 3

Exercice 3 — Équation de Van der Pol modifiée : Lindstedt et moyennisation

#ode#perturbation-methods#lindstedt#averaging-method#van-der-pol

Considérons l'équation différentielle non linéaire suivante :

d2x(t,ε)dt2+ε(x2(t,ε)1)dx(t,ε)dt+(1+aε)x(t,ε)=0,x(0,ε)=A,dx(0,ε)dt=0,(1)\frac{d^2x(t,\varepsilon)}{dt^2} + \varepsilon(x^2(t,\varepsilon) - 1)\frac{dx(t,\varepsilon)}{dt} + (1 + a\varepsilon)x(t,\varepsilon) = 0, \quad x(0,\varepsilon) = A, \quad \frac{dx(0,\varepsilon)}{dt} = 0, \quad (1)

0<ε10 \lt \varepsilon \ll 1, et λ\lambda, aa deux constantes réelles.

  1. (4 pts) Trouver la solution approximative de Lindstedt du premier ordre de l'équation (1).
  2. (4 pts) Appliquer la méthode de la moyenne pour trouver une solution approximative de l'équation (1).
الحل

1. Méthode de Lindstedt

On réécrit : x¨+x=ε((x21)x˙ax)F=(x21)x˙ax\ddot{x} + x = \varepsilon(-(x^2-1)\dot{x} - ax) \Rightarrow F = -(x^2-1)\dot{x} - ax.

On pose y(θ,ε)=y0(θ)+εy1(θ)y(\theta,\varepsilon) = y_0(\theta) + \varepsilon y_1(\theta) et ω(ε)=1+εω1+\omega(\varepsilon) = 1 + \varepsilon\omega_1 + \ldots

Ordre 0 : y¨0+y0=0\ddot{y}_0 + y_0 = 0, solution y0(θ)=Acosθy_0(\theta) = A\cos\theta.

Ordre 1 : y¨1+y1=2ω1y¨0y˙02\ddot{y}_1 + y_1 = -2\omega_1 \ddot{y}_0 - \dot{y}_0^2... Après calcul et élimination des termes séculaires (condition ω1=a/2\omega_1 = -a/2, A=2A = 2) :

y1(θ)=34sinθ14sin3θy_1(\theta) = \frac{3}{4}\sin\theta - \frac{1}{4}\sin 3\theta y(θ,ε)=2cosθ+ε(34sinθ14sin3θ),ω(ε)=1εa2+O(ε)\boxed{y(\theta,\varepsilon) = 2\cos\theta + \varepsilon\left(\frac{3}{4}\sin\theta - \frac{1}{4}\sin 3\theta\right), \quad \omega(\varepsilon) = 1 - \varepsilon\frac{a}{2} + O(\varepsilon)}

2. Méthode de la moyenne

Avec ω=1\omega = 1, F=(x21)x˙axF = -(x^2-1)\dot{x} - ax. On pose x=Acosθx = A\cos\theta, x˙=Asinθ\dot{x} = -A\sin\theta.

Les équations moyennées donnent A˙=εA8(4A2)\dot{A} = \frac{\varepsilon A}{8}(4 - A^2) et Φ˙=εa2\dot{\Phi} = \frac{\varepsilon a}{2}.

On résout A˙=εA8(4A2)\dot{A} = \frac{\varepsilon A}{8}(4-A^2) par séparation de variables : A(t)=2[(4/A021)eεt+1]1/2A(t) = \frac{2}{[(4/A_0^2-1)e^{-\varepsilon t}+1]^{1/2}}.

A(t)2A(t) \to 2 quand t+t \to +\infty (indépendamment de A0A_0). Avec Φ(t)=εat2+π2\Phi(t) = \frac{\varepsilon a t}{2} + \frac{\pi}{2} :

x(t)=2[(4/A021)eεt+1]1/2sin(t+εat2+π2)\boxed{x(t) = \frac{2}{[(4/A_0^2-1)e^{-\varepsilon t}+1]^{1/2}} \sin\left(t + \frac{\varepsilon a t}{2} + \frac{\pi}{2}\right)}