UNIVERSITÁ DEGLI STUDI DI PADOVA Dipartimento di Ingegneria Industriale TESI DI LAUREA IN INGEGNERIA MECCANICA E MECCATRONICA CURRICULUM MECCANICO LA CAVITAZIONE NELLE MACCHINE A FLUIDO INCOMPRIMIBILE Relatori: Ch.mo Prof. MIRTO MOZZON Ch.mo Prof. RINO A. MICHELIN Laureando: STEVEN TIEPPO Anno Accademico 2013/2014 Introduzione La tesi tratta uno dei limiti principali delle macchine a fluido incomprimibile e dei dispositivi idraulici come le valvole. La cavitazione è quel fenomeno che si genera quando la pressione del liquido, in locazioni specifiche, scende al di sotto della tensione di vapore creando delle bolle che, in seguito, spostandosi in una zona di pressione superiore, implodono generando onde d’urto che possono danneggiare i materiali del dispositivo. Tale fenomeno è caratterizzato inoltre da un rumore caratteristico che si vedrà essere il miglior modo per rilevare l’entità del danno potenziale. Si passeranno in rassegna le tipologie di nucleazione delle bolle di vapore, le loro dinamiche sferiche e non sferiche, nel caso della crescita e del collasso, andando infine a fare una carrellata delle principali forme di cavitazione in grande scala. Si vedranno le tipologie del danno che questi cicli di implosione possono provocare, ed alcuni metodi per arginarli ed eventualmente eliminarli. In alcune occasioni non è possibile evitare l’insorgenza della cavitazione, quindi è necessario monitorarla per riuscire a prevenire i danni più severi. 1 INDICE CAPITOLO 1 - Cambio fase, nucleazione e cavitazione 1.1 Introduzione 5 1.2 Lo stato liquido 5 1.3 Fluidità ed elasticità 6 1.4 Resistenza alla trazione 7 1.5 Cavitazione ed ebollizione 8 1.6 Tipi di nucleazione 9 1.7 Nucleazione omogenea 9 1.8 Confronto con gli esperimenti 11 1.9 Nucleazione eterogenea 12 1.10 Numero dei siti di nucleazione 14 1.11 Contaminazione di gas 15 1.12 Nucleazione nei liquidi in movimento 16 1.13 Effetti viscosi 18 1.14 Incipiente cavitazione 19 1.15 Condizione di incipiente cavitazione 20 1.16 Dimensionamento di incipiente cavitazione 21 CAPITOLO 2 - Dinamiche della bolla sferica 2.1 Introduzione 23 2.2 Equazione di Rayleigh-Plesset 23 2.3 Contenuto delle bolle 25 2.4 Assenza dell’effetto termico 26 2.5 Stabilità delle bolle di vapore/gas 28 2.6 Cavitazione dalla diffusione di massa 31 2.7 Effetti termici sulla crescita 32 2.8 Crescita termicamente controllata 34 2.9 Effetti di convezione 35 2.10 Effetti della rugosità 36 2.11 Perturbazioni non sferiche 37 2 CAPITOLO 3 - Collasso delle bolle 3.1 Introduzione 39 3.2 Collasso della bolla 39 3.3 Collasso controllato termicamente 41 3.4 Effetti termici sul collasso 42 3.5 Forma non sferica nel collasso 43 3.6 Danni della cavitazione 48 3.7 Collasso delle nuvole 52 3.8 Rumore della cavitazione 52 3.9 Luminescenza 56 CAPITOLO 4 - Dinamica delle bolle oscillanti 4.1 Introduzione 59 4.2 Frequenza naturale 59 4.3 Costante politropica effettiva 62 4.4 Termini di smorzamento 64 4.5 Effetti non lineari 65 4.6 Analisi debolmente non lineari 67 4.7 Soglia per la cavitazione transitoria 69 4.8 Rettificazione della diffusione di massa 70 4.9 Forze di Bjerknes 72 CAPITOLO 5 - Flussi cavitanti 5.1 Introduzione 74 5.2 Bolle in movimento 74 5.3 Interazioni bolle/flusso 75 5.4 Osservazioni sperimentali 76 5.5 Cavitazione in larga scala 79 5.6 Cavitazione in vortici 80 5.7 Nuvole di cavitazione 85 5.8 Fogli di cavitazione o cavitazione allegata 85 3 5.9 Lamine cavitanti 88 5.10 Chiusura di cavità 89 Conclusioni 90 Bibliografia e sitografia 92 4 CAPITOLO 1 Cambio fase, nucleazione e cavitazione 1.1 Introduzione In questo capitolo ci si focalizza sui meccanismi di formazione della miscela bifase liquido-vapore, facendo particolare attenzione al processo di creazione delle bolle di vapore in un dato liquido. Si osservano i fondamentali fisici della nucleazione di tali bolle, separando i due campi di interesse: cavitazione ed ebollizione. Una utile approssimazione per distinguere questi due processi è definire la cavitazione come il processo di nucleazione in un liquido, quando la pressione cade al di sotto della tensione di vapore del liquido stesso, mentre l’ebollizione come il processo di nucleazione in un liquido, quando la temperatura è aumentata sopra a quella di saturazione liquido/vapore. Certamente da un punto di vista fisico ci sono poche differenze tra i due processi, ma tali differenze avvengono per via delle diverse complicanze nei flussi cavitanti, da una parte, e nei gradienti di temperatura e negli effetti delle pareti dall’altra. 1.2 Lo stato liquido Ogni discussione sul cambio di fase da liquido a gas o viceversa deve necessariamente essere preceduta da una discussione sullo stato liquido. Questa discussione avviene utilizzando un diagramma di fase che può essere considerato valido per varie sostanze di interesse pratico. La Figura 1.1 mostra un tipico grafico di pressione, p, temperatura, T, e volume specifico, V, dove viene indicato ogni stato della sostanza. Il punto triplo, è quel punto nel diagramma in cui coesistono gli stati solido, liquido e di vapore, e questo indica che la sostanza ha tre configurazioni stabili alternative. La linea di saturazione liquido/vapore (“binodal” nella figura di destra) si estende da questo punto, al punto critico, che viene termodinamicamente definito come il punto al quale il potenziale chimico di due fasi coesistenti deve essere uguale (“critical point”). Guardando il diagramma p-V: considerando di partire da una condizione di liquido puro nel punto A, se si depressurizza mantenendo la temperatura costante, lo stato del liquido varia lungo un’isoterma fino al punto B che è la condizione di liquido saturo, alla pressione del vapor saturo; se sono presenti sufficienti siti di nucleazione il liquido saturo diviene vapor saturo proseguendo da B verso C. Se invece non sono presenti sufficienti siti di nucleazione, la configurazione scende lungo la “theoretical isotherm” fino al punto D che è un punto di equilibrio metastabile. Un liquido allo stato D è detto “in tensione”, e la differenza di pressione tra B e D indica quello stato di tensione. 5 Figura 1.1: Diagramma di fase. Brennen (1995) Si può arrivare allo stato D anche da D’ aumentando la temperatura; la differenza di temperatura tra D e D’ è detto “surriscaldamento”. Caso analogo con F-F’ dove la differenza di temperatura è detta “grado di sottoraffreddamento del vapore”. 1.3 Fluidità ed elasticità Prima di procedere con maggiori dettagli è necessario vedere alcune differenze importanti tra le caratteristiche dello stato liquido e compararle con quelle dello stato solido cristallino e gassoso. La prima differenza tra liquido saturo e vapor saturo sta nella densità, nel liquido rimane praticamente costante e simile a quella del solido fino alla temperatura critica, mentre nel vapore può cambiare da 2 a 5 o più ordini di grandezza, cambiando radicalmente con la temperatura. Il calore specifico del liquido saturo è circa dello stesso ordine di quello del solido, eccetto in un intervallo di temperatura prossima a quella critica, . Di conseguenza anche il movimento termico delle molecole del liquido è simile a quello del solido, e questo implica la formazione di vibrazioni di piccola ampiezza. In molti processi pratici si ha dominante fluida piuttosto che elastica, e per questo si tende a considerare un liquido molto più simile ad un gas anziché ad un solido, ma anche nel solido si ha una componente di fluidità, che è mostrata quando si raggiunge una temperatura di circa il 60-80% di quella di fusione, questa fluidità allo stato solido è nota come creep, fenomeno che avviene per migrazione delle molecole nel reticolo cristallino a causa dell’agitazione termica. Si può allora, sia nel solido che nel liquido, prevedere un certo tempo, tm , per la migrazione di una molecola da una posizione della struttura ad un’altra. Se il tempo, t, associato alle forze applicate a tale struttura, è molto piccolo in confronto a tm allora il comportamento mostrato è di tipo elastico, d’altra parte se t tm si ha fluidità. L’osservazione del tempo, t, diventa molto importante quando il fenomeno viene controllato da eventi come la diffusione delle “vacanze” nella diffusione al creep. Per un 6 periodo di tempo molto lungo è probabile che le vacanze si fondano, creando una sacca di vapore che condurrebbe all’avvio della nucleazione. Si ha anche la possibilità che il liquido possa rimanere in uno stato di tensione (pressione negativa) per un significante periodo di tempo prima che si formi quella sacca, questa abilità del liquido è una chiara manifestazione della sua elasticità. 1.4 Resistenza alla trazione Frenkel (1955) illustrò la potenziale resistenza alla trazione di un liquido puro con il semplice metodo che segue. Considerando due molecole separate da una distanza, s, l’energia potenziale, ϕ, tra le due, dipende dal tipo di legame intermolecolare presente e ha la forma di Figura 1.2. L’equilibrio avviene alla distanza, , che è tipicamente dell’ordine di 10 -10m. La forza di attrazione è la derivata spaziale dell’energia potenziale F = (∂ϕ/∂x) cioè la pendenza della tangente alla curva in un dato x. Figura 1.2: Potenziale intermolecolare. Brennen (1995) Si nota che la forza mostra un massimo in dove il rapporto ⁄ è dell’ordine di 1.1 o 1.2, e a tale distanza si ha un seguente aumento di volume, / 0, di circa un terzo del valore iniziale. L’applicazione di una tensione costante, σ, equivalente alla forza F causerebbe la rottura del legame liquido o solido. Da notare che il gradiente di F è negativo fino a , al quale la pendenza si annulla, positivo aumentando fino a , dove si presenta un flesso (massimo della pendenza), e si annulla proseguendo oltre. Di conseguenza per la forza F è insufficiente per contrastare la tensione σ. Infatti i liquidi e i solidi hanno un modulo di comprimibilità, k, di solito nell’intervallo da 1010 a 1011 kg/ms2, e poiché la pressione è data da p = -k( / 0) segue che la pressione tipica alla quale si rompe il liquido va da -3x109 a -3x1010 Pa. Su tale base si può stimare che un liquido o un solido dovrebbero resistere ad una tensione da 3x104 a 3x105 atmosfere. 2 L’energia elastica accumulata per unità di volume del sistema è data da k /2 0 oppure per sopra da p /2, quindi l’energia per dividere le molecole e vaporizzare deve essere |p |/6 cioè tra 5x108 e 5x109 Pa. 7 Questo risultato è in accordo con l’ordine di grandezza del calore latente di vaporizzazione calcolato per molti liquidi. Si può anche stimare l’ordine di grandezza della assumendo che l’energia cinetica del movimento termico per molecola, k , sia uguale all’energia richiesta per dividere le molecole (con k costante di Boltzmann 1.38x10-23 kgm2/s2K). Prendendo 1030 molecole per metro cubo, questo implica che 30 300 K, che è in accordo con gli ordini di grandezza degli attuali valori. Nonostante questo modello preveda correttamente gli ordini di grandezza di calore latente e di temperatura critica, non riesce a stimare correttamente la tensione alla quale il liquido può resistere. Questo perché la tensione, dipende dalle debolezze interne alla struttura liquida che sono difficilmente calcolabili e quantificabili. 1.5 Cavitazione ed ebollizione La resistenza alla trazione di un liquido si può manifestare nei due modi già accennati: 1. Mantenendo la temperatura costante, diminuendo la pressione, p, al di sotto della tensione di vapore, pV. Il valore della differenza (pV – p) è detto tensione p, e la grandezza alla quale avviene la rottura del legame liquido è detta resistenza alla trazione del liquido pC. Il processo di “rottura” del liquido in questo modo è detto cavitazione. 2. Mantenendo il liquido a pressione costante, aumentando la temperatura, T, ad un valore superiore a quello di saturazione, TS. Il valore della differenza (T –TS) è il grado di surriscaldamento , e il valore al quale avviene la formazione del vapore è detto surriscaldamento critico C. Questo processo è detto ebollizione. Tali valori di tensione e surriscaldamento possono essere relazionati tra loro quando le grandezze sono abbastanza piccole. Dalla relazione di Clausius-Clapeyron nelle condizioni di saturazione si ha: p 1 1 [1.1] con V e L densità di saturazione del vapore e del liquido rispettivamente, mentre L è il calore latente di vaporizzazione. Nella maggior parte dei casi, tranne vicino al punto critico, si ha che la densità del liquido è molto maggiore di quella del vapore, di conseguenza, il termine L-1 nella [1.1] diventa trascurabile, quindi si ha p/ ≈ VL/T quindi integrando: p 8 [1.2] 1.6 Tipi di nucleazione Nelle applicazioni pratiche le debolezze che possono dar luogo alla nucleazione avvengono in 2 modi. I moti termici all’interno del liquido formano dei temporanei vuoti microscopici, che possono costituire i nuclei necessari per la rottura del legame, e la formazione e crescita di bolle macroscopiche. Questa nucleazione è detta omogenea. Nelle pratiche ingegneristiche le situazioni più comuni di punti d’innesco avvengono al confine tra liquido e parete solida del contenitore, o tra liquido e piccole particelle sospese nel liquido stesso, questo tipo di nucleazione è detta eterogenea. Ulteriori debolezze nella massa liquida possono essere, microbolle di gas presenti nelle fessure (cavità) del legame solido, o nelle particelle sospese o semplicemente nel liquido. Nell’acqua, le microbolle di aria sono sempre presenti e impossibili da rimuovere completamente. 1.7 Teoria della nucleazione omogenea In un liquido puro la tensione superficiale, S, è la manifestazione macroscopica della forza intermolecolare che tende a mantenere unite le molecole, e prevenire la formazione di buchi. La pressione del liquido, p, esterna ad una bolla di raggio R, è relazionata alla pressione interna della bolla, pB, dalla seguente relazione: p 2 p [1.3] Se la temperatura è uniforme e la bolla contiene solo vapore, allora la pressione interna, pB, è la pressione di saturazione del vapore pV, quindi la pressione esterna, che dalla [1.3] diventa p = pV – 2S/R, deve essere minore della tensione di vapore per avere una condizione di equilibrio. Di conseguenza, se la pressione esterna fosse, anche solo leggermente, minore di pV – 2S/R, la bolla crescerebbe e aumentando di raggio risulterebbe avere una pressione interna maggiore di quella esterna. Segue che se la massima dimensione dei vuoti presenti è RC (raggio critico), allora la resistenza alla trazione del liquido diventa: 2 p [1.4] Questa relazione è la prima di tre equazioni che costituiscono la teoria della nucleazione omogenea. La seconda equazione che bisogna identificare è quella riguardante l’energia che si deve fornire al corpo liquido, in modo da creare un nucleo o una microbolla della dimensione critica RC. Assumendo che il nucleo critico sia termodinamicamente in equilibrio con l’ambiente circostante dopo la creazione, allora l’energia che bisogna dare al liquido è composta da due fattori: 9 1. Dell’energia per tener conto di quella presente sulla superficie della bolla. Dalla definizione di tensione superficiale quest’importo è dato da S per unità di superficie, per un totale di 4πRC2S. 2. Il liquido deve essere disposto verso l’esterno, e questo implica un lavoro dal o verso il sistema. Tale lavoro è il prodotto tra volume della bolla e differenza di pressione: 4π C3 pC/3 questo è il lavoro fatto dal liquido per ottenere lo spostamento implicato dalla creazione della bolla. L’energia netta che deve essere data per formare la bolla, è: 4π 2 4 π 3 3 p 4 π 3 2 [1.5] eliminando il raggio RC: 3 16π 3 p2 [1.6] L’ultimo passaggio della teoria della nucleazione omogenea è la valutazione del meccanismo di deposito dell’energia, e la probabilità che questa possiede di raggiungere la grandezza WCR nel tempo a disposizione. Nel corpo del liquido completamente isolato da ogni radiazione esterna, ci si riconduce alla valutazione della probabilità che il movimento termico delle molecole possa portare ad una locale perturbazione di energia di quella grandezza. Molte teorie su questa tipologia di nucleazione, relazionano questo lavoro WCR alla energia cinetica tipica delle molecole (vista in §1.4), esprimendo la formula in termini del numero di Gibbs risulta [1.7] k Di conseguenza, un dato numero di Gibbs corrisponde ad una certa probabilità di nucleazione in un dato volume durante un dato tempo a disposizione. Altre formule, invece, relazionano il tasso di nucleazione, J, definito come numero degli eventi di nucleazione in un’unità di volume per unità di tempo, con il numero di Gibbs, assumendo la forma seguente: [1.8] dove JO è il fattore di proporzionalità che per Blander e Katz (1975), assume la tipica forma successiva: 1 2 2 ( ) πm [1.9] dove N è il numero di densità del liquido, espresso come molecole/m3, ed m è la massa di una molecola. 10 1.8 Confronto con gli esperimenti Il tasso di nucleazione J varia con la temperatura, ed è quindi importante identificarne l’ordine di grandezza per poter capire le osservazioni sperimentali sul fenomeno. Dalla [1.6], [1.7] e [1.8] si può ricavare il valore di pC come: p ( 16π 3k ln ( 1 2 3 , [1.10] * Questa formula si usa per calcolare la resistenza alla trazione del liquido, data la temperatura, e conoscendo la variazione della tensione superficiale con la temperatura, e altre proprietà del fluido come la definizione di un tasso di nucleazione critico J. Si nota che l’effetto più importante causato dalla temperatura è quello della variazione di S3 al numeratore, dato che S è approssimativamente lineare con la temperatura, andando a zero al punto critico, si ha che la resistenza alla trazione è fortemente funzione della temperatura vicino al punto critico. Mentre, ogni dipendenza dalla temperatura nel termine JO è praticamente trascurabile, poiché si trova all’argomento di un logaritmo. Si divide in due la discussione dei risultati sperimentali, sopra e sotto la temperatura alla quale la pressione spinodale è circa zero. Questa temperatura è empiricamente data da T/TC = 0.9, per temperature comprese tra (0.9-1)TC, la resistenza alla trazione del fluido è abbastanza modesta, poiché il raggio critico, a tali temperature, è piuttosto grande. Ad esempio per 1bar di tensione si ha un RC di circa un micrometro, questo implica che le microcontaminazioni di particelle, o le microbolle hanno un piccolo effetto sull’esperimento in questo intervallo di temperatura. Figura 1.3: Osservazione sperimentale della vita media di etere dietilico surriscaldato a quattro pressioni diverse, in funzione della temperatura di saturazione. Skripov (1974) 11 Figura 1.4: Limite di surriscaldamento per cinque liquidi diversi, comparato con la linea del liquido spinodale da cinque equazioni di stato diverse. Brennen (1995) Nella Figura 1.3 viene riportata l’osservazione sperimentale della vita media in secondi, plottata con la temperatura di saturazione di un’unità di volume di liquido surriscaldato (nello specifico etere dietilico). Sempre nell’intervallo TC > T > 0.9TC, si nota che finché J-1 è minore di 5 secondi si ha comportamento di nucleazione omogenea, mentre la radicale deviazione dell’esperimento si ha per J-1 > 5s causata principalmente da radiazioni che inducono nucleazione ad un più piccolo valore del surriscaldamento. La figura mostra anche quanto debolmente il limite del surriscaldamento dipenda dal tasso di nucleazione critico. Finché le linee sono per lo più verticali, si può ottenere, dai risultati sperimentali, un massimo valore del surriscaldamento possibile senza la necessità di stipulare uno specifico tasso di nucleazione critico. La Figura 1.4 mostra i dati del limite di surriscaldamento per 5 tipi diversi di liquido. Per molti liquidi in questo intervallo di pressioni positive, il massimo surriscaldamento possibile è accuratamente previsto dalla teoria della nucleazione omogenea. Se si considerano temperature più basse di 0.9TC la teoria della nucleazione omogenea funziona ancora in modo soddisfacente, però aumentano le incertezze, poiché questa teoria implica tensioni pC sempre più grandi e quindi raggi critici sempre più piccoli. Di conseguenza ci deve essere qualche altro tipo di nucleazione che prende il sopravvento, causando la rottura a tensioni molto minori di quelle previste dalla omogenea. 1.9 Nucleazione eterogenea Nella nucleazione omogenea si sono considerati vuoti microscopici di raggio R, i quali crescono causando la rottura, quando la pressione del liquido si abbassa al di sotto del valore critico pV – 2S/R, con tensione data dalla [1.3]. Ora invece si considerano un numero di analoghe situazioni all’interfaccia solido/liquido come illustrato in Figura 1.5. L’angolo di contatto nell’intersezione liquido/vapore/solido è denotato da θ. Nel caso di piano a superficie idrofoba, la resistenza alla trazione è data da 2 sinθ/ con R dimensione massima del vuoto, quindi tende a zero quando θ tende a π. Nell’altro caso, con superficie idrofila, la resistenza alla trazione è simile a quella della nucleazione omogenea finché le dimensioni massime dei vuoti sono comparabili. 12 Figura 1.5: Vari modelli di nucleazione eterogenea. Brennen (1995) Perciò, si potrebbe concludere che la presenza di una superficie idrofoba, causerebbe nucleazione eterogenea, conferendo una resistenza alla trazione molto minore. In scala microscopica, le superfici non sono uniformi o piane, perciò bisogna considerare gli effetti di superfici locali di diverse geometrie. Ad esempio, la cavità conica di Figura 1.5 è usata di solito per esemplificare gli effetti della geometria di superficie. Se metà dell’angolo al vertice della cavità è chiamato α, è chiaro che la resistenza alla trazione è nulla quando θ α + π/2 piuttosto che θ → π. Inoltre, se θ > α + π/2 è chiaro che la bolla di vapore crescerà fino a riempire la cavità, a pressione superiore a quella di saturazione del vapore. Considerando quindi l’intervallo di superfici geometriche microscopiche, non sorprende che le sacche di vapore crescano all’interno di particolari cavità superficiali a pressione vicina a pV, soprattutto quando la superficie è idrofoba. Viene spontaneo porsi il quesito su come si formi la prima sacca di vapore, è plausibile che minuscole sacche di gas estraneo siano assorbite dalla superficie del solido durante la sua formazione. Quindi, la presenza di questi gas estranei e la geometria della superficie, promuovono la crescita e l’aspetto macroscopico delle bolle di vapore. Allora essi sono detti siti di nucleazione, e con più si riduce la pressione maggiori sono i siti che diventano attivi, divenendo capaci di formare e rilasciare bolle nel corpo del liquido. Questi eventi possono essere visti, in modo simile, durante l’ebollizione dell’acqua, dove all’inizio le bolle si formano solo in certi siti vicini alle pareti del contenitore, e proseguendo con l’ebollizione si aumenta il numero di siti attivi, quindi la densità dei siti di nucleazione è funzione della temperatura. 13 1.10 Numero dei siti di nucleazione In un serbatoio in ebollizione, il liquido più caldo si trova a contatto con la parete solida riscaldata, quindi tutti i siti di nucleazione importanti stanno su quella superficie. Per quantificare il processo di nucleazione è necessario definire una funzione della densità di distribuzione sulla superficie per questi siti di nucleazione, N(RP), con RP raggio del sito di nucleazione e dove N(RP)dRP è il numero di siti con grandezza tra RP e RP+dRP per unità di area. Bisogna inoltre conoscere l’intervallo di grandezze utilizzato nelle operazioni ad un dato surriscaldamento . Tutte le dimensioni maggiori di un certo raggio RP* sono allora eccitate da una tensione di β / RP* con β costante circa pari all’unità, questo corrisponde ad un surriscaldamento critico β T * Quindi, il numero di siti per unità di superficie n( specifico surriscaldamento è dato da n ∫ β [1.11] ), portato nelle operazioni con uno [1.12] La Figura 1.6 illustra questo effetto, sulla sinistra sono riportate le misure del numero di siti attivi per unità di superficie, n, per una superficie di rame lucidato, e tre liquidi diversi. Sulla destra si ha il grafico tracciato usando la [1.11] con β = 2, che dimostra la veridicità della [1.12] per la particolare superficie. Se si differenziano i nuclei chiamando, quelli sulla parete solida nuclei di superficie, e quelli all’interno della massa fluida nuclei di flusso libero, si può dire che è molto più complicato identificare il carattere dei nuclei di flusso libero, ad esempio come possano esistere microbolle piene di gas estraneo al liquido per una certo intervallo di tempo, piuttosto che il carattere dei nuclei di superficie. Non è possibile valutare separatamente il numero delle particelle solide e delle microbolle presenti, nonostante entrambi possano agire come siti di nucleazione per la cavitazione, le microbolle sono più prontamente crescenti e osservabili macroscopicamente. Un metodo per contare solo questi nuclei, consiste nel risucchiare il liquido attraverso un venturimetro. I nuclei cavitano per la bassa pressione della gola nel venturimetro, e possono essere contati a patto che la concentrazione sia abbastanza piccola, in modo che gli eventi siano separati nel tempo. Inoltre, in questo modo, si nota che la concentrazione dei nuclei può essere ottenuta come funzione del livello di pressione nella gola, conoscendo la portata del fluido. Tali accorgimenti sono conosciuti come “metri di suscettibilità alla cavitazione” e tendono ad essere limitati ad una concentrazione inferiore ai 10 cm-3. 14 Figura 1.6: Dati sperimentali sul numero dei siti di nucleazione attivi per unità di superficie, nel caso di rame lucidato. Griffith e Wallis (1960) 1.11 Effetti della contaminazione dei gas Virtualmente ogni liquido contiene gas disciolti, ed è impossibile eliminarli completamente da un dato volume di liquido. Per esempio, ci vogliono settimane di deaerazione per ridurre la concentrazione di aria in una certa quantità d’acqua sotto a 3ppm (la saturazione a pressione atmosferica è di circa 15ppm). Se le bolle di nucleazione contengono del gas, allora la pressione nelle bolle è la somma della pressione del gas, pG, e di quella del vapore, pV. Quindi l’equilibrio della pressione nel liquido risulta, dalla [1.3], p = pV + pG – 2S/R quindi, segue che la resistenza alla trazione (dalla [1.4]) è pC = 2S/R - pG, cioè tale resistenza diminuisce con la presenza dei gas disciolti. Inoltre, se la concentrazione dei gas disciolti è sufficiente da portare grandi valori di pG, allora la resistenza alla trazione può risultare negativa, e le bolle, in queste condizioni, crescono ad una pressione del liquido molto più alta della tensione di vapore. Ci si riferisce alle circostanze alle quali il liquido non è saturo con il gas alla pressione di mantenimento, in teoria non possono esistere bolle di gas in equilibrio, in un liquido insaturo di tale gas, esse dovrebbero dissolversi e sparire, causando un notevole aumento della resistenza alla trazione del liquido. Nonostante questo sia vero, che degassazione o trattamenti ad alte pressioni incrementino la resistenza, l’effetto non è così grande come ci si aspetterebbe. Per tale accorgimento ci sono tre spiegazioni plausibili: 1. Harvey (1944) disse che la bolla esiste in una particella o in una fessura della superficie e persiste perché la sua geometria è tale che la superficie libera mostra una curvatura altamente convessa vista dal fluido, cosicché la tensione superficiale sopporti l’alta pressione del liquido. 15 2. Possibilità della continua produzione di nuclei, per effetto della radiazione cosmica 3. Fox e Herzfeld (1954) ipotizzarono la presenza di una pelle organica che desse alla superficie della bolla un’elasticità sufficiente per resistere alla alta pressione. Inizialmente meno plausibile delle altre due, ma ora maggiormente accettata perché alcuni sviluppi tecnologici (microscopi elettronici), mostrano che piccole quantità di contaminazioni nel liquido possono generare grandi effetti di elasticità superficiale. 1.12 Nucleazione nei liquidi in movimento La più comune forma di cavitazione si ha nel moto fluido dei liquidi, dove gli effetti idrodinamici si hanno in regioni del fluido dove la pressione cade sotto alla tensione di vapore del fluido stesso. Reynolds (1873), fu il primo che tentò di spiegare il comportamento insolito delle eliche navali a velocità angolari elevate. Ipotizzò la possibilità del trascinamento di aria nelle scie delle pale dell’elica, fenomeno che oggi viene chiamato ventilazione. Sembra però che non avesse considerato la passibilità di avere vapore nelle scie, e lasciò a Parsons (1906) lo studio di questo fenomeno, e lo sviluppo del primo esperimento sulla cavitazione. Si considerano ora delle ipotesi: la singola fase di flusso di un liquido Newtoniano di densità costante, L, un campo di velocità ui(xi), e un campo di pressione p(xi), in tutti questi campi di flusso è conveniente definire una velocità di riferimento U∞ ed una pressione di riferimento p∞. Nel flusso esterno attorno al corpo solido queste due sono convenzionalmente la velocità e pressione del flusso uniforme a monte del corpo. Definendo un coefficiente di pressione come xi 2 p xi – p∞ ∞ 2 [1.13] ci sono delle posizioni, xi*, nel flusso, dove CP e p sono minime, e quel valore del coefficiente di pressione viene denotato per convenienza con CPmin che è un numero negativo. Gli effetti viscosi nel flusso sono caratterizzati da LU∞l/μL = U∞ /νL (numero di Reynolds) dove μL è la viscosità dinamica e νL è la viscosità cinematica ed l è la lunghezza caratteristica. Si tratta ora di stabilire cosa accade quando in un dato flusso la pressione totale diminuisce o la velocità del flusso aumenta, in modo che ad un certo punto la pressione si avvicini alla tensione di vapore del liquido alla temperatura di riferimento T∞ . Si definisce, per tale proposito, il numero caratteristico di cavitazione, σ come σ p∞ p ( 1 2 2 ∞ ∞) [1.14] Ogni portata di flusso, in cavitazione o no che sia, ha un certo valore di σ. L’avvio della cavitazione avviene ad un valore detto numero caratteristico di incipiente cavitazione 16 σi. Per il momento si trascurano le difficoltà pratiche dell’osservazione dell’incipiente cavitazione. In un flusso ipotetico di un liquido che non può resistere ad alcuna tensione, e nel quale le bolle di vapore appaiano istantaneamente quando la pressione raggiunge il valore pV, è chiaro che σi min [1.15] e il numero di cavitazione può essere accertato con osservazioni o misure sul flusso in singola fase. Per spiegare questo si considera la nucleazione, di nuclei a flusso libero, come il movimento lungo una linea di flusso contenente (xi*) come in Figura 1.7, nella quale per σ > -CPmin la pressione è maggiore della tensione di vapore, per σ -CPmin i nuclei trovano una pressione pari a pV solo per un momento infinitesimo, mentre ne caso in cui sia σ > -CPmin i nuclei sono ad una pressione inferiore a quella di vapore per un certo periodo di tempo finito. Ci sono due fattori per cui σi è diverso da -CPmin: - il primo motivo è che la nucleazione potrebbe non avvenire a p = pV; in un liquido degassato la nucleazione potrebbe aver bisogno di una tensione positiva, pC, quindi la nucleazione richiederebbe un numero caratteristico di incipiente cavitazione minore rispetto all’opposto del coefficiente di pressione minimo, cioè σi = -CPmin - 2 pC/ LU∞2 . In un liquido con grandi quantità di gas contaminante, pC potrebbe essere negativo (come visto in §1.11) e dare quindi un numero caratteristico maggiore di -CPmin. - il secondo motivo è che la crescita di un nucleo fino ad una grandezza osservabile richiede un certo tempo sotto le condizioni di p < pV - pC. L’effetto del tempo di permanenza causa un’osservazione di σi < -CPmin - 2 pC/ LU∞2. Come si vedrà più avanti il numero caratteristico di cavitazione può dipendere, inoltre, anche dalla temperatura del liquido. Figura 1.7: Schematica distribuzione di pressione su una linea di flusso. Brennen (1995) 17 1.13 Effetti viscosi Fino ad ora si è parlato di condizioni ideali con portata uniforme. Quando gli effetti viscosi vengono considerati, CPmin diventa funzione del numero di Reynolds, e si vorrebbe osservare una dipendenza da questo anche nel numero caratteristico di cavitazione per un flusso costante. Fino ad ora si è assunto che il flusso e le pressioni fossero laminari e costanti, tuttavia nell’ingegneria si ha a che fare con flussi non solo turbolenti, ma anche non costanti, che verificano distacco libero e forzato di vortici. Questo aspetto ha un’importante conseguenza per l’incipiente cavitazione, perché al centro del vortice la pressione può essere significativamente minore di quella media nel flusso avendo una maggiore velocità del liquido, e quindi la misura di -CPmin susciterebbe la pressione media minore, mentre la cavitazione potrebbe avvenire prima, in un vortice transitorio nel quale la pressione fosse molto più bassa di quella media minore. Diversamente, il tempo di permanenza tenderebbe a causare un più alto valore di σi di quello che ci si aspetterebbe, causando anche la sua variazione in dipendenza con Re; per distinguerlo dall’effetto di Re su CPmin, lo si chiama effetto di turbolenza. Quindi ci sono diverse ragioni per cui σi potrebbe essere diverso da -CPmin calcolato sulla conoscenza della pressione in una portata di liquido monofase: Presenza di resistenza alla trazione → diminuisce σi Tempo di permanenza → riduce σi Presenza di gas contaminanti → aumenta σi Effetti viscosi costanti dovuti alla dipendenza di CPmin a Re, possono causare la dipendenza del numero di incipiente cavitazione come funzione di Re 5. Presenza di turbolenze → incremento di σi 1. 2. 3. 4. Quindi per questi fattori il criterio di σi = -CPmin non può essere usato, perciò ci sono dei parametri che devono essere controllati in condizione di incipiente cavitazione quali: Il numero caratteristico di cavitazione σ Il numero di Reynolds Re La temperatura media del fluido T∞ Qualità del liquido inclusi numero e natura dei nuclei di flusso libero, l’ammontare dei gas disciolti e la turbolenza del flusso libero Qualità del solido, superfici di confine, incluse la rugosità, porosità e popolazione nella cavità. 18 1.14 Incipiente cavitazione Il primo effetto che si illustra è quello dell’incertezza nella stima della resistenza alla trazione del liquido. È molto difficile caratterizzare e praticamente impossibile rimuovere da un corpo liquido tutte le particelle, microbolle e le varie contaminazioni da gas che interessano la nucleazione. Tali inclusioni possono causare delle sostanziali differenze nel numero caratteristico di cavitazione in “strutture” diverse, ma anche nelle stesse strutture trattando il liquido in modo differente. La Figura 1.8 mostra la variazione di σi in diversi esperimenti sulla stessa sagoma in diverse parti del mondo. I dati ottenuti da Keller (1974) per σi in flussi intorno a corpi emisferici sono riportati in Figura 1.9, dove l’acqua è stata trattata in modi diversi per farle acquisire differenti concentrazioni di nuclei. Come si poteva immaginare, l’acqua con maggior concentrazione di nuclei mostra un maggior numero caratteristico di cavitazione. Poiché i nuclei di cavitazione sono cruciali per la comprensione del fenomeno di incipiente cavitazione, è ormai riconosciuto che il liquido in ogni esperienza di quella condizione deve essere monitorato, misurando il numero dei nuclei presenti. Figura 1.8: Numero caratteristico di incipiente cavitazione misurato per la stessa forma assialsimmetrica in un tunnel ’a qua in iv rs parti l mon o. in gr n ohnsson (1966), ohnsson (1969). È da notare che la stessa cavitazione è fonte di nuclei, poiché l’aria dissolta nel liquido tende ad uscire dalla soluzione a basse pressioni e a contribuire con una pressione parziale all’interno delle bolle macroscopiche della cavitazione. Quando la bolla si porta in una regione con maggiore pressione e il vapore condensa, rimane una piccola bolla d’aria che si ridiscioglie molto lentamente. Si è osservato che dopo pochi minuti di funzionamento con un corpo cavitante, nella sezione di lavoro, le bolle di cavitazione crescono rapidamente in numero e si diffondono per tutta la struttura in considerazione, fino ad arrivare all’osservazione di un flusso bifase. 19 Figura 1.9: Istogrammi della popolazione dei nuclei in acqua di rubinetto trattata e non trattata, e i corrispondenti numeri di incipiente cavitazione su forme emisferiche di tre diverse diametri, 3cm, 4.5cm, e 6cm. Keller (1974) Uno degli effetti della cavitazione sulla popolazione dei nuclei in una struttura è che il numero caratteristico di cavitazione al quale la cavitazione sparisce (σd), quando la p viene aumentata, può essere diverso da quello al quale appare quando la pressione cala. Il primo termine, σd, è detto numero caratteristico di cavitazione desinente. La differenza tra σd e σi è detto fenomeno di isteresi della cavitazione. Una complicanza è quella della visione dell’apparizione della cavitazione, poiché non è sempre possibile da vedere e quando è visibile la manifestazione è particolarmente soggettiva, invece, il rumore prodotto dal fenomeno è più semplice da misurare, e le misure sono maggiormente ripetibili, e quantificabili. Se il liquido è abbastanza pulito con solo nuclei molto grandi, la tensione alla quale il liquido può resistere sta ad indicare la minima pressione che si dovrebbe avere per cadere sotto alla tensione di vapore e far iniziare la cavitazione. Così σi sarebbe molto più piccolo di -CPmin, e questo indica che la qualità dell’acqua e dei suoi nuclei potrebbe causare un numero caratteristico di cavitazione incipiente più grande o più piccolo di quel valore trovato nella [1.15]. 1.15 Condizione di incipiente cavitazione Si considerano prima le caratteristiche di incipienza di una singola aletta idrodinamica con l’angolo d’ingresso variabile. I dati di Kermeen (1956) ottenuti per una NACA 4412 sono riportati in Figura 1.10, dove si nota che per angoli d’ingresso positivi, la regione di bassa pressione in cui si verifica la condizione di incipiente cavitazione si trova sulla superficie di aspirazione, mentre per angoli d’ingresso negativi, il fenomeno si sposta 20 sulla superficie di mandata, inoltre, come l’angolo aumenta, in entrambe le direzioni, il valore di -CPmin aumenta, aumentando di conseguenza il valore del numero caratteristico di incipiente cavitazione. Dati tipici sull’incipiente cavitazione per una singola aletta, sono riportati in Figura 1.11, con dati per tre diverse dimensioni dell’aletta, a velocità differenti, plottati con il numero di Reynolds, nella speranza di ottenere una curva singola; il fatto che ciò non avvenga dimostra che è presente una dimensione con effetto separato dal numero di Reynolds. Per completare la lista dei fattori che influenzano la cavitazione, è necessario parlare di rugosità della superficie e del livello di turbolenza nella portata; questi due parametri sono legati in qualche modo, poiché la rugosità influenza il livello di turbolenza. La rugosità però può anche influenzare il flusso, ritardando la separazione della strato limite, andando a modificare la pressione e velocità in modo più globale. La turbolenza, come già detto, provoca vortici al cui interno la pressione potrebbe essere minore di quella media, promuovendo la cavitazione, ma può alterare la pressione globale cambiando la posizione di separazione del flusso. Figura 1.10: Numero di incipiente cavitazione per una NACA 4412. Kermeen (1956) 1.16 Dimensionamento di incipiente cavitazione Quando si dimensiona un dispositivo idraulico di grandi dimensioni, il metodo utilizzato è quello della similitudine, studiando le caratteristiche di un modello in scala rispetto al prototipo e successivamente trasferendo i dati al prototipo stesso. Per quanto riguarda il dimensionamento di incipiente cavitazione, cambiando le dimensioni del dispositivo, durante il trasferimento dei dati (modello-prototipo), non solo si cambia l’effetto del 21 tempo di permanenza, ma anche il numero di Reynolds, inoltre i nuclei diventano di dimensione diversa rispetto al modello. Per recuperare le condizioni del modello, bisogna cambiare il livello di pressione che potrebbe alterare i nuclei contenuti, variando inoltre anche la rugosità del prototipo rispetto al modello. Un’altra questione è quella di riuscire a prevedere il fenomeno della cavitazione in un liquido, sulla base dei dati ottenuti su un altro. Nella letteratura, ci sono grandi quantità di dati sull’acqua, mentre sono scarsi per altri tipi di liquidi, quindi passando ad altri liquidi il dimensionamento viene eseguito per tentativi. Fino ad ora si sono trattate le condizioni per l’incipienza della cavitazione, ed una volta che diventa stabile, ciò che accade è meno sensibile a fattori come il contenuto dei nuclei, e quindi il dimensionamento può procedere con maggior confidenza rispetto a quello dell’incipienza. Figura 1.11: um ro i avitazion sin nt p r tr iv rs im nsioni ll’al tta i oukowski in funzione del numero di Reynolds. Holl e Wislicenus (1961) 22 CAPITOLO 2 Dinamiche della bolla sferica 2.1 Introduzione Dopo aver considerato la formazione delle bolle, si procede con l’identificazione delle sequenze dinamiche della crescita e del collasso della bolla. Il comportamento di una bolla in un dominio di liquido fermo, con temperatura uniforme lontano dalla bolla si esamina per primo. La situazione di simmetria sferica fornisce il caso più semplice che rivela un numero importante di fenomeni. Figura 2.1: Schema di una bolla sferica in un liquido. Brennen (1995) 2.2 Equazione di Rayleigh-Plesset Considerando una bolla sferica di raggio R(t) dove t è il tempo, in un dominio infinito di liquido, le cui, temperatura e pressione, lontane dalla bolla sono T∞ e p∞(t). La temperatura T∞ si assume come una semplice costante, poiché il riscaldamento del liquido è uniforme, i gradienti di temperatura sono eliminati a priori e non si prendono in considerazione fonti di calore interne o radiazioni. Dall’altra parte la pressione p∞(t) si assume nota e controllata. Mentre densità del liquido, L, e viscosità dinamica, μL, si considerano costanti. Si assume, inoltre, che il contenuto della bolla sia omogeneo e che la temperatura, TB(t), e pressione pB(t) nella bolla siano sempre uniformi. Come indicato in Figura 2.1 la posizione radiale all’interno del liquido è denotata come la distanza r dal centro della bolla, la pressione p(r,t), la velocità radiale verso l’esterno u(r,t) e la temperatura T(r,t) all’interno del liquido sono così assunte. La conservazione della massa richiede la soddisfazione della seguente formula: 23 F t r2 u r,t [2.1] dove la F(t) è relazionata al raggio R(t), dalle condizioni cinematiche della superficie della bolla. Nel caso ideale di assenza di deflusso di massa attraverso tale superficie è chiaro che u(R,t) = dR/dt, e quindi: 2 F t [2.2] t Questa è una buona approssimazione anche quando l’evaporazione o la condensazione avvengono all’interfaccia. Per dimostrarlo si considera una bolla di vapore dove il volume del tasso di produzione di vapore deve essere uguale all’incremento della dimensione della bolla 4π 2dR/dt, in più il tasso di evaporazione della massa dev’essere V(TB)4π 2dR/dt dove V(TB) è la densità del vapor saturo alla temperatura della bolla. Questo a sua volta deve eguagliare il flusso di massa di liquido verso l’interno, relativo all’interfaccia dato da V(TB)( / t)/ L, ma poiché V ≪ L si ha che l’approssimazione [2.2] è adeguata, e quindi si continua ad usarla. Assumendo un fluido newtoniano le equazioni di Navier-Stokes per il movimento nella direzione r, dopo la sostituzione con la [2.1] e dopo integrazione danno: p p∞ 1 F r 1 F2 2 r4 [2.3] Per completare questa analisi è necessario costruire una condizione dinamica sulla superficie della bolla. A tale proposito si considera un volume di controllo che consiste in una piccola lamina come in Figura 2.2. La forza netta su questa lamina in direzione radiale uscente per unità di superficie è (σrr)r=R + pB - 2S/R e poiché σrr = -p + 2μL(∂u/∂r) si ha: σrr r p p 4μ r 2 [2.4] t In assenza di deflusso di massa attraverso la superficie questa forza deve essere nulla, e sostituendo il valore di (p)r=R dalla [2.3] con F = R2dR/dt si ottiene la generalizzazione dell’equazione di Rayleigh-Plesset per le dinamiche della bolla: (p t p∞ t ) 2 4ν t2 2 t 2 3 ( * 2 t [2.5] Dato p∞ (t) si ha un’equazione che può essere risolta per trovare R(t) con pB(t) noto. 24 Figura 2.2: Porzione di superficie della bolla sferica. Brennen (1995) 2.3 Contenuto delle bolle In aggiunta all’equazione di Rayleigh-Plesset è necessario considerare il contenuto della bolla. Per rimanere in generale si assume che la bolla contenga una certa quantità di gas con pressione parziale pG0 e grandezza R0 a temperatura T∞ . Trascurando il trasferimento di gas dal o al liquido, la pressione della bolla risulta pB(t) = pV(TB) + pG0(TB/T∞)(R0/R)3. Rimane da determinare la TB(t), questo non è sempre necessario, perché sotto alcune condizioni la differenza tra la temperatura della bolla e quella del fluido è trascurabile, ma in altre condizioni questo può condizionare le dinamiche della bolla; la differenza di temperatura porta a diversi valori della tensione di vapore pV(TB) rispetto alla situazione in assenza di tali effetti termici, alterando il tasso di crescita o di collasso della bolla. Sostituendo ora il valore della pressione della bolla nella [2.5] si ottiene la forma generale della Rayleigh-Plesset: (p ∞ -p∞ t ) + (p -p ( ∞) + p 0 ( *( ∞ 0 3 * 2 t2 + 4ν t + 2 + 2 3 ( * 2 t [2.6] Dove il primo termine a sinistra è la tensione istantanea (o termine guida {1}) determinato dalle condizioni lontane dalla bolla. Il secondo termine è riferito all’effetto termico {2}, e più avanti si vedrà che le dinamiche della bolla possono essere previste in dipendenza della grandezza di questo termine. Quando la è piccola è conveniente usare l’espansione di Taylor fino al primo ordine, per valutare L e L alla temperatura nota T∞. Il grado con il quale la temperatura della bolla, devia da quella del liquido, può avere un effetto sulle dinamiche della bolla. La determinazione di TB - T∞ necessita di 2 passi: 1. È richiesta la soluzione dell’equazione di diffusione del calore per determinare la distribuzione di temperatura, T(r,t), all’interno del liquido. 2. È richiesto un bilancio di energia per la bolla. Supponendo che il calore fornito all’interfaccia dal liquido sia utilizzato tutto per la vaporizzazione, si può valutare la produzione della massa di vapore e relazionarla all’incremento di volume della bolla. 25 La natura del problema dell’effetto termico è chiara, il termine termico nella relazione di Rayleigh-Plesset richiede una relazione tra la differenza di temperatura (TB - T∞) e il raggio della bolla. Mentre, il bilancio di energia richiede una relazione tra la variazione della temperatura con la distanza dal centro della bolla (r=R) e il raggio stesso. La relazione finale è dunque tra la differenza di temperatura e la variazione della temperatura con la distanza dal centro della bolla, e necessita della soluzione dell’equazione di diffusione del calore. Questa soluzione è difficoltosa, per via dell’evidente non linearità nell’equazione, della quale non esiste una risoluzione analitica esatta. Plesset e Zwick (1952) approssimarono una utile soluzione di questa, confinando tale soluzione ai casi in cui lo spessore dello strato limite termico, δT, attorno alla bolla, sia piccolo in confronto al raggio della bolla. Sotto tali condizioni, il termine termico linearizzato nella Rayleigh-Plesset diventa, utilizzando la soluzione Plesset-Zwick: ∞ n ∞ ∞ 1 * (n 2* t [2.7] Dove C(n) costante vicino all’unità, e R* valore del raggio a regime. Mentre, il parametro termodinamico 2 2 ( ∞) [2.8] 0.5 ∞α 2 si vedrà essere di importanza cruciale per il comportamento dinamico della bolla. 2.4 In assenza dell’effetto termico Si considerano, prima, alcune caratteristiche delle dinamiche della bolla in assenza di effetti termici significanti. Questo tipo di comportamento dinamico della bolla è detto “inerzialmente controllato”, per distinguerlo da quello “termicamente controllato”. Sotto queste circostanze, la temperatura del liquido è assunta uniforme, e il secondo termine della Rayleigh-Plesset è nullo. Inoltre, assumendo che il comportamento del gas nella bolla sia politropico, vale pG = pG0(R0/R)3K con K costante che quando vale 1 si ha comportamento isotermo, e quando vale γ si ha comportamento adiabatico. Con le precedenti assunzioni più la simmetria sferica la Rayleigh-Plesset diventa p ∞ p∞ t p 0 ( 0 3 ̈ * 3 2 ( ̇) 2 4ν ̇ 2 [2.9] Sono richieste le condizioni iniziali e nel contesto di flusso cavitante è appropriato assumere che le microbolle di raggio R0 siano in equilibrio a t = 0 nel fluido a pressione p∞(0) cosicché pG0 = p∞(0) - pV(T∞) + 2S/R0 e che la velocità dR/dt = 0 in t = 0. Una soluzione tipica per la [2.9] sotto tali condizioni e con la pressione che prima cala sotto al valore di p∞(0) e poi ritorna al suo valore originale è mostrata in Figura 2.3. Le caratteristiche generali di questa soluzione sono caratterizzate dalla risposta della bolla, che passa attraverso delle regioni di bassa pressione che riflettono la forte nonlinearità della [2.9]. 26 Figura 2.3: ipi a soluzion ll’ quazion i ayl igh-Plesset per una bolla sferica. Il nucleo entra in una regione di bassa pressione al tempo adimensionale 0 e ritorna alla sua dimensione iniziale al tempo adimensionale 500. Brennen (1995) La crescita è piuttosto regolare e la dimensione massima si verifica con la pressione minima. Il processo di collasso è molto diverso, la bolla collassa in modo catastrofico e questo viene seguito da successive ridelimitazioni (o ricostruzioni) e collassi. In assenza di effetti dissipativi come la viscosità queste ridelimitazioni (rebounds) continuerebbero senza attenuazioni. Nel caso aviscoso, dopo l’integrazione e l’applicazione delle condizioni iniziali (dR/dt=0), si ottiene: 2 ( ̇) p*∞ ) 2(p 3 3 0 ] 3 [1 2p 3 1 0 [ 3 0 3 3 0 ] 3 2 [1 2 0 ] 2 [2.10] Considerando prima il comportamento caratteristico per la crescita della bolla, che si mostra quando p∞* < p∞(0), la [2.10] mostra che il tasso di crescita per R R0 è dato da: 2p ̇ →* 3 p*∞ 1 2 + [2.11] Dopo un periodo iniziale di accelerazione la cui durata, tA, può essere stimata da t [ 2 2 0 (p 3(p∞ 0 27 p*∞ ) ] 2 p*∞ ) [2.12] la seguente velocità dell’interfaccia è relativamente costante. L’equazione 2.11 rappresenta la crescita esplosiva della bolla, il cui volume aumenta con t3. Il comportamento della bolla al collasso, avviene quando si ha l’incremento della pressione da p∞ a p∞* , in questo caso, quando R ≪ R0, la [2.10] risulta: ̇ → ( 3 0 2 * [ 2(p*∞ p ) 2p 2 3 3 0 0 ( 1 0 3 * 1 2 1 ] [2.13] La maggior parte dei movimenti nel collasso della bolla, sono così rapidi che il comportamento del gas si assume adiabatico. Per una bolla con un sostanziale contenuto di gas, la velocità di collasso data dalla [2.13] non viene raggiunta, e la bolla oscilla ad un nuovo, ma più piccolo raggio di equilibrio. D’altra parte in una bolla con poco gas, la velocità in ingresso continua ad aumentare, finché l’ultimo termine in parentesi quadra non raggiunge una grandezza comparabile con gli altri termini. La velocità di collasso, in seguito, diminuisce e si raggiunge una dimensione minima data da 1 0[ min p 1 1 p*∞ 3 0 p ] 3 1 [2.14] 0 e la bolla si ridelimita. Se la pressione del gas è bassa, anche il raggio minimo è piccolo. La pressione e la temperatura al raggio minimo risultano le massime e sono date da pmax = pG0[(K - 1)(p∞* - pV + 3S/R0)] K/(K-1) e Tmax = T0[(K - 1)(p∞* - pV + 3S/R0)]. In assenza di tensione superficiale e contenuto di gas, Rayleigh (1917), ottenne il tempo tTC richiesto per il collasso totale da R = R0 a R = 0: t 0.915 ( * p∞ 2 0 p 1 2 + [2.15] I risultati ottenuti, sono pratici per la crescita della bolla, risultano invece fuorvianti nel collasso. Si vedranno, infatti, due assunzioni tra quelle fatte fino ad ora, che possono essere violate durante il collasso, quella di incomprimibilità del liquido e quella della simmetria sferica della bolla. 2.5 Stabilità delle bolle di vapore/gas È importante sottolineare che le condizioni di equilibrio pV - p∞ + pGE - 2S/RE = 0 non sempre rappresentano un equilibrio stabile ad un raggio RE con una pressione parziale del gas pGE. 28 Si considera una piccola perturbazione nella dimensione della bolla da R = RE ad R = RE(1 + ϵ) con ϵ ≪ 1, e il responso risultante dalla Rayleigh-Plesset. Si devono distinguere i due casi possibili: i) La pressione parziale del gas durante la perturbazione rimane pressoché costante ii) La massa del gas e la temperatura della bolla rimangono le stesse Da un punto di vista pratico in i) la perturbazione è generata in un tempo sufficientemente lungo da acconsentire una adeguata diffusione di massa nel liquido, cosicché la pressione parziale del gas sia mantenuta al valore iniziale. La Rayleigh-Plesset per i due casi diventa: ̈ 4ν ̇ 3 2 ̇ 2 ϵ 2 [ 3n p ] [2.16] con n = 0 in i) e n = 1 in ii) Si nota che il termine a destra ha lo stesso segno della perturbazione ϵ nel caso in cui 2S/RE > 3nKpGE, e questo indica che l’accelerazione e/o la velocità del raggio della bolla hanno lo stesso segno della perturbazione, di conseguenza, si ha sempre equilibrio instabile finché si devia da RE. D’altra parte si ha equilibrio stabile quando la disuguaglianza diventa 2S/RE < 3nKpGE. Nel caso i) si ha n = 0 quindi si ha sempre equilibrio stabile. Nel caso ii) una bolla in equilibrio stabile richiede la soddisfazione della seguente disequazione m 4 π 3 p 2 3 [2.17] 3 dove mG è la massa del gas nella bolla, e KG è la costante del gas. Per una data massa di gas esiste una dimensione critica della bolla RC data da: 9 m [ 8π ] 1 2 [2.18] Tutte le bolle di raggio minore del raggio critico possono esistere in equilibrio stabile, mentre quelle di raggio maggiore al critico sono instabili. Questa dimensione critica potrebbe essere raggiunta da una diminuzione della pressione ambiente fino al valore critico p∞ : p∞ p 4 8π [ 3 9 m spesso chiamata pressione soglia di Blake (1949). 29 1 2 ] [2.19] Figura 2.4: Raggio di equilibrio stabile e instabile della bolla in funzione di diverse masse di gas nella bolla. Daily e Johnson (1956) Nella Figura 2.4 sopra, è rappresentato il caso isotermo (K = 1) dove le linee continue rappresentano le condizioni di equilibrio di una bolla di raggio RE in funzione della tensione (pV - p∞), per diverse masse di gas nella bolla e tensioni superficiali. Il raggio critico della massa corrisponde al massimo della curva. La linea tratteggiata è il luogo dei valori di RC = RE di equazione (pV - p∞) = 4S/3RE. La regione alla destra rappresenta l’equilibrio instabile, mentre alla sinistra quello stabile. Considerando di partire dalla condizione di pressione media maggiore di quella di vapore e assumendo che la massa di gas nella bolla sia costante, il raggio di equilibrio aumenta con la diminuzione di p∞. La bolla passa in zone di equilibrio finché non si raggiunge la pressione critica che corrisponde al massimo della curva, dove si ha raggio critico. A questo punto ogni diminuzione della pressione corrisponde ad una crescita esplosiva della cavitazione, inoltre è chiaro che la tensione critica di un liquido dovrebbe essere data da 4S/3R piuttosto che 2S/R come nel capitolo 1, poiché le condizioni di equilibrio stabile non esistono nell’intervallo 4S/3R < pV - p∞ < 2S/R. In più si nota, sempre dalla Figura 2.4, che per ogni tensione sotto a quella critica, esistono due stati di equilibrio, uno stabile e l’altro instabile. Ipotizzando la bolla nell’equilibrio stabile, essa è soggetta ad oscillazioni di pressione di grandezza sufficiente da superare momentaneamente il raggio critico; a questo punto la bolla cresce esplosivamente senza limite. È chiaro che tutti i nuclei di grandezza maggiore di quella critica diventano instabili, crescendo improvvisamente e cavitando, mentre quelli più piccoli della dimensione critica reagiscono passivamente. Sebbene la risposta della bolla sia dinamica e la pressione locale del liquido attorno alla bolla p∞ cambi continuamente, si può anticipare che la dimensione critica dei nuclei 30 è data approssimativamente da 4S/3(pV - p∞)*, dove (pV - p∞)* è una rappresentazione della tensione nella zona di bassa pressione. Più basso è il valore della pressione locale attorno alla bolla, più piccolo è il raggio critico, e maggiore è il numero di nuclei attivati. Un esempio quantitativo dell’effetto è rappresentato in Figura 2.5, dove il risultato deriva dall’integrazione della Rayleigh-Plesset per bolle di forma assialsimmetrica. In figura è mostrata la massima dimensione delle bolle come funzione delle dimensioni 2 dei nuclei originali, per un tipico numero di Weber, Wb LRHU∞ /S di 28000 dove U∞, velocità di flusso libero e RH, raggio della forma, le curve sono date per quattro diversi numeri tipici di cavitazione. Per σ < 0.5 le curve hanno un tratto verticale corrispondente ad una certa dimensione critica dei nuclei. Si nota inoltre che indipendentemente dalla dimensione iniziale dei nuclei, tutti quelli instabili crescono approssimativamente alla stessa dimensione massima. Questo perché, sia il tasso di crescita sia il tempo necessario per tale crescita, sono relativamente indipendenti dalla dimensione originale dei nuclei. Figura 2.5: Dimensione massima, RM, a cui una bolla di cavitazione cresce in accordo con la Rayleigh-Plesset, come funzione della dimensione dei nuclei originali, R 0, e del numero aratt risti o i avitazion σ, n l flusso attorno a una forma assialsimm tri a i dimensione RM. Ceccio e Brennen (1991) 2.6 Cavitazione dalla diffusione di massa In molte delle circostanze considerate si assume che gli eventi accadano molto rapidamente, per il significante trasferimento di massa, di gas contaminanti, tra la bolla e il liquido. Quindi l’assunzione della costanza della massa deve essere riconsiderata. Uno dei problemi basilari è il fatto che ogni microbolla con gas, presente nel substrato liquido, dovrebbe dissolversi se la pressione fosse sufficientemente elevata. L’esperienza è però contraria a questo comportamento, e le microbolle persistono anche quando la p∞ dell’ambiente liquido è portata a varie atmosfere per un tempo 31 prolungato. Una spiegazione possibile è perché l’interfaccia è immobilizzata da effetti di contaminazione superficiale, oppure perché la bolla è incastonata in una particella solida che inibisce i processi di solubilizzazione del gas. 2.7 Effetti termici sulla crescita Bisogna esaminare il regime di validità dell’analisi precedente, cioè l’assenza dell’effetto termico. È necessario valutare la grandezza del termine [2.7]. Si esamina prima il caso della crescita della bolla; il tasso di crescita dato dalla [2.11] è costante, quindi nel caso caratteristico di pressione costante i termini della RayleighPlesset (eccetto quello in esame {2}) sono tutti circa costanti o diminuiscono in grandezza con il progredire del tempo. Inoltre un tasso di crescita costante corrisponde al caso n = 1 e il termine termico linearizzato per piccoli valori di (T∞ - TB) risulta t rmin t rmi o 1 * 1 [2.20] t2 sotto queste condizioni, anche se il termine termico è inizialmente trascurabile, guadagna in grandezza relativamente agli altri termini e alla fine influenza la crescita in modo principale. Allo stesso modo bisognerebbe aggiungere l’assunzione di un spessore di strato limite termico, δT, piccolo rispetto al raggio, e solo in circostanze di crescita molto lenta tale assunzione sarebbe violata. Figura 2.6: alor l param tro t rmofisi o, , om funzion Brennen (1995) 32 l rapporto / C per diversi liquidi. Usando la [2.20] si può definire un tempo critico, tC1 (primo tempo critico), durante il quale il termine termico raggiunge l’ordine di grandezza degli altri termini della Rayleigh-Plesset, t p 1 p*∞ 2 [2.21] che come si vede non dipende solo dalla tensione ma anche da ( ∞), che è una quantità puramente termofisica, funzione solo della temperatura del liquido. Riferendosi alla [2.8] si nota che tale termine termofisico cambia molto in un liquido con la temperatura che varia dal punto triplo al punto critico, poiché è proporzionale a ( V/ L)4, quindi il primo tempo critico varia di molti ordini di grandezza. Alcuni valori di , per un certo numero di liquidi, sono tracciati in Figura 2.6 come funzione della riduzione di temperatura T/TC, e in Figura 2.7 come funzione della tensione di vapore. Ad esempio, un flusso cavitante (acqua) con tensione di circa 10000 Pa a 20°C ha un valore di di circa 1 m/s3/2 con tC1 = 10s, che è molto più grande del tempo di crescita della bolla, quindi la crescita avviene senza ostacoli da parte del termine termico. Se invece fosse riscaldato a circa 100°C si avrebbe 1000 m/s3/2 con un tC1 = 10μs. In questo caso tutte le crescite sarebbero termicamente controllate. 2 Figura 2.7: Valore del parametro termodinamico, , come funzione della tensione di vapore (kg/ms ) per diversi liquidi. Brennen (1995) 33 2.8 Crescita termicamente controllata Quando il primo tempo critico viene superato, chiaramente la relativa importanza dei vari termini nell’equazione di Rayleigh-Plesset cambiano, i più importanti diventano il termine guida {1} e il termine termico {2}, la cui grandezza risulta molto maggiore di quella del termine inerziale ({4} nella [2.6]). Quindi se la tensione (pV - p∞*) rimane costante, allora la soluzione di forma [2.7], per il termine termico, deve avere n = 1/2. Il comportamento asintotico risulta allora: * p*∞ p ∞ 1 ( * 2 [2.22] Conseguentemente, al progredire del tempo, i termini inerziale, viscoso, gassoso, e di tensione superficiale della [2.6] diminuiscono di importanza, e dopo il primo tempo critico il tasso di crescita diminuisce incrementando il raggio con t1/2 invece che con t. Un esempio di questa inibizione termica sulla crescita della bolla si vede in Figura 2.8. Figura 2.8: Osservazione sperimentale della crescita di tre bolle di vapore, in acqua surriscaldata a 103.1°C in funzione della crescita prevista dalla teoria di Plesset-Zwick. Dergarabedian (1953) Si nota che l’assunzione di Plesset-Zwick, dello strato limite piccolo rispetto al raggio della bolla, rimane valida. Quando la crescita è causata da una diminuzione della pressione cosicché p∞(t) cambi sostanzialmente durante la crescita, l’approssimazione della 2.22 diventa più complicata, risultando necessario risolvere: la Rayleigh-Plesset, l’equazione di diffusione del calore e quella dell’energia, simultaneamente. Per il controllo sulla crescita con il termine termico, molti termini della [2.6] diventano trascurabili, giustificando la semplificazione pV(TB) = p∞(t). 34 Un fattore che potrebbe influenzare le dinamiche di crescita controllata termicamente, può essere quello del se il liquido all’interfaccia si trova in equilibrio termico con il vapore nella bolla. In molte analisi, si assume che la temperatura del liquido all’interfaccia sia la temperatura del vapor saturo nella bolla, ma questo potrebbe non essere corretto per l’elevato tasso di evaporazione. In Figura 2.9 si ha il confronto tra equilibrio apparente dei dati analitici di Theofanous et al. (1969), Jones e Zuber (1978), Cha e Henry (1981) e il risultato sperimentale di Hewitt e Parker (1968). Figura 2.9: Dati sulla crescita di una bolla di vapore in azoto liquido di Hewitt e Parker (1968) confrontati con quelli analitici di Theofanous et al. (1969), Jones e Zuber (1978), Cha e Henry (1981). Brennen (1995) 2.9 Effetti di convezione Un altro modo in cui il tasso di calore trasferito all’interfaccia può variare è attraverso la convezione causata dal moto relativo tra la bolla e il liquido. Tale valorizzazione del trasferimento di calore è rappresentata dal numero di Nusselt, Nu, definito come rapporto tra il tasso di calore reale scambiato e il tasso di calore trasferito per conduzione. Pertanto, in questo caso, il termine termico nella RayleighPlesset è funzione anche di Nu. Si definisce il numero di Peclet come /αL dove W è la tipica velocità relativa di traslazione della bolla, tale numero è in relazione con Nu nei casi comuni dei problemi di trasferimento di calore per conduzione dove Nu = 1 per Pe ≪ 1 e Nu = Pem per Pe 1. Si considerano gli effetti qualitativi per due possibili movimenti traslatori nella crescita della bolla con R = R* t n: 35 1) Spinta idrostatica: la velocità relativa in assenza di trascinamenti viscosi è proporzionale a gt 2) Bolla che cresce su una parete solida: la velocità relativa è approssimata dalla derivata temporale di R, ottenendola proporzionale a R * t n-1 Quindi, i numeri di Peclet sono rispettivamente: * 1) 2 gt n-1 α ( * ) t 2n-1 2) α [2.23] Considerando il caso di crescita controllata inerzialmente con n = 1, l’effetto di trasferimento di calore convettivo si ha solo per ≥ 1 per t > tC2 dove t 2 1 4 α2 * (p p*∞ )g2 + [2.24] Segue che l’aumento del trasferimento convettivo di calore è presente nella crescita controllata inerzialmente se tC2 < tC1. Ciascuno dei due moti della bolla presenta una temperatura al di sotto della quale ci si aspetta che il numero di Peclet raggiunga l’unità prima del primo tempo critico. La questione è su ciò che accade dopo, poiché chiaramente l’effetto termico che avrebbe inizio a tC1 sta per essere modificato dal maggior trasferimento di calore. Quando Pe > 1 il termine termico, nella Rayleigh-Plesset, cresce come t1/2/Nu cioè, nel caso dei due movimenti come t1/2-2m e t1/2-m. Come in molti problemi pratici di trasferimento del calore per convezione, m = 1/2 e le inibizioni termiche della crescita vengono annullate e il controllo inerziale continua indefinitamente. Nell’altro caso in cui tC2 > tC1 l’effetto di trasferimento di calore convettivo può influenzare la crescita controllata termicamente. Con n = 1/2 il numero di Peclet per il movimento indotto dalla spinta idrostatica diventa unitario in un tempo t 3 * 2 3 α (p p*∞ )g + [2.25] Nel caso di crescita sulla parete solida Pe rimane al valore minore dell’unità che ha raggiunto a tC1. 2.10 Effetti della rugosità superficiale Un altro fenomeno che potrebbe influenzare lo scambio termico all’interfaccia, durante la crescita della bolla, è lo sviluppo di un’instabilità sull’interfaccia stessa. Se la superficie della bolla diventa rugosa, l’incremento dell’area di superficie e l’instabilità nel movimento del liquido vicino ad essa può portare ad un sostanziale miglioramento dello scambio termico. L’effetto si traduce in un ritardo dell’alterazione dovuta al termine termico. 36 Degli studi, condotti da Shepherd e Sturtevant (1982) e da Frost e Sturtevant (1986), esaminarono la rapida crescita di bolle di nucleazione vicino al limite di surriscaldamento, trovando tassi di crescita più grandi rispetto alla crescita controllata termicamente. Le foto in Figura 2.10 rivelano che la superficie delle bolle è rugosa e irregolare, e la valorizzazione dello scambio termico causato da questa rugosità è probabilmente responsabile del maggior tasso di crescita trovato empiricamente. É stata, inoltre attribuita una somiglianza della rugosità della bolla, all’instabilità della fiamma. Figura 2.10: Fotografie della rapida crescita di una bolla in una goccia di etere surriscaldato sospesa in glicerina. La bolla è scura mentre la goccia è chiara e trasparente. Frost e Sturtevant (1986) 2.11 Perturbazioni non sferiche Fino ad ora si è assunto che la bolla rimanga sferica durante la crescita e il collasso, cioè stabile a distorsioni non sferiche. Ci sono però circostanze in cui questo non è vero. Se l’inerzia del gas nella bolla si assume trascurabile, allora l’ampiezza a(t) di una distorsione armonica sferica di ordine n, (n > 1) è governata da 2 a t2 3 a t t [ n 1 2 t2 n 1 n 1 n 2 3 ]a 0 [2.26] nella quale bisogna conoscere R(t). Le circostanze più instabili si hanno quando dR/dt < 0 e d2R/dt2 ≥ 0, queste condizioni sono raggiunte appena prima della ridelimitazione di una cavità collassante. Dall’altra parte le circostanze più stabili si hanno con dR/dt > 0 e d2R/dt2 < 0, caso che accade quando, durante la crescita, si sta raggiungendo la dimensione massima. Il fatto che i coefficienti della [2.26] non rimangano costanti nel tempo causa una divergenza rispetto all’equivalente di instabilità studiata da Rayleigh-Taylor che non consideravano l’effetto termico. Plesset e Mitchell (1956) esaminarono il caso particolare di una bolla di vapore/gas, inizialmente in equilibrio, soggetta alla variazione della pressione, trascurando gli effetti termici e viscosi e l’effetto della massa fissa all’interno della bolla. Nel caso del collasso, una caratteristica importante è che l’accelerazione risulta essere maggiore di zero, e questo dà origine all’instabilità. In un caso reale l’accelerazione iniziale positiva è di durata limitata, quindi il problema è di stabilire se l’instabilità ha 37 tempo sufficiente per una crescita significante. Riscrivendo in termini di y = R/R0 la [2.26], si ha che a(y) deve soddisfare: 2 y a y2 y y a y n 1 y a 3 y 0 [2.27] Dove i termini A(y), B(y) e C(y) sono funzione di β1 = 2S/R0(pV - p∞) e β2 = pG0/(pV - p∞), dove il primo è l’effetto della tensione superficiale e il secondo l’effetto del contenuto di gas. Un valore positivo di (pV - p∞)(1 - β1 + β2) implica la crescita della bolla da t = 0, mentre un valore negativo implica il collasso. In Figura 2.11 si riportano tipiche integrazioni numeriche della [2.27]. In caso di crescita della bolla si nota un picco dell’ampiezza del disturbo e una successiva diminuzione. Per alcuni parametri di β1 e β2 esiste una particolare armonica sferica (n = nA) dove si raggiunge la massima ampiezza del disturbo. Ricapitolando, si può vedere che l’accelerazione positiva, all’inizio della crescita della bolla, è instabile per perturbazioni armoniche sferiche di ordine n, d’altra parte si può concludere che la restante fase di crescita, durante la quale dR/dt > 0 e dR2/dt2 < 0, è stabile a tutte alle stesse perturbazioni. Quindi, se durante la fase di accelerazione è disponibile un tempo inadeguato per la crescita di perturbazioni, allora la bolla rimane imperturbata durante tutta la crescita. Figura 2.11: s mpi i r s ita ll’ampi zza, a, i un istur o armoni o sf ri o ( i or in n) sulla superficie della bolla di cavitazione per due tipiche scelte della tensione superficiale e dei parametri l gas ont nuto β1 β2. Brennen (1995) 38 CAPITOLO 3 Collasso delle bolle 3.1 Introduzione In questo capitolo si focalizza l’attenzione sulle dinamiche della bolla in fase di collasso considerando le conseguenze delle bolle di cavitazione piene di vapore. Valutando soprattutto le tipologie dei danni che questo fenomeno provoca sui materiali del dispositivo in cui si manifesta. 3.2 Collasso della bolla Il collasso della bolla è un fenomeno particolarmente importante per il rumore e per il danno ai materiali che provoca, causati dalle elevate velocità e dalle pressioni e temperature massime che possono entrare in gioco in questa fase. Nel capitolo precedente si era tentato di valutare la grandezza di queste entità (dal §2.4). Sotto ad un certo numero di assunzioni compresa quella della rimanenza sferica della bolla, si ottengono i massimi valori di pressione, temperatura e danno materiale. La deviazione dalla sfericità, quindi, può diffondere la focalizzazione del collasso e ridurre le pressioni e le temperature massime. Quando una bolla di cavitazione cresce, da un piccolo nucleo a molte volte la sua dimensione originale, il collasso comincia una volta raggiunto il raggio massimo RM, con una pressione parziale dei gas, pGM, che è molto piccola. In un flusso cavitante il raggio massimo è spesso dell’ordine di cento volte la dimensione del nucleo originale R0. Segue che se la pressione originale del gas nei nuclei è di 1bar il valore della pressione massima all’inizio del collasso potrebbe essere di circa 10-6bar, se (p∞* - p∞(0)) ≈ 0.1bar seguirebbe una pressione massima di circa 1010bar con una temperatura massima che potrebbe essere 4x104 volte quella ambiente (T∞). Molti fattori come la diffusione del gas dal liquido alla bolla e l’effetto della comprimibilità del liquido, attenuano questi valori. In Figura 3.1 si mostrano dei dati ottenuti da Herring (1941), e Gilmore (1952), che assunsero la bolla priva di gas incondensabili, trascurarono gli effetti termici mantenendo una pressione costante e introdussero il numero di Mach per il moto del collasso. Si nota quindi che, nel caso ideale, il numero di Mach della superficie della bolla, aumenta al diminuire del raggio della bolla. La linea “incompressible” corrisponde al caso in cui la comprimibilità del liquido è stata trascurata nell’equazione del moto 2.13 , da notare che la compressibilità tende a diminuire la velocità del collasso. Quando la bolla contiene qualche gas incondensabile o quando il termine termico diventa importante, la soluzione diventa più complessa poiché la pressione nella bolla non è più costante. 39 Figura 3.1: Numero di Mach della superficie della bolla, -(dR/dt)/c, in funzione del raggio per una differenza di pressione di 0.517 bar. Per il metodo analitico di Herring (1941) e Gilmore (1952). Brennen (1995) Sotto queste circostanze è utile trovare il modo di incorporare l’effetto di comprimibilità del liquido in una versione modificata della Rayleigh-Plesset. Keller e Kolodner (1956) proposero la seguente modifica in assenza del termine termico, viscoso e di tensione superficiale: (1 2 1 t * 2 t 3 (1 2 1 3 2 t *( t * (1 1 t * 1 [p p∞ p (t *] p t Dove pc(t) rappresenta la parte variabile della pressione del liquido, nel centro della bolla in assenza della bolla. Si suggerirono altre equazioni, ma tutte erano ugualmente valide, e lineari col numero di Mach (dR/dt)/c, e si dimostrò che le modifiche alla [2.6] erano tutte accurate finché il numero di Mach è dell’ordine di 0.3, mentre a numeri superiori, le equazioni devono essere risolte numericamente. Tuttavia, finchè c’è del gas presente per rallentare il collasso, l’importanza principale della comprimibilità del liquido non è l’effetto che ha sulle dinamiche della bolla (debole), ma il ruolo nella formazione di onde d’urto durante la ridelimitazione che segue il collasso. La Figura 3.2 presenta un esempio dei risultati per la distribuzione della pressione nel liquido prima (sinistra) e dopo (destra) del momento di dimensione minima, cioè prima e dopo la ridelimitazione. Il grafico di destra mostra la propagazione di un impulso di pressione o di urti, lontani dalla bolla, seguendo la dimensione minima. Concludendo si può dire che l’impulso di pressione presenta un’attenuazione mentre si propaga lontano dalla bolla, raggiungendo un picco nell’ampiezza pp che è dato all’incirca da pp 100 40 p∞ r [3.2] Figura 3.2: Risultati tipici per la distribuzione della pressione nel liquido prima e dopo il collasso. Con valori -3 di p∞ 1 ar, γ 1.4 pr ssion inizial n lla olla i 10 bar. Brennen (1995) Questa formula dà l’ordine di grandezza dell’impulso di pressione che urta contro una superficie solida a pochi raggi di distanza; ad esempio se si ha una pressione media nel liquido di 1bar con una distanza r = RM si ha un picco dell’impulso di pressione dell’ordine di 100bar. Tutte queste analisi assumono la simmetria sferica della bolla. 3.3 Collasso controllato termicamente Bisogna ricordare che potrebbero verificarsi effetti termici nelle prime fasi del collasso delle bolle, nello stesso modo in cui avvengono nella fase di crescita. Considerando una bolla di raggio R0, inizialmente a riposo al tempo t = 0, in un liquido a pressione p∞, il collasso inizia con l’incremento della pressione dell’ambiente liquido a p∞*. Dalla Rayleigh-Plesset il moto iniziale in assenza degli effetti termici ha la forma 1 p t2 [3.3] 0 dove pc è definita come p [ p*∞ p 2 2 0 3m 3 0 8π ∞ 5 0 ] [3.4] Si definisce inoltre un tempo critico, tC4, necessario per lo sviluppo di effetti termici significanti come 41 t ( 4 2 0 3 [3.5] * Procedendo con l’analisi, se tC4 ≪ tTC dove tTC è il tempo tipico per il collasso, allora il collasso termicamente controllato inizia presto in questo processo. Tali condizioni si verificano se (p 2 1 3 0 )2 ( * [3.6] ≪1 In tal caso, il movimento iniziale è effettivamente dominato dal termine termico ed è dato della seguente forma, p 0 2 0 1 [ 4π 2 1 ∫0 y/(y 1 1 y 1 )2 ]t2 [3.7] dove tra parentesi quadra si ha un termine di ordine circa unitario. Se la disequazione 3.6 è violata, l’effetto termico comincerà a diventare importante solo più tardi nel processo di collasso. 3.4 Effetti termici sul collasso Anche se gli effetti termici sono trascurabili in molte fasi del collasso, essi giocano un ruolo molto importante nello stadio finale, quando il contenuto della bolla è altamente compresso dall’inerzia del liquido. La pressione e la temperatura che sono previste nel gas all’interno della bolla durante il collasso sferico sono molto alte. Poiché il tempo trascorso è molto piccolo (μs) sembrerebbe un’approssimazione ragionevole assumere che i gas non condensabili nella bolla si comportino adiabaticamente. Attraverso un esperimento, Benjamin (1958) ottenne una temperatura massima del gas di circa 8800K nel centro della bolla. Ma nonostante il piccolo tempo di questo picco, è stato dimostrato che lo scambio termico tra liquido e gas è importante per l’alto gradiente di temperatura e per la piccola distanza in gioco. Altri valori arrivarono intorno a 6700K e 848bar, per una frazione di microsecondi. Altri effetti che possono essere importanti sono le interdiffusioni di gas e vapore nella bolla, che possono causare un accumulo di gas non condensanti all’interfaccia, e pertanto creare una barriera attraverso la quale il vapore deve diffondersi per condensare sull’interfaccia. 42 3.5 Forma non sferica nel collasso Si considera ora il collasso di una bolla che contiene principalmente vapore. Come già fatto si distinguono due fasi importanti del moto, escludendo l’accelerazione iniziale verso l’interno. 1. La forma asintotica del collasso in cui la velocità è proporzionale ad R-3/2 che avviene a significanti compressioni del gas 2. La fase di ricostruzione (o ridelimitazione) nella quale l’accelerazione cambia segno e assume un valore positivo molto grande. Le caratteristiche stabili di queste due fasi sono molto diverse. Il calcolo di Plesset e Mitchell (1956) mostrò che la bolla è solo leggermente instabile durante la prima fase in cui l’accelerazione è negativa, d’altra parte è chiaro che la bolla può divenire molto instabile per disturbi non sferici durante la seconda fase con valori molto elevati dell’accelerazione. Tutte le bolle di vapore che collassano a dimensioni di ordini di grandezza più piccole rispetto alla dimensione massima, inevitabilmente tendono a far emergere, da quel collasso, una nuvola di bolle più piccole. Nel caso di bolle collassanti che sono predominantemente riempite di gas, questo comportamento è meno probabile, poiché i valori minori dell’accelerazione, di questi casi, rendono debole l’instabilità e a volte implicano la stabilità sferica. Quindi, bolle di cavitazione che contengono sostanzialmente gas, spesso rimangono sferiche durante la fase di ridelimitazione. In altri casi invece l’instabilità è sufficiente per causare la frammentazione. Molti esempi di frammentazione e alta distorsione delle bolle, che derivano dalla fase di ricostruzione sono mostrati in Figura 3.3, in cui gli effetti termici sono sostanziali. Figura 3.3: Fotografie di una bolla di etere in glicerina prima (sinistra) e dopo (destra) il collasso e la ridelimitazione. La nuvola sulla destra è il risultato dei successivi cicli di collasso e ridelimitazione. Frost e Sturtevant (1986) Una caratteristica sostanziale nel collasso di molte bolle di vapore è lo sviluppo di un getto rientrante, causato da un’asimmetria come nel caso di pareti solide vicine. Per causa dell’asimmetria una parte della bolla accelera verso l’interno molto più rapidamente rispetto alla parte opposta, e questo si traduce in uno sviluppo di un microgetto rientrante ad alta velocità, che penetra la bolla. 43 Figura 3.4: Formazione iniziale del microgetto rientrante a causa della parete solida trasparente segnata con la linea tratteggiata. Benjamin e Ellis (1966) Figura 3.5: Protuberanza generata dalla gravità dopo che il getti rientrante è progredito attraverso la bolla penetrando il fluido dalla parte opposta. Benjamin e Ellis (1966) La Figura 3.5 è una delle prime, se non la prima, che mostra il risultato di un getto, diretto verso l’alto, prodotto dalla gravità, che è progredito attraverso la bolla e penetrato nel fluido dall’altra parte, creando la protuberanza appuntita. La figura 3.6 invece mostra la comparazione tra l’osservazione sperimentale (Lauterborn-Bolle 1975) e quella studiata (Plesset-Chapman 1971) di un getto rientrante sviluppato in una bolla che collassa vicino ad una parete solida. 44 Figura 3.6: Collasso di una bolla molto vicina alla superficie solida in un liquido quiescente. La forma teorica di Plesset- hapman (lin a ontinua) omparata on l’oss rvazion sp rim ntal i Lauterborn-Bolle (puntini). Plesset e Prosperetti (1977) Figura 3.7: Calcolazione numerica del collasso di cinque bolle, che mostra la formazione dei getti rientranti per le quattro bolle esterne. Chahine e Duraiswami (1992) Un’altra forma di asimmetria che può causare la formazione di un getto rientrante è la vicinanza di altre bolle in una nuvola di bolle. Inoltre, le bolle sul bordo esterno di questa nuvola, tendono a sviluppare getti diretti verso il centro della nuvola, come esemplificato in Figura 3.7. Un’altra manifestazione di instabilità è dovuta al collasso di una bolla vicino alla superficie libera, che produce un getto rientrante, diretto lontano dalla superficie. Quando la bolla è inizialmente sferica, ma vicina alla parete, lo sviluppo tipico dei microgetti è illustrato in Figura 3.8, quando la bolla è più lontana dalla parete gli eventi 45 successivi sono alquanto diversi; come mostra la Figura 3.9, c’è la formazione di due vortici sulle bolla toroidale (frame 11), dopo che i microgetti hanno completato la loro penetrazione nella bolla originale. Figura 3.8: Serie di fotografie che mostrano lo sviluppo dei microgetti rientranti in una bolla collassante molto vi ina alla par t soli a. gni fram pr so opo 2μs ha una largh zza i 1.4mm. Tomita e Shima (1990) Figura 3.9: Foto simili alla precedente, ma con una maggiore distanza della bolla dalla superficie solida. Tomita e Shima (1990) 46 Inoltre, in alcune fotografie, in cui le bolle sono a circa un diametro di distanza dalla parete, si vede che il collasso iniziale è abbastanza sferico e che il getto rientrante penetra il fluido tra la bolla e la parete appena la bolla si ridelimita dopo il primo collasso. Questo stadio è molto simile alla Figura 3.5 ma con la protuberanza diretta al muro. D’altra parte, quando la bolla iniziale è molto vicina alla parete e il collasso comincia con una forma di calotta sferica, le foto mostrano una bolla che si spiattella (a forma di frittella) verso la superficie come si vede dalla Figura 3.10. In tali circostanze è difficile visualizzare il microgetto. Plesset e Chapman (1971) calcolarono numericamente la distorsione di una bolla inizialmente sferica durante il collasso vicino ad un confine solido, e come dimostra la Figura 3.6 i loro profili sono in accordo con le osservazioni sperimentali. Quando una bolla collassa in un fluido fermo, vicino ad una parete, i getti rientranti raggiungono alte velocità abbastanza presto nel processo, e molto prima che il volume del liquido raggiunga dimensioni alle quali, per esempio, la comprimibilità diventi importante. La velocità dei getti rientranti, UJ, al momento dell’impatto con la superficie opposta della bolla è data dalla 1 p 2 ξ( ) [3.8] in cui ξ è costante e p è la differenza tra la pressione remota che manterrebbe la bolla in equilibrio al suo massimo raggio, e la pressione remota presente durante il collasso. Se si ha la fissione della bolla, dovuta alla rottura causata dai microgetti o dagli effetti di instabilità della seconda fase, molte delle osservazioni sperimentali del collasso mostrano che una bolla che arriva dalla prima ricostruzione non è una singola bolla ma una nuvola di bolle più piccole. Sfortunatamente, la velocità degli ultimi momenti del collasso è così elevata che non si ha il tempo necessario per vedere i dettagli di questa fissione. Il comportamento dinamico della nuvola di bolle può essere differente da quello di una singola bolla, ad esempio, lo smorzamento dei cicli di collasso e ridelimitazione è maggiore. Per concludere, bisogna enfatizzare che le osservazioni fatte, sono pertinenti al collasso di una bolla in un fluido fermo. 47 Figura 3.10: Serie di fotografie del collasso di una bolla emisferica su una superficie solida. Benjamin e Ellis (1966) 3.6 Danni della cavitazione Il problema più comune causato dalla cavitazione è il danno ai materiali che le bolle possono provocare quando collassano vicino alla superficie solida. La questione è complicata, sia perché coinvolge fenomeni di flusso instabile sia per le reazioni del materiale di cui è fatta la superficie stessa. Per tale motivo ci sono molte regole empiriche per aiutare a valutare il tenore del danno da cavitazione in una data applicazione. In precedenza si è visto che il collasso di una bolla di cavitazione è un processo violento, che genera onde d’urto di grande ampiezza, e microgetti nel fluido. Quando tali fenomeni avvengono vicino alla superficie solida, generano su di essa degli stress altamente localizzati. Questi carichi, dovuti a ripetuti collassi, causano fenomeni di fatica localizzata sulla superficie, e il successivo distacco o sfaldamento di pezzi di materiale. Nella figura 3.11 è riportata una tipica foto del danno localizzato prodotto dalla cavitazione sulla pala di una pompa a flusso misto, esso ha l’apparenza di una frastagliatura, con rottura a fatica, facilmente distinguibile dall’erosione dovuta alle particelle solide che ha un aspetto molto più regolare e uniforme. Nel caso di acciaio, gli effetti della corrosione spesso intensificano la velocità del danno da cavitazione. 48 Figura 3.11: Danno da cavitazione sulle pale di una pompa a flusso misto. Brennen (1995) Una volta che la nuvola di cavitazione si è generata, le piccole bolle continuano a collassare collettivamente, mantenendo lo stesso comportamento dinamico qualitativo della singola bolla, includendo la possibile formazione di onde d’urto in seguito al punto di minimo volume. Per studiare la formazione della nuvola sono stati usati diversi metodi, tra cui il metodo fotografico Schlieren ad alta velocità, mostrando che l’onda d’urto sferica è realmente generata dalla nuvola residua al momento di minimo volume. La Figura 3.12 mostra una serie di foto di bolle collassanti con la corrispondente traccia di pressione; l’istante di minimo volume risulta tra il sesto e il settimo fotogramma e la traccia mostra chiaramente il picco di pressione che avviene in quell’istante. Con questo metodo fotografico il microgetto sembrerebbe avere un ruolo sussidiario. Un altro metodo è quello di utilizzare un materiale fotoelastico in modo da poter osservare gli stress nel solido e misurare le pulsazioni acustiche simultaneamente. Usando come innesco il primo collasso della bolla si impiega un intervallo di tempo variabile per fare le fotografie dello stato di stress del solido a vari istanti relativi al secondo collasso. Registrando simultaneamente la pressione nel liquido, si può confermare che lo stress impulsivo nel materiale inizia allo stesso momento dell’impulso acustico. Tuttavia in investigazioni più approfondite si osservò che l’impulso di stress era il risultato dei microgetti e anche del collasso della nuvola residua, quest’ultimo da 2 a 3 volte più grande rispetto a quello dei microgetti. Ci sono dei dispositivi che permettono di valutare queste suscettibilità comparative; il più comune consiste di un dispositivo che oscilla un campione, in un liquido, producendo periodici sviluppi e collassi di bolle di cavitazione sulla superficie del campione. I test durano molte ore, con regolari pesature per determinare le perdite di 49 peso, mostrando che la perdita di materiale non è costante, questo perché una superficie irregolare e danneggiata può produrre un’alterazione del modello di cavitazione. In queste prove si utilizza un dispositivo magnetostrittivo al fine di raggiungere frequenze standard dai 5 ai 20 kHz. Queste frequenze causano le più grandi nuvole di bolle di cavitazione sulla superficie del campione, perché sono vicine alla frequenza naturale di una significante frazione di nuclei presenti nel liquido. Nella maggior parte dei dispositivi i danni da cavitazione sono molto indesiderati, però, ci sono circostanze in cui il fenomeno è utilizzato vantaggiosamente, ad esempio nel meccanismo di taglio di rocce con alta velocità dei getti d’acqua (taglio causato almeno in parte dalla cavitazione nel getto appena si ha il contatto con la superficie rugosa della roccia) o nella cavitazione medica per l’eliminazione di adiposità localizzate. Ulteriori esempi di danni della cavitazione sono riportati in Figura 3.20 e 3.21. Figura 3.12: Bolla collassante vicino ad una parete solida in funzione della traccia di pressione caratteristica della superficie. Shima, Takayama, Tomita e Ohsawa (1983) 50 Figura 3.13: Fotografia di una nuvola di cavitazione formata acusticamente. Plesset e Ellis (1955) Figura 3.14: ist assiali ll’ingr sso lla avitazion l danno da cavitazione sul mozzo di una girante di pompa centrifuga. Soyama, Kato, and Oba (1992) 51 3.7 Collasso delle nuvole In molti dispositivi pratici il danno da cavitazione avviene in aree ben localizzate, ad esempio sulla girante di una pompa, come visto nelle figure precedenti. Questo è dovuto da periodici e coerenti collassi di una nuvola di bolle cavitanti. Una tipica nuvola di cavitazione è mostrata in Figura 3.13. In altre macchine idrauliche la periodicità può avvenire naturalmente come risultato del distacco di vortici cavitanti, o per risposta ad un disturbo periodico imposto sul flusso. In quasi tutti questi casi, il collasso coerente della nuvola può causare rumore molto più intenso e maggior potenziale di danno rispetto ad un flusso non fluttuante. Di conseguenza il deterioramento è molto severo sulla superficie del solido vicino alla locazione del collasso della nuvola. Un esempio di questo fenomeno è incluso nella Figura 3.14 in cui le nuvole di cavitazione collassano in una posizione specifica, portando al danno localizzato della foto di destra. 3.8 Rumore della cavitazione Il violento e catastrofico collasso delle bolle di cavitazione produce, oltre alla possibilità di danno materiale, anche un particolare rumore, che è una conseguenza delle momentanee alte pressioni che si vengono a generare quando il contenuto della bolla è altamente compresso. Se si denota il volume della cavità che varia in funzione del tempo con V(t), segue che la variazione della componente di pressione nel tempo, lontano dalla bolla, è data da: 2 pa [3.9] t2 4π dove pa è la pressione acustica irradiata, ed R è la distanza dal centro della cavità dal punto di misura. Poiché il rumore è direttamente proporzionale alla seconda derivata del volume rispetto al tempo, è chiaro che l’impulso del rumore, generato al collasso della bolla, avviene per il grande valore positivo di d2V/dt2, quando la bolla è vicina alla sua dimensione minima. Si usa presentare il livello del suono, usando una scarto quadratico medio di pressione, o pressione acustica, pS, definita come p2 ̅̅̅ p2a ∞ ∫ f f [3.10] 0 con frequenza f e funzione di densità spettrale G(f). Il rumore scoppiettante che accompagna la cavitazione è una delle caratteristiche più evidenti di questo fenomeno per ricercatori e ingegneri. L’insorgenza della cavitazione è spesso individuata prima da questo rumore piuttosto che dall’osservazione delle bolle, e quindi per gli ingegneri è il primo metodo di rilevazione in dispositivi come pompe e valvole. Ci sono infatti, molti metodi empirici per stimare il tasso del danno materiale misurando il rumore generato. 52 Figura 3.15: Spettri di energia del rumore da un modello di valvola che opera sotto condizioni di non avitazion (σ 0.523) i avitazion (σ 0.452 0.342). artin t al. (1981) Il rumore dovuto alla cavitazione tipica nell’orifizio di un controllo idraulico a valvola è presentato in Figura 3.15. La curva più bassa a σ 0.523 rappresenta il rumore turbolento di un flusso non cavitante; sotto il numero tipico di incipiente cavitazione si ha un drammatico incremento nel livello di rumore alla frequenza di 5 o più kHz. Il picco dello spettro tra 5 e 10 kHz corrisponde alla frequenza naturale dei nuclei presenti nel flusso. Figura 3.16: Frequenza dei picchi acustici in funzione del numero tipico di cavitazione. Dati per getti cavitanti da Franklin e McMillan (1984) 53 Nella pratica, molti esperimenti sembrano esibire un comportamento in frequenza che varia come f- -1 o f- -2, mostrando un debole decadimento. Molti autori hanno anche analizzato gli effetti della comprimibilità del liquido. I picchi in molte curve del rumore, da flussi cavitanti (esempio in Figura 3.15), tendono ad abbassare le frequenze al diminuire del numero tipico di cavitazione. Questa tendenza è ulteriormente illustrata dai dati per getti cavitanti di Figura 3.16. Ceccio e Brennen (1991) registrarono il rumore della cavitazione di bolle in una portata, riportando un tipico segnale acustico come in Figura 3.17. Il grande impulso positivo a circa 450μs corrisponde al primo collasso della bolla, questo è seguito da qualche oscillazione dipendente dall’impianto, e da un secondo impulso a circa 1100μs, corrispondente al secondo collasso dopo la ridelimitazione del primo. Figura 3.17: Segnale acustico da una singola bolla collassante. Ceccio e Brennen (1991) Una buona misura della grandezza dell’impulso di collasso è l’impulso acustico I definito come l’area sotto l’impulso t2 ∫ pa t [3.11] t1 con t1 e t2 i tempi prima e dopo l’impulso al quale pa è zero. Si definisce anche un impulso adimensionale I* come * 4π ∞ 2 [3.12] con U∞ e RH velocità di riferimento e lunghezza di riferimento nel flusso. La media degli impulsi acustici per una bolla che collassa su due forme assialsimmetriche si riporta in Figura 3.18, con gli impulsi previsti dalla Rayleigh-Plesset. Poiché queste calcolazioni teoriche assumono che la bolla rimanga sferica, la 54 discrepanza tra teoria e pratica non è così sorprendente, infatti una interpretazione della figura è che la teoria può fornire una stima dell’ordine di grandezza e un limite massimo del rumore prodotto da una singola bolla. In realtà la deviazione dalla sfericità produce un collasso meno focalizzato e quindi meno rumore. Figura 3.18: omparazion ll’impulso a usti o , pro otto al ollasso i una singola olla i avitazion su due forme assialsimmetriche in funzione del massimo volume prima del collasso. Simbolo vuoto per il corpo Schiebe, simbolo pieno per corpi ITTC. Ceccio e Brennen (1991) Nel prossimo passo si considera la sintesi del rumore da cavitazione, dal rumore prodotto da singoli eventi. Se l’impulso prodotto da ogni evento è denotato da I ed il numero di eventi per unità di tempo da NE, il livello di pressione sonora risulta p [3.13] omettendo alcuni fattori di proporzionalità per motivi di chiarezza, quindi i risultati sono intesi solo come guida qualitativa. Sia risultati sperimentali, sia analisi basate sulla Rayleigh-Plesset indicano che l’impulso adimensionale prodotto da un singolo evento di cavitazione è fortemente correlato con il massimo volume delle bolla prima del collasso, ed è anche indipendente dagli altri parametri di flusso. Segue dalla [3.10] e [3.12] che * 1 ∞ 2 *( t 55 * t2 ( t * + t1 [3.14] Se il raggio della bolla al tempo t1 è denotato da RX e il coefficiente di pressione nel liquido con CPX, allora * 2 8π ( σ * 1 2 [3.15] L’integrazione numerica della Rayleigh-Plesset per un intervallo tipico di circostanze rende RX/RM ≈ 0.62 dove RM è il raggio di massimo volume e (CPX - σ) è proporzionale a RM/RH cosicché: * β( 5 2 [3.16] * Si nota una relazione abbastanza simile tra I* e RM/RH dell’onda di suono generata dal collasso di una bolla di gas in un liquido comprimibile. Da sopra segue che 2.5 ∞ 0.5 β 12 [3.17] Quindi la valutazione dell’impulso da un singolo evento è completa con la stima di RM. Se tutti i nuclei fluenti attraverso un certo tubo di flusso noto, di area circolare AN, che cavitano similarmente, allora [3.18] ∞ con N concentrazione di nuclei nel flusso entrante, segue che il livello di pressione sonora risultante dalla sostituzione di [3.19], [3.18] nella [3.14], ha la forma: p β 3 2 ∞ 2 σ min 2.5 [3.19] Leggi diverse si applicano quando la cavitazione è generata da fluttuazioni turbolente, e le tipiche tensioni dei nuclei che si muovono lungo un percorso in una portata turbolenta sono più difficili da stimare. Inoltre, quando la quantità delle bolle è elevata, la radiazione del rumore inizia ad essere influenzata dall’interazione tra le bolle. 3.9 Luminescenza Le altissime pressioni e temperature che possono verificarsi in un gas incondensabile durante il collasso, sono da ritenersi responsabili per il fenomeno noto come luminescenza, cioè l’emissione di luce osservata durante il collasso della bolla di cavitazione. Il fatto che la luce fosse emessa al collasso fu dimostrato da Meyer e Kuttruff (1959), che osservarono la cavitazione su di un’asta oscillante magnetostrittiva e correlarono la luce al punto di collasso nel ciclo di cavitazione. 56 Come già visto nel paragrafo 3.4, le temperature possono essere stimate sulla base della compressione uniforme dei gas non condensabili, e raggiunte per frazioni di microsecondi. Queste condizioni spiegherebbero l’emissione della luce. Alcuni esperimenti individuarono temperature più alte, rispetto a quelle trovate prima, con durate dei picchi dell’ordine di picosecondi. Alcune osservazioni su questi esperimenti si concentrano sul fatto che le bolle collassanti formano una propagazione sferica di un urto nel gas contenuto dalla bolla, e che la concentrazione dell’urto al centro della bolla è un’importante ragione per l’alta temperatura associata alla radiazione sonoluminescente. Quando accade nel contesto di cavitazione acustica la luminescenza è chiamata sonoluminescenza, per evidenziare il fatto che sia la cavitazione piuttosto che il suono a causare l’emissione di luce. Come ci si aspetterebbe dalla Rayleigh-Plesset, la tensione superficiale e la tensione di vapore del liquido sono importanti per la determinazione della sonoluminescenza del flusso, come è chiaramente mostrato in Figura 3.19 La sonoluminescenza è anche fortemente dipendente dalla conduttività termica del gas, e questo è molto evidente con gas il cui coefficiente di conduzione termica è basso. Chiaramente la conduzione del gas gioca un ruolo fondamentale nel fenomeno, e la rottura della bolla per collasso completo potrebbe portare all’eliminazione dello stesso. L’emissione di luce in un flusso cavitante fu osservata per prima in un tubo venturi, identificando la fonte come la regione del collasso della bolla, trovando anche che un impulso di pressione acustica è associato ad ogni lampo di luce. Si studiarono soprattutto i comportamenti dell’acqua a questo fenomeno, con diversi tipi di gas incondensabili disciolti. Ci furono però altri esperimenti, che trovarono molte difficoltà nell’osservare luminescenze in portate cavitanti, il perché fu associato al fatto che solo il collasso di bolle con una simmetria significativamente sferica producessero luminescenze, e questo accade raramente in molti flussi. Figura 3.19: Correlazione del flusso sonoluminescente per diversi tipi di liquidi. Jarman (1959) 57 Figura 3.20: Elica navale con danni provocati dalla cavitazione in vicinanza del bordo. www.wikipedia.org Figura 3.21: Girante di pompa centrifuga danneggiata dalla cavitazione con rottura a fatica. www.wikipedia.org 58 CAPITOLO 4 Dinamica delle bolle oscillanti 4.1 Introduzione In questo capitolo ci si focalizza sulla risposta di una bolla a continue oscillazioni di pressione. Una utile classificazione usa la grandezza del raggio oscillante della bolla in risposta al campo di pressione fluttuante imposto. Si possono identificare tre regimi: 1. Per piccole ampiezze di oscillazione della pressione la risposta è lineare 2. A causa della non linearità nell’equazioni governanti, particolarmente nella Rayleigh-Plesset, la risposta di una bolla comincia ad essere influenzata da queste non linearità all’aumento dell’ampiezza di oscillazione. Tuttavia, la bolla può continuare ad oscillare stabilmente. Sotto queste circostanze ci si riferisce alla “cavitazione acustica stabile”, molti fenomeni non lineari possono influenzarla in modo importante, alcuni di questi sono: la produzione di subarmoniche, il fenomeno di diffusione rettificata e la generazione delle forze di Bjerknes. 3. Sotto altre circostanze il cambio della dimensione della bolla durante un singolo ciclo di oscillazione può diventare così grande che la bolla subisce un ciclo di crescita di cavitazione esplosivo e un collasso violento, simile a quello descritto in precedenza. Questo tipo di risposta è detta “cavitazione acustica transitoria”, ed è distinta dalla stabile per il fatto che il raggio della bolla cambia di molti ordini di grandezza durante il ciclo. È necessario notare che altri fattori sono importanti nella determinazione della risposta per un dato campo di oscillazione della pressione. Uno di questi è la relazione tra la frequenza, ω, dell’oscillazione imposta e la frequenza naturale, ωN, della bolla. Qualche volta è caratterizzato dalla relazione tra il raggio di equilibrio della bolla, RE, in assenza di oscillazioni di pressione e la dimensione della bolla ipotetica, RR, che risuonerebbe alla frequenza imposta, ω. Un altro fattore è la relazione tra l’ampiezza di oscillazione della pressione, p̃ , e la pressione media p∞. Ad esempio se p̃ < p∞ la bolla non risulta mai sotto tensione e non caviterà. Anche la predominante del riempimento della bolla è un fattore importante, se con vapore si ha “cavitazione acustica transitoria”, se con gas si ha “cavitazione acustica stabile”. 4.2 Frequenza naturale Inizialmente si sono trascurati gli effetti termici e l’influenza della comprimibilità del liquido. Come si vedrà più avanti, entrambi aumentano diventando più importanti dei termini viscosi che si erano considerati. 59 Tenendo presente la soluzione linearizzata della 2.9 quando la pressione all’infinito ̅ ∞ , su cui è stata sovrapposta una piccola oscillazione coincide col valore della media p di pressione di ampiezza, p̃ , e una pulsazione ω cosicché: p∞ ̅∞ p {p̃ jωt } [4.1] jωt [4.2] La risposta dinamica lineare della bolla è allora [1 {ϕ }] dove RE è la dimensione di equilibrio alla pressione media, e la risposta del raggio della bolla ϕ è in generale un numero complesso in modo che RE|ϕ| sia l’ampiezza dell’oscillazione del raggio della bolla. La fase, ϕ, rappresenta la sfasatura tra p∞ e R. Si assume che la massa di gas nella bolla rimanga costante (mG). Sostituendo la [4.1] e [4.2] nella [2.9], trascurando tutti i termini di ordine ϕ2 e usando la condizione di equilibrio pV - p∞ + pGE - 2S/RE = 0, si ottiene: ω2 jω 4ν 1 2 2 [ 2 3 p p̃ ] 2 ϕ [4.3] dove, come prima: p ̅∞ p 2 p 3m [4.4] 3 4π Segue che per una data ampiezza, il picco di risposta d’ampiezza avviene alla frequenza, ωP, dato dal minimo valore del raggio della parte sinistra nella [4.3] ω [ 3 p 8ν2 2 / 2 4 1 2 [4.5] ] oppure in termini di differenza di pressione: ω [ ̅∞ 3 (p p ) 2 3 2 8ν2 1 3 4 1 2 ] [4.6] A questo picco di frequenza l’ampiezza di risposta è inversamente proporzionale allo smorzamento: |ϕ|ω p̃ ω 4μ *ω2 4ν2 1/2 4+ [4.7] È utile definire la frequenza naturale dell’oscillazione della bolla come valore del picco di frequenza per uno smorzamento nullo: 60 ω { 1 2 1 2 ̅∞ [3 (p p ) 2 3 1 [4.8] ]} La connessione con il criterio di stabilità di paragrafo 2.5, è chiaro quando si osserva che non esistono frequenze naturali per tensioni maggiori di 4S/3RE, oscillazioni stabili possono esistere solo per equilibri stabili. La frequenza di picco è una quantità importante da considerare in ogni problema di dinamica della bolla. Figura 4.1: Frequenza di risonanza di una bolla in acqua a 300K, ( 0.0717, μL 0.000863 come funzione del raggio della bolla per varie pressioni. Brennen (1995) L = 996.3) Si nota dalla [4.6] che la frequenza di picco, ωP, è funzione solo della tensione, di RE, e dalle proprietà del liquido, inoltre si vede che il secondo e terzo termine dominano a raggi di equilibrio molto piccoli, e la frequenza è quasi indipendente dalla tensione, (vedi Figura 4.1 e 4.2). Per bolle grandi il termine viscoso diventa trascurabile e il picco di frequenza dipende dalla tensione. È importante prendere nota che per i nuclei comunemente trovati in acqua, nell’intervallo di grandezza da 1 a 100μm, le frequenze naturali sono dell’ordine dai 5 ai 25kHz. Inoltre si nota che, come ogni oscillatore, un nucleo eccitato alla sua frequenza di risonanza, ωP, esibisce una risposta la cui ampiezza è principalmente funzione dello smorzamento. Poiché lo smorzamento viscoso è piuttosto piccolo in molte circostanze pratiche, l’ampiezza data dalla 4.7 può essere molto grande a causa del fattore μL al 61 Figura 4.2: Frequenza di risonanza della bolla nel sodio a 800K in funzione del raggio della bolla per varie pressioni. Brennen (1995) denominatore. Questo potrebbe causare che il nucleo ecceda la sua dimensione critica, RC, e che risulti un comportamento altamente non lineare con ampiezze molto grandi. L’ampiezza di pressione critica risulta: p̃ 4μ *ω2 4ν2 4 1 2 + {[1 ω2 2 1 3 2 ] 1} [4.9] e in molte circostanze questa è approssimativamente uguale a 4μLωN. Per esempio: 10μm nell’acqua a 300K frequenza naturale circa 10kHz con l’ampiezza di pressione critica circa 0.002bar che è molto bassa. 4.3 Costante politropica effettiva Si valuta l’assunzione che il gas nella bolla si comporti politropicamente in accordo con pG = pG0(R0/R)3K. L’assunzione politropica è di solito accettabile per le seguenti ragioni: Durante la crescita di una bolla di vapore, il gas gioca un ruolo relativamente piccolo, e la preponderanza di vapore tende a determinare la temperatura della bolla. Durante la fase di collasso, quando il gas predomina, le velocità sono così alte che sembra adeguata un’assunzione adiabatica γ. 62 Poiché una bolla collassante perde la sua simmetria sferica, il movimento interno del gas genera mescolamento, che tende a negare ogni altro modello sofisticato basato sulla simmetria sferica. Il problema dell’appropriata costante politropica è direttamente accoppiato con la valutazione dello smorzamento termico effettivo della bolla, che fu analizzato in dettaglio da Prosperetti (1977), con particolare attenzione alla diffusione termica nel gas, e riuscendo a predire l’effettivo esponente politropico mostrato in Figura 4.3, dove è plottato in funzione della frequenza ridotta, ω E2/αG, per diversi valori di diffusività termica adimensionale nel gas, definita come α* α ω 2 [4.10] con αG e cG rispettivamente la diffusività termica e la velocità di propagazione del suono nel gas. Si nota che per basse frequenze il comportamento tende a diventare isotermico con K = 1, d’altra parte, con frequenze alte, il comportamento tende a diventare isoentropico. Inoltre, ad alte frequenze, l’esponente può assumere valori fuori dall’intervallo 1 < < γ. La Figura 4.4 mostra i dati di esperimenti tipici di bolle d’aria in acqua. Risultati che sono coerenti con la teoria per frequenze minori a quella di risonanza. Figura 4.3: Esponente politropico effettivo, K, in funzione della frequenza di riduzione per una determinata diffusività termica del gas. Prosperetti (1977) 63 Riassumendo 1 3 {δ} [4.11] dove {ζ} è una funzione complessa. Finchè l’uso di un’esponente effettivo politropico provvede un consistente avvicinamento per le oscillazioni lineari, si è dimostrato che può causare significativi errori quando l’oscillazione diventa non lineare. Sotto queste circostanze il comportamento del gas può divergere da quello che è il comportamento con un esponente politropico effettivo. 4.4 Termini di smorzamento Chapman e Plesset (1971) hanno presentato un utile riassunto dei tre principali contributi allo smorzamento delle oscillazioni della bolla, cioè a causa della viscosità del liquido, della comprimibilità del liquido attraverso la radiazione acustica e della conduttività termica. È particolarmente conveniente rappresentare i tre componenti di smorzamento come contributi additivi di una viscosità effettiva del liquido μE, che si può sostituire nella Rayleigh-Plesset al posto della viscosità attuale μL: μ μ μ μ [4.12] dove la viscosità “acustica” μ ω2 4 2 [4.13] con cL velocità del suono nel liquido. Figura 4.4: isur sp rim ntali ll’ spon nt politropi o, , p r oll ’aria in a qua in funzion rapporto tra raggio di equilibrio e raggio di risonanza della bolla. La linea solida rappresenta i risultati teorici. Crum (1983) 64 l La viscosità “termica”, μT, segue dalle stesse analisi usate per ottenere l’esponente politropico effettivo μ ̅∞ p 2 4ω m{δ} [4.14] La grandezza relativa dei tre componenti dello smorzamento può essere piuttosto differente per diverse grandezze di RE. In Figura 4.5 si illustrano i dati per bolle d’aria in acqua a 20°C, dove si nota che la componente viscosa è importante per bolle molto piccole, la componente acustica è predominante per bolle più grandi di un centimetro, mentre la componente termica è dominante per molte bolle di dimensioni considerevoli. Figura 4.5: Termini di smorzamento della bolla e smorzamento totale in funzione del raggio di equilibrio, per acqua. Chapman e Plesset (1971) 4.5 Effetti non lineari Nel precedente paragrafo si è assunto che la perturbazione nel raggio della bolla, ϕ, fosse sufficientemente piccolo cosicché si possa mantenere l’approssimazione lineare, però, le bolle singole, esibiscono un numero di interessanti e importanti fenomeni di non linearità. Quando un liquido che contiene inevitabilmente microbolle, viene irradiato da un suono ad una data frequenza, ω, la risposta non lineare risulta come una dispersione armonica, che non produce solo armoniche con frequenze multiple di ω (superarmoniche) ma più insolitamente subarmoniche con frequenza minore di ω della forma mω/n con m, n interi. Sia le super che le subarmoniche diventano importanti quando l’ampiezza di eccitazione aumenta. 65 Ci furono molte investigazioni analitiche e numeriche per la produzione delle subarmoniche, utilizzando soluzioni della Rayleigh-Plesset per esplorare le caratteristiche non lineari di una singola bolla eccitata da una pressione oscillante con una frequenza ω. Come ci si poteva aspettare, diversi tipi di risposta arrivano, a seconda che ω sia più grande o più piccola della frequenza naturale della bolla ωN. Figura 4.6: Risposta per frequenza minore alla natural in alto on ω/ωN = 0.8, e maggiore alla naturale in asso on ω/ωN 1.8, s mpr a un’ampi zza i 0.33 ar. Flynn (1964) La Figura 4.6 mostra due esempi del tipo di risposta, una per ω < ωN e l’altra per ω > ωN. Notare la presenza di subarmoniche in entrambi i casi. Si esaminarono soluzioni numeriche per un gran numero di frequenze di eccitazione diverse e si costruirono curve di risposta a tali frequenze come in Figura 4.7, in cui si nota il progressivo sviluppo delle risposte di picco a frequenze subarmoniche all’aumentare dell’ampiezza di eccitazione. Gli effetti non lineari non solo creano i picchi sub-armonici, ma anche picchi di risonanza distorti (a sinistra), creando delle discontinuità indicate con la linea verticale tratteggiata in figura. Queste corrispondono a biforcazioni o transizioni improvvise tra due soluzioni valide, una con un’ampiezza più grande dell’altra. 66 Figura 4.7: mpi zza ll os illazioni ra iali i una olla i raggio 1μm in a qua a 1 ar i pr ssione in funzion l rapporto ω/ωN. Lauterborn (1976) 4.6 Analisi debolmente non lineari Finchè le ampiezze di oscillazione sono ancora abbastanza piccole è valido usare una tecnica espansiva per investigare gli effetti delle piccole non linearità. Si devono mantenere solo i termini che sono quadratici nell’ampiezza di oscillazione, i termini cubici, o di ordine maggiore, sono trascurati. Il risultato è la seguente forma non lineare della [4.3]: n 1 p̃ n ω2 2 β0 n ϕn n ∑ β1 n, m ϕm ϕn m 1 m ∑ β2 n,m ϕ̂ m ϕn m [4.15] m 1 dove β0(n), β1(n,m) e β2(n,m) sono funzioni di ν/ωNRE2, / LωN2RE3, K e δ (discretizzazione del dominio di frequenza). Date le caratteristiche del fluido e della bolla la [4.15] può essere risolta iterativamente per trovare ϕn. Il valore di N dovrebbe essere abbastanza grande da comprendere tutte le armoniche con ampiezza significante. Esaminando prima le caratteristiche delle oscillazioni radiali che sono date da una singola frequenza di eccitazione, è chiaro che, dalla 4.15 , l’unica ϕn non nulla avviene a frequenze che sono multipli interi della frequenza di eccitazione. Di conseguenza per questa classe di problemi si sceglierebbe δ per rimanere alla frequenza di eccitazione. La Figura 4.8 mette a confronto la soluzione della non linearità e l’integrazione numerica della Rayleigh-Plesset che chiaramente divergono all’aumentare dell’ampiezza di oscillazione; comunque la figura mostra che le soluzioni debolmente non lineari sono qualitativamente buone. 67 Figura 4.8: Confronto tra soluzioni debolmente non lineari e integrazione numerica della Rayleigh-Plesset. Brennen (1995) La Figura 4.9 mostra esempi dei valori di |ϕn| per tre diverse ampiezze di eccitazione e dimostra come le armoniche diventano più importanti all’aumentare dell’ampiezza. Le soluzioni debolmente non lineari possono anche essere usate per costruire spettri di risposta in frequenza, come si vede in Figura 4.10, dove è tracciata la massima deviazione possibile del raggio di equilibrio (Rmax - RE)/RE, in funzione della frequenza di eccitazione dove per convenienza si stima: max ∑ ϕn [4.16] n 1 Chiaramente le soluzioni debolmente non lineari esibiscono risonanze subarmoniche simili a quelle viste in molti casi più comuni. I vantaggi delle analisi debolmente non lineari diventano più visibili quando ci sono problemi di geometria più complessa o frequenze multiple di eccitazione. 68 Figura 4.9: s mpi i gran zza ll’ampi zza ll armoni h ϕn|. Brennen (1995) Figura 4.10: Spettri di risposta in frequenza a tre diverse ampiezze di eccitazione: 0.1 a puntini, 0.2 tratteggiata e 0.3 solida. Brennen (1995) 4.7 Soglia per la cavitazione transitoria Ci si sposta ora sulle circostanze in cui avviene la cavitazione acustica transitoria. Il primo paso nello stabilire tale criterio è compiuto dall’analisi di stabilità statica di paragrafo 2.5, dove si era esplorata la stabilità di una bolla quando la pressione lontana dalla bolla era variata, ed identificando una dimensione critica del raggio [2.18], 69 ed una soglia di pressione critica dalla [2.19] che se raggiunta guiderebbe ad una crescita instabile e quindi ad una cavitazione transitoria. La complicazione che si aggiunge in questo contesto, è che ci sia solo un tempo finito durante ogni ciclo per il quale la crescita può avvenire, quindi bisogna capire se c’è o non c’è sufficiente tempo per una significante crescita instabile. Il problema è determinato dalla relazione tra la pulsazione ω dell’oscillazione imposta e la frequenza naturale, ωN della bolla. Se ω ≪ ωN, allora l’inerzia ha poca importanza nelle dinamiche della bolla, e la bolla risponde quasi staticamente. L’ampiezza critica alla quale la cavitazione transitoria avviene risulta dalla 2.18 p̃ ̅∞ p p 4 8π [ 3 9m ] 1 2 [4.17] Dall’altra parte se ω ωN il problema coinvolge le dinamiche della crescita della bolla poiché l’inerzia determina la dimensione delle perturbazioni. Si nota che una condizione necessaria ma non sufficiente per la cavitazione transitoria è che la pressione dell’ambiente cada sotto alla tensione di vapore del liquido, per parte del ciclo di oscillazione. Se assumendo che la soglia venga superata, il tasso di crescita della bolla è dato approssimativamente dal tasso di crescita asintotica della [2.11]; combinandola con il tempo disponibile per la crescita il raggio massimo della bolla è dato da: π p̃ f β * ω ̅∞ p 1 2 + [4.18] dove f(β), data da Apfel (1981), tiene conto di alcuni dettagli come la funzione di mezzo periodo per cui la pressione è negativa. L’ultimo passaggio nella costruzione del criterio per ω ωN è di argomentare che la cavitazione transitoria avviene quando RM → 2 E, e quindi la pressione critica diventa p̃ ̅∞ p 4 2 π2 f ω2 2 [4.19] 4.8 Rettificazione della diffusione di massa Si passa l’attenzione all’effetto non lineare che coinvolge il trasferimento di gas disciolti tra liquido e bolla. Questo importante effetto avviene in presenza di un campo acustico ed è conosciuto come “rettificazione della diffusione di massa”. Si considera una bolla di gas disciolto e si aggiunge una oscillazione alla pressione ambiente. Il gas tenderà ad uscire dalla soluzione della bolla durante quella parte del ciclo di oscillazione quando la bolla è più grande della media, perché la pressione parziale del gas nella bolla si annulla. Il gas si ridiscioglie durante l’altra metà del ciclo, quando la bolla è più piccola della media. 70 Tuttavia ci sono due effetti che tendono ad aumentare la massa di gas nella bolla: Il primo di questi è dovuto al fatto che il rilascio di gas dal liquido avviene durante quella parte di ciclo quando l’area di superficie è grande, e quindi l’afflusso è leggermente più grande dell’efflusso durante l’altra parte del ciclo. Il secondo è dovuto alla diffusione dello strato limite, che tende ad essere allungato e più sottile quando la bolla è più grande e questo migliora il flusso nella bolla. Tale effetto contribuisce con un secondo termine quadratico al flusso netto di gas nella bolla. Analisi recenti che includono tutti i termini di contributi non lineari, portano al seguente risultato della diffusione di massa t * ∞ ( 3 1 ̅ ∞ ))/ (1 2 / p ̅∞ ) 4 /(3 p 2 +* 1 2 π t 1/2 1/2 + [4.20] Il segno del tasso di crescita della bolla è determinato dal segno del termine ∞ 2 (1 ̅∞ p )( 3 2 * [4.21] In assenza di oscillazioni e tensioni superficiali la bolla cresce quando c∞ > cS e si dissolve quando l’opposto è vero. Tuttavia, con presenza di oscillazioni, il termine ( 3/ 2) diminuisce sotto l’unità come l’ampiezza aumenta, questo causa un incremento positivo nel tasso di crescita. Anche in un liquido subsaturo, per cui c∞ < cS, questo incremento potrebbe causare il cambiamento del segno e divenire positivo. Se si applica una pressione oscillante ad un fluido con microbolle di raggio RE, allora la [4.21] dimostra che esiste una certa soglia di ampiezza sopra la quale le microbolle cominciano a crescere per diffusione rettificata. Figura 4.11: s mpi lla r s ita i oll ’aria in a qua satura a ausa lla iffusion r ttifi ata. quattro diverse ampiezze di pressione. La linea solida indica le previsioni teoriche. Crum (1980) 71 r Figura 4.12: ati ll’ampi zza lla pr ssion i soglia p r oll in a qua istillata satura ’aria. a frequenza del suono è di 22.1kHz. Crum (1984) Questa soglia di ampiezza è sufficientemente grande da far si che il rapporto ( 3/ 2) sia così piccolo da far svanire la [4.21]. Misure sperimentali del tasso di crescita e della soglia di ampiezza di pressione sono mostrate in Figura 4.11 e 4.12. Infine, si nota ancora che molte delle teorie assumono la simmetria sferica e che la divergenza dalla sfericità potrebbe alterare la diffusione dello strato limite in modo tale che si potrebbe influenzare radicalmente il processo di trasferimento di massa. La crescita delle bolle causata è molto lenta a confronto con gli altri processi di crescita, tuttavia, questo è un metodo attraverso cui, nuclei molto piccoli e stabili possono crescere sufficientemente per diventare nuclei di cavitazione. Come osservò Blake (1949), la soglia di pressione aumenta con l’aumentare di mG, pertanto, appena la massa di gas aumenta è necessaria una minore riduzione di pressione per creare una bolla instabile. 4.9 Forze di Bjerknes Un effetto non lineare diverso, è la forza competente di una bolla in un campo acustico causato da onde sonore. Il numero d’onda spaziale è denotato da ω/ L, la presenza di queste onde implica un gradiente di pressione istantaneo nel liquido. Sostituendo p̃ p̃ * sin 72 xi [4.22] * nella [4.1], dove la costante p̃ è l’ampiezza delle onde sonore, e xi è la direzione della propagazione d’onda. Come ogni altro gradiente di pressione, questo produce una forza istantanea nella direzione di propagazione dell’onda data da ̅i F π 2 (p̃ * ) sin 2 xi (ω2 2 [4.23] ω ) Questa è conosciuta come forza primaria di Bjerknes, tale formula produce alcuni modelli interessanti di migrazione della bolla in un campo sonoro stazionario. Se ω > ωN la forza primaria di Bjerknes causa migrazione delle bolle in un verso, se ω < ωN le bolle tendono a migrare nel verso opposto. È importante infine menzionare un altro effetto non lineare, un campo acustico può causare movimenti medi nel tempo detti “acoustic streaming”. Il termine microstreaming è usato per riferirsi a questi movimenti vicini a piccole bolle. Generalmente questi movimenti prendono la forma di modelli di circolazione che potrebbero alterare i processi di diffusione di calore o di massa, e quindi modificare fenomeni come la diffusione rettificata. 73 CAPITOLO 5 Flussi cavitanti 5.1 Introduzione Si comincia questa discussione di cavitazione in flussi descrivendo l’effetto del flusso su un singolo evento di cavitazione. Questo è il termine che viene usato nel processo che avviene quando un singolo nucleo di cavitazione è portato in una regione di bassa pressione nel flusso, cresce esplosivamente fino a dimensioni macroscopiche, e collassa quando si trova in una regione con pressione più elevata. Tutti i modelli di analisi sono basati su due assunzioni, che le bolle rimangano sferiche e che i singoli eventi non interagiscano gli uni con gli altri. Le osservazioni su flussi reali mostrano che le bolle di cavitazione sono lontane dalla sfericità, e di solito non ci sono singole bolle, ma piuttosto una nuvola di bolle più piccole. La deviazione dalla sfericità è spesso il risultato dell’interazione della bolla con il gradiente di pressione e le forze nel flusso, o l’interazione con una superficie solida. Quando la frequenza degli eventi di cavitazione aumenta nello spazio o nel tempo, essi cominciano ad interagire gli uni con gli altri, e si può manifestare un fenomeno nuovo. Gli eventi di cavitazione interagendo idrodinamicamente possono portare a conseguenze gravi come quelle viste nel capitolo 3. Aumentando la densità degli eventi si possono manifestare un numero di diversi tipi di cavitazione in larga scala. 5.2 Bolle in movimento Plesset (1948) dimostrò alcune valide approssimazioni per gli eventi di cavitazione usando la relazione [2.6]. Parkin (1952) aggiunse un modello più dettagliato per la crescita di bolle di cavitazione in movimento nel flusso attorno al corpo. Questo assumeva che le bolle iniziassero come micronuclei nel flusso del liquido in arrivo, e che la bolla si muovesse con la velocità del liquido lungo una linea di flusso vicino alla superficie solida. L’incipienza della cavitazione si riteneva avvenisse quando le bolle avessero raggiunto una dimensione osservabile dell’ordine di 1mm. Parkin, credeva che la mancanza di accordo tra questa teoria e le osservazioni sperimentali, fosse a causa del fatto che si era trascurato lo strato limite. Esperimenti seguenti (Kermeen, McGraw e Parkin 1955) rivelarono che la cavitazione poteva risultare da flussi di nuclei liberi all’interno del liquido, o da nuclei originati dalle imperfezioni della superficie di contatto. La maggior parte delle osservazioni sperimentali indicano che la grande maggioranza dei nuclei è già presente nel corpo liquido. Johnson e Hsieh (1966) mostrarono che, poiché le linee di flusso che incontrano la regione di bassa pressione sono vicine alla superficie libera, e quindi vicine alle linee di 74 flusso di ristagno, i nuclei subiscono grandi accelerazioni fluide e gradienti di pressione appena passano vicino a queste zone di ristagno. Quindi l’ideale sarebbe forzare i nuclei ad uscire da queste aree. Johnson e Hsieh chiamarono questo fenomeno come “effetto screening”. Questo è solo uno degli effetti che le accelerazioni e i gradienti di pressione nel flusso possono avere sui nuclei e sulla crescita e collasso di bolle cavitanti. 5.3 Interazioni bolla/flusso Il teorema del massimo modulo afferma che il massimo di una funzione armonica deve essere sul contorno e non all’interno della regione di soluzione di quella funzione. Segue che, la pressione minima di un flusso potenziale, stabile ed aviscoso, deve giacere sul contorno di quel flusso. Inoltre, gli effetti dei fluidi reali in molte portate non alterano il fatto che la pressione minima avvenga sulla, o vicino alla, superficie solida. L’eccezione più comune di questa regola si trova nella cavitazione di vortici, dove gli effetti stabili e/o viscosi associati al distacco dei vortici o alle turbolenze, causano la divergenza dal teorema del massimo modulo. In molti flussi in cui la pressione minima avviene sul confine, si ha che le bolle di cavitazione che si formano nelle vicinanze di quel punto, sono probabilmente influenzate da quel confine (superficie solida) interagendo con esso. Si osserva inoltre che ogni curvatura della superficie solida, o meglio, delle linee di flusso nelle vicinanze del punto di minima pressione, causano gradienti di pressione normali alla superfice. Questi gradienti di pressione forzano la bolla verso la superficie, e possono causare sostanziali deviazioni dalla sfericità. Prima ancora che gli effetti dello strato limite siano presi in considerazione nel quadro generale, è evidente che le dinamiche delle bolle di cavitazione possono essere alterate significativamente da interazioni con la superficie solida vicina, e dal flusso vicino a tale superficie. Prima di descrivere alcune osservazioni sperimentali dell’interazione bolla/flusso, si considerano le dimensioni relative delle bolle di cavitazione e dello strato limite viscoso. Nella portata di un flusso uniforme di velocità, U, attorno ad un oggetto come un’aletta idrodinamica, con dimensione tipica, l, lo spessore dello strato limite laminare vicino al punto di minima pressione è dato qualitativamente da δ (νLl/U)1/2. Tra parentesi, si nota che la transizione di turbolenza, di solito avviene a valle del punto di minima pressione e di conseguenza l’appropriato spessore dello strato limite per limitare la cavitazione, confinandola nelle immediate vicinanze della regione di bassa pressione, è lo spessore dello strato limite laminare. L’analisi approssimata, del paragrafo 2.5, rende un raggio massimo di bolla dato da 2l Segue che il rapporto δ/ M σ [5.1] min è approssimativamente dato da δ 1 1 2 σ min 75 ν 2 { } l [5.2] Quindi, poiché (- σ - CPmin) è dell’ordine di 0.1 o maggiore, si ha che per alti numeri di Reynolds, Ul/νL, che sono tipici in molti flussi in cui la cavitazione è un problema, lo strato limite è di solito più sottile della dimensione tipica delle bolle. Si può anticipare che quelle parti della bolla di cavitazione più lontane dalla superficie solida, interagiscono con il flusso principalmente fuori dallo strato limite, mentre quelle più vicine alla superficie solida sono influenzate dallo strato limite. 5.4 Osservazioni sperimentali Nelle osservazioni iniziali si notò che quasi tutte le bolle di cavitazione erano vicine alla forma semisferica piuttosto che sferica, e che apparivano essere separate dalla superficie solida da un piccolo strato di film di liquido. Ricerche più recenti si focalizzarono su queste interazioni bolla/flusso, estendendole al caso di bolle di cavitazione nel flusso attorno a corpi assialsimmetrici. Si usarono due corpi di tale forma nelle investigazioni, un corpo Schiebe (Schiebe 1972) e un corpo ITTC, una caratteristica importante della forma Schiebe è che lo strato limite non separa la zona di bassa pressione in cui si formano le bolle di cavitazione, l’altro invece esibisce una separazione laminare in quella regione. Per entrambe le forme le isobare nelle vicinanze del punto di minima pressione mostrano un largo gradiente di pressione normale alla superficie, come si vede in Figura 5.1. Questo gradiente è associato alla curvatura del corpo, e quindi alle linee di flusso nelle vicinanze del punto di minima pressione. Segue, a un dato numero di cavitazione, σ, che la regione a pressione sotto alla tensione di vapore, che è racchiusa tra la superficie solida e la superficie isobara, è lunga e fina comparata con la dimensione della forma. Solo i nuclei che passano attraverso questo piccolo volume cavitano. Figura 5.1: Isobare nella vicinanza del punto di minima pressione sulla forma assialsimmetrica Schiebe, con i valori indicati per i coefficienti di pressione. Schiebe (1972) e Kuhn de Chizelle et al. (1992) 76 Si descrivono per prime le osservazioni condotte da Ceccio e Brennen (1991). Foto tipiche di bolle sulla forma Schiebe di diametro di 5.08cm durante il ciclo di crescita e collasso delle bolle sono mostrate in Figura 5.2. Figura 5.2: Serie di fotografie mostranti la crescita e il collasso di bolle di cavitazione in movimento in un flusso attorno ad un orpo hi i iam tro i 5.08 m a σ 0.45 v lo ità i 9m/s. l flusso va a destra a sinistra. Ceccio e Brennen (1991) 77 La bolla, durante la fase di crescita, ha la forma di calotta sferica, ed è separata dalla parete da un sottile strato di liquido, dello stesso ordine di grandezza dello spessore dello strato limite. Gli sviluppi successivi dipendono dalla geometria del corpo e dal numero di Reynolds. Con riferimento alla Figura 5.2, si nota che come la bolla comincia ad entrare nella regione in cui il gradiente di pressione è avverso, la superficie esterna comincia ad essere spinta verso l’interno causando il profilo a cuneo della bolla. Il collasso inizia sulla superficie esterna della bolla, e questo spesso conduce ad uno schiacciamento in avanti, come si vede in figura. Altri due processi avvengono nello stesso momento: Lo spessore della bolla in direzione parallela al flusso diminuisce più velocemente della sua “apertura alare”, cioè dello spessore della bolla in direzione parallela alla superficie della forma e normale al flusso. Quindi la dimensione maggiore della bolla risulta l’apertura alare. (processo di fissione) La bolla acquisisce significante vorticità nella direzione dell’apertura alare, attraverso le interazioni con lo strato limite durante la fase di crescita. Di conseguenza, come il collasso procede, questa vorticità si concentra, e la bolla evolve in uno o più vortici cavitanti. Altri due fenomeni addizionali furono osservati nelle forme ITTC: Si osservò che lo strato di liquido sotto alla bolla diventa perturbato da qualche instabilità. Come si vede in Figura 5.3 questo si traduce in uno strato pieno di bolle che viene rilasciato dietro alla bolla principale. Quindi l’instabilità dello strato liquido conduce ad un altro processo di fissione della bolla. Per via della separazione fisica lo strato pieno di bolle collassa dopo della bolla principale. Meno frequentemente, quando la bolla passa la zona di separazione laminare, innesca la formazione di una “cavitazione allegata” all’estremità laterali in direzione dell’apertura alare della bolla, come in Figura 5.4. Appena la bolla principale procede verso valle, le “striature” o “code” di cavitazione allegata, escono dietro alla bolla principale, poiché le estremità delle code sono attaccate alla superficie solida. Successivamente la bolla principale collassa, lasciando le code che persistono per una lunga frazione di tempo come si nota dalla foto in basso di Figura 5.4. Un’osservazione importante è la presenza di una “fossetta” sulla superficie esterna di tutte le bolle in movimento, come in Figura 5.5. Queste non sono le precorritrici dei getti rientranti, poiché sembrano rimanere stabili durante la maggior parte del processo di collasso. Un’altra nota importante è l’osservazione che ad alti numeri di Reynolds, le code allegate avvengono anche sui corpi di Schiebe, che non esibiscono normalmente separazioni laminari. Inoltre la probabilità di avvenimento di code allegate aumenta con Re e la cavitazione allegata diventa più estesa. Con l’aumento ulteriore di Re, le bolle tendono ad innescare cavità allegate sull’intera scia della bolla, come si vede nelle due foto in basso di Figura 5.5. Inoltre, le cavitazioni allegate tendono a rimanere a lungo tempo dopo la sparizione della bolla principale. 78 Figura 5.3: s mpi i profili h illustrano l’insta ilità llo strato liqui o on una olla i avitazion in movim nto. sp rim nti on una forma i iam tro i 5.08 m a σ 0.45 una velocità di 8.7m/s. Il flusso è da destra a sinistra. Ceccio e Brennen (1991) Questo aumento del numero di eventi di cavitazione con la variazione del numero di Reynolds implica una ricca complessità nella microfluidomeccanica delle bolle di cavitazione. L’importanza di questi diversi tipi di eventi, risiede nel fatto che causano differenze nel processo di collasso, che a sua volta altera la rumorosità e il danno potenziale. Ad esempio, si è notato che gli eventi di cavitazione con code allegate producono meno rumore degli eventi senza code. Suggerendo alcune strategie per la riduzione del rumore e/o del danno potenziale. 5.5 Cavitazione in larga scala Quando la densità degli eventi di cavitazione diventa grande abbastanza, essi iniziano ad interagire alterando il flusso in modo significativo. Questo aumento di densità può avvenire a causa della diminuzione del numero tipico di cavitazione, che causa l’attivazione di un maggior numero di nuclei più piccoli, oppure dall’aumento della popolazione di nuclei per il flusso che sta arrivando. 79 Un esempio dello spostamento di bolle di cavitazione di elevata densità si vede in Figura 5.6, dove sembra che le bolle si fondano formando una scia unica di vapore verso il bordo d’uscita dell’aletta idrodinamica. La cavitazione in grande scala che si forma in uno dei due modi sopra citati, può assumere diverse forme. Figura 5.4: Foto che mostrano la formazione delle code attaccate dietro alla bolla di cavitazione in movimento. La foto in basso mostra la persistenza delle code dopo il collasso della bolla principale. Esperimenti con una forma ITTC di diametro di 5.08cm a σ 0.42 una v lo ità i 9m/s. Il flusso è da destra a sinistra. Ceccio e Brennen (1991) 5.6 Cavitazione in vortici Molti flussi con alti numeri di Reynolds di importanza pratica contengono una regione di concentrata vorticità, dove la pressione nel cuore del vortice è spesso significantemente piccola rispetto al resto del flusso. Questo è il caso ad esempio del vortice di un’elica navale, o di una girante di pompa o nel flusso vorticoso presente nel tubo di scarico di una turbina idraulica. Si ha che l’incipiente cavitazione spesso avviene in questi vortici, e con ulteriore riduzione del numero tipico di cavitazione, il cuore del vortice può diventare pieno di vapore. Nelle Figure dalla 5.7 alla 5.12 sono riportati alcuni esempi di questi tipi particolari di cavitazione in larga scala. La Figura 5.7 mostra le foto di un flusso vorticoso cavitante su una aletta idrodinamica con un certo angolo di attacco. In questi esperimenti l’incipiente cavitazione avviene nel vortice ad un σ 1.4, con ulteriore diminuzione della pressione la cavitazione nel cuore diventa continua come sulla foto di destra. Questa transizione è probabilmente innescata da un accumulo di bolle nel centro del vortice che tendono a migrare verso il cuore per il gradiente di pressione centrifuga. Con ulteriore diminuzione del numero di 80 cavitazione, possono apparire cavitazione in bolle e/o in fogli sulla superficie dell’aletta, come si vede in figura. Figura 5.5: Eventi di cavitazione tipici che mostrano: In alto a sinistra la “foss tta”, in alto a stra l o all gat , al ntro in asso, la avità all gata sull’int ra s ia. sp rim nti su un orpo hi on un iam tro i 50.8 m on σ 0.605 v lo ità i 15m/s. l flusso è a stra a sinistra. Kuhn de Chizelle et al. (1992) Figura 5.6: Denso movimento di bolle di cavitazione sulla superficie di una NACA 4412 con angolo di incidenza nullo a velocità di 13.7m/s e numero di cavitazione di 0.3. Il flusso è da sinistra a destra. Kermeen (1956) 81 Figura 5.7: foto di un flusso vorticoso cavitante su una aletta idrodinamica ellissoidale con un certo angolo i atta o. ulla sinistra si ha la formazion i una ontinua avità a flusso vorti oso, a σ 1.15, on angolo ’attacco di 7.5°. Sulla stra appar a σ 0.43 on angolo ’atta o i 9.5°. Higuchi, Rogers e Arndt (1986) Il flusso vorticoso cavitante è abbastanza visibile nella girante non rivestita della turbopompa in Figura 5.8. Quando si formano vortici cavitanti continui alle estremità delle pale di un elica, essi creano una sorprendente struttura di flusso stabile, come in Figura 5.9, dove si vede che persistono per una lunga distanza dopo il passaggio dell’elica. Chiaramente la cavitazione può avvenire in ogni vortice, e le Figure 5.10 e 5.11 sono un esempio eclatante. La Figura 5.10 mostra la cavitazione nel flusso vorticoso nel tubo di scarico di una turbina Francis. Lo spargimento tridimensionale di vortici da un dispositivo può spesso portare alla formazione e propagazione di un anello di vortici con un nucleo di vapore/gas. La Figura 5.12 mostra questo anello di vortice cavitante che è appena emerso dalla regione della chiusura di una cavità allegata da una lamina oscillante. Altro esempio è riportato in Figura 5.13 dove il vortice è causato dall’oscillazione naturale di una lamina. L’anello di vortice cavitante ha la propria velocità di propagazione relativa al fluido circostante, e ha quindi spostato al di sopra il resto della scia nel momento in cui la foto è stata scattata. Figura 5.8: Flusso vorticoso cavitante su di un modello in scala della girante della turbopompa nello Space Shuttle Main ngin . n a qua on σ 0.42. raist (1979) 82 Figura 5.9: Cavitazione nel flusso vorticoso di un modello di elica. Nethelands Maritime Research Institute and Lips B.V. Figura 5.10: Vortice cavitante nel tubo di scarico di una turbina Francis. Brennen (1995) 83 Figura 5.11: Vortice cavitante nella separazione della scia di un disco piatto con un lembo. Il flusso è da destra a sinistra. A.J. Acosta (1973) Figura 5.12: Formazione di un anello di vortice cavitante nella regione di chiusura della cavità allegata, su un’al tta os illant on or a i 152mm. ’angolo i in i nza os illa tra 5° e 9° alla frequenza di 10Hz. Flusso da sinistra a destra alla velocità di 8.5m/s ed un numero di cavitazione medio di 0.5. D.P. Hart Figura 5.13: Separazion ll’an llo i vorti alla avitazion parzial p r os illazion idrodinamica. Flusso da destra a sinistra. A.J. Acosta (1973) 84 i un’al tta 5.7 Nuvole di cavitazione In molti flussi di interesse pratico si osserva la periodica formazione e collasso di una “nuvola” di bolle di cavitazione. Questa struttura, come già visto in precedenza, è detta “nuvola di cavitazione”. La periodicità temporale può avvenire in modo naturale come risultato del distacco di vortici cavitanti (vedi Figura 5.11) o come risposta ad un periodico disturbo imposto sul flusso. Esempi comuni di fluttuazioni imposte sono le interazioni tra lo statore e il rotore delle pale di turbine o pompe, o l’interazione tra l’elica e la scia non uniforme creata dallo scafo di una nave. In molti di questi casi il collasso coerente delle nuvole (Figura 3.14) può causare rumore più intenso e danni potenziali. Si sono studiati i modelli di flusso su di un’aletta oscillante. Questi studi sono esemplificati nelle foto di Figura 5.14, che mostra la formazione, separazione e collasso di una nuvola di cavitazione su un’aletta oscillante in verticale. Tutti questi studi enfatizzano che si sente un “botto” al momento del collasso della nuvola; in figura questo accade tra la foto centrale e quella di destra. Le nuvole di cavitazione continuano ad essere una preoccupazione principale per eliche e pompe, e sono continuamente soggetto di ricerca. Figura 5.14: Fotografie che mostrano formazione, separazione e collasso di una nuvola di cavitazione sulla sup rfi i i aspirazion i un’al tta (152mm i or a) oscillante alla frequenza di 5.8Hz, e angolo di incidenza di 5°. Il flusso è da sinistra a destra con velocità di 7.5m/s e il numero di cavitazione medio di 1.1. E. McKenny 5.8 Fogli di cavitazione e cavitazione allegata Un’altra classe di strutture di cavitazione in larga scala è quella che avviene quando una scia o una regione di flussi separati si riempie di vapore. Con riferimento alla Figura 5.6 si osservano densi movimenti di bolle quando l’angolo d’ingresso è piccolo. Ad angoli più grandi di 10° (o minori di -2°) la cavitazione avviene come singola zona di separazione piena di vapore come illustrato in Figura 5.15. Questa forma di cavitazione su un’aletta idrodinamica è di solito detta “foglio di cavitazione” e nel contesto di pompe è conosciuta come “pala di cavitazione”. I corpi tozzi esibiscono spesso un’improvvisa transizione da bolla di cavitazione in movimento a una singola scia piena di vapore, come il numero di cavitazione diminuisce. Un esempio è mostrato in Figura 5.16 che mostra due foto di sfere cavitanti; la transizione avviene quando la scia piena di bolle (sulla sinistra) improv85 visamente diventa un singolo vuoto pieno di vapore come mostrato nella foto di destra. Nel contesto di corpi tozzi la scia piena di vapore è spesso chiamata cavità completamente sviluppata o cavità allegata. Chiaramente fogli, pale, cavità completamente sviluppate o allegate sono termini per la stessa struttura di cavitazione in larga scala. Figura 5.15: Foglio di cavitazione sulla superficie di una NACA 4412 ad un angolo ’ingr sso i 12°, con numero di cavitazione di 1.05. Il flusso è da sinistra a destra, con velocità di 10.7m/s. Kermeen (1956) Quando il bordo d’attacco è definito da uno spigolo vivo, la cavità completamente sviluppata è spesso simile ad un vetro liscio, poiché la separazione dello strato limite è laminare. Questa superficie liscia può essere vista nella foto a destra di Figura 5.16 e nelle foto di Figura 5.17. La forma del bordo d’ingresso del corpo può far si che lo strato limite interfacciale passi velocemente a strato turbolento, come nel caso della foto della cavitazione ogiva in Figura 5.17 e della sfera cavitante di Figura 5.16. Figura 5.16: Sfera cavitante con diametro di 7.62cm. La fotografia di sinistra mostra la cavitazione di bolle e la formazione di una scia piena di bolle prima della transizione alla cavità allegata completamente sviluppata mostrata sulla sinistra. Il flusso è da destra a sinistra con velocità rispettivamente di 5.6m/s e di 10.7m/s. Brennen (1970) 86 Per altre forme la transizione può essere ritardata indefinitamente come nel caso del disco cavitante della foto a destra di Figura 5.17. Quando non c’è lo spigolo vivo per formare una cavità completamente sviluppata possono accadere molti fenomeni differenti. La separazione della cavitazione può avvenire lungo una linea ben definita e stabile sulla superficie del corpo, come visto a destra di Figura 5.16, oppure la linea può essere interrotta come in Figura 5.18. Possono esserci altre forme di cavitazione diverse da quelle illustrate, come accade nel caso di Figura 5.19, dove le cavità si manifestano come striature, anche se non sono molto chiare le ragioni per la periodicità trasversale di queste ultime figure. Figura 5.17: Due cavità completamente sviluppate su una forma ogiva di diametro 5.95cm sulla destra e su un disco di diametro di 7.62cm sulla destra. Il flusso è da destra a sinistra con velocità rispettivamente di 7.62m/s e 10.7m/s. Brennen (1970) Figura 5.18: Foglio di cavitazione su una forma ITTC. Il flusso è da sinistra a destra con velocità di 12.2m/s, e numero di cavitazione di 0.424. A.J. Acosta (1973) 87 Figura 5.19: triatur i avitazion su un’al tta i onv ssa alla v lo ità i 15.5m/s da sinistra a destra. V.H. Arakeri (1975) σ 0.11. l flusso è 5.9 Lamine cavitanti Su una lamina idrodinamica la cavitazione allegata può assumere un diverso numero di forme. Quando, come visto in Figura 5.20, la cavità allegata è vicina alla superficie di aspirazione della lamina, la condizione è di cavitazione parziale. Questa è la cavitazione allegata più comunemente osservata su eliche e pompe. A più bassi numeri di cavitazione, la cavità può essere molto vicina al bordo d’uscita della lamina, come in Figura 5.20 in basso. Questa configurazione è detta supercavitazione, ed eliche per navi ad alta velocità sono spesso designate per operare sotto queste condizioni. Tra questi regimi gli esperimenti hanno mostrato che quando la lunghezza della cavità è vicina alla lunghezza della lamina, il flusso diventa instabile e la dimensione della cavità fluttua abbastanza violentemente tra questi limiti. Figura 5.20: Tipi di cavitazione allegata su una lamina inclinata. (a) cavitazione parziale, (b) supercavitazione. Brennen (1995) 88 Durante questo ciclo di fluttuazione la cavità si allunga piuttosto dolcemente, mentre si accorcia con un processo di “pinching-off” formando una grande nuvola di bolle dal dietro della cavità, collassando violentemente come descritto in precedenza. Nelle pompe e in altri dispositivi, queste condizioni di cavitazione chiaramente necessitano di essere evitate, pervia del danno potenziale che possono causare ai materiali. 5.10 Chiusura di cavità Il flusso nelle vicinanze della chiusura della cavità necessita di alcuni commenti. Il flusso in quella ragione è invariabilmente turbolento, poiché lo strato limite che stacca dal corpo produce uno strato limite interfacciale con la superficie libera. Il livello di turbolenza in questo strato, cresce velocemente all’avvicinarsi della chiusura della regione, così il flusso in quella zona appare come un movimento schiumoso di miscelazione turbolenta. Dove le due superfici libere di flusso collidono, una parte del flusso ritorna indietro nella cavità. I cambiamenti della struttura del flusso nella regione di chiusura possono avvenire in flussi orizzontali quando la spinta idrostatica diventa significante, e questi sono i casi in cui il numero di Froude basato sulla lunghezza della cavità, l, Fr = U∞/(gl)1/2, risulta minore di qualche valore critico FrC. Quando Fr < FrC le strutture con getti rientranti non si verificano. Quando invece nella regione di chiusura si forma una coppia di vortici controrotanti con gas/vapore si ha una maggiore stabilità e minore turbolenza rispetto al caso dei getti rientranti. Tornando alla formazione dei getti rientranti sulla chiusura, si nota che questa tipologia di flussi possono esibire delle fluttuazioni significanti, causate dai distacchi ciclici di vortici dal dietro della cavità. Una volta che il cuore dei vortici condensa, solo piccole bolle di gas con componenti incondensabili residui vengono portate via dalla scia. Ci sono molte instabilità che possono innescare o promuovere questi processi ciclici di distacco. Il fenomeno di pulsazione si ha quando grandi cavità si creano dalla fornitura di gas incondensabile dalla scia del corpo, queste cavità che sono indistinguibili dalla controparte di quelle piene di vapore, sono conosciute come “cavità ventilate”; quando il gas aumenta al punto in cui il trascinamento nella regione di chiusura non è più in grado di portar via quel volume di gas, la cavità comincia a fluttuare, creando un processo di pinching-off e ricreando la cavità. Infine, un altro processo che si può trovare nella regione di chiusura nel caso di cavità con predominanza di vapore, è la condensazione istantanea che provoca un meccanismo di chiusura della cavità. Sia le fluttuazioni in larga scala sia le turbolenze in piccola scala nella regione di chiusura, agiscono trascinando le bolle e rimuovendole dalla cavità. È chiaro che il meccanismo preciso di trascinamento può differire considerevolmente da una configurazione ad un’altra. Misure del volume di trascinamento per grandi vuoti con il tipo di chiusura a getto rientrante, suggeriscono che il volume aumenta in relazione alla velocità del flusso, nella forma U∞n, dove n è poco più grande dell’unità. 89 Usando forme assialsimmetriche di diversa dimensione, b, Billet e Weir (1975) mostrarono che il volume trascinato era in scala con U∞b2, che è significantemente variabile con σ, il volume aumenta con la diminuzione del numero tipico di cavitazione e la cavità diventa più grande. Sotto condizioni stabili, la rimozione di vapore e gas non condensabile per trascinamento nella regione di chiusura è bilanciato dal processo di evaporazione e dal rilascio di gas dalla soluzione lungo la superficie libera. Questi processi forniti, sono influenzati dallo stato dello strato limite interfacciale. Uno strato turbolento valorizza il processo di diffusione di calore e di massa che produce l’evaporazione e il rilascio di gas dalla soluzione. Una delle conseguenze del bilancio tra la fornitura di gas incondensabile (aria) e la sua rimozione da trascinamento è l’inerente regolazione della pressione parziale del gas non condensabile nella cavità. Brennen (1969) mise insieme un modello semplificato di questi processi, e mostrò che i risultati per la pressione parziale dell’aria era in accordo approssimativamente con le misure sperimentali della pressione parziale. È presente, inoltre, un bilancio analogo di calore, in cui il calore latente rimosso dal processo di trascinamento deve essere bilanciato dalla diffusione di calore alla cavità attraverso lo strato interfacciale. Questo richiede una temperatura della cavità più bassa di quella del liquido circostante, essa è solitamente piccola in acqua a temperature normali, e può essere significante a più alte temperature o in altri liquidi a temperature simili a quelle in cui le bolle singole sperimentano significanti effetti termici sulla crescita. Conclusioni Come accennato al termine del paragrafo 3.6, non sempre la cavitazione è un fenomeno negativo, ad esempio nella medicina viene utilizzato vantaggiosamente. La cavitazione medica è un trattamento non chirurgico e non invasivo, finalizzato a ridurre lo spessore delle adiposità localizzate attraverso l’utilizzo di ultrasuoni a bassa frequenza. Essa rappresenta una valida alternativa non invasiva alla liposuzione, che è a tutti gli effetti un intervento chirurgico, basato su un sistema ad ultrasuoni progettato per fare implodere, e successivamente eliminare, le cellule adipose. Purtroppo però, come si è visto, la cavitazione è uno dei limiti principali delle macchine a fluido incomprimibile e deve essere tenuta in considerazione in fase di progetto. A tale scopo si utilizza il valore dell’NPSH (Net Positive Suction Head) che è definito come l’energia netta alla bocca d’aspirazione quando la macchina è nelle condizioni di incipiente cavitazione. Questa definizione è data nel caso di macchine operatrici, ma l’NPSH si utilizza anche nella progettazione delle turbine, e si usa definirlo come l’altezza totale al netto della tensione di vapore (Ventrone 2006). I valori dell’NPSH sono dati dalla seguente formula (1+λm)c22/2g + λww22/2g, dove λm e λw sono fattori di carico palare, c2 è la velocità assiale della corrente nella zona di minor pressione e w2 è la velocità relativa della corrente sempre nella zona di 90 minor pressione. Con il termine minor pressione si intende ad esempio, nel caso di pompe centrifughe, la zona all’ingresso della girante, e nel caso di turbine centripete la zona all’uscita della girante o nel tubo di scarico. Il valore dell’NPSH impone un’altezza massima sopra battente alla quale porre la macchina affinché non caviti. Quindi, per ovviare al fenomeno della cavitazione si potrebbe abbassare l’altezza del dispositivo, ma in alcuni sistemi in cui la portata continua a variare, come ad esempio nelle valvole o in impianti in cui le macchine operanti non elaborano sempre la stessa quantità di flusso, la cavitazione può verificarsi in determinate condizioni, e non potendo agire continuamente sull’altezza della macchina è necessario monitorarla. A tale proposito l’ENEA (agenzia nazionale per le nuove tecnologie, l’energia e lo sviluppo economico sostenibile) ha depositato, nel 2012, un nuovo brevetto che consiste in un dispositivo in grado di rilevare i segnali provenienti da trasduttori elettrici sensibili alle vibrazioni prodotte dai fenomeni di cavitazione che avvengono in un fluido e di effettuarne l’analisi qualitativa e quantitativa della frequenza e dell’andamento nel tempo. Il dispositivo, denominato CASBA 2012, rileva la formazione di bolle di gas in un fluido in rapido movimento, misura la frequenza con cui esse si formano e l’andamento nel tempo del fenomeno, quindi può essere utilizzato per la diagnostica e la protezione degli impianti industriali soggetti a danneggiamento derivante dalla formazione e implosione di bolle costituite dal fluido stesso in fase gassosa e da eventuali altri gas disciolti nel fluido, che creando dei microgetti ad altissima velocità di liquido possono danneggiare in maniera significativa gli organi della macchina. La cavitazione, ad esempio, diventa molto difficile da gestire per un sommergibile, dove la velocità dell’elica oltre la quale si innescano fenomeni di cavitazione è funzione della profondità, della salinità e della temperatura dell’acqua, tutti fattori che possono variare molto rapidamente. Per un sommergibile la cavitazione non è tanto pericolosa per gli eventuali danni all’elica ma perché il suono delle implosioni delle bolle di gas è facilmente individuabile, facendo così perdere al sommergibile la sua caratteristica fondamentale, quella di essere “invisibile”. Il CASBA 2012 misura e analizza i segnali elettrici provenienti da sensori piezoelettrici, installati in prossimità della zona dove si sospetta l’insorgenza del fenomeno, sensibili alle vibrazioni prodotte durante la cavitazione del liquido. Il sensore trasforma le vibrazioni in segnali elettrici che, attraverso un cavo coassiale, vengono elaborati e rappresentati su di un display. Il CASBA 2012 può collegare fino a 4 canali di acquisizione dati, dove in ognuno il segnale è campionato da una coppia di circuiti funzionanti in parallelo con un preciso sfasamento temporale, consentendo di aumentare l’accuratezza delle misure e di individuare con maggiore precisione la sorgente primaria e i meccanismi di propagazione dei fenomeni. La cavitazione è quindi un punto fondamentale nella progettazione di macchine e dispositivi idraulici, e deve essere accuratamente prevista per evitare danni ai materiale, che possono essere molto gravosi, e nel caso della mancanza di previsione del fenomeno è necessaria la sua monitorazione con dei dispositivi, come il CASBA 2012, che sono oggigiorno molto sviluppati tecnologicamente, e altamente precisi. 91 BIBLIOGRAFIA e SITOGRAFIA Acosta A.J., 1973, “Hydrofoils and Hydrofoil Craft”, Ann. Rev. Fluid Mech., 5, 161–184 Blake F.G., 1949, “The Onset of Cavitation in Liquids”, I. Acoustics Res. Lab., Harvard Univ., Tech. Memo. No. 12 Brennen C.E., 1995, “Cavitation and Bubble Dynamics”, Oxford University Press, Oxford Ceccio S.L. and Brennen C.E., 1991, “Observations of the Dynamics and Acoustics of Travelling Bubble Cavitation”, J. Fluid Mech., 233, 633–660. Frenkel J., 1955, “Kinetic Theory of Liquids”, Dover, New York Frost D. and Sturtevant B., 1986, “Effects of Ambient Pressure on the Instability of a Liquid Boiling Explosively at the Superheat Limit”, ASME J. Heat Transfer, 108, 418–424 Kermeen R.W., 1956, “Water Tunnel Tests of NACA 4412 and Walchner Profile 7 Hydrofoils in Non-Cavitating and Cavitating Flows”, Calif. Inst. of Tech. Hydro. Lab. Rep. 47-5 Lauterborn W. and Bolle H., 1975, “Experimental Investigations of Cavitation Bubble Collapse in the Neighborhood of a Solid Boundary”, J. Fluid Mech., 72, 391–399 Plesset M.S. and Ellis A.T., 1955, “On the Mechanism of Cavitation Damage”, Trans. ASME, 1055–1064 Plesset M.S. and Mitchell T.P., 1956, “On the Stability of the Spherical Shape of a Vapor Cavity in a Liquid”, Quart. Appl. Math., 13, No. 4, 419–430 Schiebe F.R., 1972, “Measurements of the Cavitation Susceptibility Using Standard Bodies”, St. Anthony Falls Hydr. Lab., Univ. of Minnesota, Rep. No. 118 Ventrone G., 2006, “Macchine per Allievi Ingegneri”, seconda edizione, Libreria Cortina, Padova www.ENEA.it www.lacavitazionemedica.it www.wikipedia.org 92
© Copyright 2024 Paperzz