Estudo Teórico e Simulação Numérica do Escoamento ...

32
Estudo Fenomenológico e Numérico do Escoamento Estratificado Óleo - água Ondulado e com Mistura na Interface Ricardo Pereira de Ávila Orientador: Prof. Assoc. Oscar M. H. Rodriguez Novembro 2015 NETeF (Núcleo de Engenharia Térmica e Fluidos) – EESC – USP

Transcript of Estudo Teórico e Simulação Numérica do Escoamento ...

Page 1: Estudo Teórico e Simulação Numérica do Escoamento ...

Estudo Fenomenológico e Numérico do Escoamento Estratificado Óleo-água Ondulado e com Mistura na

Interface

Ricardo Pereira de Ávila

Orientador: Prof. Assoc. Oscar M. H.

Rodriguez

Novembro 2015

NETeF (Núcleo de Engenharia Térmica e Fluidos) – EESC – USP

Page 2: Estudo Teórico e Simulação Numérica do Escoamento ...

Entranhamento de gotas: Fatores

de Entranhamento

Escoamento Estratificado com

Mistura na Interface (ST&MI)

NETeF-EESC-USP 2

Page 3: Estudo Teórico e Simulação Numérica do Escoamento ...

Esquema do modelo do escoamento estratificado com o entranhamento de gotículas - Extraída de Hadziabdic e Oliemans

(2007).

3NETeF-EESC-USP

Fatores de Entranhamento

Page 4: Estudo Teórico e Simulação Numérica do Escoamento ...

, o dropletso

o

AE

A

Fatores de entranhamento :( ) para a fase óleo, e, ( ) para a fase água.

4NETeF-EESC-USP

Fatores de Entranhamento

oE

, w dropletsw

w

AE

A

wE

Page 5: Estudo Teórico e Simulação Numérica do Escoamento ...

Método de identificação das alturas das interfaces, correspondentes aos pontos de inflexão nas oscilações das curvas de distribuições de fases: 34D, Um = 0,20m/s,

(vermelho) e água (azul) – Extraída de Yusoff (2012) (modificada).

5NETeF-EESC-USP

Fatores de Entranhamento - Yusoff (2012)

Page 6: Estudo Teórico e Simulação Numérica do Escoamento ...

Representação esquemática para o cálculo do fator de empacotamento de gotas bidimensional (FEG-2D).

6NETeF-EESC-USP

Fatores de Entranhamento - Yusoff (2012)

Page 7: Estudo Teórico e Simulação Numérica do Escoamento ...

Comparação entre as alturas de interfaces para diferentes modelos fenomenológicos (RO) Rodriguez e Oliemans (2006) modificado, (RB) Rodriguez e Baldani (2012),

(RO-FE) Rodriguez e Oliemans (2006) com Fator de Entranhamento, (RB-FE) Rodriguez e Baldani (2012) com Fator de Entranhamento.

7NETeF-EESC-USP

Fatores de Entranhamento - Yusoff (2012)

Page 8: Estudo Teórico e Simulação Numérica do Escoamento ...

Comparação entre os fatores de entranhamento da água:(a) sem FEG-2D, (b) com FEG-2D .

8NETeF-EESC-USP

Fatores de Entranhamento - Yusoff (2012)

(a) (b)

Page 9: Estudo Teórico e Simulação Numérica do Escoamento ...

Comparação entre os fatores de entranhamento do óleo:(a) sem FEG-2D, (b) com FEG-2D .

9NETeF-EESC-USP

Fatores de Entranhamento - Yusoff (2012)

(a) (b)

Page 10: Estudo Teórico e Simulação Numérica do Escoamento ...

Comparação entre os gradientes de pressão calculados pelos modelos em relação aos dados experimentais – Rodriguez e Oliemans (2006):

(a) RO-FE e RB-FE; (b) RB e RB-FE.10NETeF-EESC-USP

Fatores de Entranham. – Rod.Oliem.(2006)

(a) (b)

Page 11: Estudo Teórico e Simulação Numérica do Escoamento ...

Fatores de entranhamento de óleo e de água (RB-FE) em função da fração volumétrica de água (dados experimentais) – Rodriguez e

Oliemans (2006).11NETeF-EESC-USP

Fatores de Entranham. – Rod.Oliem.(2006)

(a) (b)

Page 12: Estudo Teórico e Simulação Numérica do Escoamento ...

12NETeF-EESC-USP

Fatores de Entranham. – Rod.Oliem.(2006)

(a) (b)

Áreas de entranhamento de óleo e de água (RB-FE) em função da fração volumétrica de água (dados experimentais) – Rodriguez e Oliemans

(2006).

Page 13: Estudo Teórico e Simulação Numérica do Escoamento ...

Simulações Numéricas em CFD:

escoamento estratificado ondulado

óleo-água

NETeF-EESC-USP 13

Page 14: Estudo Teórico e Simulação Numérica do Escoamento ...

• Escoamento estratificado ondulado óleo-água;

• Teste de malha (em n° de nós):

200.000: holdup e gradiente de pressão (-dp/dz);

• Malhas utilizadas (em n° de nós):

200.000: holdup e gradiente de pressão (-dp/dz);

800.000: transição de malhas;

3.200.000: geometria e parâmetros de onda;

NETeF-EESC-USP 14

Simulações Numéricas em CFD

Page 15: Estudo Teórico e Simulação Numérica do Escoamento ...

Volume de controle e pontos de tomada de pressão (cruzesamarelas em z = 1,5m e z = 1,75m) - óleo (vermelho) e água (azul) Uws = 0,15m/s e Uos = 0,10m/s, horizontal, comprimento do

tubo: 2m (CFX-Post®).

Volume de Controle e Pontos de Pressão

15NETeF-EESC-USP

Page 16: Estudo Teórico e Simulação Numérica do Escoamento ...

NETeF-EESC-USP 16

Resultados: Holdup (óleo)

Fração volumétrica do óleo - Pereira (2011): horizontal, Uws = 0,15m/s , Uos = 0,03m/s , 0,05m/s , 0,07m/s , 0,10m/s , 0,13m/s e 0,15m/s ,(a) previsões, (b) erros.

(a) (b)

Page 17: Estudo Teórico e Simulação Numérica do Escoamento ...

NETeF-EESC-USP 17

Resultados: Holdup (óleo)

Fração volumétrica do óleo - Pereira (2011): -10°,Uws = 0,15m/s , Uos = 0,03m/s , 0,05m/s , 0,07m/s , 0,10m/s , 0,13m/s e 0,15m/s ,(a) previsões, (b) erros.

(a) (b)

Page 18: Estudo Teórico e Simulação Numérica do Escoamento ...

NETeF-EESC-USP 18

Resultados: Holdup (óleo)

Fração volumétrica do óleo - Pereira (2011): +10°, Uws = 0,15m/s , Uos = 0,03m/s , 0,05m/s , 0,07m/s , 0,10m/s , 0,13m/s e 0,15m/s ,(a) previsões, (b) erros.

(b)(a)

Page 19: Estudo Teórico e Simulação Numérica do Escoamento ...

NETeF-EESC-USP 19

Resultados: Gradiente de Pressão

Fração volumétrica do óleo - Pereira (2011): horizontal, Uws = 0,15m/s , Uos = 0,03m/s , 0,05m/s , 0,07m/s , 0,10m/s , 0,13m/s e 0,15m/s ,(a) previsões, (b) erros.

(b)(a)

Page 20: Estudo Teórico e Simulação Numérica do Escoamento ...

NETeF-EESC-USP 20

Resultados: Holdup (óleo) - Erros

Erros médios globais para as previsões da fração volumétrica de óleo - Pereira (2011).

Inclinação [°] Homogêneo Rodriguez e

Oliemans

(2006)

Rodriguez e

Baldani (2012) CFD

0 38,29 48,00 3,87 15,37

-10 39,48 50,94 8,57 15,90

-20 34,50 36,54 14,38 13,18

+10 37,31 43,05 13,62 15,68

+20 30,28 24,30 26,46 5,40

Page 21: Estudo Teórico e Simulação Numérica do Escoamento ...

Inclinação [°] Homogêneo Rodriguez e

Oliemans

(2006)

Rodriguez e

Baldani (2012) CFD

0 283,90 216,50 11,70 142,00

-10 70,47 60,76 1,33 27,40

-20 31,08 29,51 2,23 10,90

+10 39,68 22,93 6,44 19,03

NETeF-EESC-USP 21

Resultados: Gradiente de Pressão - Erros

Erros médios globais para as previsões do gradiente de pressão -Pereira (2011).

Page 22: Estudo Teórico e Simulação Numérica do Escoamento ...

Malhas (em n° de nós): (a) 200.000 , (b) 800.000 , (c) 3.200.000 . (ICEM-CFD®).

Geometria e Parâmetros de Onda

22NETeF-EESC-USP

(a) (b) (c)

Page 23: Estudo Teórico e Simulação Numérica do Escoamento ...

Interface ondulada na região de simetria,seções tranversais (interface convexa e côncava), Óleo (vermelho) e água (azul),Uws = 0,15m/s e Uos = 0,10m/s, horizontal, (CFX-Post®).

23NETeF-EESC-USP

Geometria e Parâmetros de Onda

Page 24: Estudo Teórico e Simulação Numérica do Escoamento ...

Seções tranversais: (a) côncava, (b) convexa, Óleo (vermelho) e água (azul),Uws = 0,15m/s e Uos = 0,10m/s, horizontal, (CFX-Post®).

24NETeF-EESC-USP

Geometria e Parâmetros de Onda

(a) (b)

Page 25: Estudo Teórico e Simulação Numérica do Escoamento ...

Uws = 0,16 m/s Amplitude Desvio

padrão

[mm]

Comprim. Desvio

padrão

[mm]

Celeridade Desvio

padrão

[mm] Uos = 0,02 m/s [mm] [mm] [m/s]

Experimental 4,00 1,93 29,45 11,33 0,11 0,03

CFD

1° 4,21 3,01 138,26 63,99 1,18 0,46

2° 3,15 1,31 110,21 37,75 2,03 0,97

3° 3,21 1,40 106,43 35,64 2,08 0,85

Erros para os parâmetros de onda: Uws = 0,16m/s, Uos = 0,02m/s (LabVIEW®).

Geometria e Parâmetros de Onda - Erros

25NETeF-EESC-USP

Page 26: Estudo Teórico e Simulação Numérica do Escoamento ...

Agradecimentos

• Prof. Assoc. Oscar M. H. Rodriguez;

• Prof. Dr. Marcelo Souza de Castro;

• Marcel Barbosa Cavallini;

• colegas;

26NETeF-EESC-USP

Page 27: Estudo Teórico e Simulação Numérica do Escoamento ...

Simulações Numéricas em CFD:

escoamento estratificado ondulado

óleo-água

NETeF-EESC-USP 27

Page 28: Estudo Teórico e Simulação Numérica do Escoamento ...

NETeF-EESC-USP 28

Resultados: Holdup (óleo)

Fração volumétrica do óleo - Pereira (2011): -20°, Uws = 0,10m/s , Uos = 0,03m/s , 0,05m/s , 0,07m/s , 0,10m/s , 0,13m/s e 0,15m/s ,(a) previsões, (b) erros.

(b)(a)

Page 29: Estudo Teórico e Simulação Numérica do Escoamento ...

NETeF-EESC-USP 29

Resultados: Holdup (óleo)

Fração volumétrica do óleo - Pereira (2011): +20°, Uws = 0,10m/s , Uos = 0,03m/s , 0,05m/s , 0,07m/s , 0,10m/s , 0,13m/s e 0,15m/s ,(a) previsões, (b) erros.

(b)(a)

Page 30: Estudo Teórico e Simulação Numérica do Escoamento ...

NETeF-EESC-USP 30

Resultados: Gradiente de Pressão

Fração volumétrica do óleo - Castro (2013): horizontal, Uos = 0,02m/s , 0,04m/s , 0,06m/s , 0,08m/s , 0,10m/s , 0,12m/s e 0,14m/s ,(a) previsões p/ Uws = 0,13m/s , (b) previsões p/ Uws = 0,16m/s .

(b)(a)

Page 31: Estudo Teórico e Simulação Numérica do Escoamento ...

Uws [m/s] Homogêneo Rodriguez e

Oliemans

(2006)

Rodriguez e

Baldani (2012) CFD

0,13 210,18 112,69 34,89 147,53

0,16 203,48 130,64 18,68 178,03

NETeF-EESC-USP 31

Resultados: Gradiente de Pressão - Erros

Erros médios globais para as previsões do gradiente de pressão -Castro (2013), horizontal.

Page 32: Estudo Teórico e Simulação Numérica do Escoamento ...

Malhas (em n° de nós): (a) 200.000 , (b) 800.000 , (c) 3.200.000 . (ICEM-CFD®).

Geometria e Parâmetros de Onda

32NETeF-EESC-USP

Uws = 0,12 m/s Amplitude Desvio

padrão

[mm]

Comprim. Desvio

padrão

[mm]

Celeridade Desvio

padrão

[mm] Uos = 0,06 m/s [mm] [mm] [m/s]

Experimental 3,02 1,31 25,72 9,85 0,17 0,03

CFD 0,99 0,62 94,91 39,06 - -