📚 الرئيسية

مسابقة دكتوراه 2018Université Mohammed Seddik Benyahia - Jijel — الموضوع 01

مسابقة تخصص · EDP

Concours de Doctorat LMD 2018/2019, Option : EDP et Applications, Épreuve : Analyse variationnelle des EDP, Université de Jijel, Faculté des Sciences Exactes et Informatique, Département de Mathématiques.

التمرين 1

Problème elliptique mixte Dirichlet-Robin : formulation variationnelle

#elliptic-pde#mixed-boundary-conditions#robin-condition#variational-formulation#lax-milgram#trace-theorem

Soit Ω\Omega un ouvert borné connexe de RN\mathbb{R}^N de frontière Ω\partial\Omega assez régulière. On suppose que Ω=Γ0Γ1\partial\Omega=\overline{\Gamma_0}\cup\overline{\Gamma_1} avec Γ0Γ1=\Gamma_0\cap\Gamma_1=\emptyset, et que la mesure (N1)(N-1)-dimensionnelle de Γ0\Gamma_0 est non nulle. On considère le problème (1){div(p(x)u(x))+q(x)u(x)=f(x),xΩ,u(x)=0,xΓ0,p(x)uν(x)+σu(x)=g1(x),xΓ1,(1)\quad\begin{cases}-\operatorname{div}\big(p(x)\nabla u(x)\big)+q(x)u(x)=f(x),& x\in\Omega,\\ u(x)=0,& x\in\Gamma_0,\\ p(x)\dfrac{\partial u}{\partial\nu}(x)+\sigma u(x)=g_1(x),& x\in\Gamma_1,\end{cases}fL2(Ω)f\in L^2(\Omega), pL(Ω)p\in L^\infty(\Omega) avec p(x)α>0p(x)\ge\alpha>0 p.p., qL(Ω)q\in L^\infty(\Omega), q0q\ge0, σ0\sigma\ge0, g1L2(Γ1)g_1\in L^2(\Gamma_1) et ν\nu la normale unitaire extérieure.

  1. Donner une formulation variationnelle du problème aux limites (1).
  2. Montrer que cette formulation admet une solution unique.
  3. Montrer que la solution variationnelle est solution du problème aux limites (1).
الحل

Espace de travail

V={vH1(Ω): v=0 sur Γ0},V=\{v\in H^1(\Omega):\ v=0\ \text{sur }\Gamma_0\}, sous-espace fermé de H1(Ω)H^1(\Omega). Comme Γ0>0|\Gamma_0|>0, l'inégalité de Poincaré vaut sur VV : vL2CvL2\|v\|_{L^2}\le C\|\nabla v\|_{L^2}, donc vL2\|\nabla v\|_{L^2} est une norme équivalente sur VV.

1. Formulation variationnelle

On multiplie par vVv\in V, on intègre et on applique la formule de Green : ΩpuvΩpuνv+Ωquv=Ωfv.\int_\Omega p\,\nabla u\cdot\nabla v-\int_{\partial\Omega}p\frac{\partial u}{\partial\nu}v+\int_\Omega quv=\int_\Omega fv. Sur Γ0\Gamma_0, v=0v=0 ; sur Γ1\Gamma_1, puν=g1σup\frac{\partial u}{\partial\nu}=g_1-\sigma u. D'où : trouver uVu\in V tel que pour tout vVv\in V, a(u,v)=Ω ⁣(puv+quv)dx+Γ1 ⁣σuvdσ=Ωfvdx+Γ1g1vdσ=L(v).a(u,v)=\int_\Omega\!\big(p\,\nabla u\cdot\nabla v+quv\big)dx+\int_{\Gamma_1}\!\sigma uv\,d\sigma=\int_\Omega fv\,dx+\int_{\Gamma_1}g_1v\,d\sigma=L(v).

2. Existence et unicité (Lax-Milgram)

Continuité de aa : p,qLp,q\in L^\infty, le terme de bord se contrôle par le théorème de trace vL2(Γ1)CvH1\|v\|_{L^2(\Gamma_1)}\le C\|v\|_{H^1} ; donc a(u,v)CuVvV|a(u,v)|\le C\|u\|_V\|v\|_V. Coercivité : avec pαp\ge\alpha, q0q\ge0, σ0\sigma\ge0 : a(v,v)αΩv2αCPvV2.a(v,v)\ge\alpha\int_\Omega|\nabla v|^2\ge\alpha C_P\|v\|_V^2. Continuité de LL : L(v)fL2vL2+g1L2(Γ1)vL2(Γ1)CvV|L(v)|\le\|f\|_{L^2}\|v\|_{L^2}+\|g_1\|_{L^2(\Gamma_1)}\|v\|_{L^2(\Gamma_1)}\le C\|v\|_V. Par Lax-Milgram, il existe une unique uVu\in V solution de la formulation variationnelle.

3. Retour au problème fort

En prenant v=φD(Ω)v=\varphi\in\mathcal{D}(\Omega) (nulle au bord), la formulation donne div(pu)+quf,φ=0\langle-\operatorname{div}(p\nabla u)+qu-f,\varphi\rangle=0, donc div(pu)+qu=f dans Ω (au sens D).-\operatorname{div}(p\nabla u)+qu=f\ \text{dans }\Omega\ (\text{au sens }\mathcal{D}'). En supposant uu assez régulière, on réinjecte dans la formulation et la formule de Green : le terme de volume s'annule, il reste Γ1(puν+σug1)vdσ=0\int_{\Gamma_1}\big(p\frac{\partial u}{\partial\nu}+\sigma u-g_1\big)v\,d\sigma=0 pour tout vVv\in V, d'où la condition de Robin puν+σu=g1p\frac{\partial u}{\partial\nu}+\sigma u=g_1 sur Γ1\Gamma_1. La condition u=0u=0 sur Γ0\Gamma_0 vient de uVu\in V. Donc uu résout (1).

التمرين 2

Problème variationnel (P_λ) avec estimation uniforme en λ

#variational-problem#lax-milgram#a-priori-estimate#uniform-bound#neumann-type#coercivity

Soit Ω\Omega un ouvert borné connexe de R2\mathbb{R}^2 à frontière régulière et soit fL2(Ω)f\in L^2(\Omega). On considère le problème variationnel (Pλ){Chercher uH1(Ω) tel queΩuvdx+λ(Ωudx)(Ωvdx)=ΩfvdxvH1(Ω),(P_\lambda)\quad\begin{cases}\text{Chercher }u\in H^1(\Omega)\text{ tel que}\\ \displaystyle\int_\Omega\nabla u\cdot\nabla v\,dx+\lambda\Big(\int_\Omega u\,dx\Big)\Big(\int_\Omega v\,dx\Big)=\int_\Omega fv\,dx\quad\forall v\in H^1(\Omega),\end{cases}λ1\lambda\ge1. On admet l'inégalité suivante : vH1(Ω)2C1(vL2(Ω)2+Ωvdx2)vH1(Ω).\|v\|_{H^1(\Omega)}^2\le C_1\Big(\|\nabla v\|_{L^2(\Omega)}^2+\Big|\int_\Omega v\,dx\Big|^2\Big)\qquad\forall v\in H^1(\Omega).

  1. Montrer que (Pλ)(P_\lambda) admet une solution unique uu et que cette solution vérifie uH1(Ω)C2fL2(Ω)\|u\|_{H^1(\Omega)}\le C_2\|f\|_{L^2(\Omega)}, où C2>0C_2>0 est indépendante de λ\lambda.
  2. Montrer qu'il existe une constante C3>0C_3>0, indépendante de λ\lambda, telle que ΩudxC3λfL2(Ω)\Big|\int_\Omega u\,dx\Big|\le\dfrac{C_3}{\lambda}\|f\|_{L^2(\Omega)}.
الحل

1. Existence, unicité et estimation uniforme

La forme bilinéaire est aλ(u,v)=Ωuv+λ(Ωu)(Ωv)a_\lambda(u,v)=\int_\Omega\nabla u\cdot\nabla v+\lambda\big(\int_\Omega u\big)\big(\int_\Omega v\big).

Continuité : aλ(u,v)uv+λΩuL2vL2|a_\lambda(u,v)|\le\|\nabla u\|\|\nabla v\|+\lambda|\Omega|\,\|u\|_{L^2}\|v\|_{L^2} (Cauchy-Schwarz sur u\int u), donc continue sur H1H^1 (pour λ\lambda fixé).

Coercivité (uniforme en λ1\lambda\ge1) : avec l'inégalité admise et λ1\lambda\ge1, aλ(v,v)=vL22+λΩv2vL22+Ωv21C1vH12.a_\lambda(v,v)=\|\nabla v\|_{L^2}^2+\lambda\Big|\int_\Omega v\Big|^2\ge\|\nabla v\|_{L^2}^2+\Big|\int_\Omega v\Big|^2\ge\frac1{C_1}\|v\|_{H^1}^2. La constante de coercivité 1/C11/C_1 ne dépend pas de λ\lambda. Par Lax-Milgram, (Pλ)(P_\lambda) a une unique solution uH1(Ω)u\in H^1(\Omega).

Estimation : prenons v=uv=u. Alors aλ(u,u)=ΩfufL2uL2fL2uH1a_\lambda(u,u)=\int_\Omega fu\le\|f\|_{L^2}\|u\|_{L^2}\le\|f\|_{L^2}\|u\|_{H^1}. Avec la coercivité, 1C1uH12fL2uH1  uH1C1fL2=:C2fL2,\frac1{C_1}\|u\|_{H^1}^2\le\|f\|_{L^2}\|u\|_{H^1}\ \Rightarrow\ \|u\|_{H^1}\le C_1\|f\|_{L^2}=:C_2\|f\|_{L^2}, avec C2=C1C_2=C_1 indépendante de λ\lambda.

2. Décroissance en 1/λ1/\lambda de la moyenne

Prenons la fonction test constante v1H1(Ω)v\equiv1\in H^1(\Omega). Alors v=0\nabla v=0 et Ωv=Ω\int_\Omega v=|\Omega|, donc la formulation donne λ(Ωu)Ω=Ωf1dx=Ωf.\lambda\Big(\int_\Omega u\Big)|\Omega|=\int_\Omega f\cdot1\,dx=\int_\Omega f. D'où Ωu=1λΩΩf1λΩfL21L2=Ω1/2λΩfL2=1λΩ1/2fL2.\Big|\int_\Omega u\Big|=\frac{1}{\lambda|\Omega|}\Big|\int_\Omega f\Big|\le\frac{1}{\lambda|\Omega|}\|f\|_{L^2}\|1\|_{L^2}=\frac{|\Omega|^{1/2}}{\lambda|\Omega|}\|f\|_{L^2}=\frac{1}{\lambda|\Omega|^{1/2}}\|f\|_{L^2}. En posant C3=Ω1/2C_3=|\Omega|^{-1/2} (indépendante de λ\lambda) : ΩudxC3λfL2(Ω).\boxed{\Big|\int_\Omega u\,dx\Big|\le\frac{C_3}{\lambda}\|f\|_{L^2(\Omega)}.} La moyenne de uu tend vers 00 quand λ\lambda\to\infty : à la limite, uu tend vers la solution de Neumann de moyenne nulle.