Equação de Águas Rasas: mudanças entre as edições
| Linha 14: | Linha 14: | ||
<math>\dfrac{hu)^{t + \Delta t}_{i, j} - (hu)^{t}_{i, j}}{\Delta t} + \left [ \dfrac{(hu^2 + \cfrac{1}{2}gh^2)^t_{i+1, j} - (hu^2 + \cfrac{1}{2}gh^2)^t_{i-1, j}}{2 \Delta x} \right ] + \left [ \dfrac{(huv)^t_{i, j+1} - (huv)^t_{i, j-1}}{2 \Delta y} \right ] = -g h^{t}_{i, j} b_{x. i, j}</math> | <math>\dfrac{hu)^{t + \Delta t}_{i, j} - (hu)^{t}_{i, j}}{\Delta t} + \left [ \dfrac{(hu^2 + \cfrac{1}{2}gh^2)^t_{i+1, j} - (hu^2 + \cfrac{1}{2}gh^2)^t_{i-1, j}}{2 \Delta x} \right ] + \left [ \dfrac{(huv)^t_{i, j+1} - (huv)^t_{i, j-1}}{2 \Delta y} \right ] = -g h^{t}_{i, j} b_{x. i, j}</math> | ||
<math>\dfrac | <math>\dfrac{(hv)^{t + \Delta t}_{i, j} - (hv)^{t}_{i, j}}{\Delta t} + \left [ \dfrac{(huv)^t_{i+1, j} - (huv)^t_{i-1, j}}{2 \Delta x} \right ] + \left [ \dfrac{(hv^2 + 1/2 gh^2)^t_{i, j+1} - (hv^2 + 1/2 gh^2)^t_{i, j-1}}{2 \Delta y} \right ] = -g h^{t}_{i, j} b_{y. i, j} </math> | ||
Resolvendo pelo método de FTCS (para frente no tempo) e ajustando aos limites de estabilidade, temos como resultado: | Resolvendo pelo método de FTCS (para frente no tempo) e ajustando aos limites de estabilidade, temos como resultado: | ||
Edição das 22h22min de 7 de outubro de 2021
(EM EDIÇÃO) Grupo: Gabriel Schmökel, Julia Remus e Pedro Inocêncio Rodrigues Terra
Forma Conservativa
A partir das equações de conservação de momento e de massa, pode ser obtida as equações de águas rasas na forma conservativa. A forma conservativa da equação de águas rasas desconsidera a viscosidade do fluido e as tensões de cisalhamento aplicadas nele. Ela pode ser descrita como:
Falhou ao verificar gramática (erro de sintaxe): {\displaystyle Inserir fórmula aqui}
Ao aproximar por diferenças finitas obtemos o sistema de equações discretizadas a seguir.
![{\displaystyle {\dfrac {h_{i,j}^{t+\Delta t}-h_{i,j}^{t}}{\Delta t}}+\left[{\dfrac {(hu)_{i+1,j}^{t}-(hu)_{i-1,j}^{t}}{2\Delta x}}\right]+\left[{\dfrac {(hv)_{i,j+1}^{t}-(hv)_{i,j-1}^{t}}{2\Delta y}}\right]=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b833d4e6a0776ccf53f290b26e4cab8bef845ecf)
![{\displaystyle {\dfrac {hu)_{i,j}^{t+\Delta t}-(hu)_{i,j}^{t}}{\Delta t}}+\left[{\dfrac {(hu^{2}+{\cfrac {1}{2}}gh^{2})_{i+1,j}^{t}-(hu^{2}+{\cfrac {1}{2}}gh^{2})_{i-1,j}^{t}}{2\Delta x}}\right]+\left[{\dfrac {(huv)_{i,j+1}^{t}-(huv)_{i,j-1}^{t}}{2\Delta y}}\right]=-gh_{i,j}^{t}b_{x.i,j}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/520dd8da0bc59958cc15ca8033999e9705b7f49f)
![{\displaystyle {\dfrac {(hv)_{i,j}^{t+\Delta t}-(hv)_{i,j}^{t}}{\Delta t}}+\left[{\dfrac {(huv)_{i+1,j}^{t}-(huv)_{i-1,j}^{t}}{2\Delta x}}\right]+\left[{\dfrac {(hv^{2}+1/2gh^{2})_{i,j+1}^{t}-(hv^{2}+1/2gh^{2})_{i,j-1}^{t}}{2\Delta y}}\right]=-gh_{i,j}^{t}b_{y.i,j}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4030dec5c129aefa2aca422747c0456b63ce9844)
Resolvendo pelo método de FTCS (para frente no tempo) e ajustando aos limites de estabilidade, temos como resultado:
.... aqui gráfico ....
Para esse desenvolvimento encontramos algumas dificuldades para resolução do sistema de equações.