370 likes | 540 Views
Aula 9. análise das equações de conservação. sistema físico . modelo de camadas múltiplas. transporte em suspensão. L a. transporte por arrastamento. camada de mistura. substrato. h = h s + h b : profundidade do escoamento. h s : espessura da camada de transporte em suspensão.
E N D
Aula 9 análise das equações de conservação
sistema físico • modelo de camadas múltiplas transporte em suspensão La transporte por arrastamento camada de mistura substrato h = hs + hb: profundidade do escoamento hs : espessura da camada de transporte em suspensão hb : espessura da camada de transporte por arrastamento La : espessura da camada de mistura Yb : cota do fundo
modelo conceptual • granulometria uniforme, velocidade média, baixas concentrações, equilíbrio, regime permanente massa total conservação da massa de sedimentos quantidade de movimento da mistura equações de fecho: variáveis dependentes: : profundidade do escoamento : concentração de sedimentos : velocidade média do escoamento : declive da linha de energia : cota do fundo
C x análise das equações • motivação considere-se a equação de advecção pura de uma grandeza C: (A) seja, para simplificar o cálculo, e G constantes. em particular, seja: considere-se o problema de valores iniciais e de fronteira representado por (A) e pelas condições de fronteira e iniciais: comportamento esperado:
n+1 n i i+1 i1 análise das equações • motivação discretize-se a derivada espacial por diferenças “upwind” e a derivada temporal por diferenças de 1ª ordem: discretize-se o domínio de cálculo em 10 secções de cálculo (9 trechos) e proceda-se ao cálculo numérico considerando
n+1 n i i+1 i1 análise das equações • motivação discretize-se o domínio de cálculo em 10 secções de cálculo (9 trechos) e proceda-se ao cálculo numérico
n+1 n i i+1 i1 análise das equações • motivação análise: note-se que a derivada material de uma grandeza C se escreve, num referencial Eulereano: comparando com a equação (A) conclui-se que i.e., tem o significado físico de uma velocidade de propagação conclusão: existe uma velocidade física para a propagação de informação relativa a um fenómeno essencialmente advectivo e, também, uma velocidade numérica de propagação; num modelo numérico, terão que ser compatíveis!
n+1 n i i+1 i1 análise das equações • motivação direcção de propagação física: jusante para montante direcção de propagação numérica: montante para jusante se se colocar a condição de fronteira na secção correcta (jusante) e se discretizar o termo equação convectivo da equação (A) por diferenças finitas de 1ª ordem “downwind” obtém-se
análise das equações • objectivos da análise matemática dos modelos de transporte de sedimentos - determinar as velocidades (magnitude e sentido) de propagação de informação inerentes ao modelo conceptual - determinar a natureza da informação propagada • aplicações - determinação do número e natureza das condições de fronteira e iniciais - escolha dos esquemas numéricos em face da correcta propagação da informação no domínio de cálculo
t P t P x x análise das equações • exemplo, modelo #5 conservação da massa de sedimentos velocidade de propagação termos de fonte forma canónica não-conservativa de uma pde (equação diferencial parcial) derivada temporal gradiente velocidade de propagação: (perturbações na cota do fundo) grandeza transportada:
t J M I x análise das equações • exemplo, modelo #5 condições no contorno: t x 1 condição de fronteira a montante 1 condição inicial 1 condição de fronteira a jusante 1 condição inicial
análise das equações • problema unidimensional, caso geral: sistema de n pdes de 1ª ordem forma canónica não-conservativa de um sistema de n pdes V: vector das variáveis dependentes primitivas G: vector dos termos de fonte (irrelevante para as velocidades de propagação) A e B : por analogia com a equação diferencial única do exemplo anterior, as matrizes A e B dão conta das direcções e velocidades de propagação de informação física quanto à determinação das velocidades de propagação, os conceitos fundamentais são: linearidade/não linearidade hiperbolicidade
análise das equações • problema unidimensional, caso geral: sistema de n pdes de 1ª ordem linearidade/não linearidade os sistemas de equações de conservação associados a processos fluviais são sistemas de pdes de 1ª ordem quasi-lineares, i.e., as matrizes A e B são função das variáveis dependentes primitivas mas não das suas derivadas. exemplo:
análise das equações • hiperbolicidade do sistema (ver acetato para a noção de fase, velocidades de fase = características do sistema) definição: . o sistema é hiperbólico se tiver n direcções de propagação independentes . o sistema é hiperbólico se o polinómio característico (ex: ) tiver n raízes reais distintas . o sistema é hiperbólico se a matriz A-1B admitir n vectores próprios independentes nota: as características do sistema permitem conhecer a velocidade e a direcção de propagação da informação no domínio; falta conhecer as grandezas efectivamente propagadas!
análise das equações • hiperbolicidade do sistema exemplo: equações de Saint-Venant forma conservativa forma não-conservativa (TPC) polinómio característico
t t P P x x escoamento lento escoamento rápido análise das equações • hiperbolicidade do sistema exemplo: equações de Saint-Venant características do sistema: propagação da informação para jusante (cheias) propagação da informação para montante ou jusante (efeitos de regolfo)
análise das equações • hiperbolicidade do sistema exemplo: equações de Saint-Venant que informação é propagada ao longo das linhas características? ou... pode um sistema ser escrito na forma ? variáveis características: informação propagada ao longo das linhas características sim, pode, desde que, sendo S uma matriz de mudança de base, seja possível proceder à transformação nesse caso com
análise das equações • hiperbolicidade do sistema exemplo: equações de Saint-Venant a matriz de mudança de base é, por definição composta pelos vectores próprios de M. em rigor, as linhas de S são os vectores próprios esquerdos de M. para as equações de Saint-Venant, B = M e os vectores próprios são:
análise das equações • hiperbolicidade do sistema exemplo: equações de Saint-Venant
análise das equações • hiperbolicidade do sistema exemplo: equações de Saint-Venant informação propagada ao longo de informação propagada ao longo de
t P x análise das equações • hiperbolicidade do sistema exemplo: equações de Saint-Venant ao longo de ao longo de exemplo: escoamento lento
análise das equações • hiperbolicidade do sistema exemplo: modelo #4 – equilíbrio, granulometria uniforme, velocidade média massa total conservação da massa de sedimentos quantidade de movimento da mistura o polinómio característico, , é . tem três raízes reais distintas (três vectores próprios independentes). é portanto um sistema hiperbólico!
análise das equações • hiperbolicidade do sistema exemplo: modelo #4 – equilíbrio, granulometria uniforme, velocidade média linhas características, comparação com as equações de Saint-Venant para água limpa concentrações calculadas por: linha vermelha ( ) fórmula de Meyer-Peter & Muller; linha azul ( ) fórmula de Bagnold. notas: - (1) é, fundamentalmente, idêntica à água limpa; - (2) e (3) são afectadas pelo transporte de sedimentos; - se Fr < 0.7 o sistema exibe duas escalas distintas.
t t P P x análise das equações • hiperbolicidade do sistema exemplo: modelo #4 – equilíbrio, granulometria uniforme, velocidade média aspecto das linhas características (notar que as características não mudam de sinal com Fr... como se define o regime crítico?) Fr > 1.0 Fr < 0.7 x
análise das equações • estudo da hiperbolicidade; conclusões separação de escalas: se Fr < 0.7, (2) é aproximadamente igual a e (3) é aproximadamente igual a s (justifica o modelo #5) massa total massa de sedimentos quantidade de movimento da mistura
t t x x análise das equações • estudo da hiperbolicidade; conclusões separação de escalas: se Fr < 0.7, o modelo #5 representa uma boa aproximação. resolução desacoplada (I) do sistema de equações porque - propagação na fase líquida resolvida como uma sucessão de regimes permanentes; velocidades de propagação l(1) e l(2) infinitas (l(2) determina o sinal, consoante o número de Froude); - num dado Dt, a equação da fase líquida é resolvida antes da equação relativa à evolução morfológica (problemas: ver acetatos). P P
t x análise das equações • estudo da hiperbolicidade; conclusões separação de escalas: se Fr < 0.7, o modelo #5 representa uma boa aproximação. resolução desacoplada (II) do sistema de equações porque - num dado Dt, as equações (dinâmicas completas, ie. com termos de inércia local) da fase líquida são resolvidas antes das equações relativas à conservação da massa de sedimentos e do leito (problemas: ver acetatos). - localmente: M
t M J I x análise das equações • estudo da hiperbolicidade; conclusões condições iniciais e de fronteira; exemplo: modelos de transporte em equilíbrio: jusante: curva de vazão (Fr < 1) ou evolução temporal da cota do fundo (Fr > 1) CFM 1 CFM 2 montante: hidrograma de caudais sólidos e líquidos ou hidrograma de caudais sólidos e alturas do escoamento CFJ CIs iniciais: h, u e Yb
t I x análise das equações • estudo da hiperbolicidade; conclusões condições iniciais e de fronteira; exemplo: modelos de transporte em equilíbrio: montante: notar que a variável dependente é a cota do fundo, Yb, mas a condição de fronteira relativa a ls = l(2) é expressa em termos de Qs. há que converter, na vizinhança da fronteira, o caudal sólido em equilíbrio em cotas do fundo; M J CFM 1 CFM 2 CFJ em modelos desacoplados este procedimento pode levar ao mau-condicionamento do problema (oscilações que crescem a partir da fronteira). CIs
P J M t CFJ CFM 1 CFM 2 x análise das equações • estudo da hiperbolicidade; conclusões condições iniciais e de fronteira; exemplo: modelos de transporte em desequilíbrio (Cb é variável dependente): notas: i) Qs é facilmente introduzido na equação de conservação da massa de sedimentos na camada de transporte; ii) não se pode prescrever a cota do fundo nas fronteiras sob pena de provocar o mau condicionamento do problema. I CIs
análise das equações • estudo da hiperbolicidade; conclusões condições iniciais e de fronteira: - em geral, o número de condições independentes a especificar numa dada superfície de contorno é igual ao número de linhas características que “entram” por essa superfíce; - simbolicamente: em que C(k) é a expressão vectorial da linha caracterísitca cuja velocidade de fase é l(k).
análise das equações • estudo da hiperbolicidade; conclusões problemas descontínuos, soluções fracas. forma conservativa: forma não-conservativa: com: o aparecimento de soluções descontínuas é inevitável!
análise das equações • estudo da hiperbolicidade; conclusões problemas descontínuos, soluções fracas. t1+dt U+ caminho do choque dt U- t1 x1 x1+dx dx
análise das equações • estudo da hiperbolicidade; conclusões soluções fracas – solução descontínua num número contável de pontos. exemplo: solução de Stoker do o problema de Riemann para as shallow-water equations. teorema de Lax (1957): se i) o sistema de equações é estritamente hiperbólico, se os fluxos são funções contínuas e diferenciáveis e iii) se a amplitude da descontinuidade inicial é finita, então a solução do problema de Riemann consiste em n ondas (choques ou ondas de expansão), em que n é a dimensão da matriz jacobiana do sistema, mediados por n+1 estados constantes.
análise das equações • estudo da hiperbolicidade; conclusões soluções fracas – solução descontínua num nuúmero contável de pontos. exemplo: solução de Stoker do o problema de Riemann para as shallow-water equations. onda de expansão associada a l- t através do choque - condições de Rankine-Hugoniot: choque associado a l+ x através da onda de expansão – quasi-invariantes de Riemann:
análise das equações • estudo da hiperbolicidade; conclusões soluções fracas – solução descontínua num nuúmero contável de pontos. exemplo: solução de Stoker do o problema de Riemann para as shallow-water equations. solução para a onda de expansão: solução para o estado constante e para o choque: incógnitas: h*, u* e S
análise das equações • estudo da hiperbolicidade; conclusões soluções fracas – solução descontínua num nuúmero contável de pontos. exemplo: solução de Stoker do o problema de Riemann para as shallow-water equations.