Tese Duque
Tese Duque
Tese Duque
Ciências
ii
Método dos Elementos Finitos para problemas com Fronteiras Livres
iii
Método dos Elementos Finitos para problemas com Fronteiras Livres
iv
Método dos Elementos Finitos para problemas com Fronteiras Livres
Agradecimentos
Finda esta etapa, desejo manifestar o meu sincero agradecimento a todas as pessoas que,
de algum modo, contribuíram para a sua realização, nomeadamente:
- à minha esposa Ana e à minha lha Matilde pela compreensão das minhas ausências e tam-
bém por sempre me apoiarem e encorajarem, especialmente nos momentos mais difíceis;
- aos meus pais por me incentivarem a seguir esta área e por me terem dado condições para
iniciar os meus estudos académicos superiores;
- ao meu co-orientador, Prof. Dr. Rui Almeida pela orientação, permanente disponibilidade
e pelos valiosos conselhos, críticas e correções;
- ao meu orientador, Prof. Dr. Stanislav Antontsev pela orientação, ensinamentos, enorme
disponibilidade e, principalmente, pela imensa paciência perante as minhas lacunas;
- ao Prof. Dr. Jorge Ferreira, da Universidade de Pernambuco, pela colaboração e por ter
dado um novo ânimo ao trabalho numa fase crítica.
A todos um Bem-Haja.
v
Método dos Elementos Finitos para problemas com Fronteiras Livres
vi
Método dos Elementos Finitos para problemas com Fronteiras Livres
Resumo
Neste trabalho pretende-se fazer um estudo sobre a aplicação do método dos elementos
nitos a diversos problemas de reação-difusão com fronteiras livres. Para se obter estimativas
do erro e simulações representativas é necessário obter alguns resultados teóricos sobre regu-
laridade e algumas propriedades físicas das soluções. Neste sentido, o outro objetivo deste
trabalho é o estabelecimento de resultados teóricos relativos a problemas em aberto.
O primeiro problema a ser estudado é a equação parabólica seguinte:
Como o problema pode ser degenerado, utiliza-se um problema aproximado, regularizado atra-
vés da introdução de um parâmetro ε. Demonstra-se, sob algumas condições em γ e f que
a solução fraca do problema aproximado converge para a solução fraca do problema inicial,
quando o parâmetro ε tende para zero. São calculadas soluções discretas utilizando o método
dos elementos nitos contínuo no espaço e descontínuo no tempo e é provada a convergência
destas soluções para a solução fraca do problema inicial.
Estuda-se também a aplicação do método da malha móvel a esta equação considerando-a
R
num domínio livre em 2 . É desenvolvido um conjunto de subrotinas em Matlab que permitem
calcular e representar gra camente soluções aproximadas de vários problemas de reacção-difu-
são com fronteiras livres. A discretização espacial é de nida por uma partição do domínio em
triângulos. Em cada elemento nito, a solução é aproximada por uma função seccionalmente
polinomial de grau r ≥ 1 utilizando polinómios interpoladores de Lagrange em coordenadas de
área. Os vértices dos triângulos podem mover-se segundo um sistema de equações diferenciais
parciais que é adicionado ao problema. Posteriormente, deduz-se uma equação para mover os
vértices da fronteira. O sistema resultante é convertido num sistema de equações diferenciais
ordinárias no tempo, que é resolvido utilizando um integrador apropriado. Os integrais que
surgem são calculados utilizando a quadratura de Gauss. Finalmente, são apresentados alguns
resultados de aplicação.
O outro problema estudado é o sistema não linear da forma
(
ut = a1 (l1 (u), l2 (v))∆u + λ1 |u|p−2 u + f1 (x, t)
, x ∈ Ω ⊂ Rd , t ∈]0, T ],
vt = a2 (l1 (u), l2 (v))∆v + λ2 |v|p−2 v + f2 (x, t)
Palavras-chave
Método dos elementos nitos, equação de meios porosos, expoente variável, malha móvel,
fronteira livre, termo difusivo não local, efeitos de localização, comportamento assintótico.
vii
Método dos Elementos Finitos para problemas com Fronteiras Livres
viii
Método dos Elementos Finitos para problemas com Fronteiras Livres
Abstract
The aim of this work is to study the application of the nite element method to various
reaction-diffusion problems with free boundaries. To obtain error estimates and representative
simulations, it is necessary to obtain some theoretical results on regularity and some physical
properties of the solutions. In this sense, the other goal of this work is to establish theoretical
results concerning open problems.
The rst problem to be studied is the following parabolic equation:
Since the problem may be of degenerate type, an approximate problem is used, regularized by
introducing a parameter ε. It is proved, under certain conditions on γ and f , that the weak
solution of the approximate problem converges to the weak solution of the initial problem,
when the parameter ε tends to zero. Discrete solutions are built using the continuous in space
and discontinuous in time nite element method, and the convergence of these to the weak
solution of the initial problem is proved.
The application of the moving mesh method to this problem, considering it a free boundary
problem is also studied. A set of subroutines in Matlab, that can calculate and graphically
represent approximate solutions of various reaction-diffusion problems with free boundaries, is
developed. The spatial discretization is de ned by a triangulation of the domain. In each nite
element, the solution is approximated by a piecewise polynomial function of degree r ≥ 1,
using Lagrange interpolating polynomials in area coordinates. The vertices of the triangles
move according to a system of differential equations which is added to the equations of the
problem. The resulting system is converted into a system of ordinary differential equations in
time variable, which is solved using a suitable integrator. The integrals that arise in the system
of ordinary differential equations are calculated using the Gaussian quadrature. Finally, some
numerical results of application are presented.
The other problem to be studied is the nonlinear nonlocal reaction-diffusion coupled system
(
ut = a1 (l1 (u), l2 (v))∆u + λ1 |u|p−2 u + f1 (x, t)
, x ∈ Ω ⊂ Rd , t ∈]0, T ],
vt = a2 (l1 (u), l2 (v))∆v + λ2 |v|p−2 v + f2 (x, t)
where a1 and a2 are Lipschitz-continuous positive functions, l1 and l2 are continuous linear
forms, λ1 , λ2 ≥ 0 and p ≥ 2. The existence and uniqueness of weak and strong solutions to
these systems and localization properties of the solutions, including the waiting time effect
and stable localization, are proved. Moreover, important results on polynomial and exponential
decay and the vanishing of the solutions in nite time, are also presented. The convergence
of a linearized Euler-Galerkin nite element method is proved, and the optimal order of con-
vergence is obtained. Finally, the presented method is implemented and simulated in Matlab
environment.
Keywords
Finite element method, porous media equation, variable exponent, moving mesh, free boun-
dary, non local diffusion term, localization effects, asymptotic behavior.
ix
Método dos Elementos Finitos para problemas com Fronteiras Livres
x
Método dos Elementos Finitos para problemas com Fronteiras Livres
Índice
1 Introdução 1
Bibliogra a 107
xi
Método dos Elementos Finitos para problemas com Fronteiras Livres
xii
Método dos Elementos Finitos para problemas com Fronteiras Livres
Lista de Figuras
3.1 Malha espacial para h = 0.1 usada no problema com o segundo membro arti cial. 56
3.2 Evolução no tempo da aproximação obtida no problema com o segundo membro
arti cial. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.3 Erro da aproximação obtida no problema com o segundo membro arti cial em
função de h e de δ. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.4 Erro da aproximação obtida no problema com o segundo membro arti cial em
função de ε. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
3.5 Evolução no tempo da aproximação obtida no problema com suporte não convexo. 58
3.6 Malha espacial para h = 1/15 usada no problema com absorção. . . . . . . . . . 59
3.7 Evolução no tempo da aproximação obtida no problema com absorção. . . . . . 59
3.8 Cálculo do vetor normal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66
3.9 Domínio computacional e malha inicial usados no problema com expoente cons-
tante. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
3.10 Evolução no tempo da aproximação obtida no problema com expoente constante. 67
3.11 Evolução no tempo da malha obtida no problema com expoente constante. . . . 67
3.12 Erro das aproximações da solução e da fronteira, obtidas no problema com expo-
ente constante, em função do número de triângulos. . . . . . . . . . . . . . . . 68
3.13 Domínio computacional e malha inicial usados no problema com suporte não con-
vexo. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
3.14 Evolução no tempo da aproximação obtida no problema com suporte não convexo. 69
3.15 Evolução no tempo da malha obtida no problema com suporte não convexo. . . . 69
3.16 Evolução no tempo da aproximação obtida no problema com absorção e expoente
constante. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
3.17 Evolução no tempo da malha obtida no problema com absorção e expoente cons-
tante. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
3.18 Erro das aproximações da solução e da fronteira, obtidas no problema com absor-
ção e expoente constante, em função do número de triângulos. . . . . . . . . . 70
3.19 Domínio computacional e malha inicial usados no problema com absorção e expo-
ente variável. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
3.20 Evolução no tempo da aproximação obtida no problema com absorção e expoente
variável. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
3.21 Evolução no tempo da malha obtida no problema com absorção e expoente variável. 72
xiii
Método dos Elementos Finitos para problemas com Fronteiras Livres
xiv
Método dos Elementos Finitos para problemas com Fronteiras Livres
Capítulo 1
Introdução
Nas últimas décadas, foram utilizadas várias equações para descrever alguns sistemas físicos,
químicos, biológicos e ecológicos. Uma das mais bem sucedidas e cruciais é a equação de
reação-difusão que é uma equação diferencial parcial parabólica do tipo:
∂u
− div(a∇u) − f = 0, (1.1)
∂t
J = −a∇u(x),
1
Método dos Elementos Finitos para problemas com Fronteiras Livres
tn ∈]0, T ]. A solução nal obtida pode então ser Un (x) = U (x, tn ) ou a interpolação linear
t−tn−1
de nida, em cada intervalo de tempo ]tn−1 , tn ], por U (x, t) = tn −tn−1 Un (x) − t−tn
tn −tn−1 Un−1 (x).
No caso mais simples, todos os intervalos têm a mesma amplitude δ. Se a solução exata tiver
pouca regularidade em t, um método viável para integrar no tempo é o método dos elementos
nitos descontínuo que aproxima a solução em cada intervalo de tempo por um polinómio, mas
não obriga a que a solução seja contínua nos instantes tn ∈]0, T ]. Para problemas mais difí-
ceis recorre-se frequentemente a integradores mais e cientes. Geralmente, estes integradores
utilizam diversas fórmulas e intervalos de amplitude variável, para poderem controlar o erro
das estimativas, em cada instante tn . Um exemplo deste tipo de integradores é o ode15s ,
implementado em ambiente Matlab.
O primeiro problema de fronteira livre a tratar neste trabalho surge ao estudar o movi-
mento de um gás ideal barotrópico num meio poroso. Se se supuser que a pressão p depende
explicitamente da densidade u e da temperatura θ da forma
p = p(u, θ) = k̃|u|γ̃(θ(x,t)) ,
então, pela lei da conservação da massa e pela lei de Darcy, obtém-se uma equação do tipo
(1.1) com um termo de difusão
a = k|u|γ(x,t) .
O comportamento das equações da forma (1.3) pode ser bem diferente do caso não degene-
rado do tipo
∂u
∂t = div(a(x, t, u)∇u) + f (x, t, u) em ΩT
u(x, t) = 0 em ΓT , (1.5)
u(x, 0) = u0 (x) em Ω
com
R.
− +
∂a ∂f
0 < a ≤ a ≤ a ≤ ∞, + ≤ C,
∂u ∂u ∀(x, t) ∈ Ω̄×]0, T ], ∀u ∈
Se os dados iniciais forem su cientemente regulares, então o problema (1.5) admite soluções
fracas únicas e estas mantêm-se regulares (ver [52]). É conhecido que dados iniciais em L2 ou
distribuições produzem soluções regulares.
Aplicando o MEFG a problemas do tipo (1.5) conseguem-se erros de ordem ótima. Em [72],
2
Método dos Elementos Finitos para problemas com Fronteiras Livres
1
(Un − Un−1 , χ) + (a(Un )∇Un , ∇χ) = (fn , χ), χ ∈ Sh
δ
prova-se que:
kUn − u(tn )kL2 (Ω) ≤ C(h2 + δ),
n ≥ 1.
com
1
c= e m = γ + 1.
γ+1
O coe ciente c não é essencial e, sem perda de generalidade, pode considerar-se c = 1. Se se
considerar f = f (x, t), o problema (1.6), em geral, não admite soluções clássicas. A existên-
cia e unicidade de soluções fracas para este tipo de problemas foi demonstrada em dimensão
um por Ole nik, Kala²inkov e ƒºou ([61]) e para dimensões superiores a um, por Lions ([55]).
Lions provou que se f ∈ L m+1 (ΩT ) e u0 ∈ L2 (Ω), então existe uma única solução do pro-
m
m−1
blema (1.6) pertencente a L∞ (0, T ; L2 (Ω)). Provou ainda que |u| 2 u ∈ L2 (0, T ; H01 (Ω)) e
−1, m+1
ut ∈ L m+1 (0, T ; W m (Ω)).
m
Uma diferença deste tipo de equações para as não degeneradas é que soluções iniciais
suaves, ou até mesmo analíticas, não produzem necessariamente soluções suaves. Foi demons-
trado em [46] que, para m ≥ 2, dados iniciais suaves e com suporte compacto nunca produzem
soluções em C 1 . A derivada espacial torna-se descontínua em algum instante nito.
Aronson ([6]) provou que se um−1
0 é Lipschitz-contínua em x, então um−1 também é Lips-
chitz-contínua e u Hölder-contínua em x com expoente min{1, m−1
1
} e em t com expoente
1
min{ 13 , 2m−1 }. Provou ainda que (um )x existe e é contínua em x e se 1 < m < 2, então ux
também existe e é contínua em x.
Ole nik, Kala²inkov e ƒºou provaram também que se u0 tem suporte compacto, então u tem
suporte compacto para qualquer tempo positivo. Esta propriedade é usualmente denotada por
propagação com velocidade nita ( nite speed of propagation ). Kala²inkov [46] estudou
então o caso em que u0 satisfaz as condições:
é limitado por duas curvas x = Γ1 (t) e x = Γ2 (t), onde Γ1 (t) é contínua, monótona não crescente
e contém o ponto (a1 , 0) e Γ2 (t) é contínua, monótona não decrescente e contém o ponto (a2 , 0).
Às curvas Γ1 (t) e Γ2 (t), geralmente, chamamos interfaces ou fronteiras livres.
As questões interessantes sobre a regularidade aparecem junto à interface que é melhor
3
Método dos Elementos Finitos para problemas com Fronteiras Livres
Mais tarde, Knerr ([49]) provou que Γ1 e Γ2 são de Lipschitz e que a derivada lateral de Γi
existe e satisfaz
Γ0i (t) = −νx (Γi (t), t), t ∈ [0, T ]. (1.7)
A equação (1.7) deve ser entendida como o limite quando x → Γi (t). Knerr provou também que
se ν0 (x) ≥ (a2 − x)$ , para a2 − ≤ x ≤ a2 com 0 < < 1 e $ < 2, então ν(a2 , t) > 0, t > 0 e
que se ν0 (x) ≤ (x − a2 )2 , para a2 − < x < a2 e > 0, então existe uma constante tw tal que
Γ2 (t) = a2 , para 0 ≤ t ≤ tw . Este resultado cou conhecido por tempo de espera ( waiting time
effect ) e, grosso modo, indica que a interface pode começar instantaneamente a mover-se
ou pode haver um intervalo de tempo nito [0, tw ] no qual a fronteira permanece estacionária
e passado esse tempo começa a mover-se com velocidade nita.
Em 1987, Aronson e Vazquez ([8]) provaram que Γi é uma curva C ∞ no intervalo ]tw , +∞[.
No caso tw = 0, a interface é C ∞ ( R+). Se tw > 0, esta regularidade só não é veri cada em
t = tw , pois Γi (t) = ai , t ∈ [0, tw ].
Angenant ([2]) melhorou estes resultados mostrando que no intervalo ]tw , +∞[, Γi é uma
função analítica.
Para dimensões superiores a um, Caffarelli e Friedman ([14]) provaram que a superfície
Γ(t) = ∂Ω(t) é Hölder-contínua e Caffarelli e Wolanski ([15]) provaram que, impondo algumas
condições em ν, existe 0 < κ < 1 tal que Γ(t) é de classe C 1,κ .
Daskalopoulos e Hamilton ([23]) provaram que existe um intervalo ]0, T [ no qual Γ é uma
superfície de classe C ∞ , mas as condições são: ν0 , Dν0 e dD2 ν0 se puderem estender conti-
nuamente até à fronteira por funções Hölder-contínuas de classe C κ , κ > 0, Dν0 6= 0 em Γ e
ν0 > 0. d denota a distância à fronteira do domínio. A condição Dν0 6= 0 em Γ garante que a
interface começa a mover-se em t = 0.
Ko ([51]) estendeu este resultado para t > T no caso de a solução ter simetria radial.
Para f = f (u) e γ constante os casos em que f (u) = −up são os mais usuais, obtendo-se o
problema:
∂u
∂t = c∆(um ) − up em ΩT
u(x, t) = 0 em ΓT . (1.8)
u(x, 0) = u0 (x) ≥ 0 em Ω
0 é Lipschitz-contí-
Kala²inkov ([47]) provou que existe solução fraca e que é única se p ≥ 1 e um
nua e que (um−1 )x é limitada para t > 0. Destes resultados, conclui-se que u é Hölder-contínua
em ΩT relativamente a x com expoente 1
m e, consequentemente, relativamente a t com expo-
ente 2m+1 .
1
Tal como no caso da PME, a solução tem as derivadas ut , ux e uxx , limitadas no
conjunto Ω(t), e as questões principais de regularidade surgem na fronteira desse conjunto.
Em [47], provou-se que no caso p > 1 a solução possui a propriedade de propagação com
velocidade nita. Se 1 < p < m, então Ω(t) é limitado para t ≥ 0 e no caso p > m, Ω(t) expande
4
Método dos Elementos Finitos para problemas com Fronteiras Livres
in nitamente.
Knerr ([50]) provou que Ω(t) é um intervalo ]Γ1 (t), Γ2 (t)[ e que Γi (t) ∈ C 0 ([0, +∞[) ∩
C 0,1
(]0, +∞[), no caso m > 1 e p ≥ 1. Com base nestes resultados, Herrero e Vazquez ([42])
deram condições para a existência de tempo de espera tw e provaram que Γi (t) ∈ C 1 (]tw , +∞[).
Demonstraram também que Γi (t) satisfaz a equação (1.7), no caso p ≥ m. Shmarev e Vazquez
([69]) provaram que a equação (1.7) também é válida para 1 < p ≤ m.
Se 0 < p < 1, a equação (1.7) deixa de ser válida, pois o suporte da solução encolhe com
o tempo e a solução extingue-se num intervalo de tempo nito, ou seja, existe te > 0 tal que
Ω(t) 6= ∅, t < te e Ω(t) = ∅, t ≥ te . A esta propriedade chamamos extinção em tempo nito
( nite time extinction ). Em ([38]), provou-se que para p < 1 e m + p > 2, Γi (t) satisfaz (1.7)
quando o domínio expande e satisfaz a equação:
1−p
Γ0i (t) = , t ∈ [0, T ] (1.9)
(u1−p )x (Γi (t), t)
1−p
Γ0i (t) = −νx (Γi (t), t) +
(u1−p )x (Γi (t), t)
e que o domínio expande até um determinado instante e depois contrai até se extinguir num
tempo nito. Neste caso, provou-se que a fronteira é analítica.
O caso mais curioso surge quando m > 1, 0 < p < 1 e 0 < m + p < 2 onde se prova ([39]) que
a fronteira segue a equação:
1 β
Γ0i (t) = (α − (u(m−p)/2 )x )1/β (Γi (t), t),
κ(1 + β) x
com as constantes
2−m−p
(m − p)2
m−p m−p m−p
2−m−p
α= p , κ= α , β= .
2m(m + p) 2m(1 − p)(m + p + 2) 2m m−p
Para dimensões superiores a um, Shmarev ([67, 68]) provou que se m > 1, p > 0 e m + p ≥ 2,
a velocidade da fronteira Γ(t) é dada por
0 m m−1
ν = Γ (t) = lim − ∇u + ∇Υ ,
x→Γ m−1
div(u∇Υ) = up , Υ = 0 em Γ
e provou que u e Γ são funções analíticas em t e preservam a regularidade inicial nas variáveis
espaciais. A regularidade de ν em relação às variáveis espaciais é melhor que no instante
inicial.
A procura de soluções na forma explícita é uma área importante no estudo de equações
diferenciais. Em 1952, Barenblat ([13]) deduziu e estudou uma solução de (1.6), com c = 1 e
5
Método dos Elementos Finitos para problemas com Fronteiras Livres
com propriedades diferentes. Esta solução não é limitada, tende para in nito quando t → T − .
Para o problema (1.8) em dimensão um, Kersner ([48]) encontrou uma solução explícita para
m > 1, 0 < p < 1 e m + p = 2 da forma:
1
m−1
m−1 h 2
i
u(x, t) = Ct m+1 − (m + 1)2 t2 − x2 .
2m(m + 1)t +
Um problema do tipo (1.3) foi primeiro estudado por Antontsev e Shmarev [4], para γ depen-
dente das variáveis espaciais e do tempo. Nesse artigo foi tratado o problema com f = f (x, t)
e −1 < γ − ≤ γ(x, t) ≤ γ + < ∞. Provaram que se
Z T
ku0 kL∞ (Ω) + kf kL∞ (Ω) dt ≤ ∞,
0
então o problema (1.3) tem pelo menos uma solução pertencente a L∞ (0, T ; L∞ (Ω)) tal que
γ
|u| 2 ∇u ∈ L2 (0, T ; L2 (Ω)) e ut ∈ L2 (0, T ; H0−1 (Ω)). No caso γ − > 0, os autores conseguiram
provar que se supx∈Ω̄ |∇γ| ∈ L2 (]0, T ]), então a solução é única e no caso γ + ≤ 0, tem-se que
k∇ukL2 (ΩT ) ≤ C. Antontsev e Shmarev estudaram também as propriedade físicas das soluções
deste tipo de problemas. Concluíram que se γ é contínuo e positivo, então u possui a propri-
edade de propagação com velocidade nita. Estabeleceram ainda condições a impor a u0 e f
su cientes para que exista tempo de espera. Eles provaram que se γ + < 0, então u tem a
propriedade de extinção em tempo nito. Mais tarde, em [5], Antontsev e Shmarev demons-
traram a existência e unicidade de solução fraca para um problema mais geral considerando
f = b|u|γ(x,t)/2 ∇u − c|u|σ(x,t)−2 u + d.
6
Método dos Elementos Finitos para problemas com Fronteiras Livres
Em 2006, Henriques e Urbano ([41]) provaram que as soluções de (1.3), com γ − > 0, são lo-
calmente contínuas demonstrando que é possível de nir uma sequência de cilindros encaixados
Qn e uma sequência wn convergente para zero tais que:
ut = div(|u|m ∇u), m ≥ 1,
foi obtida por Rose ([64]) que utilizou um esquema regularizado, substituindo o coe ciente de
difusão |u|m por um novo coe ciente da forma aε (u) que satisfaz as condições:
R
aε (α) ∈ C 4 ( ), ε ∈]0, 1]
m
aε (α) = |α| , α≥ε
aε (α) ≥ ε/2, α≥0 .
a0ε (α) ≥ 0, α≥0
R
aε (−α) = aε (α), α∈
Rose aplicou, ao problema regularizado, o método dos elementos nitos de Galerkin contínuo
com aproximações lineares no espaço e o método de Euler regressivo no tempo. Escolhendo
2
0 < ε ≤ Ch m+1 , conseguiu provar uma estimativa da forma
1
! m+2
X 1 κ 2
δ ku(tn ) − Un km+2
Lm+2 (Ω) ≤ C(δ m+1 + (ln(1/h)) (m+1)(m+2) h m+1 ),
n
para dimensões 1, 2 e 3, supondo que ut ∈ L m+2 (0, T ; L m+2 (Ω)) e que a triangulação é quase
m+1 m+1
uniforme. Para dimensão 1, κ = 0, para dimensão 2 e 3, κ = 1.
Nochetto e Verdi ([59]) estudaram a equação:
1
ut = ∆v, u ∈ k(v) = v m
1
(Un − Un−1 , χ) + (∇Vn , ∇χ) = 0 χ ∈ Sh , n ≥ 1,
δ
7
Método dos Elementos Finitos para problemas com Fronteiras Livres
min(s/ε, k(s)), s > 0
Un = Πh kε (Vn ), kε (s) = 0, s=0 .
max(s/ε, k(s)), s < 0
4m 2m−2
Para δ = Ch 3m−1 e ε = Ch 3m−1 provaram a estimativa:
4m
ku − U kLm+2 (ΩT ) ≤ Ch (m+1)(3m−1) ,
A equação:
1
ut = ∆v + f, u ∈ k(v) = v m
1
(Un − Un−1 , χ) + (∇Vn , ∇χ) = (f n/2 , χ) χ ∈ Sh ,
δ
Un = Πh kε (Vn ), n ≥ 1.
ku − U kL∞ (0,T ;H −1 (Ω)) + kv − V kL2 (0,T ;L2 (Ω)) ≤ C(δ + (ln(1/h))α h).
Utilizando o mesmo esquema que Nochetto e Verdi e usando uma propriedade de não dege-
nerescência na fronteira, Ebmeyer [34] provou uma estimativa do tipo
m2 +6m+8
ku − U kL2 (0,T ;L2 (Ω)) ≤ Ch 6m2 +2m ,
5m2 +4m−1
escolhendo δ ≤ Ch 3m2 +m . É de salientar que esta estimativa é limitada inferiormente por 16 ,
ao contrário das anteriores que tendiam para zero quando m tendia para in nito.
1 1
(Un − Un−1 , χ) + (∇(|Un |m−1 Un ), ∇χ) = (fn , χ) χ ∈ Sh , n ≥ 1,
δ m
quando m ≥ 1.
Mais recentemente, Ebmeyer e Liu ([35]) de niram uma nova seminorma e utilizaram o
esquema:
1
(Un − Un−1 , χ) + (∇Vn , ∇χ) = 0 χ ∈ Sh , n ≥ 1.
δ
8
Método dos Elementos Finitos para problemas com Fronteiras Livres
onde v = |u|m−1 u.
Contudo, ainda não se fez um estudo sobre a convergência, análise do erro e efeitos nu-
méricos para o método dos elementos nitos, quando aplicado a equações do tipo (1.3) com
expoente variável.
Antontsev e Shmarev [4] provaram que as soluções do problema (1.3) têm a propriedade de
propagação com velocidade nita. Se a solução inicial tiver suporte compacto, então origina-se
uma fronteira livre que delimita o suporte da solução. Assim, podemos tratar o problema como
um problema de fronteira livre. Surge então o problema de calcular o domínio Ω(t) e a função
u.
Para abordar o problema em domínio livre é necessário mover a fronteira de um modo
adequado. Desta forma, é particularmente útil utilizar um método onde a malha se possa
mover com a fronteira, mas também que se adapte à solução obtida no interior.
O Método dos Elementos Finitos com malha adaptativa é amplamente utilizado para resolver
equações diferenciais, em domínios xos, que têm soluções com grandes variações no domínio.
Ficou já evidente que podem ser obtidos resultados mais precisos se a malha concentrar pontos
em regiões onde a solução precisa de maior de nição. Se houver a garantia de que o número de
elementos nitos é su ciente para descrever bem as variações da solução e o movimento dos
nós é adequado, então obtém-se um método extremamente e ciente. A principal di culdade
reside na escolha do movimento dos nós de modo a que não ocorram singularidades.
O Método dos Elementos Finitos Móveis (MEFM) de Miller e Miller [58] é uma extensão natural
do método das linhas para malhas xas. O método foi aplicado a equações do tipo:
ut = L(u),
com L um operador onde apenas surgem as derivadas espaciais. A solução procurada é da forma
X
U = U (x(t), t) = Uj (t)ψj (x(t)),
j
onde x(t) é o interpolador linear da posição dos nós da malha {Xj (t)} e ψj é a função básica de
Lagrange de grau um associada ao nó Xj . A solução e a posição dos nós são obtidos através da
minimização da norma L2 (Ω) do resíduo
Λ = Ut − L(U ).
∂Uj ∂Xj
Após aplicar a regra da cadeia e da minimização em relação a ∂t e ∂t , resultam as
equações: Z Z
∂U ∂X
ψi − · ∇U dx = ψi L(U ) dx
Ω ∂t ∂t Ω
e Z Z
∂U ∂X
(ψi ∇U ) − · ∇U dx = ψi ∇U L(U ) dx.
Ω ∂t ∂t Ω
Estas equações formam um sistema de equações diferenciais ordinárias que tem de ser inte-
grado no tempo considerando como valores iniciais U (0) e a posição inicial da malha. Obtemos,
9
Método dos Elementos Finitos para problemas com Fronteiras Livres
deste modo, em cada instante, a solução do problema e a posição da malha. Em princípio, este
método é ótimo na norma L2 , pois minimiza o resíduo. No entanto, o problema pode tornar-se
degenerado se alguma das componentes de ∇U for contínua em algum nó e também é possível
que a malha mude de topologia, pois a medida de algum dos elementos nitos pode tornar-se
zero ou negativa. Para evitar isto, é necessário adicionar um termo de penalização que vai
tornar o sistema de equações diferenciais ordinárias bastante rígido ( stiff ) e obrigar ao uso
de um integrador mais robusto para manter a e ciência.
Mais tarde, Carlson e Miller ([16, 17]) introduziram uma função peso da forma
1
W= ,
1 + |∇U |2
e Z Z
∂U ∂X
(ψi ∇U ) − · ∇U W dx = ψi ∇U L(U )W dx,
Ω ∂t ∂t Ω
e consideram x como uma aplicação entre dois domínios. Os autores provaram que a solução
10
Método dos Elementos Finitos para problemas com Fronteiras Livres
existe e é regular desde que os dados iniciais e a fronteira sejam regulares. A aplicação é não
singular, mas a malha pode tornar-se bastante distorcida ( skewed ). A solução do problema
inicial e a posição dos nós são calculados separadamente em cada intervalo ]tn−1 − tn ]. Foram
simulados diversos problemas com diferentes funções monitor, mas revelou-se difícil encontrar
uma função, no caso geral, que satisfaça todas as necessidade de adaptação.
O método de Baines, Hubbard e Jimack ([11, 10, 12]) utiliza a lei da conservação geomé-
trica, na sua forma integral, para calcular as velocidades. A malha é acoplada com a solução
e é movida, em cada instante, de forma a procurar conservar no tempo o integral de uma
função monitor igualmente distribuído pelos elementos nitos. A equação diferencial parcial é
utilizada com a versão euleriana da equação da conservação para gerar o campo de velocida-
des, que depois é usado para mover a malha na versão lagrangeana. Nos casos mais simples,
a solução da equação diferencial parcial pode ser reconstruída, à posteriori, através da ver-
são lagrangeana da lei de conservação. Nos casos mais complexos, pode-se utilizar a versão
lagrangeana da equação diferencial parcial. Este método foi aplicado com sucesso a equações
com fronteiras móveis e livres obtendo-se resultados precisos tanto para a solução como para
a posição da fronteira. Recentemente, Marlow, Hubbard e Jimack ([57]) aplicaram o método
a problemas com fronteira móvel e utilizaram, para mover os pontos na fronteira, uma função
monitor que podia ser diferente da função monitor utilizada para mover os pontos interiores.
Neste trabalho, estendem-se os resultados de Rose para o caso do expoente depender das
variáveis espaciais e os resultados de Baines, Hubbard e Jimack para o caso de problemas em
domínios com fronteira livre.
O segundo problema com fronteira livre a ser estudado aqui surge da biologia. Pretende-se
estudar a difusão de uma população (por exemplo de bactérias). Neste caso, é razoável supor
que a mobilidade no meio depende do quanto populoso ele está. Este problema pode ser
modelado pela equação (1.1) onde a depende do total da população nesse domínio Ω, isto é,
Z
a = a(l(u)), com l(u) = u(x, t) dx,
Ω
Considere-se um sistema não linear do tipo reação-difusão, num domínio limitado e com
termos de difusão não lineares aplicados a funcionais lineares li , da forma
(
ut − a1 (l1 (u), l2 (v))∆u + λ1 |u|p−2 u = f1 (x, t) em Ω×]0, T ]
. (1.11)
vt − a2 (l1 (u), l2 (v))∆v + λ2 |v|p−2 v = f2 (x, t) em Ω×]0, T ]
Neste caso, u e v podem descrever a densidade de duas populações que interagem através das
funções a1 e a2 . Supõe-se que a taxa de mortalidade na espécie u é proporcional a |u|p−2 u com
o factor λ1 ≥ 0 e, na espécie v, é proporcional a |v|p−2 v com o fator λ2 ≥ 0. É assumido que
p ≥ 1. f1 e f2 representam as in uências por fontes externas.
Este tipo de termo de difusão foi proposto por Chipot e Rodrigues em 1992 [20] para pro-
blemas elípticos. Nesse trabalho os autores provaram a existência de soluções fracas e que a
unicidade não se veri ca no caso geral.
Mais tarde, Chipot e Lovat [18, 19] estudaram a existência e unicidade de soluções para
11
Método dos Elementos Finitos para problemas com Fronteiras Livres
Em 2004, Corrêa, Menezes e Ferreira [22] estenderam os resultados obtidos por Chipot e
Lovat, considerando a = a(l(u)) e f = f (x, u) funções contínuas. De facto, em [22], os autores
melhoraram os resultados dos trabalhos [20, 18, 19] considerando tanto o caso estacionário
como o caso evolutivo onde a não linearidade aparece, não só no operador u → a(l(u))∆u, mas
também na função não linear f do segundo membro.
No caso dos sistemas, a maioria dos autores assume que a matriz a de difusão é diagonal,
cando o acoplamento entre as equações apenas no termo reativo f . No entanto, a difusão
cruzada pode existir, veja-se por exemplo Oliveira [25] e as respectivas referências. Oliveira
considerou um sistema de reação-difusão onde a é uma matriz real ne × ne não diagonizável e
f: Rne → Rne é uma função de classe C 2.
Raposo et al. [62], em 2008, estudaram o sistema de reação-difusão da forma
(
ut − a(l(u))∆u + f (u − v) = α(u − v) em Ω×]0, T ]
,
vt − a(l(v))∆v − f (u − v) = α(v − u) em Ω×]0, T ]
com a(l) > 0, f funções Lipschitz-contínuas, l uma forma linear contínua e α uma constante
positiva e provaram a existência, unicidade e decaimento exponencial das soluções fracas.
Há poucos estudos sobre a aproximação numérica de problemas não locais e os que existem
são bastante restritos a condições de fronteira não locais ou termos de reação não locais.
Yin e Xu [76] aplicaram o método dos volumes nitos para aproximar soluções de um pro-
blema não local em escoamentos reativos em meios porosos e deduziram ordens de convergên-
cia ótimas na norma L2 .
Sidi Ammi e Torres [70] propuseram a aplicação do método dos elementos nitos no espaço
e os métodos de Euler ou Crank Nickolson no tempo para uma discretização completa de um
problema não local resultante do estudo da temperatura num termistor. Os autores provaram
ordens de convergência ótimas na norma L2 .
12
Método dos Elementos Finitos para problemas com Fronteiras Livres
com a(l) > 0 para todo l 6= 0, a(0) ≥ 0 e f Lipschitz-contínua satisfazendo f (0) = 0 e provaram
a existência e unicidade de soluções fracas. Deram também condições em u0 para a extinção
em tempo nito ou persistência no tempo e propuseram um esquema de diferenças nitas para
estudarem o comportamento das soluções.
Em [62], os autores também realizaram simulações numéricas para comparar o comporta-
mento das soluções do problema com coe ciente de difusão constante, com o comportamento
das soluções do problema não local, e este último com o de Ackleh e Ke [1]. Utilizaram um es-
quema implícito de diferenças nitas em dimensão um e de volumes nitos ([37]) em dimensão
dois, nas variáveis espaciais.
Até à data ainda não foi feito um estudo da convergência do método dos elementos nitos
de Euler-Galerkin quando aplicado ao sistema (1.11). Uma vez que o problema ainda não foi es-
tudado por outros autores, também não foram provadas a existência, unicidade e propriedades
físicas das soluções globais deste tipo de sistemas.
Os resultados apresentados neste trabalho estendem os de Chipot e Lovart, Corrêa, Menezes
e Ferreira, Raposo et al. e Simsen e Ferreira para sistemas não locais e provam novos resultados
sobre as propriedades das soluções.
Este trabalho está organizado da seguinte forma: no capítulo 2, são de nidos alguns con-
ceitos básicos e notações e enunciados alguns resultados conhecidos que irão ser utilizados nos
capítulos seguintes. É também relembrado o método dos elementos nitos de Galerkin e as suas
principais propriedades. No capítulo 3, estuda-se a equação em meios porosos com expoente
variável. Em domínio xo, de ne-se um problema regularizado introduzindo um parâmetro ε.
Demonstram-se estimativas para a solução do problema regularizado em função de ε e prova-se
a convergência deste para o problema inicial quando ε tende para zero. Neste capítulo, tam-
bém é deduzida a ordem de convergência do método dos elementos nitos de Galerkin quando
aplicado ao problema regularizado. São ainda apresentados exemplos de implementação, em
Matlab, do método. No domínio livre, o problema é abordado numa nova formulação onde
se deduz uma equação para o movimento da fronteira. Para a malha evoluir com o domínio,
constrói-se uma malha adaptável e discretizam-se a equação do problema e a equação do mo-
vimento da fronteira nessa malha. Para testar o método, são apresentados exemplos da sua
aplicação a problemas concretos. O sistema não linear com termo de difusão não local é tra-
tado no capítulo 4 onde se formula o problema e as hipóteses sobre os dados do problema; se
prova a existência e unicidade de soluções fortes globais no tempo; se estudam também as pro-
priedades de localização e comportamento assintótico das soluções; se discretiza o problema
no espaço e tempo e se provam estimativas do erro e ordem de convergência para as soluções
discretas. O capítulo termina com a apresentação de exemplos de implementação do método
em Matlab.
Alguns dos resultados obtidos, no estudo feito para a elaboração desta tese, foram apre-
sentados pelo autor em vários seminários e workshops na Universidade da Beira Interior e no
Centro de Matemática e Aplicações Fundamentais da Universidade de Lisboa. O autor partici-
pou também nos congressos internacionais: International Conference on Engineering UBI2011
que teve lugar na Covilhã, Portugal, em 2011; 6th Workshop on Statistics, Mathematics and
13
Método dos Elementos Finitos para problemas com Fronteiras Livres
14
Método dos Elementos Finitos para problemas com Fronteiras Livres
Capítulo 2
a2 b2
ab ≤ + .
2 2
Lema 2.2 (Young) Sejam a, b ≥ 0, 1 < p, p0 < +∞ e > 0 números reais, então
0 0
ap 1−p bp 1 1
ab ≤ + , + = 1.
p p0 p p0
Demonstração
Suponha-se que a > b, então
a−b n Z a−b
Z a Z
η + a + b dη ≥ (n + 1)
2 2
|a|n a − |b|n b = (n + 1) |ξ|n dξ = (n + 1)
|η|n dη =
b − a−b
2
2 a−b
− 2
a−b
a − b n a − b
Z
2 n+1 n
= 2(n + 1) |η|n dη = 2(n + 1) = |a − b| (a − b).
0 2 2 2n
15
Método dos Elementos Finitos para problemas com Fronteiras Livres
n+1 n
|a|n a − |b|n b ≤ |a − b| (a − b).
2n
Lema 2.4 Sejam a, b ∈ R e p ≥ 2, então existe uma constante C > 0, apenas dependente de p,
tal que
p−2 p−2 p−2
|x| x − |y| y ≤ C (|x| + |y|) |x − y| .
Demonstração
Para p ≥ 2, a função g(z) = |z|p−2 z é contínua e diferenciável em R, com g0(z) = (p −
1)|z|p−2 . Pelo teorema do valor médio de Lagrange, existe ϑ entre a e b tal que
16
Método dos Elementos Finitos para problemas com Fronteiras Livres
O conjunto M é fechado se M = M̄ .
∂M = M̄ ∩ (X \ M )
é chamado a fronteira de M .
De nição 2.11 (Conjunto conexo) Um conjunto aberto Ω ∈ Rd é dito conexo se dois quaisquer
pontos de Ω poderem ser unidos por uma linha poligonal totalmente contida em Ω.
De nição 2.15 (Derivada forte (clássica)) Seja f : Rn → R uma função contínua em x, então
de ne-se derivada parcial de f em relação à coordenada i:
∂f f (x + hei ) − f (x)
(x) = lim .
∂xi h→0 h
De nição 2.16 (C k (Ω)) Para k ∈ N e Ω ⊂ Rn, C k (Ω) representa o conjunto de todas as funções
de nidas em Ω e que têm todas as derivadas até à ordem k contínuas em Ω.
De nição 2.17 (Conjunto convexo) O conjunto M é dito convexo se dois quaisquer pontos do
conjunto poderem ser unidos por um segmento totalmente contido no conjunto, ou seja,
a, b ∈ M, λ ∈ [0, 1] ⇒ λa + (1 − λ)b ∈ M.
De nição 2.18 (Função convexa) Seja Ω um conjunto convexo, uma função f de nida em Ω é
dita convexa se
17
Método dos Elementos Finitos para problemas com Fronteiras Livres
Um funcional linear contínuo de nido num espaço linear normado X é uma aplicação linear
contínua f : X → X. Para denotar o valor de f em u ∈ X utiliza-se a notação f (u) ou hf, ui.
|hf, ui|
kf kX ∗ = sup , f ∈ X ∗,
u∈X kukX
u6=0
é um espaço de Banach.
De nição 2.22 (Convergência fraca) Seja X um espaço linear normado e {un } ∈ X, diz-se
que {un } converge fracamente para u ∈ X se
Observação 2.26 Um operador compacto também pode ser denominado operador linear com-
pletamente contínuo.
Observação 2.29 Por vezes, refere-se um domínio com fronteira Lipschitz-contínua apenas
por domínio Lipschitz.
Se gk ∈ C ϑ , para todo k = 1, . . . , k, então diz-se que ∂Ω ∈ C ϑ .
18
Método dos Elementos Finitos para problemas com Fronteiras Livres
2.3 Distribuições
é um conjunto compacto de Ω.
Observação 2.34 As derivadas distributivas, até à ordem k, de uma função f ∈ C k (Ω̄) são
iguais às derivadas clássicas de f , por isso denotar-se-á por
∂ |α| f
,
∂xα
1
1
. . . xα
n
n
O estudo costuma ser feito num conjunto Ω(t) e num intervalo de tempo ]t0 , T ], T > t0 . A
função u está então de nida em
∂u ∂u
∇u = ( ,..., )
∂x1 ∂xd
d
X ∂2f
∆u =
i=1
∂x2i
19
Método dos Elementos Finitos para problemas com Fronteiras Livres
∂2f
denota o laplaciano de f , onde ∂x2i
é a derivada de segunda ordem de f em relação a xi .
Se f = (f1 , . . . , fd ), então o divergente de f é de nido por
d
X ∂fi
div f = .
i=1
∂xi
De nição 2.35 (Espaço Lp ) Para p ∈ [1, +∞[, Lp (A) é o espaço linear de todas as ( classes de
equivalência das ) funções u mensuráveis em A para as quais
Z
|u|p dx < ∞.
A
Z p1
p
kukLp (A) = |u(x)| dx , 1 ≤ p < ∞,
A
é um espaço de Banach.
De nição 2.37 (Espaço L∞ ) L∞ (A) é o espaço linear de todas as (classes de equivalência das)
funções u mensuráveis em A para as quais:
é um espaço de Banach.
2. (Desigualdade de Minkowski)
3. (Desigualdade de Hölder)
Z
1 1
|u(x)v(x)| dx ≤ kukLp (X) kvkLp0 (X) , 1 ≤ p, p0 ≤ +∞, + = 1, ∀u ∈ Lp (X), v ∈ Lp0 (X).
G p p0
20
Método dos Elementos Finitos para problemas com Fronteiras Livres
1 1
kukLp (X) ≤ (meas(X)) p − q kukLq (X) , ∀u ∈ Lp (X).
6. Sejam 1 ≤ p ≤ r ≤ q ≤ +∞ e 1
r = α
p + q ,
1−α
0 ≤ α ≤ 1, e u ∈ Lp (X) ∩ Lq (X) então
u ∈ Lr (X) e
1−α
kukLr (X) ≤ kukα
Lp (X) kukLq (X) .
Lema 2.42 Seja A um aberto limitado de Rd, gk e g funções de Lq (A), 1 < q < ∞ tais que:
kgk kLq (A) ≤ C, gk → g, qtp em A,
então gk * g em Lq (A).
e
k
X
kukW k,∞ (Ω) = ess sup |Dα u|.
Ω
|α|=0
21
Método dos Elementos Finitos para problemas com Fronteiras Livres
Z k
X
(u, v) = Dα uDα v dx.
Ω |α|=0
Neste trabalho, interpreta-se W0k,p (Ω) como o conjunto das funções u ∈ W k,p (Ω) tais que
Dα u = 0 em ∂Ω para todo |α| ≤ k − 1. Denota-se também H0k (Ω) = W0k,2 (Ω)
Observação 2.46 |.|W k,p,k (Ω) é uma norma em W0k,p (Ω). Em particular, denota-se
Z 12
kukH01 (Ω) = |∇u|2 dx .
Ω
Observação 2.47
Lp (Ω) = W 0,p (Ω) ⊃ W 1,p (Ω) ⊃ W 2,p (Ω) ⊃ . . . .
0 0
De nição 2.48 (W −k,p (Ω)) Se 1 ≤ p < ∞, denota-se por W −k,p (Ω) o espaço dual de W0k,p (Ω),
1
p + 1
p0 = 1, e por H0−k (Ω) o dual de H k (Ω).
2.6 Imersões
De nição 2.49 Sejam X ⊂ Y dois espaços lineares normados, com as normas k · kX e k · kY
respetivamente. De ne-se o operador identidade por:
I:X → Y
.
u 7→ I(u) = u
De nição 2.50 (Imersão contínua) Se I é contínuo, então diz-se que é uma imersão contínua
de X em Y e escreve-se X ,→ Y.
De nição 2.52 (Imersão compacta) Se I é completamente contínuo, então diz-se que I é uma
imersão compacta e escreve-se X ,→,→ Y .
X ,→ Y ⇒ Y ∗ ,→ X ∗
e
X ,→,→ Y ⇒ Y ∗ ,→,→ X ∗ .
22
Método dos Elementos Finitos para problemas com Fronteiras Livres
Se d
p <k< d
p + 1, então W k,p (Ω) ,→ C 0,k−d/p (Ω̄).
Se k = d
p + 1, então W k,p (Ω) ,→ C 0,α (Ω̄) ∀α ∈]0, 1[.
Se k > d
p + 1, então W k,p (Ω) ,→ C 0,1 (Ω̄).
De nição 2.57 (Lp (a, b; X)) Seja X um espaço de Banach. Para 1 ≤ p ≤ ∞ de ne-se Lp (a, b; X)
como sendo o espaço das funções u : [a, b] → X, tais que:
! p1
Z b
kukLp (a,b;X) = ku(t)kpX dt < ∞ para 1 ≤ p < ∞
a
e
kukL∞ (a,b;X) = ess sup ku(t)kX < ∞.
t∈[a,b]
De nição 2.59 Sejam a, b ∈ R, a < b e X um espaço de Banach com norma k · kX . Diz-se que a
função u : [a, b] → X é contínua num ponto t0 ∈ [a, b] se
Denota-se por C 0 ([a, b], X) o conjunto de todas as funções contínuas no intervalo [a, b] e com
valores em X.
23
Método dos Elementos Finitos para problemas com Fronteiras Livres
é um espaço de Banach.
De nição 2.61 Diz-se que a função u : [a, b] → X é diferenciável no ponto t0 ∈ [a, b] se existir
w ∈ X tal que:
u(t0 + h) − u(t0 )
lim k − wkX = 0.
h→0 h
Neste caso, du
dt (t0 ) = ut (t0 ) é chamada derivada forte de u em t0 .
De nição 2.62 Uma função f ∈ L1 (a, b; X) tem derivada generalizada se existir uma função
g ∈ L1 (a, b, X) tal que:
Z b Z b
dϕ
(t)f (t) dt = − ϕ(t)g(t) dt, ∀ϕ ∈ D(a, b).
a dt a
De nição 2.63 (W k,p (a, b; X)) Seja X um espaço de Banach. Para k = 1, 2, . . . e p ∈ [1, ∞],
de ne-se
dj f
W k,p (a, b; X) = {f ∈ Lp (a, b; X) : ∈ Lp (a, b; X), j = 1, . . . , k}.
dtj
Estes espaços têm as mesmas propriedades que foram enunciadas para os espaços Lp .
De uma forma semelhante, podem de nir-se os espaços C k ([a, b]; X) das funções contínuas
e diferenciáveis, em t, até à ordem k, no intervalo [a, b] e com valores em X.
De nição 2.64
dj f
C k ([a, b]; X) = {f ∈ C 0 ([a, b]; X) : ∈ C 0 ([a, b]; X), j = 1, . . . , k}.
dtj
Teorema 2.65 Se u ∈ L2 (0, T ; H01 (Ω)) e ut ∈ L2 (0, T ; H −1 (Ω)) então u ∈ C 0 ([0, T ]; L2 (Ω)).
Nesta secção enunciam-se alguns teoremas que são necessários nas demonstrações, mas não
se enquadram nas secções anteriores.
u ∈ X tais que:
unk * u.
24
Método dos Elementos Finitos para problemas com Fronteiras Livres
f ∈ X ∗ tais que:
∗
fnk * f.
Teorema 2.70 (Schauder) Seja X um espaço de Banach, M ⊂ X um conjunto não vazio, limi-
tado, fechado e convexo e seja F : M → M um operador compacto. Então existe, pelo menos,
um ponto xo u ∈ M de F .
Teorema 2.71 (Green) Seja Ω um conjunto conexo limitado com fronteira Lipschitz-contínua
e sejam u, v ∈ C 2 (Ω̄). Então
Z Z Z
v∆u dx = v∇u · n ds − ∇v · ∇u dx,
Ω ∂Ω Ω
Forma diferencial Suponha-se que h e r são integráveis em ]a, b[ e não negativas qtp em ]a, b[.
Se y ∈ C([a, b]), yt ∈ L1 (]a, b[) e a desigualdade seguinte é satisfeita:
então: Z t Rs
R
t
− r(τ ) dτ
y(t) ≤ y(a) + h(s)e a ds e a r(s) ds
.
a
Forma integral Suponha-se que h é contínua em [a, b], r integrável em ]a, b[ e h, r ≥ 0 q.t.p.
em ]a, b[. Se y é uma função contínua em [a, b] e satisfaz a inequação:
Z t
y(t) ≤ h(t) + r(s)y(s) ds, t ∈ [a, b],
a
então Z t Rt
r(τ ) dτ
y(t) ≤ h(t) + h(s)r(s)e s ds, t ∈ [a, b].
a
Lema 2.73 (Aubin-Lions) Sejam B0 , B e B1 três espaços de Banach tais que B0 e B1 são re e-
xivos de modo que B0 ,→,→ B ,→ B1 . Seja
25
Método dos Elementos Finitos para problemas com Fronteiras Livres
• Ω = ∪nt
k=1 Tk ;
• int(Tk ) 6= ∅, ∀Tk ∈ Th ;
• int(Ti ) ∩ int(Tj ) = ∅, ∀Ti , Tj ∈ Th , j 6= i;
• cada lado de Tk ou pertence à fronteira de Ω ou é lado de outro Ti ∈ Th ;
• cada Tk tem fronteira Lipchitz-contínua .
h = max{hk , k = 1, . . . , nt}.
De nição 2.74 Uma família Th é dita regular se forem satisfeitas simultaneamente as duas
condições:
2. h → 0.
Uma vez estabelecida a malha Th , de ne-se o espaço dos elementos nitos Shr .
Primeiro, de ne-se um espaço linear de dimensão nita Qkr da forma
Neste trabalho, considera-se Qkr o conjunto dos polinómios de grau menor ou igual a r em x e
y e exige-se que vh seja contínua em Ω e se anule na fronteira de Ω. Resumindo, de ne-se
Este teorema permite que, na resolução dos problemas, se possa considerar Shr , pois Shr ⊂ S.
A inclusão anterior não é veri cada no caso da fronteira de Ω ser uma curva que não pode ser
exatamente descrita por triângulos.
Finalmente, tem de se de encontrar um conjunto Σk de condições linearmente independen-
tes que permitam determinar, de forma unívoca, cada elemento de Qkr . Um polinómio de grau
menor ou igual a r em x e y ca perfeitamente de nido pelos valores em nk = (r + 1)(r + 2)/2
pontos distintos. Sejam a1 , . . . , ank pontos distintos igualmente espaçados em Tk , então pode-se
26
Método dos Elementos Finitos para problemas com Fronteiras Livres
O espaço Shr descrito desta forma tem uma base canónica ϕ1 , . . . , ϕnk de nida por ϕj (ai ) = δij ,
sendo δij a função delta de Kronecker. Assim, cada elemento nito é de nido pelo terno
(Tk , Qkr , Σkr ).
Suponha-se que (Tk , Qkr , Σkr ) é um elemento nito onde Tk é um triângulo de vértices a1 ,
a2 e a3 e seja T̂k um outro triângulo de vértices â1 , â2 e â3 . Existe uma única aplicação a m:
Φk : T̂k → Tk
x̂ → x = Bk x̂ + bk ,
De nição 2.76 Dois elementos nitos (Tk , Qkr , Σkr ) e (T̂k , Q̂kr , Σ̂kr ) dizem-se a m-equivalen-
tes se existir uma aplicação a m invertível Φk , tal que:
Tk = Φk (T̂k ) (2.2)
Qkr = {pr : Tk → R| pr = p̂r ◦ Φ−1
k , p̂r ∈ Q̂kr } (2.3)
Σkr = {pr (Φk (âi ))}. (2.4)
Uma família Th de elementos nitos é dita família a m se todos os elementos nitos que
a constituem são a m-equivalentes a um único elemento nito, que se irá chamar elemento
nito de referência (TR , QRr , ΣRr ).
Proposição 2.77 Qualquer elemento de uma família a m Th , ou da sua imagem por uma apli-
cação a m invertível, pode ser escolhido como elemento de referência.
Pela proposição, pode veri car-se que TR não tem de pertencer à família e que até pode
ser escolhido de uma forma bastante livre. Neste trabalho, considera-se
TR = {(x, y) ∈ R2 | x + y ≤ 1 ∧ x ≥ 0 ∧ y ≥ 0}.
Seja ΠPr (Tk ) o operador interpolador que a cada função v de nida em Tk faz corresponder o
polinómio, pr , interpolador de Lagrange, de grau menor ou igual a r, que interpola v nos pontos
ai , i = 1, . . . , nk. É óbvio que pr ∈ Qkr e pela de nição da base ϕ1 , . . . , ϕnk pode de nir-se
ΠPr (Tk ) por
nk
ΠPr (Tk ) (v) =
X
ϕi v(ai ).
i=1
Sejam Ph e Eh o conjunto dos vértices e dos lados dos triângulos, respetivamente. Retirando
as repetições e introduzindo uma ordenação em Ph e Eh ca-se com Ph = {P1 , . . . , Pnv } e
27
Método dos Elementos Finitos para problemas com Fronteiras Livres
Eh = {e1 , . . . , enl }.
Se ek ∈ ∂Ω, denota-se por ~nk o vetor normal exterior a ek .
De nição 2.78 Seja v ∈ L2 (Ω). ṽ ∈ Shr é a projeção ortogonal de v sobre Shr denotado por
ṽ = ΠShr v se satisfaz a condição
Z Z
ṽχ dx = vχ dx, ∀χ ∈ Shr .
Ω Ω
para 0 ≤ m ≤ r e 1
p + 1
p0 = 1.
De nição 2.80 (Ordem de convergência) Suponha-se que limh→0+ vh = v. Diz-se que a ordem
de convergência é p > 0 se
kvh − vkX ≤ Chp , ∀h > 0, (2.5)
onde C é uma constante positiva que não depende de h. No caso de limn→∞ vn = v, diz-se que
a ordem de convergência é p > 0 se
Quanto maior for p, maior é a velocidade com que a aproximação converge para o valor exato.
Aplicando o logaritmo em (2.5), obtém-se
o que signi ca que se se representa gra camente log(kvh − vkX ) em função de log(h), o grá co
deverá estar abaixo da reta que passa por (0, log(C)) e tem declive p. Este processo revela-se
útil na estimação numérica da ordem de convergência. Para (2.6) o raciocínio é o mesmo.
De nição 2.82 No espaço Shr de ne-se o laplaciano discreto ∆h : Shr → Shr da forma
Z Z
(∆h u)χ dx = − ∇u · ∇χ dx, ∀u, χ ∈ Shr .
Ω Ω
De nição 2.83 (Projeção de Ritz com peso) A projeção de Ritz com peso ρh = Rh ρ ∈ Shr é
de nida por Z Z
A∇ρh · ∇χ dx = A∇ρ · ∇χ dx, ∀χ ∈ Shr .
Ω Ω
28
Método dos Elementos Finitos para problemas com Fronteiras Livres
ni−1
[
]0, T ] = In , In =]tn , tn+1 ], tn+1 = tn + δ.
n=0
No caso da solução exata ter pouca regularidade em t, um método mais adequado é o método
dos elementos nitos descontínuo que aproxima a solução em cada intervalo de tempo por um
polinómio. Para isso, de ne-se o espaço das funções seccionalmente polinomiais em relação ao
tempo com coe cientes em Shr como sendo
s
X
S δs = {W : [0, +∞[→ H | W|In = tj χj , χj ∈ H}.
j=0
As funções de S δs podem ser descontínuas nos instantes tn , mas são contínuas à esquerda de
tn . Denota-se Wn = W (tn ) e Wn+ = limt→t+
n
W (t).
Para deduzir estimativas de erros, utilizam-se algumas projeções. De ne-se a projeção no
tempo ΠS δs . ṽ = ΠS δs v é obtido pelas condições:
ṽ(tn+1 ) = v(tn+1 ), n = 0, . . . , ni − 1,
Z
(ṽ − v, W ) dt = 0, ∀W ∈ S δs−1 , n = 0, . . . , ni − 1.
In
Para provar que está bem de nido, ver [72] na página 207. No caso H = Shr obtém-se o espaço
δs
Shr .
De ne-se também a projeção no espaço e tempo ΠShr
δs . Ṽ = ΠS δs v é solução de
hr
Z Z
Ṽ χ dx = ṽχ dx, ∀χ ∈ Shr , ∀t ∈ In ,
Ω Ω
Z Z
(Ṽ − v, W ) dt = (ṽ, W ) − (v, W ) dt = 0, δs−1
∀W ∈ Shr
In In
29
Método dos Elementos Finitos para problemas com Fronteiras Livres
e
(Ṽ (tn ) − v(tn ), χ) = (ṽ(tn ), χ) − (v(tn ), χ) = 0, ∀χ ∈ Shr .
e
|v(t) − ṽ(t)|H 1 (Ω) ≤ Cδ s+1 |v|W s+1,∞ (In ;H 1 (Ω)) , ∀t ∈ In .
Demonstração
Pelo teorema 12.1 em [72], tem-se que
Z
|v(t) − ṽ(t)|2H k (Ω) ≤ Cδ 2s+1 |v (s+1) (t)|2H k (Ω) dt, k = 0, 1.
In
Utilizando a desigualdade
Z
kv (s+1) (t)k2H k (Ω) dt ≤ δ|v|2W s+1,∞ (In ;H k (Ω)) ,
In
obtém-se o pretendido.
Teorema 2.86 Seja v ∈ L∞ (In ; H r+1 (Ω)) ∩ W s+1,∞ (In ; L2 (Ω)) e Ṽ a projeção de v em Shr
δs
,
então:
kv(t) − Ṽ (t)kL2 (Ω) ≤ Chr+1 |v|L∞ (In ;H r+1 (Ω)) + Cδ s+1 |v|W s+1,∞ (In ;L2 (Ω)) .
Demonstração
Pela desigualdade triangular, obtém-se
kv(t) − Ṽ (t)kL2 (Ω) ≤ kv(t) − ΠShr v(t)k + kΠShr v(t) − ΠShr ṽ(t)k
≤ Chr+1 |v|L∞ (In ;H r+1 (Ω)) + Cδ s+1 |v|W s+1,∞ (In ;L2 (Ω)) .
Para problemas mais complicados, existem integradores mais robustos, por exemplo, o
ode15s [66] que é uma subrotina do Matlab e que aproxima numericamente a solução de equa-
ções e sistemas EDO da forma M (t, y)y 0 = F (t, y). Este integrador utiliza um método de ordem
variável (de 1 a 5) baseado nas fórmulas de diferenciação numérica (NDFs). Opcionalmente,
utiliza o método de Gear (BDFs). Começa no instante inicial e calcula a aproximação no instante
seguinte, usando um determinado passo. Se o erro da aproximação satis zer a tolerância dese-
jada, então é considerado um sucesso e passa para o instante seguinte, com o mesmo passo. Se
não satis zer a tolerância, é considerado um passo falhado e é calculada uma nova estimativa
com o tamanho do passo reduzido a metade. Este processo continua desta forma até conseguir
um sucesso ou o tamanho do passo ser muito pequeno.
30
Método dos Elementos Finitos para problemas com Fronteiras Livres
Nesta secção, calculam-se explicitamente as funções ϕ1 , . . . , ϕnk . Uma vez que as apro-
ximações pretendidas são polinómios de grau superior a 1 e os vértices do triângulos podem
mover-se, é vantajoso de nir as funções ϕ em cada Tk utilizando coordenadas triangulares.
Para se obter um polinómio interpolador de grau r em cada triângulo, é necessário de nir
nk pontos nesse triângulo. Segundo [9], estes pontos são de nidos criando uma malha uniforme
como mostra a gura 2.1. Seja Tk ∈ T um elemento qualquer. Considere-se que as coordenadas
X2 Y3 − X3 Y2 + (Y2 − Y3 )x + (X3 − X2 )y
ψ1 (x, y) = ,
24
X3 Y1 − X1 Y3 + (Y3 − Y1 )x + (X1 − X3 )y
ψ2 (x, y) = ,
24
X1 Y2 − X2 Y1 + (Y1 − Y2 )x + (X2 − X1 )y
ψ3 (x, y) = ,
24
onde
X1 (Y2 − Y3 ) + X2 (Y3 − Y1 ) + X3 (Y1 − Y2 )
4= .
2
Por vezes, às funções ψ chamam-se funções pirâmide, pois ψj é linear e toma valor 1 em
(Xj , Yj ) e zero nos outros pontos. Para a aproximação com polinómios de grau r, o polinómio
interpolador associado ao ponto Pi é construído à custa das funções interpoladoras de Lagrange
de grau 1 da forma:
ϕi (x, y) = f¯a (ψ1 (x, y))f¯b (ψ2 (x, y))f¯c (ψ3 (x, y)),
31
Método dos Elementos Finitos para problemas com Fronteiras Livres
onde k
Y1
(rψj − i + 1), 1 ≤ k ≤ r
f¯k (ψj ) = i=1
i
1, k=0
com
k k
X r Y 1
(rψj − i + 1), 2 ≤ k ≤ r
∂ f¯k
m=1
m i=1 i
= i6=m
∂ψj
r, k=1
0, k=0
e
∂ψ1 Y2 − Y3 ∂ψ2 Y3 − Y1 ∂ψ3 Y1 − Y2
= , = , = ,
∂x 24 ∂x 24 ∂x 24
∂ψ1 X3 − X2 ∂ψ2 X1 − X3 ∂ψ3 X2 − X1
= , = , = .
∂y 24 ∂y 24 ∂y 24
Embora os integrais envolvendo apenas as funções interpoladoras possam ser calculados exa-
tamente, não é prático o cálculo exato dos integrais que envolvem outras funções do problema,
possivelmente não lineares.
Neste trabalho, todos os integrais são calculados numericamente utilizando a quadratura
de Gauss-Legendre [56]. Os pontos e os pesos de quadratura são calculados num triângulo de
referência TR e considera-se Φk : TR → Tk onde cada elemento Tk da partição é imagem
do triângulo de referência TR . Sejam (α, β) as coordenadas locais no triângulo Tk e (x, y) as
coordenadas locais do triângulo TR , então
Z Z
I= h(α, β) dα dβ = (h ◦ Φk )(x, y)|JΦk | dx dy.
Tk TR
Se se considerar TR o triângulo de vértices (0, 0), (1, 0) e (0, 1), então pretende-se calcular
Z Z 1 Z 1−x
I= f (x, y) dx dy = f (x, y) dy dx.
TR 0 0
32
Método dos Elementos Finitos para problemas com Fronteiras Livres
1 1
1 + ξ (1 − ξ)(1 − η) 1−ξ
Z Z
I= f , dξ dη.
−1 −1 2 4 8
ni X
ni
X 1 − ξi 1 + ξi (1 − ξi )(1 − ηj )
I≈ wi wj f , ,
i=1 j=1
8 2 4
Logo, " #
X2 − X1 X3 − X1
JΦk = e |JΦk | = 24.
Y2 − Y1 Y3 − Y1
(1−ξi )(1−ηj )
Sejam W xi = 1−ξi
8 wi , W yj = wj , Cxij = 1−ξi
2 e Cyij = 4 , i, j = 1 . . . m, então
m X
X m
I= W xi W yi (h ◦ Φk )(Cxij , Cyi,j ) 24.
i=1 j=1
33
Método dos Elementos Finitos para problemas com Fronteiras Livres
34
Método dos Elementos Finitos para problemas com Fronteiras Livres
Capítulo 3
De nição 3.1 (Solução fraca) Uma função u(x, t) é dita solução fraca do problema (3.1) se
γ
(i) u ∈ L∞ (0, T ; L∞ (Ω)), |u| 2 ∇u ∈ L2 (0, T ; L2 (Ω)), ut ∈ L2 (0, T ; W −1,2 (Ω)),
(ii) u = 0 em ΓT ,
Em [4], Antontsev e Shmarev provaram que o problema (3.1) tem solução fraca segundo a
de nição 3.1.
35
Método dos Elementos Finitos para problemas com Fronteiras Livres
+ γ(x)
εγ ≤ aε (x, v) = (v 2 + ε2 ) 2 , 0 < ε < 1. (3.4)
De nição 3.2 (Solução fraca) Uma função v(x, t) é dita solução fraca do problema (3.3) se
1
(i) v ∈ L∞ (0, T ; L∞ (Ω)), aε (x, v) 2 ∇v ∈ L2 (0, T ; L2 (Ω)), vt ∈ L2 (0, T ; W −1,2 (Ω)),
(ii) v = 0 em ΓT ,
Segundo a teoria clássica, para cada ε > 0, o problema (3.3) tem solução fraca que satisfaz
esta de nição e o sistema (3.3). Além disso, se ∇γ ∈ L2 (ΩT ), então vt , ∆v ∈ L2 (ΩT ).
Uma vez que a in uência do parâmetro ε na regularidade das soluções do problema regularizado
não é evidente, demonstram-se alguns resultados sobre a regularidade das soluções.
Começa-se por obter estimativas para as soluções do problema aproximado.
Teorema 3.3 Seja γ uma função mensurável em Ω que satisfaz a condição (3.2). Se
Z T
ku0 kL∞ (Ω) + kf kL∞ (Ω) dt < C, (3.5)
0
Demonstração
Multiplicando (3.3) por v 2k−1 e integrando em Ω, obtém-se
Z Z Z
vt v 2k−1 dx + div(aε (x, v)∇v)v 2k−1 dx = f v 2k−1 dx.
Ω Ω Ω
36
Método dos Elementos Finitos para problemas com Fronteiras Livres
O que resulta em
Z
1 d
kvk2k
L2k (Ω) + (2k − 1) v 2k−2 aε (x, v)|∇v|2 dx ≤ kf kL2k (Ω) kvk2k−1
L2k (Ω) .
2k dt Ω
d
kvkL2k (Ω) ≤ kf kL2k (Ω) .
dt
Integrando em t ca
Z T
kv(x, t)kL2k (Ω) ≤ kv(x, 0)kL2k (Ω) + kf kL2k (Ω) dt.
0
Lema 3.4 Se γ satisfaz (3.2) e a é de nido por (3.4) com v solução de (3.3), então
aε (x, v) ≤ C,
Demonstração
Segue imediatamente de (3.4) e de (3.6).
Teorema 3.5 Suponha-se que v é solução de (3.3), com γ ∈ H 1 (Ω) e a condição (3.5) é válida,
então Z T Z
+
|∇v|2 dxdt ≤ Cε−γ ,
0 Ω
Demonstração
Da mesma forma como a demonstração anterior, multiplicando (3.3) por v e integrando em
Ω, obtém-se Z Z Z
vt v dx + div(aε (x, v)∇v)v dx = f v dx.
Ω Ω Ω
O que resulta em
Z
1 d
kvk2L2 (Ω) + aε (x, v)|∇v|2 dx ≤ kf kL2 (Ω) kvkL2 (Ω) ,
2 dt Ω
37
Método dos Elementos Finitos para problemas com Fronteiras Livres
Atendendo a (3.4), ca
Z T Z
+
|∇v|2 dxdt ≤ Cε−γ .
0 Ω
Lema 3.6 Se γ ∈ H 1 (Ω) satisfaz (3.2) e a é de nido por (3.4) com v solução de (3.3), então
+
k∇aε (x, v)kL2 (Ω) < Cε−γ ,
Demonstração
Segue imediatamente do teorema anterior, de (3.6) e de
γ γ ∇γ
∇aε (x, v) = γv(v 2 + ε2 ) 2 −1 ∇v + (v 2 + ε2 ) 2 ln(v 2 + ε2 ).
2
Teorema 3.7 Suponha-se que u0 ∈ H 1 (Ω), γ ∈ H 2 (Ω) e a condição (3.5) é válida. Se v é solução
de (3.3), então
Z T Z
+
(∆v)2 dxdt ≤ Cε−8γ ,
0 Ω
Demonstração
Considere-se a seguinte função:
Z v(x,t) Z v
γ(x)
w(x, t) = (ς 2 + ε2 ) 2 dς = aε (x, ς) dς,
0 0
ou seja, Z v
1
vt = wt e aε (x, v)∇v = ∇w − ∇aε (x, ς) dς.
aε (x, v) 0
ou de forma equivalente,
wt = aε (x, v)∆w − aε (x, v)g1 ,
com
v v
|∇γ|2
Z Z
∆γ
g1 = (∇aε (x, ς)|ς=v ) · ∇v + aε (x, ς) ln(ξ 2 + ε2 ) dς + aε (x, ς) ln2 (ς 2 + ε2 ) dς − f.
0 2 0 4
38
Método dos Elementos Finitos para problemas com Fronteiras Livres
Em particular,
Z Z
−γ + +
|∇w(x, T )| dx ≤ Cε 2
e (∆w)2 dxdt ≤ Cε−2γ .
Ω ΩT
Por outro lado ([53], lema 4.5, pág. 63) tem-se que:
Z Z
+
|∇w|4 dxdt ≤ C ess sup |w|2 (∆w)2 dxdt ≤ Cε−2γ .
ΩT (x,t)∈ΩT ΩT
então
+
|∇v| ≤ Cε−γ (|∇w| + 1).
γ
∆w = aε (x, v)∆v + aε (x, v) ln(v 2 + ε2 )∇γ · ∇v + γ(v 2 + ε2 ) 2 −1 v|∇v|2 + g2 ,
onde,
v v
|∇γ|2
Z Z
∆γ
g2 = aε (x, ς) ln(ς 2 + ε2 ) dς + aε (x, ς) ln2 (ς 2 + ε2 ) dς.
0 2 0 4
Finalmente, Z
+
|∇v|4 dxdt ≤ Cε−6γ
ΩT
e
Z Z Z 2
γ
(aε (x, v)∆v)2 dxdt ≤ (∆w)2 dxdt + γ(v 2 + ε2 ) 2 −1 v|∇v|2 dxdt+
ΩT ΩT ΩT
Z Z
+ a2ε (x, v) ln2 (v 2 2
+ ε )(∇γ · ∇v) dxdt + 2
g22 dxdt
ΩT ΩT
+ + + +
≤ Cε−2γ + Cε−6γ + Cε−γ + C ≤ Cε−6γ .
39
Método dos Elementos Finitos para problemas com Fronteiras Livres
Logo, Z
+
(∆v)2 dxdt ≤ Cε−8γ .
ΩT
Teorema 3.8 Suponha-se que γ ∈ H 2 (Ω), u0 ∈ H 1 (Ω), v é solução de (3.3) e a condição (3.5) é
válida, então
Z T Z
+
vt2 dxdt ≤ Cε−2γ ,
0 Ω
Demonstração
Da mesma forma que o teorema 3.7 seja
Z v
w= aε (x, ς) dς,
0
então Z v
∇w = aε (x, v)∇v + ∇aε (x, ς) dς
0
e
wt = aε (x, v)vt .
ou seja,
Z Z Z Z Z v
1
(wt )2 dx + ∇w∇wt dx = wt f dx + ∇aε (x, ς) dς∇wt dx.
Ω aε (x, v) Ω Ω Ω 0
onde: v v
|∇γ|2
Z Z
∆γ
g= aε (x, ς) ln(ς 2 + ε2 ) dς + aε (x, ς) ln2 (ς 2 + ε2 ) dς.
0 2 0 4
Integrando esta equação em t resulta
Z Z Z
1 1 1
(wt )2 dxdt + |∇w(x, T )| dx = 2
|∇w(x, 0)|2 dx+
ΩT aε (x, v) 2 Ω 2 Ω
Z Z Z
+ f wt dxdt − wt (∇aε (x, ς)|ς=v ) · ∇v dxdt − wt g dxdt.
ΩT ΩT ΩT
Z Z
+C aε (x, v)|(∇aε (x, ς)|ς=v ) · ∇v|2 dxdt + C aε (x, v)g 2 dxdt.
ΩT ΩT
40
Método dos Elementos Finitos para problemas com Fronteiras Livres
Como
1
|vt | = |wt |,
aε (x, v)
segue-se que Z Z
+ 1 +
(vt )2 dxdt ≤ ε−γ (wt )2 dxdt ≤ ε−2γ
ΩT ΩT aε (x, v)
e o resultado está provado.
Lema 3.9 Se γ satisfaz (3.2) e a é de nido por (3.4) com v solução de (3.3), então
+
k(aε (x, v))t kL2 (Ω) < Cε−γ ,
Demonstração
Segue imediatamente do teorema anterior e de:
γ
(aε (x, v))t = γvt (v 2 + ε2 ) 2 −1 .
Lema 3.10 Se γ satisfaz (3.2) e v é limitado, então para 0 < ε < 1 a desigualdade seguinte é
válida:
γ(x)
||v|γ(x) − (v 2 + ε2 ) 2 | ≤ Cε,
Demonstração
γ
Aplica-se o teorema do valor médio de Lagrange à função g(y) = (v 2 + y 2 ) 2 no intervalo
[0, ε] e conclui-se que, para cada x ∈ Ω, existe c ∈ [0, ε] ⊂ [0, 1], tal que:
γ γ−1
g(ε) − g(0) = γc(v 2 + c2 ) 2 −1 ε ≤ γ(v 2 + c2 ) 2 ε,
logo
γ(x)
||v|γ(x) − (v 2 + ε2 ) 2 | ≤ Cε.
Para demonstrar a convergência, seguem-se as ideias da demonstração do teorema de unici-
dade em [4] e suas referências. Primeiro, prova-se que a diferença entre as soluções dos dois
problemas satisfaz um determinado problema de Cauchy. De seguida, obtêm-se estimativas
para as soluções do problema auxiliar. Finalmente, estas estimativas são usadas para concluir a
estimativa pretendida.
41
Método dos Elementos Finitos para problemas com Fronteiras Livres
Lema 3.11 Seja u solução fraca de (3.1) e v solução fraca de (3.3). A função w = u − v satisfaz
a equação
Z Z
w(−ζt − A∆ζ + B∇ζ) dxdt = (F [v] − Fε [v])∆ζ + (G[v] − Gε [v])∇ζ dxdt (3.7)
ΩT ΩT
com
ζ ∈ L2 (0, T ; H01 (Ω)) ∩ L∞ (0, T ; L∞ (Ω)), ζt ∈ L2 (0, T ; L2 (Ω)), ζ(x, T ) = 0,
γ γ
1 u|u| − v|v|
A= , B = A∇ln(γ + 1) − D, (3.8)
γ+1 u−v
u|u|γ ln(u2 ) − v|v|γ ln(v 2 )
∇γ
D= ,
2(γ + 1) u−v
u|u|γ
F (x, t) = F [u(x, t)] = ,
γ+1
Z v
γ
Fε (x, t) = Fε [v(x, t)] = (ς 2 + ε2 ) 2 dς,
0
ln(u2 )
∇γ γ 1
G[v] = u|u| −
γ+1 2 γ+1
e v
∇γ 2
Z
γ
Gε [v] = (ς + ε2 ) 2 ln(ς 2 + ε2 ) dς.
0 2
Demonstração
As funções F e Fε satisfazem
ln(u2 )
γ ∇γ γ 1
∇F = |u| ∇u + u|u| −
γ+1 2 γ+1
e v
∇γ 2
Z
γ γ
2 2
∇Fε = (v + ε ) ∇v + 2 (ς + ε2 ) 2 ln(ς 2 + ε2 ) dς.
0 2
Logo,
ln(u2 )
γ ∇γ γ 1
|u| ∇u = ∇F [u] − u|u| − = ∇F [u] − G[u]
γ+1 2 γ+1
e v
∇γ 2
Z
γ γ
(v 2 + ε2 ) 2 ∇v = ∇Fε − (ς + ε2 ) 2 ln(ς 2 + ε2 ) dς = ∇Fε [v] − Gε [v].
0 2
Z Z Z Z T
− vζt dxdt + aε (x, v)∇v · ∇ζ dxdt − f ζ dxdt = − vζ dx
ΩT ΩT ΩT Ω t=0
da forma
Z Z Z Z Z T
− uζt dxdt + ∇F [u]∇ζ dxdt − G[u]∇ζ dxdt − f ζ dxdt = − uζ dx ,
ΩT ΩT ΩT ΩT Ω t=0
Z Z Z Z Z T
− vζt dxdt + ∇Fε [v]∇ζ dxdt − Gε [v]∇ζ dxdt − f ζ dxdt = − vζ dx .
ΩT ΩT ΩT ΩT Ω t=0
42
Método dos Elementos Finitos para problemas com Fronteiras Livres
Z Z Z Z Z T
− vζt dxdt − Fε [v]∆ζ dxdt − Gε [v]∇ζ dxdt − f ζ dxdt = − vζ dx .
ΩT ΩT ΩT ΩT Ω t=0
Z T Z Z T Z Z T Z
− wζt dxdt − (F [u] − Fε [v])∆ζ dxdt = (G[u] − Gε [v])∇ζ dxdt. (3.9)
0 Ω 0 Ω 0 Ω
Considere-se
F [u] − Fε [v] = (F [u] − F [v]) + (F [v] − Fε [v]) e G[u] − Gε [v] = (G[u] − G[v]) + (G[v] − Gε [v]).
Pode escrever-se F [u] − F [v] = Aw e G[u] − G[v] = −Bw, com A e B de nidos em (3.8). Então,
a equação (3.9) pode ser reescrita da forma (3.7), como pretendido.
onde > 0 é um parâmetro pequeno, φ ∈ L2 (ΩT ) uma função arbitrária e A e B estão de nidas
em (3.8). Se
φ2
≤ C, (3.11)
A+
então Z Z tZ
2
|∇η| dx + (A + )(∆η)2 dxdt ≤ C.
Ω 0 Ω
Demonstração
Demonstra-se facilmente que
|B|2
0 ≤ A ≤ C, |D| ≤ C, |B| ≤ C, ≤ C,
A
43
Método dos Elementos Finitos para problemas com Fronteiras Livres
Tem-se que
|B|2
Z Z
1
|I1 | ≤ (A + )(∆η)2 dx + |∇η|2 dx
4 Ω Ω A+
e
φ2
Z Z
1 2
|I2 | ≤ (A + )(∆η) dx + dx.
4 Ω Ω A+
Juntando as estimativas calculadas, obtém-se
|B|2 φ2
Z Z Z Z
1 d 1
|∇η|2 dx + (A + )(∆η)2 dx ≤ (∇η)2 dx + dx.
2 dt Ω 2 Ω Ω A+ Ω A+
Já se provou que
|B|2 |B|2
≤ ≤ C|∇γ|2 ≤ C.
A+ A
Pela hipótese (3.11) e pelo teorema de Gronwall, obtém-se o pretendido.
então
γ+
+1 1
ku − vkL2 (ΩT ) ≤ Cε 2 ,
γ+
+1
2
Demonstração
Seja w = u − v, então satisfaz a equação (3.7). Escolhendo ζ(x, t) = η(x, T − t) solução de
(3.10) vem que
Z Z
wφ dxdt = (w ∆η + (F [v] − Fε [v])∆η + (G[v] − Gε [v])∇η) dxdt =
ΩT ΩT
= I3 + I4 + I5 .
Tem-se que
v v v v
v|v|γ
Z Z Z Z
γ γ γ
F [v] − Fε [v] = − (ς 2 + ε2 ) 2 dς = |ς|γ dς − (ς 2 + ε2 ) 2 dς = |ς|γ − (ς 2 + ε2 ) 2 dς,
γ+1 0 0 0 0
então
|F [v] − Fε [v]| ≤ C sup |v|ε
e Z v
ln(u2 )
∇γ
γ 1 ∇γ 2 γ
G[v] − Gε [v] = u|u| − − (ς + ε2 ) 2 ln(ς 2 + ε2 ) dς
Z v γ+1 2 Zγ v+ 1 0 2
γ ∇γ 2 2 γ2
= |ς| ∇γln(|ς|) dς − (ς + ε ) ln(ς 2 + ε2 ) dς
0
Z v 0 2
∇γ 2 γ
= |ς|γ ∇γln(|ς|) − (ς + ε2 ) 2 ln(ς 2 + ε2 ) dς,
0 2
logo
|G[v] − Gε [v]| ≤ C sup |v|ε.
Pela desigualdade de Schwarz e pelas estimativas calculadas no lema anterior, demonstra-se
que
Z 21
1 p 1
2
|I3 | ≤ sup(|w|) 2 T meas(Ω) (∆η) dxdt ≤ C 2 ,
ΩT
Z 21
− 12 1
p
|I4 | ≤ sup(|F [v] − Fε [v]|) T meas(Ω) 2
(∆η) dxdt ≤ Cε− 2
ΩT
44
Método dos Elementos Finitos para problemas com Fronteiras Livres
e Z 21
p
|I5 | ≤ sup(|G[v] − Gε [v]|) T meas(Ω) (∇η)2 dxdt ≤ Cε.
ΩT
Logo, Z
1 1
wφ dxdt ≤ C( 2 + ε− 2 + ε).
ΩT
Se se considerar = ε, então
Z
1
wφ dxdt ≤ Cε 2 .
ΩT
Ao considerar φ = |w|κ sign(w), com κ uma constante positiva, pode deduzir-se que
1 1
φ2 φ2 |u − v|2κ
≤ = u|u|γ −v|v|γ ≤ C,
A+ A
(γ+1)(u−v)
Nesta secção, segue-se [72] e aplica-se o método de Galerkin, contínuo no espaço e descon-
tínuo no tempo, para obter soluções aproximadas do problema regularizado. Sempre que não
houver perigo de confusão, deixar-se-á de explicitar a dependência de a em relação a x.
Seja h > 0, Th uma partição de Ω e Shr o espaço de elementos nitos
s
X
δs
Shr = {W = Whδ : [0, +∞[→ Sh | W|In = tn χn (x), χn ∈ Shr }.
n=0
ou seja,
ni−1
XZ Z ni−1
XZ Z ni−1
XZ Z Z T
−V ζt dxdt + aε (x, V )∇V · ∇ζ dxdt = f ζ dxdt − V ζdx .
n=0 In Ω n=0 In Ω n=0 In Ω Ω 0
45
Método dos Elementos Finitos para problemas com Fronteiras Livres
ni−1
XZ Z ni−1
X Z tn+1 ni−1
XZ Z
Vt ζ dxdt − V ζ dx + aε (x, V )∇V · ∇ζ dxdt =
In Ω Ω + In Ω
n=0 n=0 tn n=0
ni−1
XZ Z Z
= f ζ dxdt + V0+ ζ0 dx,
n=0 In Ω Ω
ni−1
XZ Z ni−1
XZ ni−1
XZ Z ni−1
XZ Z
Vt ζ dxdt + [V ]n ζ dx + aε (x, V )∇V · ∇ζ dxdt = f ζ dxdt.
n=0 In Ω n=1 Ω n=0 In Ω n=0 In Ω
Escolhendo Shr
δs
para espaço das funções teste, obtém-se
ni−1
XZ Z ni−1
XZ ni−1
XZ Z
Vt W dxdt + [V ]n Wn dx + aε (x, V )∇V · ∇W dxdt =
n=0 In Ω n=1 Ω n=0 In Ω
ni−1
XZ Z
= f W dxdt. (3.12)
n=0 In Ω
Como W pode não ser contínuo nos ti , i = 0, . . . , ni pode escolher-se W de modo que se
anule fora de In , n = 0, . . . , ni − 1. Neste caso, a equação (3.12) transforma-se em ni equações,
uma para cada In . Assim, o problema discreto pode ser encontrar V ∈ Shr
δs
de modo que
Z Z Z Z Z
Vt W dxdt + +
Vn−1 Wn−1
+
dx + aε (x, V )∇V · ∇W dxdt =
In Ω Ω In Ω
Z Z Z
= f W dxdt + Vn−1 Wn−1
+
dx, (3.13)
In Ω Ω
δs
∀W ∈ Shr , ∀n ∈ {0, . . . , ni}.
Tem agora de se provar que este problema tem solução única. Note-se que para demonstrar
a existência e unicidade de solução para o problema (3.12), é su ciente provar que, em cada
intervalo In e com o valor de Vn−1 xo, existe uma única solução para a equação (3.13).
Teorema 3.15 Se Vn−1 ∈ L2 (Ω) e f ∈ L2 (Ω × In ), então o problema (3.13) tem uma solução V .
Demonstração
Z Seja n ≥ 1 xo. ParaZ cada
Z h, δ > 0, de Zne-se a aplicação contínua F : Shr → Shr da forma
δs δs
Z
F (V ) W dxdt = Vt W dxdt + +
Vn−1 Wn−1
+
dx+
In Ω In Ω Ω
Z Z Z Z Z
+ aε (x, V )∇V · ∇W dxdt − f W dxdt − Vn−1 Wn−1
+
dx, δs
∀W ∈ Shr .
In Ω In Ω Ω
escolhendo W = V ,
46
Método dos Elementos Finitos para problemas com Fronteiras Livres
Z Z Z Z Z
+
F (V ) V dxdt = Vt V dxdt + (Vn−1 )2 dx+
In Ω In Ω Ω
Z Z Z Z Z
2 +
+ aε (x, V )|∇V | dxdt − f V dxdt − Vn−1 Vn−1 dx.
In Ω In Ω Ω
Se V pertencer a B = {W ∈ Shr
δs
: kW kL2 (Ω×In ) ≤ , > Cε1γ + kVn−1 kL2 (Ω) + (Cεγ1+ )2 kfn kL2 (Ω×In ) },
então Ω F (V ) V dx ≥ 0, ∀V ∈ ∂B. O corolário do teorema de ponto xo de Brouwer implica a
R
Teorema 3.16 Se Vn−1 ∈ L2 (Ω) e f ∈ L2 (Ω × In ), então o problema (3.13) tem uma única
solução V .
Demonstração
Sejam Vn−1 e f duas funções conhecidas e V1 e V2 duas soluções de (3.13), então, para
W ∈ Shr
δs
, tem-se que
Z Z Z Z Z
(V1 )t W dxdt + n−1 Wn−1 dx +
(V1 )+ +
aε (x, V1 )∇V1 · ∇W dxdt =
In Ω Ω In Ω
Z Z Z
= f W dxdt + Vn−1 Wn−1
+
dx
In Ω Ω
e Z Z Z Z Z
(V2 )t W dxdt + n−1 Wn−1 dx +
(V2 )+ +
aε (x, V2 )∇V2 · ∇W dxdt =
In Ω Ω In Ω
Z Z Z
= f W dxdt + Vn−1 Wn−1
+
dx.
In Ω Ω
47
Método dos Elementos Finitos para problemas com Fronteiras Livres
Escrevendo
Z 1
d
aε (x, V1 )∇V1 − aε (x, V2 )∇V2 = (aε (x, ςV1 + (1 − ς)V2 )∇(ςV1 + (1 − ς)V2 )) dς =
0 dς
Z 1 Z 1
γ
= aε (x, ςV1 +(1−ς)V2 )∇V dς+ γ((ςV1 +(1−ς)V2 )2 +ε2 ) 2 −1 (ςV1 +(1−ς)V2 )(ς∇V1 +(1−ς)∇V2 )V dς
0 0
= A∇V + BV,
com Z 1
A= aε (x, ςV1 + (1 − ξ)V2 ) dς
0
e Z 1
γ
B= γ((ςV1 + (1 − ς)V2 )2 + ε2 ) 2 −1 (ςV1 + (1 − ς)V2 )(ς∇V1 + (1 − ς)∇V2 ) dς,
0
+
Neste problema φ é uma função arbitrária de L2 (Ω × int(In )). Uma vez que εγ ≤ A ≤ C e
+
|B| ≤ Cεγ , o problema tem uma solução fraca em Shr
δs
. Então a equação (3.15) transforma-se
em Z Z
V φ dxdt = 0, (3.17)
In Ω
pois Wn = W (x, tn ) = η(x, tn−1 ) = 0. Como (3.17) é válida para qualquer φ ∈ L2 (Ω × int(In )),
conclui-se que V = 0 e, consequentemente, a solução é única.
O resto desta secção é dedicada à demonstração de uma estimativa do erro em termos
dos parâmetros da discretização para a solução totalmente discreta do problema (3.12). Esta
demonstração é muito semelhante à demonstração do teorema 3.13 e está também repartida
em dois lemas e um teorema.
Lema 3.17 Seja v solução fraca do problema aproximado e V a solução discreta, então a
função e = V − v satisfaz, para todo W ∈ Shr
δs
, a equação
Z Z ni−1
XZ Z
(−eWt + A∇e · ∇W + eB · ∇W ) dxdt = en [W ]n dx − eni Wni dx, ∀W ∈ Shr
δs
(3.18)
I Ω n=1 Ω Ω
48
Método dos Elementos Finitos para problemas com Fronteiras Livres
com Z 1
A= aε (x, ςV + (1 − ξ)v) dς
0
e Z 1
γ
B= γ((ςV + (1 − ς)v)2 + ε2 ) 2 −1 (ςV + (1 − ς)v)(ς∇V + (1 − ς)∇v) dς.
0
Demonstração
A equação (3.12) pode ser reescrita da forma
ni−1
XZ Z ni−1
XZ ni−1
XZ Z
− V Wt dxdt − Vn [W ]n dx + aε (x, V )∇V · ∇W dxdt =
n=0 In Ω n=1 Ω n=0 In Ω
ni−1
XZ Z Z
= f W dxdt − Vni Wni dx.
n=0 In Ω Ω
Somando
Z Z ni−1
XZ Z Z Z Z Z
− V Wt dxdt− Vn [W ]n dx+ aε (x, V )∇V ·∇W dxdt = f W dxdt− Vni Wni dx.
I Ω n=1 Ω I Ω I Ω Ω
Z Z ni−1
XZ Z Z Z Z Z
− vWt dxdt− vn [W ]n dx+ aε (x, v)∇v·∇W dxdt = f W dxdt− vni Wni dx.
I Ω n=1 Ω I Ω I Ω Ω
Z Z ni−1
XZ Z Z Z
− eWt dxdt− en [W ]n dx+ (aε (x, V )∇V −aε (x, v)∇v)·∇W dxdt+ eni Wni dx = 0,
I Ω n=1 Ω I Ω Ω
= A∇e + Be.
Então,
Z Z ni−1
XZ Z Z Z Z Z
− eWt dxdt− en [W ]n dx+ A∇e·∇W dxdt+ eB·∇W dxdt+ eni Wni dx = 0.
I Ω n=1 Ω I Ω I Ω Ω
49
Método dos Elementos Finitos para problemas com Fronteiras Livres
com Z 1
A= aε (x, ςV + (1 − ς)v) dς,
0
e Z 1
γ
B= γ((ςV + (1 − ς)v)2 + ε2 ) 2 −1 (ςV + (1 − ς)v)(ς∇V + (1 − ς)∇v) dς.
0
Se
φ2
≤ C, (3.20)
A
então Z Z tZ
|∇η|2 dx + A(∆η)2 dxdt ≤ C. (3.21)
Ω 0 Ω
Demonstração
+ +
Tem-se que εγ ≤ A ≤ C e |B| ≤ Cεγ , logo pela teoria clássica [52], para qualquer ε > 0
e φ ∈ L2 (ΩT ) este problema tem uma única solução fraca contínua η tal que ηt , ∆η ∈ L2 (ΩT ).
A seguir encontram-se estimativas para η e respetivas derivadas.
Multiplicando (3.19) por ∆η e integrando, obtém-se
Z Z Z Z
ηt ∆η dx − A(∆η)2 dx = − B̃∇η∆η dx + φ∆η dx,
Ω Ω Ω Ω
com
1
∇γ
Z
B̃ = aε (x, ςV + (1 − ς)v) ln((ςV + (1 − ς)v)2 + ε2 ) dς.
0 2
Tem-se que
|B̃|2
Z Z
1
|I1 | ≤ A(∆η)2 dx + |∇η|2 dx
4 Ω Ω A
e
φ2
Z Z
1 2
|I2 | ≤ A(∆η) dx + dx.
4 Ω Ω A
|B̃|2 φ2
Z Z Z Z
1 d 1
|∇η|2 dx + A(∆η)2 dx ≤ |∇η|2 dx + dx.
2 dt Ω 2 Ω Ω A Ω A
R 2
1
C 0
aε (x, ςV + (1 − ς)v) dς
≤ R1 ≤ C.
0
aε (x, ςV + (1 − ς)v) dς
50
Método dos Elementos Finitos para problemas com Fronteiras Livres
Teorema 3.19 Seja v solução fraca do problema aproximado e V a solução discreta, então
γ+
+1 +
kV − vkL2 (ΩT ) ≤ Cε−γ (hr+1 kvkL∞ (I,H r+1 (Ω)) + δ s+1 |v (s+1) |W s+1,∞ (I,L2 (Ω)) ),
γ+
+1
2
Demonstração
Primeiro, de ne-se Ṽ ∈ Shr
δs
em cada In por
Z Z
Ṽ χ dx = ṽχ dx, ∀t ∈ In , χ ∈ Shr ,
Ω Ω
onde
ṽ(tn ) = v(tn ), n = 1, . . . , ni
e Z Z Z Z
δ,s−1
ṽW dxdt = vW dxdt, ∀W ∈ Shr .
In Ω In Ω
Para veri car que Ṽ está bem de nido, pode consultar-se [72].
Pode então escrever-se e = θ + ρ, onde θ = V − Ṽ e ρ = Ṽ − v.
Pelo teorema 2.86, obtêm-se as estimativas seguintes:
kρkL2 (Ω) ≤ C(hr+1 kvkL∞ (In ;H r+1 (Ω)) + δ s+1 |v|W s+1,∞ (In ;L2 (Ω)) ), ∀t ∈]0, T ].
Z Z ni−1
XZ Z Z Z Z Z
− θWt dxdt− θn [W ]n dx− A∇θ·∇W dxdt+ θB·∇W dxdt+ θni Wni dx =
I Ω n=1 Ω I Ω I Ω Ω
Z Z ni−1
XZ Z Z Z Z Z
= ρWt dxdt+ ρn [W ]n dx+ A∇ρ·∇W dxdt+ ρB ·∇W dxdt+ ρni Wni dx.
I Ω n=1 Ω I Ω I Ω Ω
Z Z ni−1
XZ Z
ρWt dxdt = 0, ρn [W ]n dx = 0 e ρni Wni dx = 0.
I Ω n=1 Ω Ω
51
Método dos Elementos Finitos para problemas com Fronteiras Livres
Z Z 12
−γ +
≤ Cε (hr+1
kvkL∞ (In ;H r+1 (Ω)) + δ s+1
|v|W s+1,∞ (In ;L2 (Ω)) ) A(∆h W ) dxdt2
I Ω
+
≤ Cε−γ (hr+1 kvkL∞ (In ;H r+1 (Ω)) + δ s+1 |v|W s+1,∞ (In ;L2 (Ω)) ),
Z Z 12 Z Z 21
2 2
|I2 | ≤ ρ dxdt (B · ∇W ) dxdt
I Ω I Ω
Z Z 21
r+1 s+1 −γ + 2
≤ C(h kvkL∞ (In ;H r+1 (Ω)) +δ |v|
W s+1,∞ (I n ;L2 (Ω))
)Cε (∇W ) dxdt
I Ω
+
≤ Cε−γ (hr+1 kvkL∞ (In ;H r+1 (Ω)) + δ s+1 |v|W s+1,∞ (In ;L2 (Ω)) ).
+
Considerando φ = |θ|γ /2
sign(θ), obtém-se o resultado pretendido.
No caso de s = 0, as modi cações são evidentes.
Juntando este teorema com o teorema 3.13 e as estimativas para as derivadas de v, pode-se
estimar o erro entre u e V da forma
γ+
+1 + + 1
kV − ukL2 (ΩT ) ≤ Cε−αγ hr+1 + Cε−βγ δ s+1 + Cε 2 .
γ+
+1
2
γ+
+1 1 + +
ku − V kL2 (ΩT ) ≤ C(ε 2 + h2 ε−9γ + δε−3γ ).
γ+
+1
2
1 9 + 1 +
Escolhendo h = ε 4 + 2 γ e δ = ε 2 +3γ , obtém-se
γ+
+1 1
ku − V kL2 (ΩT ) ≤ Cε 2 .
γ+
+1
2
52
Método dos Elementos Finitos para problemas com Fronteiras Livres
1) Para k ≥ 1, Wk ∈ Shr
δ0
é calculado de modo que para todo W ∈ Shr
δ0
Z Z Z Z
Wk W dx + δ aε (x, Wk−1 )∇Wk · ∇W dx = gn W dx + W0 W dx, ∀W ∈ Shr .
Ω Ω Ω Ω
(3.23)
2) Se
kWk − Wk−1 kL2 (Ω) ≥ tol,
3) Se
kWk − Wk−1 kL2 (Ω) < tol,
de ne-se Vn = Wk e termina.
Prova-se, antes de mais, que todas as estimativas obtidas pelo algoritmo são limitadas.
Lema 3.21 Supondo que kVn−1 kL2 (Ω) ≤ C, onde C não depende de n e que gn ∈ L2 (Ω), então
existe 0 < C 0 < ∞ tal que
kWk kL2 (Ω) ≤ C 0 , ∀k > 0,
Demonstração
Seja W = Wk em (3.23), então
Z Z Z Z
Wk2 dx + δ aε (x, Wk−1 )(∇Wk ) dx =
2
gn Wk dx + W0 Wk dx,
Ω Ω Ω Ω
ou seja,
Z Z Z
γ
kWk k2L2 (Ω) + δ 2
(Wk−1 + ε2 ) 2 ∇Wk · ∇Wk dx = gn Wk dx + W0 Wk dx.
Ω Ω Ω
γ
Como δ + ε2 ) 2 ∇Wk · ∇Wk dx é não negativo,
R 2
Ω
(Wk−1
Z Z
kWk k2L2 (Ω) ≤ gn Wk dx + W0 Wk dx.
Ω Ω
kWk k2L2 (Ω) ≤ kgn kL2 (Ω) kWk kL2 (Ω) + kW0 kL2 (Ω) kWk kL2 (Ω) .
Simpli cando
kWk kL2 (Ω) ≤ kgn kL2 (Ω) + kW0 kL2 (Ω) .
Teorema 3.22 Seja tol > 0. Existe C > 0 tal que: δ < 1
C(h−2 ε−γ + +εγ + )
implica a existência de
∗
k ∈ N de modo que
kWk − Wk−1 kL2 (Ω) ≤ tol, ∀k > k ∗ .
Demonstração
53
Método dos Elementos Finitos para problemas com Fronteiras Livres
kWk − Wk−1 kL2 (Ω) < C1 kWk−1 − Wk−2 kL2 (Ω) , C1 < 1, ∀k ≥ 2.
e
Z Z Z Z
Wk−1 W dx + δ aε (x, Wk−2 )∇Wk−1 · ∇W dx = gn W dx + W0 W dx. (3.25)
Ω Ω Ω Ω
Considere-se
Ek = Wk − Wk−1 ∈ Sh e W = Ek ,
então Z Z
Ek2 dx + δ (aε (x, Wk−1 )∇Wk − aε (x, Wk−2 )∇Wk−1 ) · ∇Ek dx = 0.
Ω Ω
Escreve-se
aε (x, Wk−1 )∇Wk − aε (x, Wk−2 )∇Wk−1 =
= aε (x, Wk−1 )∇Wk−1 − aε (x, Wk−1 )∇Wk−1 + aε (x, Wk−1 )∇Wk − aε (x, Wk−2 )∇Wk−2 +
+aε (x, Wk−2 )∇Wk−2 − aε (x, Wk−2 )∇Wk−1 =
= aε (x, Wk−1 )∇Ek + aε (x, Wk−1 )∇Wk−1 − aε (x, Wk−2 )∇Ek−1 − aε (x, Wk−2 )∇Wk−2 .
Logo,
Z Z Z
Ek2 dx + δ aε (x, Wk−1 )|∇Ek | dx + δ 2
(aε (x, Wk−1 )∇Wk−1 − aε (x, Wk−2 )∇Wk−2 ) · ∇Ek dx =
Ω Ω Ω
Z
=δ aε (x, Wk−2 )∇Ek−1 · ∇Ek dx.
Ω
Volta a escrever-se
com Z 1
A= aε (x, ςWk−1 + (1 − ς)Wk−2 ) dς
0
e
Z 1
γ
B= γ((ςWk−1 + (1 − ς)Wk−2 )2 + ε2 ) 2 −1 (ςWk−1 + (1 − ς)Wk−2 )(ς∇Wk−1 + (1 − ς)∇Wk−2 ) dς.
0
+ +
Tem-se que εγ ≤ A ≤ C e |B| ≤ Cεγ . Então,
Z Z Z Z
Ek2 dx + δ aε (x, Wk−1 )|∇Ek |2 dx + δ A∇Ek−1 · ∇Ek dx + δ BEk−1 · ∇Ek dx =
Ω Ω Ω Ω
54
Método dos Elementos Finitos para problemas com Fronteiras Livres
Z
=δ aε (x, Wk−2 )∇Ek−1 · ∇Ek dx,
Ω
ou seja,
Z Z Z Z
Ek2 dx+δ aε (x, Wk−1 )|∇Ek |2 dx = δ (aε (x, Wk−2 )−A)∇Ek−1 ·∇Ek dx−δ BEk−1 ·∇Ek dx.
Ω Ω Ω Ω
+ + + +
ε−γ εγ εγ εγ
≤ δ(C k∇Ek−1 k2L2 (Ω) + k∇Ek k2L2 (Ω) ) + δ(C k∇Ek−1 k2L2 (Ω) + k∇Ek k2L2 (Ω) ).
2 2 2 2
Simpli cando,
+ +
Cδε−γ Cδεγ
Z
Ek2 dx ≤ k∇Ek−1 k2L2 (Ω) + kEk−1 k2L2 (Ω) .
Ω 2 2
Atendendo à desigualdade inversa k∇Ek−1 kL2 (Ω) ≤ ChkEk−1 kL2 (Ω) , obtém-se
+ +
kEk k2L2 (Ω) ≤ Cδ(h−2 ε−γ + εγ )kEk−1 k2L2 (Ω)
+ + 1
kEk kL2 (Ω) ≤ (Cδ(h−2 ε−γ + εγ )) 2 kEk−1 kL2 (Ω) .
Se
1
δ< ,
C(h−2 ε−γ + + εγ + )
+ + 1
então obtém-se o pretendido, considerando C1 = (Cδ(h−2 ε−γ + εγ )) 2 < 1.
Para testar o método, considera-se primeiro um problema com o segundo membro calculado
de modo que o problema tenha solução exata.
Seja ΩT =] − 1, 1[2 ×]0, T ] e considere-se o problema (3.1), onde
x 2 y 2
γ =2− −
2 2
seja solução do problema. Utilizou-se uma malha estruturada como mostra a gura 3.1.
A tolerância usada foi tol = 10−10 . Para aproximar a solução, foram utilizados polinómios de
Lagrange de grau 1.
A gura 3.2, mostra a evolução da solução obtida com ε = 10−8 , δ = 10−3 e h = 0.1, em
alguns valores de t. Observa-se que o comportamento é idêntico ao da solução exata.
55
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 3.1: Malha espacial para h = 0.1 usada no problema com o segundo membro arti cial.
Figura 3.2: Evolução no tempo da aproximação obtida no problema com o segundo membro arti cial.
Com vista a analisar a convergência, zeram-se várias simulações para diversas combinações
de ε, h e δ. Na primeira imagem da gura 3.3, apresenta-se o logaritmo do erro da aproximação,
na norma de L2 (Ω), em função do logaritmo de h para alguns valores de δ. Conclui-se que o
erro diminui ao reduzir-se h e δ e que a ordem de convergência de h é aproximadamente 2.
A segunda imagem da gura 3.3 tem representado o logaritmo do erro da aproximação, na
norma de L2 (Ω), em função do logaritmo de δ para alguns valores de h. Pode observar-se que
o erro diminui ao reduzir-se h e δ e que a ordem de convergência de δ é aproximadamente 1.
A gura 3.4 tem representado o logaritmo do erro da aproximação, na norma de L2 (Ω), em
função do logaritmo de ε para alguns valores de h e com δ = 10−3 . O grá co evidencia que o
erro diminui ao reduzir-se ε e que a ordem de convergência é aproximadamente 1, melhor do
que é demonstrado no teorema 3.13.
56
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 3.3: Erro da aproximação obtida no problema com o segundo membro arti cial em função de h e
de δ.
Figura 3.4: Erro da aproximação obtida no problema com o segundo membro arti cial em função de ε.
Considere-se agora um problema onde a solução inicial tem um buraco no suporte. Seja
ΩT =] − 1.5, 1.5[2 ×]0, T ] e um problema do tipo (3.1) com
x 2 y 2
γ= + + 1.1.
2 2
57
Método dos Elementos Finitos para problemas com Fronteiras Livres
desaparece ao m de um instante de tempo nito t∗ > 0. Após esse instante, o suporte torna-se
compacto e a solução comporta-se de um modo semelhante às soluções de Barenblat.
Figura 3.5: Evolução no tempo da aproximação obtida no problema com suporte não convexo.
Para terminar, simula-se um problema com um termo de absorção com expoente variável da
forma
∂u
∂t = div(uγ(x,t) ∇u) − uσ(x,t) em ΩT
u=0 em ΓT .
u(x, 0) = u0 (x) em Ω
Não foi provada a convergência para este tipo de problemas, mas as modi cações não são
muito signi cativas. Considere-se então os expoentes como sendo
x2 + y 2
γ= e σ = x2 + y 2 + 1 + e−t .
t2 + 1
O domínio escolhido foi Ω =] − 1, 1[2 e usou-se uma malha estruturada composta por 3600
triângulos, como mostra a gura 3.6. A solução foi aproximada por polinómios de grau 1 no
espaço e grau zero no tempo. Terminou-se a simulação em T = 0.1. Na gura 3.7, está
58
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 3.6: Malha espacial para h = 1/15 usada no problema com absorção.
O objetivo desta secção é implementar, em Matlab, um método que permita encontrar so-
luções aproximadas para o problema (3.1), em domínios espaciais livres de duas dimensões. A
técnica utiliza o método dos elementos nitos com uma malha móvel e aproximações de grau
maior que um. O uso de funções base seccionalmente polinomiais de grau r ≥ 1 vai permitir
a redução do número de elementos nitos mantendo a precisão. É aplicado um sistema de
EDP parabólicas que envolve a velocidade dos nós para mover a malha para certas regiões do
domínio. As equações da malha e do problema físico são resolvidas simultaneamente.
A abordagem do problema num domínio livre exige uma nova formulação que contemple o
movimento da fronteira. Nesta secção formula-se, de um modo exato, o problema num domínio
livre de R2 e de ne-se o que se entende por solução fraca neste contexto.
Seja Ω(t) ⊂ R2 e ]0, T ] um intervalo de tempo nito. Considere-se o seguinte problema:
encontrar a função u e o domínio Ω(t) que satisfaçam as condições
∂u
∂t − div(|u|γ(x,t) ∇u) = f (x, t, u) em Ω(t)×]0, T ]
u(x, t) = 0 em ∂Ω(t) × [0, T ] , (3.26)
u(x, 0) = u0 (x) em Ω(t)
59
Método dos Elementos Finitos para problemas com Fronteiras Livres
Para o problema car bem de nido, tem de se obter uma equação para a velocidade da fron-
teira, ν = ( ∂x ∂y
∂t , ∂t ).
Pela lei da conservação da massa,
∂u
+ div(ν · u) = 0, (3.28)
∂t
A equação (3.29) indica a velocidade em cada x ∈ Ω(t). Uma vez que a solução tem a
propriedade de propagação com velocidade nita, é plausível assumir que limx→Γ ν existe e é
nito. Mesmo assim, a equação (3.29) não se adequa muito ao cálculo numérico da velocidade
da fronteira, pois na fronteira u = 0, logo, para que ela se mova, ∇u deve ser in nito, o que não
é viável em termos numéricos. Tem então de se procurar uma nova equação para a velocidade
normal da fronteira que se possa utilizar com o método dos elementos nitos. Como
então
∇|u|γ ∇γ γ
ν=− + |u| ln(|u|).
γ γ
Na fronteira, u = 0. Se γ > 1 e ∇γ é limitado em Ω, então u → 0 implica que
∇γ γ
|u| ln(|u|) → 0.
γ
Logo, a velocidade na fronteira pode ser dada por ν = − γ1 ∇|u|γ e a velocidade normal por
1 ∂|u|γ
νn = − . (3.30)
γ ∂n
De ne-se agora o que se entende por solução fraca para este problema. Seja ζ uma função
teste. Se se multiplicar a equação (3.31) por ζ e se for integrada em Ω(t), obtém-se
Z Z
∂u
ζ dx = div(uγ ∇u) ζ dx.
Ω(t) ∂t Ω(t)
60
Método dos Elementos Finitos para problemas com Fronteiras Livres
Então, obtém-se Z Z
∂u
ζ dx = − uγ 5 u . 5 ζ dx. (3.32)
Ω(t) ∂t Ω(t)
De nição 3.23 Uma função u(x, t) é dita solução fraca do problema (3.31) se:
γ
(i) u ∈ L∞ (0, T ; L∞ (Ω(t))), u 2 ∇u ∈ L2 (0, T ; L2 (Ω(t))), ∂u
∂t ∈ L2 (0, T ; W −1,2 (Ω(t))),
(ii) u = 0 em ∂Ω(t),
(iii) para qualquer função de teste ζ ∈ H u e para todo 0 < t ≤ T , a equação (3.32) é válida,
A condição (v) não é comum na de nição de solução fraca, pois surge implicitamente na equação
(3.32), mas não existe nenhuma demonstração exata de que (3.29) é válido na fronteira. Por
essa razão, decidiu-se impor (v) como uma condição na de nição.
61
Método dos Elementos Finitos para problemas com Fronteiras Livres
onde m é uma matriz quadrada ou um escalar que será chamada função monitor que, geral-
mente, depende de u e que será de nida mais à frente.
As direções de descida máxima de F são dadas pelos simétricos das derivadas parciais,
então pode de nir-se
∂x ∂F ∂y ∂F
=− , =− .
∂t ∂x ∂t ∂y
Na prática, por vezes, é preferível utilizar direções descendentes diferentes da mais rápida, mas
que permitam ajustar a descida. Então, utilizar-se-á uma modi cação das equações anteriores
∂x 1 ∂E ∂y 1 ∂E
=− , =− ,
∂t τ ∂x ∂t τ ∂y
com τ > 0 (um parâmetro a de nir pelo utilizador), que permite ajustar a velocidade com que
a malha se ajusta. Desta forma, as aplicações x(ξ, η, t) e y(ξ, η, t) de Ωc para Ωf são calculadas
resolvendo as equações
∂x 1 ∂y 1
= div(m∇x), = div(m∇y).
∂t τ ∂t τ
Para m adequado é possível adicionar condições iniciais e na fronteira de modo a que exista
solução única.
onde f x (ξ, t) está relacionada com o movimento da fronteira. Se a fronteira é xa, então
f (ξ, t) = x0 (ξ). Se a fronteira é móvel, então f (ξ, t) é de nida pela linha poligonal obtida
x x
pela solução do sistema formado pelas equações (3.40) e (3.41). τ é um parâmetro positivo e
m é a função monitor. Para o problema (3.33) ser do tipo parabólico, m deve ser uma função
escalar positiva ou matriz de funções de nida positiva. Foram realizadas diferentes simulações
com diferentes funções monitor e escolheu-se
s 2 2
∂u ∂u
m= 1 + 100 + 100 ,
∂x ∂y
pois é conhecida por concentrar pontos onde o gradiente de u é elevado e por saber-se que o
gradiente da solução é máximo próximo da fronteira, onde se quer mais de nição. Para que a
transformação x seja aplicável à geração da malha, têm de garantir-se algumas propriedades.
A existência, unicidade e invertibilidade são as indispensáveis, além disso, deve garantir-se
adaptabilidade e suavidade.
62
Método dos Elementos Finitos para problemas com Fronteiras Livres
(3.33), obtém-se Z Z
∂x 1
ζ dξ = div(m∇x) ζ dξ.
Ωc ∂t Ωc τ
Aplicando o teorema de Green e supondo que ζ se anula na fronteira e tem regularidade
su ciente, obtém-se
Z Z
∂x 1
ζ dξ = − m∇x · ∇ζ dξ.
Ωc ∂t τ Ωc
De nição 3.24 Uma função x = (x(ξ, η, t), y(ξ, η, t)) é dita solução fraca do problema (3.33) se
√ √
(i) x, y ∈ L∞ (0, T ; W 1,2 (Ωc )), m∇x, m∇y ∈ L2 (0, T ; L2 (Ωc )), ∂x ∂y
∂t , ∂t ∈ L2 (0, T ; W −1,2 (Ωv )),
(ii) x = f x em ∂Ωc ,
Z Z
∂y 1
ζ dξ = − m∇y · ∇ζ dξ, ∀ζ ∈ H0x , (3.35)
Ωc ∂t τ Ωc
u
Shr = {uh ∈ C00 (Ω̄)|uh|Tk é um polinómio de grau r ∀Tk ∈ Th }
e
Shx = {xh ∈ C00 (Ω̄)|xh|Tk é um polinómio de grau 1 ∀Tk ∈ Th }.
np
X np
X
u≈U = Uj (t)ϕj , ζ= ζi (t)ϕj ,
j=1 i=1
onde (Xk , Yk ) = Pk são as coordenadas do vértice k. Alguns valores de Uj são utilizados para
satisfazer as condições da fronteira e os restantes valores são as incógnitas a determinar. Os
valores de χi são conhecidos ( arbitrários ).
Calculando a derivada em relação a t e aplicando a regra da cadeia obtemos, obtém-se
63
Método dos Elementos Finitos para problemas com Fronteiras Livres
np nv nv
∂U X ∂Uj X ∂ϕj ∂Xk X ∂ϕj ∂Yk
= ( ϕj + Uj + Uj )
∂t j=1
∂t ∂Xk ∂t ∂Yk ∂t
k=1 k=1
np nv nv
X ∂Uj ∂ϕj X ∂Xk ∂ϕj X ∂Yk
= ( ϕj − Uj ψk − Uj ψk ),
j=1
∂t ∂x ∂t ∂y ∂t
k=1 k=1
np np
Z nv nv
!
X ∂Uj ∂ϕj X ∂Xk ∂ϕj X ∂Yk X
ϕj − Uj ψk − Uj ψk ζi ϕi dx =
Ω(t) j=1 ∂t ∂x ∂t ∂y ∂t i=1
k=1 k=1
Z np
XX np
=− U γ Uj ζi ∇ϕj · ∇ϕi dx.
Ω(t) j=1 i=1
np Z nv Z nv Z
X ∂Uj X ∂U ∂Xk X ∂U ∂Yk
ϕj ϕi dx − ϕi ψk dx − ϕi ψk dx =
j=1 Ω(t) ∂t Ω(t) ∂x ∂t Ω(t) ∂y ∂t
k=1 k=1
np Z
X
=− U γ ∇ϕj · ∇ϕi dx Uj , i = 1, . . . , np. (3.36)
j=1 Ω(t)
nv
X nv
X nv
X
x≈X= Xj (t)ψj (ξ, η), y≈Y = Yj (t)ψj (ξ, η), ζ= ζi (t)ψi (ξ, η).
j=1 j=1 i=1
e
nv Z nv Z
X ∂Yj 1X
ψj ψi dξ = − Yj m∇ψj · ∇ψi dξ , i = 1, . . . , nv. (3.38)
j=1
∂t Ωc τ j=1 Ωc
Juntando as equações (3.36), (3.37) e (3.38), tem-se o sistema ODE nal com as incógnitas
∂Uj ∂Xk
∂t , ∂t e ∂t . Este sistema pode ser escrito da forma matricial com
∂Yk
∂U
M 11 M 12 M 13 ∂t F1
M 22 ∂X (3.39)
2
0 0 = F ,
∂t
∂Y
0 0 M 33 ∂t F3
onde
64
Método dos Elementos Finitos para problemas com Fronteiras Livres
Z nt Z
X
11
Mij = ϕi ϕj dx = ϕi ϕj dx,
Ωf k=1 Tk∗
Z np nt Z np
12
X ∂ϕl X X ∂ϕl
Mij =− ϕi ψj Ul dx = − ϕi ψj Ul dx,
Ωf ∂x ∗ ∂x
l=1 k=1 Tk l=1
Z np nt Z np
13 ∂ϕl X X X ∂ϕl
Mij =− Ul dx = −
ϕi ψj ϕi ψj Ul dx,
Ωf ∂y Tk∗ ∂y
l=1 k=1 l=1
Z nt Z
X
22
Mij = ψi ψj dξ = ψi ψj dξ,
Ωc k=1 Tk
Z Xnt Z
33
Mij = ψi ψj dξ = ψi ψj dξ,
Ωc k=1 Tk
np
X Z Xnp Z
Fi1 = − Uj ( Ul ϕl )γ ∇ϕj · ∇ϕi dx + f ϕi dx =,
j=1 Ωf l=1 Ωf
nt
X np Z
X np
X Z
= − ( Ul ϕl )γ ∇ϕj · ∇ϕi dxUj + f ϕi dx,
k=1 j=1 Tk∗ l=1 Tk∗
nv Z nt nv Z
1 X 1 X X
Fi2 = − m∇ψj · ∇ψi dξXj = − m∇ψj · ∇ψi dξXj ,
τ j=1 Ωc
τ j=1 Tk
k=1
nv Z nt nv Z
1 X 1 X X
Fi3 = − m∇ψj · ∇ψi dξYj = − m∇ψj · ∇ψi dξYj .
τ j=1 Ωc
τ j=1 T kk=1
De forma a satisfazer as condições especí cas do problema, alguns destes elementos devem
ser alterados. O sistema (3.39) vai ser resolvido utilizando o integrador ode15s [66], no intervalo
]0, T ].
Para aplicar o método dos elementos nitos é necessário discretizar a equação (3.30). A
fronteira é aproximada por uma linha poligonal. É necessário então calcular-se as derivadas
parciais de uγ e o vetor normal à curva φ(x, y, t) = 0 em cada ponto da malha que pertença
γ y 1 ∂uγ
à fronteira móvel. Para obter as aproximações νhx ≈ γ1 ∂u∂x e νh ≈ γ ∂y , primeiro obtém-se o
polinómio interpolador, Zh , de uγ no espaço Shr
u
. Então, pode escrever-se
np
X
Zh = Zj ϕj .
j=1
O valor das derivadas parciais espaciais em cada ponto vértice é o resultado da média
aritmética das derivadas em cada triângulo que contém este ponto, isto é,
∂ϕj 1 X ∂ϕj
(Pm ) = (Pm ),
∂x ]Im ∂x Ti
i∈Im
Finalmente,
np np
1 X ∂ϕj 1 X ∂ϕj
νhx (Pm ) = − Zj (Pm ) e νhy (Pm ) = − Zj (Pm ).
γ j=1 ∂x γ j=1 ∂y
65
Método dos Elementos Finitos para problemas com Fronteiras Livres
O vetor normal em cada ponto da fronteira é a soma normalizada dos dois vectores normais
aos dois segmentos que concorrem nesse ponto, como mostra a gura 3.8.
Suponha-se que no ponto Pm , pertencente à fronteira móvel, o vetor normal é n = (nx , ny ),
então a equação para a velocidade da fronteira em (3.31) pode ser discretizada da forma
∂Xm ∂Ym
nx + ny = νhx (Pm )nx + νhy (Pm )ny . (3.40)
∂t ∂t
Esta equação é adicionada às equações (3.36), (3.37) e (3.38) substituindo a equação relativa
ao ponto Xm .
Como é preciso mais uma equação para o movimento da fronteira, impõe-se que os pontos
fronteiros se movam segundo retas que passam na origem. Para o ponto Pm se mover segundo
uma recta que passa na origem tem de satisfazer a equação
∂Xm ∂Ym
Ym − Xm = 0. (3.41)
∂t ∂t
Esta equação também é adicionada às equações (3.36), (3.37) e (3.38) substituindo a equação
relativa ao ponto Ym .
Para testar o novo método, simula-se primeiro o caso em que γ é constante. Neste caso, são
conhecidas soluções exatas, o que permite medir o erro. No caso γ = 2 e f = 0 a solução do
problema (3.1) é
" 2 2 # 21
1 x + y2
u(x, y, t) = 2 1 − ,
ϑ %0 ϑ
+
onde 61
%20
t
ϑ= , t0 = .
t0 6
Iniciou-se a simulação com o raio inicial %0 = 0.5 e de niu-se o tempo t = t − t0 .
O domínio computacional é o suporte da solução inicial, ou seja, Ωc = Ω(0) = supp(u0 ) =
círculo de centro na origem e raio 0.5. Utilizaram-se aproximações de u de grau 3 e de x de
grau 1. O parâmetro de suavização foi τ = 0.1 A malha inicial é uma malha estruturada com
512 triângulos, como mostra a gura 3.9.
66
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 3.9: Domínio computacional e malha inicial usados no problema com expoente constante.
Na gura 3.10, está representada a solução obtida, para alguns valores de t. Pode observar-
-se que evolui da mesma forma que a solução exata.
Figura 3.10: Evolução no tempo da aproximação obtida no problema com expoente constante.
A gura 3.11 tem representada a malha em alguns instantes. A tracejado tem-se a posição
exata da fronteira no respetivo instante. Veri ca-se que a malha se adapta ao movimento da
fronteira.
Figura 3.11: Evolução no tempo da malha obtida no problema com expoente constante.
Foram realizadas várias simulações para malhas iniciais diferentes. Em cada simulação, o
número de triângulos foi dividido por quatro. Foi medido o erro da solução na norma L2 (Ω) e o
erro da fronteira na norma L∞ . Os resultados estão apresentados na gura 3.12.
A gura tem representado o logaritmo do erro em função do logaritmo do número de triân-
gulos. Pode observar-se que tanto a solução como a fronteira convergem, mas a convergência
é lenta.
67
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 3.12: Erro das aproximações da solução e da fronteira, obtidas no problema com expoente
constante, em função do número de triângulos.
Este exemplo é o mesmo do segundo exemplo simulado para a malha xa, ou seja, com
x 2 y 2
γ= + + 1.1.
2 2
Neste caso, o domínio computacional é Ωc = Ω(0) = supp(u0 ) = coroa circular de raios 0.5 e 1.
Utilizaram-se aproximações de u de grau 2 e de x de grau 1. O parâmetro de suavização é o
mesmo do exemplo anterior.
A malha inicial foi calculada utilizando o programa initmesh do Matlab, que utiliza a trian-
gulação de Delaunay e está representa na gura 3.13. A malha é composta por 212 triângulos.
Figura 3.13: Domínio computacional e malha inicial usados no problema com suporte não convexo.
Na gura 3.14, pode observar-se que a solução obtida evolui da mesma forma que a solução
obtida no caso de malha xa. Não é possível simular a solução após o desaparecimento do
buraco interior, pois, neste caso, todos os pontos da fronteira interior convergem para um
ponto e o método torna-se degenerado. Para evitar que isso aconteça, para-se o algoritmo
quando os pontos da fronteira estão muito próximos.
Na gura 3.15, está representada a evolução da malha. Nota-se que se adapta bem ao
movimento da fronteira da solução obtida.
68
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 3.14: Evolução no tempo da aproximação obtida no problema com suporte não convexo.
Figura 3.15: Evolução no tempo da malha obtida no problema com suporte não convexo.
Agora, simulam-se dois casos de problemas de equações em meios porosos mas com absor-
ção.
Considere-se, então, o problema
∂u
∂t = div(uγ(x,t) ∇u) − uσ(x,t) em ΩT
u=0 em ΓT .
u(x, 0) = u0 (x) em Ω
Para estes valores dos expoentes, a velocidade da fronteira é conhecida e a expressão coincide
com a expressão para o caso sem absorção. Iniciou-se a simulação em t0 = 0.01 e normalizou-se
o tempo da forma t = t − t0 . O domínio computacional é o suporte da solução inicial, isto é,
Ωc = Ω(0) = supp(u0 ) = círculo centrado na origem de raio 0.7.
Utilizou-se aproximações de U de grau 4 no espaço e grau zero no tempo. A malha compu-
tacional usada foi uma malha estruturada como mostra a gura 3.11.
Na gura 3.16, mostra-se a evolução da solução obtida com uma malha de 512 triângulos.
Pode observar-se que o comportamento da solução aproximada é semelhante ao comportamento
da solução exata.
Na gura 3.17, estão representadas as malhas obtidas em alguns instantes. A posição exata
da fronteira é indicada por uma linha a tracejado. Observa-se que, mais uma vez, a malha se
69
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 3.16: Evolução no tempo da aproximação obtida no problema com absorção e expoente constante.
Figura 3.17: Evolução no tempo da malha obtida no problema com absorção e expoente constante.
O problema foi simulado para diferentes malhas com diferente número de triângulos e cal-
culou-se o erro da solução e da fronteira em cada simulação. Os resultados obtidos do erro
na norma L2 (Ω) para a solução e na norma L∞ (Ω) na fronteira foram representados na gura
3.18. Para analisar a ordem de convergência, representou-se o logaritmo do erro em função do
logaritmo do número de triângulos. Os resultados indicam que existe convergência tanto para
solução como para a fronteira, mas essa convergência é lenta.
Figura 3.18: Erro das aproximações da solução e da fronteira, obtidas no problema com absorção e
expoente constante, em função do número de triângulos.
Por m, simula-se um problema com absorção onde ambos os expoentes são variáveis. Neste
caso, não existe nenhum resultado teórico que indique o movimento da fronteira. Para expo-
entes constantes e para σ ≥ 1 em dimensão um, o movimento da fronteira é dado pela mesma
equação que o problema sem absorção. Em dimensões superiores a um, surge um termo que
70
Método dos Elementos Finitos para problemas com Fronteiras Livres
pensa-se não ter muita in uência para estes valores de σ e que pode ser desprezado. Neste
sentido, supõe-se que a mesma equação é válida para expoentes variáveis e domínios em R2.
Para se poder validar as suposições, testa-se o mesmo exemplo simulado para malha xa, ou
seja, com os expoentes
x2 + y 2
γ= e σ = x2 + y 2 + 1 + e−t ,
t2 + 1
e valores iniciais u0 de nidos por
( p p
cos(2 π x2 + y 2 ), x2 + y 2 < 0.5
u0 = .
0 , restante domínio
Figura 3.19: Domínio computacional e malha inicial usados no problema com absorção e expoente
variável.
Na gura 3.20, pode observar-se que a solução evolui da mesma forma que a solução obtida
em malha xa.
Figura 3.20: Evolução no tempo da aproximação obtida no problema com absorção e expoente variável.
Na gura 3.21, está representada a evolução da malha. Note-se que ela se comporta da
mesma forma que a malha obtida para o caso de expoentes constantes.
71
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 3.21: Evolução no tempo da malha obtida no problema com absorção e expoente variável.
72
Método dos Elementos Finitos para problemas com Fronteiras Livres
Capítulo 4
As funções u e v podem descrever as densidades de duas populações que interagem através das
funções a1 e a2 .
Este capítulo, diz respeito à prova da existência, unicidade e propriedades físicas das so-
luções globais regulares do sistema (4.1). O outro objetivo é provar a convergência de uma
solução totalmente discreta com o método dos elementos nitos de Euler-Galerkin.
Na secção 1, formula-se o problema e as hipóteses sobre os dados. Na secção 2, prova-se a
existência e unicidade de soluções globais fortes. Na secção 3, investigam-se as propriedades de
localização de soluções locais. Na secção 4, apresenta-se um resultado sobre o comportamento
assintótico da solução global do sistema. Na secção 5, discretiza-se o problema no espaço e
prova-se a convergência da solução semidiscreta. A secção 6 é dedicada à discretização no
tempo e à prova de convergência da solução totalmente discreta. Por m, apresentam-se al-
guns exemplos da aplicação desta teoria, na secção 7.
De uma forma exata, seja Ω um domínio de Rd com fronteira regular ∂Ω. Considere-se o
seguinte problema: encontrar o par (u, v) que satisfaz as seguintes condições:
ut − a1 (l1 (u), l2 (v))∆u + λ1 |u|p−2 u = f1 (x, t) em Ω×]0, T ]
v − a (l (u), l (v))∆v + λ |v|p−2 v = f (x, t) em
t 2 1 2 2 2 Ω×]0, T ]
, (4.2)
u(x, t) = v(x, t) = 0 em ∂Ω×]0, T ]
em
u(x, 0) = u (x), v(x, 0) = v (x)
0 0 Ω
H1: v0 , u0 ∈ L2 (Ω),
73
Método dos Elementos Finitos para problemas com Fronteiras Livres
H4: li : L2 (Ω) → R é uma forma linear contínua, isto é, existe gi ∈ L2(Ω), tal que:
Z
li (u) = lgi (u) = gi (x)u(x)dx, para todo u ∈ L2 (Ω), i = 1, 2,
Ω
De nição 4.1 (Solução Fraca) Diz-se que o par (u, v) é uma solução fraca do problema (4.2) se
(ii) para todo ζ ∈ H01 (Ω), as seguintes igualdades em D0 (0, T ) são verdadeiras:
Z Z Z Z
d p−2
uζ dx + a1 (l1 (u), l2 (v)) ∇u.∇ζ dx + λ1 |u| uζ dx = f1 ζ dx, (4.3)
dt Ω Ω Ω Ω
e
Z Z Z Z
d p−2
vζ dx + a2 (l1 (u), l2 (v)) ∇v.∇ζ dx + λ2 |v| vζ dx = f2 ζ dx, (4.4)
dt Ω Ω Ω Ω
Teorema 4.2 (Existência) Seja p > 1 e 0 < T < +∞. Se as hipóteses H1-H5 forem satisfeitas,
então existe uma solução fraca (u, v) do sistema (4.2) no sentido da de nição 4.1.
Demonstração
Sejam {χn (x)}n∈N a base Hilbertiana de H01 (Ω) e Sn o espaço gerado por χ1, χ2 , . . . , χn ,
n = 1, 2, . . . .
Considere-se
n
X n
X
un (x, t) = Uin (t)χi (x) e vn (x, t) = Vin (t)χi (x),
i=1 i=1
as soluções fracas do sistema aproximado correspondente a (4.2), onde Uin (t) e Vin (t) são as
soluções do sistema ODE, não linear, na variável t, seguinte:
Z Z Z Z
p−2
(un )t χ dx+a1 (l1 (un ), l2 (vn )) ∇un .∇χ dx+λ1 |u| uχ dx = f1 χ dx ∀χ ∈ Sn , (4.5)
Ω Ω Ω Ω
Z Z Z Z
p−2
(vn )t χ dx+a2 (l1 (un ), l2 (vn )) ∇vn .∇χ dx+λ2 |v| vχ dx = f2 χ dx ∀χ ∈ Sn , (4.6)
Ω Ω Ω Ω
74
Método dos Elementos Finitos para problemas com Fronteiras Livres
n
X
vn (0) = v0n = V0in χi → v0 em L2 (Ω). (4.8)
i=1
Como é bem sabido, o sistema (4.5)-(4.8) tem uma solução local (un (t), vn (t)) em algum
intervalo [0, tn [, 0 < tn < T . Para provar que essa solução pode ser estendida a todo o intervalo
]0, T ], para qualquer T > 0, usa-se a primeira estimativa à priori calculada abaixo.
Fazendo χ = un (x, t) em (4.5), obtém-se
Z
1 d
kun k2L2 (Ω) + a1 (l1 (un ), l2 (vn ))kun k2H 1 (Ω) + λ1 kun kpLp (Ω) = f1 un dx. (4.9)
2 dt 0
Ω
a−
Z
|f1 un | dx ≤ Ckf1 k2L2 (Ω) + kun k2H 1 (Ω) .
Ω 2 0
1 d a−
kun k2L2 (Ω) + kun k2H 1 (Ω) + λ1 kun kpLp (Ω) ≤ Ckf1 k2L2 (Ω) .
2 dt 2 0
d
kun k2L2 (Ω) + a− kun k2H 1 (Ω) + 2λ1 kun kpLp (Ω) ≤ Ckf1 k2L2 (Ω) . (4.10)
dt 0
Então, pode estender-se a solução ao intervalo ]0, T ]. De seguida, passa-se ao limite quando
n → ∞. Tem-se que
p−2
(un )t = a1 (l1 (un ), l2 (vn ))∆un − λ1 |un | un + f1 ∈ H −1 (Ω).
Note-se que −a1 (l1 (un ), l2 (vn ))∆un de ne um elemento de H −1 (Ω), dado por
Z
h−a1 (l1 (un ), l2 (vn ))∆un , ζi = a1 (l1 (un ), l2 (vn )) ∇un · ∇ζdx, ∀ζ ∈ H01 (Ω).
Ω
75
Método dos Elementos Finitos para problemas com Fronteiras Livres
Então,
(un )t e (vn )t são limitadas em L2 (0, T ; H −1 (Ω)). (4.16)
Por outro lado, H01 (Ω) ,→,→ L2 (Ω) ,→ H −1 (Ω). Juntando as convergências (4.18) e (4.19) e
o lema de compacidade de Aubin-Lions, obtém-se
logo, passando, se necessário a uma subsucessão, ainda denotada por (un ), tem-se
un → u qtp em Ω×]0, T ].
p−2
Uma vez que z 7→ |z| z é uma função contínua, deduz-se
p−2 p−2
|un | un → |u| u qtp em Ω×]0, T ]. (4.21)
p−2
(|un | un ) é limitada em Lp0 (0, T ; Lp0 (Ω)). (4.22)
p−2 p−2
|un | un * |u| u em Lp0 (0, T ; Lp0 (Ω)),
76
Método dos Elementos Finitos para problemas com Fronteiras Livres
u ∈ C 0 ([0, T ]; L2 (Ω)).
Deste modo, faz sentido calcular u(0). Seja ϑ ∈ C 1 (0, T ; R), com ϑ(0) = 1 e ϑ(T ) = 0.
Devido à convergência (4.19) tem-se que
Z T Z T
(u0n , ζ)ϑdt → (u0 , ζ)ϑdt, ζ ∈ L2 (Ω). (4.23)
0 0
Utilizando a convergência (4.18) em (4.24), obtém-se que (un (0), ζ) → (u(0), ζ), para todo
ζ ∈ H01 (Ω). Mas un (0) converge fortemente para u0 em L2 (Ω), consequentemente, fracamente
em L2 (Ω). Além disso, (un (0), ζ) → (u0 , ζ), para todo ζ ∈ H01 (Ω). Pela unicidade do limite,
tem-se que (u(0), ζ) → (u0 , ζ), para todo ζ ∈ H01 (Ω). Logo u(0) = u0 . Da mesma forma,
conclui-se que v(0) = v0 . Então, o problema (4.2) tem uma solução.
Teorema 4.3 (Unicidade) Seja p > 1 e 0 < T < +∞. Se as hipóteses H1-H5 forem satisfeitas,
então o problema (4.2) tem, no máximo, uma solução.
Demonstração
Z Z Z Z
d
p−2
dt u2 ζ dx + a1 (l1 (u2 ), l2 (v2 )) ∇u2 .∇ζ dx + λ1 |u2 | u2 ζ dx = f1 ζ dx
ZΩ ZΩ ZΩ ZΩ .
d p−2
v2 ζ dx + a2 (l1 (u2 ), l2 (v2 )) ∇v2 .∇ζ dx + λ2 |v2 | v2 ζ dx = f2 ζ dx
dt Ω Ω Ω Ω
77
Método dos Elementos Finitos para problemas com Fronteiras Livres
Usou-se o facto de
Z Z
a1 (l1 (u1 ), l2 (v1 )) ∇u1 · ∇ζ1 dx − a1 (l1 (u2 ), l2 (v2 )) ∇u2 · ∇ζ1 dx =
Ω Ω
Z Z
= a1 (l1 (u1 ), l2 (v1 )) ∇u1 · ∇ζ1 dx − a1 (l1 (u2 ), l2 (v2 )) ∇u2 · ∇ζ1 dx+
Ω Ω
Z Z
+a1 (l1 (u1 ), l2 (v1 )) ∇u2 · ∇ζ1 dx − a1 (l1 (u1 ), l2 (v1 )) ∇u2 · ∇ζ1 dx =
Ω Ω
Z Z
= a1 (l1 (u1 ), l2 (v1 )) ∇ζ1 · ∇ζ1 dx − (a1 (l1 (u2 ), l2 (v2 )) − a1 (l1 (u1 ), l2 (v1 ))) ∇u2 · ∇ζ1 dx.
Ω Ω
Somando e usando o facto de (|u1 |p−2 u1 − |u2 |p−2 u2 , u1 − u2 ) ≥ 0 e (|v1 |p−2 v1 − |v2 |p−2 v2 , v1 − v2 ) ≥ 0,
deduz-se que
1 d
(kζ1 k2L2 (Ω) + kζ2 k2L2 (Ω) ) + a1 (l1 (u1 ), l2 (v1 ))kζ1 k2H 1 (Ω) + a2 (l1 (u1 ), l2 (v1 ))kζ2 k2H 1 (Ω) ≤
2 dt 0 0
Z
≤ |a1 (l1 (u2 ), l2 (v2 )) − a1 (l1 (u1 ), l2 (v1 ))| ∇u2 · ∇ζ1 dx+
Ω
Z
+|a2 (l1 (u2 ), l2 (v2 )) − a2 (l1 (u1 ), l2 (v1 ))| ∇v2 · ∇ζ2 dx.
Ω
1 d
(kζ1 k2L2 (Ω) + kζ2 k2L2 (Ω) ) + a− (kζ1 k2H 1 (Ω) + kζ2 k2H 1 (Ω) ) ≤
2 dt 0 0
≤ (A1 |l1 (u2 ) − l1 (u1 )| + B1 |l2 (v2 ) − l2 (v1 )|)ku2 kH01 (Ω) kζ1 kH01 (Ω) +
+(A2 |l1 (u2 ) − l1 (u1 )| + B2 |l2 (v2 ) − l2 (v1 )|)kv2 kH01 (Ω) kζ2 kH01 (Ω)
≤ (A1 |l1 (u2 − u1 )| + B1 |l2 (v2 − v1 )|)ku2 kH01 (Ω) kζ1 kH01 (Ω) +
+(A2 |l1 (u2 − u1 )| + B2 |l2 (v2 − v1 )|)kv2 kH01 (Ω) kζ2 kH01 (Ω)
≤ C1 (ku2 − u1 kL2 (Ω) + kv2 − v1 kL2 (Ω) )ku2 kH01 (Ω) kζ1 kH01 (Ω) +
+C2 (ku2 − u1 kL2 (Ω) + kv2 − v1 kL2 (Ω) )kv2 kH01 (Ω) kζ2 kH01 (Ω) .
Logo,
1 d a− M (t)
(kζ1 k2L2 (Ω) + kζ2 k2L2 (Ω) ) + (kζ1 k2H 1 (Ω) + kζ2 k2H 1 (Ω) ) ≤ (kζ1 k2L2 (Ω) + kζ2 k2L2 (Ω) ),
2 dt 2 0 0 2
78
Método dos Elementos Finitos para problemas com Fronteiras Livres
com
C12 C2
M (t) = −
ku2 k2H 1 (Ω) + 2− kv2 k2H 1 (Ω) ∈ L1 (0, T ).
4a 0 4a 0
De nindo z(t) = (kζ1 k2L2 (Ω) + kζ2 k2L2 (Ω) ), tem-se que
dz
≤ M (t)z, z(0) = 0.
dt
Demonstração
Seja k > 0, se se multiplicar a primeira equação do sistema (4.2) por u2k−1 e se for integrada
em Ω, obtém-se
Z Z Z Z
1 d 2k 2 2(k−1) p+2k−2
u dx + (2k − 1)a1 (l1 (u), l2 (v)) |∇u| u dx + λ1 |u| dx = f1 u2k−1 dx.
2k dt Ω Ω Ω Ω
Então,
1 d 2k−1
kuk2k
L2k (Ω) ≤ kf kL2k (Ω) kukL2k (Ω) ,
2k dt
ou, de modo equivalente,
d
kuk2k−1
L2k (Ω) kukL2k (Ω) ≤ kf kL2k (Ω) kuk2k−1
L2k (Ω) .
dt
Z t
kukL2k (Ω) ≤ ku0 kL2k (Ω) + kf kL2k (Ω) dt.
0
Demonstração
Multiplicando a primeira equação do sistema (4.2) por ∆u e integrando em Ω, resulta
Z Z Z Z
ut ∆u dx − a1 (l1 (u), l2 (v))(∆u)2 dx + λ1 |u|p−2 u∆u dx = f1 ∆u dx.
Ω Ω Ω Ω
79
Método dos Elementos Finitos para problemas com Fronteiras Livres
por isso,
Z Z Z Z
1 d
|∇u|2 dx + a1 (l1 (u), l2 (v))(∆u)2 dx + λ1 (p − 1) |u|p−2 |∇u|2 dx = − f1 ∆u dx.
2 dt Ω Ω Ω Ω
Simpli cando
!
Z Z T Z Z Z T Z
2 − 2 2
|∇u(x, T )| dx + a (∆u) dxdt ≤ C |∇u0 | dx + f12 dxdt .
Ω 0 Ω Ω 0 Ω
e
vt = a2 (l1 (u), l2 (v))∆v − λ2 |v|p−2 v + f2 ∈ L2 (0, T ; L2 (Ω)).
De seguida, estudam-se alguns efeitos de localização das soluções deste tipo de sistemas.
No decorrer desta secção, supõe-se sempre que as hipóteses dos teoremas de existência e
unicidade são satisfeitas. Sejam x0 ∈ Ω e %0 ∈]0, dist(x0 , ∂Ω)], então de ne-se
B% ≡ B% (x0 ) = x ∈ Rd : |x − x0| < % ⊂ Ω, e
80
Método dos Elementos Finitos para problemas com Fronteiras Livres
Assume-se que
u(x, 0) = 0 e v(x, 0) = 0 qtp em B%0 . (4.25)
De nição 4.6 ([3],119) A função %(t) : [0, t∗ [→ [0, +∞[, %(0) ≤ %0 , é dita uma taxa no ponto
x0 se
∀t ∈ [0, t∗ [, u(x, t) = 0 qtp em B%(t) (x0 ) = {x : |x − x0 | < %(t)} ⊂ Ω.
De nição 4.7 ([3],120) Diz-se que a função u(x, t) possui a propriedade de localização estável
se para algum x0 ∈ Ω, existe uma taxa estritamente positiva %(t), no ponto x0 , de nida no
conjunto [0, ∞[ de tal modo que
lim inf %(t) > 0.
t→∞
Z tZ
Eu (%, t) = |∇u|2 dxdς, E u (%) = sup Eu (%, ς),
0 B% 0≤ς≤T
Z tZ
p
Du (%, t) = |u| dxdς, e Du (%) = sup Du (%, ς).
0 B% 0≤ς≤T
As funções de energia locais para v são de nidas de forma análoga. Por sup pretende dizer-se
ess sup.
De ne-se
b = bu + bv , E = Eu + Ev e D = Du + Dv .
Teorema 4.8 (Localização estável) Seja (u, v) solução fraca de (4.2) com 1 < p < 2 em B%0 ×]0, T ],
(B%0 ⊂ Ω) e suponha-se que
Então,
u(x, t) = 0 e v(x, t) = 0 qtp em (x, t) ∈ B% × R+,
onde % é de nido pela equação (4.32).
81
Método dos Elementos Finitos para problemas com Fronteiras Livres
Demonstração
1
b(%) + a− E(%, t) + λ1 D(%, t) ≤ I. (4.28)
2
Z tZ
|I| ≤ a+ |∇u| |u| dsdς
0 Γ%
Z tZ ! 12 Z tZ ! 12
+ 2 2
≤ a |∇u| dsdς |u| dsdς
0 Γ% 0 Γ%
Z tZ ! 21
1
2
≤ a+ (E% ) 2
|u| dsdς . (4.29)
0 Γ%
com
d(2 − q) + q 2−q
ϑ= < 1, 1 < q < 2, $ = − 1 + d ,
d(2 − q) + 2q 2q
e usa-se o facto de
2(q−p) p(2−q) 2−p p 4
kukLq (B% ) ≤kukLq(2−p)
2 (B% )
kukLq(2−p)
p (B% )
= kukL22(B% ) kukL2 p (B% ) , p < q = < 2.
4−p
e (4.30) ca da forma
1 2−p 1
ϑ 2−p 1 1−ϑ
kukL2 (Γ% ) ≤ C Et2 +%−$ b 4 Dt2 b 4 Dt2 . (4.31)
82
Método dos Elementos Finitos para problemas com Fronteiras Livres
Z tZ ! 21
1
2
|I| ≤ a+ (E% ) 2
|u| dsdς
0 Γ%
1
Z t
1 2−p
2ϑ 2−p 1 2(1−ϑ) 12
1
−$
≤ C (E% ) 2
Et +%
2
b 4 Dt 2
b 4 Dt2 dς
0
Z t 12
1 (2−p)(1−ϑ)
−$ϑ ϑ (1−ϑ)
≤ % K1 (E% ) 2
(b) 4 (Et +Dt ) (Dt ) dς
0
Z t 12
1 (2−p)(1−ϑ)
−$ϑ
≤ % K1 (E% ) 2
(b) 4 (Et + Dt ) dς
0
1 (2−p)(1−ϑ) 1
≤ %−$ϑ K1 (E% ) 2 ( b ) 4 (E+ D) 2
1 12 + (2−p)(1−ϑ)
≤ %−$ϑ K1 (E% ) 2 b+E+ D 4
,
(2−p)ϑ
com K1 = Cmax(1, ( b ) 4
0 ). Substituindo a última desigualdade em (4.28),
(%0 )) max(1, %$ϑ
obtém-se
21 12 + (2−p)(1−ϑ)
b + E + D ≤ %−$ϑ K1 E % b+E + D 4
,
12 − (2−p)(1−ϑ) 21
b+E+D 4
≤ %−$ϑ K1 E % .
Se se de nir
(2 − p)(1 − ϑ)
℘= < 1,
2
então tem-se que
( E )1−℘ ≤ %−2$ϑ K12 E % ,
De modo semelhante,
(E v )1−℘ ≤ %−2$ϑ K22 (E v )% .
Kω
%ω ω
(E(%0 ))℘ − (E(%))℘ , ω = 1 + 2$ϑ.
0 −% ≤
℘
Kω
%ω (t) = %ω
0 − (E(%0 ))℘ . (4.32)
℘
83
Método dos Elementos Finitos para problemas com Fronteiras Livres
não for muito elevada, então %(t) > 0 e existe um cilindro B% × R+ onde u e v são nulas qtp.
Seguidamente, investiga-se o efeito de tempo de espera no sentido da de nição seguinte.
Diz-se que u(x, t) possui a propriedade de tempo de espera generalizada se para algum tw > 0,
a função %(t) ≡ %0 é uma taxa no ponto x0 no intervalo [0, tw ].
Assume-se que
Z Z
2 2
|u0 | dx = 0, |v0 | dx = 0, f1 (x, t) = 0, f2 (x, t) = 0, (x, t) ∈ B%0 × [0, T ], (4.33)
B%0 B% 0
e os parâmetros CR , ℘, R , %0 satisfazem
1
1−℘
CR ℘
− 1−℘ 1
G = CR − C5 (R − %0 ) − C6 (R − %0 ) ℘ = 0, (4.36)
℘
Teorema 4.10 (Efeito de tempo de espera) Seja (u, v) solução fraca de 4.2 com 1 ≤ p < 2 em
BR ×]0, T ]. Se (4.33)-(4.36) são satisfeitas então
Demonstração
Repetindo os primeiros passos da demonstração do teorema 4.8 para %0 ≤ % ≤ R chega-se
às desigualdades
1 1
bu + Eu (%, t) + Du (%, t) ≤ C1 (Eu )%1−℘ + C2 [% − %0 ]+℘ (4.38)
e
1 1
bv + Ev (%, t) + Dv (%, t) ≤ C3 (Ev )%1−℘ + C4 [% − %0 ]+℘ . (4.39)
84
Método dos Elementos Finitos para problemas com Fronteiras Livres
Utilizando as propriedades
Z tZ
∇u + ∇v + |u|p + |v|p dxdς,
2 2
E% ≤ Ψ% , Ψ=E+D=
0 B%
Z tZ !
∂ ∇u + ∇v + |u|p + |v|p dsdς
2 2
sup Ψ% (%, ς) = sup =
ς∈[0,t] ∂% ς∈[0,t] 0 Γ%
! !
Z T Z
∇u + ∇v + |u| + |v| dsdt = ∂
p p
2 2
= sup Ψ(%, ς) ,
0 Γ% ∂% ς∈[0,t]
onde,
Z T Z
∇u + ∇v + |u|p + |v|p dxdς.
2 2
Ξ(%) = sup Ψ(%, t) =
t∈[0,T ] 0 B%
Considere-se o problema:
1 1
φ(%) = C5 φ%1−℘ + C6 [% − %0 ]+℘ , %0 ≤ % ≤ R. (4.42)
Nesta secção, estuda-se o comportamento assintótico das soluções quando t tende para
in nito. Introduz-se a função de energia global
Z
1
b(t) = u2 + v 2 dx,
2 Ω
85
Método dos Elementos Finitos para problemas com Fronteiras Livres
Demonstração
Suponha-se que f1 ≡ f2 ≡ 0 e λ1 = λ2 = 0. Multiplicando a primeira equação do sistema
(4.2) por u, a segunda por v e integrando em Ω, obtém-se
Z Z
1 d
u2 dx + a1 (l1 (u), l2 (v))|∇u|2 dx = 0
2 dt Ω Ω
e Z Z
1 d
v 2 dx + a2 (l1 (u), l2 (v))|∇v|2 dx = 0.
2 dt Ω Ω
d
b(t) + 2Ca− b(t) ≤ 0
dt
e
−
b(t) ≤ b(0)e−2Ca t .
Teorema 4.12 (Extinção em tempo nito) Suponha-se que 1 < p < 2 e (u, v) é solução de
(4.2).
2℘−1
t 2−2℘
kfi kL2 (Ω) ≤ 1− , i = 1, 2,
tf +
com tf ≥ t0 , ℘ de nidos por (4.45) e 0 < < 0 , implica que (u, v) extingue-se num tempo
nito.
Demonstração
Multiplicando a primeira equação de (4.2) por u e integrando em Ω, obtém-se
Z Z Z Z
1 d
u2 dx + a1 (l1 (u), l2 (v))|∇u|2 dx + λ1 |u|p dx = f1 u dx,
2 dt Ω Ω Ω Ω
Z Z ϑ Z 2(1−ϑ)
p
2 2 p
u dx ≤ |∇u| dx |u| dx ≤
Ω Ω Ω
86
Método dos Elementos Finitos para problemas com Fronteiras Livres
Z ϑ+ 2(1−ϑ)
p
2 p
≤ |∇u| + |u| dx , ϑ ∈]0, 1[,
Ω
onde,
p
℘= < 1. (4.45)
pϑ + 2(1 − ϑ)
Substituindo esta última desigualdade em (4.44), obtém-se
Z Z ℘ Z
1 d
u2 dx + C u2 dx ≤ f1 u dx. (4.46)
2 dt Ω Ω Ω
Integrando em t, obtém-se
Logo,
b(0)1−℘
b = 0 para t ≥ te = <∞
C(1 − ℘)
e o mesmo acontece para u e v.
Suponha-se agora que f1 6≡ 0 ou f2 6≡ 0. Neste caso, escreve-se (4.46) e (4.47) da forma
Z Z ℘ Z 21 Z 12
1 d 2
u dx + C 2
u dx ≤ f12 dx 2
u dx (4.48)
2 dt Ω Ω Ω Ω
e Z ℘ Z 21 Z 12
Z
1 d 2
v dx + C 2
v dx ≤ 2
f2 dx 2
v dx . (4.49)
2 dt Ω Ω Ω Ω
1
b0 + C1 b℘ ≤ C2 (t)b 2 ,
1
onde C2 (t) = C max{kf1 kL2 (Ω) , kf2 kL2 (Ω) }. Introduzindo a função g(t) = b 2 , pode escrever-se
esta desigualdade da forma
2℘−1
0 C1 2℘−1 t 2−2℘
g + g ≤ C2 (t) ≤ 1− .
2 2 tf +
2℘−1
0 C1 2℘−1 t 2−2℘
h + h = 1− . (4.50)
2 2 tf +
87
Método dos Elementos Finitos para problemas com Fronteiras Livres
As soluções desta equação são majorantes das funções g. Veri ca-se facilmente que se h(0),
e tf satisfazem
C1 h(0)
h(0)2℘−1 − − ≥ 0, (4.51)
2 2(1 − ℘)tf 2
então as funções
1
t 2−2℘
h = h(0) 1 −
tf +
Teorema 4.13 (Decaimento polinomial) Suponha-se agora que p > 2 e (u, v) é solução do pro-
blema (4.2). Se f1 6≡ 0 ou f2 6≡ 0, então existem K > 0 e B > 0, dependentes de d, p e Ω, tais
que:
K
kfi kL2 (Ω) ≤ 2℘−1 , i = 1, 2, (4.52)
(Bt + 1) 2(℘−1)
implica que b satisfaça
b(0)
b(t) ≤ 1 .
(Bt + 1) ℘−1
Demonstração
Multiplicando a primeira equação de (4.2) por u, e integrando em Ω, obtém-se
Z Z Z Z
1 d
u2 dx + a1 (l1 (u), l2 (v))|∇u|2 dx + λ1 |u|p dx = f1 u dx. (4.53)
2 dt Ω Ω Ω Ω
Então, Z Z ℘
p 2
|u| dx ≥ C u dx ,
Ω Ω
onde ℘ = p
2 > 1.
Substituindo esta desigualdade em (4.53) e ignorando o segundo termo resulta
Z Z ℘ Z
1 d 2 2
(4.54)
u dx + C u dx ≤ f1 u dx .
2 dt Ω Ω Ω
e Z ℘ Z 21 Z 21
Z
1 d 2
v dx + C v dx2
≤ f22 dx 2
v dx . (4.57)
2 dt Ω Ω Ω Ω
88
Método dos Elementos Finitos para problemas com Fronteiras Livres
1
b0 + C1 b℘ ≤ C2 (t)b 2 ,
C1 2℘−1 K
g0 + g ≤ 2℘−1 .
2 (Bt + 1) 2(℘−1)
C1 2℘−1 K
h0 + h = 2℘−1 . (4.58)
2 (Bt + 1) 2(℘−1)
Se as constantes satisfazem
2℘−1 2℘−1 1
K < Cb(0) 2 , B = (b(0) 2 − K)(℘ − 1)b(0)− 2 ,
então 1
b02
h= 1
(Bt + 1) 2(℘−1)
são soluções desta equação e são majorantes das funções g. Voltando à função b, conclui-se a
estimativa desejada.
Nesta secção, discretiza-se o problema no espaço e obtêm-se estimativas do erro para essa
discretização. Neste problema, considera-se apenas polinómios de grau 1. Então, Sh irá denotar
o conjunto das funções contínuas em Ω̄ que se anulam na fronteira ∂Ω e que são polinómio de
grau 1, em cada triângulo, ou seja,
89
Método dos Elementos Finitos para problemas com Fronteiras Livres
np
X np
X
uh (x, t) = αj (t)ψj (x) e vh (x, t) = βj (t)ψj (x),
j=1 j=1
tais que:
np Z np np np Z
X X X X
0
αj (t) ψi ψj dx + a1 (l1 ( αk ψk ), l2 ( βk ψk )) ∇ψi · ∇ψj dx αj =
j=1 Ω j=1 k=1 k=1 Ω
Xnp Z np
X Xnp Z
p−2
= −λ 1 | α j ψj | α j ψ j ψ i dx = ψi f1 dx, i = 1, . . . , np
Ω j=1 Ω
j=1 j=1
np np np np
(4.60)
Z Z
X
0
X X X .
β j (t) ψ ψ
i j dx + a (l
2 1 ( α ψ
k k ), l2 ( β ψ
k k )) ∇ψi · ∇ψj dx βj =
Ω Ω
j=1 j=1 k=1 k=1
np Z np np Z
X X X
p−2
= −λ2 | βj ψj | βj ψj ψi dx = ψi f2 dx, i = 1, . . . , np
j=1 Ω j=1 Ω
j=1
α (0) = α0 = u (P ), β (0) = β 0 = v (P ), j = 1, . . . , np
j j 0 j j j 0 j
Z np
X np
X Z
Fiv = −λ2 | βj ψj |p−2 βj ψj ψi dx + ψi f2 dx.
Ω j=1 j=1 Ω
90
Método dos Elementos Finitos para problemas com Fronteiras Livres
Teorema 4.14 Seja p ≥ 2. Se (u, v) é solução do problema (4.2) e (uh , vh ) é solução do pro-
blema (4.59), supondo que u e v são su cientemente regulares, então
Demonstração
Antes de mais, saliente-se que a solução (u, v) de (4.2) também satisfaz (4.59).
Escrevendo uh − u = uh − Rh u + Rh u − u = θ + η e vh − v = vh − Rh v + Rh v − v = ρ + µ e
pelo lema 2.84, pode concluir-se que
kηkL2 (Ω) ≤ Ch2 |u|H 2 (Ω) e kµkL2 (Ω) ≤ Ch2 |v|H 2 (Ω) .
≤ C(kθk2L2 (Ω) +kηk2L2 (Ω) )+Ckθk2L2 (Ω) +Ckηt k2L2 (Ω) +Ckθk2L2 (Ω) +
+C(kθk2L2 (Ω) +kηk2L2 (Ω) + kρk2L2 (Ω) +kµk2L2 (Ω) ) + a− k∇θk2L2 (Ω) .
e
|a1 − a1h | ≤ A1 |l1 (u) − l1 (uh )| + B1 |l2 (v) − l2 (vh )| = A1 |l1 (u − uh )| + B1 |l2 (v − vh )|
91
Método dos Elementos Finitos para problemas com Fronteiras Livres
d
kθk2L2 (Ω) ≤ C(kθk2L2 (Ω) + kηk2L2 (Ω) + kηt k2L2 (Ω) + kρk2L2 (Ω) + kµk2L2 (Ω) ).
dt
d
kρk2L2 (Ω) ≤ C(kρk2L2 (Ω) + kµk2L2 (Ω) + kµt k2L2 (Ω) + kθk2L2 (Ω) + kηk2L2 (Ω) ).
dt
d
(kθk2L2 (Ω) + kρk2L2 (Ω) ) ≤ C(kθk2L2 (Ω) + kρk2L2 (Ω) + kηk2L2 (Ω) + kµk2L2 (Ω) + kηt k2L2 (Ω) + kµt k2L2 (Ω) ).
dt
Z T
≤ kθ(x, 0)k2L2 (Ω) + kρ(x, 0)k2L2 (Ω) +C kηk2L2 (Ω) + kµk2L2 (Ω) + kηt k2L2 (Ω) + kµt k2L2 (Ω) dt.
0
Pelo Lema 2.79, tem-se que kθ(x, 0)kL2 (Ω) ≤ Ch2 |u0 |H 2 (Ω) e kρ(x, 0)kL2 (Ω) ≤ Ch2 |v0 |H 2 (Ω) .
Recordando as estimativas de η e µ e como kηt kL2 (Ω) = kΠSh ut − ut kL2 (Ω) ≤ Ch2 |ut |H 2 (Ω) e
kµt kL2 (Ω) = kRh vt − vt kL2 (Ω) ≤ Ch2 |vt |H 2 , pode inferir-se que
Como se considera que não há possibilidade de confusão, para simpli car a notação, omi-
te-se a dependência de h nesta secção. Seja δ > 0 e considere-se a partição ]0, T ] = ∪ni
n=1 In ,
In =]tn−1 , tn ], com tn = tn−1 + δ. A discretização no tempo é feita utilizando o método de Euler
regressivo. Neste método, a derivada temporal, ut (x, tn ), é aproximada em cada instante tn
pela diferença nita
¯ u(x, tn ) − u(x, tn−1 )
∂u(x, tn ) = .
δ
92
Método dos Elementos Finitos para problemas com Fronteiras Livres
Pnp Pnp
Utilizando outra vez a base {ψj }np
j=1 e de nindo Un = j=1 αjn ψj e Vn = j=1 βjn ψj , o
sistema (4.62) pode ser escrito em forma matricial como
O sistema de equações algébricas (4.63) é não linear e tem de ser resolvido em cada passo
no tempo. A não linearidade está na presença de ωn nos termos A(ωn ) e F(ωn ). Para se evitar
a aplicação de um método iterativo, e consequente aumento do custo computacional e do erro,
considera-se uma linearização do sistema substituindo ωn por ωn−1 , nesses dois termos. Então,
o problema totalmente discreto é calcular Un e Vn , tais que:
Z Z Z
¯
∂Un χ dx + a1 (l1 (Un−1 ), l2 (Vn−1 )) ∇Un · ∇χ dx + λ1 |Un−1 |p−2 Un−1 χ dx =
Ω Ω Ω
Z
= f1 (x, tn )χ dx, ∀χ ∈ Sh , n = 1, . . . , ni
Ω
Z Z Z
¯
∂Vn χ dx + a2 (l1 (Un−1 ), l2 (Vn−1 )) ∇Vn · ∇χ dx + λ2 |Vn−1 |p−2 Vn−1 χ dx = ,
Ω Ω Ω
Z
= f2 (x, tn )χ dx, ∀χ ∈ Sh , n = 1, . . . , ni
Ω
P1 (Sh ) u0 , V0 = ΠP1 (Sh ) v0
U =Π
0
(4.64)
o que é equivalente a calcular o vetor ωn = (α1 , . . . , αnp , β1 , . . . , βnp ) solução do seguinte sis-
n n n n
tema linear:
(M + δA(ωn−1 ))ωn = δF(ωn−1 ) + Mωn−1 , n = 1, . . . , ni. (4.65)
Este sistema possui sempre solução, uma vez que M + δA(ωn−1 ) é de nida positiva. Esta
secção termina com a estimativa do erro para a aproximação totalmente discreta da solução do
problema inicial.
Teorema 4.15 Seja p ≥ 2. Se (u, v) é solução do problema (4.2) e (Un , Vn ) solução de (4.64),
93
Método dos Elementos Finitos para problemas com Fronteiras Livres
kUn (x) − u(x, tn )k2L2 (Ω) + kVn (x) − v(x, tn )k2L2 (Ω) ≤ C(h2 + δ)2 , n = 1, . . . , ni,
Demonstração
Denote-se por un = u(x, tn ) e vn = v(x, tn ). Da mesma forma como anteriormente, escreve-
-se
Un − un = (Un − Rh un ) + (Rh un − un ) = θn + ηn
e
Vn − vn = (Vn − Rh vn ) + (Rh vn − vn ) = ρn + µn .
De seguida, estima-se θn e ρn .
Aplicando as propriedades de Rh , tem-se sucessivamente
Z Z
¯ n χ dx + a1 (l1 (Un−1 ), l2 (Vn−1 ))
∂θ ∇θn · ∇χ dx =
Ω Ω
Z Z
= ¯ n χ dx + a1 (l1 (Un−1 ), l2 (Vn−1 ))
∂U ∇Un · ∇χ dx−
Ω Ω
Z Z
− ¯ h un χ dx − a1 (l1 (Un−1 ), l2 (Vn−1 ))
∂R ∇Rh un · ∇χ dx =
Ω Ω
Z Z Z
p−2
= −λ1 |Un−1 | Un−1 χ dx + f1 (x, tn )χ dx − a1 (l1 (un ), l2 (vn )) ∇Rh un · ∇χ dx−
Ω Ω Ω
Z Z
− (un )t χ dx + (a1 (l1 (Un−1 ), l2 (Vn−1 )) − a1 (l1 (un ), l2 (vn ))) ∇Rh un · ∇χ dx−
Ω Ω
Z
− ¯ h un − (un )t ) χ dx =
(∂R
Ω
Z Z Z
= λ1 (|un |p−2 un − |Un−1 |p−2 Un−1 ) χ dx − ¯ n χ dx −
∂η ¯ n − (un )t ) χ dx+
(∂u
Ω Ω Ω
Z
+(a1 (l1 (Un−1 ), l2 (Vn−1 )) − a1 (l1 (un ), l2 (vn ))) ∇Rh un · ∇χ dx.
Ω
¯ n kL (Ω) + k∂u
+k∂η ¯ n − (un )t kL (Ω) )k∇θn kL (Ω) .
2 2 2
Na desigualdade acima, usa-se a desigualdade de Poincaré kθn kL2 (Ω) ≤ Ck∇θn kL2 (Ω) e as duas
desigualdades seguintes:
¯ n kL (Ω) ),
kUn−1 − un kL2 (Ω) ≤ C(kθn−1 kL2 (Ω) + kηn−1 kL2 (Ω) + δk∂u 2
¯ n kL (Ω) ).
kVn−1 − vn kL2 (Ω) ≤ C(kρn−1 kL2 (Ω) + kµn−1 kL2 (Ω) + δk∂v 2
94
Método dos Elementos Finitos para problemas com Fronteiras Livres
Então, após simpli car o termo a− k∇θn k2L2 (Ω) com a desigualdade de Cauchy, chega-se a
¯ n k2 2 ¯
∂kθ L2 (Ω) ≤ Ckθn−1 kL2 (Ω) + C(kηn−1 kL2 (Ω) + δk∂un kL2 (Ω) + kρn−1 kL2 (Ω) + kµn−1 kL2 (Ω) +
¯ n kL (Ω) + k∂η
+δk∂v ¯ n kL (Ω) + k∂u
¯ n − (un )t kL (Ω) )2 .
2 2 2
Utilizando o lema 2.79, conclui-se que kηn kL2 (Ω) ≤ C(u)h2 e kµn kL2 (Ω) ≤ C(v)h2 . Tem-se
também que
Z tj
¯
−1
k∂ηn kL2 (Ω) =
δ (η(x, ς))t dς
≤ C(u)h2
tj−1
L2 (Ω)
e
Z tj
¯
−1
k∂µn kL2 (Ω) =
δ (µ(x, ς))t dς
≤ C(v)h2 .
tj−1
L2 (Ω)
e
Z tj
¯
−1
k∂vn − (vn )t kL2 (Ω) =
δ (s − tj−1 )vtt (ς) dς
≤ C(v)δ.
tj−1
L2 (Ω)
Então, obtém-se
¯ n k2 2 2 2
∂kθ L2 (Ω) ≤ Ckθn−1 kL2 (Ω) + C(h + δ) .
¯ n k2 2 2 2
∂kρ L2 (Ω) ≤ Ckρn−1 kL2 (Ω) + C(h + δ) .
¯ n k2 2 2 2 2 2
∂(kθ L2 (Ω) + kρn kL2 (Ω) ) ≤ C(kθn−1 kL2 (Ω) + kρn−1 kL2 (Ω) ) + C(h + δ) .
Mostrou-se que
kθn k2L2 (Ω) + kρn k2L2 (Ω) ≤ (1 + Cδ)(kθn−1 k2L2 (Ω) + kρn−1 k2L2 (Ω) ) + Cδ(h2 + δ)2 .
kθn k2L2 (Ω) + kρn k2L2 (Ω) ≤ C(kθ0 k2L2 (Ω) + kρ0 k2L2 (Ω) ) + C(h2 + δ)2 .
Adicionado as estimativas de kηn kL2 (Ω) e kµn kL2 (Ω) , termina-se a demonstração.
95
Método dos Elementos Finitos para problemas com Fronteiras Livres
Para testar o algoritmo numérico, simula-se um exemplo em que a solução inicial está
localizada apenas em parte do domínio. Considere-se o sistema (4.2) em Ω =]0, 1[2 com p = 3 e
Os coe cientes λ1 e λ2 são ambos iguais a um. As fontes externas são dadas por
1 −x−y 1
f1 (x, y, t) = 0.1 e e f2 (x, y, t) = 0.1(x − y) .
(t + 1)10 (t + 1)7
Os parâmetros do método dos elementos nitos são h = 0.05, δ = 10−3 e a malha está
representada na gura 4.3.
Na gura 4.1, mostra-se a evolução no tempo para estimativa de u obtida e, na gura 4.2,
mostra-se a evolução no tempo para a estimativa de v. A solução expande-se a todo o domínio
e tende para zero. O decaimento é mais evidente na gura 4.4 onde se representa a função
b(t).
96
Método dos Elementos Finitos para problemas com Fronteiras Livres
97
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 4.3: Malha espacial para h = 0.05 usada no Figura 4.4: Função de energia b(t) no problema da
problema da localização da solução. localização da solução.
98
Método dos Elementos Finitos para problemas com Fronteiras Livres
Numa secção anterior, foi demonstrado que as soluções do sistema (4.2) com 1 < p < 2 têm
a propriedade da extinção em tempo nito se as funções fi satis zerem algumas condições e
têm a propriedade de decaimento polinomial se p ≥ 2 e fi decaírem de forma adequada. No
sentido de estudar estas propriedades, simula-se um exemplo com p = 1.6 e depois com p = 3
mantendo iguais todos os outros parâmetros. Embora não tenha sido provada a convergência
para p = 1.6, os resultados estão de acordo com a teoria.
Considere-se o sistema (4.2) em Ω =]0, 1[2 com
Figura 4.5: Evolução no tempo da aproximação de u obtida no problema com extinção em tempo nito e
decaimento polinomial.
Na gura 4.5, mostra-se a evolução no tempo para a estimativa de u e na gura 4.6 mostra-se
a evolução no tempo para a estimativa de v obtidas para p = 1.6. Como era de esperar, as
funções u e v são zero para t ≥ 0.065. O comportamento da solução para p = 3 é semelhante ao
comportamento da solução com p = 1.6, excepto o facto de o valor de u e v ser pequeno para
t ≥ 0.065, mas não ser zero.
De forma a veri car o comportamento assintótico, representou-se na gura 4.8 a função
de energia b, com uma escala logarítmica no eixo das ordenadas. Pela observação da gura
torna-se evidente que b se extingue em t ≈ 0.065 para p = 1.6 e, no caso p = 3, b tende para
zero mas não se extingue.
99
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 4.6: Evolução no tempo da aproximação de v obtida no problema com extinção em tempo nito e
decaimento polinomial.
Figura 4.7: Malha espacial para h = 0.025 no Figura 4.8: Função de energia b(t) no problema com
problema com extinção em tempo nito e extinção em tempo nito e decaimento polinomial.
decaimento polinomial.
100
Método dos Elementos Finitos para problemas com Fronteiras Livres
Para último exemplo, considere-se o sistema (4.2) em Ω =] − 1, 1[2 com p = 3, mas, neste
caso, os termos de difusão são da forma
1 1 2 1
a1 (l1 , l2 ) = 2 + − , a2 (l1 , l2 ) = 3 + − .
1 + l22 1 + l12 1 + l12 1 + l22
Estas funções não satisfazem as condições dos teoremas, mas é interessante ver o comporta-
mento das soluções. Os parâmetros do método dos elementos nitos e malha são os mesmos do
exemplo anterior.
Da mesma forma e como no exemplo anterior, na gura 4.9, mostra-se a evolução, no
tempo, da estimativa obtida para u e, na gura 4.10, mostra-se a evolução da estimativa
obtida para v nos mesmos instantes. Veri ca-se que a solução u decresce e, uma vez que já
estava inicialmente de nida em todo o domínio, não se veri ca nenhuma expansão. A solução v
tem um comportamento completamente diferente. Inicialmente, expande-se a todo o domínio
e decresce, mas, a partir de certo instante, começa a crescer inde nidamente. Este facto
pode ser justi cado pela expressão de f2 que é crescente no tempo e o seu efeito começa ser
dominante a partir de certo instante.
101
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 4.9: Evolução no tempo da aproximação de u obtida no problema com fontes externas ilimitadas.
102
Método dos Elementos Finitos para problemas com Fronteiras Livres
Figura 4.10: Evolução no tempo da aproximação de v obtida no problema com fontes externas ilimitadas.
103
Método dos Elementos Finitos para problemas com Fronteiras Livres
104
Método dos Elementos Finitos para problemas com Fronteiras Livres
Capítulo 5
105
Método dos Elementos Finitos para problemas com Fronteiras Livres
106
Método dos Elementos Finitos para problemas com Fronteiras Livres
Bibliogra a
[1] A. S. Ackleh and L. Ke. Existence-uniqueness and long time behavior for a class of non-
local nonlinear parabolic evolution equations. Proc. Amer. Math. Soc., 128(12):3483 3492
(electronic), 2000. 13
[2] S. Angenent. Analyticity of the interface of the porous media equation after the waiting
time. Proc. Amer. Math. Soc., 102(2):329 336, 1988. 4
[3] S. N. Antontsev, J. I. Díaz, and S. I. Shmarev. Energy methods for free boundary problems.
Progress in Nonlinear Differential Equations and their Applications, 48. Birkhäuser Boston
Inc., Boston, MA, 2002. Applications to nonlinear PDEs and uid mechanics. 81, 84
[4] S. N. Antontsev and S. I. Shmarev. A model porous medium equation with variable exponent
of nonlinearity: existence, uniqueness and localization properties of solutions. Nonlinear
Anal., 60(3):515 545, 2005. 6, 7, 9, 35, 41
[5] S. N. Antontsev and S. I. Shmarev. Existence and uniqueness of solutions of degenerate pa-
rabolic equations with variable exponents of nonlinearity. Fundam. Prikl. Mat., 12(4):3 19,
2006. 6
[6] D. G. Aronson. Regularity propeties of ows through porous media. SIAM J. Appl. Math.,
17:461 467, 1969. 3
[7] D. G. Aronson. Regularity properties of ows through porous media: The interface. Arch.
Rational Mech. Anal., 37:1 10, 1970. 4
[8] D. G. Aronson and J. L. Vázquez. Eventual C ∞ -regularity and concavity for ows in one-
dimensional porous media. Arch. Rational Mech. Anal., 99(4):329 348, 1987. 4
[9] M. J. Baines. Moving nite elements. Monographs on Numerical Analysis. The Clarendon
Press Oxford University Press, New York, 1994. Oxford Science Publications. 31
[10] M. J. Baines, M. E. Hubbard, and P. K. Jimack. A moving mesh nite element algorithm
for uid ow problems with moving boundaries. Internat. J. Numer. Methods Fluids,
47(10-11):1077 1083, 2005. 8th ICFD Conference on Numerical Methods for Fluid Dyna-
mics. Part 2. 11
[11] M. J. Baines, M. E. Hubbard, and P. K. Jimack. A moving mesh nite element algorithm
for the adaptive solution of time-dependent partial differential equations with moving
boundaries. Appl. Numer. Math., 54(3-4):450 469, 2005. 11
[13] G. I. Barenblatt. On some unsteady motions of a liquid and gas in a porous medium. Akad.
Nauk SSSR. Prikl. Mat. Meh., 16:67 78, 1952. 5
107
Método dos Elementos Finitos para problemas com Fronteiras Livres
[15] L. A. Caffarelli and N. I. Wolanski. C 1,α regularity of the free boundary for the
N -dimensional porous media equation. Comm. Pure Appl. Math., 43(7):885 902, 1990.
4
[16] N. N. Carlson and K. Miller. Design and application of a gradient-weighted moving nite
element code. I. In one dimension. SIAM J. Sci. Comput., 19(3):728 765, 1998. 10
[17] N. N. Carlson and K. Miller. Design and application of a gradient-weighted moving nite
element code. II. In two dimensions. SIAM J. Sci. Comput., 19(3):766 798, 1998. 10
[18] M. Chipot and B. Lovat. Some remarks on nonlocal elliptic and parabolic problems. In
Proceedings of the Second World Congress of Nonlinear Analysts, Part 7 (Athens, 1996),
volume 30, pages 4619 4627, 1997. 11, 12
[19] M. Chipot and B. Lovat. On the asymptotic behaviour of some nonlocal problems. Positi-
vity, 3(1):65 81, 1999. 11, 12
[20] M. Chipot and J.-F. Rodrigues. On a class of nonlocal nonlinear elliptic problems. RAIRO
Modél. Math. Anal. Numér., 26(3):447 467, 1992. 11, 12
[21] P. G. Ciarlet. The nite element method for elliptic problems, volume 40 of Classics
in Applied Mathematics. Society for Industrial and Applied Mathematics (SIAM), Phila-
delphia, PA, 2002. Reprint of the 1978 original [North-Holland, Amsterdam; MR0520174
(58 #25001)]. 15
[23] P. Daskalopoulos and R. Hamilton. Regularity of the free boundary for the porous medium
equation. J. Amer. Math. Soc., 11(4):899 965, 1998. 4
[24] C. de Boor. Good approximation by splines with variable knots. II. In Conference on the
Numerical Solution of Differential Equations (Univ. Dundee, Dundee, 1973), pages 12 20.
Lecture Notes in Math., Vol. 363. Springer, Berlin, 1974. 10
108
Método dos Elementos Finitos para problemas com Fronteiras Livres
[34] C. Ebmeyer. Error estimates for a class of degenerate parabolic equations. SIAM J. Numer.
Anal., 35(3):1095 1112 (electronic), 1998. 8
[35] C. Ebmeyer and W. B. Liu. Finite element approximation of the fast diffusion and the
porous medium equations. SIAM J. Numer. Anal., 46(5):2393 2410, 2008. 8
[37] R. Eymard, T. Gallouët, R. Herbin, and A. Michel. Convergence of a nite volume scheme
for nonlinear degenerate parabolic equations. Numer. Math., 92(1):41 82, 2002. 13
[41] E. Henriques and J. M. Urbano. Intrinsic scaling for PDE's with an exponential nonlinearity.
Indiana Univ. Math. J., 55(5):1701 1721, 2006. 7
[42] M. A. Herrero and J. L. Vázquez. The one-dimensional nonlinear heat equation with ab-
sorption: regularity of solutions and interfaces. SIAM J. Math. Anal., 18(1):149 167, 1987.
5
[43] W. Huang, Y. Ren, and R. D. Russell. Moving mesh partial differential equations (MMPDES)
based on the equidistribution principle. SIAM J. Numer. Anal., 31(3):709 730, 1994. 10,
61
109
Método dos Elementos Finitos para problemas com Fronteiras Livres
[44] W. Huang and R. D. Russell. Moving mesh strategy based on a gradient ow equation for
two-dimensional problems. SIAM J. Sci. Comput., 20(3):998 1015 (electronic), 1999. 10,
62
[48] R. Kershner. The behavior of temperature fronts in media with nonlinear heat conductivity
under absorption. Vestnik Moskov. Univ. Ser. I Mat. Mekh., (5):44 51, 1978. 6
[49] B. F. Knerr. The porous medium equation in one dimension. Trans. Amer. Math. Soc.,
234(2):381 415, 1977. 4
[50] B. F. Knerr. The behavior of the support of solutions of the equation of nonlinear heat
conduction with absorption in one dimension. Trans. Amer. Math. Soc., 249(2):409 424,
1979. 5
[51] Y. Ko. Regularity of the interface for the porous medium equation. Electron. J. Differential
Equations, pages No. 68, 12 pp. (electronic), 2000. 4
[52] O. A. Ladyzhenskaya, V. A. Solonnikov, and N. N. Ural0 tseva. Linear and Quasi-Linear Equa-
tions of Parabolic Type. Translations of Mathematical Monographs. American Mathematical
Society, Providence, 1968. 2, 43, 50
[53] O. A. Ladyzhenskaya and N. N. Ural0 tseva. Linear and quasilinear elliptic equations. Trans-
lated from the Russian by Scripta Technica, Inc. Translation editor: Leon Ehrenpreis. Aca-
demic Press, New York, 1968. 39
[54] S. Lian, W. Gao, C. Cao, and H. Yuan. Study of the solutions to a model porous medium
equation with variable exponent of nonlinearity. J. Math. Anal. Appl., 342(1):27 38, 2008.
7
[55] J.-L. Lions. Quelques méthodes de résolution des problèmes aux limites non linéaires.
Dunod, Paris, 1969. 3
[56] J. N. Lyness and R. Cools. A survey of numerical cubature over triangles. In Mathematics
of Computation 1943 1993: a half-century of computational mathematics (Vancouver,
BC, 1993), volume 48 of Proc. Sympos. Appl. Math., pages 127 150. Amer. Math. Soc.,
Providence, RI, 1994. 32
[57] R. Marlow, M. E. Hubbard, and P. K. Jimack. Moving mesh methods for solving parabolic
partial differential equations. Comput. & Fluids, 46:353 361, 2011. 11
[58] K. Miller and R. N. Miller. Moving nite elements. I. SIAM J. Numer. Anal., 18(6):1019 1032,
1981. 9
[59] R. H. Nochetto and C. Verdi. Approximation of degenerate parabolic problems using nume-
rical integration. SIAM J. Numer. Anal., 25(4):784 814, 1988. 7
110
Método dos Elementos Finitos para problemas com Fronteiras Livres
[61] O. A. Ole nik, A. S. Kala²inkov, and Yu -Lin0 ƒºou. The Cauchy problem and boundary
problems for equations of the type of non-stationary ltration. Izv. Akad. Nauk SSSR. Ser.
Mat., 22:667 704, 1958. 3
[64] M. E. Rose. Numerical methods for ows through porous media. I. Math. Comp.,
40(162):435 467, 1983. 7
[65] J. Rulla and N. J. Walkington. Optimal rates of convergence for degenerate parabolic
problems in two dimensions. SIAM J. Numer. Anal., 33(1):56 67, 1996. 8
[66] L. F. Shampine and M. W. Reichelt. The MATLAB ODE suite. SIAM J. Sci. Comput.,
18(1):1 22, 1997. Dedicated to C. William Gear on the occasion of his 60th birthday.
30, 65
[70] M. R. Sidi Ammi and D. F. M. Torres. Numerical analysis of a nonlocal parabolic problem
resulting from thermistor problem. Math. Comput. Simulation, 77(2-3):291 300, 2008. 12
[71] J. Simsen and J. Ferreira. A global attractor for a nonlocal parabolic problem.
http://ptmat.fc.ul.pt/arquivo/docs/preprints/pdf/2013/preprint_Jorge_Ferreira.pdf,
March 2013. 12
[72] V. Thomée. Galerkin nite element methods for parabolic problems, volume 25 of Sprin-
ger Series in Computational Mathematics. Springer-Verlag, Berlin, second edition, 2006.
2, 15, 29, 30, 45, 51
[73] J. L. Vázquez. The porous medium equation. Oxford Mathematical Monographs. The
Clarendon Press Oxford University Press, Oxford, 2007. Mathematical theory. 6
111
Método dos Elementos Finitos para problemas com Fronteiras Livres
[74] A. Wacher and I. Sobey. String gradient weighted moving nite elements in multiple di-
mensions with applications in two dimensions. SIAM J. Sci. Comput., 29(2):459 480 (elec-
tronic), 2007. 10
[75] D. Wei and L. Lefton. A priori Lρ error estimates for Galerkin approximations to porous
medium and fast diffusion equations. Math. Comp., 68(227):971 989, 1999. 8
[76] Z. Yin and Q. Xu. A fully discrete symmetric nite volume element approximation of
nonlocal reactive ows in porous media. Math. Probl. Eng., pages Art. ID 175904, 7, 2013.
12
112
Método dos Elementos Finitos para problemas com Fronteiras Livres
Índice de notações
Geral
n denota o vetor normal unitário exterior à fronteira
∗
* denota uma convergência fraca estrela
p0 = p
p−1 é o conjugado de p > 1
|α| =
Pn
i=1 αi para α= (α1 , . . . , αn ) ∈ Nn
sign(a) denota o sinal de a (sign(a) = 1, se a > 0 e sign(a) = −1, se a < 0).
Acrónimos
113
Método dos Elementos Finitos para problemas com Fronteiras Livres
Conjuntos
Ω̄ denota o fecho de Ω
ΩT = Ω×]0, T ]
Γ = ∂Ω denota a fronteira de Ω
ΓT = ∂Ω × [0, T ]
Operadores
d
X ∂ui
div u = é o divergente de u
i=1
∂xi
∂u ∂u
∇u = ,..., denota o vetor gradiente de u
∂x1 ∂xd
d
X ∂2u
∆u = é o laplaciano de u
i=1
∂x2i
∂u
ut = denota a derivada ( fraca ou forte ) de u em relação à variável temporal
∂t
114
Método dos Elementos Finitos para problemas com Fronteiras Livres
∂ |α| u
denota a derivada forte de ordem α de u
∂xα
1
1
. . . xα
n
n
Rh u é a projeção de Ritz de u
Normas e seminormas
p1
k
X
kukW k,p (Ω) = kDα ukpLp (Ω) , para 1 ≤ p < ∞
|α|=0
k
X
kukW k,∞ (Ω) = ess sup |Dα u|
Ω
|α|=0
p1
X
|u|W k,p,j (Ω) = kDα ukpLp (Ω) , para 1 ≤ p < ∞ e 0 ≤ j ≤ k.
|α|=j
R p1
b
kukLp (a,b;X) = a
ku(t)kpX dt < ∞, para 1 ≤ p < ∞
p1
k j
X d u
kukW k,p (a,b;X) = k j kLp (a,b;X)
j=1
dt
Espaços de funções
115
Método dos Elementos Finitos para problemas com Fronteiras Livres
Rb
Lp (a, b; X) = {u :]a, b[→ X| u é mensurável e a
kukpX dt < ∞}
Rb
W k,p (a, b; X) = {u ∈ W k,p (Ω)| u é mensurável e a
kukpX dt < ∞}
0
W −k,p (Ω) é o espaço dual de W0k,p (Ω)
s
X
δs
Shr = {W : [0, +∞[→ Shr | W|In = tj χj , χj ∈ Shr }
j=0
116
Método dos Elementos Finitos para problemas com Fronteiras Livres
Índice Remissivo
convergência método
forte, 16 dos elementos nitos de Galerkin, 1
fraca, 18 elementos nitos descontínuo, 29
fraca estrela, 18 Euler, 29
corolário Euler regressivo, 92
Brouwer, 25 Faedo-Galerkin, 74
ponto xo, 52
desigualdade malha, 26
Cauchy, 15 regular, 26
Gronwall, 25
Hölder, 20 ordem de convergência, 28
inversa, 55
projeção
Minkowski, 20
de Ritz, 28
multiplicativa, 82, 86
no espaço e tempo, 29
Poincaré, 23
no tempo, 29
Young, 15
ortogonal, 28
domínio, 17
propriedade
equação de extinção em tempo nito, 5
em meios porosos, 3 de localização estável, 81
reação-difusão, 1 de propagação com velocidade nita, 3
espaço de tempo de espera, 4
de Bochner, 23
sistema
de elementos nitos, 26, 45, 63
não local, 11
de Lebesgue, 20
reação-difusão, 1
de Sobolev, 21
solução
dual, 18
forte, 79
fronteira fraca, 35, 36, 61, 63, 74
Lipschitz-contínua, 18 semidiscreta, 89
livre, 3 totalmente discreta, 46, 92
função
taxa, 81
de energia global, 85
teorema
de energia local, 81
Banach-Alaoglu, 24
Hölder-contínua, 17
Brouwer, 25
Lipschitz-contínua, 17
do valor médio, 41
monitor, 62
Green, 25
integrador, 29 Gronwall, 25
interface, 3
velocidade
laplaciano, 20 normal, 60
discreto, 28 vetor normal, 28
lema
Aubin-Lions, 25
117
Método dos Elementos Finitos para problemas com Fronteiras Livres
118