- Em quais casos não é possível calcular o limitante superior para o erro de truncamento no método das diferenças finitas?
Passo 1
O enunciado pergunta claramente: “Em quais casos não é possível calcular o limitante superior para o erro de truncamento no método das diferenças finitas?”\n\nInformações relevantes que o enunciado implica e que vou usar depois:
- O “erro de truncamento” é um termo local que aparece quando uma aproximação por diferenças finitas é obtida truncando uma expansão (por ex. uma série de Taylor).
- Para calcular um limitante superior desse erro normalmente é preciso conhecer a expressão do termo de resto (como função de derivadas da solução) e ter um limite superior (uma constante) para as derivadas envolvidas na expressão.
- O tamanho do passo
e a ordem do esquema (expoente de no erro) também entram na fórmula do limitante.
Passo 2
Explicando como normalmente se obtém um limitante superior (aula rápida):\n\nQuando derivamos um esquema de diferenças finitas usamos a expansão de Taylor com resto na forma de Lagrange. A forma geral do resto é
onde
é algum ponto entre
e
. Aqui:
-
é a derivada de ordem da função que estamos aproximando; -
traz potências do passo quando .\n\nExemplo concreto muito usado: aproximação central para a primeira derivada:
Para obter um limitante superior do termo de erro basta ter um
tal que
em todo o intervalo de interesse; então
Repare que tudo funcionou porque conhecemos (ou assumimos) que a derivada de ordem necessária existe e é limitada.
Passo 3
Agora, quando isso deixa de ser possível? Vou listar os casos e explicar por que em cada um não dá para calcular um limitante superior teórico do erro de truncamento.
- A função não tem derivada da ordem exigida (falta de regularidade)
- Se o esquema exige
e a função só tem derivivadas até ordem menor que ou essas derivadas não existem em pontos do domínio, a fórmula do resto com não é válida. Exemplo: não tem , então um esquema que precisa de numa vizinhança de não tem limitante superior clássico.
- A derivada existe mas é não-limitada (singularidade)
- Se
existe porém é não limitada ou tende a infinito em algum ponto do domínio (p.ex. comportamento perto de ), não se consegue encontrar uma constante finita para fazer a estimativa , logo não há limitante finito do tipo .
- A solução é desconhecida e não há estimativas a priori
- Em muitos problemas (por ex. resolver uma EDP numérica) a função exata
não é conhecida em forma fechada. Para obter um limitante do tipo é preciso um limite para . Se não existem teoremas de regularidade que garantam esse limite (ou se não se pode provar que é limitada), então não dá para calcular o limitante superior teórico.
- Dados discretos sem modelo contínuo
- Quando só se tem valores
em pontos (dados experimentais) sem saber a função contínua subjacente nem sua regularidade, não é possível garantir uma constante para as derivadas; assim não há como calcular um limitante teórico.
- Presença de descontinuidades ou singularidades internas (choques, fronteiras de material, cantos)
- Se a solução apresenta descontinuidade ou esquina na derivada em algum lugar do domínio (p.ex. solução de conservação com choque), as hipóteses para a expansão de Taylor falham localmente e não se pode dar um bound uniforme do erro de truncamento baseado em derivadas.
- Esquemas não lineares ou termos de erro dependentes da solução de modo intratável
- Para alguns esquemas não lineares o termo de truncamento envolve combinações de derivadas e da própria solução de forma que não se consegue isolar um
conhecido; sem estimativas da solução, o limitante não é computável.
- Mal condicionamento / falta de estabilidade para transformar erro local em global
- Mesmo quando existe um limitante local do erro de truncamento, para o erro global é preciso também estabilidade do método. Se o sistema/passi não é conhecido como estável (ou a constante de estabilidade é desconhecida e possivelmente grande), então não se consegue um limitante útil do erro total — formalmente um limitante para o erro de truncamento local pode existir, mas não há como transformá-lo em um limitante do erro global.
Passo 4
Exemplos práticos e curtos para fixar a ideia (sou teu amigo explicando na boa):\n\n- Exemplo 1 (não dá): usar esquema de segunda derivada central perto de
para
. A quarta derivada de
diverge em
, então não existe
finito para
e o limitante do tipo
não pode ser calculado.\n\n- Exemplo 2 (dá): para
em
, todas as derivadas existem e são limitadas por
. Para o esquema central da segunda derivada (erro
) podemos tomar
e obter um limitante explícito
.\n\nQuando não se consegue o bound teórico, recomendações práticas:
- Fazer estudo numérico de convergência (rodar com
e estimar ordem); - Usar estimadores a posteriori (residual-based) que não exigem conhecer derivadas da solução exata;
- Procurar teoremas de regularidade (por ex. para EDP elípticas) que garantam bounds para derivadas da solução; se não existir, aceitar que só estimativas numéricas são possíveis.
Passo 5
Recapitulando de forma objetiva (antes da resposta final): para calcular um limitante superior do erro de truncamento é necessário poder escrever o termo de resto (geralmente em função de derivadas de certa ordem) e possuir um limite superior conhecido para essas derivadas. Sempre que essa informação faltar — por falta de regularidade, por singularidades, por desconhecimento da solução ou por dependência não tratável da solução no termo de erro — não é possível obter o limitante superior teórico.
Resposta
- Quando a função não tem derivadas da ordem exigida na expansão (falta de regularidade).
- Quando as derivadas necessárias existem mas são não-limitadas (singularidades), impossibilitando encontrar uma constante
finita. - Quando a solução exata é desconhecida e não existem estimativas a priori (teoremas de regularidade) que forneçam bounds para as derivadas.
- Quando só há dados discretos sem modelo contínuo ou sem informação sobre a regularidade da função.
- Quando há descontinuidades ou choques na solução que invalidam a expansão de Taylor local.
- Quando o termo de truncamento depende de forma intratável da própria solução (esquemas não lineares) e não se pode isolar uma constante limitante.
- Quando, embora exista um bound local, falta informação de estabilidade (ou a constante de estabilidade é desconhecida), impedindo um limitante útil para o erro global.