CFD中文网

@李东岳 李老师，问题算是解决了，正如楼上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文章说这个参数很重要，给了三种设置的方法，但是我一直没引起重视。。。。您这个参数是如何来设置的呢？谢谢！

@李东岳 入口速度是在平均速度基础上加的扰动，这些扰动就是通过那个算法实现的，因此在每个时间步入口速度场的扰动时一直存在的。

@xjwang 同时，我在DSRFG方法中采用的参数都是根据自己计算工况进行的设置，比如x,y,z坐标数据、湍流积分尺度l等，这是否可以理解为您指出的“因为参数选错了，而引起了divergence-free问题”？之前一直在网格和时间步长设置方面找原因，没有考虑过本身生成的速度场是否问题。呜呜0.0！

@xjwang 辛苦辛苦^ ^，这是计算得到的瞬时速度等值线图，可以看出脉动速度在下游是衰减的（因为一直通过瞬时速度波动情况来观察下游湍流情况，不知这样判断是否有问题？）。另外关于您提到的湍流积分尺度参数，是指huang文章中说的turbulence integral length scales Lu/Lv/Lw么？我在想我是否也是这个问题导致的？

@xjwang 感谢回复，我的湍流入口方法也是采用DSRFG方法(Huang et al,2010)，速度入口边界，也是把脉动速度场定义在入口边界面上了，计算过程中通过监测不同轴向截面速度的等值线图，发现在入口附近速度场有小幅值脉动，在远离入口处速度场的脉动几乎消失（感觉和你说的沿着顺风向，脉动功率谱密度减小的问题类似）。而理想的结果是否应该是由于入口扰动的持续作用，在下游能一直维持较强脉动的湍流场。不知您当时是如何解决的？

@李东岳 李老师好，计算域是完全的矩形，即三维槽道流。

@xjwang 朋友你好，我也遇到了同样的问题，采用LES计算，时变的速度入口数据导入后，随着非稳态进行，脉动速度在向下游传递过程中很快减弱至消失，不知是否已经解决？

@cccrrryyy 好的，谢谢

@luoluo 感谢回复，我有一个疑问，由于脉动速度的长时间相关总是等于零，经过足够长的时间后初始场的随机状态对湍流场的发展几乎没有影响，因此对初始场添加扰动的原因还有些疑问，为什么要添加这个扰动，如果不加扰动，直接用RANS计算得到的时间结果作为初场，进出口也给周期边条，采用大涡模拟能否自发生成湍流？之前算过，好像是产生不了湍流的，只是其中的原因一直不清楚。