有限差分法を使用してPDEを解くときに境界条件を選択する方法を説明するのに役立つリソースを見つけようとしています。
私が現在アクセスできるすべての本とメモは、同様のことを言っています:
境界が存在する場合の安定性を管理する一般規則は、入門テキストとしては複雑すぎます。彼らは洗練された数学的な機械を必要とします
(A.イセル微分方程式の数値解析の最初のコース)
たとえば、移流方程式に2段階のリープフロッグ法を実装しようとする場合:
MATLABを使用して
M = 100; N = 100;
mu = 0.5;
c = [mu 0 -mu];
f = @(x)(exp(-100*(x-0.5).^2));
u = zeros (M, N);
x = 1/(M+1) * (1:M);
u(:,1) = f(x);
u(:,2) = f(x + mu/(M+1));
for i = 3:N
hold off;
u(:,i) = conv(u(:,i-1),c,'same') + u(:,i-2);
plot(x, u(:,i));
axis( [ 0 1 0 2] )
drawnow;
end
このソリューションは、境界に到達するまでうまく動作し、突然突然動作が異常になります。
このような境界条件の処理方法はどこで学べますか?