Processos de mistura em águas rasas
Storyboard
Os mecanismos de mistura em áreas rasas são gerados por diversos tipos de ondas. Entre eles estão as ondas internas, as ondas superficiais, a interação entre ondas e correntes, as marés e a quebra das ondas na costa.
ID:(1629, 0)
Mecanismos
Iframe
Mecanismos
ID:(15614, 0)
Número de Strouhal em função do número de Reynold
Imagem
O número de Strouhal ($St$)9480 está relacionado empiricamente com o número de Reynolds ($Re$)9107. O número de Strouhal ($St$)9480 está associado com la frequência de geração de vórtice ($\omega$)9496, la velocidade de fricção ($U_d$)9466 e la profundidade total ($H$)9465 é
$ St \equiv \displaystyle\frac{ \omega H }{ U_d }$ |
Isso permite estimar, através de o número de Reynolds ($Re$)9107, a frequência com que a concentração pode trocar os componentes a serem difundidos. No entanto, deve-se ter em mente que o processo pode ser interrompido se a frequência for menor do que a das marés.
ID:(12199, 0)
Estresse cinemático
Descrição
Se assumirmos que não há vento na superfície, pode-se supor que não há tensão na superfície. Portanto, haverá apenas a tensão da água no fundo. Esta tensão diminuirá linearmente do fundo para a superfície. Para simplificar a modelagem, pode-se usar a proporção entre la profundidade ($z$)9470 e la profundidade total ($H$)9465, o que nos dá um fator adimensional la profundidade relativa ($\xi$)10349. La estresse cinemático ($\tau_x$)10351 será, portanto, proporcional a
$\tau_x \propto 1-\xi$
Como la estresse cinemático ($\tau_x$)10351 é equivalente à densidade de energia dividida pela densidade, o valor no fundo deve ser proporcional ao quadrado da velocidade no fundo. Isso é descrito no modelo com la velocidade de fricção ($U_d$)9466 e significa que
$\tau_x \propto U_d^2$
Finalmente, há o efeito de la rugosidade ($k$)9463 do fundo do mar, ou seja, a proporção de o desigualdade ($d$)9464 e la profundidade total ($H$)9465. Isso significa que la estresse cinemático ($\tau_x$)10351 deve ser corrigido por um fator análogo à profundidade:
$\tau_x \propto \displaystyle\frac{1-\xi}{1-k}$
Assim, obtém-se um modelo da seguinte forma:
$ \tau_x \equiv \displaystyle\frac{ U_d ^2}{1- k }(1- \xi )$ |
que é representado graficamente a seguir:
ID:(15630, 0)
Comprimento de mistura
Descrição
La comprimento de mistura ($l$)9489 corresponde ao que poderia ser descrito como o tamanho dos vórtices. Próximo à parede, esses vórtices só podem ter um tamanho máximo igual à distância até a parede, que é mínima. À medida que nos aproximamos da superfície, esses vórtices podem se tornar cada vez maiores, então a função deve atingir um máximo nesse ponto.
Para simplificar a modelagem, pode-se usar a proporção entre la profundidade ($z$)9470 e la profundidade total ($H$)9465, o que nos proporciona um fator adimensional la profundidade relativa ($\xi$)10349. Assim, uma função simples que atende a essa descrição é:
$l \propto \xi\left(1-\displaystyle\frac{1}{2}\xi\right)$
Por outro lado, o modelo de camada limite de Prandtl mostra que estas são uma fração do fluxo com uma largura igual a la profundidade total ($H$)9465 e uma proporção de la constante de Karman ($\kappa$)9477, então:
$l \propto \kappa H$
Finalmente, devemos corrigir pelo efeito da rugosidade da mesma forma que para o estresse cinemático:
$l \propto \displaystyle\frac{\kappa H}{1-k}$
Portanto, la comprimento de mistura ($l$)9489 pode ser modelado da seguinte forma:
$ l \equiv \displaystyle\frac{ \kappa H }{1 - k } \xi \left(1 - \displaystyle\frac{1}{2} \xi \right)$ |
ID:(12201, 0)
Viscosidade do redemoinho
Top
Quando Prandtl modela a formação de turbilhões perto das paredes, ele estabelece a relação entre la viscosidade turbulenta ($A$)9469, la comprimento de mistura ($l$)9489 e o gradiente de o perfil da velocidade ($u_z$)9473 em la profundidade ($z$)9470 da seguinte forma:
$ A = l ^2\displaystyle\frac{\partial u_z }{\partial z }$ |
Por outro lado, a força viscosa típica, que é modelada como a viscosidade multiplicada pela superfície de contato e o gradiente de velocidade, corresponde a la estresse cinemático ($\tau_x$)10351 no caso das turbulências:
$ \tau_x = A \displaystyle\frac{\partial u_z }{\partial z }$ |
De ambas as equações, surge a relação:
$ \tau_x = \displaystyle\frac{ A ^2 }{ l ^2 }$ |
Essa relação permite calcular la viscosidade turbulenta ($A$)9469 em função de la estresse cinemático ($\tau_x$)10351 e la comprimento de mistura ($l$)9489, que são modelados neste caso. Assim, com la profundidade total ($H$)9465, la velocidade de fricção ($U_d$)9466, la rugosidade ($k$)9463, la profundidade relativa ($\xi$)10349 e la constante de Karman ($\kappa$)9477:
$ A = \displaystyle\frac{ \kappa H U_d }{(1- k )^{3/2}} \xi \left(1- \displaystyle\frac{1}{2} \xi \right)\sqrt{1- \xi }$ |
que é representado a seguir:
O resultado é que a viscosidade turbulenta é máxima na profundidade média e se reduz a valores mínimos tanto próximo ao fundo quanto próximo à superfície. Ou seja, nessas zonas a mistura e a perda de momento são menores.
ID:(15624, 0)
Perfil de velocidade
Top
Como la estresse cinemático ($\tau_x$)10351 é igual a la viscosidade turbulenta ($A$)9469 e ao gradiente de o perfil da velocidade ($u_z$)9473 em la profundidade ($z$)9470, a equação pode ser integrada para obter o perfil de velocidade:
$u_z = \displaystyle\int_d^z \frac{\tau_x}{A} , dz'$
Após integrar essa expressão, com la velocidade de fricção ($U_d$)9466, la constante de Karman ($\kappa$)9477, la rugosidade ($k$)9463 e la profundidade relativa ($\xi$)10349, obtemos:
$ u_z = \displaystyle\frac{ U_d }{ \kappa }\sqrt{1-k}\ln\left(\displaystyle\frac{ \xi }{ k }\right)$ |
que corresponde à famosa lei logarítmica desenvolvida por Prandtl e Schlichting.
O perfil é mostrado no seguinte gráfico:
O perfil também permite relacionar tanto ($$)9467 com la velocidade de fricção ($U_d$)9466 em função de la rugosidade ($k$)9463 e la constante de Karman ($\kappa$)9477, o que, por sua vez, permite definir um o arrastar para baixo ($C_D$)9468 com:
$ U ^2 = \displaystyle\frac{ U_d ^2}{ C_D }$ |
e
$ C_D = \displaystyle\frac{ \kappa ^2 }{(1- k ) \ln^2(1/ k )}$ |
ID:(15623, 0)
Concentração de sedimentos
Top
Se considerarmos o comportamento do material suspenso, observamos dois principais fatores. Primeiramente, há uma tendência de que ele se sedimente com uma velocidade la taxa de sedimentação ($\omega_s$)9476, gerando um fluxo que depende de ($$)9478, expresso como:
$\omega_s c_z$
Por outro lado, os redemoinhos tendem a misturar a água, gerando uma difusão que transporta os sedimentos em direção à superfície. Esse fluxo, representado por la viscosidade turbulenta ($A$)9469, é dado pelo gradiente de ($$)9478 em la profundidade ($z$)9470, igual a:
$A\displaystyle\frac{\partial c_z}{\partial z}$
A distribuição se forma quando os sedimentos alcançam o equilíbrio, onde o fluxo de sedimentação é igual à difusão gerada pelos redemoinhos em direção à superfície. Integrando ambos os termos da equação com o taxa de erosão ($E$)10350 e o desigualdade ($d$)9464, obtemos a distribuição:
$c_z=\displaystyle\frac{E}{\omega_s}e^{\displaystyle\int_d^z \omega_s/A dz'}$
Após empregar a expressão obtida para la viscosidade turbulenta ($A$)9469 com o fator de Rouse ($R_s$)10352, la rugosidade ($k$)9463 e la profundidade relativa ($\xi$)10349, derivamos a expressão:
$ c_z = \displaystyle\frac{ E }{ \omega_s }\left(\displaystyle\frac{ k }{ \xi }\right)^{ R_s }$ |
que pode ser representada graficamente como:
ID:(15631, 0)
Modelo
Top
Parâmetros
Variáveis
Cálculos
Cálculos
Cálculos
Equações
$ A = \displaystyle\frac{ \kappa H U_d }{(1- k )^{3/2}} \xi \left(1- \displaystyle\frac{1}{2} \xi \right)\sqrt{1- \xi }$
A = kappa * U_d * H * xi * (1- xi /2)*sqrt(1- xi )/(1- k )^(3/2)
$ A = l ^2\displaystyle\frac{\partial u_z }{\partial z }$
A = l ^2*@DIF( u_z , z )
$ C_D = \displaystyle\frac{ \kappa ^2 }{(1- k ) \ln^2(1/ k )}$
C_D = kappa ^2/((1- k )* log(1/ k )^2)
$ c_z = \displaystyle\frac{ E }{ \omega_s }\left(\displaystyle\frac{ k }{ \xi }\right)^{ R_s }$
c_z =( E / omega_s )*( k / xi )^ R_s
$ k \equiv \displaystyle\frac{ d }{ H }$
k = d / H
$ l \equiv \displaystyle\frac{ \kappa H }{1 - k } \xi \left(1 - \displaystyle\frac{1}{2} \xi \right)$
l = kappa * H * xi *(1 - xi /2)/(1 - k )
$ R_0 \equiv \displaystyle\frac{ \omega_s }{ \kappa U_d }$
R_0 = omega_s /( kappa * U_d )
$ R_s \equiv R_0 (1 - k )^{3/2}$
R_s = R_0 (1 - k )^(3/2)
$ St \equiv \displaystyle\frac{ \omega H }{ U_d }$
St = omega * H / U_d
$ \tau_x = A \displaystyle\frac{\partial u_z }{\partial z }$
tau_x = A *@DIF( u_z , z )
$ \tau_x = \displaystyle\frac{ A ^2 }{ l ^2 }$
tau_x = A ^2 / l ^2
$ \tau_x \equiv \displaystyle\frac{ U_d ^2}{1- k }(1- \xi )$
tau_x = U_d ^2 *(1- xi )/(1- k )
$ U ^2 = \displaystyle\frac{ U_d ^2}{ C_D }$
U ^2 = U_d ^2/ C_D
$ u_z = \displaystyle\frac{ U_d }{ \kappa }\sqrt{1-k}\ln\left(\displaystyle\frac{ \xi }{ k }\right)$
u_z = U_d sqrt(1- k )*log( xi / k ))/ kappa
$ \xi \equiv \displaystyle\frac{ z }{ H }$
xi = z / H
ID:(15618, 0)
Número de Strouhal
Equação
O número de Strouhal ($St$)9480 caracteriza la frequência de geração de vórtice ($\omega$)9496. Compara a velocidade associada a la frequência de geração de vórtice ($\omega$)9496 e seu tamanho com o do fluxo representado por la profundidade total ($H$)9465.
Portanto, com isso, temos
$ St \equiv \displaystyle\frac{ \omega H }{ U_d }$ |
ID:(12198, 0)
Rugosidade inferior
Equação
O comportamento da corrente e a turbulência a ser gerada ou amortecida dependem de la rugosidade ($k$)9463 do leito marinho. Isso é definido comparando o perfil médio de o desigualdade ($d$)9464 com o perfil de la profundidade total ($H$)9465 onde ele está localizado.
Portanto, define-se que la rugosidade ($k$)9463 é
$ k \equiv \displaystyle\frac{ d }{ H }$ |
ID:(12183, 0)
Profundidade relativa
Equação
La profundidade relativa ($\xi$)10349 é definido em termos de la profundidade ($z$)9470 e la profundidade total ($H$)9465, expresso da seguinte forma:
$ \xi \equiv \displaystyle\frac{ z }{ H }$ |
ID:(12182, 0)
Estresse cinemático
Equação
Para um fluxo laminar, la força viscosa ($F_v$)4979 pode ser calculado a partir de os superfícies paralelas ($S$)10119, la viscosidade ($\eta$)5422, la diferença de velocidade entre superfícies ($\Delta v$)5556 e la distância entre superfícies ($\Delta z$)5436 utilizando a seguinte equação:
$ F_v =- S \eta \displaystyle\frac{ \Delta v }{ \Delta z }$ |
No caso de fluxo turbulento, pode-se estabelecer uma analogia definindo uma viscosidade turbulenta ($A$)9469,1 como a viscosidade dividida pela densidade, associada à força por área e densidade, o que denominaremos la estresse cinemático ($\tau_x$)10351, calculado em função de o perfil da velocidade ($u_z$)9473 e la profundidade ($z$)9470 da seguinte forma:
$ \tau_x = A \displaystyle\frac{\partial u_z }{\partial z }$ |
ID:(12191, 0)
Lei de Prandtl
Equação
Em 1925, Prandtl introduziu o conceito de uma camada limite na qual os turbilhões misturam o fluido e transferem momento de forma semelhante à modelagem da transferência em nível molecular, gerando comportamento viscoso. O tamanho desta zona é definido como la comprimento de mistura ($l$)9489 e o efeito é descrito com um análogo à viscosidade que corresponde a la viscosidade turbulenta ($A$)9469. Isso pode ser estimado com o gradiente de o perfil da velocidade ($u_z$)9473 em la profundidade ($z$)9470 utilizando:
$ A = l ^2\displaystyle\frac{\partial u_z }{\partial z }$ |
[1] Prandtl, Ludwig (1925). "Bericht über Untersuchungen zur ausgebildeten Turbulenz" (Relatório sobre investigações em turbulência desenvolvida). Z. Angew. Math. Mech. 5 (2): 136.
ID:(12186, 0)
Relação de tensão cinemática
Equação
La estresse cinemático ($\tau_x$)10351 pode ser calculado a partir de la viscosidade turbulenta ($A$)9469 e la comprimento de mistura ($l$)9489 usando o seguinte método:
$ \tau_x = \displaystyle\frac{ A ^2 }{ l ^2 }$ |
Assim como la viscosidade turbulenta ($A$)9469 se relaciona com la comprimento de mistura ($l$)9489, o perfil da velocidade ($u_z$)9473 e la profundidade ($z$)9470 são definidos como
$ A = l ^2\displaystyle\frac{\partial u_z }{\partial z }$ |
e, dado que la estresse cinemático ($\tau_x$)10351 é
$ \tau_x = A \displaystyle\frac{\partial u_z }{\partial z }$ |
eliminar o gradiente resulta em
$ \tau_x = \displaystyle\frac{ A ^2 }{ l ^2 }$ |
ID:(15633, 0)
Difusividade de vórtice oceânico mais profundo
Equação
La estresse cinemático ($\tau_x$)10351 será máximo próximo ao fundo do oceano e nulo na superfície, desde que não haja vento na superfície do oceano. Como está associado a la velocidade de fricção ($U_d$)9466 no fundo, mas precisa ser corrigido pelo efeito de la rugosidade ($k$)9463, pode ser modelado com base em la profundidade relativa ($\xi$)10349 da seguinte forma:
$ \tau_x \equiv \displaystyle\frac{ U_d ^2}{1- k }(1- \xi )$ |
ID:(12202, 0)
Comprimento de mistura
Equação
A zona de mistura introduzida por Prandtl, de tamanho la comprimento de mistura ($l$)9489, é estimada como uma fração da ordem de la constante de Karman ($\kappa$)9477 de la profundidade total ($H$)9465. Além disso, deve-se ajustar pelo efeito de la rugosidade ($k$)9463, e considerar que la comprimento de mistura ($l$)9489 depende de la profundidade relativa ($\xi$)10349, sendo nulo no fundo e aproximadamente constante e máximo perto da superfície. Portanto, pode ser modelado da seguinte forma:
$ l \equiv \displaystyle\frac{ \kappa H }{1 - k } \xi \left(1 - \displaystyle\frac{1}{2} \xi \right)$ |
[1] Prandtl, Ludwig (1925). "Bericht über Untersuchungen zur ausgebildeten Turbulenz" (Relatório sobre investigações em turbulência desenvolvida). Z. Angew. Math. Mech. 5 (2): 136.
ID:(12194, 0)
Viscosidade do redemoinho
Equação
A partir de medições, podemos modelar la viscosidade turbulenta ($A$)9469 com la profundidade relativa ($\xi$)10349, la profundidade total ($H$)9465, la rugosidade ($k$)9463, la velocidade de fricção ($U_d$)9466 e la constante de Karman ($\kappa$)9477 usando a expressão:
$ A = \displaystyle\frac{ \kappa H U_d }{(1- k )^{3/2}} \xi \left(1- \displaystyle\frac{1}{2} \xi \right)\sqrt{1- \xi }$ |
Assim como la estresse cinemático ($\tau_x$)10351 se relaciona com la viscosidade turbulenta ($A$)9469 e la comprimento de mistura ($l$)9489, tem-se que:
$ \tau_x = \displaystyle\frac{ A ^2 }{ l ^2 }$ |
Se forem usados la constante de Karman ($\kappa$)9477, la profundidade total ($H$)9465 e la rugosidade ($k$)9463:
$ l \equiv \displaystyle\frac{ \kappa H }{1 - k } \xi \left(1 - \displaystyle\frac{1}{2} \xi \right)$ |
e com la velocidade de fricção ($U_d$)9466:
$ \tau_x \equiv \displaystyle\frac{ U_d ^2}{1- k }(1- \xi )$ |
obtém-se:
$ A = \displaystyle\frac{ \kappa H U_d }{(1- k )^{3/2}} \xi \left(1- \displaystyle\frac{1}{2} \xi \right)\sqrt{1- \xi }$ |
ID:(12185, 0)
Perfil de velocidade
Equação
O perfil da velocidade ($u_z$)9473 é uma função de la profundidade relativa ($\xi$)10349 e dos parâmetros la rugosidade ($k$)9463, la velocidade de fricção ($U_d$)9466 e la constante de Karman ($\kappa$)9477, representada da seguinte forma:
$ u_z = \displaystyle\frac{ U_d }{ \kappa }\sqrt{1-k}\ln\left(\displaystyle\frac{ \xi }{ k }\right)$ |
Assim como la estresse cinemático ($\tau_x$)10351 se relaciona com la viscosidade turbulenta ($A$)9469, o perfil da velocidade ($u_z$)9473 e la profundidade ($z$)9470, ele é definido por
$ \tau_x = A \displaystyle\frac{\partial u_z }{\partial z }$ |
e pode ser integrado de o desigualdade ($d$)9464 a la profundidade ($z$)9470 para obter a velocidade usando a seguinte expressão:
$u_z=\displaystyle\int_d^z\displaystyle\frac{\tau_x}{A}dz'$
Com a formulação de la viscosidade turbulenta ($A$)9469 em termos de la profundidade relativa ($\xi$)10349 juntamente com la profundidade total ($H$)9465, la rugosidade ($k$)9463 e la velocidade de fricção ($U_d$)9466, e considerando que
a seguinte equação para velocidade é derivada:
$u_z=\displaystyle\frac{U_d\sqrt{1-k}}{\kappa}(\ln(z/d) + \Phi(\xi,k))$
onde
$\Phi=2[\arctan(\lambda)-\arctan(\lambda_0)]-\ln\left(\displaystyle\frac{1+\lambda}{1+\lambda_0}\right)$
é definido com
$\lambda=\sqrt{1-\xi}$
e
$\lambda=\sqrt{1-k}$
Visto que em grande parte da profundidade
$\ln(z/d) \gg \Phi(\xi,k)$
o perfil de velocidade pode ser simplificado para
$ u_z = \displaystyle\frac{ U_d }{ \kappa }\sqrt{1-k}\ln\left(\displaystyle\frac{ \xi }{ k }\right)$ |
ID:(12187, 0)
Coeficiente de amortecimento
Equação
($$)9467 é proporcional a la velocidade de fricção ($U_d$)9466, com uma constante de proporcionalidade que depende de la constante de Karman ($\kappa$)9477 e la rugosidade ($k$)9463, conforme:
$ C_D = \displaystyle\frac{ \kappa ^2 }{(1- k ) \ln^2(1/ k )}$ |
ID:(12184, 0)
Velocidade de superfície
Equação
La velocidade de fricção ($U_d$)9466 é proporcional a ($$)9467, sendo a constante de proporcionalidade o arrastar para baixo ($C_D$)9468, que representa a relação entre as respectivas energias cinéticas:
$ U ^2 = \displaystyle\frac{ U_d ^2}{ C_D }$ |
ID:(12188, 0)
Difusão de uma concentração de sedimentos
Equação
($$)9478 é uma função de la profundidade relativa ($\xi$)10349 que depende de o taxa de erosão ($E$)10350, la taxa de sedimentação ($\omega_s$)9476, la rugosidade ($k$)9463 e o fator de Rouse ($R_s$)10352, e é calculada da seguinte forma:
$ c_z = \displaystyle\frac{ E }{ \omega_s }\left(\displaystyle\frac{ k }{ \xi }\right)^{ R_s }$ |
Os sedimentos tendem a cair para o fundo com uma taxa de sedimentação ($\omega_s$)9476,1, enquanto a difusão, que neste caso corresponde à mistura gerada por turbilhões, induz um fluxo igual a la viscosidade turbulenta ($A$)9469 e o gradiente de ($$)9478 em la profundidade ($z$)9470, da seguinte forma:
$A\displaystyle\frac{\partial c_z}{\partial z}+\omega_s c_z= 0$
Integrando essa expressão, obtemos:
$c_z = \displaystyle\frac{E}{\omega_s}e^{-\displaystyle\int_d^z \omega_s/A dz'}$
com la comprimento de mistura ($l$)9489:
$ \tau_x = \displaystyle\frac{ A ^2 }{ l ^2 }$ |
temos:
$c_z=\displaystyle\frac{E}{\omega_s}e^{-\displaystyle\int_d^z \omega_s/l\sqrt{\tau_x} dz'}$
o que resulta em:
$c_z=\displaystyle\frac{E}{\omega_s}\left(\displaystyle\frac{z}{d}\right)^{R_s}\Phi_c(\xi,k)$
com o fator de Rouse ($R_s$)10352 e o número de despertar ($R_0$)9475:
$ R_s \equiv R_0 (1 - k )^{3/2}$ |
onde:
$\Phi=\left(\displaystyle\frac{1+\lambda}{1+\lambda_0}\right)^{2R_s}e^{2R_s[\arctan(\lambda)-\arctan(\lambda_0)]}$
com:
$\lambda=\sqrt{1-\xi}$
e:
$\lambda=\sqrt{1-k}$
Como em grande parte da profundidade:
$\Phi\sim 1$
temos a distribuição de concentração:
$ c_z = \displaystyle\frac{ E }{ \omega_s }\left(\displaystyle\frac{ k }{ \xi }\right)^{ R_s }$ |
ID:(12193, 0)
Número de Rousse
Equação
O valor o número de despertar ($R_0$)9475 compara a velocidade de sedimentação, que compete com a difusão associada à corrente no fundo. Quando combinado com la taxa de sedimentação ($\omega_s$)9476, la velocidade de fricção ($U_d$)9466 e la constante de Karman ($\kappa$)9477, resulta em:
$ R_0 \equiv \displaystyle\frac{ \omega_s }{ \kappa U_d }$ |
ID:(12195, 0)
O fator Rouse
Equação
O valor o número de despertar ($R_0$)9475 avalia a velocidade de sedimentação, que compete com a difusão associada à corrente no fundo. Combinado com la taxa de sedimentação ($\omega_s$)9476, la velocidade de fricção ($U_d$)9466 e la constante de Karman ($\kappa$)9477, resulta em:
$ R_0 \equiv \displaystyle\frac{ \omega_s }{ \kappa U_d }$ |
Nos casos em que o fundo do mar não é liso, existe la rugosidade ($k$)9463, o que leva a uma correção do número de Rouse, que chamamos de fator de Rouse:
$ R_s \equiv R_0 (1 - k )^{3/2}$ |
ID:(15632, 0)