Aller au contenu principal

Cours · Bac+3 (ingénieur)

Modélisation

Fiches de cours et sous-notions liées à ce chapitre.

Connecte-toi pour t'entraîner

Pas assez de questions pour un entraînement ciblé sur cette notion.

Équations aux dérivées partielles

Non commencé

Définitions

Une EDP relie une fonction inconnue u(x,t)u(x,t) (ou u(x,y)u(x,y)) à ses dérivées partielles. Classification du second ordre : L[u]=auxx+buxy+cuyy+\mathcal{L}[u] = a u_{xx} + b u_{xy} + c u_{yy} + \cdots ; discriminant Δ=b24ac\Delta = b^2 - 4ac : elliptique (Δ<0\Delta < 0, Laplace), parabolique (Δ=0\Delta = 0, chaleur), hyperbolique (Δ>0\Delta > 0, ondes).

Conditions aux limites : Dirichlet (uΩ=gu|_{\partial\Omega} = g), Neumann (u/nΩ=h\partial u/\partial n|_{\partial\Omega} = h), Robin. Problème bien posé (Hadamard) : existence, unicité, stabilité.

EDP : classification elliptique / parabolique / hyperbolique via Δ=b24ac\Delta=b^2-4ac. Chaleur, Laplace, ondes.

CL Dirichlet / Neumann / Robin. Problème bien posé (Hadamard).

Pourquoi c’est central : cette notion structure les preuves et calculs du programme ; maîtriser définitions et hypothèses évite les applications hors cadre.

Formules

  • Chaleur : ut=αΔu(parabolique)\text{Chaleur : } u_t = \alpha \Delta u \quad \text{(parabolique)}
  • Laplace : Δu=0(elliptique)\text{Laplace : } \Delta u = 0 \quad \text{(elliptique)}
  • Ondes : utt=c2Δu(hyperbolique)\text{Ondes : } u_{tt} = c^2 \Delta u \quad \text{(hyperbolique)}
  • Δu=2ux2+2uy2\Delta u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2}

Exemples

Exemple 1

Type de ut=αuxxu_t=\alpha u_{xx}.

Méthode

(1) Lire l'énoncé et repérer les hypothèses / la forme utile. (2) Parabolique (chaleur). (3) Vérifier le résultat (ordre de grandeur, cas particulier, ou dérivation/substitution).

Résultat

Parabolique.

Exemple 2

Type de utt=c2Δuu_{tt}=c^2\Delta u.

Méthode

(1) Lire l'énoncé et repérer les hypothèses / la forme utile. (2) Hyperbolique (ondes). (3) Vérifier le résultat (ordre de grandeur, cas particulier, ou dérivation/substitution).

Résultat

Hyperbolique.

À retenir

Méthodes numériques

Non commencé

Définitions

Discrétisation des EDP : différences finies (DF) sur maillage régulier, éléments finis (EF) sur maillage triangulaire/tétraédrique, volumes finis (conservation locale). DF : u(xi)ui+12ui+ui1h2u''(x_i) \approx \frac{u_{i+1} - 2u_i + u_{i-1}}{h^2}.

Stabilité CFL pour explicite : ΔtCh2/α\Delta t \leq C h^2/\alpha (chaleur). EF : formulation faible uv=fv\int \nabla u \cdot \nabla v = \int fv sur espace VhV_h. Erreur O(hp)O(h^p) selon ordre éléments. Maillage adaptatif raffine zones de forte variation.

Discrétisation : DF, EF, VF. Schéma u(ui+12ui+ui1)/h2u''\approx(u_{i+1}-2u_i+u_{i-1})/h^2. CFL explicite chaleur : Δth2/(2α)\Delta t\leq h^2/(2\alpha).

EF : formulation faible. Erreur O(hp)O(h^p).

Pourquoi c’est central : cette notion structure les preuves et calculs du programme ; maîtriser définitions et hypothèses évite les applications hors cadre.

Formules

  • u(xi)ui+12ui+ui1h2+O(h2)u''(x_i) \approx \frac{u_{i+1} - 2u_i + u_{i-1}}{h^2} + O(h^2)
  • utnun+1unΔtu_t^n \approx \frac{u^{n+1} - u^n}{\Delta t}
  • CFL chaleur 1D : Δth22α\text{CFL chaleur 1D : } \Delta t \leq \frac{h^2}{2\alpha}
  • EF : trouver uhVh, a(uh,vh)=(vh) vhVh\text{EF : trouver } u_h \in V_h,\ a(u_h, v_h) = \ell(v_h)\ \forall v_h \in V_h

Exemples

Exemple 1

Condition CFL pour le schéma explicite de la chaleur 1D.

Méthode

(1) Lire l'énoncé et repérer les hypothèses / la forme utile. (2) Stabilité : Δth22α\Delta t \leq \frac{h^2}{2\alpha}. Si h=102h=10^{-2} et α=105\alpha=10^{-5}, alors Δt5\Delta t\leq 5. (3) Vérifier le résultat (ordre de grandeur, cas particulier, ou dérivation/substitution).

Résultat

Δth2/(2α)\Delta t \leq h^2/(2\alpha) (ex. Δt5\Delta t\leq 5 s pour ces valeurs).

Exemple 2

Ordre DF centrée ordre 2.

Méthode

(1) Lire l'énoncé et repérer les hypothèses / la forme utile. (2) Erreur de consistance O(h2)O(h^2) pour uu''. (3) Vérifier le résultat (ordre de grandeur, cas particulier, ou dérivation/substitution).

Résultat

Consistance O(h2)O(h^2).

À retenir