15
Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa [email protected] Taciana Toledo de Almeida Albuquerque [email protected] Universidade de São Paulo stituto de Astronomia Geofísica e Ciências Atmosfér Departamento de Meteorologia

Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa [email protected]

Embed Size (px)

Citation preview

Page 1: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Descrição do comportamento da camada de mistura através de um

modelo numérico de fechamento de primeira ordem.

América Murguía [email protected]

Taciana Toledo de Almeida [email protected]

Universidade de São PauloInstituto de Astronomia Geofísica e Ciências Atmosféricas

Departamento de Meteorologia

Page 2: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

ResumoIntroduçãoMetodologiaDadosResultadosConclusões

Page 3: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

O que é Camada Limite Planetária?

.Introdução

-É a região da atmosfera onde os fenômenos meteorológicos de pequena escala são importantes;

-Está sob a influência direta da superfície e responde as forçantes superficiais numa escala de tempo menor que 1 hora.

Estrutura Espacial da CLPDiurno

Noturno

Page 4: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Forçantes Superficiais:- Aquecimento Solar;

- Resfriamento Radiativo;

- Vento Horizontal.

De acordo com a importância dessas forçantes a CLP pode ser classificada como:

-Convectiva: aquecimento solar da superfície é suficiente para manter a convecção térmica;

-Neutra: quando nem o aquecimento solar e nem o resfriamento radiativo são suficientes para alterar as características da turbulência de origem mecânica.

-Estável: o resfriamento radiativo da superfície é suficiente para manter uma estratificação térmica através da qual a turbulência terá que realizar trabalho (Produção Mecânica);

Page 5: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Fluxo de Calor Sensível (H):

- É proporcional a covariância da velocidade vertical (w) e da temperatura potencial do ar ();

wcH P

- Radiação solar aquece a superfície

Energia é transferida para

atmosfera através do fluxo

turbulento de calor sensível.

>0 – existe transferência da superfície para atmosfera w w <0 – existe transferência de calor da atmosfera para superfície

Em Condições Convectivas a CLP se forma:

- Aquecimento da superfície pela radiação solar;

- Transferência da energia disponível na superfície para a atmosfera através dos fluxos turbulentos de calor sensível;

Page 6: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Caracterizar alguns parâmetros que influenciam a estrutura da CLC (6-18h), dentro da Camada de Mistura:

- Temperatura Potencial;

- Fluxo de Calor Sensível;

- Altura da Camada de Mistura;

- Intensidade da Inversão de Temperatura no topo da CM.

OBJETIVO

Camada de Mistura

- Região onde o gradiente vertical da temperatura potencial médio vale zero, devido a mistura turbulenta intensa;

- A mistura pode ser de origem térmica (convecção térmica), ou de origem mecânica, associada ao cisalhamento vertical do vento.

Page 7: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Descrição da Turbulência

- Não existe nenhuma descrição determinística dos escoamentos turbulentos, pois não existe até hoje solução das equações do movimento que descrevam um escoamento turbulento;- Os escoamentos turbulentos são descritos do ponto de vista estatístico:

iii uuu Flutuação estatística em torno da média

Comportamento médio

- Problema de fechamento: número de incógnitas é maior do que o de equações.

- Aplica-se então as propriedades da média de Reynolds

1) Fechamento de 1ª ordem “bulk” Teoria Fluxo Gradiente (K)

Teoria do Comprimento de Mistura

Parametrização

- O surgimento das covariâncias conduz a um problema de fechamento de 1ª ordem;

Page 8: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Metodologia

Aspectos do Modelo

1) Aspectos dinâmicos: Equação da Termodinâmica

2) Aspectos numéricos: Esquema Avançado no Tempo (n, n+1)

Utiliza o método de Diferenças Finitas

3) Aspectos físicos:-Condições dentro e sobre a CLC é considerada horizontalmente homogênea (advecção horizontal não é considerada);

- Céu limpo: convecção livre --- PT >>PM

Page 9: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

z

w

t

Descrição do Modelo

Eq. da Termodinâmica

h

ww

tiM 0

Integrando com à altura Temperatura Potencia na CM

iw

t

h ´)´(Altura da CLP

h

ww

t

h

ti 0)''()''(

Intensidade da inversão de temperatura no topo da CM

Page 10: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Parâmetros Caso1 Caso2 Caso3 Dia 23 23 23 Mês 7 7 7 Latitude (º) -20,5 -20,5 -20,5 Temperatura da superfície do solo (K)

294 300 305

Temperatura potencial da CM (K) 294 300 300

Intensidade da inversão de altitude (o K) 2 5 2

Gradiente vertical de temperatura potencial na atmosfera livre (o K/m)

0,0075 0,003 0,0075

Altura inicial da CM (m) 500 100 100

Passo de tempo de integração (s) 25 25 25 Taxa de entranhamento 0,1 0,1 0,1 Velocidade vertical sinótica (m/s) -0,005 -0,005 -0,005 Parâmetro de rugosidade (m) 0,02 0,02 0,1 Velocidade média na camada de mistura (m/s) 5,01 5,01 5,01 Pressão de vapor (mb) 12 12 12

Dados

Page 11: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Caso 1

293

294

295

296

297

298

299

300

4 6 8 10 12 14 16 18 20

Tempo Local (h)

0,0

0,5

1,0

1,5

2,0

2,5

4 6 8 10 12 14 16 18 20

Tempo Local (h)

400

500

600

700

800

900

1000

4 6 8 10 12 14 16 18 20

Tempo Local (h)

-50

0

50

100

150

200

4 6 8 10 12 14 16 18 20

Tempo Local (h)

w

Resultados

Page 12: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Caso 2

298300

302304

306308

310312

4 6 8 10 12 14 16 18 20

Tempo Local (h)

0

1

2

3

4

5

6

4 6 8 10 12 14 16 18 20

Tempo Local (h)

0

500

1000

1500

2000

2500

4 6 8 10 12 14 16 18 20

Tempo Local (h)

-50

50

150

250

350

450

4 6 8 10 12 14 16 18 20

Tempo Local (h)

w

Resultados

Page 13: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Caso 3

299

300

301

302

303

304

305

4 6 8 10 12 14 16 18 20

Tempo Local (h)

0,0

0,5

1,0

1,5

2,0

2,5

4 6 8 10 12 14 16 18 20

Tempo Local (h)

0

100

200

300

400

500

2 4 6 8 10 12 14 16 18 20

Tempo Local (h)

-10

0

10

20

30

40

50

4 6 8 10 12 14 16 18 20

Tempo Local (h)

w

Resultados

Page 14: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Conclusões

Boa representação do comportamento da CM no período convectivo.O fluxo turbulento de calor, nos três casos, mostrou que a entrada de ar mais quente numa camada mais fria fez que o ar comece a resfriar e com isso aumenta a altura da CM.Estes resultados podem ser comparados com o trabalho de Tennekes (1973), onde as distribuições são similares no comportamento da altura e a intensidade da inversão no topo.O modelo possui limitações, por ser simples para fazer simulações muito complexas, restringindo a utilização dele para ser aplicado em condições atmosféricas diferentes.

Page 15: Descrição do comportamento da camada de mistura através de um modelo numérico de fechamento de primeira ordem. América Murguía Espinosa americam@model.iag.usp.br

Obrigada!