@李东岳 李老师，问题算是解决了，正如楼上xjwang所说的那样，不应该通过监测瞬时速度云图的方式来判断，后来我通过数次调整积分尺度的值，并在下游检测不同位置点的能谱情况，发现是可以得到湍流场的，虽然能量有些衰减，但就目前来看，结果是可以接受的。

@sunss 感谢分享。我做的是槽道流，能谱是通过槽道流实验得到的（沿轴向的一维谱），通过DSRFG方法生成了入口速度场，得到的入口速度场能谱分布与实验能谱拟合的较好，但是在计算中会出现我上面提到下游脉动减弱的问题。下图是我还原Huang文章中的卡门谱结果，整体看起来方法的实现应该是没有问题。

@李东岳 还没有，还在找原因@@

@sunss 感谢分享，但是根据huang2010论文中指出该方法严格保证了无散度条件，并且给出了证明过程。

@xjwang 你好，我的研究不是ABL，属于槽道流，但是想用DSRFG这个方法生成湍流入口用于LES计算。

@xjwang 朋友，是用fluent画的图。我后来更新了结果，在前面的帖子里放上了结果，你看下结果如何？

@李东岳 如李老师所言：“像是这种速度波动越往后越少的，类似开始有湍流后来无湍流”。我有个疑惑，采用中心差分格式的扩散项可以将扰动向四周传递，迎风格式的对流项具有迁移性，即可以将扰动向下游传递。从上面来看，入口的扰动应该是可以持续向下游传递的。如果几何轴向长度非常长的话，在很远的下游是不是也会保持一定湍流状态？还是由于壁面摩擦力及流体粘性对扰动有削减作用，在下游一直衰减甚至脉动全无，就是这种入口给定方法本身就会出现的现象呢？

@xjwang 根据您的建议，我修改了湍流积分尺度，由于没有实验，这个量我是通过湍动能和耗散率来估算的。重新进行计算，目前能够看出下游流场表现出湍流状态。这个结果和您之前的计算发展趋势相似么，以说明我目前入口给定结果的正确性。

@李东岳 李老师，麻烦能说一下internalField的作用么，我不是很了解，谢谢

@xjwang 感谢回复，我的计算域如图所示，几何尺寸：X=1.2m、Y=0.2m、Z=0.4m。展向和轴向为均匀网格。空气沿x轴正向进入，平均速度是1.2m/s。入口采用速度入口边界条件，出口为压力出口边界条件，展向是周期性边界条件（即Z方向），上部和底部为壁面，采用无滑移边界条件。三个方向网格节点：64X64X128，垂直壁面方向采用非均匀网格，首层网格高度△y=0.0001m，保证Yplus<1。目前程序中湍流积分尺度设置为Ls=0.014（大致估算的，目前看来是有问题的）。Huang文章说这个参数很重要，给了三种设置的方法，但是我一直没引起重视。。。。您这个参数是如何来设置的呢？谢谢！