1. 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 h e a ordem do esquema (expoente de h 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 é
R_n=\frac{f^{(n+1)}(\xi)}{(n+1)!}(x-x_0)^{n+1},
onde \xi é algum ponto entre x e x_0. Aqui:
  • f^{(n+1)} é a derivada de ordem n+1 da função que estamos aproximando;
  • (x-x_0)^{n+1} traz potências do passo h quando x=x_0\pm h.\n\nExemplo concreto muito usado: aproximação central para a primeira derivada:
\frac{f(x+h)-f(x-h)}{2h}=f'(x)-\frac{h^2}{6}f^{(3)}(\xi).
Para obter um limitante superior do termo de erro basta ter um M tal que |f^{(3)}(t)|\le M em todo o intervalo de interesse; então
\left|\text{erro}\right|\le \frac{h^2}{6}M.
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.
  1. A função não tem derivada da ordem exigida (falta de regularidade)
  • Se o esquema exige f^{(m)} e a função só tem derivivadas até ordem menor que m ou essas derivadas não existem em pontos do domínio, a fórmula do resto com f^{(m)}(\xi) não é válida. Exemplo: f(x)=|x| não tem f''(0), então um esquema que precisa de f'' numa vizinhança de 0 não tem limitante superior clássico.
  1. A derivada existe mas é não-limitada (singularidade)
  • Se f^{(m)} existe porém é não limitada ou tende a infinito em algum ponto do domínio (p.ex. comportamento f^{(m)}(x)\sim 1/\sqrt{x} perto de 0), não se consegue encontrar uma constante M finita para fazer a estimativa |f^{(m)}|\le M, logo não há limitante finito do tipo C h^p.
  1. 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 u não é conhecida em forma fechada. Para obter um limitante do tipo C h^p é preciso um limite para u^{(m)}. Se não existem teoremas de regularidade que garantam esse limite (ou se não se pode provar que u^{(m)} é limitada), então não dá para calcular o limitante superior teórico.
  1. Dados discretos sem modelo contínuo
  • Quando só se tem valores f(x_i) em pontos (dados experimentais) sem saber a função contínua subjacente nem sua regularidade, não é possível garantir uma constante M para as derivadas; assim não há como calcular um limitante teórico.
  1. 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.
  1. 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 M conhecido; sem estimativas da solução, o limitante não é computável.
  1. 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 x=0 para f(x)=\sqrt{x}. A quarta derivada de \sqrt{x} diverge em 0, então não existe M finito para |f^{(4)}| e o limitante do tipo C h^2 não pode ser calculado.\n\n- Exemplo 2 (dá): para f(x)=\sin x em [0,\pi], todas as derivadas existem e são limitadas por 1. Para o esquema central da segunda derivada (erro \propto h^2 f^{(4)}(\xi)/12) podemos tomar M=1 e obter um limitante explícito \le h^2/12.\n\nQuando não se consegue o bound teórico, recomendações práticas:
  • Fazer estudo numérico de convergência (rodar com h,h/2,h/4 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 M 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.

Exercícios relacionados