我正在努力将广义加权残差方法与 RK4 结合起来。GWRM 部分将 PDE 分解到光谱空间，其中未知数是 Chebyshev 系数 a_k。但是，我很难看到在这种情况下如何包含边界条件。在其他光谱方法中，包括物理网格，因此边界条件可以显式设置或仅包含在 Chebyshev 差分矩阵中。
这里是我正在解决的 ODE 的简短推导。有人对如何包含边界条件有任何想法吗？
请记住，A、b和c都是向量。素数表示第一个和项除以 2。
P.s 得到的方程是可以用 rk4 离散化的 ODE。
这是我目前对 BCs 如何实现的理解，但是随着每个时间步骤，解决方案越来越远离真正的边界条件。
可以用最高模式 K 和 K-1 处的切比雪夫系数代替边界方程，
0
答案（75 ％ 确定？）是，由于频谱空间中没有明确的边界条件，因此不可能采用明确的时间积分方案。基函数必须满足边界条件，或者需要明确设置边界条件。
为了使用 GWRM 求解 PDE，您需要在频谱分解中包含时域并求解一组线性 / 非线性代数方程，或者使用隐式时间积分方案，如 Backward Euler 或隐式 RK4。
隐式方法而不是显式方法起作用的原因是，在隐式方法中，下一个时间步的切比雪夫系数出现在方程的两侧。因此，您可以用最高模式代替边界条件，并迭代直到切比雪夫系数的下一步满足 PDE 和边界条件。
本站系公益性非盈利分享网址，本文来自用户投稿，不代表边看边学立场，如若转载，请注明出处
评论列表（73条）