如何处理 PDE 过实线

计算科学 pde 边界条件
2021-12-04 21:33:17

我定义了一个 PDER,对此我没有确切的解决方案,我要用有限差分来近似它,所以我需要输入一些 BC。

任何人都可以提出任何很好的参考来处理无限间隔作为 BC 吗?

可能的尝试:

  1. 输入[M,M]作为大 M 的域并尝试猜测确切的解x=±M.

  2. 因为我的IC的对称性(logcoshx),识别左右边界,即u(M,t)=u(M,t).

  3. 尝试一些 Neumann 条件。

  4. 尝试周期性边界。

语境

感兴趣的 PDE 是

ut=a1uxx+a2Φ(t)ux2
u(0,x)=logcoshx, 定义在(x,t)R×[0,1],Φ(t)是标准高斯的 CDF 和ai是常数。

这没有明确的分析解决方案。

1个回答

除了 PDE 本身之外,您还必须有边界条件。例如,对于您的情况,常见的边界条件是u(t,x)0为了x±对于任何t. 我见过其他情况是这种形式的情况x2nu(t,x)0为了x±.

我不知道可以直接求解此类方程的有限差分方法(尽管有限元方法有“无限元”公式)。但是您可以引入坐标变换,例如(1,1)yx(,)它将您的方程转换为在实线的有限子集上提出的方程。

为了避免一些复杂性,让我假设您最初的问题是在半线上提出的x(0,)并且我们想在区间上将其转换为 1y(0,1). 然后我们可以(在许多其他选项中)选择y=1ex,产生x=yxy=exy=(1y)y.

然后,您将替换所有关于x通过关于衍生品y使用这个表达式。这会产生具有空间可变系数的 PDE(取决于y),您可以对其应用有限差分法。