|||
生成湍流速度场的时候通常要使速度场满足某种分布(http://blog.sciencenet.cn/home.php?mod=space&uid=117333&do=blog&id=693433)。要做到这一点,最方便的做法就是在傅里叶空间完成。
实际实现的时候可以采取这样的方法:
1. 生成实空间中的三维随机场。
2. 傅里叶变换。
3. 注意到傅里叶变换得到的场只有八分之一(一个卦限)独立,对于numpy,三维数组的八个顶角都是原点,计算k矢量的模的时候各卦限计算到相应顶角的距离。
4. 每个点乘以k矢量模的某个幂次,使得速度场满足某种分布。
5. 傅里叶反变换。
6. 将实空间速度场归一化。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-7-18 14:25
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社