ACOUSTIQUE

Physique des Ondes et des Vibrations : De l'équation de d'Alembert à l'Imagerie Médicale.

MÉCANIQUE DES FLUIDES
ONDES
SIGNAL

Introduction : La Nature du Son

1. Qu'est-ce qu'une onde acoustique ?

Le son n'est pas un objet qui se déplace, mais une perturbation mécanique qui voyage à travers un milieu matériel (gaz, liquide ou solide). Sans matière (dans le vide), il n'y a pas de son (contrairement à la lumière).

"Les particules du milieu ne voyagent pas avec l'onde. Elles oscillent autour d'une position d'équilibre, transmettent leur énergie à la voisine, et reviennent à leur place."

On peut imaginer le milieu comme une chaîne infinie de masses reliées par des ressorts. Si vous poussez la première masse, la compression se propage de proche en proche. C'est un transport d'énergie sans transport de matière.

2. Propagation vs Diffusion : Le duel Physique

Ce cours couple Acoustique et Thermique car elles utilisent des outils mathématiques similaires, mais elles décrivent deux réalités opposées.

La Diffusion (Chaleur)

C'est un phénomène de "lissage". L'information se dilue. Si vous chauffez un point d'une barre de fer, la chaleur s'étale lentement. Si vous arrêtez de chauffer, la température s'uniformise. C'est un phénomène irréversible.

$$ \frac{\partial^2 T}{\partial x^2} = A \frac{\partial T}{\partial t} $$

Dérivée 1ère en temps (Lent)

La Propagation (Son)

C'est un phénomène de "transport". L'information voyage intacte (ou presque). Si vous criez, votre voix arrive à 340 m/s à l'autre bout. L'onde possède une "inertie" qui la fait avancer. C'est un phénomène réversible (en théorie).

$$ \frac{\partial^2 p}{\partial x^2} = \frac{1}{c^2} \frac{\partial^2 p}{\partial t^2} $$

Dérivée 2nde en temps (Rapide)

3. Les Modes de Propagation

La manière dont la matière se déforme définit le type d'onde. C'est une distinction cruciale pour comprendre pourquoi le son se comporte différemment dans l'eau et dans l'acier.

A. ONDE LONGITUDINALE (P-Wave)
Mécanisme : Compression / Dilatation.
Direction : La déformation est parallèle à la direction de propagation.
Milieux : Existe dans TOUS les milieux (Gaz, Liquides, Solides). C'est le mode principal du son aérien.
B. ONDE TRANSVERSE (S-Wave)
Mécanisme : Cisaillement (Shear).
Direction : La déformation est perpendiculaire à la propagation (comme une vague ou une corde de guitare).
Milieux : Existe UNIQUEMENT dans les SOLIDES (et fluides très visqueux).
Pourquoi pas de transverses dans l'eau ?

Pour qu'une onde transverse existe, il faut que le milieu résiste au cisaillement (si je tire une particule vers le haut, elle doit tirer sa voisine vers le haut). Les fluides parfaits (eau, air) n'ont pas de cohésion transversale : les couches glissent les unes sur les autres. Une onde de cisaillement s'y amortit instantanément.

4. La Double Périodicité

Une onde acoustique harmonique (sinusoïdale) possède une double identité. Elle se répète dans le temps, mais aussi dans l'espace.

TEMPS ($t$)

On regarde une particule fixe vibrer.

$T = \frac{1}{f}$

Période (s) & Fréquence (Hz)

$\times c =$
ESPACE ($x$)

On prend une photo de l'onde à un instant $t$.

$\lambda = c \cdot T$

Longueur d'onde (m)

Relation fondamentale : Une onde à haute fréquence ($f \nearrow$) a une courte longueur d'onde ($\lambda \searrow$).
Ex: À 1 kHz dans l'air, $\lambda \approx 34$ cm. En échographie (5 MHz) dans l'eau, $\lambda \approx 0.3$ mm.

Le Spectre Acoustique

Visualisation Fréquentielle

Passez votre souris sur les différentes zones pour visualiser le comportement de l'onde. Notez comment la période ($T$) diminue quand la fréquence ($f$) augmente.

20 Hz - 20 kHz AUDIBLE
1 Hz 100 Hz 10 kHz 1 MHz 1 GHz
INFRASONS
< 20 Hz
Séismes, Éléphants, Vibrations des ponts.
AUDIBLE
20 Hz - 20 kHz
Voix humaine, Musique. Zone sensible.
ULTRASONS
20 kHz - 1 GHz
Échographie, Sonar, CND.
HYPERSONS
> 1 GHz
Physique quantique (Phonons), Microscopie.

Milieux Limités & Guides

Quand l'onde rencontre des bords

Dans un milieu limité (surface libre, plaque), les réflexions multiples créent des modes guidés. Observez attentivement le mouvement des particules ci-dessous.

Onde de Rayleigh

Surface

Mouvement : Elliptique Rétrograde.
Regardez une particule en surface : elle tourne à l'envers de la vague.

Onde de Love

Couche/Substrat
VUE DE DESSUS

Mouvement : Cisaillement Horizontal.
Mouvement de "serpent" amorti en profondeur.

Onde de Lamb

Plaque

Mouvement : Flexion ($A_0$).
Toute l'épaisseur de la plaque ondule.

Phase vs Groupe (Dispersion)

Le Paquet d'Ondes et la Déformation

Dans la réalité, une impulsion sonore n'est jamais une fréquence unique, c'est un "paquet" d'ondes.

Vitesse de Phase ($V_{\phi}$)
$$ V_{\phi} = \frac{\omega}{k} $$

Vitesse des crêtes (vaguelettes)

Vitesse de Groupe ($V_g$)
$$ V_g = \frac{d\omega}{dk} $$

Vitesse de l'enveloppe (Énergie)

SIMULATION : NON-DISPERSIF (Air)

Milieu Non-Dispersif (ex: Air, Eau) : $V_{\phi} = V_g$.
L'enveloppe et les vaguelettes avancent à la même vitesse. Le signal ne se déforme pas.

Conséquence : Dans un milieu dispersif (solides, guides d'onde), $V_{\phi} \neq V_g$. Les différentes fréquences voyagent à des vitesses différentes.
Résultat : Le signal "s'étale" dans le temps et perd de l'amplitude.

🚴

L'Analogie du Peloton (Tour de France)

Imaginez un peloton de cyclistes vu d'hélicoptère.

Vitesse de Groupe ($V_g$)

LE PELOTON

C'est la vitesse à laquelle l'ensemble du groupe avance vers la ligne d'arrivée.
C'est la vitesse réelle de l'énergie (l'information).

Vitesse de Phase ($V_{\phi}$)

LE CYCLISTE

C'est la vitesse d'un cycliste individuel à l'intérieur du groupe. Il peut rouler plus vite que le groupe (remonter) ou moins vite (se laisser glisser).

En résumé : Dans un milieu dispersif, les cyclistes (les ondes) traversent le peloton. Ils apparaissent à la queue, remontent le groupe, et disparaissent en tête.
Le cycliste va vite ($V_{\phi}$), mais le message (le groupe) avance doucement ($V_g$).

Bref Historique

1877 - Théorie
Lord Rayleigh
Publie "The Theory of Sound". Il pose les bases mathématiques de l'acoustique moderne.
1880 - Matériaux
Pierre & Jacques Curie
Découverte de la Piézoélectricité (Quartz). C'est la capacité de transformer une pression mécanique en électricité (et inversement). La base de tous les transducteurs modernes.
1912/1917 - La Guerre & Le Titanic
Paul Langevin
Suite au naufrage du Titanic et à la guerre sous-marine, Langevin invente le premier SONAR actif (Sound Navigation and Ranging) utilisant des céramiques piézoélectriques pour détecter les sous-marins (l'ASDIC).
Années 1950 - Médical & Industrie
L'Essor des Ultrasons
Développement de l'échographie médicale (Wild & Reid) et du Contrôle Non Destructif (Firestone) pour vérifier les matériaux industriels sans les détruire.
DÉBUT DU COURS : PHYSIQUE ONDULATOIRE

I.1 L'Équation d'Onde

1. Le Cadre de Validité (Hypothèses)

Avant d'écrire la moindre équation, nous devons définir notre "monde" acoustique. Nous supposons un Fluide Idéal défini par 4 piliers :

1. Milieu Continu & Homogène

Le milieu est le même partout ($\rho_0$ constant). On ignore la structure atomique (l'échelle macroscopique $\gg$ distance inter-atomique).

2. Isotrope

Les propriétés physiques sont les mêmes dans toutes les directions. Le son se propage de la même façon vers le haut, le bas, la gauche...

3. Fluide Parfait

Pas de viscosité. On néglige les frottements internes et les pertes thermiques (transformations adiabatiques). L'énergie se conserve localement.

4. Petits Mouvements

"L'Acoustique Linéaire". La pression acoustique $p$ est négligeable devant la pression atmosphérique $P_0$.
Conséquence math : On néglige les termes d'ordre 2 (ex: $v^2 \approx 0$).

2. Définition des Grandeurs (Notation)

Avant de démarrer la démonstration, définissons précisément les termes utilisés pour une tranche de fluide au repos située en $x$.

$u(x,t)$

DÉPLACEMENT (m)

Position par rapport à l'équilibre.

$\rho_0$

MASSE VOLUMIQUE (kg/m³)

Constante (milieu homogène).

$\theta$

DILATATION CUBIQUE

Variation relative de volume ($dV/V$).

$\chi_s$

COMPRESSIBILITÉ ($Pa^{-1}$)

Coefficient adiabatique du milieu.

3. Démonstration : Tranche de Fluide (3D)

On considère une tranche de fluide élémentaire d'épaisseur $dx$ et de section $S$. Masse : $dm = \rho_0 S dx$.

A. Équation de la Dynamique NEWTON

La force nette est la différence de pression entre $x$ et $x+dx$. On effectue un développement limité au 1er ordre :

$$ \sum F = [P(x) - P(x+dx)] \cdot S $$ $$ \approx [P(x) - (P(x) + \frac{\partial p}{\partial x} dx)] \cdot S = - \frac{\partial p}{\partial x} dx S $$
RÉSULTAT (1) : $\rho_0 \frac{\partial^2 u}{\partial t^2} = - \frac{\partial p}{\partial x}$
$P(x)S$ $P(x+dx)S$ $dm$

B. Conservation de la Masse DILATATION

Le volume varie car le déplacement $u$ n'est pas uniforme.
• $V_0 = S \cdot dx$
• $V = S \cdot (dx + du) = S(dx + \frac{\partial u}{\partial x}dx)$

$$ \theta = \frac{V - V_0}{V_0} = \frac{S \frac{\partial u}{\partial x} dx}{S dx} = \frac{\partial u}{\partial x} $$
RÉSULTAT (2) : $\theta = \partial_x u$
REPOS ($dx$) DILATÉ $u(x)$ $u(x+dx)$

C. Loi de Comportement

On utilise le coefficient de compressibilité adiabatique $\chi_s$ (définition thermodynamique) :

$$ \chi_s = - \frac{1}{V} \left( \frac{\partial V}{\partial P} \right)_s $$

On intègre pour les petites variations ($dV/V = \theta$) :

$$ p = - \frac{1}{\chi_s} \theta = - \frac{1}{\chi_s} \frac{\partial u}{\partial x} \quad \text{(3)} $$

D. ÉQUATION FINALE

On dérive (3) par rapport à $x$ et on injecte dans (1) :
$\frac{\partial p}{\partial x} = - \frac{1}{\chi_s} \frac{\partial^2 u}{\partial x^2}$

$$ \frac{\partial^2 u}{\partial x^2} - \rho_0 \chi_s \frac{\partial^2 u}{\partial t^2} = 0 $$

4. Généralisation (3D & Variables)

Ce que nous avons démontré pour le déplacement $u$ en 1D est en réalité une loi universelle de l'acoustique linéaire.

Universalité des Variables

La Pression $p$, la Vitesse $v$, et la Masse Volumique $\rho$ obéissent toutes à la même équation.

$$ \frac{\partial^2 (\dots)}{\partial x^2} - \frac{1}{c^2} \frac{\partial^2 (\dots)}{\partial t^2} = 0 $$

Remplacez (...) par $p, v, \rho$ ou $u$.

Passage en 3D (Laplacien)

Dans l'espace, la dérivée seconde $\partial^2/\partial x^2$ est remplacée par l'opérateur Laplacien $\Delta$.

$$ \Delta p - \frac{1}{c^2} \frac{\partial^2 p}{\partial t^2} = 0 $$

$\Delta = \partial_x^2 + \partial_y^2 + \partial_z^2$

5. La Solution Mathématique

D'Alembert a montré que toute solution de cette équation s'écrit comme la somme de deux ondes voyageant en sens opposé.

$$ p(x,t) = \underbrace{f(t - \frac{x}{c})}_{\text{Onde Progressive } (\rightarrow)} + \underbrace{g(t + \frac{x}{c})}_{\text{Onde Régressive } (\leftarrow)} $$
Onde Progressive ($x \nearrow$)

Terme en $(t - x/c)$. Le profil avance vers les x positifs.

Onde Régressive ($x \searrow$)

Terme en $(t + x/c)$. Le profil recule (ex: réflexion sur un mur).

6. L'Outil Mathématique : Notation Complexe

En physique, manipuler des sinus et cosinus est pénible (formules triguinometriques). On préfère utiliser l'exponentielle complexe pour transformer les équations différentielles en équations algébriques simples.

L'Onde Plane Harmonique

On note la solution sous forme :

$$ \underline{p}(x,t) = A \cdot e^{i(\omega t - kx)} $$
  • $\omega$ : Pulsation temporelle ($2\pi f$)
  • $k$ : Vecteur d'onde ($2\pi / \lambda$)
⚠️ Règle d'Or

1. On fait tous les calculs avec les complexes ($e^{i\dots}$).
2. À la toute fin, pour avoir le sens physique, on ne garde que la Partie Réelle.

$$ p_{reelle} = \text{Re}(\underline{p}) = A \cos(\omega t - kx) $$
Relation de Dispersion

En réinjectant cette solution dans l'équation d'onde, on trouve le lien fondamental entre temps et espace :

$$ k = \frac{\omega}{c} $$

6. Physique de la Célérité (Vitesse)

La vitesse du son n'est pas une constante universelle. C'est une compétition entre deux propriétés du milieu :

$$ c = \sqrt{\frac{\text{Propriété Élastique (Ressort)}}{\text{Propriété Inertielle (Masse)}}} = \frac{1}{\sqrt{\rho_0 \chi_s}} $$

PLUS C'EST RAIDE, PLUS C'EST VITE. PLUS C'EST LOURD, PLUS C'EST LENT.

A. Le Cas des Gaz (Newton vs Laplace)

Le débat historique portait sur la nature thermodynamique de la transformation.

1687
Isaac Newton

Hypothèse : ISOTHERME ($T = cte$)
Il pensait que la chaleur avait le temps de s'équilibrer entre les zones comprimées et dilatées.

Compressibilité isotherme :

$\chi_T = 1 / P_0$
$c = \sqrt{P_0 / \rho_0} \approx 280 \text{ m/s}$

❌ FAUX (-20%)

1816
Pierre-Simon de Laplace

Hypothèse : ADIABATIQUE ($Q = 0$)
Le son va trop vite. La chaleur reste "piégée" dans les compressions, ce qui augmente la raideur du gaz.

Compressibilité adiabatique :

$\chi_s = 1 / (\gamma P_0)$
$$ c = \sqrt{\frac{\gamma P_0}{\rho_0}} \approx 340 \text{ m/s} $$

✅ CORRECT

De la Pression à la Température

En utilisant la loi des Gaz Parfaits ($\frac{P_0}{\rho_0} = \frac{RT}{M}$), on obtient la forme usuelle qui montre que $c$ ne dépend que de la température :

$$ c_{gaz} = \sqrt{\frac{\gamma R T}{M}} $$

B. Comparaison : Liquides et Solides

Eau (Liquide) ~1500 m/s

L'eau est 800x plus lourde que l'air ($\rho \nearrow$), ce qui devrait ralentir le son.
MAIS, elle est incompressible ($\chi \searrow \searrow$).

$$ c = \frac{1}{\sqrt{\rho_0 \chi}} $$

La raideur l'emporte sur l'inertie.

Acier (Solide) ~5000 m/s

On remplace la compressibilité par l'inverse du Module d'Young ($E$). Les solides transmettent l'énergie très efficacement.

$$ c_{barre} = \sqrt{\frac{E}{\rho_0}} $$

Onde de compression (Longitudinale).

🎮 Simulateur : Vitesse dans l'Air

Voyez l'influence de la température sur la célérité.

  • $\gamma = 1.4$ (Air diatomique)
  • $M = 29$ g/mol
343.2 m/s
❄️ 🔥
T = 20°C

IV. Champs Réels & Diffraction

1. L'Hypothèse de l'Onde Plane vs Réalité

Dans les exercices, on aime l'Onde Plane (qui avance comme un mur infini). Dans la réalité, les sources (haut-parleur, piston) ont une taille finie. Elles créent des ondes complexes.

Principe de Huygens-Fresnel

"Toute surface émettrice peut être considérée comme une infinité de petites sources ponctuelles émettant des ondes sphériques."

$$ P_{tot} = \sum \text{Ondes Sphériques} $$
Interférences

2. Structure du Faisceau (Zone de Fresnel & Fraunhofer)

Si on regarde le champ de pression devant un émetteur (ex: Sonde échographique), on distingue deux zones (voir Slide 5).

Source
Zone de Fresnel (Proche)
Limite $a^2/\lambda$
Zone de Fraunhofer (Lointain)
Champ Proche (Fresnel)

Zone de turbulences et d'interférences. L'amplitude varie énormément. On ne peut pas y faire de mesures fiables.
Limite : $Y \approx a^2 / \lambda$

Champ Lointain (Fraunhofer)

L'onde s'organise et devient sphérique. L'énergie décroît régulièrement en $1/r$. C'est ici qu'on utilise le sonar/écho.

I.2 Le Potentiel des Vitesses ($\Phi$)

1. Définition et Démonstration (PFD)

Lorsque seules des ondes longitudinales sont considérées, on peut simplifier l'étude en remplaçant le vecteur vitesse $\vec{v}$ par un scalaire $\Phi$.

Démonstration : Du PFD au Potentiel

Reprenons le Principe Fondamental de la Dynamique (PFD) en 3D pour un fluide parfait (sans forces externes) :

$$ \rho_0 \frac{\partial \vec{v}}{\partial t} = - \vec{\nabla} p $$

On introduit le potentiel tel que $\vec{v} = - \vec{\nabla} \Phi$. L'équation devient :

$$ \rho_0 \frac{\partial}{\partial t} (- \vec{\nabla} \Phi) = - \vec{\nabla} p \implies \vec{\nabla} (\rho_0 \frac{\partial \Phi}{\partial t}) = \vec{\nabla} p $$

En intégrant spatialement (à une constante près négligeable), on obtient la relation fondamentale pression-potentiel :

$$ p(x,t) = \rho_0 \frac{\partial \Phi}{\partial t} $$

Note : En cas d'onde transverse (solides), il faudrait introduire un potentiel vectoriel via un rotationnel ($\vec{v} = \vec{\nabla} \wedge \vec{\Psi}$), mais c'est hors programme ici.

2. Stratégie de Résolution

Au lieu de chercher $p$ ou $\vec{v}$ séparément, la méthode universelle en acoustique est la suivante :

1
Trouver $\Phi$

On résout d'abord l'équation scalaire :

$$ \Delta \Phi - \frac{1}{c^2} \frac{\partial^2 \Phi}{\partial t^2} = 0 $$
2
Déduire le reste

Tout découle de $\Phi$ par simple dérivation :

  • $\vec{v} = - \vec{\nabla} \Phi$
  • $p = \rho_0 \partial_t \Phi$

3. Pourquoi l'Onde Sphérique ? (Réalité vs Modèle)

L'hypothèse de l'onde plane (milieu infini) est souvent fausse en pratique. Pourquoi ?

Exemple Acoustique Audible
  • Célérité $c \approx 340$ m/s
  • Fréquence $f = 400$ Hz (La)
  • $\lambda = c/f \approx 0.85$ m

La longueur d'onde est de l'ordre du mètre ! Les sources réelles (HP, cordes vocales) sont rarement beaucoup plus grandes que $\lambda$.

Source Réelle DIFFRACTION

"Les ondes émises ne sont pas planes."

Conséquence (Principe de Huygens / Intégrale de Rayleigh) :
Toute source quelconque est vue comme une superposition de sources ponctuelles émettant des ondes sphériques. D'où l'importance de la solution harmonique sphérique :

$$ \Phi(r, t) = \frac{A}{r} e^{i(\omega t - kr)} $$

Notez le terme en $1/r$ : l'amplitude décroît avec la distance (dilution de l'énergie).

III. Aspects Énergétiques

1. Énergies mises en jeu (Cinétique & Potentielle)

Comme un ressort qui oscille, l'onde acoustique transporte deux formes d'énergie qui s'échangent en permanence. Pour un volume unitaire, on parle de densité d'énergie :

Énergie Cinétique ($E_c$)

Liée à la Vitesse des particules.

$$ e_c = \frac{1}{2} \rho_0 v^2 $$
Énergie Potentielle ($E_p$)

Liée à la Compression (Pression).

$$ e_p = \frac{p^2}{2 \rho_0 c^2} $$

Pour une onde plane progressive, on montre que $e_c = e_p$. L'énergie est équi-répartie.

2. Intensité Acoustique ($I$) et Valeur Efficace

L'intensité instantanée est le flux de puissance : $\vec{I}(t) = p(t) \cdot \vec{v}(t)$.
Problème : Pour un signal sinusoïdal, la moyenne est nulle ! (Ça oscille autour de 0).

Solution : On utilise la valeur efficace (RMS - Root Mean Square).

Formule Pratique (Onde Plane)
$$ I = \frac{p_{eff}^2}{\rho_0 c} \quad \text{avec} \quad p_{eff} = \frac{p_{max}}{\sqrt{2}} $$

C'est cette valeur que mesurent les sonomètres.

Représentation Graphique

Signal $p(t)$ (Bleu) vs $p_{eff}$ (Rouge)

3. L'Échelle des Décibels (Logarithmique)

L'oreille gère une dynamique énorme ($10^6$ en pression, $10^{12}$ en intensité). On utilise donc une échelle log.

Niveau de Pression ($L_p$)

$$ L_p = 20 \log_{10} \left( \frac{p_{eff}}{20 \mu Pa} \right) $$

Niveau d'Intensité ($L_I$)

$$ L_I = 10 \log_{10} \left( \frac{I}{10^{-12} W/m^2} \right) $$
Remarque : La limite physique du son

La pression acoustique est une fluctuation autour de la pression atmosphérique $P_{atm}$ ($10^5$ Pa).
Comme on ne peut pas descendre en dessous du vide (pression nulle), l'amplitude maximale théorique d'une onde sinusoïdale dans l'air est $P_{max} = P_{atm}$.

$$ L_{max} = 20 \log \left( \frac{10^5 / \sqrt{2}}{20 \cdot 10^{-6}} \right) \approx 191 \text{ dB} $$

Au-delà, l'air ne se comporte plus linéairement (ondes de choc).

⚠️ L'Arithmétique des Décibels

Les décibels ne s'additionnent pas directement ! $50 \text{ dB} + 50 \text{ dB} \neq 100 \text{ dB}$.
Il faut repasser par les puissances (échelle linéaire), additionner, puis revenir aux log.

Formule Générale (Sources Incohérentes)
$$ L_{total} = 10 \log_{10} \left( \sum_{i=1}^{N} 10^{\frac{L_i}{10}} \right) $$
Cas : Doublement de source

2 enceintes identiques

+ 3 dB

$10 \log(2) \approx 3$

Cas : Décuplement

10 violons au lieu d'1

+ 10 dB

$10 \log(10) = 10$

IV. Réflexion & Transmission

1. Introduction & Hypothèses

En pratique, les milieux ne sont jamais infinis. L'onde rencontre des interfaces. Pour simplifier l'étude mathématique (comme vu en Slide 55), on se place dans un cadre précis :

1. Ondes Planes

On suppose des sources étendues et une observation loin de la source (diffraction négligée).

2. Milieux Non Dissipatifs

L'impédance $Z = \rho c$ est réelle (pas de partie imaginaire, pas d'absorption). Le vecteur d'onde $k$ est réel.

2. Conditions de Continuité à l'Interface

Au niveau de la séparation (en $x=0$), la matière ne se déchire pas et il n'y a pas de vide qui se crée. Cela impose la continuité de deux grandeurs :

  • Pression ($p$) : Les forces de part et d'autre sont égales.
  • Vitesse Normale ($v_n$) : Ou déplacement normal. Le fluide 1 reste collé au fluide 2.
⚠️ MISE EN GARDE (Slide 56)

Ne surtout pas confondre :
$c$ (Célérité) : Vitesse de l'onde (~340 m/s). Dépend du milieu.
$v$ (Vitesse Particulaire) : Vitesse de vibration des molécules (~mm/s). Dépend de l'amplitude.
La continuité s'applique sur $v$, pas sur $c$ !

3. Incidence Oblique : Démonstrations

Considérons une onde incidente arrivant avec un angle $\theta_i$. Elle génère une réfléchie ($\theta_r$) et une transmise ($\theta_t$).

$\vec{k}_i$ $\vec{k}_r$ $\vec{k}_t$ Interface (x=0)
DÉMONSTRATION : Pourquoi Snell-Descartes ? (Slide 59)

Écrivons les champs de pression à l'interface ($x=0$, et quel que soit $y$) :

$$ \underbrace{A e^{i(\omega t - k_i y \sin \theta_i)}}_{\text{Incidente}} + \underbrace{B e^{i(\omega t - k_r y \sin \theta_r)}}_{\text{Réfléchie}} = \underbrace{C e^{i(\omega t - k_t y \sin \theta_t)}}_{\text{Transmise}} $$

Pour que cette égalité soit vraie pour tout $y$ et tout $t$, il faut impérativement que les phases soient identiques :

$$ k_i \sin \theta_i = k_r \sin \theta_r = k_t \sin \theta_t $$

Sachant que $k = \omega / c$ et que la fréquence $\omega$ est conservée :

LOI DE SNELL-DESCARTES
$$ \frac{\sin \theta_i}{c_1} = \frac{\sin \theta_r}{c_1} = \frac{\sin \theta_t}{c_2} $$

Implique $\theta_i = \theta_r$ (Réflexion) et $\sin \theta_t = \frac{c_2}{c_1} \sin \theta_i$ (Réfraction).

Relation des Amplitudes (Slide 59) :
En simplifiant les exponentielles (qui sont égales), l'équation de continuité des pressions donne simplement :

$$ 1 + r = t $$

4. Interface Liquide / Solide (Conversion de Mode)

C'est le cas général le plus riche. Un solide supporte deux types d'ondes : Longitudinales ($L$) et Transverses ($T$).
Une onde incidente pure ($L$) dans l'eau va générer deux ondes réfractées dans le solide.

Onde L (Longitudinale)

Célérité rapide ($C_{2L}$)

$\sin \theta_{tL} = \frac{C_{2L}}{C_{1L}} \sin \theta_i$
Onde T (Transverse)

Célérité plus lente ($C_{2T} < C_{2L}$)

$\sin \theta_{tT} = \frac{C_{2T}}{C_{1L}} \sin \theta_i$

Les Angles Critiques (Réflexion Totale)

Comme la célérité dans le solide est généralement plus grande que dans le liquide ($C_2 > C_1$), l'angle transmis $\theta_t$ atteint $90^\circ$ ($\pi/2$) avant $\theta_i$. L'onde devient "rasante" puis disparaît (évanescente).

Hiérarchie des Angles Critiques

Comme $C_{2L} > C_{2T}$, l'onde Longitudinale disparaît en premier !

1. Premier Angle Critique ($\theta_{icL}$)

$$ \theta_{icL} = \arcsin(C_1 / C_{2L}) $$

$\theta > \theta_{icL}$ : L'onde L transmise disparaît.

2. Second Angle Critique ($\theta_{icT}$)

$$ \theta_{icT} = \arcsin(C_1 / C_{2T}) $$

$\theta > \theta_{icT}$ : L'onde T transmise disparaît aussi. Réflexion totale.

Phénomène Onde de Surface

Juste après ces angles critiques, des ondes de surface complexes (type Rayleigh) sont générées à l'interface (combinaison des ondes évanescentes).

VII. Systèmes Viscoélastiques

Du modèle idéal à la réalité : Frottements, Viscosité et Absorption.

1. Pourquoi ce chapitre ? (Objectifs)

Jusqu'ici, nous avons supposé que l'énergie se conservait indéfiniment. Dans la réalité (et surtout dans le corps humain pour l'imagerie médicale), ce n'est pas le cas. Le milieu résiste et absorbe de l'énergie.

Les Sources de Pertes (Atténuation Globale)

L'énergie de l'onde diminue à cause de deux phénomènes majeurs :

  • L'Absorption (Viscosité/Frottements) : Conversion de l'énergie acoustique en chaleur (friction solide ou viscosité fluide).
  • La Diffusion (Scattering) : Hétérogénéités du milieu (bulles, fibres, grains, nœuds dans le bois) qui éparpillent l'onde.
$$ \alpha_{globale} = \alpha_{friction} + \alpha_{inclusion} + \alpha_{viscosité} + \dots $$

2. La Notion de Viscosité de Cisaillement

La viscosité traduit les frottements entre les couches de fluide qui glissent les unes sur les autres (Écoulement de Couette).

Contrainte de Cisaillement ($\tau$)

Force $F$ appliquée sur une surface $S$.

$$ \tau = \frac{F}{S} $$

Unité : Pascal (Pa)

Taux de Cisaillement ($\dot{\gamma}$)

Gradient de vitesse ($V_{max}$ sur hauteur $H$).

$$ \dot{\gamma} = \frac{\partial V}{\partial z} \approx \frac{V_{max}}{H} $$

Unité : $s^{-1}$

Loi de Comportement (Fluide Newtonien)

La viscosité dynamique $\eta$ est le coefficient de proportionnalité qui relie la cause ($\dot{\gamma}$) à l'effet ($\tau$) :

$$ \tau = \eta \cdot \dot{\gamma} $$

Unité de $\eta$ : Pa.s (Pascal seconde)

  • Si $\eta$ est constant (indépendant de $\dot{\gamma}$) : Fluide Newtonien (Eau, Glycérine).
  • Si $\eta$ dépend de la vitesse : Fluide Non-Newtonien (Sang).

3. Conséquences : L'Atténuation Exponentielle

L'ajout du terme de viscosité dans l'équation d'Euler (Navier-Stokes) entraîne une perte d'amplitude progressive de l'onde.

Onde Plane
$$ p(x,t) = A_0 e^{-\alpha x} e^{i(\omega t - kx)} $$
Onde Sphérique
$$ p(r,t) = \frac{A_0}{r} e^{-\alpha r} e^{i(\omega t - kr)} $$

$\alpha$ : Coefficient d'atténuation acoustique
Unité : $m^{-1}$ ou $Neper.m^{-1}$

Vecteur d'Onde Complexe ($k^*$)

Pour une onde transverse se propageant à la célérité $c_T$, on définit un vecteur d'onde complexe qui inclut l'atténuation :

$$ k^* = \frac{\omega}{c_T} - i\alpha $$
  • $\frac{\omega}{c_T}$ : Partie réelle (Propagation classique)
  • $\alpha = -Im(k^*)$ : Partie imaginaire (Atténuation transverse)
LA CÉLÉRITÉ COMPLEXE ($c_T^*$)

On introduit alors aussi une célérité complexe définie par la relation :

$$ k^* = \frac{\omega}{c_T^*} $$

Attention à ne pas confondre $c_T$ (célérité réelle de phase) et $c_T^*$ (célérité complexe incluant les pertes).

LOI EN FRÉQUENCE (Important)

Pour un fluide Newtonien (viscosité constante), l'atténuation $\alpha$ dépend du carré de la fréquence :

$$ \alpha \propto f^2 $$

4. Viscoélasticité (Modules Complexes)

Dans les solides mous ou visqueux, on définit un module de cisaillement complexe $G^*$ qui combine élasticité (réel) et pertes (imaginaire).

$$ G^* = G' + iG'' $$
G' (Stockage)

Partie Élastique (Ressort)

$$ G' = \rho c_T^2 \frac{1-(\alpha c_T/\omega)^2}{(1+(\alpha c_T/\omega)^2)^2} $$
G'' (Perte)

Partie Visqueuse (Amortisseur)

$$ G'' \approx \omega \eta $$
Relation Viscosité - Atténuation

Si l'amortissement est faible, on peut calculer la viscosité $\eta$ à partir de l'atténuation mesurée $\alpha$ et la vitesse $c$ :

$$ \eta \approx \frac{2 \rho \alpha c^3}{\omega^2} $$

Valable pour ondes transverses ($c_T$, $\alpha_T$) et longitudinales ($c_L$, $\alpha_L$).

⚠️ CHAÎNE DE CONSÉQUENCES (Slide 73)

Tout part de la friction interne (Viscosité). Cette perte d'énergie impose une diminution de l'amplitude (Atténuation).

Pour modéliser cette perte, les mathématiques nous obligent à utiliser des grandeurs complexes. Cela se propage jusqu'à l'impédance acoustique.
Conséquence physique concrète : Contrairement à un milieu parfait, une onde se réfléchissant sur un milieu visqueux subira un décalage temporel (Déphasage) dû à cette partie imaginaire.

Viscosité ($\eta$) [Pa.s]
Atténuation ($\alpha$) [$m^{-1}$]
Vecteur d'onde complexe ($k^*$)
Célérité complexe ($c^*$) & Impédance ($Z^*$)
DÉPHASAGE À LA RÉFLEXION

V. Applications & Ingénierie

De la Piézoélectricité à l'Imagerie Médicale et la Thérapie.

1. L'Effet Piézoélectrique ($\pi\iota\epsilon\zeta\omega$)

Slide 77-78

Pour faire de l'échographie, il faut un convertisseur Électro-Mécanique. On utilise des matériaux piézoélectriques (découverts par Pierre & Jacques Curie, 1880).

Émission (Effet Inverse)

Tension Électrique $\to$ Déformation Mécanique

On envoie une impulsion électrique, le cristal vibre et crée l'onde.

Réception (Effet Direct)

Pression Mécanique $\to$ Tension Électrique

L'écho revient, comprime le cristal qui génère un courant.

Autres méthodes (Slide 76) : Laser (optique), EMAT (magnétique), cMUT (condensateur).

2. Le Principe Pulse-Echo

C'est la base du SONAR (Langevin, 1915) et de l'échographie (Wild & Reid, 1950). On mesure deux grandeurs fondamentales sur le signal retour (Slide 79) :

Temps de Vol ($\Delta t$)
$\to$ Donne la Distance et la Rigidité (via la vitesse $c$).
Amplitude ($A$)
$\to$ Donne l'Atténuation (Viscosité, Pathologie).
$$ d = \frac{c \cdot \Delta t}{2} $$ (On divise par 2 car l'onde fait l'aller-retour).

3. La Sonde : Focalisation Électronique

Slide 85

Une sonde moderne n'est pas un mono-élément. C'est une barrette de centaines de petits cristaux. En jouant sur les retards d'excitation entre les éléments, on peut manipuler le faisceau sans bouger la sonde.

Balayage (Steering)

Retard linéaire.

Le faisceau part de biais. Permet de scanner une zone large (secteur).

Focalisation (Focusing)

Retard parabolique.

L'onde converge en un point précis pour une résolution maximale.

4. L'Effet Doppler (Flux)

Si la cible bouge (globules rouges), la fréquence de l'écho change.

$$ f_r - f_i = \frac{2v}{c} f_i \cos \theta $$
  • Rouge : Le sang s'approche de la sonde ($f_r > f_i$).
  • Bleu : Le sang s'éloigne ($f_r < f_i$).

5. Thérapie (HIFU)

High Intensity Focused Ultrasound.
On ne cherche plus à voir, mais à brûler. En focalisant une forte puissance, on détruit des tumeurs ou on coagule des tissus sans chirurgie.

Problème majeur : La localisation précise (guidage par IRM ou écho).

6. Elastographie ("Palpation à Distance")

Slide 92

Technique de rupture. En échographie standard, la vitesse longitudinale ($V_L$) varie peu entre un tissu sain et une tumeur.
L'idée : On génère une onde de cisaillement (Onde T). Sa vitesse ($V_T$) dépend directement de la rigidité (dureté) du tissu.

👉
Application : Cancer du sein/prostate

Une tumeur est plus dure qu'un tissu sain $\to$ $V_T$ augmente localement.

7. Applications Matériaux & Agroalimentaire

Slide 93-97

Les ultrasons sondent la matière opaque là où la lumière s'arrête.

🪵
Le Bois & Patrimoine

Le bois est anisotrope (la vitesse dépend du sens des fibres).
Exemple : Analyse de la dégradation du navire Vasa (1626) sans prélèvement destructif.

🍯
Le Miel

Suivi de la maturation et de la cristallisation. La vitesse et l'atténuation changent drastiquement selon la teneur en sucre et la viscosité.

🌡️
Thermométrie

Comme vu au Chap. I, $c \propto \sqrt{T}$.
Mesurer la vitesse permet de connaître la température à cœur d'un organe (ex: ablation cardiaque) de façon non-invasive.

VI. Exercices & Applications

Démonstrations mathématiques et calculs basés sur les supports de cours.

Exercice 1 : Propagation Harmonique

Slide 27

Montrer que les fonctions $u(x,t) = A e^{i(\omega t - kx)}$ et $u(x,t) = B e^{i(\omega t + kx)}$ sont solutions de l'équation de propagation : $$ \frac{\partial^2 u}{\partial x^2} - \frac{1}{c^2} \frac{\partial^2 u}{\partial t^2} = 0 $$ Condition : On pose $k = \omega / c$.

VOIR LA DÉMONSTRATION DÉTAILLÉE ▼

Pour la solution progressive $u = A e^{i(\omega t - kx)}$ :

Dérivées spatiales $$ \frac{\partial u}{\partial x} = -ik \cdot u $$ $$ \frac{\partial^2 u}{\partial x^2} = (-ik)^2 \cdot u = -k^2 u $$
Dérivées temporelles $$ \frac{\partial u}{\partial t} = i\omega \cdot u $$ $$ \frac{\partial^2 u}{\partial t^2} = (i\omega)^2 \cdot u = -\omega^2 u $$

En remplaçant dans l'équation d'onde :

$$ -k^2 u - \frac{1}{c^2} (-\omega^2 u) = 0 \implies -k^2 u + \frac{\omega^2}{c^2} u = 0 $$

L'équation est vérifiée si et seulement si $k^2 = \frac{\omega^2}{c^2}$, soit $k = \frac{\omega}{c}$.
Note : Le raisonnement est identique pour le terme en $+kx$ (onde régressive).

Exercice 2 : Équation de Pression

Slide 27

Sachant que l'on a obtenu une équation de propagation pour le déplacement $u(x,t)$, montrer que la pression $p(x,t)$ vérifie la même équation en utilisant la relation : $$ p = -\frac{1}{\chi} \frac{\partial u}{\partial x} $$

VOIR LES MÉTHODES DE RÉSOLUTION ▼

Méthode 1 (Algébrique) :

On remplace $p$ par sa définition dans l'équation d'onde. Comme $\chi$ est une constante, les opérateurs de dérivation s'appliquent directement sur $\partial u / \partial x$.

Méthode 2 (Notation complexe) :

Si $u = A e^{i(\omega t - kx)}$, alors :

$$ p = -\frac{1}{\chi} (-ik u) = \frac{ik}{\chi} u $$

En utilisant la forme $i = e^{i\pi/2}$ :

$$ p = \frac{k}{\chi} A e^{i(\omega t - kx + \pi/2)} $$

Conclusion : $p$ et $u$ sont proportionnels et liés par une relation linéaire. Puisque $u$ suit l'équation de d'Alembert, $p$ la suit également. On remarque un déphasage de $\pi/2$ (quadrature).

Exercice 3 : Célérité & Température

Slide 29

On suppose le passage de l'onde adiabatique réversible. La vitesse du son dans un gaz parfait est $c = (\gamma RT / M)^{0.5}$.
1. Montrer que pour l'air ($\gamma=1.4$, $M=29$ g/mol), $c \approx 20 \sqrt{T}$.
2. Calculer $c$ pour $T = 300$ K.

VOIR LA RÉSOLUTION NUMÉRIQUE ▼

1. Simplification de la formule :

$$ c = \sqrt{\frac{1.4 \times 8.314}{0.029}} \cdot \sqrt{T} \approx \sqrt{401.3} \cdot \sqrt{T} \approx 20 \sqrt{T} $$

2. Application numérique à 300K :

$$ c = 20 \cdot \sqrt{300} \approx 20 \cdot 17.32 \approx \mathbf{346.4 \, m/s} $$

Note : À 20°C (293K), on retrouve la valeur usuelle de 340 m/s.

Questions de Réflexion

Conceptuel

"C'est quoi un gaz parfait ? Que veut dire adiabatique ? Réversible ? Est-ce réaliste ?"

RÉPONSES SYNTHÉTIQUES ▼
  • Gaz Parfait : Modèle sans interaction entre molécules (hormis chocs) et sans volume propre. Excellent pour l'air à pression atmosphérique.
  • Adiabatique : Évolution sans échange thermique. Justifié ici car la propagation du son est trop rapide pour permettre des transferts de chaleur significatifs.
  • Réversible : Transformation sans pertes d'énergie (frottements/viscosité négligés).
  • Réalisme : Oui, pour les petits déplacements (acoustique linéaire). Les limites apparaissent avec les "gaz comprimés" ou les ondes de choc.

Exercice 4 : L'Onde Sphérique

Slide 36-37

En réalité, les sources émettent des ondes sphériques. Le Laplacien en coordonnées sphériques (cas radial pur) s'écrit : $\Delta = \frac{1}{r} \frac{\partial^2 (r \cdot)}{\partial r^2}$.

1. Montrer que $p(r,t) = \frac{f(t - r/c)}{r}$ est solution de l'équation de d'Alembert.
2. Quelle est la conséquence physique du terme en $1/r$ sur l'amplitude et l'énergie ?

DÉMONSTRATION & INTERPRÉTATION ▼
1. Vérification de la solution

On injecte la forme proposée dans l'équation : $\Delta p - \frac{1}{c^2} \frac{\partial^2 p}{\partial t^2} = 0$.

Terme Spatial (Laplacien)

$$ \Delta p = \frac{1}{r} \frac{\partial^2}{\partial r^2} \left( r \cdot \frac{f(t-r/c)}{r} \right) $$ $$ = \frac{1}{r} \frac{\partial^2 f}{\partial r^2} $$

On pose $u = t - r/c$, donc $\frac{\partial f}{\partial r} = -\frac{1}{c} f'$.

$$ \Delta p = \frac{1}{r} \left( \frac{1}{c^2} f'' \right) $$

Terme Temporel

$$ \frac{1}{c^2} \frac{\partial^2 p}{\partial t^2} = \frac{1}{c^2} \frac{\partial^2}{\partial t^2} \left( \frac{f}{r} \right) $$

Ici $r$ est constant par rapport au temps.

$$ = \frac{1}{c^2} \frac{1}{r} f'' $$

Les deux termes sont identiques : l'équation est vérifiée.

2. Conséquence Physique (Atténuation Géométrique)

Contrairement à l'onde plane (amplitude constante), l'onde sphérique voit son amplitude décroître en $1/r$.

L'énergie est proportionnelle au carré de la pression ($E \propto p^2$).
Donc l'énergie décroît en $1/r^2$.

C'est logique : l'énergie se répartit sur la surface d'une sphère ($4\pi r^2$) qui grandit. Même sans frottement, le son s'affaiblit avec la distance.

Exercice 5 : Application Onde Plane (Relation p/v)

Slide 42

On considère une onde plane harmonique.
1. Montrer que pour une onde se propageant vers les $x$ croissants : $p = \rho_0 c \cdot v$.
2. Montrer que pour une onde se propageant vers les $x$ décroissants : $p = -\rho_0 c \cdot v$.
3. En déduire que l'intensité acoustique s'écrit $I = p^2 / \rho_0 c$.

DÉMONSTRATION DÉTAILLÉE ▼

1. Cas de l'onde progressive ($x \nearrow$)

On part du potentiel harmonique progressif : $\Phi = A e^{i(\omega t - kx)}$.

Calcul de la Vitesse

$$ v = - \frac{\partial \Phi}{\partial x} = -(-ik)\Phi = ik\Phi $$

Calcul de la Pression

$$ p = \rho_0 \frac{\partial \Phi}{\partial t} = \rho_0 (i\omega)\Phi $$

On fait le rapport $p/v$ :

$$ \frac{p}{v} = \frac{\rho_0 i \omega \Phi}{i k \Phi} = \rho_0 \frac{\omega}{k} $$

Or $c = \omega / k$, donc $p = \rho_0 c \cdot v$

2. Cas de l'onde régressive ($x \searrow$)

Le potentiel est $\Phi = B e^{i(\omega t + kx)}$.

  • La dérivée temporelle reste inchangée ($i\omega$).
  • La dérivée spatiale change de signe ($\partial_x e^{ikx} = ik e^{ikx}$).
$$ v = - (ik)\Phi \quad \text{et} \quad p = \rho_0(i\omega)\Phi $$

Donc $\frac{p}{v} = -\rho_0 \frac{\omega}{k} \implies$ $p = -\rho_0 c \cdot v$

3. Déduction de l'Intensité

L'intensité est le produit $I = p \cdot v$.

En remplaçant $v$ par $p / \rho_0 c$ (en valeur absolue pour la norme) :

$$ I = p \cdot \left( \frac{p}{\rho_0 c} \right) = \frac{p^2}{\rho_0 c} $$

Exercice 6 : Impédance de l'Onde Sphérique

Slide 53

On considère une onde sphérique divergente harmonique définie par le potentiel $\Phi(r,t) = \frac{A}{r} e^{i(\omega t - kr)}$.
1. Exprimer la Pression $p$ et la Vitesse $v$.
2. En déduire l'impédance acoustique spécifique $Z(r) = p/v$.
3. Que se passe-t-il "très loin" de la source ($r \to \infty$) ?

SOLUTION DÉTAILLÉE ▼
1. Calculs de p et v

On rappelle les relations : $p = \rho_0 \frac{\partial \Phi}{\partial t}$ et $v = - \frac{\partial \Phi}{\partial r}$.

A. Pression (Temporel)

Dérivée par rapport à $t$ :

$$ p = \rho_0 \frac{A}{r} \frac{\partial}{\partial t} \left[ e^{i(\omega t - kr)} \right] $$ $$ = \rho_0 \frac{A}{r} (i\omega) e^{i(\omega t - kr)} $$

On reconnaît $\Phi = \frac{A}{r} e^{i(\omega t - kr)}$, donc :

$$ p = i \omega \rho_0 \Phi $$
B. Vitesse (Spatial)

Dérivée par rapport à $r$ (Produit $f \cdot g$) :

On pose $f = \frac{1}{r}$ et $g = e^{i(\omega t - kr)}$.
$\frac{\partial f}{\partial r} = -\frac{1}{r^2}$
$\frac{\partial g}{\partial r} = -ik e^{i(\omega t - kr)}$
$$ \frac{\partial \Phi}{\partial r} = A \left[ \left(-\frac{1}{r^2}\right) e^{i(\omega t - kr)} + \left(\frac{1}{r}\right) (-ik) e^{i(\omega t - kr)} \right] $$

On factorise par $\Phi = \frac{A}{r} e^{i(\omega t - kr)}$ :

$$ \frac{\partial \Phi}{\partial r} = \Phi \left[ -\frac{1}{r} - ik \right] $$

Comme $v = - \frac{\partial \Phi}{\partial r}$, les signes s'inversent :

$$ v = \Phi \left[ \frac{1}{r} + ik \right] $$
2. Expression de l'Impédance Z

On fait le rapport $Z = p / v$ :

$$ Z(r) = \frac{i \omega \rho_0 \Phi}{\Phi \left( \frac{1}{r} + ik \right)} = \frac{i \omega \rho_0}{\frac{1}{r} + ik} $$

Pour simplifier, on multiplie numérateur et dénominateur par $r$ :

$$ Z(r) = \rho_0 c \cdot \frac{ikr}{1 + ikr} $$

(Rappel : on a utilisé $\omega = c \cdot k$)

3. Limite à l'infini (Champ Lointain)

Si on s'éloigne de la source ($r \to \infty$), le terme $ikr$ devient prédominant devant 1 ($ikr \gg 1$).

$$ \lim_{r \to \infty} Z(r) \approx \rho_0 c \cdot \frac{ikr}{ikr} = \rho_0 c $$

Conclusion physique : Près de la source, l'impédance est complexe (déphasage Pression/Vitesse). Loin de la source, la courbure du front d'onde devient négligeable : l'onde sphérique se comporte localement comme une onde plane ($Z \to \rho c$).

Exercice 7 : Réflexion/Transmission en Incidence Normale (1D)

Slide 63 + Notes

On considère une interface plane en $x=0$ séparant deux milieux caractérisés par leurs impédances $Z_1$ et $Z_2$.
Objectif : Retrouver les coefficients de réflexion/transmission en Vitesse ($r_v, t_v$) et en Pression ($r_p, t_p$), puis vérifier la conservation de l'énergie.

DÉMONSTRATION COMPLÈTE ▼

1. Pose du système (Continuités)

À l'interface ($x=0$), il y a continuité de la Pression et de la Vitesse :

$$ \begin{cases} p_i + p_r = p_t & \text{(1) Pression} \\ v_i + v_r = v_t & \text{(2) Vitesse} \end{cases} $$
⚠️ Point Clé : Le signe de la Vitesse

Lien Impédance/Vitesse selon le sens de propagation (voir notes) :

  • Onde Incidente ($\rightarrow$) : $p_i = Z_1 v_i$
  • Onde Transmise ($\rightarrow$) : $p_t = Z_2 v_t$
  • Onde Réfléchie ($\leftarrow$) : $p_r = - Z_1 v_r$ (Car elle recule !)

2. Résolution : Vitesse vs Pression

A. Coefficients en Vitesse

On remplace les $p$ par des $v$ dans (1) :

$$ Z_1 v_i - Z_1 v_r = Z_2 v_t $$

On utilise (2) $v_t = v_i + v_r$ :

$$ Z_1(v_i - v_r) = Z_2(v_i + v_r) $$

On regroupe les termes :

$$ r_v = \frac{v_r}{v_i} = \frac{Z_1 - Z_2}{Z_1 + Z_2} $$

De même pour la transmission :

$$ t_v = \frac{v_t}{v_i} = \frac{2 Z_1}{Z_1 + Z_2} $$
B. Coefficients en Pression

On remplace les $v$ par des $p$ dans (2) :

$$ \frac{p_i}{Z_1} - \frac{p_r}{Z_1} = \frac{p_t}{Z_2} $$

On utilise (1) $p_t = p_i + p_r$ :

$$ \frac{p_i - p_r}{Z_1} = \frac{p_i + p_r}{Z_2} $$

On regroupe les termes :

$$ r_p = \frac{p_r}{p_i} = \frac{Z_2 - Z_1}{Z_2 + Z_1} $$

De même pour la transmission :

$$ t_p = \frac{p_t}{p_i} = \frac{2 Z_2}{Z_2 + Z_1} $$
Remarque fondamentale : $r_v = - r_p$

3. Bilan de Puissance (Intensité)

On définit les coefficients de réflexion ($R$) et transmission ($T$) en énergie.

$$ R = - r_v \cdot r_p = r_p^2 = \left( \frac{Z_2 - Z_1}{Z_2 + Z_1} \right)^2 $$
$$ T = t_v \cdot t_p = \frac{2 Z_1}{Z_1+Z_2} \cdot \frac{2 Z_2}{Z_1+Z_2} = \frac{4 Z_1 Z_2}{(Z_1 + Z_2)^2} $$

Vérification de la conservation de l'énergie ($R+T=1$) :

$$ R + T = \frac{(Z_2-Z_1)^2 + 4Z_1Z_2}{(Z_1+Z_2)^2} = \frac{(Z_2+Z_1)^2}{(Z_1+Z_2)^2} = 1 $$

CQFD.

Exercice 8 : Célérité Complexe ($C_T^*$)

Slide 68

On considère une onde transverse avec atténuation. On définit le vecteur d'onde complexe $k^*$ de deux manières :
1. Par les paramètres physiques : $k^* = \frac{\omega}{C_T} - i\alpha$ (où $C_T$ est la célérité de phase réelle).
2. Par la définition de la célérité complexe : $k^* = \frac{\omega}{C_T^*}$.

Montrer que la partie réelle et imaginaire de $C_T^*$ s'écrivent : $$ C_T^* = \frac{C_T \omega^2}{C_T^2 \alpha^2 + \omega^2} + i \frac{C_T^2 \omega \alpha}{C_T^2 \alpha^2 + \omega^2} $$

DÉMONSTRATION ÉTAPE PAR ÉTAPE ▼

1. Égalité des vecteurs d'onde

On égale les deux définitions de $k^*$ :

$$ \frac{\omega}{C_T^*} = \frac{\omega}{C_T} - i\alpha $$

2. Inversion pour isoler $C_T^*$

On divise par $\omega$ :

$$ \frac{1}{C_T^*} = \frac{1}{C_T} - i\frac{\alpha}{\omega} $$

On met au même dénominateur le terme de droite :

$$ \frac{1}{C_T^*} = \frac{\omega - i \alpha C_T}{\omega C_T} $$

On inverse l'expression :

$$ C_T^* = \frac{\omega C_T}{\omega - i \alpha C_T} $$

3. Séparation Réelle / Imaginaire (Conjugué)

Pour supprimer le $i$ au dénominateur, on multiplie le numérateur et le dénominateur par le conjugué $(\omega + i \alpha C_T)$ :

$$ C_T^* = \frac{\omega C_T (\omega + i \alpha C_T)}{(\omega - i \alpha C_T)(\omega + i \alpha C_T)} $$

Dénominateur : Identité remarquable $(A-iB)(A+iB) = A^2 + B^2$.

$$ D = \omega^2 + (\alpha C_T)^2 $$

Numérateur : On développe.

$$ N = \omega^2 C_T + i (\omega C_T)(\alpha C_T) = \omega^2 C_T + i \omega \alpha C_T^2 $$

4. Résultat Final

On sépare la fraction en deux parties :

$$ C_T^* = \underbrace{\frac{C_T \omega^2}{\omega^2 + \alpha^2 C_T^2}}_{Re(C_T^*)} + i \underbrace{\frac{C_T^2 \omega \alpha}{\omega^2 + \alpha^2 C_T^2}}_{Im(C_T^*)} $$
Note Importante (Slide 68) :

On remarque que $Re(C_T^*) \neq C_T$. La célérité réelle est modifiée par l'atténuation.
Toutefois, si l'atténuation est très faible ($\alpha \to 0$), alors $Re(C_T^*) \approx C_T$.

Exercice 9 : Module de Cisaillement Complexe ($G^*$)

Slide 70

Dans un solide viscoélastique, on définit le module de cisaillement complexe par la relation $G^* = \rho (C_T^*)^2$.
On pose $G^* = G' + iG''$, où :

En utilisant la relation de dispersion $k^* = \frac{\omega}{C_T^*} = \frac{\omega}{C_T} - i\alpha$, démontrer les expressions exactes de $G'$ et $G''$.

DÉMONSTRATION DÉTAILLÉE ▼

1. Expression de $C_T^*$

On part de l'égalité des vecteurs d'onde :

$$ \frac{\omega}{C_T^*} = \frac{\omega}{C_T} - i\alpha = \frac{\omega}{C_T} \left( 1 - i \frac{\alpha C_T}{\omega} \right) $$

On inverse pour isoler la célérité complexe :

$$ C_T^* = C_T \cdot \frac{1}{1 - i X} \quad \text{avec } X = \frac{\alpha C_T}{\omega} $$

2. Calcul de $G^*$

On applique la définition $G^* = \rho (C_T^*)^2$ :

$$ G^* = \rho C_T^2 \left( \frac{1}{1 - i X} \right)^2 = \rho C_T^2 \frac{1}{(1 - i X)^2} $$

On développe le carré au dénominateur :

$$ (1 - iX)^2 = 1 - 2iX + (iX)^2 = (1 - X^2) - 2iX $$

Donc :

$$ G^* = \rho C_T^2 \frac{1}{(1 - X^2) - 2iX} $$

3. Séparation Réelle / Imaginaire

Pour éliminer la partie imaginaire du dénominateur, on multiplie par le conjugué $D^* = (1 - X^2) + 2iX$ :

$$ G^* = \rho C_T^2 \frac{(1 - X^2) + 2iX}{ [ (1 - X^2) - 2iX ] \cdot [ (1 - X^2) + 2iX ] } $$

Le dénominateur devient : $(1-X^2)^2 + (2X)^2 = 1 - 2X^2 + X^4 + 4X^2 = 1 + 2X^2 + X^4 = (1+X^2)^2$.

On obtient l'expression finale séparée :

$$ G^* = \rho C_T^2 \left[ \frac{1 - X^2}{(1+X^2)^2} + i \frac{2X}{(1+X^2)^2} \right] $$

4. Identification

Partie Réelle ($G'$)
$$ G' = \rho C_T^2 \frac{1 - (\frac{\alpha C_T}{\omega})^2}{(1 + (\frac{\alpha C_T}{\omega})^2)^2} $$
Partie Imaginaire ($G''$)
$$ G'' = \rho C_T^2 \frac{2 \frac{\alpha C_T}{\omega}}{(1 + (\frac{\alpha C_T}{\omega})^2)^2} $$