Passa ai contenuti principali

Analisi matematica \\ Considerazioni sulla divergenza


La divergenza è un operatore vettoriale largamente utilizzato non solo in fisica ma anche in dimostrazioni di teoremi puramente matematici. Si definisce come divergenza di un campo F la struttura div(F)=F. L' operatore nabla è invece definito in un riferimento cartesiano tridimensionale i,j,k come =i(∂/∂x)+j(∂/∂x)+k(∂/∂x). Per il calcolo della divergenza si ricorre al procedimento di calcolo del prodotto scalare (se si opera su una matrice si segue la regola riga per colonna). E' interessante notare che la divergenza applicata ad una struttura matematica tende a "scalarizzarla", ovvero se la applichiamo ad un vettore otteniamo uno scalare,se la applichiamo a un tensore otteniamo un vettore.

Considerazioni fisico-geometriche:
Si è visto che la divergenza rappresenta un prodotto scalare e dunque se lo applichiamo ad un campo vettoriale otterremo tre funzioni scalari che descrivono gli incrementi che un ipotetico campo subisce nello spazio lungo ogni direzione considerata. Le tre funzioni scalari ovviamente sono sommate nell'unica funzione scalare che la divergenza rappresenta.
Nel complesso (come per le semplici derivate) la divergenza rappresenta un incremento,ma proiettato sugli assi di riferimento.
Immaginiamo questi due campi:


(potremmo associare il primo a quello di un fluido in moto laminare che accelera in una condotta e il secondo al moto di un fluido su una lastra che genera taglio oppure un tornado che si avvita su un asse)
Le linee del primo campo sono tutte parallele tra di loro, per evidenziare l'incremento lungo un solo asse e semplificare l'analisi. Se applichiamo una divergenza al primo campo si otterrebbe un valore diverso da zero e positivo (il campo cresce lungo l'asse concorde)
Se invece prendiamo ad esame il secondo campo la variazione di intensità ortogonale non darà contributo alla divergenza bensì al rotore. Ortogonalmente non si ha nemmeno una componente, ma avrei potuto metterla e mantenerla costante per ottenere lo stesso risultato.

Approfondimento: divergenza di un gradiente
Approfondimento: divergenza di un prodotto scalare-vettore
Approfondimento: teorema della divergenza(solo enunciato)

Commenti

Post popolari in questo blog

Matlab \\ Metodo Montecarlo

Il codice richiede di inserire soltanto il numero dei valori in cui si desidera dividere l'intervallo e la funzione. Maggiori saranno gli intervalli del dominio e più preciso sarà il calcolo dell' integrale definito.E' comunque dimostrato che il metodo Montecarlo ha una convergenza molto lenta e richiede, in generale, un numero elevato di intervalli per ottenere una stima precisa. %calcolo di integrali di funzioni in X=[0,1] function[]=metodo_montecarlo() disp('inserire numero di elementi del dominio') t=input('t='); X=linspace(0,1,t); disp('definire la funzione f(X),(x maiuscola)') f=input('f(X)='); y=f'; s=length(X); E=sum(y(:))/s; disp('integrale definito della funzione in X=[0,1],E='); disp(E);

Analisi matematica \\ Teorema degli zeri

Enunciato: Sia una funzione f continua in [a,b] e sia f(a)f(b)<0 allora esiste un punto c∈ (a,b) tale che f(c)=0 . Graficamente si comprende meglio il significato del problema, cioè l'esistenza di almeno un punto di intersezione di f con l'asse delle x se la funzione passa per due punti opposti rispetto all'asse. Dimostrazione: La dimostrazione vale per la funzione raffigurata nell'immagine (crescente con f(a) < f(b) ) ma può essere estesa a qualsiasi caso che rispetti le condizioni del teorema. Utilizzando il metodo di bisezione si divide l'intervallo [a,b] scrivendo il punto medio c 1 =(b-a)/2 . Se con il punto appena trovato vale la relazione f(c 1 )=0 il teorema è dimostrato. Altrimenti bisogna reiterare il processo prendendo ora il punto c 1 e sostituirlo con la scrittura a 1 se f(c 1 ) <0 o b 1 se f(c 1 )>0 (questa è una semplice formalità per poter reiterare il processo similmente al primo passaggio...

Matlab \\ Equazioni concatenate di Allievi

Programma per il calcolo delle sovrappressioni del colpo d'ariete. Questo codice prevede l'inserimento manuale dei dati all avvio del programma. Per creare il vettore dei tempi,come indicato nel programma, è consigliabile usare la funzione "linspace(inizio, fine, lunghezza del vettore)". %h0=carico piezometrico iniziale %n0=valore della funzione di variazione sezione all otturatore al tempo 0 %a=celerità %t_off=tempo della manovra di chiusura %vettore_t=vettore contenenti i tempi a intervallo costante function[]=allievi() disp('definire il carico piezometrico iniziale h0(metri)') h0=input('h0='); disp('definire la funzione apertura otturatore n0(inizio manovra)') n0=input('n0='); disp('definire la celerità di propagazione') a=input('a='); disp('definire tempo della manovra di chiusura(secondi)') t_off=input('t_off='); disp('vettore dei tempi [funzione "linspace(inizio,fine,num...