14
Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório Nacional de Computação Científica

Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

Embed Size (px)

Citation preview

Page 1: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

Simulação de Grandes Escalas com Métodos Estabilizados de Elementos

Finitos

R. C. Almeida e M. JansenCoordenação de Mecânica Computacional

Laboratório Nacional de Computação Científica

Page 2: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório
Page 3: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

Hierarquia dos métodos de solução:

Page 4: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório
Page 5: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

Problemas Transientes

• As Médias de Reynolds não distinguem, e não podem, os efeitos da turbulência e da não estacionaridade

• nenhum tipo de média pode fazê-lo a menos que o efeito transiente seja periódico (ou de alguma forma estruturado) E ocupe um freqüência, ou comprimento de onda, diferente da turbulência

Page 6: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

Problemas Transientes

• Um escoamento separado é muito sensível às condições de contorno e iniciais;

• ‘Neste caso, não é promissor o uso de modelos baseados em médias (de Reynolds)’

Peter Bradshaw

CTR / Stanford

Page 7: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

• Escoamentos que são muito sensíveis às condições de contorno necessitam– simulações 3D – dependência no tempo– DNS, LES, DES ou …

• DES (detached eddy simulation), alimentado por um modelo de médias de Reynolds, é preciso somente se a camada limite for tão fina que detalhes da turbulência não são relevantes

(x,y,z,t)

Page 8: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

Elementos Finitos + Turbulência

• Até recentemente:– restrito à médias de Reynolds– muito caro– pouco preciso

• Mais competitivo:– formulações estáveis de EF– malhas não estruturadas– adaptatividade h e p

Page 9: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

Métodos Estabilizados

Page 10: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

LES + Métodos Estabilizados

• efeito da difusão artificial

• malhas não-estruturadas

• efeito da compressibilidade

• modelagem sub-malha

Page 11: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

0).().(

ut

uDt

D

z

u

x

w

zx

v

y

u

yu

x

u

xx

pX

Dt

Du ).(3

22

z

v

y

w

zx

v

y

u

xu

y

v

yy

pY

Dt

Dv ).(3

22

z

v

y

w

yx

w

z

u

xu

z

w

zz

pZ

Dt

Dw ).(3

22

TkupDt

De).(

Page 12: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

ydyxGtyutxuw

)(),(),(

0).(

ut

i

j

j

i

ji

ji

j

i

x

u

x

u

xx

puu

xt

u 1

jijiijuuuuT

contrário caso 02

x se 1

)( k3xG box

Page 13: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

'~FFF

FF ~

ijllijtijTST

3

12

llTpP

3

1

ijt

jij

i

j

i Sxx

P

x

uu

t

u

21

SCSt

2)( 21

2ijijSSS

Page 14: Simulação de Grandes Escalas com Métodos Estabilizados de Elementos Finitos R. C. Almeida e M. Jansen Coordenação de Mecânica Computacional Laboratório

jijiijtijuuuuS 2

ijijjijiij

TuuuuL

jijiijuuuuT