Resumos
Neste artigo discutimos a propagação de ondas de pequena amplitude num líquido de profundidade constante com base em uma equação de ondas. A solução geral mostra que a velocidade de propagação é dependente do comprimento de onda de modo que ocorre o fenômeno da dispersão. Soluções para os casos de pequena e grande profundidades são também obtidas.
ondas de gravidade; dispersão de ondas; ondas num líquido
In this article we discuss the propagation of small amplitude waves in a liquid of constant depth based on a wave equation. The general solution shows that propagation speed depends on the wavelength and consequently there is waves dispersion. Solutions for shallow and the deep water are shown.
gravity waves; wave dispersion; waves on the liquid
ARTIGOS GERAIS
Ondas superficiais de gravidade
Surface gravity waves
Alcione S. Fernandes1 1 E-mail: alcionesf@ig.com.br ; Giselle M. Alves
Departamento de Física, Universidade Federal do Paraná, Curitiba, PR, Brasil Escola Técnica da Universidade Federal do Paraná, Curitiba, PR, Brasil
RESUMO
Neste artigo discutimos a propagação de ondas de pequena amplitude num líquido de profundidade constante com base em uma equação de ondas. A solução geral mostra que a velocidade de propagação é dependente do comprimento de onda de modo que ocorre o fenômeno da dispersão. Soluções para os casos de pequena e grande profundidades são também obtidas.
Palavras-chave: ondas de gravidade, dispersão de ondas, ondas num líquido.
ABSTRACT
In this article we discuss the propagation of small amplitude waves in a liquid of constant depth based on a wave equation. The general solution shows that propagation speed depends on the wavelength and consequently there is waves dispersion. Solutions for shallow and the deep water are shown.
Keywords: gravity waves, wave dispersion, waves on the liquid.
1. Introdução
O estudo da propagação de ondas é fundamental para a formação dos alunos dos cursos de física, matemática e das engenharias. No seu desenvolvimento é possivel aplicar varias técnicas matemáticas como o emprego de uma teoria linearizada, ferramenta simples e extremamente útil na resolução de problemas complexos, bem como a obtenção da solução de equações diferenciais com condições de contorno.
Neste artigo desenvolvemos, com base em conceitos acessíveis aos alunos dos cursos iniciais de cálculo e física, alguns aspectos da teoria da propagação de perturbações harmônicas na superfície de um líquido de profundidade constante e inicialmente em repouso. Uma abordagem computacional deste tema foi desenvolvida recentemente [1]. Informações descritivas [2] ou mais avançadas [3, 4] são encontradas na literatura.
Inicialmente, consideremos o escoamento de um fluido não viscoso e sujeito a uma forca gravitacional. O movimento deste fluido ideal é descrito pela equação de Euler [5]
e pela equação da continuidade
Nestas equações 
(r, t) representa a densidade do fluido, v(r, t) a velocidade, p(r, t) a pressão e f a força gravitacional específica (forca por unidade de massa).
Suponhamos que o fluido, inicialmente em equilíbrio ( = o, v = 0 e p = p0), experimenta uma perturbação de modo a adquirir a velocidade de escoamento 
(r, t), com densidade e pressão dadas por
com 
.
Vamos considerar que as perturbações das grandezas de equilíbrio e seus gradientes são pequenos de modo que seus produtos podem ser desconsiderados e somente os termos lineares são mantidos. Além disso, da Eq. (1) ∇po = 0f no equilíbrio. Substituindo as Eqs. (3) e (4) nas Eqs. (1) e (2), obtemos
onde admitimos que a densidade de equilíbrio é aproximadamente uniforme de modo que pudemos descartar o termo de segunda ordem · ∇0.
Se tomarmos o rotacional de ambos os membros da Eq. (5) e tendo em conta que a força gravitacional é conservativa, obtemos
Como ∇ x é proporcional à velocidade angular em cada ponto do fluido e esta não muda com o tempo, então a solução possível para a Eq. (7) é
de modo que o fluido é irrotacional em todos os pontos. Desta forma, podemos relacionar o campo de velocidades a um potencial escalar Φ(r, t), através de
com componentes
Se considerarmos que o fluido é incompressível, isto é, a densidade permanece constante, obtemos da Eq. (6) que
Substituindo a Eq. (9) na anterior, resulta a equação de Laplace ∇2Φ = 0 ou, em coordenadas cartesianas
Portanto, um fluido irrotacional e incompressível pode ser descrito por uma função potencial escalar que e solução da equação de Laplace e com componentes da velocidade dadas pelas Eqs. (10).
2. Ondas superficiais de gravidade
Consideremos em seguida, a propagação de uma perturbação provocada na superfície livre de um líquido com profundidade uniforme e sujeito a um campo gravitacional uniforme. Estas ondas progressivas longitudinais são denominadas ondas superficiais de gravidade [3]. As amplitudes perturbativas são consideradas pequenas em relação às grandezas de equilíbrio, de modo que podemos reter apenas os termos lineares. Uma situação adequada a estas condições é a de uma brisa suave ou a queda de uma pequena pedra na superfície de mar calmo.
Para estudar as ondas superficiais de gravidade num líquido de profundidade uniforme h, vamos considerar que:
(a) o líquido é perturbado a partir do repouso, por exemplo, pela ação de um vento suave; assim as velocidades e os deslocamentos das partículas são pequenos o que nos permite desconsiderar os seus produtos;
(b) o fluxo é irrotacional de modo que o campo de velocidades do líquido deriva de um potencial escalar;
(c) a amplitude ςo das ondas é muito menor que o comprimento de onda (ς0/λ « 1);
(d) não são considerados os efeitos de tensão superficial, viscosidade e condutividade térmica.
Para ondas progressivas na direção x do plano xz (veja Fig. 1) não há dependência da coordenada y de modo que o nosso problema é bidimensional. Desta forma, temos da Eq. (12) que
com as componentes das velocidades das partículas do líquido dadas por
O primeiro passo é obter o potencial escalar Φ(x, t) da Eq. (13) e que está sujeito às seguintes condições de contorno:
(a) no fundo a componente vertical da velocidade do líquido é nula, ou seja
(b) A superfície livre do líquido é descrita por z(x,t) = ς(x, t) e a velocidade para as partículas da superfície é dada por
Se expandirmos o primeiro membro desta equação em série de Taylor em torno de z = 0, concluimos que
onde a aproximação linear foi usada. Esta é uma condição de contorno cinemática pois envolve a velocidade da superfície livre do líquido.
(c) Como a tensão superficial não é levada em conta então as pressões em ambos os lados da superfície livre do líquido são iguais; além disso, como a pressão atmosférica é constante ao longo da superfície livre podemos considerá-la nula sem perda de generalidade, ou seja
e a equação de Bernoulli (A.8) torna-se
Em resumo, o problema consiste em determinar o potencial escalar Φ (x, z, t) para a equação de Laplace
sujeito às condições de contorno
Após obtermos a função Φ(x, z, t) podemos determinar as componentes da velocidade do líquido através das Eqs. (14). Para isso, entretanto, é necessário explicitar a forma para a superfície livre ς(x, t).
3. A função Φ (x, z, t) para perturbações periódicas
Consideremos as perturbações periódicas da superfície livre do líquido na forma de uma onda progressiva na direção x
onde ςo é a amplitude constante da onda. Esta escolha pode ser justificada pelo fato de que uma perturbação arbitrária pode ser construída adicionando componentes sinusoidais através de uma série de Fourier.
Comparando a equação anterior com as Eqs. (22) e (23) verificamos que o potencial escalar deve ser da forma
onde a amplitude Φ(z) e a relacão de dispersão ω = ω(k) devem ser determinadas.
Substituindo a Eq. (25) na (20), obtemos a equação diferencial linear de segunda ordem para a amplitude
com solução geral
onde a e b são constantes de integração. Da Eq. (25) temos
e utilizando as condições de contorno (21) e (22) resulta
Portanto, o potencial escalar (25) é dado por
Após determinarmos o potencial escalar podemos passar para a obtenção das demais variáveis que caracterizam o movimento do líquido como velocidades e posições das partículas, freqüência angular e velocidade de propagação da onda.
4. Resultados
4.1. Velocidades e trajetórias das partículas, freqüência angular e velocidade de fase da onda
As componentes da velocidade das partículas são dadas, através das Eqs. (14) e (30), por
Consideremos que as partículas de fluido oscilam com pequenas amplitudes em torno da posição de equilíbrio (xo, zo) e que diferem pouco das coordenadas (x, z); assim, podemos integrar no tempo as equações (31) para obtermos
de modo que as trajetórias das partículas são elipses descritas por
com semi-eixos maior e menor que dependem da profundidade de acordo com
Da condição de contorno (23) e Eq. (24) obtemos a freqüência angular ω como função de k, denominada relação de dispersão, como
Para a velocidade de fase ω/k, temos
onde utilizamos k = 2π/λ. A Eq. (36) esta representada pela curva (a) da Fig. 2.
O movimento ondulatório apresenta, então, as seguintes características:
(a) a Eq. (36) representa uma relação dispersiva pois a cada comprimento de onda corresponde uma velocidade de propagação diferente, de modo que as ondas se dispersam no líquido;
(b) na curva (a) da Fig. 2 observamos que, para uma dada profundidade h, a velocidade de fase cresce rapidamente com o comprimento de onda λ e tende para o valor limite (gh)1/2, como será mostrado na seção 5.2;
(c) as Eqs. (32) mostram que as fases dos movimentos (argumentos das funções sinusoidais) independem da profundidade z0; assim, as partículas executam movimentos oscilatórios em fase, isto é, em cada instante, se para um dado x0 uma delas esta no topo da trajetória então todas estarão;
(d) as Eqs. (34) indicam que os semi-eixos das trajetórias elípticas diminuem com a profundidade e, no fundo (z0 = -h), o movimento das partículas é horizontal pois b(z0) = 0.
Ao desprezarmos o termo ∇ consideramos que 
na equação de movimento (1); com as Eqs. (31) verificamos então, que os nossos resultados são válidos para
isto é, para situações em que a amplitude da perturbação é muito menor que o comprimento da onda que propaga, como foi observado no item (c) da seção 2.
4.2. Variações de pressão
A pressão hidrostatica na profundidade -z2
2
As coordenadas z e zo só precisam ser distinguidas quando tratamos do movimento das partículas, como na seção 4.1.
e p0(z) = -z desde que consideramos a pressão atmosférica nula e z < O abaixo da superfície livre do líquido. Sendo δp(x, z, t) = p(x, z, t) - p0 a perturbação na pressão, temos que
Utilizando a equação de Bernoulli (A.8), obtemos
que, juntamente com as Eqs. (30) e (35) resulta na seguinte expressão para as variações na pressão do líquido que ocorrem pela passagem da onda
Desta forma, com um sensor colocado na profundidade -z medimos a pressão local p(x, z, t) e, através da Eq. (38), determinamos as variações de pressão δp(x, z, t) que podem ser confrontadas com a Eq. (40) para obtermos, por exemplo, as amplitudes ς0 das ondas que se propagam na superfície.
5. Aproximação dos resultados para diferentes profundidades
Nesta seção, vamos adaptar os resultados gerais obtidos para locais com grandes ou pequenas profundidades h relativas ao comprimento de onda λ da onda que se propaga.
5.1. Ondas em águas profundas
Inicialmente, vamos analisar quais são as características das ondas que se propagam em regiões em que a profundidade é muito maior que o comprimento de onda, isto é, em que kh = 2πh/λ » 1.
Na Fig. 3 verificamos a função tghx → 1 para grandes valores de x. No entanto não é necessário que x seja muito grande para que isto ocorra; por exemplo, para x = kh = 1,75, que corresponde a h/λ = 0,28, já obtemos tghx = 0,94. Com esta aproximação, podemos escrever a equação para a velocidade de fase (36)
Esta equação estão representada pela curva (b) da Fig. 2 e exibe o fato que as ondas com maiores comprimentos de onda se deslocam mais rapidamente, isto é, apresentam o fenômeno da dispersão.
Se estamos sentados à beira da praia e percebemos a chegada de cristas de ondas em intervalos regulares, por exemplo, de T = 8 s, então, da Eq. (41) e de c = λ/T, concluimos que o o comprimento de onda correspondente é
e velocidade de fase c ~ 45 km/h. Observemos que este comprimento de onda se refere às regiões distantes da praia, onde devemos ter a profundidade h > 28 m; próximo da praia o comprimento de onda será muito menor (veja seção 5.2), embora o período da onda permaneça inalterado, de modo análogo ao de ondas eletromagnéticas que passam através de meios com índices de refração diferentes.
Para águas profundas as trajetórias das partículas de líquido são dadas, através das Eqs. (32) e (34) e com senh x ≈ cosh x ≈ ex, por
as quais representam circunferências centradas em (x0, z0) e raios r(z) = ς0 eKZ que decrescem com a profundidade, como mostrado na Fig. 4(a). Para uma profundidade igual à metade do comprimento de onda, temos que
de modo que o raio da trajetória é apenas 4 por cento do seu valor na superfície.
Para a variação de pressão nesta situação, temos das Eqs. (40) e (24),
e, da mesma maneira que para o raio, para a profundidade z = -λ/2 seu valor cai a 4 por cento do valor na superfície. Assim, um detector colocado no fundo de um lago não poderá detectar variações na pressão devido às ondas superficiais de gravidade cujos comprimentos de onda sejam muito menores que o dobro da profundidade do lago.
5.2. Ondas em águas rasas
Em seguida, vamos analisar as características das ondas que se propagam em regiões em que a profundidade é muito menor que o comprimento, isto é, k h = 2πh/λ « 1.
Para pequenos valores de x, tgh x ~ x e, da Fig. 3 que, para x < 0,45 esta já é uma boa aproximação. Assim, a velocidade de fase (36) pode ser aproximada por
a qual independe do comprimento de onda, aumenta com a profundidade do líquido e está representada pela curva (c) da Fig. 2. Portanto, para águas rasas devemos ter, efetivamente, λ » h enquanto que no caso de águas profundas não necessariamente precisamos ter λ « h (veja Eq. (41) para a validade das aproximações utilizadas.
Para ondas que chegam em uma praia com profundidade h = 1 metro com intervalos regulares de 8 s, a velocidade de propagação é de c = 11 km/h é o comprimento de onda é λ = 25 m. Neste exemplo, as ondas apresentam uma redução de 75 por cento na velocidade de propagação e no comprimento de onda quando passam de águas profundas para regiões com pouca profundidade (veja seção 5.1).
Para pequenos valores de x, senh x ~ x e cosh x ~ 1; assim, as componentes de velocidade (31), são aproximadas por
sendo a componente vertical nula para z = -h.
Integrando as Eqs. (47) em relacão ao tempo, obtemos para as trajetórias das partículas
que representam elipses com centro em (x0, z0) em que o semi-eixo maior é constante e o semi-eixo menor diminui com a profundidade até se anular para z0 = -h, como ilustrado na curva (b) da Fig. 4.
Para esta região, a variação linearizada para a pressão (40) torna-se
que independe da profundidade. Com a Eq. (24) reescrevemos a equação anterior como
cuja forma é a mesma da pressão hidrostática gz. Neste caso, ao contrário do que ocorre quando vale a aproximação de águas rasas Eq. (45), um sensor de pressão colocado em z = -h pode, em princípio, detectar as flutuações na pressão, uma vez que a forma de δp não é atenuada pela funçao eKZ.
6. Considerações finais
No desenvolvimento deste trabalho pudemos mostrar que, através de um potencial escalar que é solução de uma equação de onda, é possível determinar as velocidades e trajetórias das partículas de um líquido incialmente em repouso e com profundidade constante, percorrido por uma perturbação sinusoidal na superfície líquida bem como, as flutuações de pressão nos diferentes pontos do fluido. Verificamos, através de uma teoria linearizada, que as ondas com pequenas amplitudes que se propagam no líquido, o fazem com diferentes velocidades que dependem da profundidade local. Este fato permite-nos compreender, por exemplo, porque as frentes de onda chegam à praia praticamente paralelas à linha da praia. Quando as ondas provenientes das regiões mais profundas, portanto com com maiores velocidades de fase, se aproximam das regiões mais rasas, experimentam uma diminuição de velocidade, de modo que as ondas sofrem seguidas refrações e passam a propagar-se numa direção que é perpendicular à linha da praia, como mostrado na Fig. 5.
Recebido em 16/7/2008
Revisado em 22/1/2009
Aceito em 6/2/2009
Publicado em 26/6/2009
Apêndice: Equação de Bernoulli
A equação de conservação do momento linear para um fluido não viscoso com densidade g(r, t), pressão p(r, t) e submetido à forca externa específica f(r), também denominada equação de Euler, é
onde v( r, t) representa as componentes da velocidade do fluido. Se a força externa específica é conservativa podemos relacioná-la ao potencial escalar φ(r, t) por
Para um fluxo irrotacional, a identidade vetorial v x (∇ x v) = ∇(v2/2) - (v · ∇)v = 0, permite escrever que
de modo que, com as Eqs.(A.2) e (A.3), temos para a Eq. (A.1)
Podemos introduzir o potencial escalar Φ(x,t) para velocidade através da Eq. (9) e, para um fluido incompressível (g = const.), temos para a equação anterior
ou seja,
onde a função F(t) não depende da posição e pode ser absorvida pela redefinição do potencial Φ (r, t) de modo que o campo de velocidades da Eq. (9) permanece inalterado.
Para o campo gravitacional
é a energia potencial gravitacional por unidade de massa e a Eq. (A.6) pode ser escrita como
que é a forma da equação de Bernoulli dependente do tempo, aplicável a um fluido incompressível, não viscoso situado num campo gravitacional e que experimenta um fluxo irrotacional.
- [1] S. Meirelles e N. Violante Carvalho, Revista Brasileira de Ensino de Física 29, 4 (2007).
- [2] M. Alonso e E.J. Finn, Física (Addison-Wesley Longman, São Paulo, 1992).
- [3] J. Lighthill, Waves in Fluids (Cambridge University Press, Cambridge, 1978).
- [4] M.S.D. Cattani, Elementos de Mecânica dos Fluidos (Edgard Blucher, São Paulo, 1990).
- [5] K.R. Symon, Mecânica (Editora Campus, Rio de Janeiro, 1982).
Datas de Publicação
-
Publicação nesta coleção
08 Out 2009 -
Data do Fascículo
Jun 2009
Histórico
-
Aceito
06 Fev 2009 -
Revisado
22 Jan 2009 -
Recebido
16 Jul 2008
