@xjwang
CFDRFG的附录code还有个问题,求解$K$,不是按照论文中的原理计算的。他先用$PQ$没有代入$S(f)$的值,算了一组$K$,然后所有坐标点都用了同一组$K$。而根据理论,应该是要根据每个坐标点的非均匀$PQ$,求对应该点的$K$。DSRFG如果生成非均匀,理论上也要这么计算。
739858e5-ac4f-4736-94bd-9f4ec9f07a9d-image.png
关于CDRFG方法,他们的code里是对每一个点都生成一组PQK
最近看到一篇JFM文章《Guo, H., et al., An efficient and low-divergence method for generating inhomogeneous and anisotropic turbulence with arbitrary spectra. Journal of Fluid Mechanics, 2023. 970: p. A2.》,用的就是三维能谱,不是频谱,生成非均匀湍流场也是low-divergence。