Formulação Diferencial das Equações de Transporte

67
IM250 Prof. Eugênio Rosa Formulação Diferencial das Equações de Transporte

description

Formulação Diferencial das Equações de Transporte. Forma Integral das Equações de Transporte. O TTR permite escrever as Equações de Transporte a partir do conceito de V olume de C ontrole:. J e f são fontes genéricos associados a SC e ao VC. Teorema de Gauss. - PowerPoint PPT Presentation

Transcript of Formulação Diferencial das Equações de Transporte

Page 1: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Formulação Diferencial

das Equações de Transporte

Page 2: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Forma Integral das Equações de Transporte • O TTR permite escrever as Equações de Transporte

a partir do conceito de Volume de Controle:

b( / )B M

Source

Massa 1 0

Movimento V

1a Lei e

2a Lei s

S Source

VCSCVC SCr dfdAJ dAVnd

dt

d

CVCS

dgndA

T

VCSCSC

k dqdAVndAnq

T

SC VC

k PsdT

qdAn

T

q

J e f são fontes genéricos associados a SC e ao VC

Page 3: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Teorema de Gauss

• O Teorema de Gauss transforma a avaliação de uma integral de superfície em integral de volume.

• Ele aplica-se a grandezas escalares, vetoriais e tensorias:

VCSC

VCSC

VCSC

d Adn

dVAdVn

d Adn

TT

é o operador nabla, f é o gradiente de um escalar (vetor); xV é o rotacional de um vetor (vetor) e .T é o divergente de um tensor (vetor)

Page 4: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Aplicação do Teorema de Gauss

• Aplicando o Teorema de Gauss à Equação de Transporte vamos transformar os termos de superfície em volume:

VC

S Source

VCSCVC SCr

0dfJVdt

d

Gauss de Teorema

dfdAJ dAVnddt

d

A transformação é válida para V.C. não deformáveis, isto é, seu volume não varia com o tempo.

Page 5: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Forma Diferencial

• Como representação Integral acima o tamanho do VC é arbitrário, para a identidade ser válida para qualquer volume é necessário que seu argumento seja nulo!

VC 0dfJV

t

Volume

de fonteSuperfície

de fonte convectivotransiente

fJVt

Page 6: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Diferencial da Massa

• A equação da Massa é obtida fazendo-se b = 1 e J = f = 0,

• Note que para fluidos incompressíveis, isto é, r constante, ela se reduz para:

0Vt

0V

Page 7: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Diferencial da Massa

• Desmembrando o segundo termo da equação vamos encontrar:

• Para regime permanente e um fluido incompressível, a sua densidade não varia ao longo de uma linha de corrente, logo Dr/dt =0 portanto:

0Vt

0V

0VDt

D ou 0VV

tDtD

Veja discussão sobre escoamento estratificado no material do curso

Page 8: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Forma Conservativa e Não-Conservativa

• A equação de transporte acima está na sua forma Conservativa. Os termos transiente e convectivos podem ser desdobrados :

• Nota-se que a forma Conservativa mantinha implicitamente a equação da massa. Após a simplificação chega-se a forma Não-Conservativa

fJVt

fJVt

Vt

0

fJVt

Page 9: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Derivada Substantiva ou Total

• Em cinemática o termo acima tem um significado especial.

• Ele coincide com a taxa de variação de uma propriedade seguindo uma partícula, isto é, a partir de um referencial Lagrangeano.

Vt

VtDt

D

Page 10: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Diferencial da Q. Movimento

• A equação da Q. Movimento é obtida fazendo-se b =V, J = T e f = rg,

• A Equação da Q. Movimento é vetorial, possui 3 componentes,

• Todos os termos possuem unidades de Força/Volume (N/m3)

• O termo rVV é um produto diádico, possui natureza tensorial e representa o fluxo de Q. movimento que cruza a S.C.

gVVt

V

T

Page 11: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Diferencial da Q. Movimento Forma Não-Conservativa

• Desmembrando os termos de transporte e eliminando a equação da massa encontra-se:

• A derivada total da velocidade DV/Dt dá a aceleração seguindo uma partícula!

• Note que a derivada total resgata o conceito da análise de Sistemas pois ele segue uma partícula infinitesimal com identidade fixa!

gVVt

V

T

g

Dt

VD ou gVV

t

V

TT

Page 12: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

PARTIÇÃO DO TENSOR DE TENSÕES

• O Tensor de tensões é decomposto em duas parcelas: uma correspondendo a tensão hidrostática e outra à tensão viscosa.

• A primeira parcela o tensor é isotrópico e refere-se a pressão estática do fluido.

• A segunda parce o tensor é simétrico e deve-se ao movimento relativo das partículas de fluido e a viscosidade; assim

• ou em notação indicial:

P T T

ij ij ijP T T

Page 13: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Transporte da Energia Cinética, K

• Multiplicando-se ambos os lados da Eq. NS por V vamos encontrar:

DVV V P V V g

Dt

T

KDt

DVV

2

1

Dt

DV

Dt

DV

DK V P V V g

Dt

T

• A energia cinética por unidade de massa K é:

• E sua equação de transporte é:

Page 14: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação da Diferencial da Energia ‘e’

• A equação da Energia é obtida fazendo-se b = e, J = -qk + T.V e f = q’’’;

• O lado esquerdo representa o transporte da energia.• O lado direito representa os termos de calor e trabalho

(1a lei) e também um fonte de energia volumétrico.• Na forma não conservativa:

k

eVe q V q

t

T

k

Deq V q

Dt

T

Page 15: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Transporte de ‘e’

• Neste estágio é conveniente substituir T = -P+T’ e expandir os termos na eq. da energia. Vamos utilizar a forma não-conservativa:

qqVVVPPV

Dt

De

qqVVPDt

De

k

VVP

k

T

:TT

T

T’:V é o produto ‘escalar’ entre o tensor desvio da tensão e o tensor deformação do fluido, seu resultado é um escalar. Veja definições no material impresso do curso

Page 16: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Modos de Energia ‘e’

• Vamos considerar três modos de energia: interna, cinética e potencial:

• onde û é a energia interna, g a aceleração da gravidade e r o vetor posição

rgVV2

1ue

ˆ

VgrDt

Dg rg

VDt

DVVV

2

1

Dt

DVV

2

1

Dt

uD u

ˆ

ˆ

VgDt

VDV

Dt

uD

Dt

De

ˆ

• A derivada total em termos das parcelas de ‘e’ fica sendo:

Page 17: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Transporte da Energia Interna, û

• Subtraindo a Equação da Energia Cinética da Equação de ‘e’ vamos ter:

qqV VP Dt

uD

gVV PV Dt

VDV-

qqVVVPPVgVDt

VDV

Dt

uD

k

k

VVP

:Tˆ

T

:TTˆ

T

q VPVq Dt

uDk

:T

ˆ

Page 18: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Diferencial da 2a Lei

• A 2a Lei é obtida fazendo-se b = s, J = -qk/T e f = q’’’/T,

• Os primeiro e segundo termos (lado direito) referem-se à produção ou à destruição de s devido a transferência de calor na fronteira e devido a geração de energia internamente ao volume.

• O último termo refere-se a produção de entropia devido as irreversibilidades do sistema.

• Na forma não-conservativa:

PsT

q

T

qsV

t

s k

kqDs q Ps

Dt T T

Page 19: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

1a e 2a Leis Forma Não-Conservativa

• De maneira similar a equação da massa e Q. de movimento, os termos transiente e convectivos podem ser desmembrados , a equação da massa eliminada e gerando a forma não conservativa da 1a e 2a leis:

PsT

q

T

q

Dt

Ds

qVqDt

De

k

k

T

Page 20: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Função Corrente e Vorticidade

Page 21: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Função Corrente e Eq. Massa (Lagrange 1781)

• A função corrente é um conceito matemático, tal que o escalar deve sempre satisfazer a equação da massa.

• A equação da massa 2D e incompressível, para um sistema cartesiano ou polar, reduz para:

A definição da função corrente é:

Page 22: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

As Linhas de e as Linhas de Corrente

• Para coordenadas cartesianas,

• Para = constante, d = 0 e

• Note que a sua definição coincide com a definição de linha de corrente; somente para regime permanente.

dyudxvdyy

dxx

d

u

v

dx

dyou

dy

v

dx

u

000

Page 23: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Tubo de Corrente

• Duas linhas de corrente definem um tubo de corrente, pq não há velocidade normal às linhas por definição!

• Considere dois tubos espaçados por uma distância Dx,Dy, a vazão volumétrica que passa por eles é:

Dx<0

Dy>0

x

y dxvyudAnVQ

12Q

• A diferença entre duas linhas de corrente define a vazão no tubo!

• Deve-se atribuir um valor a uma única linha de corrente, os valores das demais vem da integração.

2

1

Page 24: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

A Função Corrente y e a Vorticidade

• Para escoamentos 2D a vorticidade se reduz a apenas uma componente:

• Então para um escoamento 2D e irrotacional a função corrente é determinada satisfazendo a equação de Laplace:

2z

z

xxyyx

v

y

u

02

Page 25: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Notas Finais da Parte I

• As equações de transporte, especificamente a Quantidade de Movimento, Energia e 2a Lei estão expressas em função do campo de tensões T.

• Não é possível resolvê-las nesta forma porque não se conhece como o campo de tensão se comporta com o campo de velocidades.

• É necessário estabelecer as equações constitutivas para o fluido onde será modelado como a tensão varia com o campo de velocidades, nosso próximo tópico.

Page 26: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Parte II

Equações Constitutivas

Page 27: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Introdução

• Por equação constitutiva entende-se ‘modelos’ que expressam uma variável em função de outra.

• Por exemplo, a tensão em função da taxa de deformação do fluido.

• Estes ‘modelos’ não são leis físicas mas podem representar sob condições estabelecidas o comportamento físico do fluido.

• Nesta seção serão desenvolvidas equações constitutivas para a – Tensão T no fluido , – Taxa de Calor por condução no fluido, qk.

• Das duas equações a mais envolvente é a equação constitutiva para tensão, vamos começar por ela.

Page 28: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Sobre a Natureza da Tensão T

• As tensões que agem no fluido podem ser Normais ou Cisalhantes;

• Além disto, no estado estático (sem movimento relativo) só agem tensões normais enquanto que para fluido em movimento surgem tensões normais e cisalhantes devido ao atrito (deslizamento) entre as camadas de fluido.

• A tensão T é divida em duas partes, uma devido a pressão P (forças normais) e outra denominada por desvio da tensão, T’ associada ao movimento relativo das partículas no fluido:

TPT

Page 29: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

A Pressão

• A pressão é um tensor isotrópico, isto é, ela não depende da orientação, seus elementos da diagonal são iguais e fora da diagonal são nulos, por isto o tensor pode ser representado por um único escalar:

P000P000P

PPDy

PDx

PDA

Page 30: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Propriedades do Tensor Desvio das Tensões, T’

• O tensor desvio das tensões existe somente se houver movimento relativo entre as partículas de fluido.

• Ele possui tensões normais e cisalhantes,

• Ele é simétrico, isto é, os elementos fora da diagonal são idênticos, T’ij = T’ji

Page 31: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Similaridades Sólido - Fluido• Uma tensão aplicada a um corpo sólido causa uma

deformação, lei de R. Hooke (1635-1703)

• Fluido se deforma continuamente quando sujeito a uma tensão. Newton propôs, por similaridade, que a tensão é proporcional a taxa de deformação

dy

dG

dy

du

dy

dtd

CoeficienteLamé (N/m2) Deformação

viscosidade (N.s/m2)

TaxaDeformação

Page 32: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Viscosidade Dinâmica (Absoluta)

• Fluidos Newtonianos (água, todos os gases e maioria dos líquidos) são aqueles que apresentam uma relação linear entre a tensão e a taxa de deformação.

txy

(N/m2)

du/dy (1/s)

sm

gkou

m

sN

dy/du 2

• A viscosidade m é uma propriedade do fluido e tem natureza escalar.

Page 33: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Extensão para Escoamentos 3D

• A lei de Newton pode ser

estendida para escoamentos

3D a partir do conhecimento

da taxa de deformação

Page 34: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Tensor Deformação, Dij

Em notação indicial, o tensor deformação, Dij, é definido por

Em notação vetorial,

z

w

y

w

x

wz

v

y

v

x

vz

u

y

u

x

u

DDDDDDDDD

x

u

333231

232221

131211

j

iji,D

TT Vou Vgrad

DD

Page 35: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Operação com Tensores

Qualquer tensor pode ser decomposto em uma parte simétrica e outra anti-simétrica:

Simétrico- AntiTensorSimétrico Tensor

2

1

2

1ij,ji,ij,ji,ji, DDDDD

Como T’ é um tensor simétrico ele é proporcional a parte simétrica do tensor Deformação (paralelo a lei de Newton)

Page 36: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Decomposição do Tensor Deformação

1. A diagonal do tensor simétrico está associada a dilatação linear do elemento

2. Os elementos fora da diagonal do tensor simétrico estão associados a deformação angular

3. Os elementos do tensor anti-simétrico estão associados a rotação do elemento fluido.

SIMÉTRICO-ANTI TENSORSIMÉTRICO TENSOR

0z

v

y

w

2

1

z

u

x

w

2

1

y

w

z

v

2

10

y

u

x

v

2

1

x

w

z

u

2

1

x

v

y

u

2

10

z

w

z

v

y

w

2

1

z

u

x

w

2

1

y

w

z

v

2

1

y

v

y

u

x

v

2

1

x

w

z

u

2

1

x

v

y

u

2

1

x

u

Page 37: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

O Tensor, Sij

• O tensor S é a parte simétrica do tensor deformação D.

• Ele existe devido ao movimento relativo do fluido que causa deformações normais e angulares ao elemento de fluido.

TVV2

1S

são tensores que representam o gradiente de velocidades e seu transposto TV e V

Page 38: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Constitutiva para Fluido Newtoniano

• Para fluidos incompressíveis (r constante)

• Para fluidos compressíveis

• Onde I é o tensor identidade

SIT 2P

SIT

SIIT

T

2 V3

2 P

ou 2 V3

2 P

100010001

I

Page 39: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Porque Tensão e Deformação são Linearmente Dependentes?

• A relação t = mdu/dy é um modelo! Portanto não há razão alguma que na natureza os fluidos devam seguir este modelo.

• Entretanto, os gases seguem este modelo;

• Água, óleos em geral e uma grande maioria de líquidos podem ser bem representados por este modelo;

• Mas há líquidos que não são representados: tintas, fluidos biológicos, emulsões em geral.

Page 40: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Fluidos Newtonianos Generalizados

• Eles descrevem fluidos com comportamento não-linear tensão x deformação mas não reproduzem efeitos de:

– tensão normal,– efeitos dependentes do tempo,– ou efeitos elásticos

Page 41: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Fluidos Newtonianos Generalizados

• A relação ‘mais’ geral entre tensão e deformação:

• n – índice de comportamento do escoamento.• k – índice de consistência.

n = 1, fluido newtoniano, k =m

n > 1, fluido dilatante

n < 1 fluido pseudo plástico

Page 42: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

• Shear thickening fluids of this sort are being researched for bullet resistant body armor[2], useful for their ability to absorb the energy of a high velocity projectile impact but remain soft and flexible while worn. Some shear thickening fluids are also used in all wheel drive systems utilizing a viscous coupling unit for power transmission.

• A familiar example of the opposite, a shear thinning fluid, or pseudoplastic fluid, is paint: one wants the paint to flow readily off the brush when it is being applied to the surface being painted, but not to drip excessively.

• There are fluids which have a linear shear stress, shear strain relationship, that requires a finite yield stress before they begin to flow. That is the shear stress, shear strain curve doesn't pass through the origin. These fluids are called Bingham plastics. Several examples are clay suspensions, drilling mud, toothpaste, mayonnaise, chocolate, and mustard. The classic case is ketchup which will not come out of the bottle until you stress it by shaking.

• There are also fluids whose strain rate is a function of time. Fluids that require a gradually increasing shear stress to maintain a constant strain rate are referred to as rheopectic. An opposite case of this, is a fluid that thins out with time and requires a decreasing stress to maintain a constant strain rate (thixotropic).

Page 43: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Effects of Rheology on Non-Newtonian Impact• At impact, the drop experiences a

large shear rate which rapidly decreases as the drop deforms. For non-Newtonian fluids, the effective viscosity of the fluid depends on the shear rate. Therefore one would expect that the effective viscosity of a shear thickening fluid such as cornstarch in water (images in the middle) would behave like a highly viscous fluid at shortly after impact and would gradually spread as a less viscous fluid as time progressed. Sure enough, the cornstarch / water mixture behaves more similar to the highly viscous glycerol (images on the right) than the lower viscosity water (images on the left) at short times. However at later times, one sees that the cornstarch / water fluid spreads out significantly faster than its glycerol counterpart - suggesting a lower effective viscosity. Such experiments provide another method to investigate the rheology of complex fluids.

Page 44: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Viscosidade Aparente, h• É uma conveniência matemática para ajustar a

forma de modelos lineares. • Desmembrando a tensão em um termo linear e

outro com potência (n-1):

• A viscosidade aparente é h = k(du/dy)^(n-1). • Note que ela não é mais propriedade do fluido

mas depende do campo de velocidades.• Ela pode variar ponto a ponto dentro do campo

do escoamento

Page 45: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Constitutiva para Fluido Newtoniano Generalizado

• Para fluidos incompressíveis (r constante)

SIT S2P

• onde S é um escalar com dimensão de (1/s) e é definido pelo produto escalar do tensor S

S:S21S

• e h é uma função tipo lei de potência de S,

1nkS

S:S é o produto escalar entre dois tensores, veja definição em Bird, Stewart and Lightfoot

Page 46: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

REOLOGIA

Fluidos comportamento Não-linear

tensão x deformação

Sólidos comportamento Não-linear

tensão x deformação

Materiais comportamento Visco-elástico

Fluido Newtoniano Comportamento

Puramente Viscoso Linear

Mecânica dos Fluidos

Sólido Hookeano Comportamento

Puramente Elástico Linear

Mecânica dos Sólidos

du/dy

tan =

tan = G

Campo da Reologia

Page 47: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Difusão de Calor, Lei de Fourier

• A condução ou difusão de calor tem natureza vetorial e é dada pela Lei de Fourier:

• onde k é o coeficiente de condução ou difusão térmica, W/moC.

2km

W Tkq

Page 48: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Difusão de Massa, Lei de Fick• O fluxo de massa por difusão de uma espécie

química em outra é proporcional ao gradiente de concentração mássica da espécie :

• onde m’’ é o vetor fluxo de massa (kg/(s.m2);• r é a densidade da mistura;• Dj é o coef. Difusão de massa, (m2/s);

• e wj é a fração mássica ou concentração do componente j, wj = mj/M.

2jjjs.m

kg wDm

Page 49: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Notas Finais

• As equações constitutivas para tensão, calor e massa permitem que as equações de transporte de Q. Movimento e Energia sejam escritas em termos das variáveis básicas: Velocidades, Pressão e Temperatura.

• Na Parte IV desta aula vamos retornar às Equações de Transporte para fazermos esta substituição e chegarmos a sua forma final!

Page 50: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Retorno às Equações Diferencias de Transporte

Page 51: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Navier Stokes

• Substituindo a Eq. constitutiva da Tensão para fluido Newtoniano vamos chegar às Equações de Navier-Stokes (NS):

V

VV P gt

T'

g2V3

2PVV

t

V

S

• onde • A Eq. acima é válida para escoamentos compressíveis,

com viscosidade variável (regime laminar ou turbulento?)• Filmes: (1), (2) e (3).

TVV2

1 S

• A Eq. Transporte de Q. Movimento é:

Page 52: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Navier Stokes Compressível

• vamos chegar às Equações de Navier-Stokes (NS) para um fluido compressível com m constante:

gVV3

1PVV

t

V 2

• Para m constante e considerando a identidade:

VVVV2 2T S

Page 53: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Navier Stokes Incompressível

• Para r e m constantes temos que, .V =0, logo:

• Esta é a forma mais popular das Equações de Navier Stokes: fluido incompressível e com viscosidade constante.

2

2

VVV P V g

t ou

DV P V g

Dt

Page 54: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

A Função Dissipação, • O trabalho realizado pelas tensões para ‘deformar’ o fluido

converte ‘energia mecânica’ do escoamento em ‘energia térmica’ .

• O nome dissipação sugere que em mecânica ‘dissipada’ em térmica, portanto é um termo que introduz irreversibilidades no escoamento.

• Para um fluido Newtoniano ela é definida:

• ou em notação indicial:

• f é a função dissipação, sempre positiva para atender 2a lei.

02V3

2V

2

S:S:T

0x

V

x

V

2

1

x

V

3

22

i

j

j

i2

i

i

222222

2

y

W

z

V

z

U

x

W

x

V

y

U

z

W

y

V

x

U2V

3

2

a função dissipação para coordenadas cartesianas, veja mais detalhes na brochura ‘Forma Dif. Eq. Transporte’.

Page 55: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Transporte da Energia Interna, û

• Substituindo as equações constitutivas para o tensor desvio da tensão e da condução vamos ter:

• Dû/Dt é o transporte de energia interna;

• kT é fluxo calor líquido por condução na S.C.;

• -P.V é trabalho de compressão, fluidos compr.;

• f é a função dissipação, converte trabalho de deformação em energia interna (veja próx slide);

• q’’’ representa geração volumétrica de energia dentro do volume (reação química, radiação outras fontes)

qVPTkDt

uD ˆ

Page 56: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Transporte da Entalpia, h

• O termo do trabalho de pressão pode ser re-escrito em função da equação da massa:

• Substituindo a definição: h = û+P/r na equação de û, chega-se a forma não-conservativa da Equação de Transporte da Entalpia:

• ou a sua forma conservativa:

Dt

DPP

Dt

D

Dt

D1PVP

qDt

DPTk

Dt

Dh

qDt

DPTkhV

t

h

veja mais detalhes na brochura ‘Forma Dif. Eq. Transporte’.

Page 57: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Transporte da Entalpia Total, h0

• A entalpia específica e a entalpia total de um fluido compressível são definidas por:

• Somando à equação da entalpia a energia cinética:

VV21hh e Puh 0

ˆ

gVqTk 2V3

2V

t

P

Dt

Dh

gV 2V3

2VPVVV

2

1

Dt

D

qTk Dt

DP

Dt

Dh

Viscosos Termos

0

S

S

Page 58: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação Transporte da Entalpia Total, h0

• Em geral a entalpia total é empregada para escoamentos compressíveis onde o termo de trabalho das forças de campo é desprezível, neste caso:

• Para tornar sua representação mais compacta é freqüente agrupar os termos viscosos num único operador:

qTk 2V3

2V

t

P

Dt

Dh

Viscosos Termos

0

S

qTk VVt

P

Dt

Dh

Viscosos Termos

0

T

Page 59: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Transporte da Temperatura

• A partir da Equação de transporte da Entalpia e da relação termodinâmica para uma substância pura:

• onde b é o coef expansão volumétrica, • Pode-se mostrar que a forma não-conservativa da

Equação de Transporte para Temperatura é:

• e a sua forma conservativa:

qDt

DPTTk

Dt

DTCP

dP

T1dTC

P

h

T

hdh p

TP

PT

1

qDt

DPTTkTVC

t

TC PP

veja mais detalhes na brochura ‘Forma Dif. Eq. Transporte’.

Page 60: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Equação de Transporte da Entropia• A equação de transporte de S é:

• o termo de produção, Os, é determinado a partir da relação termodinâmica para uma substância pura:

• substituindo as eqs. para h e s na relação acima vamos encontrar:

• As irreversibilidades estão associadas a uma troca térmica com diferença de temperatura ou ao trabalho viscoso realizado pelo fluido

veja mais detalhes na brochura ‘Forma Dif. Eq. Transporte’.

PsT

q

T

Tk

Dt

Ds

Dt

DP1

Dt

DsT

Dt

Dh

dPTdsdh

0

TT

TkPs 2

2

Page 61: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Notas Finais

• Estas são as formas finais de algumas das equações de transporte.

• Há diversas outras que não foram abordadas neste aula, entre elas: transporte de um escalar, e transporte de vorticidade.

• As duas últimas estão na brochura anexa para referência.

• O desafio da próxima aula será simplificar algumas equações e procurar expressá-las numa única Equação Geral de Transporte.

Page 62: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Forma Geral das Equações de Transporte

• O método dos volumes finitos parte da forma conservativa das Eq. Transporte. Considere uma variável escalar f genérica:

SVt

• onde G é o coeficiente difusivo definido por:

T

T

L

L

PrPr

• O fonte S tem natureza diversa: i) representam as condições de contorno do fenômeno;

ii) modelam a ação de forças ou energia de novos mecanismos físicos ou ;

iii) representam todos os outros termos da eq. particular que se quer representar e que não são representados pelo lado esquerdo da equação!

Page 63: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Notação Indicial Eq. Geral de Transporte

SVt

S

xV

xt jj

j

• onde j pode variar de 1 a 3 representando cada uma das direções ortogonais.

• é uma variável escalar genérica e

• A Eq. de Transporte em Notação vetorial

• também pode ser representada em notação indicial pelos operadores

Page 64: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Eq. Geral Escalar e Termos Fontes

onde S

xV

xt jj

j

e k ,w ,h T, h,

PrPr

m0

T

T

qx

PV

t

P

ii

PiiPPi

iP C

1

xx

Tk

C

q

Cx

PV

t

P

C

T

qx

K

Cp

k

x2

x

V

3

2

xV

t

P

jiij

j

j

ii

S

m

kCP

kC

2

1k1

kP

k

CPr

ph

ScD

Prw

1Prk

1Pr

h

T

h0

wm

k

e

k

CPr

pT

k

CPr

ph0

Pr S

Page 65: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Tabela Equações de Transporte

• Obtenha no link abaixo uma tabela extraída do Bird, Stewart and Lightfoot que descreve as componentes das equações de transporte nos sistemas de coordenadas cartesiano, cilindrico-polar e esférico.

• TABELAS EQUAÇÕES DE TRANSPORTE

Page 66: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

Referências[1] White, F.M.; "Viscous Fluid Flow", McGraw Hill (1974)[2] Moore, F.K.; "Theory of Laminar Flows", Princeton Un. Press

(1964)[3] Rosenhead, L.; "Laminar Boundary Layers", Oxford (1963)[4] Warsi, Z.U.A., "Fluid Dynamics: Theoretical and Computational

Approaches", CRC (1993)[5] Panton, R. “Incompressible Flow”, John Wiley (1984)[6] Tennekes, H. and Lumley, J.L., “A First Course in Turbulence”, MIT

Press, 1972,[7] Reynolds W.C. and Perkins, H.C., “Engineering

Thermodynamics”, Mc Graw Hill, (1977)[8] Hinze, J.O., “Turbulence”, McGraw Hill, (1959)[9] Townsend, A.A., “The Strucuture of Turbulent Shear Flow”,

Cambridge Un. Press, 2nd ed., (1976).[10] Wilcox, D.C., “Turbulence Modeling for CFD”, 2nd ed., DCW

Industries, (1998).[11] Astarita, G. and Marrucci, G., “Principles of Non-Newtonian Fluid

Mechanics” , McGraw Hill(1974)

Page 67: Formulação Diferencial das Equações de Transporte

IM250 Prof. Eugênio Rosa

FIM