Esquemas Numéricos Eficientes para as Equações de
Navier-Stokes Bidimensionais Incompressı́veis em Notação
Complexa
Fabiano Barbosa Mendes da Silva∗
Jorge Ferreira
Universidade Federal Rural de Pernambuco, UFRPE - UAG,
55296-901, Garanhuns - PE
E-mail: [email protected], [email protected],
Sı́vio Marques de Almeida Gama
Universidade do Porto - Faculdade de Ciências - DMA,
4169-007, Porto, Portugal
E-mail: [email protected],
Ana Carla Percontini da Paixão
Universidade Estadual de Feira de Santana - UEFS
44036-900, Feira de Santana - BA
E-mail: [email protected].
RESUMO
Os fluidos incompressı́veis definem-se como sendo as soluções das equações de Navier-Stokes
que, em notação Euleriana [4], escrevem-se:

2

 ∂t u + (u.∇)u = −∇p + ν∇ u + f,
(1)

 ∇.u = 0,
suplementadas com uma condição inicial e condições fronteiras dadas (por exemplo, condições
fronteiras periódicas ou anulação de u nas vizinhanças de infinito). Consideramos V o volume
ocupado pelo fluido, um subconjunto compacto e conexo d-dimensional de Rd independente do
tempo, supondo que a força é incompressı́vel, isto é, ∇.f = 0.
Escrevemos as equações de Navier-Stokes bidimensionais incompressı́veis (NS2D) considerando
as componentes u1 e u2 do campo de velocidades bidimensional como a parte real e imaginária
de um campo de velocidades complexo. Para isto, introduzimos o campo complexo de velocidades
χ = u1 + iu2 e finalmente escrevemos as equações NS2D em notação complexa [3], isto é, à custa
de χ:

1
2
0
2

 ∂t χ + 2 (∂1 − i∂2 )χ = −(∂1 + i∂2 )p + ν∇ χ,
(2)

 Re [(∂ − i∂ ) χ] = 0,
1
2
onde p0 designa-se por pressão modificada. A vantagem de termos NS2D escritas em notação
complexa reside na particularidade de os termos não-lineares (T N L) serem simplesmente o
quadrado do campo complexo χ. Isto possibilita o uso de métodos pseudo-espectrais que requerem somente duas FFT’s (transfomadas de Fourier rápidas) complexas-complexas por passo
∗
Doutorando em Matemática Computacional - UFPE
1173
de tempo para o cálculo dos termos não-lineares. Se considerássemos NS2D forçadas, a equação
(2) teria o termo adicional F = f1 + if2 , onde (f1 , f2 ) são as componentes da força. Assumindo periodicidade L = 2π em cada eixo de coordenadas e supondo que o campo complexo de
P
bk (t)eik·x , no espaço de Fourier, a equação (2) escreve-se:
velocidades χ = χ(x1 , x2 ; t) = k∈Z 2 χ
bk = −νk 2 χ
bk + Td
∂t χ
N Lk ,
(3)
sendo
1
c2 − 1 (ik − k )3 χ
d
2∗ ,
(4)
Td
N Lk = − (ik1 + k2 )χ
1
2
k
k
4
4k 2
onde k = (k1 , k2 ) é o vetor de onda e k 2 = k.k = k12 + k22 é o quadrado do número de onda k.
Saliente-se que NS2D, quando escritas em termos de função corrente, requerem o uso de quatro
FFT’s reais-complexas [1], o que é menos performante que duas FFT’s complexas-complexas.
Com efeito, uma FFT complexa-complexa é equivalente a duas FFT’s reais-complexas mais um
interface que prepara os dados reais em complexos, já que toda a transformada de Fourier rápida
real-complexa passa necessariamente pelo uso de transformadas rápidas complexas-complexas
[5]. O esquema numérico para a evolução temporal de (3) é o slaved-frog [2]:
1 − e−2αδt
fn ,
(5)
n > 0,
α
tal que δt é o passo de tempo, qn = q (nδt) e fn = f (nδt) . Este esquema numérico apresenta um
erro de truncatura local de 3a ordem no tempo. Aplicando-o a (3), e escrevendo adequadamente
seus termos através de quantidades pré-calculadas, obtemos:
qn+1 = e−2αδt qn−1 +
c2 (n) + B χ
d
2∗ (n).
bk (n + 1) = M χ
bk (n − 1) + B1 χ
χ
2
k
k
(6)
c2 e χ
d
2∗ não são complexos conjugados entre si. Porém, relacionam-se do seguinte
Notar que χ
modo:
d
c2 (−k , −k ) ∗ .
2∗ (k , k ) = χ
(7)
χ
1 2
1
2
Essa igualdade decorre da simetria Hermitiana. Notar que, se k = 0, ou seja k12 + k22 = 0, em
(6) tem-se uma indeterminação no segundo membro da igualdade. Calculando o limite k → 0,
desenvolvendo o segundo termo do lado direito da igualdade em série de Taylor, eliminando
assim a indeterminação, concluı́mos que:
b 0 (n + 1) = χ
b 0 (n − 1).
χ
(8)
b 0 é constante.
Desse modo, χ
Palavras-chave: Navier-Stokes, Notação complexa, Esquemas numéricos.
Referências
[1] C. Basdevant, Thecnical improvements for direct numerical simulations of homogeneous
thee-dimensional turbulence, J. Comp. Phys., 50 (1983) 209-214.
[2] U. Frisch, Z.S. She & O. Thual, Viscoelastic Behavior of cellular solutions to the KuramotoSivashinsky model, J. Fluid. Mech., 168 (1986) 221-240.
[3] S. Gama, U. Frisch, H. Scholl, The 2-D Navier-Stokes equations with a large scale instability
of Kuramoto-Sivashinsky type: numerical exploration on the Connexion Machine, J. Sci.
Comp., 6 (1991) 425-452.
[4] S.B. Pope, “Turbulent Flows”, Cambridge University Press, 2000.
[5] W.H. Press, S.A. Teukolsky, W.T. Vetterlig & B.P. Flannery, “Numerical Recipes in Fortran: The Art of Scientific Computing”, Cambridge University Press, 1992.
1174
Download

2010Métodos Numéricos e Aplicações P