In analisi matematica l'equazione delle onde, conosciuta anche come equazione di d'Alembert[1], è un'equazione differenziale alle derivate parziali iperbolica di grande importanza in diversi campi della fisica, tra cui acustica, elettromagnetismo e fluidodinamica (varianti dell'equazione si trovano anche in meccanica quantistica e relatività generale), descrivendo solitamente la propagazione di un'onda, lineare e non dispersiva, nelle variabili spaziali e temporali, tra cui le onde sonore ed elettromagnetiche. Storicamente il primo problema in cui è stata derivata è stato quello della corda vibrante di uno strumento musicale, studiato da Jean le Rond d'Alembert, Eulero, Daniel Bernoulli e Joseph-Louis Lagrange.
Indice
1L'equazione
2Soluzione
2.1Equazione della corda vibrante
3Equazione scalare in più dimensioni
3.1Equazione in tre dimensioni
3.2Equazione in due dimensioni
3.3Equazione in dimensione arbitraria
3.3.1Dimensioni dispari
3.3.2Dimensioni pari
4Equazione delle onde non omogenea
5Esempi
6Note
7Bibliografia
8Voci correlate
9Altri progetti
10Collegamenti esterni
L'equazione
La forma generale dell'equazione delle onde riguarda una funzione (in generale ) della posizione (in generale ) e del tempo . Si tratta di un'equazione alle derivate parziali iperbolica la cui espressione generale è[2]:
dove rappresenta la velocità di propagazione dell'onda. Per un'onda sonora che si propaga nell'aria la velocità vale all'incirca 330 metri al secondo, mentre per una corda vibrante può assumere valori assai diversi. L'equazione può essere scritta con l'operatore di d'Alembert:
Se l'onda si propaga in un mezzo dispersivo la velocità è dipendente dalla frequenza, e deve essere rimpiazzata dalla velocità di fase, dove indica il numero d'onda:
La funzione incognita esprime l'intensità dell'onda in una particolare posizione al tempo . Per un'onda sonora che si propaga nell'aria, ad esempio, esprime la pressione dell'aria nei diversi punti dello spazio. Per una corda vibrante, invece, esprime lo spostamento fisico della corda dalla sua posizione di riposo. Il simbolo denota l'operatore di Laplace rispetto alla variabile posizione , in generale un vettore. Anche può consistere in una grandezza scalare o vettoriale.
L'equazione pone inoltre in evidenza la diretta proporzionalità tra la concavità della funzione incognita con la sua accelerazione. Un'onda può sovrapporsi ad un altro movimento, ed in tale caso la funzione scalare contiene un fattore di Mach (che ha valore positivo per l'onda che si muove lungo il flusso, e valore negativo per l'onda riflessa).
Soluzione
Lo stesso argomento in dettaglio: Onda (fisica).
L'equazione delle onde può essere scritta come[3]:
e quindi:
Si tratta della somma di due onde che si propagano in verso opposto, come fu mostrato da Jean le Rond d'Alembert[4][5][6].
Definendo ed si ha:
da cui si ricava:
In questo modo l'equazione delle onde assume la forma[7]:
la cui soluzione è:
Si tratta di due onde che si propagano in verso opposto con velocità .
Le funzioni e si determinano a partire dalle condizioni iniziali:
ottenendo la formula di d'Alembert[8][9]:
Se classicamente e , allora , sebbene e possano essere distribuzioni. Ad esempio, se si tratta di funzioni deltiformi la soluzione può essere vista come un impulso che si propaga in una direzione.
Un modo equivalente per giungere alla soluzione si ottiene definendo le variabili:
e considerando l'equazione delle onde:
Si calcolano quindi le derivate:
e le derivate di espresse in funzione di ed :
Inserendo tali espressioni nell'equazione delle onde, tutti i termini si semplificano tranne la derivata mista:
L'ultima equazione implica che:
Allora la soluzione è data dalla somma di e , e ritornando nelle variabili originali si ha:
dove:
e per determinare e si devono imporre le due condizioni iniziali.
L'equazione delle onde nel caso monodimensionale si può derivare nel modo che segue. Si immagini una fila di corpuscoli di massa che sono interconnessi mediante minuscole barrette limitatamente flessibili ciascuna di lunghezza . Le barrette sono caratterizzate da massa trascurabile da una rigidezza (flessionale), cioè di una resistenza alle forze che tendono a fletterla, che viene misurata da . Per questo modello misura la distanza della posizione d'equilibrio della corpuscolo posto in al tempo . L'equazione di moto per il corpuscolo nella posizione è:
Si supponga vi siano di questi oggetti distribuiti in modo uniforme sulla lunghezza . Essi hanno complessivamente massa , mentre la rigidezza totale della fila è . Si può allora scrivere l'equazione precedente nella forma:
Passando al limite per e e si ottiene:
dove è il quadrato della velocità di propagazione in questo caso particolare.
Equazione scalare in più dimensioni
La soluzione per il problema dei valori iniziali dell'equazione in tre dimensioni può essere ottenuta dalla soluzione per un'onda sferica, e tale risultato può essere utilizzato per ricavare la soluzione in due dimensioni. Per spazi di generica dimensione si considerano separatamente il caso di dimensione dispari e pari.
L'equazione delle onde rimane inalterata in seguito alla rotazione delle coordinate spaziali, essendo il laplaciano invariante sotto rotazione: si vuole sfruttare tale simmetria per ottenere una soluzione che dipenda solo dalla distanza radiale dal punto di osservazione. Soluzioni di questo tipo devono soddisfare la relazione[10]:
dove e sono rispettivamente la derivata parziale prima e seconda rispetto a , e lo stesso vale per e . L'espressione può essere scritta come:
dove la quantità soddisfa l'equazione monodimensionale. Vi sono pertanto soluzioni che hanno la forma:
dove e sono funzioni arbitrarie, corrispondenti a due onde che si propagano sfericamente in verso opposto alla velocità .
Si consideri una sorgente che emette a frequenza fissata costante con fase nulla per e con un'ampiezza picco-picco pari a . Detta la distanza dalla sorgente, l'ampiezza dell'onda è data da[11]:
Un'onda di questo tipo, caratterizzata da una sola frequenza di propagazione, è detta monocromatica.
Una somma di onde sferiche è ancora soluzione dell'equazione d'onda, ed in questo modo si può costruire un numero arbitrario di soluzioni. Sia una funzione arbitraria e si supponga che la forma dell'onda sia una Delta di Dirac. Si consideri una famiglia di onde sferiche con centro e sia la distanza radiale da tale punto. Si ha:
e se è una sovrapposizione di onde di questo tipo pesata dalla funzione allora:
Dalla definizione di funzione deltiforme:
dove , e sono le coordinate sulla sfera unitaria . Risulta che è t-volte il valor medio di su una sfera di raggio centrata in :
da cui segue che:
Il valor medio è una funzione pari di , e quindi se:
allora:
che fornisce la soluzione per il problema relativo al valore iniziale.
In accordo con il principio di Huygens-Fresnel, ogni elemento di un fronte d'onda si può considerare formalmente come una sorgente secondaria di onde sferiche in fase con la sorgente primaria e di ampiezza proporzionale a quella dell'onda primaria e all'area dell'elemento. La perturbazione prodotta in un punto dello spazio si può sempre ottenere come sovrapposizione di tutte le onde sferiche secondarie che raggiungono quel punto.
Equazione in due dimensioni
In uno spazio a due dimensioni l'equazione d'onda ha la forma:
Se si tratta come una funzione ambientata in uno spazio di tre dimensioni che è indipendente dalla terza dimensione:
la formula risolutiva in tre dimensioni diventa:
dove e sono le prime due coordinate sulla sfera unitaria, e è l'elemento di superficie sulla sfera. L'integrale può essere scritto come un integrale sul disco con centro e raggio :
Equazione in dimensione arbitraria
Si vuole ricavare la soluzione dell'equazione:
per , con:
Dimensioni dispari
Sia un intero dispari e sia .[12] Si supponga e per . Definendo attraverso la relazione:
si ha che ed in vale la relazione[13]:
Inoltre:
Dimensioni pari
Sia un intero pari e sia .[14] Si supponga e per . Definendo attraverso la relazione:
si ha che ed in vale la relazione[15]:
Inoltre:
Equazione delle onde non omogenea
L'equazione delle onde non omogenea in una dimensione ha la forma:
con condizioni iniziali:
La funzione è detta sorgente poiché descrive l'effetto delle sorgenti di onde sul mezzo nel quali si propagano. Per esempio, in elettromagnetismo la radiazione elettromagnetica ha come termine sorgente una densità di carica e/o di corrente.
Per ottenere la soluzione dell'equazione con condizioni iniziali date si può sfruttare il fatto che essa obbedisce al principio di causalità, ovvero in ogni punto il valore di dipende soltanto da e e dal valore della funzione tra e . Tali quantità sono infatti le uniche presenti nella formula risolutiva di d'Alembert, e la condizione è fisicamente dovuta al fatto che la velocità della luce è la massima velocità di propagazione possibile: e questo implica che l'ampiezza dell'onda in un punto dello spazio e in un certo istante di tempo è relazionata all'ampiezza dall'onda in un punto distante dal primo ad un tempo diverso, non istantaneamente. In termini di calcolo della soluzione questo si traduce nel fatto che in un certo arco di tempo per ogni punto si deve considerare un'area corrispondente che ne è causalmente relazionata. Integrando pertanto l'equazione non omogenea in tale regione:
si ottiene la somma di tre integrali di linea lungo i limiti della regione causalmente connessa. Si ha:
mentre il termine dell'espressione precedente integrato rispetto al tempo si annulla in quanto l'intervallo di tempo è nullo, sicché .
Per i restanti due limiti della regione si nota che è costante, da cui si ottiene . Si ha nuovamente:
ed in modo simile:
Sommando i tre risultati ottenuti ed inserendoli nell'integrale iniziale:
dove si è esplicitato il contorno dell'integrale sulla funzione sorgente. Tale soluzione è valida per ogni scelta di compatibile con l'equazione d'onda, ed i primi due termini sono la formula di d'Alembert soluzione dell'equazione omogenea. La differenza risiede quindi nel terzo termine, l'integrale sulla sorgente.
Esempi
In generale la velocità della propagazione ondosa varia con la frequenza dell'onda, fenomeno chiamato dispersione. Un'altra comune correzione del modello base consiste nel fatto che, nei sistemi realistici, la velocità può dipendere dall'ampiezza dell'onda, cosa che conduce a un'equazione nonlineare:
La presenza di termini non lineari e dispersivi può dar luogo a comportamenti ondulatori particolari, come i solitoni, o la turbolenza d'onda.
In tre dimensioni, ad esempio per studiare la propagazione del suono nello spazio, si considera:
L'equazione delle onde elastiche in tre dimensioni descrive la propagazione delle onde in un mezzo isotropoomogeneo elastico. I materiali solidi in gran parte sono elastici, quindi questa equazione descrive fenomeni come le onde sismiche della Terra e le onde ultrasoniche usate per rivelare difetti nei materiali. Questa equazione è ancora lineare, ma ha forma più complessa di quella delle equazioni presentate in precedenza, in quanto deve rendere conto sia del moto longitudinale che del trasversale:
dove e sono i cosiddetti moduli di Lamé che descrivono le proprietà elastiche del mezzo, esprime la densità, è la funzione sorgente che esprimente la forza che causa il moto e è lo spostamento. Si noti che in questa equazione sia la forza che lo spostamento sono grandezze vettoriali, e pertanto essa viene anche chiamata equazione vettoriale delle onde.
^(FR) Jean le Rond d'Alembert (1747) "Recherches sur la courbe que forme une corde tenduë mise en vibration", Histoire de l'académie royale des sciences et belles lettres de Berlin, vol. 3, pp. 214-219.
^(FR) Jean le Rond d'Alembert (1747) "Suite des recherches sur la courbe que forme une corde tenduë mise en vibration", Histoire de l'académie royale des sciences et belles lettres de Berlin, vol. 3, pp. 220-249.
^(FR) Jean le Rond d'Alembert (1750) "Addition au mémoire sur la courbe que forme une corde tenduë mise en vibration," Histoire de l'académie royale des sciences et belles lettres de Berlin, vol. 6, pp. 355-360.