Universit degli Studi di Perugia Facolt di Ingegneria

  • Slides: 41
Download presentation
Università degli Studi di Perugia Facoltà di Ingegneria Lezioni del corso di TERMOFLUIDODINAMICA E

Università degli Studi di Perugia Facoltà di Ingegneria Lezioni del corso di TERMOFLUIDODINAMICA E IMPIANTI TERMOTECNICI a. a. 2013/2014

n OBIETTIVI Fornire agli allievi le conoscenze in materia di complementi di trasmissione del

n OBIETTIVI Fornire agli allievi le conoscenze in materia di complementi di trasmissione del calore, di termofluidodinamica applicata e di impianti tecnici ai problemi dell’ingegneria meccanica. Fornire agli allievi conoscenze in materia di termofluidodinamica computazionale e di modelli di dispersione di inquinanti in atmosfera.

CONTENUTI n Conduzione: proprietà termofisiche; casi non stazionari; approssimazione di corpo sottile; problemi non

CONTENUTI n Conduzione: proprietà termofisiche; casi non stazionari; approssimazione di corpo sottile; problemi non lineari: integrale di conducibilità; transitori in sistemi a temperatura non uniforme; superfici alettate. n Trasmissione di calore per irraggiamento; metodo della radiosità. n Caratteri della convezione; equazioni di Navier-Stokes; equazione dell'energia nei fluidi, forma adimensionale delle equazioni della convezione; approssimazione di strato limite; valutazione dello spessore degli strati limite meccanico e termico; equazioni dello strato limite; soluzione di similitudine su lastra piana; strato limite in geometrie non piane: separazione. n Flusso laminare in tubi; calcolo delle perdite di carico; convezione laminare nei flussi interni; convezione naturale laminare. n Caratteri della turbolenza; transizione alla turbolenza; struttura della turbolenza; sforzi di Reynolds; lunghezza di rimescolamento; profili di velocità; perdite di carico in flussi turbolenti; tubi scabri; diffusività termica turbolenta; analogia di Reynolds; analogia di Prandtl-Taylor; relazioni di scambio termico in flussi interni.

Segue contenuti n n n Trasmissione del calore per condensazione; condensazione a film; effetti

Segue contenuti n n n Trasmissione del calore per condensazione; condensazione a film; effetti di turbolenza; correlazioni della condensazione a film; condensazione a gocce; effetto degli incondensabili. Ebollizione; curva di Nukiyama; surriscaldamento; crescita delle bolle; flusso critico; regimi di ebollizione in convezione forzata; correlazioni per l'ebollizione in convezione forzata. Scambiatori di calore; metodi della DTLM e dell’efficienza; scambiatori a più passaggi; tipologie di scambiatori di calore. Termofluidodinamica computazionale: metodo delle differenze finite; metodi alle differenze finite nella conduzione; applicazione delle differenze finite a problemi di conduzione; metodo degli elementi finiti; equazioni per il metodo degli elementi finiti nella conduzione stazionaria; applicazione del metodo agli elementi finiti ad un caso di conduzione non stazionaria; metodo agli elementi finiti nella conduzione non a regime; integrazione delle equazioni del metodo agli elementi finiti non stazionario. Applicazioni e casi di studio. Modelli matematici per il calcolo della dispersione di inquinanti in atmosfera: modelli generali, modelli gaussiani, puff model, box model, street Canyon, modelli lagrangiani; modelli specifici per traffico stradale e per sorgenti puntiformi; criteri generali di scelta dei modelli di diffusione di inquinanti in atmosfera. Applicazioni e casi di studio

Segue contenuti n n n n n Benessere termoigrometrico, indici globali, cause di discomfort

Segue contenuti n n n n n Benessere termoigrometrico, indici globali, cause di discomfort locale, Sindrome da edifici insalubri. Processi psicrometrici, analisi del miscuglio aria vapore, caratterizzazione degli stati termodinamici, trasformaszioni del fluido connesse ai trattamenti dell’aria. Unità di Trattamento Aria, elementi fondamentali, criteri di progettazione e dimensionamento, schemi impiantistici. Classificazione degli impianti di climatizzazione e condizionamento, elementi degli impianti a tutt’aria e misti aria acqua. Generatori di calore: caratteristiche principali, elementi di selezione, rendimenti. Gruppi termici tradizionali, a temperatura scorrevole, a condensazione. Fluidi frigorigeni, requisiti caratteristiche ed applicazioni. Classificazione e tipologie delle pompe di calore. Caratteristiche tecniche macchine frigorifere a compressione e ad assorbimento Elementi terminali degli impianti di riscaldamento e condizionamento, radiatori, fan coil, pannelli radianti, termostrisce, aerotermi. Organi di spinta, canali e tubazioni, elementi per la progettazione e il dimensionamento.

Segue contenuti n n n Sistemi di smaltimento del calore (a secco, a umido,

Segue contenuti n n n Sistemi di smaltimento del calore (a secco, a umido, misti) Regolazione degli impianti, dispositivi per la regolazione, algoritmi matematici basilari per il funzionamento degli stessi. Collaudo degli impianti, normativa, requisiti e figure professionali coinvolte. Sistemi di cogenerazione, caratteristiche e criteri di dimensionamento Cenni di Prevenzione incendi, chimica del fuoco, normativa, sistemi di protezione attiva e passiva.

Segue contenuti PREREQUISITI: n Fisica tecnica 1 e 2. n TESTI CONSIGLIATI: n G.

Segue contenuti PREREQUISITI: n Fisica tecnica 1 e 2. n TESTI CONSIGLIATI: n G. Guglielmini, C. Pisoni, Elementi di trasmissione del calore, Ed. Veschi n C. Buratti: Impianti di climatizzazione e condizionamento , Ed. Morlacchi n Saranno inoltre distribuite dispense da parte dei docenti MODALITÀ DI VERIFICA DEL PROFITTO: n La verifica del profitto consiste in una prova scritta e in un colloquio orale della durata di circa 30’.

Modalità di trasmissione del calore CONDUZIONE CONVEZIONE IRRAGGIAMENTO

Modalità di trasmissione del calore CONDUZIONE CONVEZIONE IRRAGGIAMENTO

Conduzione Il manico scotta! La conduzione è il principale sistema di trasmissione di calore

Conduzione Il manico scotta! La conduzione è il principale sistema di trasmissione di calore nei solidi. Afferrando il manico di una pentola riscaldata il calore sarà condotto attraverso il metallo verso la mano.

Convezione L’aria più calda è meno densa, pertanto sale attraverso gli strati più freddi.

Convezione L’aria più calda è meno densa, pertanto sale attraverso gli strati più freddi. La convezione è il principale sistema di trasmissione di calore nei liquidi e gas. L’aria calda sopra la pentola sale verso l’alto poichè è più leggera dell’aria fredda che la sovrasta.

Irraggiamento L’emissione di calore per irraggiamento è generata da tutti gli oggetti che si

Irraggiamento L’emissione di calore per irraggiamento è generata da tutti gli oggetti che si trovano ad una temperatura al di sopra dello zero assoluto. L’irraggiamento è la sola modalità di trasmissione del calore che non richiede materia come mezzo di trasporto. E’ quindi l’unica possibilità di trasferire calore attraverso il vuoto.

LA CONDUZIONE

LA CONDUZIONE

I meccanismi della conduzione 1/2 Lato caldo Calore Lato freddo Guardiamo nel dettaglio cosa

I meccanismi della conduzione 1/2 Lato caldo Calore Lato freddo Guardiamo nel dettaglio cosa avviene alle particelle di materia quando un corpo è riscaldato ad una estremità. Il calore fa vibrare le particelle all’interno del corpo; tali vibrazioni sono trasferite da una particella all’adiacente ed in tal modo il calore è trasmesso attraverso tutto il corpo.

I meccanismi della conduzione 2/2 Lato caldo Lato freddo elettrone In tutti i solidi,

I meccanismi della conduzione 2/2 Lato caldo Lato freddo elettrone In tutti i solidi, la trasmissione del calore per conduzione avviene attraverso due meccanismi: Calore 1. Il calore fa vibrare le particelle, tale movimento è trasferito da una particella all’altra 2. Il “mare di elettroni” esterno che i corpi possiedono (in modo particolare i metalli) acquista energia cinetica all’atto del riscaldamento; nei metalli sono proprio gli elettroni che conducono la maggior parte del calore.

Analisi della conduzione Evidenze sperimentali su pareti piane di spessore << altezza permettono di

Analisi della conduzione Evidenze sperimentali su pareti piane di spessore << altezza permettono di ricavare il calore q” scambiato per unità di tempo e superficie. T Distribuzione della temperatura nel solido T T 0 x

Il postulato di Fourier (1768— 1830) Definisce la quantità di calore che attraversa una

Il postulato di Fourier (1768— 1830) Definisce la quantità di calore che attraversa una superficie infinitesima comunque orientata in direzione ad essa normale n La conducibilità termica, k, è caratteristica del materiale ed è una delle sue proprietà termofisiche.

Alcuni valori di conducibilità termica N. B. : in condizioni normali di temperatura e

Alcuni valori di conducibilità termica N. B. : in condizioni normali di temperatura e pressione.

Variabilità della conducibilità termica con la temperatura SOLIDI

Variabilità della conducibilità termica con la temperatura SOLIDI

Variabilità della conducibilità termica con la temperatura Liquidi non metallici in condizioni di saturazione

Variabilità della conducibilità termica con la temperatura Liquidi non metallici in condizioni di saturazione

Variabilità della conducibilità termica con la temperatura Gas a pressione normale

Variabilità della conducibilità termica con la temperatura Gas a pressione normale

L’equazione di Fourier 1/4 z (x, y, z) L’ elemento di volume infinitesimo d.

L’equazione di Fourier 1/4 z (x, y, z) L’ elemento di volume infinitesimo d. V è centrato nel punto (x, y, z). y x Superficie del sistema z y x

L’equazione di Fourier 2/4 Variazione di energia nel volume V = Somma dei flussi

L’equazione di Fourier 2/4 Variazione di energia nel volume V = Somma dei flussi di calore attraverso la superficie di V + Flusso di calore per generazione interna in V

L’equazione di Fourier 3/4 Questa espressione rappresenta l’EQUAZIONE GENERALE DELLA CONDUZIONE in coordinate cartesiane

L’equazione di Fourier 3/4 Questa espressione rappresenta l’EQUAZIONE GENERALE DELLA CONDUZIONE in coordinate cartesiane per sistemi tridimensionali ρ = densità locale del mezzo C = calore specifico locale del mezzo = generazione interna di calore

L’equazione di Fourier 4/4 Se il mezzo è omogeneo ed isotropo, introducendo la diffusività

L’equazione di Fourier 4/4 Se il mezzo è omogeneo ed isotropo, introducendo la diffusività termica a=k (ρC)-1: in cui il termine ∇2 T rappresenta l’operatore laplaciano della temperatura: Caso di regime stazionario e assenza di generazione interna di calore…. n Coordinate cilindriche n Coordinate sferiche

Parete a simmetria cilindrica 1/9 Cilindro cavo di lunghezza L e raggi r 1

Parete a simmetria cilindrica 1/9 Cilindro cavo di lunghezza L e raggi r 1 ed r 2 r 1 r 2 L Sezione trasversale r r

Parete a simmetria cilindrica 2/9 Ipotesi n Flusso unidimensionale: T = T(r) n Assenza

Parete a simmetria cilindrica 2/9 Ipotesi n Flusso unidimensionale: T = T(r) n Assenza di generazione interna di calore n Regime stazionario r n Mezzo omogeneo ed isotropo Importanti applicazioni n n n Tubi isolati Isolamento di cavi elettrici Scambiatori di calore r

Parete a simmetria cilindrica 3/9 L’equazione generale della conduzione in coordinate cilindriche: si trasforma,

Parete a simmetria cilindrica 3/9 L’equazione generale della conduzione in coordinate cilindriche: si trasforma, con le ipotesi enunciate, come segue: L’integrale generale si esprime come: T(r) = C 1+C 2 ln(r)

Parete a simmetria cilindrica 4/9 Le condizioni al contorno si scrivono: T(r 1) =

Parete a simmetria cilindrica 4/9 Le condizioni al contorno si scrivono: T(r 1) = T 1 T(r 2) = T 2 che, applicate all’integrale generale: dove:

Parete a simmetria cilindrica 5/9 Il profilo di temperatura La soluzione T(r 1) =

Parete a simmetria cilindrica 5/9 Il profilo di temperatura La soluzione T(r 1) = T 1 Profilo logaritmico T(r) T(r 2) = T 2 r

Parete a simmetria cilindrica 6/9 Il flusso di calore per unità di superficie si

Parete a simmetria cilindrica 6/9 Il flusso di calore per unità di superficie si valuta attraverso l’espressione di Fourier: Il flusso di calore che attraverso la generica isoterma è pari a:

Parete a simmetria cilindrica 7/9 Con condizioni al contorno di tipo convettivo all’interno e

Parete a simmetria cilindrica 7/9 Con condizioni al contorno di tipo convettivo all’interno e all’esterno: r 1 Fluido che scorre all’interno a Tf 1 e con coefficiente medio di convezione h 1 r 2 Fluido che scorre all’esterno a Tf 2 e con coefficiente medio di convezione h 2

Parete a simmetria cilindrica 8/9 Rconv 1 Tf 1 Rconv 2 Rcond T 1

Parete a simmetria cilindrica 8/9 Rconv 1 Tf 1 Rconv 2 Rcond T 1 T 2 Tf 2

Parete a simmetria cilindrica 9/9 L’espressione del flusso termico in forma “apparentemente” semplificata è:

Parete a simmetria cilindrica 9/9 L’espressione del flusso termico in forma “apparentemente” semplificata è: La semplificazione scompare nell’esplicitare il coefficiente globale di trasmissione U (caso di parete cilindrica con n strati):

Raggio critico di isolamento 1/3 Ricoprire una tubazione con materiale isolante non porta necessariamente

Raggio critico di isolamento 1/3 Ricoprire una tubazione con materiale isolante non porta necessariamente ad un aumento della resistenza termica complessiva del sistema. Se da un lato si ha una crescita della resistenza per conduzione nell’attraversamento dello spessore, dall’altro, l’aumento della superficie esterna disperdente, fa diminuire la resistenza per convezione. Se si raggruppano in SR le resistenze del cilindro fino allo strato di isolante, si può riscrivere l’espressione del calore scambiato fra il fluido e l’esterno. Pertanto la resistenza termica totale è:

Raggio critico di isolamento 2/3 Da questa relazione si evince che, fissati r 1,

Raggio critico di isolamento 2/3 Da questa relazione si evince che, fissati r 1, r 2, L, k ed h, Rcond aumenta logaritmicamente con r mentre Rconv diminuisce con r secondo una curva iperbolica. L’annullamento della derivata prima della funzione Rt (r) porta alle relazioni: La derivata seconda è: che, calcolata per r = rc: ovvero in rc si ha un punto di minimo.

Raggio critico di isolamento 3/3 Se r 2 > rc l’adozione dell’isolante comporta comunque

Raggio critico di isolamento 3/3 Se r 2 > rc l’adozione dell’isolante comporta comunque un aumento della resistenza termica Se r 2 < rc l’aggiunta di spessore di isolante riduce la resistenza totale, a meno di aumentare lo spessore oltre r 3.

Mezzi a conducibilità dipendente dalla temperatura 1/5 L’equazione generale della conduzione si esprime come:

Mezzi a conducibilità dipendente dalla temperatura 1/5 L’equazione generale della conduzione si esprime come: ovvero, in forma compatta: per la sua integrazione si ricorre alla trasformazione di KIRCHHOFF: T 0 = temperatura di riferimento, k 0 = k(T 0) e ad una funzione f tale che:

Mezzi a conducibilità dipendente dalla temperatura 2/5 e poichè: inoltre: L’equazione generale diventa dunque:

Mezzi a conducibilità dipendente dalla temperatura 2/5 e poichè: inoltre: L’equazione generale diventa dunque: La diffusività termica a non dipende dalla temperatura in molti casi, quindi l’equazione differenziale diventa lineare

Mezzi a conducibilità dipendente dalla temperatura 3/5 CASO MONODIMENSIONALE T 1 T 2 L

Mezzi a conducibilità dipendente dalla temperatura 3/5 CASO MONODIMENSIONALE T 1 T 2 L • Regime stazionario • Assenza di sorgenti di calore • Geometria piana L’equazione della conduzione diventa: x La prima integrazione porta a: la seconda integrazione porta a: (q” di Fourier)

Mezzi a conducibilità dipendente dalla temperatura 4/5 Introducendo la conduttività media km: si può

Mezzi a conducibilità dipendente dalla temperatura 4/5 Introducendo la conduttività media km: si può scrivere che: ed integrando fino allo spessore x: Se è nota k(T) si ottiene l’andamento della temperatura T=T(x) Spesso la dipendenza della conducibilità con la temperatura è di tipo lineare:

Mezzi a conducibilità dipendente dalla temperatura 5/5 La conduttività media diventa: La distribuzione della

Mezzi a conducibilità dipendente dalla temperatura 5/5 La conduttività media diventa: La distribuzione della temperatura risulta del secondo ordine: • γ = 0 distribuzione lineare T 1 • γ>0 T 2 L x • γ<0