Formulação Diferencial das Equações de Transporte
-
Upload
teagan-rich -
Category
Documents
-
view
22 -
download
0
description
Transcript of 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
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)
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.
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
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
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
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
qVqeVt
ek
T
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.
PsT
q
T
qsV
t
s k
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
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
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
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 é:
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
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
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
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
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
IM250 Prof. Eugênio Rosa
Notas Finais da Parte II
• 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.
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.
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
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
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
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
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.
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
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
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)
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
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
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
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.
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
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
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).
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.
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
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
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
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
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
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!
IM250 Prof. Eugênio Rosa
Equação Diferencial da Massa
• Note que para fluidos incompressíveis, isto é, r constante, ela se reduz para:
0Vt
0V
0VVt
ou
DtD
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):
gVVt
V
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 é:
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
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.
gVPVVt
V 2
IM250 Prof. Eugênio Rosa
Equação de Transporte de ‘e’
• A equação de transporte da Energia ‘e’, na sua forma não-conservativa é:
• Neste estágio é conveniente substituir T = -P+T’ e expandir os termos:
qVqDt
Dek
T
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
IM250 Prof. Eugênio Rosa
Equação de Transporte de ‘e’
• Para se chegar a forma final da Equação da Energia é necessário definir:
1. As formas de energia que ‘e’ representa;
2. A difusão do calor, qk
3. O tensor das tensões no fluido e seus produtos
• Estas tarefas serão feitas na seqüência.
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:
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:
gV2V3
2VPV
Dt
VDV
S
KDt
DVV
2
1
Dt
DV
Dt
DV
gV2V3
2VPVK
Dt
D
S
• A energia cinética K é:
• E sua equação 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
ˆ
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’.
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 ˆ
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’.
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
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
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’.
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
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.
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!
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
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
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
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
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)