本文对傅里叶赝谱的Dirichlet BC强加

计算科学 流体动力学 边界条件 谱法
2021-12-05 15:45:35

我试图实现论文“将傅里叶伪谱方法应用于瑞利-贝纳德对流的狄利克雷边界条件”中的算法。

我很难理解施加边界条件的方式。

作者将无滑移(上下边界)边界条件改写为

qf~,pq=0p

在哪里pq是水平和垂直波数。

我们如何强加它?

1个回答

我不是该领域的专家,但使用 DCT(或 DST)是使用 FFT 施加 Dirichlet(或 Neumann)边界条件的标准方法。

在 FFTW 包的文档中有很好的解释

http://www.fftw.org/fftw3_doc/Real-even_002fodd-DFTs-_0028cosine_002fsine-transforms_0029.html#Real-even_002fodd-DFTs-_0028cosine_002fsine-transforms_0029

我希望它有所帮助。